• <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ì)
    村上凉子中文字幕在线| 波多野结衣高清无吗| 寂寞人妻少妇视频99o| 亚洲自拍偷在线| 99热只有精品国产| av天堂在线播放| 国产精品一及| 国模一区二区三区四区视频| 亚洲国产精品久久男人天堂| 免费看光身美女| 亚洲乱码一区二区免费版| 国产黄a三级三级三级人| 亚洲婷婷狠狠爱综合网| 美女高潮的动态| 久久这里有精品视频免费| 欧美日韩综合久久久久久| 少妇人妻一区二区三区视频| 欧美日韩乱码在线| 久久草成人影院| 亚洲丝袜综合中文字幕| 成人毛片60女人毛片免费| 热99在线观看视频| 日韩欧美精品免费久久| 内地一区二区视频在线| 波野结衣二区三区在线| 久久久成人免费电影| 99国产精品一区二区蜜桃av| 免费黄网站久久成人精品| 天天一区二区日本电影三级| 欧美日韩在线观看h| 一级毛片电影观看 | 最近中文字幕高清免费大全6| 在现免费观看毛片| 国产色婷婷99| ponron亚洲| 长腿黑丝高跟| 乱系列少妇在线播放| 99热这里只有是精品50| 99久久久亚洲精品蜜臀av| 99热只有精品国产| 一进一出抽搐动态| 日韩成人伦理影院| 中文欧美无线码| 我要看日韩黄色一级片| 婷婷色av中文字幕| 亚洲av二区三区四区| 一区福利在线观看| 日日干狠狠操夜夜爽| 久久精品久久久久久久性| 天堂网av新在线| 两性午夜刺激爽爽歪歪视频在线观看| 黄色欧美视频在线观看| av视频在线观看入口| 观看美女的网站| 18禁在线无遮挡免费观看视频| 免费人成视频x8x8入口观看| 日韩成人av中文字幕在线观看| 一夜夜www| 免费看美女性在线毛片视频| 变态另类丝袜制服| 久久久久久久午夜电影| 一级二级三级毛片免费看| 九九热线精品视视频播放| 欧美日韩乱码在线| 美女内射精品一级片tv| 禁无遮挡网站| 久久人人爽人人片av| 国产日本99.免费观看| 午夜精品一区二区三区免费看| 黄色一级大片看看| 亚洲一级一片aⅴ在线观看| 亚洲精品久久国产高清桃花| 免费人成在线观看视频色| 1000部很黄的大片| 青春草亚洲视频在线观看| 三级毛片av免费| 啦啦啦韩国在线观看视频| 两个人的视频大全免费| 嫩草影院精品99| 久久精品国产亚洲网站| 亚洲性久久影院| 最好的美女福利视频网| 69人妻影院| 91精品国产九色| 亚洲欧美日韩东京热| 3wmmmm亚洲av在线观看| 亚洲欧美清纯卡通| 欧美潮喷喷水| 97人妻精品一区二区三区麻豆| 搡老妇女老女人老熟妇| 亚洲第一电影网av| 精品一区二区三区人妻视频| 99九九线精品视频在线观看视频| 日本三级黄在线观看| 国产精品一区二区在线观看99 | 亚洲国产日韩欧美精品在线观看| 日韩欧美精品v在线| 狠狠狠狠99中文字幕| 亚洲久久久久久中文字幕| h日本视频在线播放| 久久久成人免费电影| 春色校园在线视频观看| 日韩在线高清观看一区二区三区| 国产一区二区三区av在线 | 欧美日韩一区二区视频在线观看视频在线 | 我的女老师完整版在线观看| 亚洲四区av| 99热精品在线国产| 亚洲欧美精品专区久久| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久电影中文字幕| 日韩欧美 国产精品| 国产精品三级大全| 一卡2卡三卡四卡精品乱码亚洲| 人人妻人人澡人人爽人人夜夜 | 免费观看在线日韩| 男女边吃奶边做爰视频| 午夜老司机福利剧场| 少妇丰满av| 五月伊人婷婷丁香| 亚洲自偷自拍三级| 成人特级av手机在线观看| 九九热线精品视视频播放| 日本爱情动作片www.在线观看| 久久久久网色| 在线免费十八禁| 亚洲欧美精品综合久久99| 国产 一区 欧美 日韩| av专区在线播放| 日韩欧美在线乱码| 在线观看66精品国产| 日本与韩国留学比较| 国产成年人精品一区二区| 欧美一区二区亚洲| 国产精品久久久久久精品电影小说 | 日本av手机在线免费观看| 18禁在线无遮挡免费观看视频| 久久精品夜色国产| 久久草成人影院| 国产日韩欧美在线精品| 亚洲经典国产精华液单| 午夜福利在线在线| 男人舔女人下体高潮全视频| 九九热线精品视视频播放| 99久久无色码亚洲精品果冻| 国产亚洲精品久久久com| 欧美另类亚洲清纯唯美| 久久久久久久久久成人| 不卡视频在线观看欧美| 免费无遮挡裸体视频| 欧美精品一区二区大全| 麻豆久久精品国产亚洲av| 欧美日本亚洲视频在线播放| 成人国产麻豆网| 99久久精品热视频| av在线天堂中文字幕| 男女视频在线观看网站免费| 久久精品夜色国产| 亚洲国产精品久久男人天堂| 国产精品99久久久久久久久| 国产乱人视频| 黄片无遮挡物在线观看| 亚洲精品456在线播放app| 九九爱精品视频在线观看| 亚洲在久久综合| 草草在线视频免费看| 国产亚洲av嫩草精品影院| 18+在线观看网站| 午夜a级毛片| 国产老妇伦熟女老妇高清| 亚洲无线观看免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 一本精品99久久精品77| 久久久午夜欧美精品| 日日摸夜夜添夜夜爱| 欧美一级a爱片免费观看看| 久久欧美精品欧美久久欧美| 午夜福利视频1000在线观看| 久久精品国产亚洲av涩爱 | 欧美变态另类bdsm刘玥| 久久精品国产亚洲av香蕉五月| 欧美色视频一区免费| 亚洲无线观看免费| 一边亲一边摸免费视频| 国内久久婷婷六月综合欲色啪| 我的老师免费观看完整版| 国产精品精品国产色婷婷| 日本三级黄在线观看| 日韩人妻高清精品专区| 一区二区三区免费毛片| 国产伦在线观看视频一区| 精品一区二区免费观看| 免费看av在线观看网站| 国产免费一级a男人的天堂| 久久综合国产亚洲精品| 亚洲人成网站高清观看| 欧美日韩在线观看h| 一本久久精品| 久久鲁丝午夜福利片| 免费不卡的大黄色大毛片视频在线观看 | 国产一区二区激情短视频| 国产一级毛片在线| 此物有八面人人有两片| 亚洲国产精品成人综合色| av在线老鸭窝| 一区福利在线观看| 女人被狂操c到高潮| 女同久久另类99精品国产91| 国产黄片视频在线免费观看| 国产av麻豆久久久久久久| 久久亚洲精品不卡| 免费看av在线观看网站| 久久精品人妻少妇| 亚洲四区av| 不卡视频在线观看欧美| 亚洲内射少妇av| 晚上一个人看的免费电影| 国产一级毛片七仙女欲春2| 一夜夜www| 欧美最黄视频在线播放免费| 91麻豆精品激情在线观看国产| 免费在线观看成人毛片| 两个人的视频大全免费| 国产成人影院久久av| 麻豆成人av视频| 在线观看一区二区三区| av又黄又爽大尺度在线免费看 | 国产成人aa在线观看| 亚洲天堂国产精品一区在线| 久久鲁丝午夜福利片| 日本免费一区二区三区高清不卡| 久久久久久久久久黄片| 美女内射精品一级片tv| 国产一级毛片在线| 中文资源天堂在线| 亚洲精品久久久久久婷婷小说 | 国产乱人视频| 97在线视频观看| 久久久久网色| 成人二区视频| 伦精品一区二区三区| 久久久久久久久久久免费av| 成人毛片60女人毛片免费| 波多野结衣高清无吗| 久久精品影院6| 嫩草影院新地址| 春色校园在线视频观看| 国产精品一区二区在线观看99 | 亚洲乱码一区二区免费版| 国产亚洲av嫩草精品影院| 亚洲国产精品合色在线| 成人漫画全彩无遮挡| 99热这里只有是精品50| 91麻豆精品激情在线观看国产| 性欧美人与动物交配| 国内精品美女久久久久久| 国产高清有码在线观看视频| 国产午夜精品论理片| 精品午夜福利在线看| 99热这里只有是精品在线观看| 国产成人精品久久久久久| 69人妻影院| 天堂av国产一区二区熟女人妻| 亚洲av男天堂| 成人亚洲精品av一区二区| 神马国产精品三级电影在线观看| 好男人在线观看高清免费视频| 久久鲁丝午夜福利片| 成人特级av手机在线观看| 观看免费一级毛片| 日日撸夜夜添| 国产久久久一区二区三区| av福利片在线观看| 男女视频在线观看网站免费| 久久精品国产亚洲av涩爱 | 18+在线观看网站| 亚洲自偷自拍三级| 观看免费一级毛片| 国产私拍福利视频在线观看| 又爽又黄无遮挡网站| 我的老师免费观看完整版| 久久午夜福利片| 两性午夜刺激爽爽歪歪视频在线观看| 美女xxoo啪啪120秒动态图| 黄色欧美视频在线观看| 日韩三级伦理在线观看| 精品99又大又爽又粗少妇毛片| av.在线天堂| 夜夜爽天天搞| 老熟妇乱子伦视频在线观看| 成人午夜高清在线视频| 欧美日本视频| 九九热线精品视视频播放| 欧美一级a爱片免费观看看| 国产蜜桃级精品一区二区三区| 蜜桃亚洲精品一区二区三区| 麻豆久久精品国产亚洲av| 免费看美女性在线毛片视频| 日韩视频在线欧美| 国产精品久久久久久久电影| 精品国内亚洲2022精品成人| 看免费成人av毛片| 国产精品免费一区二区三区在线| 国产精品一区www在线观看| 国产久久久一区二区三区| 九草在线视频观看| 亚洲欧美日韩无卡精品| 国产高清有码在线观看视频| 国产亚洲欧美98| 成人综合一区亚洲| 日本黄色视频三级网站网址| av国产免费在线观看| 欧美性感艳星| 欧美3d第一页| 精品久久久久久久久久免费视频| 人妻夜夜爽99麻豆av| 麻豆国产97在线/欧美| 午夜爱爱视频在线播放| 亚洲欧美成人精品一区二区| 色吧在线观看| 久久99热6这里只有精品| 看十八女毛片水多多多| 女人十人毛片免费观看3o分钟| 在线国产一区二区在线| 亚洲无线在线观看| 亚洲精华国产精华液的使用体验 | av在线老鸭窝| 亚洲av.av天堂| 精品99又大又爽又粗少妇毛片| 少妇人妻一区二区三区视频| 极品教师在线视频| 少妇丰满av| 国产精品久久久久久久久免| eeuss影院久久| 青春草国产在线视频 | 久久国内精品自在自线图片| 国产精品久久久久久久电影| 男人和女人高潮做爰伦理| 国产老妇女一区| 一个人观看的视频www高清免费观看| 亚洲欧美日韩东京热| 一级毛片久久久久久久久女| 久久精品国产亚洲av涩爱 | 国产熟女欧美一区二区| 国产亚洲av嫩草精品影院| 国产熟女欧美一区二区| 亚洲四区av| 免费观看在线日韩| 日韩亚洲欧美综合| 麻豆久久精品国产亚洲av| 国内精品宾馆在线| 亚洲一区二区三区色噜噜| 给我免费播放毛片高清在线观看| 成人特级黄色片久久久久久久| 亚洲高清免费不卡视频| av天堂中文字幕网| 久久久久九九精品影院| 国产高清不卡午夜福利| 精品久久久久久久末码| 三级经典国产精品| 日韩av在线大香蕉| 99久久九九国产精品国产免费| 日日摸夜夜添夜夜爱| 国产精品国产高清国产av| 成人午夜高清在线视频| .国产精品久久| 久久韩国三级中文字幕| 三级男女做爰猛烈吃奶摸视频| 好男人视频免费观看在线| 一个人看的www免费观看视频| 日日摸夜夜添夜夜爱| 国产精品一区二区三区四区免费观看| 亚洲五月天丁香| 国产精品一区二区三区四区免费观看| 亚洲婷婷狠狠爱综合网| av.在线天堂| 国产黄a三级三级三级人| 欧美最黄视频在线播放免费| 九九热线精品视视频播放| 波多野结衣巨乳人妻| 成人特级av手机在线观看| 亚洲人成网站在线观看播放| 嫩草影院精品99| 久久久久免费精品人妻一区二区| 69人妻影院| 黄色一级大片看看| 成人无遮挡网站| 在线天堂最新版资源| 国产亚洲av嫩草精品影院| 不卡一级毛片| 午夜激情欧美在线| 少妇人妻一区二区三区视频| 内地一区二区视频在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 中文欧美无线码| 久久精品国产99精品国产亚洲性色| 国产一区二区三区av在线 | 五月玫瑰六月丁香| 精品久久久久久久人妻蜜臀av| 久久久精品94久久精品| 久久久久免费精品人妻一区二区| 成人午夜高清在线视频| 女同久久另类99精品国产91| 国产精品嫩草影院av在线观看| 亚洲精品日韩av片在线观看| 国产精品蜜桃在线观看 | 青春草国产在线视频 | 精品无人区乱码1区二区| 日本五十路高清| 日本黄色片子视频| 日韩中字成人| 午夜久久久久精精品| 我的女老师完整版在线观看| 精品国内亚洲2022精品成人| 舔av片在线| 欧美日韩国产亚洲二区| 一个人观看的视频www高清免费观看| 在线a可以看的网站| 18禁裸乳无遮挡免费网站照片| 九九在线视频观看精品| 国产精品一区www在线观看| 2022亚洲国产成人精品| 亚洲av不卡在线观看| 国内揄拍国产精品人妻在线| 中文在线观看免费www的网站| 国产高清视频在线观看网站| 久久久欧美国产精品| 中国国产av一级| 亚洲性久久影院| 久久6这里有精品| 一级黄色大片毛片| 在线免费十八禁| 国产av在哪里看| 女的被弄到高潮叫床怎么办| 色吧在线观看| 赤兔流量卡办理| 久久综合国产亚洲精品| 人人妻人人澡人人爽人人夜夜 | 亚洲国产欧美在线一区| 亚洲av中文字字幕乱码综合| 日本av手机在线免费观看| 婷婷精品国产亚洲av| 99riav亚洲国产免费| 悠悠久久av| 日本撒尿小便嘘嘘汇集6| 亚洲精品成人久久久久久| 内地一区二区视频在线| 亚洲精品色激情综合| 日韩一区二区视频免费看| 噜噜噜噜噜久久久久久91| 日韩,欧美,国产一区二区三区 | 精品人妻一区二区三区麻豆| 欧美精品国产亚洲| 欧美日韩国产亚洲二区| videossex国产| 成人毛片60女人毛片免费| 免费在线观看成人毛片| 亚洲欧美成人综合另类久久久 | 日韩欧美 国产精品| a级毛色黄片| eeuss影院久久| 国产 一区 欧美 日韩| 蜜臀久久99精品久久宅男| 一个人免费在线观看电影| 亚洲人成网站高清观看| 日韩制服骚丝袜av| 极品教师在线视频| 久久99蜜桃精品久久| 国产成人aa在线观看| 欧美精品一区二区大全| 日本欧美国产在线视频| 你懂的网址亚洲精品在线观看 | 老司机影院成人| 插逼视频在线观看| 国产成人a∨麻豆精品| 看免费成人av毛片| 日韩成人av中文字幕在线观看| 久久久久久久久久成人| 两个人视频免费观看高清| 国产熟女欧美一区二区| 少妇熟女欧美另类| 国产精品av视频在线免费观看| 又黄又爽又刺激的免费视频.| 亚洲国产欧美人成| 亚洲在久久综合| 亚洲va在线va天堂va国产| 又粗又爽又猛毛片免费看| 国产v大片淫在线免费观看| 免费电影在线观看免费观看| 看片在线看免费视频| 综合色丁香网| 在线天堂最新版资源| 中文字幕精品亚洲无线码一区| 色噜噜av男人的天堂激情| 精品久久久噜噜| 久久久精品欧美日韩精品| 国语自产精品视频在线第100页| 99在线视频只有这里精品首页| 村上凉子中文字幕在线| 国产亚洲av片在线观看秒播厂 | 亚洲成人精品中文字幕电影| 精品人妻一区二区三区麻豆| 日韩强制内射视频| 白带黄色成豆腐渣| 精品人妻熟女av久视频| 五月伊人婷婷丁香| 青春草视频在线免费观看| 日本av手机在线免费观看| 亚洲国产高清在线一区二区三| 韩国av在线不卡| 草草在线视频免费看| 99久久精品一区二区三区| 欧美zozozo另类| 全区人妻精品视频| 国产亚洲精品久久久com| 精品一区二区三区人妻视频| 国产精品一及| 亚洲av免费高清在线观看| 亚洲性久久影院| 搡女人真爽免费视频火全软件| 午夜精品国产一区二区电影 | 欧美最新免费一区二区三区| 日韩av在线大香蕉| 免费无遮挡裸体视频| 美女内射精品一级片tv| 亚洲成人av在线免费| 18禁裸乳无遮挡免费网站照片| 熟女电影av网| 99久久精品国产国产毛片| 两个人的视频大全免费| 午夜爱爱视频在线播放| 又粗又爽又猛毛片免费看| 麻豆精品久久久久久蜜桃| 亚洲自拍偷在线| 一卡2卡三卡四卡精品乱码亚洲| 精品国产三级普通话版| 国产一级毛片在线| 在线观看一区二区三区| 久久精品国产亚洲av天美| 午夜a级毛片| av黄色大香蕉| 五月玫瑰六月丁香| 美女被艹到高潮喷水动态| a级毛片a级免费在线| 国产探花在线观看一区二区| 成年女人永久免费观看视频| 亚洲无线观看免费| 校园人妻丝袜中文字幕| 国产亚洲精品久久久久久毛片| 日本黄色片子视频| 精品午夜福利在线看| 成人高潮视频无遮挡免费网站| 久久久色成人| 神马国产精品三级电影在线观看| 国产精品伦人一区二区| 热99在线观看视频| 免费av观看视频| 日韩高清综合在线| 日韩欧美 国产精品| 国产日本99.免费观看| 在线观看一区二区三区| 中文字幕av在线有码专区| 18禁裸乳无遮挡免费网站照片| 黄片无遮挡物在线观看| 亚洲三级黄色毛片| 免费搜索国产男女视频| 国产三级中文精品| av黄色大香蕉| 日韩成人av中文字幕在线观看| 国产免费一级a男人的天堂| 国产91av在线免费观看| 淫秽高清视频在线观看| 成人特级黄色片久久久久久久| 亚洲精品456在线播放app| 九色成人免费人妻av| 亚洲经典国产精华液单| 青春草国产在线视频 | 在线播放国产精品三级| 免费看美女性在线毛片视频| 1000部很黄的大片| 色尼玛亚洲综合影院| 欧美一区二区亚洲| 亚洲成人久久性| 99久久人妻综合| 青春草亚洲视频在线观看| 寂寞人妻少妇视频99o| 成年av动漫网址| 黄色视频,在线免费观看| av.在线天堂| 99在线人妻在线中文字幕| 波野结衣二区三区在线| 欧美丝袜亚洲另类| 国产私拍福利视频在线观看| 日本熟妇午夜| 51国产日韩欧美| 97超视频在线观看视频| 99久久九九国产精品国产免费| 久久精品国产亚洲网站| 国产精品一区二区三区四区久久| 看十八女毛片水多多多| 成年免费大片在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲国产欧洲综合997久久,| 久久久久久久久久久免费av| 久久99热6这里只有精品| 哪里可以看免费的av片| 日本成人三级电影网站| 一级黄片播放器| 亚洲精品色激情综合| 国产黄片美女视频| 日韩一区二区三区影片| 国产av不卡久久| 小说图片视频综合网站| 欧美成人免费av一区二区三区| 亚洲最大成人av|