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

    平板和靜葉表面氣流-水膜耦合流動(dòng)特性的數(shù)值研究

    2016-12-23 01:27:08范小軍李亮李森張翔
    關(guān)鍵詞:葉柵水膜剪切力

    范小軍,李亮,李森,張翔

    (1.西安交通大學(xué)能源與動(dòng)力工程學(xué)院,710049,西安;2.中國(guó)航發(fā)商用發(fā)動(dòng)機(jī)有限責(zé)任公司,200241,上海)

    ?

    平板和靜葉表面氣流-水膜耦合流動(dòng)特性的數(shù)值研究

    范小軍1,李亮1,李森2,張翔2

    (1.西安交通大學(xué)能源與動(dòng)力工程學(xué)院,710049,西安;2.中國(guó)航發(fā)商用發(fā)動(dòng)機(jī)有限責(zé)任公司,200241,上海)

    針對(duì)平板表面的空氣水膜和透平靜葉柵中的水蒸氣水膜耦合流動(dòng)特性提出了分析氣流與壁面水膜耦合作用的數(shù)值方法,即氣相主流和液相水膜視為相對(duì)獨(dú)立的開(kāi)口系,通過(guò)在兩相各自的控制方程中添加考慮相間動(dòng)量和能量交換的源項(xiàng),實(shí)現(xiàn)水膜和氣流的雙向耦合計(jì)算。研究表明:平板表面水膜厚度和進(jìn)口水膜雷諾數(shù)近似呈1/2指數(shù)冪關(guān)系,與來(lái)流馬赫數(shù)呈反比關(guān)系,與出口背壓呈弱負(fù)指數(shù)次冪關(guān)系;水膜流速與出口背壓近似呈線性關(guān)系;水膜附加損失對(duì)水膜流量的變化較為敏感,對(duì)來(lái)流馬赫數(shù)并不敏感;透平靜葉壓力面?zhèn)人ず穸确植驾^為均勻,在葉頂角區(qū)的葉片表面出現(xiàn)水膜聚積,吸力面?zhèn)人ず穸茸兓^為劇烈,在葉頂和葉根角區(qū)的葉片表面出現(xiàn)水膜聚積。從而得出影響透平葉柵通道中水膜分布的機(jī)理是,壁面曲率通過(guò)影響相間剪切力來(lái)影響壁面水膜厚度和速度的分布,端壁上從壓力面到吸力面的二次流引起水膜向吸力面?zhèn)染鄯e,通道渦和角渦引起的徑向壓力梯度導(dǎo)致葉片表面水膜沿葉高重新分布,水膜聚積對(duì)葉片表面的壓力分布和出口氣流角產(chǎn)生顯著影響。

    氣流;水膜;耦合流動(dòng)特性;數(shù)值研究

    葉輪機(jī)械中航空發(fā)動(dòng)機(jī)吸雨[1]、燃?xì)廨啓C(jī)濕壓縮[2]以及蒸汽透平中的水蒸氣發(fā)生非平衡凝結(jié)[3]都會(huì)形成不同尺寸的液態(tài)水滴,這些水滴在慣性力[4]和湍流擴(kuò)散作用[5]下沉積在葉輪機(jī)械的葉片和端壁表面且形成水膜。如圖1所示,水膜在相間剪切力、表面張力、固體表面摩擦力以及壓力梯度的作用下運(yùn)動(dòng),部分水膜在葉片尾緣破碎形成大水滴。對(duì)于壓氣機(jī),葉片表面水膜的存在顯著改變了氣流速度和落后角的徑向分布,使得級(jí)的氣動(dòng)效率降低,工作穩(wěn)定裕度下降[6];對(duì)于蒸汽透平,葉片和汽缸表面水膜的運(yùn)動(dòng)則引起葉片型面損失的增加和葉片的水蝕[3]。

    圖1 壁面水膜流動(dòng)示意圖

    水膜和氣流之間的耦合作用是引起葉輪機(jī)械氣動(dòng)性能變化的根本原因,為了深刻認(rèn)識(shí)水膜出現(xiàn)時(shí)葉輪機(jī)械性能變化的原因與規(guī)律,必須深入了解葉輪機(jī)械中水膜和高速氣流之間耦合作用的機(jī)理。由于葉輪機(jī)械中水膜的厚度在數(shù)十微米的量級(jí),目前水膜特性的實(shí)驗(yàn)測(cè)量還存在很大困難,因此數(shù)值計(jì)算便成為了研究這一問(wèn)題的有效手段。

    Kirillov將葉片簡(jiǎn)化為平板,最早建立了平板葉片表面水膜運(yùn)動(dòng)的方程,計(jì)算了不同安裝角下平板葉片表面水膜的運(yùn)動(dòng)軌跡[7]。Williams等將Kirillov的方法推廣到真實(shí)三維葉片表面水膜軌跡的計(jì)算中,在計(jì)算方法上取得突破[8]。但是,上述兩種方法都只能給出水膜的速度和軌跡,而不能給出水膜的厚度分布。Nikolaidis等基于三維計(jì)算流體動(dòng)力學(xué)(CFD)氣動(dòng)計(jì)算結(jié)果發(fā)展了壓氣機(jī)葉片表面水膜特性的計(jì)算程序[1],Fendler等對(duì)于透平發(fā)展了基于子午通流計(jì)算結(jié)果的水膜特性計(jì)算程序[9],這兩種方法能夠反映葉輪機(jī)械實(shí)際三維流動(dòng)條件下葉片表面水膜速度和厚度分布的更多細(xì)節(jié),但僅考慮了氣流對(duì)水膜的作用,仍然屬于單向耦合方法。

    本文發(fā)展了葉輪機(jī)械中水膜和氣流耦合作用的數(shù)值計(jì)算方法,對(duì)平板表面的空氣-水膜流動(dòng)以及葉柵通道中的水蒸氣-水膜流動(dòng)進(jìn)行了研究。

    1 數(shù)值方法

    假設(shè)氣相主流與壁面的液相水膜之間無(wú)相變引起的傳質(zhì)過(guò)程,僅存在相間的動(dòng)量和能量交換。將氣相主流和液相水膜視為相對(duì)獨(dú)立的開(kāi)口系,通過(guò)在兩相各自的控制方程中添加考慮相間動(dòng)量和能量交換的源項(xiàng),實(shí)現(xiàn)水膜和氣流的雙向耦合計(jì)算。

    1.1 流動(dòng)控制方程

    在葉輪機(jī)械中,氣流作用下壁面的水膜厚度遠(yuǎn)小于壁面的曲率半徑,因此可以不考慮水膜沿壁面法向的運(yùn)動(dòng)。這樣,水膜在葉片和端壁三維曲面表面進(jìn)行二維流動(dòng),流動(dòng)控制方程為

    (1)

    A=[h,hVf,hTf]T

    B=[hVf,hVfVf,hVfTf]T

    式中:h為水膜厚度;Vf為水膜速度;Tf為水膜平均溫度;s為壁面梯度算子;A為虛擬時(shí)間推進(jìn)項(xiàng);B為平均通量項(xiàng);Sf為源項(xiàng);K為內(nèi)部黏性項(xiàng)。Sf、K表達(dá)式分別為

    (2)

    (3)

    動(dòng)量源第一項(xiàng)-hsPL/ρf中的PL包含了氣相壓力項(xiàng)Pv、法向重力項(xiàng)Pg和表面張力項(xiàng)Pσ;動(dòng)量源第二項(xiàng)gτh代表平行于水膜流動(dòng)方向τ的重力分量;第三項(xiàng)3τfs/2ρf代表相界面處的氣流剪切力;內(nèi)部黏性項(xiàng)中νf為水的黏性系數(shù)。相應(yīng)表達(dá)式為

    PL=Pv+Ph+Pσ

    Pg=-ρfh(ng)

    (4)

    Pσ=-σs·(sh)

    式中:ρf為水的密度;g為重力加速度;n為法向矢量;σ為表面張力。

    能量方程源項(xiàng)中,Ts為相界面上的溫度,Tw為壁面溫度。水膜內(nèi)的溫度可以認(rèn)為沿壁面法向線性分布,在靠近壁面的h/2區(qū)域內(nèi)水膜溫度由Tw變?yōu)門(mén)f,在靠近相界面的h/2區(qū)域內(nèi)水膜溫度由Tf變?yōu)門(mén)s。

    氣相的控制方程采用Navier-Stokes方程。對(duì)于笛卡爾坐標(biāo)系下的三維黏性流動(dòng)有

    ·D=F+Sv

    (5)

    式中:p、ρv、U、Ev分別為氣相的壓力、密度、速度和總能;Π為應(yīng)力張量;F為黏性項(xiàng);Sv為氣相源項(xiàng)。方程(1)、(5)和狀態(tài)方程構(gòu)成了求解氣流與水膜雙向耦合流動(dòng)的全部控制方程。

    數(shù)值求解時(shí)借助商用軟件Fluent進(jìn)行,其中氣相流動(dòng)控制方程(5)由軟件直接求解,源項(xiàng)Sv通過(guò)軟件提供的用戶(hù)接口UDF添加。水膜流動(dòng)控制方程(1)通過(guò)用戶(hù)接口UDS求解。空間離散采用二階迎風(fēng)格式,氣相流動(dòng)求解采用k-ω湍流模型。

    1.2 數(shù)值方法驗(yàn)證

    Hammitt在水蒸氣風(fēng)洞中對(duì)平板表面的水膜厚度進(jìn)行了實(shí)驗(yàn)研究[10],實(shí)驗(yàn)段如圖2所示。平板前緣沿垂直流動(dòng)方向布置了寬度為1 mm的注水縫隙,流出注水縫隙的水在氣流剪切力作用下的平板表面形成水膜;水膜厚度由布置在平板中心線上的電導(dǎo)傳感器1進(jìn)行測(cè)量。

    圖2 平板表面水膜實(shí)驗(yàn)[10]

    按照實(shí)驗(yàn)段參數(shù),建立了用于數(shù)值模擬的矩形通道模型。通道流動(dòng)截面為30 mm×30 mm,軸向弦長(zhǎng)c=100 mm,通道底面前緣設(shè)置了寬度為1 mm的注水縫隙。保持來(lái)流水蒸氣的溫度T0v=325.45 K、注水溫度T0f=325.45 K以及出口背壓P1=20 kPa不變,改變進(jìn)口水蒸氣的流速和注水流量m對(duì)通道底面的水膜厚度進(jìn)行了計(jì)算,并與文獻(xiàn)[10]中實(shí)驗(yàn)數(shù)據(jù)擬合曲線進(jìn)行了比較,結(jié)果如圖3所示??梢钥吹?本文數(shù)值計(jì)算的平板表面水膜厚度在水蒸氣流速約從50~200 m/s、水膜流量從10~70 cm3/min變化的范圍內(nèi)與實(shí)驗(yàn)測(cè)量結(jié)果一致。

    圖3 平板表面水膜厚度分布

    2 平板表面的空氣-水膜耦合流動(dòng)

    研究平板表面的空氣-水膜耦合流動(dòng)是進(jìn)一步分析復(fù)雜葉柵通道內(nèi)水膜和氣流耦合流動(dòng)特性的基礎(chǔ)。本節(jié)采用與上節(jié)相同的幾何模型,對(duì)不同來(lái)流馬赫數(shù)、出口背壓和水膜流量條件下的平板表面空氣-水膜耦合流動(dòng)特性進(jìn)行了研究。

    2.1 氣流剪切力作用下的平板表面水膜流動(dòng)特性

    保持空氣進(jìn)口溫度T0a=300 K,注水溫度T0f=300 K,對(duì)來(lái)流馬赫數(shù)Ma為0.1~0.6、進(jìn)口水膜雷諾數(shù)Ref為33~188、出口背壓為0.05~0.25 MPa共48個(gè)工況條件下的平板水膜進(jìn)行了計(jì)算。水膜雷諾數(shù)Ref定義為

    Ref=ρfud/μf=Q/μf

    (6)

    式中:u為注水縫隙出口處液相的速度;d為注水縫隙的寬度;Q為水膜流量;μf為水的運(yùn)動(dòng)黏性系數(shù)。圖4給出了來(lái)流馬赫數(shù)Ma=0.4時(shí)不同水膜雷諾數(shù)下沿流向x/c(x為以進(jìn)氣邊為原點(diǎn)的軸向坐標(biāo))水膜厚度和相間剪切力Lfs的變化。

    圖4 平板表面水膜厚度和相間剪切力沿流向的分布

    在注水縫隙處,由于水膜剛剛形成,其流向速度近似為0,此時(shí)相界面上的速度梯度最大,水膜厚度和相間剪切力均達(dá)到最大值;離開(kāi)注水縫隙后,受氣流的剪切力作用,水膜流速增加而厚度減小,同時(shí)導(dǎo)致相界面上的速度梯度減小,相間剪切力隨之迅速減小;在遠(yuǎn)離注水縫隙的下游,氣流處于充分發(fā)展的湍流區(qū),相間剪切作用力逐漸下降并趨于穩(wěn)定,因而水膜厚度緩慢增加并趨于穩(wěn)定。隨著進(jìn)口水膜雷諾數(shù)的增加,水膜厚度依次增加,但分布趨勢(shì)類(lèi)似;相間剪切力的大小和分布基本不變。這表明相間剪切力主要受氣流馬赫數(shù)的影響,與水膜厚度基本無(wú)關(guān)。

    為了分析水膜的厚度和平均速度與各影響因素之間的關(guān)系,圖5給出了x/c=0.8位置處水膜參數(shù)隨邊界條件的變化趨勢(shì)。在給定出口背壓P和水膜進(jìn)口雷諾數(shù)時(shí),水膜厚度和來(lái)流馬赫數(shù)近似呈反比關(guān)系(見(jiàn)圖5a),而水膜平均流速Vr與來(lái)流馬赫數(shù)呈線性關(guān)系(見(jiàn)圖5d)。在給定來(lái)流馬赫數(shù)和出口背壓時(shí),水膜厚度和水膜平均流速與水膜進(jìn)口雷諾數(shù)近似呈1/2指數(shù)冪關(guān)系(見(jiàn)圖5b、5e)。另外,

    間剪切力與進(jìn)口水膜雷諾數(shù)無(wú)關(guān),而僅隨來(lái)流馬赫數(shù)變化(見(jiàn)圖5d)。在給定馬赫數(shù)和水膜雷諾數(shù)時(shí),水膜厚度隨背壓的增加而減小,水膜速度則呈現(xiàn)增加的趨勢(shì);背壓通過(guò)影響氣水分界面上的溫度來(lái)影響水膜的黏性系數(shù)μf,最終影響水膜的厚度和速度分布。

    采用線性回歸方法對(duì)Ref=20~200范圍內(nèi)水膜厚度與馬赫數(shù)、水膜雷諾數(shù)的關(guān)系進(jìn)行了擬合,結(jié)果為

    (7)

    (8)

    2.2 水膜對(duì)氣流總壓損失特性的影響

    邊界層損失是內(nèi)流損失的主要組成部分之一。當(dāng)壁面存在水膜時(shí),由于湍流邊界層的厚度遠(yuǎn)遠(yuǎn)大于當(dāng)?shù)厮さ暮穸?所以水膜對(duì)主流的影響可以忽略,但其顯著影響邊界層內(nèi)的流動(dòng)。如圖6所示,當(dāng)壁面存在水膜時(shí),氣水分界面上水膜頂層受氣流的剪切力作用發(fā)生運(yùn)動(dòng),同時(shí)消耗氣流的部分動(dòng)能,因而可能影響氣流的總壓損失特性。

    為了定量描述水膜對(duì)總壓損失特性的影響,定義有、無(wú)水膜時(shí)通道總壓損失系數(shù)比

    (a)水膜厚度隨Ma的變化趨勢(shì) (b)水膜厚度隨Ref的變化趨勢(shì) (c)水膜厚度隨出口背壓的變化趨勢(shì)

    (d)水膜平均流速隨Ma的變化趨勢(shì) (e)水膜平均流速隨Ref的變化趨勢(shì) (f)水膜平均流速隨出口背壓的變化趨勢(shì)圖5 x/c=0.8位置處水膜參數(shù)隨邊界條件的變化

    圖6 有、無(wú)水膜時(shí)的邊界層流動(dòng)

    (9)

    圖7 總壓損失特性隨馬赫數(shù)的變化規(guī)律

    3 透平葉柵中的水蒸氣-水膜耦合流動(dòng)

    對(duì)某核電汽輪機(jī)末級(jí)靜葉柵中水蒸氣-水膜耦合流動(dòng)進(jìn)行研究。計(jì)算時(shí)采用結(jié)構(gòu)化網(wǎng)格,經(jīng)網(wǎng)格無(wú)關(guān)性驗(yàn)證確定網(wǎng)格總數(shù)為230萬(wàn)。邊界條件從葉柵所在的低壓透平全三維計(jì)算結(jié)果[12]中獲得。進(jìn)口蒸汽流量為3.96 kg/s,總溫為337.971 K,進(jìn)口湍流度為5%,出口壓力為16 618.2 Pa。葉片和上端壁表面水滴沉積形成水膜的流量分別為0.045 kg/(m2·s)和0.019 kg/(m2·s);葉片和上端壁表面水膜的初始平均厚度為45.7 μm;假設(shè)水滴撞擊壁面形成的水膜保留原始水滴速度的20%[1],則水膜的初始速度為47.62 m/s。

    3.1 葉片和上端壁表面的水膜分布

    圖8給出了主流水蒸氣作用下葉片表面水膜的分布。壓力面水膜厚度較為均勻,沿葉高方向大部分區(qū)域水膜厚度維持在25~35 μm范圍內(nèi);在葉頂與端壁交匯的角區(qū),壓力面出現(xiàn)了帶狀的水膜聚積區(qū),其厚度達(dá)到90 μm左右。與壓力面相比,吸力面水膜厚度的變化較為劇烈:在20%軸向弦長(zhǎng)之前,吸力面水膜的厚度小于20 μm,且沿葉高方向基本不變,沿流動(dòng)方向略有增加;在20%軸向弦長(zhǎng)之后,水膜厚度迅速增大,并在靠近葉頂和葉根的區(qū)域形成了兩個(gè)水膜聚集區(qū),這兩個(gè)區(qū)域水膜的厚度均達(dá)到了90 μm以上。

    圖8 葉片表面水膜厚度的分布

    葉片表面水膜厚度的分布與相間剪切力、水膜黏性力、表面曲率、表面張力和壁面剪切力作用下葉片表面的水膜流動(dòng)行為有關(guān)。圖9給出了葉片表面的水膜流速分布。與水膜厚度分布相對(duì)應(yīng),壓力面水膜流速在80%葉高以下沿徑向基本均勻,但在80%葉高以上沿徑向出現(xiàn)波動(dòng)。在吸力面,葉根和葉頂兩個(gè)角區(qū)水膜出現(xiàn)低速區(qū),其他區(qū)域的水膜速度梯度與流向大體相向。

    圖9 葉片表面水膜速度的分布

    葉片表面水膜厚度和速度的上述特征與葉柵通道內(nèi)的二次流動(dòng)密切相關(guān)。圖10給出了葉柵通道內(nèi)及尾跡區(qū)二次渦的發(fā)展過(guò)程,其中c為軸向弦長(zhǎng),x為以進(jìn)氣邊為原點(diǎn)的軸向坐標(biāo)。結(jié)合圖8~圖10可以看到,在吸力面與上端壁形成的角區(qū)附近,主流通道渦和角渦所誘導(dǎo)的徑向速度通過(guò)相間剪切力帶動(dòng)葉頂區(qū)域的水膜向葉展中部遷移,形成了吸力面頂部區(qū)域的水膜分布特征。造成吸力面葉根角區(qū)以及壓力面水膜沿葉高分布特征的原因與此類(lèi)似。

    圖10 葉柵通道中二次渦的發(fā)展過(guò)程

    圖11給出了上端壁的水膜厚度和速度的分布。圖11a中:在葉片前緣?、?處,子午型線突然擴(kuò)張使得上端壁邊界層區(qū)域的氣流產(chǎn)生分離,形成回流區(qū),上端壁表面的水膜在氣流的剪切作用下產(chǎn)生堆積;由?和?區(qū)域水膜的流線分布可以觀察到,上端壁處水膜從壓力面向吸力面遷移,這使得上端壁表面靠近吸力面的水膜不斷增厚,形成水膜的帶狀聚積。根據(jù)連續(xù)性方程,帶狀聚積區(qū)水膜流量的增大也使得此區(qū)域的水膜流速明顯增加,如圖11b所示。

    (a)水膜厚度(b)水膜速度 圖11 上端壁表面水膜特性參數(shù)的分布

    在上端壁靠近葉片尾緣處,來(lái)自吸力面和壓力面的主氣流以及水膜均各自交匯并相互影響。主氣流流經(jīng)葉片尾緣時(shí)形成低速的尾跡區(qū),導(dǎo)致上端壁圖11a中 ?區(qū)的水膜速度比較低;在尾緣下游壓力面和吸力面氣流的壓差作用下,上端壁表面尾緣下游的水膜向吸力面偏移。圖11a中 ?區(qū)和 ?區(qū)的水膜聚積是主流通道渦在尾跡區(qū)遷移引起的,如圖12所示。通道渦在葉柵出口下游截面得到擴(kuò)展,并導(dǎo)致上端壁表面的水膜厚度重新分布。

    (a)x/c=1.55 (b)x/c=1.6 (c)x/c=1.8圖12 葉柵下游二次渦的發(fā)展過(guò)程

    總的來(lái)看,影響水膜分布的主要因素有3點(diǎn):一是壁面曲率,它通過(guò)影響壁面壓力分布改變了主流場(chǎng)的速度分布,速度分布又進(jìn)一步影響相間的剪切力分布,從而影響表面水膜的速度和厚度分布;二是端壁表面附近橫向二次流對(duì)水膜的遷移作用,引起端壁表面水膜向葉片吸力面?zhèn)染鄯e;三是通道渦和角渦的影響,引起水膜沿葉高重新分布。

    3.2 水膜對(duì)葉柵流場(chǎng)的影響

    葉片和端壁表面的水膜與主流相互作用,對(duì)葉柵的氣動(dòng)特性產(chǎn)生了影響。圖13給出了85%~95%葉高處葉片表面壓力的分布??梢钥吹?在吸力面靠近尾緣區(qū)域,沿相對(duì)葉高方向水膜引起的表面壓力差異越來(lái)越明顯,在95%葉高處表面壓力變化幅值達(dá)到±300 Pa,約為該截面進(jìn)、出口壓降的13%。該位置對(duì)應(yīng)圖8中吸力面水膜聚積的區(qū)域。可以看到,吸力面靠近葉頂區(qū)域的水膜聚積是由角區(qū)氣流的二次流引起的,而水膜的聚積又反過(guò)來(lái)對(duì)葉柵中的流動(dòng)產(chǎn)生了影響。

    圖13 85%~95%葉高處葉片表面壓力分布

    水膜的存在也使得葉片的出口氣流角發(fā)生了變化。如圖14所示,與無(wú)水膜時(shí)相比,進(jìn)口氣流角相同,而水膜使葉柵的出口氣流角有所增加,在約85%葉高處增加了0.6°。

    圖14 進(jìn)、出口氣流角沿葉高的分布

    對(duì)于本文分析的透平靜葉,由于水膜引起的氣流特性變化僅存在于局部區(qū)域,因此對(duì)葉柵的氣動(dòng)特性沒(méi)有實(shí)質(zhì)性改變。但是,考慮到透平葉片較厚,而水膜的厚度僅在數(shù)十微米量級(jí),因此水膜引起葉片表面壓力和出口氣流角的變化相當(dāng)顯著。對(duì)于航空發(fā)動(dòng)機(jī)吸雨的情況,壓氣機(jī)葉片要薄得多,水膜對(duì)葉柵氣動(dòng)特性可能產(chǎn)生更為顯著的影響。

    4 結(jié) 論

    本文對(duì)平板表面和透平靜葉中氣流和水膜的耦合作用進(jìn)行了研究,主要結(jié)論如下。

    (1)平板表面水膜厚度與進(jìn)口水膜雷諾數(shù)呈1/2指數(shù)冪關(guān)系,與來(lái)流馬赫數(shù)呈反比關(guān)系,與出口背壓呈弱負(fù)指數(shù)次冪關(guān)系。水膜平均流速隨來(lái)流馬赫數(shù)的增加線性增大,與進(jìn)口水膜雷諾數(shù)呈1/2指數(shù)冪關(guān)系,并隨出口背壓的增加逐步增大,但增速逐漸降低。直通道總壓損失系數(shù)比隨來(lái)流馬赫數(shù)、水膜流量的增大而增大,水膜附加損失對(duì)水膜流量和出口背壓的變化較為敏感。

    (2)葉柵通道中影響壁面水膜分布的機(jī)理為:壁面曲率通過(guò)影響相間剪切力來(lái)影響水膜的厚度和速度;端壁上從壓力面到吸力面的橫向二次流引起水膜在吸力面附近聚積;通道渦和角渦引起水膜沿葉高重新分布。葉柵通道中水膜對(duì)氣流產(chǎn)生影響。沿葉高水膜聚集區(qū)的葉片表面壓力分布發(fā)生顯著變化,葉柵出口氣流角增加。

    [1] NIKOLAIDIS T, PILIDIS P. The effect of water ingestion on an axial flow compressor performance [J]. Journal of Aerospace Engineering, 2014, 228(3): 411-423.

    [2] BHARGAVE R, MEHERHOMJI C, CHAKER C, et al. Gas turbine fogging technology: a state-of-the-art review part I Inlet evaporative fogging analytical and experimental aspects [J]. Journal of Engineering for Gas Turbines and Power, 2007, 129(2): 443-453.

    [3] MOORE M J, SIEVERDING C H. Aerothermodynamic of low pressure steam turbines and condensers [M]. Carlsbad, USA: Hemisphere Publishing Corporation, 1986: 1-190.

    [4] YOUNG J B, YAU K K. The inertial deposition of fog droplets on steam turbine blades [J]. ASME Journal of Turbomachinery, 1988, 110(2): 155-162.

    [5] YAU K K, YOUNG J B. The deposition of fog droplets on steam turbine blades by turbulent diffusion [J]. ASME Journal of Turbomachinery, 1987, 109(3): 429-435.

    [6] 劉波, 曹志鵬, 高嵩, 等. 來(lái)流含水對(duì)航空發(fā)動(dòng)機(jī)風(fēng)扇/壓氣機(jī)特性的影響 [J]. 航空動(dòng)力學(xué)報(bào), 2006, 20(6): 1041-1047. LIU Bo, CAO Zhipeng, GAO Song, et al. Influence of inlet water ingestion on aero-engine fan-compressor performance [J]. Journal of Aerospace Power, 2006, 20(6): 1041-1047.

    [7] KIRILLOV I I, IABLONIK R M. Fundamentals of the theory of turbines operating on wet steam [M]. Washington, DC, USA: NASA, 1970: 1-36.

    [8] WILLIAMS J, YOUNG J B. Movement of deposited water on turbomachinery rotor blade surfaces [J]. ASME Journal of Turbomachinery, 2007, 129(2): 394-403.

    [9] FENDLER Y, DOREY J M, STANCIU M, et al. Developments for modeling of droplets deposition and liquid film flow in a throughflow code for steam turbines: GT2012-68968 [R]. New York, USA: ASME, 2012.

    [10]HAMMITT F G. Liquid film and droplet stability consideration as applied to wet steam flow [J]. Forschung im Ingenieurwesen: A, 1981, 47(1): 1-14.

    [11]SCHLICHTING H, GERSTEN K. Boundary-layer theory [M]. Heidelberg, Berlin: Springer-Verlag, 2003: 497-518.

    [12]LI Liang, YANG Jiandao, YOU Wei, et al. Investigation of the vapour-liquid two-phase flow in the low-pressure cylinder of a 1 000 MW nuclear power steam turbine [J]. Proc IMechE: Part A Journal of Power and Energy, 2014, 228(2): 178-185.

    (編輯 苗凌)

    Numerical Investigation for Coupled Flow Behavior Between Gas Flow and Wall Film of Flat and Turbine Stator

    FAN Xiaojun1,LI Liang1,LI Sen2,ZHANG Xiang2

    (1. School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China; 2. AECC Commercial Aircraft Engine Co. Ltd., Shanghai 200241, China)

    A numerical algorithm was employed to investigate the coupled flow behavior between gas flow and wall film. The air flow and wall film on a flat, and the water vapor flow and wall film in a turbine cascade were focused on. It is concluded that the film thickness is approximately proportional to the root of the inlet wall film Reynolds number, and inversely proportional to the mainstream flow Mach number. Film thickness and the outlet pressure have a weak negative exponential order relation, and the interface velocity is linearly related to the outlet pressure. The additional losses due to water film are sensitive to the change of water flow rate, but not sensitive to the flow Mach number. As for the turbine stator, the pressure surface film thickness is relative even, except that a film accumulation region appears at the tip corner region. In contrast, the suction side thickness changes obviously, and the film accumulation is observed at both the tip and root corner regions. The wall film distribution in a turbine cascade is affected by three factors. The wall curvature influences the film thickness and velocity by influencing the interphase shear force. The endwall secondary flow from the pressure side to the suction side leads to film accumulation near the suction surface. The radial pressure gradient caused by passage vortex and corner vortex results in film redistribution along the blade height. On the other hand, the film accumulation on the blade exerts a significant impact on the surface pressure and outlet flow angle.

    gas flow; wall film; coupling flow behavior; numerical investigation

    2016-05-23。 作者簡(jiǎn)介:范小軍(1992—),男,碩士生;李亮(通信作者),男,副教授,博士生導(dǎo)師。

    時(shí)間:2016-09-08

    10.7652/xjtuxb201611002

    TK263.3

    A

    0253-987X(2016)11-0007-07

    網(wǎng)絡(luò)出版地址:http:∥www.cnki.net/kcms/detail/61.1069.T.20160908.1101.002.html

    猜你喜歡
    葉柵水膜剪切力
    一種汽車(chē)發(fā)動(dòng)機(jī)用橡膠減震器
    巧測(cè)水膜張力
    少兒科技(2022年4期)2022-04-14 23:48:10
    基于Fluent的不同彎道剪切力分布特性數(shù)值研究
    亞聲速壓氣機(jī)平面葉柵及其改型的吹風(fēng)試驗(yàn)
    濕滑跑道飛機(jī)著陸輪胎-水膜-道面相互作用
    水流剪切力對(duì)供水管道管壁生物膜生長(zhǎng)的影響
    非能動(dòng)核電站安全殼外壁下降水膜的穩(wěn)定性分析
    水潤(rùn)滑軸承水膜壓力無(wú)線測(cè)試系統(tǒng)研究
    超音速透平葉柵造型設(shè)計(jì)及分析
    針對(duì)軸流壓氣機(jī)的非軸對(duì)稱(chēng)端壁造型優(yōu)化設(shè)計(jì)
    精品国内亚洲2022精品成人 | 丁香六月天网| 91精品三级在线观看| 1024视频免费在线观看| 老司机午夜福利在线观看视频 | 1024香蕉在线观看| 一级片免费观看大全| 久久久久久亚洲精品国产蜜桃av| 亚洲伊人色综图| 亚洲第一青青草原| 大香蕉久久网| 国产熟女午夜一区二区三区| 国产一区二区 视频在线| 在线十欧美十亚洲十日本专区| 精品人妻1区二区| 中文字幕人妻丝袜一区二区| 免费人妻精品一区二区三区视频| 多毛熟女@视频| 欧美精品啪啪一区二区三区| 91成人精品电影| 日韩熟女老妇一区二区性免费视频| 久久精品亚洲av国产电影网| 国产精品一区二区在线不卡| av免费在线观看网站| 欧美在线黄色| 亚洲欧美一区二区三区黑人| 两人在一起打扑克的视频| 国产1区2区3区精品| 999精品在线视频| 久久中文看片网| 一边摸一边抽搐一进一小说 | 久热爱精品视频在线9| 精品国产一区二区三区四区第35| 性色av乱码一区二区三区2| 99久久精品国产亚洲精品| 青青草视频在线视频观看| 一边摸一边做爽爽视频免费| 少妇被粗大的猛进出69影院| 久久精品91无色码中文字幕| 午夜成年电影在线免费观看| 成年人午夜在线观看视频| 亚洲欧美色中文字幕在线| 丝袜美腿诱惑在线| 一级片'在线观看视频| 99国产精品一区二区三区| 精品亚洲成a人片在线观看| 午夜福利在线免费观看网站| 日韩有码中文字幕| 亚洲成人手机| 亚洲国产毛片av蜜桃av| 纵有疾风起免费观看全集完整版| 久久久国产一区二区| 老司机深夜福利视频在线观看| 亚洲美女黄片视频| 欧美精品人与动牲交sv欧美| 亚洲国产中文字幕在线视频| 女性被躁到高潮视频| 在线观看一区二区三区激情| 久久人妻熟女aⅴ| 久久久久视频综合| 国产精品久久久久久精品电影小说| 国产又爽黄色视频| 午夜福利在线观看吧| svipshipincom国产片| 欧美精品亚洲一区二区| 91麻豆精品激情在线观看国产 | 亚洲国产精品一区二区三区在线| 精品久久蜜臀av无| 在线观看66精品国产| 国产不卡av网站在线观看| 国产成人影院久久av| 在线播放国产精品三级| 午夜视频精品福利| 久久中文字幕一级| 成人18禁在线播放| 欧美日韩亚洲高清精品| 少妇裸体淫交视频免费看高清 | 天天操日日干夜夜撸| 免费久久久久久久精品成人欧美视频| 亚洲精品成人av观看孕妇| 久久天堂一区二区三区四区| 超色免费av| 两个人看的免费小视频| 女人爽到高潮嗷嗷叫在线视频| 99久久99久久久精品蜜桃| 中亚洲国语对白在线视频| 精品福利永久在线观看| 中国美女看黄片| 中亚洲国语对白在线视频| 如日韩欧美国产精品一区二区三区| 免费观看人在逋| 黄色毛片三级朝国网站| 色婷婷av一区二区三区视频| 99香蕉大伊视频| 国产精品偷伦视频观看了| 夜夜夜夜夜久久久久| 国产高清videossex| 亚洲av片天天在线观看| 黄网站色视频无遮挡免费观看| 久久精品国产a三级三级三级| 国产一卡二卡三卡精品| 欧美激情高清一区二区三区| 如日韩欧美国产精品一区二区三区| 精品午夜福利视频在线观看一区 | 黄色视频在线播放观看不卡| 国产在线一区二区三区精| 人人妻人人爽人人添夜夜欢视频| 51午夜福利影视在线观看| av网站在线播放免费| 搡老熟女国产l中国老女人| 日本五十路高清| 黑丝袜美女国产一区| 亚洲成国产人片在线观看| 99国产精品一区二区三区| 午夜福利乱码中文字幕| 桃红色精品国产亚洲av| 亚洲成人手机| 国产精品久久久人人做人人爽| 免费人妻精品一区二区三区视频| 日日夜夜操网爽| 纯流量卡能插随身wifi吗| www.999成人在线观看| 黄色片一级片一级黄色片| 国产成人免费无遮挡视频| 97在线人人人人妻| 日韩大码丰满熟妇| 丰满少妇做爰视频| 一进一出好大好爽视频| 99久久国产精品久久久| 精品人妻熟女毛片av久久网站| 久久久国产成人免费| 精品少妇内射三级| h视频一区二区三区| 日韩大码丰满熟妇| 香蕉久久夜色| 亚洲第一青青草原| 菩萨蛮人人尽说江南好唐韦庄| 一本综合久久免费| 欧美老熟妇乱子伦牲交| 亚洲伊人久久精品综合| 九色亚洲精品在线播放| 啦啦啦 在线观看视频| 免费不卡黄色视频| 国产精品熟女久久久久浪| 一区二区日韩欧美中文字幕| 国产av又大| 王馨瑶露胸无遮挡在线观看| av在线播放免费不卡| 黄片大片在线免费观看| 搡老熟女国产l中国老女人| 国产精品香港三级国产av潘金莲| 一级,二级,三级黄色视频| 久久影院123| 老熟妇仑乱视频hdxx| 俄罗斯特黄特色一大片| 免费av中文字幕在线| 不卡一级毛片| 中文字幕人妻丝袜制服| 成人18禁在线播放| 成年人午夜在线观看视频| tocl精华| 久久亚洲真实| 桃红色精品国产亚洲av| 黑人猛操日本美女一级片| 亚洲欧美一区二区三区久久| 91九色精品人成在线观看| 国产精品一区二区精品视频观看| 精品亚洲成a人片在线观看| 亚洲中文av在线| 亚洲人成电影免费在线| 啦啦啦在线免费观看视频4| 午夜视频精品福利| videosex国产| 欧美日韩一级在线毛片| 黄色怎么调成土黄色| 亚洲国产欧美日韩在线播放| 久久久精品免费免费高清| 亚洲男人天堂网一区| 日韩三级视频一区二区三区| 久久天躁狠狠躁夜夜2o2o| 在线av久久热| 午夜福利视频精品| 亚洲欧美日韩另类电影网站| 欧美在线黄色| 一夜夜www| 十分钟在线观看高清视频www| 在线观看免费高清a一片| 精品国产超薄肉色丝袜足j| 欧美变态另类bdsm刘玥| 满18在线观看网站| 国产激情久久老熟女| 欧美日韩国产mv在线观看视频| 国产av一区二区精品久久| 午夜精品久久久久久毛片777| 国产亚洲精品一区二区www | 国产三级黄色录像| 久久ye,这里只有精品| 自线自在国产av| 欧美精品一区二区大全| 大片免费播放器 马上看| 日韩欧美一区视频在线观看| av福利片在线| 又紧又爽又黄一区二区| 国产有黄有色有爽视频| 国产精品电影一区二区三区 | 国内毛片毛片毛片毛片毛片| 中文字幕最新亚洲高清| 亚洲第一av免费看| 三级毛片av免费| 9热在线视频观看99| 亚洲精品国产精品久久久不卡| 精品午夜福利视频在线观看一区 | 久久久久久久精品吃奶| 夜夜夜夜夜久久久久| 在线 av 中文字幕| 人人妻人人爽人人添夜夜欢视频| 桃红色精品国产亚洲av| 色综合婷婷激情| 在线永久观看黄色视频| aaaaa片日本免费| 国产视频一区二区在线看| 国产精品自产拍在线观看55亚洲 | 丰满饥渴人妻一区二区三| 大陆偷拍与自拍| 国产精品一区二区在线不卡| 久久久国产欧美日韩av| 欧美日韩一级在线毛片| 建设人人有责人人尽责人人享有的| 99国产综合亚洲精品| 精品第一国产精品| 亚洲第一av免费看| 日韩一卡2卡3卡4卡2021年| 国产日韩欧美视频二区| 亚洲一码二码三码区别大吗| 在线观看免费高清a一片| 久久久精品免费免费高清| 精品国产一区二区久久| 丰满饥渴人妻一区二区三| 亚洲av第一区精品v没综合| 黄片播放在线免费| 午夜福利一区二区在线看| 精品熟女少妇八av免费久了| 水蜜桃什么品种好| 欧美激情高清一区二区三区| 久久青草综合色| 巨乳人妻的诱惑在线观看| 香蕉国产在线看| 五月天丁香电影| 五月开心婷婷网| 国产片内射在线| 免费看a级黄色片| 在线观看免费日韩欧美大片| 在线 av 中文字幕| 国产伦人伦偷精品视频| 精品国产乱子伦一区二区三区| 欧美在线黄色| 首页视频小说图片口味搜索| 女人被躁到高潮嗷嗷叫费观| 色精品久久人妻99蜜桃| 99国产精品一区二区三区| 啦啦啦 在线观看视频| 亚洲国产中文字幕在线视频| 男女午夜视频在线观看| 黄色片一级片一级黄色片| 一边摸一边做爽爽视频免费| 51午夜福利影视在线观看| 国产1区2区3区精品| 蜜桃在线观看..| 一本—道久久a久久精品蜜桃钙片| 国产成人欧美| 大香蕉久久成人网| 曰老女人黄片| 男女免费视频国产| 99九九在线精品视频| 天堂动漫精品| 亚洲精品av麻豆狂野| 精品视频人人做人人爽| 99久久精品国产亚洲精品| 国产无遮挡羞羞视频在线观看| 在线看a的网站| 午夜免费鲁丝| 极品少妇高潮喷水抽搐| 亚洲九九香蕉| 亚洲精品国产一区二区精华液| 欧美日韩av久久| 1024香蕉在线观看| 男女边摸边吃奶| 久久人人爽av亚洲精品天堂| 亚洲 国产 在线| 日韩视频一区二区在线观看| 欧美一级毛片孕妇| 动漫黄色视频在线观看| 老司机影院毛片| 一本色道久久久久久精品综合| 不卡av一区二区三区| 波多野结衣av一区二区av| 51午夜福利影视在线观看| 亚洲国产中文字幕在线视频| 9热在线视频观看99| 久久毛片免费看一区二区三区| 在线观看免费日韩欧美大片| 黄色毛片三级朝国网站| 免费看a级黄色片| 中文字幕人妻丝袜制服| 黄色丝袜av网址大全| 一区二区三区激情视频| 岛国毛片在线播放| 可以免费在线观看a视频的电影网站| 欧美激情久久久久久爽电影 | 精品免费久久久久久久清纯 | 在线av久久热| 久热爱精品视频在线9| 亚洲av美国av| 一级片'在线观看视频| 十八禁高潮呻吟视频| 久久香蕉激情| 成人特级黄色片久久久久久久 | 亚洲av片天天在线观看| 女性生殖器流出的白浆| 亚洲色图综合在线观看| 老汉色av国产亚洲站长工具| 久久午夜综合久久蜜桃| 91麻豆精品激情在线观看国产 | 无人区码免费观看不卡 | 一级片免费观看大全| 免费一级毛片在线播放高清视频 | 丝袜人妻中文字幕| 国产极品粉嫩免费观看在线| 欧美大码av| 亚洲欧美激情在线| 大香蕉久久网| 一本—道久久a久久精品蜜桃钙片| 精品少妇一区二区三区视频日本电影| 天堂动漫精品| 亚洲免费av在线视频| 亚洲午夜精品一区,二区,三区| 人人妻人人爽人人添夜夜欢视频| 国产精品秋霞免费鲁丝片| 高清欧美精品videossex| 日韩欧美一区二区三区在线观看 | 露出奶头的视频| 国产精品欧美亚洲77777| 午夜视频精品福利| 亚洲精品中文字幕一二三四区 | 亚洲男人天堂网一区| 老司机午夜十八禁免费视频| 国产精品一区二区精品视频观看| 国产一区二区激情短视频| 久久久久国产一级毛片高清牌| 国产视频一区二区在线看| 亚洲精品国产精品久久久不卡| 国产97色在线日韩免费| 国产亚洲一区二区精品| 男女下面插进去视频免费观看| 成人永久免费在线观看视频 | 久久青草综合色| 午夜福利视频在线观看免费| 亚洲国产av新网站| 亚洲欧美日韩另类电影网站| 亚洲精品一二三| 欧美+亚洲+日韩+国产| 亚洲人成77777在线视频| 最新在线观看一区二区三区| 一边摸一边抽搐一进一出视频| 国产91精品成人一区二区三区 | 激情在线观看视频在线高清 | 亚洲熟妇熟女久久| 久久久久国产一级毛片高清牌| 精品一区二区三区视频在线观看免费 | 久久天堂一区二区三区四区| 国产一区二区在线观看av| 久热爱精品视频在线9| 一个人免费在线观看的高清视频| 黑人猛操日本美女一级片| 国产深夜福利视频在线观看| 男女边摸边吃奶| 国产在视频线精品| 国产精品二区激情视频| 精品国内亚洲2022精品成人 | 国产成人精品在线电影| 搡老熟女国产l中国老女人| 老司机影院毛片| 久久精品国产综合久久久| 久久国产精品大桥未久av| 精品第一国产精品| 精品乱码久久久久久99久播| 精品一区二区三区av网在线观看 | 国产男女内射视频| 久久毛片免费看一区二区三区| 伦理电影免费视频| 亚洲男人天堂网一区| 精品一区二区三区av网在线观看 | 色综合欧美亚洲国产小说| 精品一区二区三卡| 国产精品 欧美亚洲| 午夜视频精品福利| 麻豆av在线久日| 啦啦啦视频在线资源免费观看| 午夜免费鲁丝| 精品视频人人做人人爽| 女人精品久久久久毛片| 飞空精品影院首页| 国产一区二区激情短视频| 18禁黄网站禁片午夜丰满| 男男h啪啪无遮挡| 国产精品影院久久| 国产欧美亚洲国产| 99国产精品一区二区蜜桃av | 久久天躁狠狠躁夜夜2o2o| 50天的宝宝边吃奶边哭怎么回事| 成人国产一区最新在线观看| 天天操日日干夜夜撸| 搡老乐熟女国产| 成人国产一区最新在线观看| 狠狠婷婷综合久久久久久88av| 免费av中文字幕在线| 妹子高潮喷水视频| 成在线人永久免费视频| 少妇的丰满在线观看| 麻豆成人av在线观看| 午夜福利免费观看在线| 婷婷丁香在线五月| 亚洲 国产 在线| bbb黄色大片| 亚洲国产成人一精品久久久| 麻豆国产av国片精品| 动漫黄色视频在线观看| 一级,二级,三级黄色视频| 99国产精品一区二区三区| 90打野战视频偷拍视频| 黄色a级毛片大全视频| 人妻 亚洲 视频| 91麻豆av在线| 国产成人系列免费观看| 18禁裸乳无遮挡动漫免费视频| 悠悠久久av| 日韩大码丰满熟妇| 美女国产高潮福利片在线看| 岛国毛片在线播放| 色尼玛亚洲综合影院| 国产主播在线观看一区二区| 最近最新中文字幕大全电影3 | 亚洲九九香蕉| 老熟妇乱子伦视频在线观看| 十八禁人妻一区二区| 亚洲精品一二三| 久久久水蜜桃国产精品网| 伦理电影免费视频| 99re在线观看精品视频| 欧美黄色淫秽网站| 乱人伦中国视频| 午夜精品久久久久久毛片777| 亚洲中文字幕日韩| 视频在线观看一区二区三区| 手机成人av网站| 国产免费福利视频在线观看| 老司机深夜福利视频在线观看| 国产伦人伦偷精品视频| 不卡av一区二区三区| 99国产精品一区二区蜜桃av | 不卡av一区二区三区| 18禁美女被吸乳视频| 一本综合久久免费| 亚洲av第一区精品v没综合| 一区二区三区乱码不卡18| 99在线人妻在线中文字幕 | www.精华液| 99国产极品粉嫩在线观看| 菩萨蛮人人尽说江南好唐韦庄| av线在线观看网站| 国产欧美日韩一区二区精品| 亚洲精品中文字幕在线视频| 两性夫妻黄色片| 女人被躁到高潮嗷嗷叫费观| 欧美日韩一级在线毛片| 久久人人97超碰香蕉20202| 国产极品粉嫩免费观看在线| www.自偷自拍.com| 久久这里只有精品19| 极品少妇高潮喷水抽搐| 老司机午夜十八禁免费视频| 成人国产一区最新在线观看| 一进一出抽搐动态| 久久久久久人人人人人| 日韩有码中文字幕| 婷婷成人精品国产| 成人黄色视频免费在线看| 久久国产精品影院| 久久精品成人免费网站| 水蜜桃什么品种好| 亚洲黑人精品在线| 精品福利观看| cao死你这个sao货| 午夜激情av网站| 欧美精品人与动牲交sv欧美| 老司机在亚洲福利影院| 亚洲一码二码三码区别大吗| 亚洲va日本ⅴa欧美va伊人久久| 极品教师在线免费播放| 午夜福利视频精品| 性高湖久久久久久久久免费观看| 免费一级毛片在线播放高清视频 | 一级,二级,三级黄色视频| 9热在线视频观看99| 精品免费久久久久久久清纯 | 五月天丁香电影| 国产福利在线免费观看视频| 老熟妇乱子伦视频在线观看| 交换朋友夫妻互换小说| 手机成人av网站| 中国美女看黄片| 成人特级黄色片久久久久久久 | 搡老乐熟女国产| 久久久精品免费免费高清| 成人国产av品久久久| 日韩熟女老妇一区二区性免费视频| 日本av手机在线免费观看| 99精品欧美一区二区三区四区| 久久影院123| 免费在线观看完整版高清| 9色porny在线观看| 18禁国产床啪视频网站| 久久人妻av系列| 国产精品免费一区二区三区在线 | 成人手机av| 999久久久精品免费观看国产| 嫁个100分男人电影在线观看| 成人特级黄色片久久久久久久 | xxxhd国产人妻xxx| 国产精品二区激情视频| 两个人看的免费小视频| 天堂动漫精品| 午夜福利视频精品| 一级,二级,三级黄色视频| 亚洲 国产 在线| 人妻久久中文字幕网| 久久久水蜜桃国产精品网| 在线十欧美十亚洲十日本专区| 亚洲色图 男人天堂 中文字幕| tube8黄色片| 色视频在线一区二区三区| 黄频高清免费视频| 婷婷成人精品国产| 自拍欧美九色日韩亚洲蝌蚪91| 精品少妇黑人巨大在线播放| 在线十欧美十亚洲十日本专区| 免费观看a级毛片全部| 啦啦啦免费观看视频1| 精品少妇久久久久久888优播| 最近最新中文字幕大全免费视频| 日韩一区二区三区影片| 亚洲色图 男人天堂 中文字幕| 一二三四在线观看免费中文在| 精品一区二区三卡| 欧美日韩视频精品一区| 好男人电影高清在线观看| 男人操女人黄网站| 在线观看66精品国产| 香蕉丝袜av| 久久久久网色| a级毛片在线看网站| 亚洲色图av天堂| 精品少妇久久久久久888优播| 精品乱码久久久久久99久播| 青青草视频在线视频观看| 亚洲成av片中文字幕在线观看| 亚洲欧美日韩高清在线视频 | 亚洲成国产人片在线观看| 欧美在线黄色| 国产亚洲午夜精品一区二区久久| 久久毛片免费看一区二区三区| 丰满迷人的少妇在线观看| 91精品三级在线观看| 黄色毛片三级朝国网站| 亚洲人成77777在线视频| 日韩免费高清中文字幕av| 一夜夜www| 一边摸一边抽搐一进一出视频| 久久精品成人免费网站| 欧美日韩黄片免| 国产精品美女特级片免费视频播放器 | 少妇的丰满在线观看| 十八禁网站网址无遮挡| 水蜜桃什么品种好| 日本vs欧美在线观看视频| 99精品欧美一区二区三区四区| 女性被躁到高潮视频| 精品熟女少妇八av免费久了| 啦啦啦免费观看视频1| 啦啦啦视频在线资源免费观看| 中国美女看黄片| 精品亚洲成国产av| 国产精品电影一区二区三区 | www.熟女人妻精品国产| 国产av精品麻豆| 精品国产亚洲在线| 高清av免费在线| 亚洲av电影在线进入| 妹子高潮喷水视频| 亚洲avbb在线观看| 99国产综合亚洲精品| 999精品在线视频| 多毛熟女@视频| 久久中文字幕一级| 在线亚洲精品国产二区图片欧美| 男女高潮啪啪啪动态图| 狠狠婷婷综合久久久久久88av| 久久久久久人人人人人| 正在播放国产对白刺激| 18禁美女被吸乳视频| 欧美+亚洲+日韩+国产| 久久热在线av| cao死你这个sao货| 精品亚洲乱码少妇综合久久| 青青草视频在线视频观看|