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

    基于CFD的潛液式液化天然氣泵導(dǎo)葉設(shè)計(jì)

    2016-12-21 03:31:38錢(qián)濤陳紅勛梁成鵬
    關(guān)鍵詞:喉部導(dǎo)葉轉(zhuǎn)角

    錢(qián)濤,陳紅勛,梁成鵬

    (上海大學(xué)上海市應(yīng)用數(shù)學(xué)和力學(xué)研究所,上?!?00072)

    基于CFD的潛液式液化天然氣泵導(dǎo)葉設(shè)計(jì)

    錢(qián)濤,陳紅勛,梁成鵬

    (上海大學(xué)上海市應(yīng)用數(shù)學(xué)和力學(xué)研究所,上海200072)

    潛液式液化天然氣(liquefied natural gas,LNG)泵工作時(shí),屏蔽電機(jī)和泵體全部浸沒(méi)在低溫液體中.為減小泵的徑向和軸向尺寸,潛液式LNG泵采用了一種特殊結(jié)構(gòu)的導(dǎo)葉.在分析潛液式LNG泵導(dǎo)葉結(jié)構(gòu)特點(diǎn)的基礎(chǔ)上,研究導(dǎo)葉進(jìn)口喉部寬度和折轉(zhuǎn)角對(duì)泵設(shè)計(jì)工況水力性能的影響.首先,設(shè)計(jì)不同幾何參數(shù)的導(dǎo)葉,并分別與同一葉輪進(jìn)行匹配;再通過(guò)ANSYS CFX軟件,采用標(biāo)準(zhǔn)k-ε湍流模型對(duì)各導(dǎo)葉分別進(jìn)行全流場(chǎng)數(shù)值計(jì)算.計(jì)算結(jié)果表明:進(jìn)口喉部寬度是潛液式LNG泵導(dǎo)葉的關(guān)鍵尺寸,設(shè)計(jì)時(shí)需重點(diǎn)考慮;進(jìn)口喉部寬度存在最優(yōu)值,且最優(yōu)值大于經(jīng)驗(yàn)值;進(jìn)口折轉(zhuǎn)角對(duì)泵揚(yáng)程和效率影響較小.因此,設(shè)計(jì)導(dǎo)葉時(shí)可優(yōu)先確定其他關(guān)鍵尺寸,再通過(guò)調(diào)節(jié)進(jìn)口折轉(zhuǎn)角改善導(dǎo)葉的結(jié)構(gòu).

    潛液式液化天然氣泵;導(dǎo)葉;喉部寬度;折轉(zhuǎn)角;數(shù)值模擬

    近年來(lái),清潔能源液化天然氣(liquefied natural gas,LNG)的使用在我國(guó)發(fā)展迅速,潛液式LNG泵(簡(jiǎn)稱LNG泵)的應(yīng)用范圍也在不斷擴(kuò)大.LNG泵的結(jié)構(gòu)如圖1所示,其結(jié)構(gòu)設(shè)計(jì)的特點(diǎn)有助于屏蔽電機(jī)并使泵體全部浸沒(méi)在低溫液體中,達(dá)到零泄漏[1].圖2為L(zhǎng)NG泵的過(guò)流部件,其基本工作過(guò)程如下:與葉輪和誘導(dǎo)輪共軸的電機(jī)帶動(dòng)泵軸高速旋轉(zhuǎn),從而帶動(dòng)誘導(dǎo)輪和葉輪轉(zhuǎn)動(dòng),導(dǎo)葉則固定在泵體上;低溫液體經(jīng)過(guò)誘導(dǎo)輪進(jìn)入葉輪后在離心力作用下從葉輪出口甩出,然后進(jìn)入導(dǎo)葉入口,減速擴(kuò)壓后進(jìn)入到下一級(jí)的葉輪,再經(jīng)過(guò)下一級(jí)導(dǎo)葉壓出.

    圖1 LNG泵的結(jié)構(gòu)Fig.1 Structure of LNG pump

    圖2 LNG泵的過(guò)流部件Fig.2 Overflow parts of LNG pump

    導(dǎo)葉是LNG泵的關(guān)鍵過(guò)流部件之一,其設(shè)計(jì)的好壞對(duì)泵的性能有著重要的影響.已有研究表明,導(dǎo)葉內(nèi)的水力損失約占泵內(nèi)水力損失的40%~50%[2-6].導(dǎo)葉的主要作用是收集葉輪中甩出的液體,并將液體的動(dòng)能轉(zhuǎn)化為壓能后輸送至下一級(jí)吸入口或排出口.根據(jù)結(jié)構(gòu)特點(diǎn)的不同,現(xiàn)有的導(dǎo)葉主要有徑向?qū)~、流道式導(dǎo)葉和空間導(dǎo)葉(見(jiàn)圖3).徑向?qū)~的正導(dǎo)葉和反導(dǎo)葉相互隔開(kāi),導(dǎo)葉出口的直徑相比葉輪直徑增大30%左右[7].但只有反導(dǎo)葉而沒(méi)有正導(dǎo)葉的徑向?qū)~雖然可以縮小徑向尺寸,但匹配高轉(zhuǎn)速葉輪時(shí),由于甩出的液體圓周分速度較大,液體會(huì)直接進(jìn)入反導(dǎo)葉從而產(chǎn)生堵塞導(dǎo)致水力性能降低.流道式導(dǎo)葉的工作效率較高,但設(shè)計(jì)和制造難度較大[8].空間導(dǎo)葉和流道式導(dǎo)葉類似,正反導(dǎo)葉合為一體,其軸向尺寸較大[9].由于LNG泵轉(zhuǎn)速較高,且其泵體浸沒(méi)在液體容器中,緊湊的結(jié)構(gòu)設(shè)計(jì)特點(diǎn)對(duì)節(jié)約泵裝置安裝空間和生產(chǎn)成本意義重大,因此LNG泵的導(dǎo)葉采用如圖2所示的特殊設(shè)計(jì).考慮到導(dǎo)葉結(jié)構(gòu)的特殊性,有必要分析其設(shè)計(jì)參數(shù)對(duì)泵性能的影響,這不僅可為L(zhǎng)NG泵導(dǎo)葉的設(shè)計(jì)優(yōu)化提供參考,也可為多級(jí)離心泵導(dǎo)葉的設(shè)計(jì)提供新思路.

    圖3 泵導(dǎo)葉的常見(jiàn)設(shè)計(jì)Fig.3 Popular design of pump diffuser

    隨著流體力學(xué)的發(fā)展和計(jì)算機(jī)性能的提高,計(jì)算流體動(dòng)力學(xué)(computational fluid dynamics,CFD)技術(shù)成為模擬泵內(nèi)部流動(dòng)規(guī)律和預(yù)測(cè)水力性能的重要工具.相對(duì)于離心泵的葉輪及其蝸殼內(nèi)部流場(chǎng)的分析計(jì)算和特性研究,國(guó)內(nèi)對(duì)于導(dǎo)葉性能的數(shù)值研究較少.徐媛暉等[10]設(shè)計(jì)出一款新型軸向?qū)~,并將同一葉輪分別與兩個(gè)導(dǎo)葉進(jìn)行匹配,再對(duì)這兩組模型分別進(jìn)行了全流場(chǎng)數(shù)值計(jì)算,結(jié)果表明該新型軸向?qū)~能夠改善導(dǎo)葉內(nèi)部流動(dòng),提高井泵的水力性能.周嶺等[11]針對(duì)深井離心泵,設(shè)計(jì)了三維曲面導(dǎo)葉與圓柱形導(dǎo)葉,并進(jìn)行了全流場(chǎng)數(shù)值模擬,為優(yōu)化導(dǎo)葉水力性能提供了參考.叢小青等[12]對(duì)一空間導(dǎo)葉式離心泵進(jìn)行全流道三維湍流數(shù)值計(jì)算,分析了設(shè)計(jì)工況下整個(gè)流道、環(huán)形空間及其空間導(dǎo)葉內(nèi)部的流場(chǎng)分布,并通過(guò)試驗(yàn)驗(yàn)證了計(jì)算結(jié)果.李家文等[13]針對(duì)液氫泵提出了一種流道式導(dǎo)葉的設(shè)計(jì)方法,降低了級(jí)間導(dǎo)葉的損耗,提高了液氫泵的工作效率.

    本工作主要借助數(shù)值模擬分析不同幾何參數(shù)對(duì)LNG泵導(dǎo)葉水力性能的影響.首先,分析LNG泵導(dǎo)葉的結(jié)構(gòu)特點(diǎn)和幾何參數(shù)選取方法;然后,在保證其他參數(shù)不變的情況下,分別改變進(jìn)口喉部寬度和轉(zhuǎn)折角,設(shè)計(jì)出不同的導(dǎo)葉,并將導(dǎo)葉與同一葉輪進(jìn)行匹配.通過(guò)ANSYS CFX軟件,在設(shè)計(jì)工況下對(duì)各模型進(jìn)行全流場(chǎng)定常數(shù)值計(jì)算;最后,根據(jù)三維數(shù)值計(jì)算結(jié)果分析進(jìn)口喉部寬度和折轉(zhuǎn)角對(duì)LNG泵導(dǎo)葉水力性能的影響.

    1 LNG泵導(dǎo)葉的幾何描述

    圖4為L(zhǎng)NG泵導(dǎo)葉的三維模型,其中R3為導(dǎo)葉進(jìn)口半徑,R4為導(dǎo)葉外圓半徑,α3和β為導(dǎo)葉葉片進(jìn)口的安裝角和折轉(zhuǎn)角,a3和b3為導(dǎo)葉進(jìn)口喉部的寬度和高度,δ為導(dǎo)葉葉片入口的厚度,a4和b4為導(dǎo)葉出口喉部的寬度和高度.由圖4可以看出,流道沿周向均勻布置,圓弧AB,BC沿著導(dǎo)葉中心軸旋轉(zhuǎn)相應(yīng)的角度得到下一個(gè)流道的進(jìn)口邊.

    圖4 LNG泵導(dǎo)葉的三維模型Fig.4 3D model of LNG pump diffuser

    LNG泵導(dǎo)葉的主要結(jié)構(gòu)為進(jìn)口擴(kuò)散段、轉(zhuǎn)向段和出口導(dǎo)流柵三部分.液體以一定的液流角沿AB段匯集到導(dǎo)葉入口,再流入導(dǎo)葉進(jìn)口擴(kuò)散段,由于β>α3,流道寬度和深度都逐漸增加,液體的速度降低,動(dòng)能轉(zhuǎn)化為壓能;接著,液體經(jīng)過(guò)導(dǎo)葉圓弧轉(zhuǎn)向段,沿圓周半徑方向流出導(dǎo)葉擴(kuò)散段;最后,不同流道的液體在空腔匯集后經(jīng)過(guò)導(dǎo)流柵沿軸向流入下一級(jí)葉輪.液體沿軸向進(jìn)入導(dǎo)葉入口時(shí),流動(dòng)類似于空間導(dǎo)葉,而進(jìn)口擴(kuò)散段又類似于徑向?qū)~的正導(dǎo)葉,因此LNG泵導(dǎo)葉具有軸向和徑向尺寸均較小的優(yōu)點(diǎn).另外,流道多為圓弧形結(jié)構(gòu),設(shè)計(jì)和加工都較為方便.

    2 水力模型

    LNG泵設(shè)計(jì)工況參數(shù)如下:流量Q=20.4 m3/h,揚(yáng)程H=200 m,轉(zhuǎn)速n=6 000 r/min,級(jí)數(shù)為兩級(jí).一級(jí)葉輪和導(dǎo)葉與二級(jí)葉輪和導(dǎo)葉的水力設(shè)計(jì)相同,因此本工作只取LNG泵的一級(jí)葉輪和導(dǎo)葉進(jìn)行數(shù)值模擬.考慮到誘導(dǎo)輪的主要作用是改變?nèi)~輪的進(jìn)口入流條件以及提高葉輪的抗空化性能,本工作計(jì)算時(shí)未考慮誘導(dǎo)輪的作用.流體域由進(jìn)口段、一級(jí)葉輪、一級(jí)導(dǎo)葉、腔體以及出口段組成.

    2.1葉輪水力模型

    LNG泵一級(jí)葉輪的幾何參數(shù)如表1所示,其三維造型如圖5所示.

    表1 LNG泵一級(jí)葉輪的幾何參數(shù)Table 1 Geometric parameters of the primary impeller in LNG pump

    圖5 LNG泵葉輪的三維模型Fig.5 3D model of LNG pump impeller

    2.2導(dǎo)葉水力模型

    借鑒徑向?qū)~的設(shè)計(jì)方法,LNG泵導(dǎo)葉主要幾何參數(shù)的計(jì)算如下.

    (1)根據(jù)葉輪出口直徑D2,選取導(dǎo)葉進(jìn)口直徑D3,有

    (2)選取導(dǎo)葉葉片數(shù)Z2.本工作中葉輪的葉片數(shù)為4,為減小泵的振動(dòng),導(dǎo)葉的葉片數(shù)選為7.(3)根據(jù)速度系數(shù)法[2],計(jì)算進(jìn)口喉部面積F3,有

    (4)確定導(dǎo)葉出口喉部面積F4,有

    式中,出口速度v4選取3~5 m/s的流速,目的是減少液體在環(huán)形空腔內(nèi)的流動(dòng)損失.

    (5)計(jì)算導(dǎo)葉進(jìn)出口的喉部寬度a3,a4以及高度b3,b4.因?yàn)?/p>

    且有數(shù)據(jù)表明喉部斷面接近正方形效果較佳[2],所以可取

    另外,導(dǎo)葉的外圓直徑為

    (6)選取導(dǎo)葉葉片的入口厚度δ,有

    (7)計(jì)算導(dǎo)葉進(jìn)口安裝角α3,有

    式中,vm3為導(dǎo)葉進(jìn)口軸面速度,vu3為導(dǎo)葉進(jìn)口圓周速度,α為進(jìn)口沖角.

    (8)根據(jù)進(jìn)口安裝角α3,選取折轉(zhuǎn)角β=α3+(6°~10°),其中β的選取應(yīng)考慮導(dǎo)葉的結(jié)構(gòu)特點(diǎn),避免導(dǎo)葉的部分面穿透或太薄.

    根據(jù)上述設(shè)計(jì)方法初步得到的LNG泵導(dǎo)葉A0的幾何參數(shù)如表2所示.可以看出,導(dǎo)葉外圓直徑僅比葉輪直徑增加12.8%.LNG泵導(dǎo)葉的進(jìn)口擴(kuò)散段的幾何參數(shù)有α3,β,a3,它們決定了導(dǎo)葉進(jìn)口喉部面積的大小,也決定了擴(kuò)散段的形狀,從而影響導(dǎo)葉與葉輪的匹配程度,其中α3可通過(guò)計(jì)算液流進(jìn)口角確定,進(jìn)口喉部寬度a3和折轉(zhuǎn)角β的確定則有待進(jìn)一步研究.

    表2 LNG泵導(dǎo)葉的幾何參數(shù)Table 2 Geometric parameters of LNG pump diffuser

    3 數(shù)值模擬

    為分析LNG泵導(dǎo)葉在設(shè)計(jì)工況下的水力性能,對(duì)LNG泵一級(jí)葉輪和一級(jí)導(dǎo)葉匹配后的流道進(jìn)行定常數(shù)值計(jì)算.

    3.1流道建模和網(wǎng)格劃分

    計(jì)算流體域分塊建模如下:通過(guò)UG軟件的求差運(yùn)算功能得到葉輪旋轉(zhuǎn)部分流道和導(dǎo)葉靜止部分流道,其中進(jìn)口和出口段適當(dāng)延長(zhǎng),考慮葉輪蓋板與腔體及導(dǎo)葉之間的間隙,而不考慮液體軸間的泄露.LNG泵的幾何模型如圖6所示.LNG泵流道的網(wǎng)格劃分在ICEM CFD軟件中完成(見(jiàn)圖7),其中進(jìn)口段、出口段以及葉輪采用全六面體結(jié)構(gòu)網(wǎng)格,導(dǎo)葉部分由于幾何結(jié)構(gòu)過(guò)于復(fù)雜,采用適應(yīng)性較好的非結(jié)構(gòu)四面體網(wǎng)格.計(jì)算時(shí)通過(guò)網(wǎng)格拼接技術(shù)將各塊網(wǎng)格耦合,其中進(jìn)口為9萬(wàn)計(jì)算單元,葉輪為152萬(wàn)計(jì)算單元,出口為13萬(wàn)計(jì)算單元.對(duì)于不同的導(dǎo)葉,采用同一套進(jìn)口、葉輪和出口網(wǎng)格.

    圖6 LNG泵的幾何模型Fig.6 Geometric model of LNG pump

    圖7 LNG泵的計(jì)算網(wǎng)格Fig.7 Computational mesh of LNG pump

    導(dǎo)葉進(jìn)行網(wǎng)格劃分時(shí),導(dǎo)葉考慮擴(kuò)散段為關(guān)鍵過(guò)流段,因此對(duì)其流道采用密度盒加密設(shè)置.為檢驗(yàn)不同網(wǎng)格數(shù)對(duì)計(jì)算結(jié)果的影響,根據(jù)A0導(dǎo)葉與葉輪匹配計(jì)算的揚(yáng)程來(lái)進(jìn)行網(wǎng)格無(wú)關(guān)性檢驗(yàn),其中導(dǎo)葉擴(kuò)散段密度盒尺寸設(shè)置為1.5~0.4 mm,導(dǎo)葉網(wǎng)格數(shù)從200萬(wàn)增加到540萬(wàn).數(shù)值計(jì)算結(jié)果如圖8所示,可以看出:當(dāng)導(dǎo)葉網(wǎng)格數(shù)從200萬(wàn)增加到540萬(wàn)時(shí),揚(yáng)程變化不超過(guò)4%;當(dāng)網(wǎng)格數(shù)大于430萬(wàn)時(shí),揚(yáng)程變化在0.5%以內(nèi).考慮到計(jì)算的效率和精度,本工作中導(dǎo)葉的網(wǎng)格設(shè)置如下:擴(kuò)散段流道密度盒尺寸為0.5 mm,全局尺寸為2 mm,導(dǎo)葉單元數(shù)為400萬(wàn)~500萬(wàn).對(duì)于不同的導(dǎo)葉,均采用相同的網(wǎng)格設(shè)置方式.

    圖8 揚(yáng)程隨導(dǎo)葉網(wǎng)格數(shù)的變化曲線Fig.8 Variation curve of head with diffuser mesh number

    3.2計(jì)算模型與定解條件

    對(duì)于離心泵內(nèi)部的不可壓縮流動(dòng),本工作采用ANSYS CFX求解守恒型雷諾時(shí)均方程,并選擇標(biāo)準(zhǔn)k-ε兩方程湍流模型進(jìn)行封閉.葉輪采用固系于旋轉(zhuǎn)葉輪上的相對(duì)參考系,導(dǎo)葉及腔體采用靜止坐標(biāo)系,動(dòng)、靜域之間采用Frozen Rotor的動(dòng)靜交接面模型.計(jì)算時(shí),對(duì)流項(xiàng)采用高分辨率格式,流場(chǎng)求解使用全隱式多重網(wǎng)格耦合方法,這不僅能夠同時(shí)求解連續(xù)方程和動(dòng)量方程,還能提高計(jì)算收斂速度和穩(wěn)定性.

    進(jìn)口條件選用總壓進(jìn)口條件,參考?jí)毫υO(shè)為101.325 kPa(1標(biāo)準(zhǔn)大氣壓),總壓設(shè)為0 Pa;出口設(shè)定流量出口;葉輪旋轉(zhuǎn)部分設(shè)定轉(zhuǎn)速6 000 r/min;時(shí)間步設(shè)為1 500步.自定義設(shè)置介質(zhì)屬性,液態(tài)天然氣密度為441.36 kg/m3,動(dòng)力黏度為0.2 mPa·s.

    3.3進(jìn)口喉部寬度的影響

    為分析進(jìn)口喉部寬度對(duì)導(dǎo)葉性能的影響,設(shè)計(jì)了不同進(jìn)口喉部寬度的導(dǎo)葉A1~A7(見(jiàn)表3).圖9為A1,A4,A7的結(jié)構(gòu)對(duì)比.可以看出,隨著進(jìn)口喉部寬度a3的增大,導(dǎo)葉擴(kuò)散段結(jié)構(gòu)由狹長(zhǎng)形變成短寬形.通過(guò)數(shù)值模擬,得到了a3不同時(shí)的單級(jí)泵在設(shè)計(jì)工況點(diǎn)下的揚(yáng)程和效率(見(jiàn)圖10).

    表3 不同LNG泵導(dǎo)葉的進(jìn)口喉部寬度Table 3 Inlet widths of different LNG pump diffusersmm

    圖9 導(dǎo)葉A1,A4和A7的對(duì)比Fig.9 Comparisons of diffusers A1,A4 and A7

    從圖10可以看出:導(dǎo)葉進(jìn)口喉部寬度對(duì)單級(jí)泵的揚(yáng)程和效率有較大的影響,隨著進(jìn)口寬度的增大,單級(jí)泵的揚(yáng)程和效率先有較大的提升,然后提升幅度逐漸收窄,達(dá)到最大值后回落;進(jìn)口喉部寬度存在最優(yōu)值,且最優(yōu)值大于速度系數(shù)法計(jì)算所得數(shù)值.

    圖10 不同進(jìn)口喉部寬度導(dǎo)葉性能的對(duì)比Fig.10 Performances comparison of diffusers with different inlet widths

    圖11為導(dǎo)葉A1,A0,A5和A7擴(kuò)散段橫截面的流線圖.可以看出,各導(dǎo)葉擴(kuò)散段流道內(nèi)都存在分離渦.導(dǎo)葉A1由于進(jìn)口喉部面積偏小,流道狹長(zhǎng),進(jìn)口流速過(guò)大,擴(kuò)散段內(nèi)的漩渦大小約占流道長(zhǎng)度的70%,流態(tài)較差,且流出擴(kuò)散段后,大部分液體仍然存在較大的切向速度,導(dǎo)葉損失較大.導(dǎo)葉A7則由于進(jìn)口喉部面積過(guò)大,且流道長(zhǎng)度過(guò)短,導(dǎo)致液體在導(dǎo)葉內(nèi)不能充分?jǐn)U散,經(jīng)過(guò)擴(kuò)散后的液體也存在較大的切向速度.導(dǎo)葉A0和A5的流態(tài)相對(duì)較好,液體流出擴(kuò)散段后切向速度較小,這有利于液體沿著軸向進(jìn)入下一級(jí)葉輪進(jìn)口端.因此,在設(shè)計(jì)LNG泵導(dǎo)葉時(shí)應(yīng)避免進(jìn)口喉部寬度取值過(guò)小,可選取比速度系數(shù)法計(jì)算的數(shù)值稍大的進(jìn)口喉部寬度.

    圖11 進(jìn)口喉部寬度不同時(shí)導(dǎo)葉擴(kuò)散段橫截面流線的對(duì)比Fig.11 Comparisons of cross section streamline for diffusers with different inlet widths

    3.4折轉(zhuǎn)角β的影響

    為分析進(jìn)口折轉(zhuǎn)角β對(duì)LNG泵導(dǎo)葉性能的影響,參照導(dǎo)葉A4的幾何結(jié)構(gòu)設(shè)計(jì)了導(dǎo)葉A8~A10(見(jiàn)表4).圖12為不同進(jìn)口折轉(zhuǎn)角導(dǎo)葉的結(jié)構(gòu)對(duì)比.可以看出,隨著β的減小,圓圈處的壁面厚度逐漸減小直至穿透,導(dǎo)葉的擴(kuò)散段長(zhǎng)度有所增加.

    表4 不同LNG泵導(dǎo)葉的折轉(zhuǎn)角Table 4 Turning angles of different LNG pump diffusers(°)

    圖12 導(dǎo)葉A4,A8和A9的對(duì)比Fig.12 Comparisons of diffusers A4,A8 and A9

    圖13為導(dǎo)葉進(jìn)口轉(zhuǎn)折角β不同時(shí)單級(jí)泵的揚(yáng)程和效率.可以看出,β對(duì)單級(jí)泵的揚(yáng)程和效率影響較小,揚(yáng)程變動(dòng)范圍為1 m以內(nèi).圖14為導(dǎo)葉A4,A8和A9擴(kuò)散段橫截面的流線圖.可以看出,導(dǎo)葉內(nèi)部流態(tài)較相似.鑒于進(jìn)口轉(zhuǎn)折角對(duì)LNG泵導(dǎo)葉的水力性能影響較小,在設(shè)計(jì)導(dǎo)葉時(shí),可以優(yōu)先確定其他關(guān)鍵尺寸,最后再通過(guò)調(diào)節(jié)進(jìn)口轉(zhuǎn)折角改善導(dǎo)葉的結(jié)構(gòu).

    圖13 不同進(jìn)口折轉(zhuǎn)角的導(dǎo)葉性能對(duì)比Fig.13 Performance comparisons of diffusers with different turning angles

    圖14 折轉(zhuǎn)角不同時(shí)導(dǎo)葉擴(kuò)散段橫截面流線的對(duì)比Fig.14 Comparisons of cross section streamline for diffusers with different turning angles

    4 結(jié)束語(yǔ)

    LNG泵的導(dǎo)葉不同于常規(guī)導(dǎo)葉,其結(jié)構(gòu)簡(jiǎn)單緊湊,徑向和軸向尺寸都相對(duì)較小,加工更加方便.本工作在分析LNG泵導(dǎo)葉結(jié)構(gòu)特點(diǎn)的基礎(chǔ)上,重點(diǎn)考慮進(jìn)口喉部寬度和折轉(zhuǎn)角對(duì)設(shè)計(jì)工況下泵水力性能的影響.通過(guò)設(shè)計(jì)不同幾何參數(shù)的LNG泵導(dǎo)葉,將設(shè)計(jì)導(dǎo)葉分別與同一葉輪進(jìn)行匹配,并借助數(shù)值模擬計(jì)算LNG泵導(dǎo)葉的性能.研究結(jié)果表明:進(jìn)口喉部寬度是LNG泵導(dǎo)葉的關(guān)鍵尺寸,設(shè)計(jì)時(shí)需重點(diǎn)考慮;進(jìn)口喉部寬度存在最優(yōu)值,且最優(yōu)值大于速度系數(shù)法計(jì)算的值;導(dǎo)葉進(jìn)口折轉(zhuǎn)角對(duì)泵揚(yáng)程影響較小,設(shè)計(jì)時(shí)可最后考慮通過(guò)改變折轉(zhuǎn)角改善導(dǎo)葉的結(jié)構(gòu).

    [1]祝勇仁,張煒,王循明.LNG汽車(chē)加氣站用潛液泵研制[J].機(jī)械科學(xué)與技術(shù),2012,31(1):163-166.

    [2]關(guān)醒凡.現(xiàn)代泵技術(shù)手冊(cè)[M].北京:宇航出版社,2010:315-321.

    [3]趙秋霞.導(dǎo)流殼幾何參數(shù)選取及其對(duì)泵性能的影響[J].太原理工大學(xué)學(xué)報(bào),2002,33(4):414-416.

    [4]孔繁余,王文廷,黃道見(jiàn),等.前置導(dǎo)葉調(diào)節(jié)混流泵性能的數(shù)值模擬[J].農(nóng)業(yè)工程學(xué)報(bào),2010,26(10):124-128.

    [5]劉建瑞,施衛(wèi)東,孔繁余,等.導(dǎo)葉渦殼組合式雙級(jí)泵的研究[J].農(nóng)業(yè)工程學(xué)報(bào),2005,21(9):80-84.

    [6]GOLcU M,PANcAR Y,SEKMEN Y.Energy saving in a deep well pump with splitter blade[J]. Energy Conversion Management,2006,47(5):638-651.

    [7]CULLEN D,RUSH S,MADISON J.Radial and axial diffusers for submerged electric motor driven pumps[J].World Pumps,2000,2000(408):24-29.

    [8]G¨ULIcH J F.Centrifugal pumps[M].Berlin:Springer-Verlag,2008:413.

    [9]ZHANG Q H,SHI W D,XU Y,et al.A new proposed return guide vane for compact multistage centrifugal pumps[J].International Journal of Rotating Machinery,2013,2013(15):1-11.

    [10]徐媛暉,張啟華,施衛(wèi)東,等.基于數(shù)值模擬的新型軸向?qū)~性能分析[J].流體機(jī)械,2014,42(9):25-30.

    [11]周嶺,施衛(wèi)東,陸偉剛.基于數(shù)值模擬的深井離心泵導(dǎo)葉性能分析[J].農(nóng)業(yè)工程學(xué)報(bào),2011,27(9):38-42.

    [12]叢小青,王光輝,袁丹青,等.空間導(dǎo)葉式離心泵的數(shù)值計(jì)算及優(yōu)化設(shè)計(jì)[J].排灌機(jī)械工程學(xué)報(bào), 2010,28(6):488-491.

    [13]李家文,朱森元,劉中祥.多級(jí)液氫泵級(jí)間導(dǎo)葉的設(shè)計(jì)與改進(jìn)[J].宇航學(xué)報(bào),2006,27(5):1072-1075.

    本文彩色版可登陸本刊網(wǎng)站查詢:http://www.journal.shu.edu.cn

    Design of liquefied natural gas submerged pump guide vane based on CFD

    QIAN Tao,CHEN Hongxun,LIANG Chengpeng
    (Shanghai Institute of Applied Mathematics and Mechanics,Shanghai University, Shanghai 200072,China)

    Liquefied natural gas(LNG)submerged pump is designed to immerse its canned motor and pump body into cryogenic liquid.To reduce radial and axial sizes,special guide vanes with different structures from popular ones are used.By analyzing the LNG pump structure,two key geometric parameters,inlet width and turning angle,were investigated to see how they affect pump hydraulic performance in the design.LNG pump guide vanes with different inlet widths and turning angles were designed to assemble the same impeller. Three-dimensional turbulent flows of different models were then numerically simulated using a standard k-ε turbulence model with ANSYS CFX.The comparative study revealed that optimum value of inlet width exited,and was larger than experience-based values. Inlet width had great influence on the pump performance,while turning angle was much less important.Therefore,designing an LNG pump guide vane,inlet width should be given high priority,and turning angle should be adjustable for structural design.

    liquefied natural gas(LNG)submerged pump;guide vane;inlet width;turning angle;numerical simulation

    TH 311

    A

    1007-2861(2016)05-0606-10

    10.3969/j.issn.1007-2861.2015.01.016

    2015-03-18

    陳紅勛(1962—),男,研究員,博士生導(dǎo)師,博士,研究方向?yàn)榱黧w機(jī)械、計(jì)算水動(dòng)力學(xué). E-mail:chenhx@shu.edu.cn

    猜你喜歡
    喉部導(dǎo)葉轉(zhuǎn)角
    軸排凝汽器喉部設(shè)計(jì)
    玩轉(zhuǎn)角的平分線
    電子喉鏡聯(lián)合窄帶成像技術(shù)對(duì)喉部早期惡性病變的診斷價(jià)值研討
    壓氣機(jī)導(dǎo)葉用GTD222合金鑄造技術(shù)研究
    模具制造(2019年4期)2019-06-24 03:36:50
    凝汽器喉部流場(chǎng)數(shù)值模擬
    三次“轉(zhuǎn)角”遇到愛(ài)
    解放軍健康(2017年5期)2017-08-01 06:27:42
    永春堂贏在轉(zhuǎn)角
    咽及喉部鱗癌放射治療技術(shù)研究進(jìn)展
    下一個(gè)轉(zhuǎn)角:邁出去 開(kāi)啟“智”造時(shí)代
    混流式模型水輪機(jī)導(dǎo)葉水力矩試驗(yàn)
    能在线免费观看的黄片| 三级国产精品欧美在线观看| 一级a爱片免费观看的视频| 长腿黑丝高跟| 精品一区二区三区视频在线| 成人美女网站在线观看视频| 三级国产精品欧美在线观看| 精品久久国产蜜桃| 午夜激情欧美在线| 国内精品久久久久精免费| 黄色日韩在线| 婷婷色综合大香蕉| 免费搜索国产男女视频| 在线播放国产精品三级| 亚洲电影在线观看av| 菩萨蛮人人尽说江南好唐韦庄 | 国产黄片美女视频| 国产伦精品一区二区三区视频9| 久久久精品欧美日韩精品| 国产人妻一区二区三区在| 国产av一区在线观看免费| 色综合色国产| 国产不卡一卡二| 欧美日韩乱码在线| 精品午夜福利视频在线观看一区| 大香蕉久久网| 一a级毛片在线观看| 一个人免费在线观看电影| 三级毛片av免费| 啦啦啦观看免费观看视频高清| 国产熟女欧美一区二区| 国产精品国产三级国产av玫瑰| 美女内射精品一级片tv| 久久久久久久久久久丰满| 国产精品一区二区三区四区免费观看 | 天美传媒精品一区二区| 国产精品久久电影中文字幕| 国产高清不卡午夜福利| 亚州av有码| 亚洲精品国产av成人精品 | 欧美精品国产亚洲| 国产黄色视频一区二区在线观看 | 99久久成人亚洲精品观看| 婷婷亚洲欧美| 日日撸夜夜添| 久久久色成人| 大又大粗又爽又黄少妇毛片口| 午夜福利在线观看免费完整高清在 | 亚洲精品粉嫩美女一区| 国产成人aa在线观看| 啦啦啦韩国在线观看视频| 我的老师免费观看完整版| 99热只有精品国产| 国产精品一区www在线观看| 99热精品在线国产| 又粗又爽又猛毛片免费看| 久久久精品94久久精品| 少妇人妻精品综合一区二区 | 99久久九九国产精品国产免费| 美女大奶头视频| 我要搜黄色片| 亚洲不卡免费看| 亚洲精品一卡2卡三卡4卡5卡| 国内揄拍国产精品人妻在线| 美女cb高潮喷水在线观看| 久久久久精品国产欧美久久久| av福利片在线观看| 免费黄网站久久成人精品| 欧美中文日本在线观看视频| 亚洲aⅴ乱码一区二区在线播放| 色视频www国产| 美女内射精品一级片tv| 亚洲熟妇中文字幕五十中出| 少妇丰满av| 内射极品少妇av片p| 尾随美女入室| 日日摸夜夜添夜夜添小说| 国产欧美日韩精品一区二区| 亚洲人成网站高清观看| 亚洲在线观看片| 亚洲三级黄色毛片| 国产大屁股一区二区在线视频| 亚洲在线观看片| 天堂av国产一区二区熟女人妻| 国产在线男女| av在线亚洲专区| 中国美白少妇内射xxxbb| 成人永久免费在线观看视频| 97超碰精品成人国产| 黄色配什么色好看| 黄色一级大片看看| 在线观看66精品国产| 在线天堂最新版资源| 国产精品一区二区性色av| 欧美成人一区二区免费高清观看| 亚洲精品粉嫩美女一区| 久久久久久久久久黄片| 丰满人妻一区二区三区视频av| 国产蜜桃级精品一区二区三区| 草草在线视频免费看| 久久这里只有精品中国| 久久亚洲精品不卡| 久久婷婷人人爽人人干人人爱| 精品国产三级普通话版| 精品一区二区免费观看| 免费电影在线观看免费观看| 欧美一区二区国产精品久久精品| 我要搜黄色片| 乱人视频在线观看| 免费看美女性在线毛片视频| 天美传媒精品一区二区| 日韩成人av中文字幕在线观看 | 国产欧美日韩精品一区二区| 午夜视频国产福利| 国产精品人妻久久久久久| 内射极品少妇av片p| 看十八女毛片水多多多| 日本成人三级电影网站| 变态另类成人亚洲欧美熟女| 精品不卡国产一区二区三区| 国产亚洲精品久久久久久毛片| 日本黄大片高清| 人人妻,人人澡人人爽秒播| 日韩精品有码人妻一区| 久久精品国产亚洲av香蕉五月| 日韩欧美三级三区| 国产一区二区亚洲精品在线观看| 成年女人永久免费观看视频| 偷拍熟女少妇极品色| 在线观看美女被高潮喷水网站| 天天躁日日操中文字幕| 一本久久中文字幕| 可以在线观看的亚洲视频| 乱系列少妇在线播放| 一本一本综合久久| 久久精品夜色国产| 精品午夜福利视频在线观看一区| 久久精品国产鲁丝片午夜精品| 欧美最黄视频在线播放免费| 成人永久免费在线观看视频| 三级经典国产精品| 国内精品美女久久久久久| 波多野结衣高清无吗| 男人舔女人下体高潮全视频| 91麻豆精品激情在线观看国产| 女生性感内裤真人,穿戴方法视频| 成人漫画全彩无遮挡| 两个人视频免费观看高清| 午夜免费激情av| 秋霞在线观看毛片| 校园春色视频在线观看| 老师上课跳d突然被开到最大视频| 国产欧美日韩精品一区二区| 伦精品一区二区三区| 一个人看的www免费观看视频| 寂寞人妻少妇视频99o| 午夜影院日韩av| 我的女老师完整版在线观看| 嫩草影院精品99| 欧美zozozo另类| 99久久精品国产国产毛片| 麻豆乱淫一区二区| 亚洲图色成人| 亚洲人成网站在线观看播放| 少妇熟女aⅴ在线视频| 欧美三级亚洲精品| 欧美又色又爽又黄视频| 午夜影院日韩av| 99热网站在线观看| 黄色视频,在线免费观看| 岛国在线免费视频观看| 亚洲四区av| 桃色一区二区三区在线观看| 99久久成人亚洲精品观看| 国产欧美日韩精品一区二区| 三级经典国产精品| 亚洲激情五月婷婷啪啪| 国产成人91sexporn| 97在线视频观看| 国模一区二区三区四区视频| 51国产日韩欧美| 蜜桃亚洲精品一区二区三区| 久久久精品大字幕| 伦精品一区二区三区| 你懂的网址亚洲精品在线观看 | 日韩欧美精品免费久久| 成人性生交大片免费视频hd| 精品久久久久久成人av| 亚洲丝袜综合中文字幕| 国产成人freesex在线 | 男插女下体视频免费在线播放| 久久精品国产亚洲网站| 日本在线视频免费播放| a级一级毛片免费在线观看| 亚洲色图av天堂| 99久久久亚洲精品蜜臀av| 少妇丰满av| 少妇猛男粗大的猛烈进出视频 | 国产视频一区二区在线看| 97碰自拍视频| 一进一出好大好爽视频| 舔av片在线| 高清毛片免费观看视频网站| aaaaa片日本免费| 成人美女网站在线观看视频| 大香蕉久久网| av国产免费在线观看| a级毛片a级免费在线| 午夜免费激情av| 久久久a久久爽久久v久久| 久久久久免费精品人妻一区二区| 白带黄色成豆腐渣| 亚洲精品一区av在线观看| 国产精品一区二区三区四区久久| 亚洲欧美日韩卡通动漫| 女人被狂操c到高潮| 真人做人爱边吃奶动态| 国产黄片美女视频| 国产 一区精品| 少妇的逼好多水| 18禁在线无遮挡免费观看视频 | 欧美日本亚洲视频在线播放| 亚洲熟妇中文字幕五十中出| 午夜激情欧美在线| 亚洲精华国产精华液的使用体验 | 一进一出抽搐gif免费好疼| 久久久成人免费电影| 午夜亚洲福利在线播放| 麻豆国产av国片精品| 亚洲精品乱码久久久v下载方式| 在线播放无遮挡| 老女人水多毛片| 看黄色毛片网站| 精品不卡国产一区二区三区| 少妇人妻一区二区三区视频| 中国国产av一级| 国产精品嫩草影院av在线观看| 直男gayav资源| 亚洲最大成人中文| 老熟妇仑乱视频hdxx| 亚洲综合色惰| 亚洲av一区综合| 久久人人爽人人片av| 免费人成在线观看视频色| 亚洲一区二区三区色噜噜| 日韩三级伦理在线观看| 在线观看午夜福利视频| 一卡2卡三卡四卡精品乱码亚洲| 国产精品99久久久久久久久| 蜜臀久久99精品久久宅男| 变态另类丝袜制服| 日韩高清综合在线| 亚洲精品久久国产高清桃花| 秋霞在线观看毛片| 亚洲精品456在线播放app| 国产成人a区在线观看| 国产黄片美女视频| 97人妻精品一区二区三区麻豆| 久久久精品欧美日韩精品| 午夜爱爱视频在线播放| 精品国内亚洲2022精品成人| 色哟哟哟哟哟哟| 十八禁网站免费在线| 久久韩国三级中文字幕| 欧美成人免费av一区二区三区| 亚洲欧美成人精品一区二区| 国产一级毛片七仙女欲春2| 我的女老师完整版在线观看| 最近手机中文字幕大全| 亚洲第一电影网av| 久久鲁丝午夜福利片| 色在线成人网| 中文字幕熟女人妻在线| 韩国av在线不卡| 欧美高清成人免费视频www| 日本免费一区二区三区高清不卡| 丰满人妻一区二区三区视频av| 一区二区三区免费毛片| 国产91av在线免费观看| 黄色配什么色好看| 国产不卡一卡二| 久久韩国三级中文字幕| 亚洲国产精品久久男人天堂| 人人妻人人澡人人爽人人夜夜 | 一个人看视频在线观看www免费| 最好的美女福利视频网| 亚洲av一区综合| 又爽又黄无遮挡网站| 人妻制服诱惑在线中文字幕| 久久精品国产自在天天线| av卡一久久| 亚洲av电影不卡..在线观看| 国内精品一区二区在线观看| 18禁黄网站禁片免费观看直播| 免费观看人在逋| 成人二区视频| 国产伦在线观看视频一区| 国产精品无大码| а√天堂www在线а√下载| 国产91av在线免费观看| 看十八女毛片水多多多| 国产精品av视频在线免费观看| 欧美zozozo另类| 精品午夜福利视频在线观看一区| 大香蕉久久网| 搡女人真爽免费视频火全软件 | 天堂影院成人在线观看| 亚洲中文字幕一区二区三区有码在线看| a级毛片免费高清观看在线播放| 亚洲国产精品成人久久小说 | 天堂网av新在线| 国产精品综合久久久久久久免费| 欧美高清成人免费视频www| 久久久午夜欧美精品| 日本免费a在线| 老司机福利观看| 小说图片视频综合网站| 久久久久久久久大av| 日本一二三区视频观看| 99热精品在线国产| 亚洲中文字幕一区二区三区有码在线看| 床上黄色一级片| 97超级碰碰碰精品色视频在线观看| 国产精品一区二区三区四区免费观看 | 日日撸夜夜添| 最近手机中文字幕大全| 一级毛片aaaaaa免费看小| 国产真实伦视频高清在线观看| 91在线精品国自产拍蜜月| 国产色爽女视频免费观看| 免费无遮挡裸体视频| 久久亚洲精品不卡| 久久久久九九精品影院| 欧美潮喷喷水| 男女啪啪激烈高潮av片| 一夜夜www| 校园春色视频在线观看| 欧美区成人在线视频| 久久精品影院6| aaaaa片日本免费| 久久久久国产精品人妻aⅴ院| 久久久久久伊人网av| 一本一本综合久久| av专区在线播放| 久久久午夜欧美精品| 欧美日韩乱码在线| 久久草成人影院| 亚洲成人av在线免费| 国产精品免费一区二区三区在线| 成人av在线播放网站| 成人精品一区二区免费| av中文乱码字幕在线| 观看美女的网站| 直男gayav资源| 人妻夜夜爽99麻豆av| 中文字幕熟女人妻在线| 人妻久久中文字幕网| 夜夜看夜夜爽夜夜摸| 有码 亚洲区| 性欧美人与动物交配| 久久久久国产精品人妻aⅴ院| 男女做爰动态图高潮gif福利片| 久99久视频精品免费| 午夜免费男女啪啪视频观看 | 国产人妻一区二区三区在| 久久久欧美国产精品| 日韩欧美精品v在线| 国产成人a∨麻豆精品| 免费看日本二区| 看十八女毛片水多多多| 十八禁网站免费在线| 最新在线观看一区二区三区| 国产亚洲精品综合一区在线观看| 少妇丰满av| 99在线视频只有这里精品首页| 搡老熟女国产l中国老女人| 在线观看一区二区三区| 毛片一级片免费看久久久久| 国产高清三级在线| 亚洲成人中文字幕在线播放| 简卡轻食公司| 伊人久久精品亚洲午夜| www.色视频.com| 久久久色成人| 国产亚洲精品av在线| 久久99热这里只有精品18| 男人和女人高潮做爰伦理| 在线播放无遮挡| 五月伊人婷婷丁香| 久久久久国内视频| 赤兔流量卡办理| 午夜a级毛片| 国产精品美女特级片免费视频播放器| 特级一级黄色大片| 免费av毛片视频| 久久精品夜色国产| 成年版毛片免费区| 亚洲精品色激情综合| 午夜激情欧美在线| 亚洲av二区三区四区| 综合色av麻豆| 一本精品99久久精品77| 亚洲成人中文字幕在线播放| 男插女下体视频免费在线播放| 99久久精品一区二区三区| 国产三级在线视频| 亚洲精品国产av成人精品 | 国产精品一区二区三区四区久久| 麻豆精品久久久久久蜜桃| 99久久九九国产精品国产免费| 男女下面进入的视频免费午夜| 久久九九热精品免费| 欧美日本视频| 久久久成人免费电影| 1024手机看黄色片| 亚洲,欧美,日韩| 欧美精品国产亚洲| 女人被狂操c到高潮| 一级黄色大片毛片| 欧美一区二区国产精品久久精品| 国内久久婷婷六月综合欲色啪| 午夜福利在线在线| 亚洲精品国产av成人精品 | 舔av片在线| 国产高清三级在线| 亚洲精品日韩av片在线观看| 日本欧美国产在线视频| 人妻制服诱惑在线中文字幕| 哪里可以看免费的av片| 久久精品国产亚洲av天美| 日日啪夜夜撸| 午夜福利视频1000在线观看| 国产探花极品一区二区| 亚洲熟妇中文字幕五十中出| 亚洲av电影不卡..在线观看| 久久热精品热| 97热精品久久久久久| 久久亚洲国产成人精品v| 欧美精品国产亚洲| 日本在线视频免费播放| 此物有八面人人有两片| 97热精品久久久久久| av在线亚洲专区| 亚洲熟妇熟女久久| 内地一区二区视频在线| 人人妻,人人澡人人爽秒播| 亚洲av免费高清在线观看| 日本熟妇午夜| 有码 亚洲区| 天堂影院成人在线观看| 亚洲美女黄片视频| 国产欧美日韩一区二区精品| 婷婷亚洲欧美| or卡值多少钱| 淫妇啪啪啪对白视频| 日本 av在线| 干丝袜人妻中文字幕| 91久久精品电影网| 国产亚洲精品综合一区在线观看| 免费在线观看影片大全网站| 亚洲va在线va天堂va国产| 亚洲最大成人中文| 日本免费a在线| 亚洲性久久影院| 精品久久国产蜜桃| 国产在视频线在精品| 国产精品久久久久久亚洲av鲁大| 久久天躁狠狠躁夜夜2o2o| 日韩国内少妇激情av| 男女边吃奶边做爰视频| 乱码一卡2卡4卡精品| 亚洲精品日韩av片在线观看| 热99re8久久精品国产| av在线蜜桃| 久久久久九九精品影院| 少妇高潮的动态图| 欧美一区二区亚洲| 国产精品野战在线观看| 国产成人91sexporn| 精品日产1卡2卡| 午夜影院日韩av| 51国产日韩欧美| 麻豆久久精品国产亚洲av| 一级毛片我不卡| 少妇高潮的动态图| 一个人免费在线观看电影| 天堂av国产一区二区熟女人妻| 我的女老师完整版在线观看| 精品99又大又爽又粗少妇毛片| 人人妻人人澡欧美一区二区| 最近2019中文字幕mv第一页| 免费观看在线日韩| 国模一区二区三区四区视频| 亚洲国产精品合色在线| 免费观看人在逋| 国产亚洲精品综合一区在线观看| 国产国拍精品亚洲av在线观看| 亚洲第一电影网av| 色噜噜av男人的天堂激情| 99久国产av精品| 国产午夜精品久久久久久一区二区三区 | 日本免费a在线| 中文字幕av在线有码专区| 99热只有精品国产| 日本在线视频免费播放| 婷婷色综合大香蕉| 97超级碰碰碰精品色视频在线观看| 日韩精品青青久久久久久| 熟女人妻精品中文字幕| 永久网站在线| 给我免费播放毛片高清在线观看| 欧美日韩综合久久久久久| 丝袜美腿在线中文| 天天一区二区日本电影三级| 久久综合国产亚洲精品| 波多野结衣巨乳人妻| 国产一区二区亚洲精品在线观看| 亚洲在线观看片| 国产高清视频在线观看网站| 亚洲av成人精品一区久久| 99久久精品热视频| 亚洲精品乱码久久久v下载方式| 18禁黄网站禁片免费观看直播| 女人被狂操c到高潮| 村上凉子中文字幕在线| 日本黄色片子视频| 日产精品乱码卡一卡2卡三| 日韩 亚洲 欧美在线| 99久久久亚洲精品蜜臀av| 国产欧美日韩一区二区精品| 久久精品91蜜桃| videossex国产| 国语自产精品视频在线第100页| 国产高清有码在线观看视频| 成人永久免费在线观看视频| or卡值多少钱| 亚洲国产日韩欧美精品在线观看| 成年女人永久免费观看视频| 美女被艹到高潮喷水动态| 国产高清视频在线观看网站| 精品久久久久久久久久免费视频| 国产精品国产三级国产av玫瑰| 精品久久久久久久久久久久久| 一进一出好大好爽视频| av黄色大香蕉| 亚洲av免费在线观看| 男女那种视频在线观看| 香蕉av资源在线| 国产又黄又爽又无遮挡在线| 黑人高潮一二区| 99热这里只有是精品50| 十八禁网站免费在线| 国产精品野战在线观看| 18禁在线播放成人免费| 日韩成人伦理影院| 俄罗斯特黄特色一大片| 国产伦在线观看视频一区| 成人漫画全彩无遮挡| 熟妇人妻久久中文字幕3abv| 国产一区二区三区在线臀色熟女| 免费一级毛片在线播放高清视频| 国产黄a三级三级三级人| 国产在线男女| 久久久精品欧美日韩精品| 免费av观看视频| 淫秽高清视频在线观看| av在线亚洲专区| 午夜亚洲福利在线播放| 亚洲欧美日韩东京热| 久久韩国三级中文字幕| 亚洲成人av在线免费| 亚洲欧美清纯卡通| 日本黄大片高清| 狂野欧美白嫩少妇大欣赏| av在线观看视频网站免费| 国产单亲对白刺激| 少妇丰满av| 亚洲天堂国产精品一区在线| 国产精品一区二区三区四区久久| 美女被艹到高潮喷水动态| 露出奶头的视频| 国产精品亚洲美女久久久| 搡老妇女老女人老熟妇| 日韩精品中文字幕看吧| 久久久久国内视频| 1000部很黄的大片| 麻豆成人午夜福利视频| 嫩草影视91久久| 免费看av在线观看网站| 国产不卡一卡二| 天堂影院成人在线观看| 久久久久久久久久黄片| 99精品在免费线老司机午夜| 露出奶头的视频| 日本熟妇午夜| 成人无遮挡网站| 国产精品一区二区性色av| 最近在线观看免费完整版| 国产高潮美女av| 在线免费十八禁| 全区人妻精品视频| 欧美成人a在线观看| 久久亚洲国产成人精品v| 国产欧美日韩精品亚洲av| 久久久久久久亚洲中文字幕| 美女大奶头视频| 国产欧美日韩一区二区精品| 国产白丝娇喘喷水9色精品| 欧美一区二区亚洲| 人人妻人人看人人澡| 久久人人爽人人爽人人片va| 天堂av国产一区二区熟女人妻| 久久久久久久午夜电影| 欧美人与善性xxx| 性插视频无遮挡在线免费观看|