• <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)
    xxx大片免费视频| 成人二区视频| 狂野欧美白嫩少妇大欣赏| 亚洲va在线va天堂va国产| 2018国产大陆天天弄谢| 97人妻精品一区二区三区麻豆| 我的老师免费观看完整版| 麻豆成人av视频| 成人午夜精彩视频在线观看| 啦啦啦啦在线视频资源| 一区二区三区乱码不卡18| 能在线免费观看的黄片| 99热这里只有是精品在线观看| 美女高潮的动态| 亚洲av中文字字幕乱码综合| 大陆偷拍与自拍| 永久免费av网站大全| 自拍偷自拍亚洲精品老妇| 国产黄频视频在线观看| 亚洲av.av天堂| 97在线视频观看| 日日干狠狠操夜夜爽| 性色avwww在线观看| 国产人妻一区二区三区在| 男人狂女人下面高潮的视频| 天堂网av新在线| 尾随美女入室| 成人毛片a级毛片在线播放| 国产三级在线视频| 日韩欧美三级三区| av在线播放精品| 2021少妇久久久久久久久久久| 日韩一区二区三区影片| 国产高清三级在线| 国产午夜精品论理片| 精品国产一区二区三区久久久樱花 | 777米奇影视久久| 深爱激情五月婷婷| 成人av在线播放网站| 99久国产av精品| 乱码一卡2卡4卡精品| 亚洲人成网站在线观看播放| 亚洲国产高清在线一区二区三| av天堂中文字幕网| 日本与韩国留学比较| 七月丁香在线播放| 在线观看美女被高潮喷水网站| 久久午夜福利片| 欧美xxxx性猛交bbbb| 亚洲国产成人一精品久久久| 美女被艹到高潮喷水动态| 九草在线视频观看| 日日啪夜夜撸| 高清av免费在线| 在现免费观看毛片| 亚洲人成网站在线播| 成人欧美大片| 黄色配什么色好看| 日日啪夜夜撸| 国产av不卡久久| 亚洲欧美日韩卡通动漫| 一级毛片我不卡| 非洲黑人性xxxx精品又粗又长| av线在线观看网站| 日韩一区二区视频免费看| 黄色一级大片看看| 婷婷色麻豆天堂久久| 亚洲av.av天堂| 国产精品精品国产色婷婷| videossex国产| 午夜激情福利司机影院| 成人美女网站在线观看视频| 国产淫语在线视频| 欧美日韩综合久久久久久| 日日啪夜夜撸| 亚洲真实伦在线观看| 精品一区二区免费观看| 欧美zozozo另类| 久久久色成人| 亚洲图色成人| 欧美激情久久久久久爽电影| 国产黄频视频在线观看| 亚洲乱码一区二区免费版| 精品一区二区三区视频在线| 精品99又大又爽又粗少妇毛片| 久久久精品免费免费高清| 老司机影院毛片| 亚洲精品视频女| 国产精品一及| 亚洲国产欧美人成| 中文乱码字字幕精品一区二区三区 | 美女被艹到高潮喷水动态| 白带黄色成豆腐渣| 2022亚洲国产成人精品| 国产精品一及| 能在线免费看毛片的网站| 日本黄色片子视频| 久久久久久久久久久丰满| 九九爱精品视频在线观看| 直男gayav资源| 久久久久久久久中文| 精品少妇黑人巨大在线播放| 国产高潮美女av| 国产伦精品一区二区三区视频9| 97在线视频观看| videossex国产| 久久草成人影院| 成人一区二区视频在线观看| 国产有黄有色有爽视频| 国内揄拍国产精品人妻在线| 中国国产av一级| 国内精品宾馆在线| 精品人妻一区二区三区麻豆| 国产成人精品福利久久| 日本午夜av视频| 日本av手机在线免费观看| 777米奇影视久久| 床上黄色一级片| 一级片'在线观看视频| 欧美性感艳星| 亚洲精品色激情综合| 亚洲精品乱久久久久久| 亚洲av福利一区| 免费观看a级毛片全部| 寂寞人妻少妇视频99o| 三级经典国产精品| 不卡视频在线观看欧美| av在线亚洲专区| kizo精华| 午夜免费观看性视频| 欧美成人精品欧美一级黄| 亚洲成人一二三区av| 国产乱来视频区| 午夜精品在线福利| 一本久久精品| 日本爱情动作片www.在线观看| 亚洲精品乱久久久久久| 2022亚洲国产成人精品| 国产免费福利视频在线观看| 国产高潮美女av| 在线观看一区二区三区| 精品一区在线观看国产| 男人和女人高潮做爰伦理| 亚洲av电影不卡..在线观看| 伦精品一区二区三区| 成人毛片60女人毛片免费| 高清日韩中文字幕在线| 卡戴珊不雅视频在线播放| 久久精品久久久久久久性| 夜夜爽夜夜爽视频| 白带黄色成豆腐渣| 97人妻精品一区二区三区麻豆| 日本猛色少妇xxxxx猛交久久| 国产精品久久视频播放| 99久国产av精品| kizo精华| 91久久精品电影网| 2018国产大陆天天弄谢| 亚洲国产色片| 国产综合懂色| .国产精品久久| 成人亚洲精品一区在线观看 | 国产黄频视频在线观看| 亚洲熟妇中文字幕五十中出| 三级毛片av免费| 国产成人a∨麻豆精品| 亚洲精品色激情综合| 青春草视频在线免费观看| 久久精品熟女亚洲av麻豆精品 | av在线亚洲专区| 成人一区二区视频在线观看| 国产综合懂色| 联通29元200g的流量卡| av在线观看视频网站免费| 国产av国产精品国产| 中文在线观看免费www的网站| 男女啪啪激烈高潮av片| 嫩草影院入口| 亚洲精品一二三| 草草在线视频免费看| 国产片特级美女逼逼视频| 亚洲,欧美,日韩| 亚洲伊人久久精品综合| 少妇人妻精品综合一区二区| 久久久精品免费免费高清| ponron亚洲| 爱豆传媒免费全集在线观看| 五月天丁香电影| 国产色婷婷99| av在线亚洲专区| 白带黄色成豆腐渣| 久久99热6这里只有精品| 久久99热6这里只有精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲天堂国产精品一区在线| 天堂影院成人在线观看| 国产色爽女视频免费观看| 男人狂女人下面高潮的视频| 国产老妇伦熟女老妇高清| 三级国产精品片| 男人爽女人下面视频在线观看| 日韩 亚洲 欧美在线| 午夜免费激情av| 久久国内精品自在自线图片| 五月天丁香电影| 18禁动态无遮挡网站| 国产色爽女视频免费观看| 日韩成人伦理影院| av天堂中文字幕网| 国产精品国产三级专区第一集| 国产女主播在线喷水免费视频网站 | 免费人成在线观看视频色| 最近最新中文字幕大全电影3| 国产淫语在线视频| 亚洲色图av天堂| 成人毛片a级毛片在线播放| 国产一区二区三区av在线| av天堂中文字幕网| 全区人妻精品视频| 网址你懂的国产日韩在线| 国产精品1区2区在线观看.| 久久久久久久久久成人| 18+在线观看网站| 国产伦在线观看视频一区| 成年女人在线观看亚洲视频 | 一级毛片黄色毛片免费观看视频| 青春草视频在线免费观看| 免费少妇av软件| kizo精华| 精品国产露脸久久av麻豆 | 熟女电影av网| 国产成人精品久久久久久| 蜜臀久久99精品久久宅男| 菩萨蛮人人尽说江南好唐韦庄| 国产成人福利小说| 91在线精品国自产拍蜜月| 亚洲欧美中文字幕日韩二区| 非洲黑人性xxxx精品又粗又长| 国产高清三级在线| 亚洲精品456在线播放app| 日韩一区二区三区影片| 99久国产av精品| 欧美激情久久久久久爽电影| 国产不卡一卡二| 91在线精品国自产拍蜜月| 中文字幕制服av| 黄色配什么色好看| 亚洲国产精品专区欧美| 好男人视频免费观看在线| 亚洲天堂国产精品一区在线| 成人毛片60女人毛片免费| 男女边吃奶边做爰视频| 成人av在线播放网站| 亚洲精华国产精华液的使用体验| 亚洲精品,欧美精品| 搡老妇女老女人老熟妇| 成年av动漫网址| 网址你懂的国产日韩在线| 亚洲精品视频女| 日本免费a在线| 91精品国产九色| 国产成人freesex在线| 国产一区二区三区av在线| 亚洲欧美日韩卡通动漫| 久久久久久久久久久丰满| 熟女电影av网| 91久久精品国产一区二区三区| 久久精品国产亚洲网站| 看十八女毛片水多多多| 色吧在线观看| 久久久久网色| 熟妇人妻久久中文字幕3abv| 亚洲av福利一区| 精品人妻熟女av久视频| 精品久久久噜噜| 久久久久久九九精品二区国产| 国产一区亚洲一区在线观看| 亚洲精品第二区| 国产又色又爽无遮挡免| 亚洲欧美精品专区久久| 久久久精品欧美日韩精品| 日韩 亚洲 欧美在线| 欧美日本视频| eeuss影院久久| 毛片一级片免费看久久久久| 亚洲精品一区蜜桃| 女人久久www免费人成看片| 少妇被粗大猛烈的视频| 丝袜美腿在线中文| 99久久人妻综合| 国产免费福利视频在线观看| 日本wwww免费看| 丰满少妇做爰视频| 久久精品国产鲁丝片午夜精品| 99久久九九国产精品国产免费| 国产黄色免费在线视频| 国产精品女同一区二区软件| 一级黄片播放器| 亚洲欧美成人精品一区二区| 精品一区二区三区人妻视频| xxx大片免费视频| 亚洲国产色片| 黄色欧美视频在线观看| 寂寞人妻少妇视频99o| 国产精品伦人一区二区| 久久久精品免费免费高清| 91av网一区二区| 性色avwww在线观看| 波多野结衣巨乳人妻| 一区二区三区高清视频在线| 国国产精品蜜臀av免费| 99热全是精品| 国产男女超爽视频在线观看| 亚洲内射少妇av| 国产69精品久久久久777片| 在线播放无遮挡| 亚洲熟女精品中文字幕| 久久久久精品性色| 波野结衣二区三区在线| 成人漫画全彩无遮挡| av黄色大香蕉| 久久精品人妻少妇| 视频中文字幕在线观看| 美女高潮的动态| 国产精品.久久久| 国国产精品蜜臀av免费| 欧美日韩亚洲高清精品| h日本视频在线播放| 免费电影在线观看免费观看| 国产一级毛片在线| 国产精品无大码| 成年人午夜在线观看视频 | 99久久中文字幕三级久久日本| 日本三级黄在线观看| 黄片wwwwww| 啦啦啦啦在线视频资源| 深夜a级毛片| 三级国产精品片| 日韩精品青青久久久久久| 日韩亚洲欧美综合| 国内精品美女久久久久久| 最后的刺客免费高清国语| 人人妻人人澡人人爽人人夜夜 | 中文资源天堂在线| 我要看日韩黄色一级片| 欧美日本视频| 综合色丁香网| 国产精品不卡视频一区二区| 国产综合精华液| 91狼人影院| 亚洲最大成人av| 午夜免费观看性视频| 少妇的逼好多水| 青春草视频在线免费观看| 亚洲av福利一区| 亚洲在线观看片| 最近最新中文字幕免费大全7| 欧美高清性xxxxhd video| 男女下面进入的视频免费午夜| 18禁动态无遮挡网站| 免费观看a级毛片全部| 男人狂女人下面高潮的视频| 3wmmmm亚洲av在线观看| 亚洲精品乱码久久久v下载方式| 男女下面进入的视频免费午夜| 中文字幕制服av| 美女主播在线视频| 中文字幕免费在线视频6| 国产不卡一卡二| 亚洲精品国产av蜜桃| 免费观看a级毛片全部| 国产精品伦人一区二区| 最新中文字幕久久久久| av黄色大香蕉| 五月玫瑰六月丁香| 久久这里有精品视频免费| 色哟哟·www| 日韩电影二区| 国产视频首页在线观看| 欧美人与善性xxx| 亚洲av福利一区| 亚洲美女搞黄在线观看| 男女啪啪激烈高潮av片| 成人午夜高清在线视频| 中文字幕人妻熟人妻熟丝袜美| 久久99精品国语久久久| 狂野欧美激情性xxxx在线观看| 成人亚洲精品av一区二区| 欧美区成人在线视频| 国产精品国产三级国产av玫瑰| 午夜免费激情av| 尾随美女入室| 国产精品精品国产色婷婷| 亚洲无线观看免费| 搡老妇女老女人老熟妇| 亚洲色图av天堂| 亚洲婷婷狠狠爱综合网| 成年女人看的毛片在线观看| 乱系列少妇在线播放| 亚洲精品亚洲一区二区| 色视频www国产| 少妇人妻一区二区三区视频| 中文字幕制服av| 亚洲av成人精品一二三区| 久久精品国产鲁丝片午夜精品| 亚洲av男天堂| 亚洲在线自拍视频| a级毛片免费高清观看在线播放| 免费看不卡的av| 免费看美女性在线毛片视频| 伦理电影大哥的女人| 午夜福利在线观看吧| 色网站视频免费| 大香蕉97超碰在线| 亚洲精品乱码久久久久久按摩| 久99久视频精品免费| 日韩av不卡免费在线播放| 青春草国产在线视频| 91av网一区二区| 99热这里只有精品一区| 高清日韩中文字幕在线| 久久精品熟女亚洲av麻豆精品 | 观看免费一级毛片| 国产精品日韩av在线免费观看| 国产精品不卡视频一区二区| 女人十人毛片免费观看3o分钟| 欧美不卡视频在线免费观看| 免费观看的影片在线观看| 三级国产精品欧美在线观看| 久久精品国产亚洲网站| 国产一区二区在线观看日韩| 老女人水多毛片| 九色成人免费人妻av| 久久久久久久久中文| 在线观看美女被高潮喷水网站| 国产毛片a区久久久久| 能在线免费看毛片的网站| 国内精品美女久久久久久| 狂野欧美白嫩少妇大欣赏| 亚洲欧美精品专区久久| 色尼玛亚洲综合影院| 五月玫瑰六月丁香| 国产精品不卡视频一区二区| 国产免费视频播放在线视频 | 国产麻豆成人av免费视频| 中文字幕av成人在线电影| 色视频www国产| 日韩欧美精品v在线| 在线免费观看不下载黄p国产| 街头女战士在线观看网站| 成人二区视频| 婷婷六月久久综合丁香| 国产亚洲av片在线观看秒播厂 | 国产在线一区二区三区精| 亚洲18禁久久av| 欧美一级a爱片免费观看看| 有码 亚洲区| 亚洲欧美成人精品一区二区| 亚洲精品国产成人久久av| 97超碰精品成人国产| 成人午夜高清在线视频| 亚洲电影在线观看av| 最近2019中文字幕mv第一页| 国产永久视频网站| 国产精品熟女久久久久浪| 精品国产一区二区三区久久久樱花 | 春色校园在线视频观看| 白带黄色成豆腐渣| 欧美成人a在线观看| 女人久久www免费人成看片| 一二三四中文在线观看免费高清| 麻豆成人av视频| 免费看美女性在线毛片视频| 亚洲精华国产精华液的使用体验| 日韩强制内射视频| 国产色婷婷99| 久久人人爽人人爽人人片va| 夜夜爽夜夜爽视频| 一级a做视频免费观看| 欧美一级a爱片免费观看看| 亚洲国产欧美在线一区| 天堂网av新在线| 777米奇影视久久| 晚上一个人看的免费电影| 国产真实伦视频高清在线观看| 人妻系列 视频| 国产成人a∨麻豆精品| 最近中文字幕2019免费版| 久久久成人免费电影| 精品久久久久久久久久久久久| 九色成人免费人妻av| 中文资源天堂在线| 久久韩国三级中文字幕| 国产免费福利视频在线观看| 欧美3d第一页| 美女国产视频在线观看| 国产激情偷乱视频一区二区| 69av精品久久久久久| 亚洲av日韩在线播放| 成年免费大片在线观看| 七月丁香在线播放| 日韩在线高清观看一区二区三区| 久久午夜福利片| 国产乱人视频| 午夜福利成人在线免费观看| 高清av免费在线| 国产探花极品一区二区| 亚洲av免费高清在线观看| 国产黄频视频在线观看| 在线观看人妻少妇| av在线天堂中文字幕| 亚洲av二区三区四区| 久99久视频精品免费| 国产一区二区在线观看日韩| 联通29元200g的流量卡| 最近中文字幕高清免费大全6| 高清午夜精品一区二区三区| 亚洲精品成人av观看孕妇| 看非洲黑人一级黄片| 亚洲人成网站高清观看| 午夜亚洲福利在线播放| 亚洲欧美一区二区三区黑人 | 国产一区二区在线观看日韩| 大又大粗又爽又黄少妇毛片口| 偷拍熟女少妇极品色| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美中文字幕日韩二区| 欧美bdsm另类| 国产伦精品一区二区三区四那| 丝袜美腿在线中文| 国产永久视频网站| 丝瓜视频免费看黄片| 亚洲内射少妇av| 国产成人精品福利久久| 熟女人妻精品中文字幕| 久久这里只有精品中国| 秋霞伦理黄片| 亚洲国产精品成人综合色| 麻豆久久精品国产亚洲av| 极品教师在线视频| 九色成人免费人妻av| 国产精品一区二区三区四区免费观看| 精品欧美国产一区二区三| 亚洲成人中文字幕在线播放| 联通29元200g的流量卡| 国产免费视频播放在线视频 | 欧美日韩亚洲高清精品| 九色成人免费人妻av| 少妇高潮的动态图| 狂野欧美白嫩少妇大欣赏| 九九久久精品国产亚洲av麻豆| 波多野结衣巨乳人妻| 国产精品久久视频播放| 亚洲真实伦在线观看| 最近最新中文字幕大全电影3| 免费大片18禁| 国产精品一区二区性色av| 亚洲av不卡在线观看| 国产成人精品婷婷| 成人高潮视频无遮挡免费网站| av女优亚洲男人天堂| 亚洲欧洲国产日韩| kizo精华| 毛片女人毛片| 少妇人妻一区二区三区视频| 亚洲三级黄色毛片| 国产成人a∨麻豆精品| 男人舔奶头视频| 国产亚洲av片在线观看秒播厂 | a级一级毛片免费在线观看| 精品久久久久久久末码| 亚洲精品,欧美精品| 国产成人精品一,二区| 成人午夜高清在线视频| 男女啪啪激烈高潮av片| 中文字幕人妻熟人妻熟丝袜美| 中国国产av一级| 天堂中文最新版在线下载 | 一级爰片在线观看| 国产av国产精品国产| 91午夜精品亚洲一区二区三区| 国产久久久一区二区三区| 舔av片在线| 国产一级毛片七仙女欲春2| 成人亚洲精品av一区二区| a级毛片免费高清观看在线播放| 精品国产一区二区三区久久久樱花 | 国产亚洲5aaaaa淫片| 亚洲av成人精品一区久久| 日本av手机在线免费观看| 精品熟女少妇av免费看| 又爽又黄无遮挡网站| 国产男人的电影天堂91| 国产白丝娇喘喷水9色精品| 麻豆精品久久久久久蜜桃| 97超视频在线观看视频| 欧美激情久久久久久爽电影| 国产成人aa在线观看| 亚洲欧美成人综合另类久久久| 国产精品不卡视频一区二区| 久久久色成人| 亚洲怡红院男人天堂| 麻豆成人av视频| 人妻少妇偷人精品九色| 亚洲av.av天堂| 日韩,欧美,国产一区二区三区| 麻豆成人午夜福利视频| 国产成人91sexporn| 高清视频免费观看一区二区 | 日韩av不卡免费在线播放| 在线观看人妻少妇| 日日摸夜夜添夜夜爱| 日产精品乱码卡一卡2卡三| 欧美一级a爱片免费观看看| 六月丁香七月|