• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    軸流式噴水推進(jìn)器在來(lái)流含氣條件下的空化特性分析

    2022-08-29 11:49:38潘中永劉月偉
    關(guān)鍵詞:含氣率含氣輪緣

    潘中永,張 帥,劉月偉

    (江蘇大學(xué) 國(guó)家水泵及系統(tǒng)工程技術(shù)研究中心,鎮(zhèn)江 212013)

    噴水推進(jìn)作為一種新型的船舶推進(jìn)技術(shù),利用噴水推進(jìn)泵的高速運(yùn)行從船底進(jìn)水流道吸入的水流,使水流從噴嘴射出,從而產(chǎn)生反作用力推動(dòng)船舶前進(jìn).噴水推進(jìn)器具有很多傳統(tǒng)螺旋槳沒(méi)有的優(yōu)點(diǎn),例如淺水區(qū)運(yùn)行效果好,高航速時(shí)推進(jìn)效率高,抗空化能力強(qiáng)等.因此,噴水推進(jìn)器在高性能船舶上的運(yùn)用變得越來(lái)越普遍.噴水推進(jìn)器實(shí)際運(yùn)行過(guò)程中的來(lái)流含氣狀況會(huì)對(duì)噴水推進(jìn)器空化特性造成影響.來(lái)流含氣狀況相比于純水,其氣液兩相流動(dòng)特性與氣泡直徑、氣體含量等因素均有不同程度的關(guān)系.噴水推進(jìn)系統(tǒng)空化會(huì)加速葉片表面剝蝕,從而影響推進(jìn)系統(tǒng)的推進(jìn)性能[1].文獻(xiàn)[2]通過(guò)試驗(yàn)方法觀察到噴水推進(jìn)泵葉輪內(nèi)垂直空泡渦(PCV)的形成過(guò)程,并發(fā)現(xiàn)PCV的脫落會(huì)間接造成泵揚(yáng)程急劇下降.文獻(xiàn)[3]在對(duì)噴水推進(jìn)滑行艇的研究中發(fā)現(xiàn)空化會(huì)破壞葉輪和導(dǎo)葉的非穩(wěn)定力,嚴(yán)重空化時(shí)空泡會(huì)侵蝕破壞泵體.文獻(xiàn)[4]基于高速成像技術(shù)觀察軸流式噴水推進(jìn)泵葉輪空化演變過(guò)程,發(fā)現(xiàn)空化最早出現(xiàn)在葉輪葉頂間隙處.文獻(xiàn)[5-6]使用高速攝影拍攝了混流式噴水推進(jìn)泵內(nèi)部空化演變規(guī)律,揭示了泵內(nèi)空化演變過(guò)程及翼型背面壓力脈動(dòng)規(guī)律.基于雷諾時(shí)均方程和均質(zhì)多相流模型,文獻(xiàn)[7]對(duì)噴水推進(jìn)軸流泵空化性能進(jìn)行仿真分析,發(fā)現(xiàn)空化會(huì)導(dǎo)致推力和扭矩陡然下降.文獻(xiàn)[8]基于數(shù)值模擬方法對(duì)噴水推進(jìn)泵葉頂泄漏渦空化進(jìn)行研究,發(fā)現(xiàn)葉頂附近監(jiān)測(cè)點(diǎn)壓力脈動(dòng)幅值隨空化程度加劇而增大.文獻(xiàn)[9]對(duì)混流式噴水推進(jìn)泵空化與渦相互作用機(jī)理進(jìn)行數(shù)值模擬研究,發(fā)現(xiàn)空化的發(fā)展加劇了泵內(nèi)渦流的產(chǎn)生和流動(dòng)不穩(wěn)定性.文獻(xiàn)[10]利用數(shù)值模擬方法對(duì)噴水推進(jìn)泵在性能陡降工況下的內(nèi)部空化流動(dòng)機(jī)理及渦與空化相互作用機(jī)制進(jìn)行深入分析.文獻(xiàn)[11]利用高速攝像技術(shù)捕捉到了噴水推進(jìn)泵葉輪在空化初生工況下的空化流動(dòng)形態(tài).文獻(xiàn)[12]對(duì)混流式噴水推進(jìn)器空化性能進(jìn)行數(shù)值模擬,發(fā)現(xiàn)空化易導(dǎo)致噴水推進(jìn)器流量、揚(yáng)程等性能參數(shù)明顯下降.針對(duì)上述文獻(xiàn)可知,國(guó)內(nèi)外對(duì)泵內(nèi)氣液兩相流對(duì)泵空化性能的影響進(jìn)行了大量的試驗(yàn)研究與數(shù)值模擬,但是對(duì)噴水推進(jìn)器內(nèi)部氣液兩相流動(dòng)研究較少,主要集中在純水介質(zhì)下的推進(jìn)泵性能和純水空化性能方面,同時(shí)大多采用直管替代彎曲進(jìn)水流道,很難真實(shí)模擬出實(shí)際航行時(shí)不同航速下的噴水推進(jìn)器性能.對(duì)噴水推進(jìn)器在實(shí)際航行時(shí)來(lái)流含氣造成的內(nèi)部氣液兩相流動(dòng)以及來(lái)流含氣后的空化流場(chǎng)狀況研究則更為匱乏.因此,文中以船用軸流式噴水推進(jìn)器為研究對(duì)象,對(duì)其在不同航速下純水介質(zhì)及不同氣相參數(shù)的氣液兩相介質(zhì)下的空化流動(dòng)定常與非定常特性展開數(shù)值模擬.

    1 數(shù)值計(jì)算模型

    1.1 模型參數(shù)

    噴水推進(jìn)器軸流泵的主要參數(shù)如表1.

    表1 噴水推進(jìn)泵參數(shù)表

    通過(guò)三維軟件Creo6.0對(duì)軸流式噴水推進(jìn)器模型進(jìn)行三維建模,如圖1.

    圖1 噴水推進(jìn)器三維模型Fig.1 3D model of intake duct and axial flow water-jet pump

    采用流場(chǎng)控制體模擬船底周圍水域并加裝在進(jìn)水流道下方,流場(chǎng)控制體尺寸為30D×10D×8D,其中D為噴水推進(jìn)泵葉輪外徑[13].裝配得到整體計(jì)算域,如圖2.

    圖2 噴水推進(jìn)器計(jì)算域Fig.2 Computational domain of water-jet propulsion

    1.2 多相流和空化模型

    目前對(duì)氣液兩相計(jì)算主要使用歐拉-歐拉法.在CFX中,歐拉-歐拉方法使用的模型有兩種:均相流模型和非均相流模型.

    均相流模型假設(shè)各相速度相同且相間無(wú)滑移,非均相流模型在考慮相間速度滑移的基礎(chǔ)上,還考慮了相間動(dòng)量傳遞和質(zhì)量傳遞等.在氣液兩相流數(shù)值模擬中,非均相流假設(shè)氣液兩相相互獨(dú)立,通過(guò)對(duì)每一相進(jìn)行單獨(dú)求解,充分考慮了氣液兩相之間的能量傳遞和速度滑移作用,更符合真實(shí)流動(dòng).因此,選用歐拉-歐拉非均相流模型進(jìn)行數(shù)值計(jì)算,相間傳遞采用粒子模型,動(dòng)量傳遞采用Schiller Naumann模型.空化模型采用CFX軟件默認(rèn)的ZGB空化模型.

    2 數(shù)值計(jì)算方法

    2.1 網(wǎng)格劃分

    利用ICEM軟件對(duì)噴水推進(jìn)器和流場(chǎng)控制體進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分,同時(shí)對(duì)進(jìn)水流道與流場(chǎng)控制體交界面處進(jìn)行局部網(wǎng)格加密,如圖3.

    圖3 進(jìn)水流道及流場(chǎng)控制體網(wǎng)格Fig.3 Mesh of intake duct and water volume

    噴水推進(jìn)泵采用TurboGrid進(jìn)行自動(dòng)結(jié)構(gòu)化網(wǎng)格劃分,并對(duì)葉片周圍壁面和輪緣間隙區(qū)進(jìn)行局部網(wǎng)格加密,保證壁面y+符合湍流模型求解要求,如圖4.

    2.2 網(wǎng)格無(wú)關(guān)性驗(yàn)證

    對(duì)噴水推進(jìn)器推力和流量進(jìn)行數(shù)值模擬,當(dāng)計(jì)算域網(wǎng)格數(shù)量在3 302 048時(shí),再增加網(wǎng)格數(shù)量,推力變化小于1%,符合網(wǎng)格無(wú)關(guān)性驗(yàn)證要求,具體方案如表2.考慮計(jì)算資源及求解精度,最終選擇網(wǎng)格劃分方案3.

    表2 網(wǎng)格無(wú)關(guān)性檢驗(yàn)

    2.3 邊界條件

    采用CFX 2019R3軟件進(jìn)行數(shù)值模擬,葉輪設(shè)為旋轉(zhuǎn)域,轉(zhuǎn)速為5 500 r/min;噴嘴、導(dǎo)葉、進(jìn)水流道和流場(chǎng)控制體設(shè)為靜止域,進(jìn)水流道內(nèi)部嵌入的驅(qū)動(dòng)軸設(shè)為相對(duì)旋轉(zhuǎn)運(yùn)動(dòng),計(jì)算域內(nèi)固體壁面條件均設(shè)為無(wú)滑移壁面.流場(chǎng)控制體模擬船底水域,控制體進(jìn)口設(shè)為速度進(jìn)口,控制體出口采用自由出流;噴嘴出口通向大氣,采用靜壓出口,壓力值為0.1 MPa.

    在材料中新建和設(shè)置空泡Vapour相物性參數(shù),純水空化流場(chǎng)模擬中,設(shè)置初始進(jìn)水流道進(jìn)口空泡、液相體積分?jǐn)?shù)分別為0和1.來(lái)流含氣空化流場(chǎng)模擬中,設(shè)置初始進(jìn)水流道進(jìn)口空泡體積分?jǐn)?shù)為0,氣液兩相體積分?jǐn)?shù)之和為1.假設(shè)流動(dòng)過(guò)程中溫度保持不變,該溫度下流體的飽和蒸汽壓力為3 574 Pa.

    以雷諾時(shí)均N-S方程作為基本控制方程,采用SSTk-ω湍流模型,殘差收斂精度設(shè)為10-5.定常計(jì)算下動(dòng)靜交界面設(shè)為Frozen Rotor,對(duì)流項(xiàng)設(shè)為High Resolution.非定常計(jì)算下將進(jìn)水流道與葉輪、葉輪與導(dǎo)葉動(dòng)靜交界面改為“Transient Rotor Stator”模式.由穩(wěn)態(tài)計(jì)算調(diào)整為瞬態(tài)計(jì)算,計(jì)算總時(shí)長(zhǎng)T=0.08 727 s,共8個(gè)周期,每3°算一步,對(duì)應(yīng)ΔT=0.000 090 91 s.

    2.4 噴水推進(jìn)器推力性能試驗(yàn)

    噴水推進(jìn)器性能實(shí)船試驗(yàn)系統(tǒng)結(jié)構(gòu)復(fù)雜,部件較多,組裝復(fù)雜,成本高,并且穩(wěn)定性差,一旦性能達(dá)不到預(yù)期效果就會(huì)導(dǎo)致試驗(yàn)資源的過(guò)度浪費(fèi).因此,靜水試驗(yàn)臺(tái)通過(guò)改變噴泵轉(zhuǎn)速n來(lái)測(cè)試噴水推進(jìn)器推力T的性能,具有結(jié)構(gòu)簡(jiǎn)單,加工成本低,組裝簡(jiǎn)單,試驗(yàn)周期短等優(yōu)勢(shì),試驗(yàn)臺(tái)和試驗(yàn)結(jié)果如圖5.由于柴油機(jī)組的問(wèn)題,試驗(yàn)過(guò)程中未能實(shí)現(xiàn)較高轉(zhuǎn)速,最高只達(dá)到額定轉(zhuǎn)速.連續(xù)進(jìn)行兩次試驗(yàn)以保證試驗(yàn)結(jié)果的準(zhǔn)確性,隨著轉(zhuǎn)速增加,流量增大,推力增加,初步分析由于仿真忽略了軸承摩擦等機(jī)械損失,導(dǎo)致各個(gè)轉(zhuǎn)速下的仿真推力值高于試驗(yàn)推力值,但試驗(yàn)和仿真推力值隨轉(zhuǎn)速變化趨勢(shì)吻合度較高,尤其是在轉(zhuǎn)速達(dá)到額定轉(zhuǎn)速5 500 r/min時(shí),試驗(yàn)推力與仿真推力誤差約為3%,證明噴水推進(jìn)器仿真模型能夠?qū)ζ湫阅苓M(jìn)行準(zhǔn)確預(yù)測(cè).

    圖5 推力試驗(yàn)及試驗(yàn)結(jié)果Fig.5 Thrust test and test results

    3 空化原理和空化特性分析

    3.1 噴水推進(jìn)器空化原理

    在噴水推進(jìn)器實(shí)際運(yùn)行過(guò)程中,對(duì)應(yīng)某一工況,進(jìn)水流道提供有效的凈正吸頭(NPSHa),噴水推進(jìn)泵具有必須的凈正吸頭(NPSHr),可以根據(jù)進(jìn)水流道提供的有效凈正吸頭與泵需要的必須凈正吸頭的關(guān)系判斷空化是否發(fā)生.

    進(jìn)水流道提供的有效凈正吸頭可表示為:

    (1)

    式中:p∞為環(huán)境壓力,Pa;pv為液體的汽化壓力,Pa;vs為航速,m/s;ω為伴流系數(shù);ε為進(jìn)水流道損失系數(shù);ρ為液體的密度,kg/m3;g為重力加速度,m/s2;hj為噴泵相對(duì)于進(jìn)水流道進(jìn)口的安裝高度,m.

    噴水推進(jìn)泵需要的必須凈正吸頭可表示為:

    (2)

    式中:v1和ω1為葉輪進(jìn)口處的絕對(duì)速度和相對(duì)速度,m/s;λ為經(jīng)驗(yàn)系數(shù),與葉輪的幾何參數(shù)有關(guān),通常取值為0.15~0.40.

    當(dāng)進(jìn)水流道提供的有效凈正吸頭(NPSHa)大于噴水推進(jìn)泵需要的必須凈正吸頭(NPSHr)時(shí),泵內(nèi)不會(huì)發(fā)生空化;反之則會(huì)發(fā)生空化.

    定義噴水推進(jìn)器空化數(shù)σwj描述噴水推進(jìn)器空化情況[14]:

    (3)

    式中:pin為泵進(jìn)流面靜壓力,Pa;pv為液體的汽化壓力,Pa;vp為泵進(jìn)流面平均速度,m/s,ρ為液體的密度,kg/m3.

    3.2 噴水推進(jìn)器空化特性分析

    以固定航速30 kn、氣泡直徑0.2 mm工況為例,分析不同進(jìn)口含氣率對(duì)葉輪流道進(jìn)流面周向壓力及葉輪流道空化分布的影響.以固定航速30 kn、進(jìn)口含氣率15%、氣泡直徑0.2 mm工況為例,分析來(lái)流含氣對(duì)噴水推進(jìn)器空化非定常流動(dòng)的影響.

    3.2.1 來(lái)流含氣對(duì)噴水推進(jìn)器空化性能影響

    圖6(a)為泵進(jìn)流面壓力和速度隨進(jìn)口含氣率a變化關(guān)系,可以發(fā)現(xiàn)泵進(jìn)流面壓力pin隨進(jìn)口含氣率增加而逐漸增大,當(dāng)含氣率大于15%時(shí),再增加含氣率時(shí)壓力pin基本保持不變.泵進(jìn)流面速度vp隨進(jìn)口含氣率增加呈線性減小趨勢(shì),這是由于進(jìn)口含氣率增加,泵進(jìn)流面流量減小導(dǎo)致的.

    將圖6(a)不同進(jìn)口含氣率下的泵進(jìn)流面壓力pin和速度vp數(shù)據(jù)轉(zhuǎn)化為空化數(shù),得到噴水推進(jìn)器空化數(shù)隨進(jìn)口含氣率變化關(guān)系,如圖6(b).由圖可知,隨著進(jìn)口含氣率增加,空化數(shù)逐漸增大,此時(shí)噴水推進(jìn)器葉輪越不容易發(fā)生空化,這是因?yàn)殡S進(jìn)口含氣率增加,泵進(jìn)流面流速vp逐漸減小,泵進(jìn)流面壓力pin逐漸增大,由式(3)可知,噴水推進(jìn)器空化數(shù)此時(shí)會(huì)逐漸增大,說(shuō)明來(lái)流含氣在一定程度上有利于增強(qiáng)噴水推進(jìn)器抗空化性能.

    圖6 進(jìn)流面壓力、流速及空化數(shù)與進(jìn)口含氣率變化關(guān)系Fig.6 Relationship between inlet surface pressure,flow rate and cavitation number and inlet air content

    3.2.2 進(jìn)口含氣率對(duì)葉輪流道進(jìn)流面周向壓力的影響

    噴水推進(jìn)器在海面上工作時(shí)經(jīng)常會(huì)從進(jìn)水流道進(jìn)水口吸入空氣,噴水推進(jìn)器來(lái)流含氣會(huì)影響葉輪各流道進(jìn)流面壓力分布,從而進(jìn)一步影響噴水推進(jìn)器的空化性能.圖7為在固定航速30 kn工況下,不同進(jìn)口含氣率時(shí)各流道進(jìn)流面周向壓力分布.為了能更清楚的展示各流道壓力分布的狀況,將空間直角坐標(biāo)系轉(zhuǎn)為圓柱坐標(biāo)系θ,d為葉片輪轂到輪緣的無(wú)量綱距離.其中,流道1處于6°~78°,流道2處于78°~150°,以此類推.

    由圖可知,在純水空化工況下,各流道進(jìn)流面上均存在明顯的低壓區(qū),流道4和流道5區(qū)間范圍內(nèi)從輪轂至輪緣均顯現(xiàn)出大面積低壓區(qū),只在靠近葉片壓力面附近才存在較小范圍的高壓區(qū)域,流道4和流道5的低壓區(qū)面積明顯大于其他流道,說(shuō)明流道4和流道5內(nèi)的空化程度較其他流道更嚴(yán)重,其余流道低壓區(qū)主要分布在輪緣處.當(dāng)進(jìn)口含氣率為5%時(shí),各流道上的低壓區(qū)域由輪轂至輪緣逐漸減少,高壓區(qū)域逐漸增多,低壓區(qū)主要集中在葉片吸力面和輪緣附近,輪緣處的低壓區(qū)明顯多于輪轂處,說(shuō)明葉片輪緣處在泄漏流的擾動(dòng)下,其空化程度要比葉片其他位置嚴(yán)重.當(dāng)進(jìn)口含氣率達(dá)到25%時(shí),高壓區(qū)范圍變得越來(lái)越大,低壓區(qū)范圍變得越來(lái)越小.由分析可得,進(jìn)口含氣率對(duì)葉輪各流道進(jìn)流面的壓力分布影響十分明顯,在固定航速下,低壓區(qū)域隨著進(jìn)口含氣率的增加變得越來(lái)越小,同時(shí)在進(jìn)流面上周向壓力分布的差異性也越來(lái)越小,說(shuō)明來(lái)流含氣在一定程度上有助于改善噴水推進(jìn)器的空化性能.

    圖7 不同進(jìn)口含氣率下葉輪各流道進(jìn)流面周向壓力分布Fig.7 Circumferential pressure distribution on inlet surface of impeller with different inlet void fraction

    對(duì)葉輪各流道進(jìn)流面輪轂側(cè),中間位置以及輪緣側(cè)3個(gè)位置上的周向壓力p進(jìn)行定量分析,如圖8.

    圖8 葉輪各流道進(jìn)流面周向壓力分布Fig.8 Circumferential pressure distribution on inlet surfaceof impeller

    圖8(a)顯示,在葉輪5個(gè)流道上輪轂位置處,純水空化時(shí)壓力最低,進(jìn)口含氣率達(dá)到25%時(shí)壓力最高,輪轂側(cè)周向壓力與進(jìn)口含氣率成正相關(guān)性.圖8(b)顯示的進(jìn)流面中間流線處周向壓力與進(jìn)口含氣率的關(guān)系與輪轂位置保持一致.與進(jìn)流面輪轂處、進(jìn)流面中間位置上的壓力與進(jìn)口含氣率關(guān)系不同,進(jìn)流面輪緣處壓力與進(jìn)口含氣率的關(guān)系表現(xiàn)出不一樣的變化趨勢(shì),在各流道的前半個(gè)流道內(nèi),周向壓力與進(jìn)口含氣率成正相關(guān)性,進(jìn)口含氣率升高,周向壓力隨之升高;在后半個(gè)流道范圍內(nèi),周向壓力與進(jìn)口含氣率成負(fù)相關(guān)性,進(jìn)口含氣率升高,周向壓力反而降低.以上分析表明,輪緣處的流動(dòng)狀況較葉片其他位置更復(fù)雜.

    3.2.3 進(jìn)口含氣率對(duì)葉輪空泡分布的影響

    圖9為T=0.087 27 s時(shí)葉輪內(nèi)空泡分布隨進(jìn)口含氣率變化規(guī)律.如圖所示,空泡主要發(fā)生在葉片吸力面進(jìn)口邊附近,同時(shí)沿著輪緣方向向葉片出口邊方向發(fā)展.進(jìn)口含氣率對(duì)固定航速工況下空化狀況影響非常明顯,在固定航速工況下,隨著進(jìn)口含氣率增加,空化區(qū)域由輪轂向輪緣方向遞減,說(shuō)明在來(lái)流含氣條件下發(fā)生空化時(shí),氣體的加入有助于抑制空泡的發(fā)展.

    圖9 進(jìn)口含氣率對(duì)葉輪流道空泡分布的影響Fig.9 Effect of inlet void fraction on cavity distribution in flow channel

    由圖9分析可知,不同進(jìn)口含氣率下,空泡主要集中在葉片吸力面的進(jìn)口邊及輪緣附近.進(jìn)一步對(duì)進(jìn)口含氣率對(duì)葉片靠近輪緣位置空泡分布的影響進(jìn)行分析.圖10為此時(shí)刻下各葉片靠近輪緣位置的空間展開.由圖可知,在純水工況時(shí),各葉片在輪緣位置的空化現(xiàn)象較為明顯,空化區(qū)域主要顯現(xiàn)在吸力面進(jìn)口邊附近,同時(shí)各個(gè)葉片上的空化分布并不均勻.隨著進(jìn)口含氣率升高,葉片輪緣附近的空泡體積分?jǐn)?shù)數(shù)值和空化范圍逐漸減小,空泡分布不對(duì)稱性減弱,葉輪域內(nèi)的空化程度得到明顯改善.

    圖10 不同進(jìn)口含氣率下葉輪葉片間空泡分布(靠近輪緣處)Fig.10 Cavitation distribution between impeller blades (near the flange) under different inlet void fraction

    3.2.4 氣相介質(zhì)對(duì)空化非定常特性的影響

    由前文分析可知,噴水推進(jìn)器來(lái)流含氣時(shí),有助于抑制葉輪流道內(nèi)的空化發(fā)展,對(duì)噴水推進(jìn)器來(lái)流含氣空化流動(dòng)的非定常特性進(jìn)行分析,同時(shí)與純水空化流場(chǎng)非定常特性進(jìn)行對(duì)比.

    圖11為純水與來(lái)流含氣空化流場(chǎng)中,葉輪5個(gè)葉片上空泡面積變化情況.由圖可知,在葉輪1個(gè)旋轉(zhuǎn)周期內(nèi),5個(gè)葉片上的空泡面積數(shù)值均呈現(xiàn)無(wú)規(guī)律變化,同時(shí)顯現(xiàn)出明顯的波峰和波谷,說(shuō)明各葉片上的空化程度變化較為明顯;各葉片空泡面積波峰和波谷出現(xiàn)的位置存在明顯的相位差,表明各葉片上空化變化規(guī)律不一致,說(shuō)明葉片上的非對(duì)稱空化程度較高.在純水空化流場(chǎng)中,葉輪各葉片上的空泡面積幅值明顯大于來(lái)流含氣空化流場(chǎng),說(shuō)明純水空化流場(chǎng)的不穩(wěn)定性強(qiáng)于來(lái)流含氣空化流場(chǎng).

    圖11 葉輪葉片空泡面積隨時(shí)間變化Fig.11 Change of cavitation area of impeller blade with time

    圖12為純水與來(lái)流含氣空化流場(chǎng)中葉輪流道內(nèi)空泡形態(tài)隨時(shí)間變化規(guī)律,在葉輪旋轉(zhuǎn)過(guò)程中,各流道內(nèi)的空泡尺度呈非對(duì)稱狀態(tài),各流道的空化區(qū)域范圍不一致,純水空化流場(chǎng)中的空泡區(qū)域明顯大于來(lái)流含氣空化流場(chǎng),說(shuō)明來(lái)流含氣對(duì)空化的發(fā)展有一定抑制作用.不同時(shí)刻下,各流道內(nèi)的空泡體積變化過(guò)程中同時(shí)存在膨脹和縮小現(xiàn)象,說(shuō)明葉輪內(nèi)部的空化情況隨時(shí)間的變化具有不穩(wěn)定性,這種不穩(wěn)定性可能會(huì)對(duì)葉輪內(nèi)的流體造成強(qiáng)烈的沖擊.

    圖12 葉輪空泡演變Fig.12 Cavitation evolution of impeller

    由上圖分析可知,在葉輪旋轉(zhuǎn)過(guò)程中,不同時(shí)刻下葉輪輪緣區(qū)的空化情況較為嚴(yán)重復(fù)雜,由于葉輪輪緣間隙流動(dòng)對(duì)葉輪內(nèi)部流場(chǎng)具有較大的影響,需要對(duì)不同時(shí)刻下葉輪輪緣區(qū)的空化流動(dòng)情況做進(jìn)一步分析.圖13為不同時(shí)刻下,純水與來(lái)流含氣空化時(shí)葉輪輪緣區(qū)的空化分布情況.由圖可知,在純水空化中,5個(gè)葉片上的空化區(qū)域均附著在葉片吸力面進(jìn)口邊附近,但分布范圍、形態(tài)及空泡體積分?jǐn)?shù)數(shù)值各不相同,同時(shí)隨著葉輪旋轉(zhuǎn),各流道輪緣側(cè)的空泡變化規(guī)律不一致,以流道A和C為例,流道A輪緣側(cè)空化區(qū)域逐漸增大,說(shuō)明空化程度逐漸增強(qiáng),流道C輪緣側(cè)空化區(qū)域逐漸減小,說(shuō)明空化程度逐漸減弱.在來(lái)流含氣空化中,5個(gè)葉片上的空化分布范圍和空泡體積分?jǐn)?shù)數(shù)值均遠(yuǎn)低于同時(shí)刻時(shí)純水空化中對(duì)應(yīng)的空化分布范圍和空泡體積分?jǐn)?shù)數(shù)值,隨著葉輪旋轉(zhuǎn),各流道輪緣側(cè)的空泡變化規(guī)律各異,說(shuō)明在葉輪輪緣側(cè),空泡的成長(zhǎng)與潰滅現(xiàn)象并存,這也表現(xiàn)出輪緣側(cè)流動(dòng)的強(qiáng)不穩(wěn)定性.在純水與來(lái)流含氣空化流場(chǎng)中,葉輪輪緣側(cè)均存在復(fù)雜的空化演變過(guò)程,純水空化中的空化范圍更大,空化程度更深,說(shuō)明來(lái)流含氣中氣體的加入對(duì)葉輪輪緣側(cè)空化的發(fā)展起到了抑制作用.

    圖13 葉輪葉片間空泡演變(靠近輪緣處)Fig.13 Cavitation evolution between impeller blades (near the flange)

    在葉輪葉片吸力面設(shè)置5個(gè)壓力監(jiān)測(cè)點(diǎn)記錄空化流場(chǎng)的壓力脈動(dòng)特性.如圖14,葉片進(jìn)口邊到葉片出口邊各監(jiān)測(cè)點(diǎn)分別為BS1~BS5.

    圖14 吸力面監(jiān)測(cè)點(diǎn)示意圖Fig.14 Schematic diagram of suction surface monitoring points

    時(shí)域圖和頻域圖如圖15.由圖可知,來(lái)流含氣空化條件下葉片吸力面上監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)系數(shù)cp的波動(dòng)峰值和波動(dòng)頻率均小于純水空化流動(dòng),這可能與來(lái)流含氣空化流動(dòng)中空泡體積的減少有關(guān);從葉片進(jìn)口邊到出口邊,各監(jiān)測(cè)點(diǎn)波動(dòng)規(guī)律性不強(qiáng),同時(shí)在葉片吸力面進(jìn)口邊附近,由于空泡體積的不穩(wěn)定變化導(dǎo)致葉片出口邊壓力脈動(dòng)幅值小于葉片進(jìn)口邊.由頻域圖可知,各監(jiān)測(cè)點(diǎn)在純水與來(lái)流含氣空化流場(chǎng)中的主頻均為1倍轉(zhuǎn)頻f/fY,純水流場(chǎng)中的次頻成分多于來(lái)流含氣空化流場(chǎng).同時(shí)在同一監(jiān)測(cè)點(diǎn)處,來(lái)流含氣空化流場(chǎng)中的主頻幅值小于純水空化流場(chǎng),這說(shuō)明氣體的加入對(duì)空化流動(dòng)不穩(wěn)定性具有一定的減緩作用.

    圖15 純水空化與來(lái)流含氣空化條件下葉輪內(nèi)壓力脈動(dòng)特性Fig.15 Characteristics of pressure fluctuation in impeller under the condition of pure water cavitation and incoming air cavitation

    4 結(jié)論

    (1) 來(lái)流含氣空化時(shí),隨進(jìn)口含氣率升高,噴水推進(jìn)器空化數(shù)增大,說(shuō)明進(jìn)口含氣率增加在一定程度上有助于增強(qiáng)噴水推進(jìn)器抗空化性能;葉輪各流道進(jìn)流面低壓區(qū)面積減小,高壓區(qū)面積增大,周向壓力分布不均勻性減弱;葉輪內(nèi)空泡體積分?jǐn)?shù)減小,葉輪空化程度減弱.

    (2) 純水與來(lái)流含氣空化流場(chǎng)中,葉輪空泡體積分?jǐn)?shù)數(shù)值呈周期性規(guī)律變化,具有較強(qiáng)的非定常特性;各葉片空泡形態(tài)在不同時(shí)刻下呈現(xiàn)不同的變化狀態(tài),這種變化可能會(huì)對(duì)葉輪內(nèi)的流體造成強(qiáng)烈的沖擊;來(lái)流含氣空化流場(chǎng)中的壓力波動(dòng)峰值和振蕩頻率均小于純水空化流場(chǎng),說(shuō)明氣體的加入對(duì)空化流動(dòng)不穩(wěn)定具有一定的減緩作用.

    猜你喜歡
    含氣率含氣輪緣
    基于電容傳感器的液體管內(nèi)含氣率測(cè)量系統(tǒng)的研究*
    煤層氣井筒氣液兩相流數(shù)值模擬
    不同含氣率對(duì)采油單螺桿泵溫度和壓力的影響
    河南科技(2022年5期)2022-04-20 06:15:47
    淺談液態(tài)和固態(tài)輪緣潤(rùn)滑裝置的差異性
    地鐵車輛輪緣厚度偏磨問(wèn)題研究
    關(guān)于優(yōu)化四方平臺(tái)動(dòng)車組輪對(duì)踏面旋修的研究
    干式輪緣潤(rùn)滑器對(duì)地鐵車輛車輪保護(hù)效果的研究
    基于黏滯—彌散理論的含氣砂巖數(shù)值模擬與分析
    井震結(jié)合砂體識(shí)別與含氣評(píng)價(jià)技術(shù)在致密氣井挖潛中的應(yīng)用
    錄井工程(2017年4期)2017-03-16 06:10:34
    含氣率對(duì)AP1000核主泵影響的非定常分析
    国产精品,欧美在线| 变态另类成人亚洲欧美熟女 | 久久久久久免费高清国产稀缺| 亚洲一区二区三区不卡视频| 亚洲av片天天在线观看| 国产精华一区二区三区| 两性夫妻黄色片| 99国产精品一区二区蜜桃av| 日韩一卡2卡3卡4卡2021年| 国产亚洲精品久久久久5区| or卡值多少钱| 久热这里只有精品99| 亚洲男人天堂网一区| 欧美日韩瑟瑟在线播放| 亚洲成国产人片在线观看| 国产精品二区激情视频| 亚洲专区字幕在线| 美国免费a级毛片| 美女 人体艺术 gogo| 女警被强在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 视频在线观看一区二区三区| 国产欧美日韩综合在线一区二区| 少妇的丰满在线观看| 19禁男女啪啪无遮挡网站| 又紧又爽又黄一区二区| 夜夜看夜夜爽夜夜摸| 国产av一区在线观看免费| 日本三级黄在线观看| svipshipincom国产片| 午夜成年电影在线免费观看| 国产精品久久电影中文字幕| 法律面前人人平等表现在哪些方面| 欧美中文综合在线视频| 看免费av毛片| 咕卡用的链子| 手机成人av网站| 国产一区二区三区视频了| 色在线成人网| 国产精品久久久久久精品电影 | 精品乱码久久久久久99久播| 久久久久久久午夜电影| 亚洲av熟女| 波多野结衣一区麻豆| 在线av久久热| 亚洲天堂国产精品一区在线| 色综合婷婷激情| 国产精品99久久99久久久不卡| 成人av一区二区三区在线看| 国产又色又爽无遮挡免费看| 91成年电影在线观看| 夜夜爽天天搞| 免费观看人在逋| 亚洲精品久久成人aⅴ小说| 欧美日韩瑟瑟在线播放| 欧美精品亚洲一区二区| 97碰自拍视频| 18禁国产床啪视频网站| 丝袜美腿诱惑在线| 国产精华一区二区三区| 欧美一级毛片孕妇| 丁香六月欧美| 一级黄色大片毛片| 国产av在哪里看| 岛国视频午夜一区免费看| cao死你这个sao货| 国产精品 国内视频| 亚洲欧美激情综合另类| 国产片内射在线| 精品卡一卡二卡四卡免费| 国产亚洲av高清不卡| 波多野结衣高清无吗| 午夜a级毛片| 51午夜福利影视在线观看| 99精品久久久久人妻精品| 校园春色视频在线观看| 一级毛片精品| 黄色成人免费大全| 黑人巨大精品欧美一区二区蜜桃| 亚洲成人国产一区在线观看| 国产av精品麻豆| 国产欧美日韩一区二区精品| 色av中文字幕| 中文字幕av电影在线播放| 97人妻天天添夜夜摸| 亚洲av成人不卡在线观看播放网| 国产精品久久久久久精品电影 | 夜夜爽天天搞| 黄色a级毛片大全视频| 欧美亚洲日本最大视频资源| 久久久久久国产a免费观看| 18禁黄网站禁片午夜丰满| 亚洲 欧美一区二区三区| 成人国产一区最新在线观看| 亚洲九九香蕉| www.熟女人妻精品国产| 精品电影一区二区在线| 性色av乱码一区二区三区2| 中出人妻视频一区二区| 岛国在线观看网站| 亚洲成人国产一区在线观看| 在线观看免费午夜福利视频| 精品久久久精品久久久| 国产蜜桃级精品一区二区三区| 日韩大码丰满熟妇| 一区二区三区国产精品乱码| 午夜影院日韩av| 亚洲欧美一区二区三区黑人| 亚洲精品国产色婷婷电影| 精品少妇一区二区三区视频日本电影| 久久精品影院6| bbb黄色大片| 亚洲av成人一区二区三| 亚洲片人在线观看| 波多野结衣一区麻豆| 看黄色毛片网站| 一边摸一边做爽爽视频免费| 人人妻人人澡欧美一区二区 | 国产精品99久久99久久久不卡| 亚洲专区国产一区二区| 在线观看一区二区三区| 久久精品国产亚洲av高清一级| 给我免费播放毛片高清在线观看| 成人国语在线视频| 可以在线观看毛片的网站| 国产视频一区二区在线看| 日本欧美视频一区| 亚洲精品粉嫩美女一区| 久久久久久国产a免费观看| 国产av精品麻豆| 99国产精品一区二区蜜桃av| 青草久久国产| 精品日产1卡2卡| 一级毛片高清免费大全| 超碰成人久久| 国产成人免费无遮挡视频| svipshipincom国产片| 两人在一起打扑克的视频| 中亚洲国语对白在线视频| 亚洲av日韩精品久久久久久密| 亚洲人成电影观看| 国产av一区在线观看免费| 亚洲最大成人中文| 成人三级黄色视频| 成在线人永久免费视频| 免费在线观看亚洲国产| 9色porny在线观看| 成人国产综合亚洲| 巨乳人妻的诱惑在线观看| 亚洲无线在线观看| 大型av网站在线播放| 桃色一区二区三区在线观看| 成人精品一区二区免费| 欧美乱码精品一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 国产国语露脸激情在线看| 久久草成人影院| 亚洲精品国产色婷婷电影| 美国免费a级毛片| 国产精品98久久久久久宅男小说| 日韩欧美国产一区二区入口| 欧美老熟妇乱子伦牲交| 日韩欧美免费精品| 中文字幕人妻熟女乱码| xxx96com| 国产1区2区3区精品| netflix在线观看网站| 国产亚洲av高清不卡| 高潮久久久久久久久久久不卡| 极品人妻少妇av视频| tocl精华| 最近最新免费中文字幕在线| 嫩草影院精品99| 18禁黄网站禁片午夜丰满| 亚洲 欧美一区二区三区| 99久久99久久久精品蜜桃| 亚洲专区国产一区二区| 欧美黑人精品巨大| 中文字幕另类日韩欧美亚洲嫩草| aaaaa片日本免费| 国产精品电影一区二区三区| 国产亚洲av高清不卡| 少妇被粗大的猛进出69影院| 亚洲aⅴ乱码一区二区在线播放 | 1024香蕉在线观看| 一级,二级,三级黄色视频| 精品久久久久久成人av| 欧美 亚洲 国产 日韩一| 男女午夜视频在线观看| 淫秽高清视频在线观看| 国产国语露脸激情在线看| 欧美最黄视频在线播放免费| 亚洲精品在线美女| 日本精品一区二区三区蜜桃| tocl精华| 一级a爱视频在线免费观看| 亚洲熟妇熟女久久| 色综合欧美亚洲国产小说| 日韩中文字幕欧美一区二区| av在线天堂中文字幕| svipshipincom国产片| 亚洲专区中文字幕在线| 一边摸一边抽搐一进一小说| 欧美日本视频| 久久精品影院6| 精品欧美一区二区三区在线| 老司机靠b影院| а√天堂www在线а√下载| 国产在线精品亚洲第一网站| 亚洲美女黄片视频| 真人一进一出gif抽搐免费| 国产伦一二天堂av在线观看| 深夜精品福利| 狂野欧美激情性xxxx| 国产一级毛片七仙女欲春2 | 亚洲国产中文字幕在线视频| 亚洲 国产 在线| 国产精品 国内视频| 麻豆一二三区av精品| 男人舔女人下体高潮全视频| 黄色片一级片一级黄色片| 欧美精品亚洲一区二区| 午夜亚洲福利在线播放| 久久国产精品人妻蜜桃| 国产不卡一卡二| 午夜精品久久久久久毛片777| 不卡av一区二区三区| 日本精品一区二区三区蜜桃| 99久久综合精品五月天人人| 免费在线观看影片大全网站| 亚洲五月天丁香| 真人做人爱边吃奶动态| 欧美激情 高清一区二区三区| 国产亚洲欧美98| 夜夜看夜夜爽夜夜摸| 制服丝袜大香蕉在线| 久久热在线av| 亚洲国产精品999在线| 嫁个100分男人电影在线观看| 色播在线永久视频| 亚洲熟妇熟女久久| 欧美久久黑人一区二区| 午夜福利一区二区在线看| 久久这里只有精品19| 国产精品免费视频内射| 亚洲狠狠婷婷综合久久图片| 国产在线精品亚洲第一网站| 一a级毛片在线观看| 国产欧美日韩综合在线一区二区| 国产精品亚洲一级av第二区| 久久中文字幕人妻熟女| 麻豆一二三区av精品| 国产欧美日韩一区二区三| 欧美色欧美亚洲另类二区 | 老司机午夜十八禁免费视频| 国产一区二区三区视频了| 国产成人精品久久二区二区91| 国产1区2区3区精品| 色综合亚洲欧美另类图片| 国产97色在线日韩免费| 日韩免费av在线播放| 免费看十八禁软件| 欧美日韩瑟瑟在线播放| 俄罗斯特黄特色一大片| 午夜日韩欧美国产| АⅤ资源中文在线天堂| xxx96com| 男女下面进入的视频免费午夜 | 午夜成年电影在线免费观看| 无遮挡黄片免费观看| 黄片大片在线免费观看| 成人国产综合亚洲| 国产精品自产拍在线观看55亚洲| 黑人欧美特级aaaaaa片| 男人舔女人下体高潮全视频| 十分钟在线观看高清视频www| 黄色 视频免费看| 一级,二级,三级黄色视频| 久久精品亚洲精品国产色婷小说| 亚洲中文av在线| 国产av一区在线观看免费| 最好的美女福利视频网| 黄色片一级片一级黄色片| 亚洲国产高清在线一区二区三 | 麻豆国产av国片精品| 亚洲色图av天堂| 黄色视频,在线免费观看| 久久国产精品人妻蜜桃| 亚洲视频免费观看视频| 国产av一区在线观看免费| 在线视频色国产色| 亚洲一区高清亚洲精品| 日韩欧美国产一区二区入口| 制服诱惑二区| 999久久久精品免费观看国产| 一二三四社区在线视频社区8| 免费在线观看亚洲国产| 韩国精品一区二区三区| 人成视频在线观看免费观看| 精品一区二区三区av网在线观看| 1024香蕉在线观看| 高清在线国产一区| 日本 av在线| 狂野欧美激情性xxxx| 黑人巨大精品欧美一区二区蜜桃| 在线观看日韩欧美| 国产99白浆流出| 视频在线观看一区二区三区| 国产乱人伦免费视频| 91av网站免费观看| 亚洲国产日韩欧美精品在线观看 | 国产成人一区二区三区免费视频网站| 久久久久久人人人人人| 真人一进一出gif抽搐免费| 又紧又爽又黄一区二区| 麻豆av在线久日| 757午夜福利合集在线观看| 亚洲性夜色夜夜综合| 丁香六月欧美| 久久久久久国产a免费观看| 国产成人av激情在线播放| av视频免费观看在线观看| 亚洲欧美精品综合久久99| 婷婷精品国产亚洲av在线| 久久久国产成人免费| 亚洲成人免费电影在线观看| a级毛片在线看网站| 不卡av一区二区三区| 黄色视频,在线免费观看| 国产成人免费无遮挡视频| 免费高清视频大片| 91大片在线观看| 精品一品国产午夜福利视频| 免费久久久久久久精品成人欧美视频| 亚洲欧美日韩无卡精品| 欧美丝袜亚洲另类 | 十分钟在线观看高清视频www| 亚洲全国av大片| 99在线视频只有这里精品首页| 中文亚洲av片在线观看爽| 免费搜索国产男女视频| 777久久人妻少妇嫩草av网站| 满18在线观看网站| 色老头精品视频在线观看| 精品国产乱码久久久久久男人| 成人精品一区二区免费| 免费高清在线观看日韩| 香蕉久久夜色| 国产伦人伦偷精品视频| 老司机深夜福利视频在线观看| av在线天堂中文字幕| 色播在线永久视频| 久久中文字幕人妻熟女| 搞女人的毛片| 亚洲一码二码三码区别大吗| 欧美日韩瑟瑟在线播放| 色播亚洲综合网| 不卡av一区二区三区| 制服丝袜大香蕉在线| 国产亚洲精品一区二区www| 久久久精品欧美日韩精品| 一级毛片高清免费大全| 天天添夜夜摸| 亚洲av成人一区二区三| 欧美激情高清一区二区三区| 亚洲专区字幕在线| 欧美成人性av电影在线观看| 成人国产综合亚洲| 亚洲五月天丁香| 国产精品一区二区精品视频观看| 久久精品91无色码中文字幕| 亚洲人成网站在线播放欧美日韩| 久久性视频一级片| 不卡一级毛片| 可以在线观看的亚洲视频| 午夜福利高清视频| x7x7x7水蜜桃| 色在线成人网| 国产成人欧美在线观看| 日韩 欧美 亚洲 中文字幕| 欧美中文日本在线观看视频| 国产一区在线观看成人免费| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利欧美成人| 午夜福利成人在线免费观看| 婷婷六月久久综合丁香| 级片在线观看| 人人妻人人爽人人添夜夜欢视频| 精品国产乱子伦一区二区三区| 人人妻人人爽人人添夜夜欢视频| 欧美最黄视频在线播放免费| 90打野战视频偷拍视频| 一本综合久久免费| 69av精品久久久久久| 国产精品国产高清国产av| 免费搜索国产男女视频| 一级黄色大片毛片| 后天国语完整版免费观看| 每晚都被弄得嗷嗷叫到高潮| 精品一区二区三区四区五区乱码| 久久人人精品亚洲av| 少妇的丰满在线观看| 美女免费视频网站| √禁漫天堂资源中文www| 一区福利在线观看| 高清在线国产一区| 亚洲熟妇中文字幕五十中出| 国产精品野战在线观看| 欧美日韩亚洲综合一区二区三区_| 女人爽到高潮嗷嗷叫在线视频| 欧美激情高清一区二区三区| 老司机在亚洲福利影院| a在线观看视频网站| 国产午夜福利久久久久久| 一区在线观看完整版| 亚洲专区中文字幕在线| 国产精品亚洲av一区麻豆| 一级a爱片免费观看的视频| www国产在线视频色| av在线播放免费不卡| 欧美国产日韩亚洲一区| 亚洲成av片中文字幕在线观看| 看免费av毛片| 俄罗斯特黄特色一大片| 性色av乱码一区二区三区2| 香蕉久久夜色| 色播亚洲综合网| 国产视频一区二区在线看| 制服诱惑二区| 亚洲伊人色综图| 国产精品亚洲av一区麻豆| 国产真人三级小视频在线观看| 亚洲国产精品久久男人天堂| 男女床上黄色一级片免费看| 丝袜美足系列| 一夜夜www| 岛国在线观看网站| a在线观看视频网站| 黄色 视频免费看| 国产成人一区二区三区免费视频网站| 国产成人系列免费观看| 久久久久精品国产欧美久久久| 99热只有精品国产| 我的亚洲天堂| or卡值多少钱| 激情在线观看视频在线高清| 禁无遮挡网站| 此物有八面人人有两片| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看免费视频网站a站| 国产片内射在线| 久久草成人影院| 亚洲精品美女久久av网站| 熟女少妇亚洲综合色aaa.| 日韩欧美国产一区二区入口| 久久香蕉激情| 国产99白浆流出| 十八禁人妻一区二区| 女人爽到高潮嗷嗷叫在线视频| 久久久久久久午夜电影| 怎么达到女性高潮| 亚洲av美国av| 在线播放国产精品三级| av欧美777| 欧美另类亚洲清纯唯美| 成人三级黄色视频| 女人精品久久久久毛片| 伊人久久大香线蕉亚洲五| 如日韩欧美国产精品一区二区三区| 一边摸一边做爽爽视频免费| 一级,二级,三级黄色视频| 国产乱人伦免费视频| 高清黄色对白视频在线免费看| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩中文字幕国产精品一区二区三区 | 多毛熟女@视频| www.999成人在线观看| 精品国产国语对白av| 最近最新中文字幕大全免费视频| 999久久久精品免费观看国产| 国内精品久久久久久久电影| 婷婷六月久久综合丁香| 日日摸夜夜添夜夜添小说| 久久人妻熟女aⅴ| 黄色成人免费大全| 日本一区二区免费在线视频| 国产精品香港三级国产av潘金莲| 日韩大尺度精品在线看网址 | 黑人巨大精品欧美一区二区mp4| 99香蕉大伊视频| 亚洲国产毛片av蜜桃av| 黑人欧美特级aaaaaa片| 国内毛片毛片毛片毛片毛片| 国产亚洲精品久久久久久毛片| 99精品久久久久人妻精品| 欧美av亚洲av综合av国产av| 国产精品一区二区免费欧美| 叶爱在线成人免费视频播放| 国产精品亚洲一级av第二区| 免费在线观看亚洲国产| 在线永久观看黄色视频| 免费一级毛片在线播放高清视频 | 精品国产超薄肉色丝袜足j| 免费在线观看亚洲国产| 国产精品久久久久久精品电影 | 亚洲少妇的诱惑av| 日本三级黄在线观看| 19禁男女啪啪无遮挡网站| 在线十欧美十亚洲十日本专区| 国产亚洲av嫩草精品影院| 国产av在哪里看| 高清在线国产一区| 露出奶头的视频| 亚洲精品一卡2卡三卡4卡5卡| 操出白浆在线播放| 日本 av在线| 一级毛片精品| 免费搜索国产男女视频| 亚洲专区中文字幕在线| 欧洲精品卡2卡3卡4卡5卡区| 国产三级黄色录像| 精品福利观看| 国产极品粉嫩免费观看在线| 国产高清激情床上av| 亚洲精品av麻豆狂野| 大型av网站在线播放| 每晚都被弄得嗷嗷叫到高潮| 老司机福利观看| 中文字幕人妻丝袜一区二区| 超碰成人久久| 国产熟女午夜一区二区三区| 99re在线观看精品视频| 亚洲精品av麻豆狂野| 日韩av在线大香蕉| 亚洲 国产 在线| 欧美国产日韩亚洲一区| 亚洲人成77777在线视频| 亚洲欧美精品综合一区二区三区| 非洲黑人性xxxx精品又粗又长| 国产蜜桃级精品一区二区三区| 热re99久久国产66热| 一个人观看的视频www高清免费观看 | АⅤ资源中文在线天堂| 欧美日韩瑟瑟在线播放| 久99久视频精品免费| 亚洲精品久久成人aⅴ小说| 真人做人爱边吃奶动态| 欧美成人一区二区免费高清观看 | 中文字幕最新亚洲高清| 免费看美女性在线毛片视频| 国产99白浆流出| 91麻豆av在线| 日韩三级视频一区二区三区| 午夜激情av网站| 99久久国产精品久久久| 午夜免费观看网址| 国产成+人综合+亚洲专区| 午夜精品在线福利| 啦啦啦免费观看视频1| 可以在线观看毛片的网站| 老鸭窝网址在线观看| 亚洲精品在线观看二区| 9色porny在线观看| 亚洲精品中文字幕一二三四区| 国产高清视频在线播放一区| 亚洲国产精品久久男人天堂| 嫩草影院精品99| 激情视频va一区二区三区| 中出人妻视频一区二区| 女生性感内裤真人,穿戴方法视频| 777久久人妻少妇嫩草av网站| 久久这里只有精品19| 久久中文看片网| 99国产精品一区二区三区| 久9热在线精品视频| 精品一品国产午夜福利视频| 亚洲午夜理论影院| 欧美成人一区二区免费高清观看 | 九色国产91popny在线| 性色av乱码一区二区三区2| 成人av一区二区三区在线看| 欧美国产日韩亚洲一区| 欧美日韩亚洲综合一区二区三区_| 亚洲一区二区三区色噜噜| 欧美日韩亚洲国产一区二区在线观看| 久久久久久久午夜电影| 亚洲三区欧美一区| www国产在线视频色| 国产午夜精品久久久久久| 麻豆成人av在线观看| 桃色一区二区三区在线观看| 长腿黑丝高跟| 一级a爱视频在线免费观看| 久久久久久亚洲精品国产蜜桃av| 男女床上黄色一级片免费看| 老熟妇乱子伦视频在线观看| 午夜老司机福利片| 在线免费观看的www视频| 日韩精品中文字幕看吧| 国内精品久久久久久久电影| 男女床上黄色一级片免费看| 老司机在亚洲福利影院| 两人在一起打扑克的视频| 久久精品影院6| 一二三四在线观看免费中文在| 首页视频小说图片口味搜索| 老汉色av国产亚洲站长工具| 欧美精品亚洲一区二区| 亚洲一区高清亚洲精品| 欧美激情极品国产一区二区三区| 一边摸一边做爽爽视频免费| 成人亚洲精品av一区二区| 国产精品自产拍在线观看55亚洲| 丝袜在线中文字幕| 久久久久精品国产欧美久久久|