• <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ì)
    可以在线观看的亚洲视频| 午夜福利欧美成人| 国产1区2区3区精品| 亚洲avbb在线观看| 美女 人体艺术 gogo| 日本一二三区视频观看| 人成视频在线观看免费观看| 最近视频中文字幕2019在线8| 久久精品影院6| 在线永久观看黄色视频| 国产三级在线视频| 国产成+人综合+亚洲专区| 久久久精品国产亚洲av高清涩受| 男女视频在线观看网站免费 | 久久精品成人免费网站| 亚洲av日韩精品久久久久久密| 国产精品一区二区三区四区久久| 十八禁网站免费在线| 亚洲精品美女久久av网站| 九色成人免费人妻av| 免费观看精品视频网站| 精品一区二区三区四区五区乱码| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久久久免费视频了| 亚洲专区字幕在线| 亚洲在线自拍视频| 国产日本99.免费观看| a级毛片在线看网站| 长腿黑丝高跟| 村上凉子中文字幕在线| 麻豆成人av在线观看| 少妇熟女aⅴ在线视频| 一卡2卡三卡四卡精品乱码亚洲| tocl精华| 免费av毛片视频| av欧美777| 亚洲美女视频黄频| 免费在线观看影片大全网站| 欧美久久黑人一区二区| 久久精品成人免费网站| 亚洲成人国产一区在线观看| 久久精品91蜜桃| 国产乱人伦免费视频| 国产亚洲精品av在线| 国内精品久久久久精免费| 亚洲精品av麻豆狂野| 男女那种视频在线观看| 日韩欧美国产一区二区入口| 久久精品亚洲精品国产色婷小说| xxx96com| 高清在线国产一区| 成年女人毛片免费观看观看9| 亚洲人成伊人成综合网2020| 丝袜人妻中文字幕| a级毛片在线看网站| 好男人在线观看高清免费视频| 亚洲五月天丁香| a级毛片在线看网站| 嫩草影院精品99| 亚洲国产欧洲综合997久久,| 巨乳人妻的诱惑在线观看| 日本免费一区二区三区高清不卡| 12—13女人毛片做爰片一| 午夜福利高清视频| 亚洲电影在线观看av| 三级国产精品欧美在线观看 | 久久伊人香网站| 国产激情久久老熟女| 亚洲国产精品久久男人天堂| 一本精品99久久精品77| 国产精品 欧美亚洲| 脱女人内裤的视频| 亚洲午夜理论影院| 99国产精品一区二区蜜桃av| 波多野结衣高清作品| 国产爱豆传媒在线观看 | 亚洲精品在线观看二区| 丁香六月欧美| 亚洲人与动物交配视频| 天堂√8在线中文| 嫩草影视91久久| 操出白浆在线播放| 欧美在线黄色| 国产一区二区三区视频了| 搞女人的毛片| 午夜日韩欧美国产| 国产aⅴ精品一区二区三区波| 在线看三级毛片| 亚洲人与动物交配视频| 欧美人与性动交α欧美精品济南到| bbb黄色大片| 一个人免费在线观看电影 | 亚洲国产精品sss在线观看| 18禁观看日本| 国内精品一区二区在线观看| 国产激情欧美一区二区| 桃色一区二区三区在线观看| 亚洲aⅴ乱码一区二区在线播放 | 成人欧美大片| 观看免费一级毛片| 桃色一区二区三区在线观看| 国产成人av教育| 两性夫妻黄色片| 最近最新免费中文字幕在线| 亚洲九九香蕉| 亚洲熟妇熟女久久| 亚洲激情在线av| 欧美在线黄色| 村上凉子中文字幕在线| 欧美+亚洲+日韩+国产| 欧美三级亚洲精品| 激情在线观看视频在线高清| 一区二区三区国产精品乱码| 男女那种视频在线观看| 欧美中文综合在线视频| 夜夜夜夜夜久久久久| 麻豆国产av国片精品| 男插女下体视频免费在线播放| 亚洲aⅴ乱码一区二区在线播放 | 国产精品 欧美亚洲| 成人国语在线视频| 国产成人av教育| 99在线视频只有这里精品首页| 在线观看午夜福利视频| 一本久久中文字幕| 国产一区二区在线av高清观看| 亚洲熟妇熟女久久| 国产一区二区激情短视频| 国产三级黄色录像| 国产真人三级小视频在线观看| 国内揄拍国产精品人妻在线| 国产精品一区二区三区四区久久| 人成视频在线观看免费观看| 国产亚洲精品av在线| 中文在线观看免费www的网站 | 日韩精品中文字幕看吧| 国产亚洲精品一区二区www| 韩国av一区二区三区四区| 亚洲成a人片在线一区二区| 久久久精品国产亚洲av高清涩受| 亚洲欧洲精品一区二区精品久久久| 精品久久久久久久人妻蜜臀av| 淫妇啪啪啪对白视频| 亚洲五月婷婷丁香| 免费在线观看完整版高清| 国产精品久久久久久久电影 | 特大巨黑吊av在线直播| 国产成人精品无人区| 久久精品国产亚洲av香蕉五月| 国产精品久久久久久久电影 | 精品一区二区三区视频在线观看免费| 国产午夜福利久久久久久| 一本综合久久免费| 人人妻人人看人人澡| 中文资源天堂在线| 黄色毛片三级朝国网站| 国产激情欧美一区二区| 99国产精品99久久久久| 国产黄片美女视频| 又粗又爽又猛毛片免费看| 成人国产一区最新在线观看| 全区人妻精品视频| 成人av一区二区三区在线看| 18禁观看日本| 亚洲自拍偷在线| 18禁黄网站禁片免费观看直播| 久久人人精品亚洲av| 亚洲片人在线观看| 人人妻人人澡欧美一区二区| ponron亚洲| 精品电影一区二区在线| 悠悠久久av| 麻豆成人av在线观看| 一本综合久久免费| 日韩精品青青久久久久久| 美女免费视频网站| 欧美日韩中文字幕国产精品一区二区三区| 久久久国产成人免费| 国产av一区在线观看免费| 午夜久久久久精精品| 亚洲中文日韩欧美视频| 麻豆成人av在线观看| 午夜福利18| 免费高清视频大片| 国产亚洲精品久久久久5区| 麻豆成人午夜福利视频| 91麻豆av在线| 真人做人爱边吃奶动态| 国产视频一区二区在线看| 日韩成人在线观看一区二区三区| 久久精品影院6| 国产精品久久视频播放| 久久久久久亚洲精品国产蜜桃av| 精品国产乱子伦一区二区三区| 老司机在亚洲福利影院| 免费看a级黄色片| 亚洲精品美女久久久久99蜜臀| 亚洲国产精品成人综合色| x7x7x7水蜜桃| 淫妇啪啪啪对白视频| 婷婷精品国产亚洲av| 精品人妻1区二区| 美女高潮喷水抽搐中文字幕| 99国产精品一区二区蜜桃av| 免费在线观看成人毛片| 91大片在线观看| 亚洲午夜精品一区,二区,三区| 日韩精品青青久久久久久| 国产成年人精品一区二区| 男女视频在线观看网站免费 | 亚洲中文字幕一区二区三区有码在线看 | 欧美丝袜亚洲另类 | av超薄肉色丝袜交足视频| 国产亚洲欧美98| 一级毛片精品| 女同久久另类99精品国产91| 久久人人精品亚洲av| 国产一区二区激情短视频| 青草久久国产| 久久九九热精品免费| 非洲黑人性xxxx精品又粗又长| 91大片在线观看| 777久久人妻少妇嫩草av网站| 久久国产乱子伦精品免费另类| 天堂√8在线中文| 免费看十八禁软件| 亚洲avbb在线观看| 亚洲国产精品999在线| 一边摸一边做爽爽视频免费| 亚洲成人中文字幕在线播放| or卡值多少钱| 久久久国产成人免费| 中文在线观看免费www的网站 | 日本免费a在线| 99久久国产精品久久久| 午夜精品一区二区三区免费看| 成年版毛片免费区| 美女免费视频网站| 亚洲一码二码三码区别大吗| av片东京热男人的天堂| 久久午夜综合久久蜜桃| 精品福利观看| 国产在线观看jvid| 琪琪午夜伦伦电影理论片6080| 在线观看免费日韩欧美大片| 午夜成年电影在线免费观看| 久久久久免费精品人妻一区二区| 18禁观看日本| 老汉色av国产亚洲站长工具| 久久天躁狠狠躁夜夜2o2o| 老司机午夜福利在线观看视频| 91九色精品人成在线观看| 久久久久久九九精品二区国产 | 99热这里只有精品一区 | 男女那种视频在线观看| 精品福利观看| 一区福利在线观看| 99国产精品一区二区蜜桃av| 人妻久久中文字幕网| 搞女人的毛片| 亚洲成人精品中文字幕电影| 俄罗斯特黄特色一大片| 很黄的视频免费| 啦啦啦观看免费观看视频高清| 妹子高潮喷水视频| 91麻豆精品激情在线观看国产| 亚洲黑人精品在线| 婷婷丁香在线五月| 日本熟妇午夜| 五月伊人婷婷丁香| 亚洲一码二码三码区别大吗| 欧美一级毛片孕妇| 韩国av一区二区三区四区| 亚洲真实伦在线观看| 国产精品av久久久久免费| 99久久国产精品久久久| 精品高清国产在线一区| 久久久久久九九精品二区国产 | 国产亚洲精品av在线| 美女午夜性视频免费| 久久久久久久精品吃奶| 亚洲欧美精品综合久久99| 日本一本二区三区精品| 亚洲人成伊人成综合网2020| 久久久久久亚洲精品国产蜜桃av| 亚洲男人的天堂狠狠| 高潮久久久久久久久久久不卡| 精品久久久久久久人妻蜜臀av| 亚洲欧美日韩东京热| 麻豆成人av在线观看| 亚洲专区中文字幕在线| av国产免费在线观看| 亚洲18禁久久av| 亚洲九九香蕉| 日韩大码丰满熟妇| 老司机午夜福利在线观看视频| 久久 成人 亚洲| 亚洲国产高清在线一区二区三| 国产一区在线观看成人免费| 18美女黄网站色大片免费观看| 成人18禁高潮啪啪吃奶动态图| 亚洲中文字幕日韩| 波多野结衣巨乳人妻| 国产一区在线观看成人免费| 男插女下体视频免费在线播放| 波多野结衣巨乳人妻| 亚洲精品在线美女| 久热爱精品视频在线9| 亚洲无线在线观看| 欧美国产日韩亚洲一区| www.自偷自拍.com| 99国产精品一区二区蜜桃av| 日本一本二区三区精品| 神马国产精品三级电影在线观看 | 中文字幕最新亚洲高清| 亚洲全国av大片| 成人午夜高清在线视频| 国产成人av激情在线播放| 欧美日本亚洲视频在线播放| 中国美女看黄片| 国产精品电影一区二区三区| 嫩草影视91久久| 欧美黑人巨大hd| 日日爽夜夜爽网站| 国产一区在线观看成人免费| 久久九九热精品免费| 国产伦一二天堂av在线观看| 国产精品av视频在线免费观看| 99精品久久久久人妻精品| 亚洲欧美日韩东京热| 男女床上黄色一级片免费看| 欧美黄色淫秽网站| 动漫黄色视频在线观看| 91av网站免费观看| 免费在线观看日本一区| 精品久久久久久,| 十八禁人妻一区二区| 真人做人爱边吃奶动态| 亚洲片人在线观看| 可以在线观看的亚洲视频| 99国产精品一区二区三区| 亚洲国产看品久久| 男人的好看免费观看在线视频 | 极品教师在线免费播放| www.熟女人妻精品国产| 成人国产一区最新在线观看| 免费人成视频x8x8入口观看| 看免费av毛片| www日本黄色视频网| 国产精华一区二区三区| 亚洲av片天天在线观看| 亚洲精品久久成人aⅴ小说| 在线观看免费午夜福利视频| 欧美日韩精品网址| 十八禁网站免费在线| 人妻久久中文字幕网| 妹子高潮喷水视频| 久久中文字幕一级| xxx96com| 成人三级黄色视频| 香蕉国产在线看| 99国产精品一区二区三区| 波多野结衣高清作品| 成年女人毛片免费观看观看9| 亚洲 欧美 日韩 在线 免费| 99国产综合亚洲精品| 久久午夜综合久久蜜桃| 此物有八面人人有两片| 午夜福利在线观看吧| 亚洲av成人一区二区三| 欧美午夜高清在线| 国产视频内射| 动漫黄色视频在线观看| 一夜夜www| 国内精品久久久久精免费| 国产一区在线观看成人免费| 久久久国产成人免费| 亚洲av五月六月丁香网| 亚洲精品粉嫩美女一区| 人妻久久中文字幕网| ponron亚洲| 中亚洲国语对白在线视频| 一级毛片高清免费大全| 一进一出好大好爽视频| 女人被狂操c到高潮| 我的老师免费观看完整版| 欧美成人性av电影在线观看| 淫妇啪啪啪对白视频| av在线播放免费不卡| 国产精品美女特级片免费视频播放器 | 婷婷丁香在线五月| 亚洲精品在线美女| 99国产精品一区二区蜜桃av| 欧美日韩乱码在线| 非洲黑人性xxxx精品又粗又长| 真人做人爱边吃奶动态| 国产亚洲精品第一综合不卡| 日本a在线网址| 欧美日韩亚洲国产一区二区在线观看| 婷婷精品国产亚洲av| 丰满人妻一区二区三区视频av | 国产精品美女特级片免费视频播放器 | 一本一本综合久久| cao死你这个sao货| 国产激情久久老熟女| 久热爱精品视频在线9| 亚洲国产欧洲综合997久久,| 又大又爽又粗| 99热这里只有是精品50| 一区二区三区高清视频在线| 舔av片在线| 首页视频小说图片口味搜索| 超碰成人久久| avwww免费| 日韩大码丰满熟妇| 国产激情欧美一区二区| 亚洲美女黄片视频| a级毛片在线看网站| 亚洲欧美日韩高清在线视频| 日韩欧美国产一区二区入口| 久久精品国产亚洲av高清一级| 日本三级黄在线观看| 人妻夜夜爽99麻豆av| av在线天堂中文字幕| 少妇的丰满在线观看| 日本成人三级电影网站| 亚洲精品粉嫩美女一区| 欧美中文日本在线观看视频| 三级男女做爰猛烈吃奶摸视频| 亚洲熟妇熟女久久| 女人被狂操c到高潮| 亚洲午夜精品一区,二区,三区| 免费av毛片视频| 毛片女人毛片| 女人爽到高潮嗷嗷叫在线视频| 巨乳人妻的诱惑在线观看| 这个男人来自地球电影免费观看| 黄色毛片三级朝国网站| 亚洲国产精品999在线| 麻豆av在线久日| 国产1区2区3区精品| 美女大奶头视频| 91国产中文字幕| 好看av亚洲va欧美ⅴa在| 国产精品久久久久久精品电影| 操出白浆在线播放| 国产精品一区二区三区四区免费观看 | 在线观看66精品国产| 久久久久久久午夜电影| 欧美另类亚洲清纯唯美| 每晚都被弄得嗷嗷叫到高潮| 丝袜美腿诱惑在线| 91大片在线观看| 国产男靠女视频免费网站| 免费看十八禁软件| 亚洲av片天天在线观看| 国产成人影院久久av| 亚洲欧美日韩高清在线视频| 亚洲国产精品999在线| 久久九九热精品免费| 国产av在哪里看| 国产精品久久视频播放| 久久久久久久久中文| 他把我摸到了高潮在线观看| 亚洲成人久久爱视频| 国语自产精品视频在线第100页| 亚洲国产欧美人成| 国产成人啪精品午夜网站| 亚洲全国av大片| 免费看十八禁软件| 91九色精品人成在线观看| 久久中文看片网| 最新在线观看一区二区三区| 日本黄色视频三级网站网址| 国产伦人伦偷精品视频| 脱女人内裤的视频| 欧美黑人精品巨大| 中文字幕精品亚洲无线码一区| 婷婷精品国产亚洲av| or卡值多少钱| 观看免费一级毛片| 一级a爱片免费观看的视频| 国产真实乱freesex| 久久久国产精品麻豆| 久久久久久久久久黄片| av福利片在线| 黄色 视频免费看| 男人的好看免费观看在线视频 | 国产高清视频在线观看网站| 女人被狂操c到高潮| 777久久人妻少妇嫩草av网站| 欧美性猛交黑人性爽| 欧美另类亚洲清纯唯美| 男的添女的下面高潮视频| 一夜夜www| 成人美女网站在线观看视频| 99视频精品全部免费 在线| 色哟哟·www| 欧美丝袜亚洲另类| 在线国产一区二区在线| 国产不卡一卡二| 精品人妻熟女av久视频| 人人妻人人澡人人爽人人夜夜 | 国产黄色视频一区二区在线观看 | 两性午夜刺激爽爽歪歪视频在线观看| 亚洲自偷自拍三级| 成熟少妇高潮喷水视频| 男女啪啪激烈高潮av片| 免费搜索国产男女视频| 亚洲精华国产精华液的使用体验 | 国产乱人偷精品视频| 亚洲欧美日韩无卡精品| 亚洲国产精品成人久久小说 | 久久久精品94久久精品| 久久草成人影院| 一区二区三区高清视频在线| 69av精品久久久久久| 国产av麻豆久久久久久久| 国产高清激情床上av| 秋霞在线观看毛片| 熟女电影av网| 日本av手机在线免费观看| 大又大粗又爽又黄少妇毛片口| 国产一区二区在线观看日韩| 国产精品久久久久久精品电影小说 | 91av网一区二区| 亚洲欧美日韩高清专用| av福利片在线观看| av视频在线观看入口| 国产精品一区二区在线观看99 | 国产伦精品一区二区三区四那| 91aial.com中文字幕在线观看| 99riav亚洲国产免费| 精品少妇黑人巨大在线播放 | 国产精品久久久久久精品电影| 精品无人区乱码1区二区| 日本黄色视频三级网站网址| 日日干狠狠操夜夜爽| 亚洲人成网站在线观看播放| 亚洲综合色惰| 男插女下体视频免费在线播放| 亚洲av中文字字幕乱码综合| 国产成人91sexporn| 日本黄色片子视频| 边亲边吃奶的免费视频| 欧美xxxx黑人xx丫x性爽| 亚洲,欧美,日韩| 久久人人爽人人爽人人片va| 亚洲图色成人| 亚洲国产精品sss在线观看| 国产精品一区www在线观看| 亚洲欧美日韩东京热| 综合色av麻豆| 青春草视频在线免费观看| 亚洲欧洲国产日韩| 干丝袜人妻中文字幕| 国产精品综合久久久久久久免费| 亚洲av免费高清在线观看| 国产老妇女一区| 一个人观看的视频www高清免费观看| 日韩,欧美,国产一区二区三区 | 久久久久久大精品| 天堂av国产一区二区熟女人妻| 变态另类丝袜制服| 久久久国产成人精品二区| 亚洲乱码一区二区免费版| 少妇人妻精品综合一区二区 | 精品久久国产蜜桃| 99热精品在线国产| 一边亲一边摸免费视频| 大又大粗又爽又黄少妇毛片口| 五月玫瑰六月丁香| 老司机福利观看| 国产真实伦视频高清在线观看| 桃色一区二区三区在线观看| 2021天堂中文幕一二区在线观| a级毛色黄片| 亚洲真实伦在线观看| 夜夜夜夜夜久久久久| 久久99蜜桃精品久久| 国产69精品久久久久777片| 亚洲真实伦在线观看| 欧美3d第一页| 亚洲成人中文字幕在线播放| 久久久久性生活片| 2021天堂中文幕一二区在线观| 在现免费观看毛片| 成年女人看的毛片在线观看| 日本欧美国产在线视频| 麻豆久久精品国产亚洲av| 亚洲国产色片| 亚洲精品亚洲一区二区| 人妻系列 视频| 深夜a级毛片| 欧美另类亚洲清纯唯美| 久久婷婷人人爽人人干人人爱| 99久久成人亚洲精品观看| www.色视频.com| 国产精品麻豆人妻色哟哟久久 | 国产激情偷乱视频一区二区| 国产精品爽爽va在线观看网站| 日日摸夜夜添夜夜添av毛片| 黄色视频,在线免费观看| 中国美白少妇内射xxxbb| 日本av手机在线免费观看| 精品人妻偷拍中文字幕| 精品无人区乱码1区二区| 国产精品福利在线免费观看| 久久草成人影院| 婷婷六月久久综合丁香| 十八禁国产超污无遮挡网站| 天堂√8在线中文| 看黄色毛片网站| 日本免费a在线| 网址你懂的国产日韩在线|