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

    發(fā)動機工況下透平級凹槽葉頂冷卻傳熱性能研究

    2021-07-13 14:37:44于金杏葉明亮晏鑫
    西安交通大學(xué)學(xué)報 2021年7期
    關(guān)鍵詞:葉頂尾緣氣膜

    于金杏,葉明亮,晏鑫

    (西安交通大學(xué)能源與動力工程學(xué)院,710049,西安)

    凹槽葉頂是目前高壓透平中應(yīng)用最為廣泛的葉頂結(jié)構(gòu)[1]。隨著透平進口溫度的升高,凹槽葉頂?shù)臒嶝摵刹粩嘣龃?需要采用葉頂冷卻技術(shù)來保證透平的運行安全。氣膜冷卻是重要的高效葉頂冷卻技術(shù)[2],對于增強熱端部件的熱防護,確保透平的運行安全起到了關(guān)鍵作用。對凹槽葉頂氣膜冷卻和傳熱性能開展研究具有十分重要的學(xué)術(shù)和工程研究價值。

    葉頂氣膜孔的布置方案是影響凹槽葉頂氣膜冷卻有效性的重要因素之一[2]。Kwak等的實驗研究發(fā)現(xiàn):當(dāng)氣膜孔布置在凹槽底部中弧線位置時,冷卻流可以對凹槽底部壓力側(cè)及尾緣形成良好的冷卻[3]。葉明亮等的數(shù)值研究發(fā)現(xiàn):在帶隔板的凹槽葉頂中,中弧線氣膜孔能夠有效降低葉頂整體熱負荷[4]。Jeong等利用一維瞬態(tài)液晶技術(shù)對凹槽底部吸力側(cè)存在氣膜孔的凹槽葉頂進行了實驗研究,發(fā)現(xiàn)靠近吸力側(cè)的冷卻流在凹槽腔室渦的作用下對凹槽底部形成了有效的覆蓋[5]。Sakaoglu等采用數(shù)值方法研究了中弧線和凹槽底部壓力側(cè)氣膜孔直徑大小的影響,發(fā)現(xiàn)當(dāng)孔徑較大時冷卻流速度較低,氣膜覆蓋范圍較大,冷卻效果較好[6]。但是,在透平級環(huán)境中,凹槽葉頂?shù)睦鋮s傳熱性能會發(fā)生顯著變化。葉明亮等對存在中弧線氣膜孔的透平級凹槽葉頂進行了數(shù)值研究,發(fā)現(xiàn)葉頂旋轉(zhuǎn)效應(yīng)增強了凹槽底部吸力側(cè)的冷卻效果[7]。

    冷卻流流量對凹槽葉頂氣膜冷卻性能有重要的影響[8]。一般情況下,冷卻流會在氣膜孔下游形成冷卻氣膜,降低氣膜孔附近及凹槽底部的熱負荷[9]。然而,當(dāng)冷卻流流量過大時,冷卻流會脫離壁面并與主流發(fā)生強烈的相互作用,導(dǎo)致氣膜孔附近出現(xiàn)高傳熱區(qū)[10]。Ma等詳細研究了冷卻流和主流的作用過程,發(fā)現(xiàn)冷卻流形成的反向旋轉(zhuǎn)渦對使泄漏流直接沖擊凹槽底部,導(dǎo)致氣膜孔之間出現(xiàn)高傳熱區(qū)[11]。因此,合理地選擇冷卻流流量,可以在減少冷氣消耗量的同時提高凹槽葉頂?shù)臍饽だ鋮s效率。

    凹槽肩壁的熱負荷較高,對肩壁進行冷卻是凹槽葉頂冷卻傳熱研究的難點[1]。壓力面?zhèn)葰饽た资抢鋮s肩壁的重要方法之一,其冷卻流不僅可以有效降低兩側(cè)肩壁的熱負荷[12],還可以增強葉頂氣膜孔對凹槽內(nèi)部的冷卻效果[13]。然而,單獨采用壓力面?zhèn)葰饽た讜r,冷卻流對凹槽內(nèi)部和肩壁的冷卻效果較小[9];采用葉頂+壓力面?zhèn)葰饽た讜r,氣膜覆蓋范圍增大,冷卻效果增強[14]。黃琰等采用數(shù)值方法對葉頂氣膜孔和葉頂+壓力面?zhèn)葰饽た椎睦鋮s傳熱性能進行了研究,結(jié)果表明葉頂+壓力面?zhèn)葰饽た啄軌蛲瑫r冷卻凹槽底面壓力側(cè)、肩壁和葉頂尾緣[13]。

    但是,黃琰等的研究是在單列葉柵、實驗低溫工況下進行的,與發(fā)動機真實工況相差較遠[13]。為了研究發(fā)動機真實工況下的凹槽葉頂?shù)睦鋮s傳熱性能,本文在黃琰等研究的基礎(chǔ)上,對發(fā)動機運行工況下透平級凹槽葉頂?shù)睦鋮s傳熱性能開展相應(yīng)的數(shù)值研究[13]。

    首先,本文基于已有的實驗數(shù)據(jù)[15]驗證了數(shù)值方法的有效性,然后對發(fā)動機運行工況下GE-E3高壓透平第一級凹槽葉頂?shù)睦鋮s傳熱性能進行了數(shù)值計算,研究了3種吹風(fēng)比(M=0.5,1.0,2.0)時2種氣膜孔布置方案(單排葉頂氣膜孔、葉頂+壓力面?zhèn)葰饽た?條件下透平級凹槽葉頂?shù)睦鋮s傳熱性能。

    1 數(shù)值計算方法

    1.1 計算模型及邊界條件

    本文的研究對象為GE-E3高壓透平的第一級,包括46個靜葉和76個動葉,圖1給出了計算模型的動靜葉流道,葉片型線和氣膜孔幾何數(shù)據(jù)均來源于NASA報告[16]。

    圖1 計算模型的動靜葉流道

    本文研究的2種氣膜孔布置方案為單排葉頂氣膜孔以及葉頂+壓力面?zhèn)葰饽た?圖2給出了氣膜孔的布置情況。葉頂氣膜孔(T1~T8)的布置方案來源于NASA報告[16]。壓力面?zhèn)葰饽た?PS1~PS9)的布置方案參考Nasir的實驗[14],孔中心與葉頂?shù)木嚯x為葉高的6.67%,第一個氣膜孔中心位于17%軸向弦長處,相鄰氣膜孔的距離為10%軸向弦長,冷卻流與壓力面間的夾角為20°。所有氣膜孔均采用冷氣腔供氣。表1給出了計算模型的具體幾何參數(shù)。

    圖2 氣膜孔布置情況

    表1 計算模型幾何參數(shù)

    采用ANSYS ICEM CFD 11.0生成多塊結(jié)構(gòu)化網(wǎng)格,圖3給出了葉頂區(qū)域的網(wǎng)格示意圖。為了保證網(wǎng)格質(zhì)量,在葉片和氣膜孔周圍生成O網(wǎng)格。網(wǎng)格的最小角度為20°。為了保證邊界層內(nèi)的流場和溫度場計算精度,對所有近壁面區(qū)域進行加密,壁面的y+控制在y+<1的范圍。

    圖3 葉頂數(shù)值計算網(wǎng)格示意圖

    選擇發(fā)動機真實運行工況為研究工況,計算邊界條件設(shè)置列于表2,其中主流進口總溫、總壓,壓比,冷卻流進口總溫與NASA報告[16]保持一致。計算傳熱特性時,壁面溫度設(shè)置為1 273 K[16];計算氣膜冷卻特性時,壁面給定絕熱邊界條件。工質(zhì)為理想空氣,其動力黏度由Sutherlands公式確定

    表2 數(shù)值計算邊界條件

    (1)

    式中:Tref為參考溫度,取273 K;μref為工質(zhì)溫度為Tref時的動力黏度,取1.171 6×10-5(N·s)/m2;T為工質(zhì)溫度;S為Sutherlands常數(shù),取110.4 K。

    1.2 參數(shù)定義

    吹風(fēng)比定義為

    (2)

    式中:ρc為冷卻流進口密度;Vc為冷卻流進口速度;ρm為主流進口密度;Vm為主流進口速度。

    當(dāng)?shù)卮碉L(fēng)比定義為

    (3)

    式中:ρc,local為冷卻流當(dāng)?shù)孛芏?Vc,local為冷卻流當(dāng)?shù)厮俣?ρm,local為主流當(dāng)?shù)孛芏?Vm,local為主流當(dāng)?shù)厮俣?i為氣膜孔編號(T1~T8,PS1~PS9)。

    傳熱系數(shù)定義為

    (4)

    式中:q為壁面熱通量;Tw為壁面溫度;Tin為主流進口總溫。

    氣膜冷卻效率定義為

    (5)

    式中:Taw,0為不通冷氣時壁面的絕熱溫度;Taw,f為通冷氣時壁面的絕熱溫度;Tt,c為冷卻流進口總溫。

    1.3 數(shù)值方法考核

    采用CFD求解器ANSYS CFX 11.0對雷諾時均N-S方程組進行數(shù)值求解?;贙wak的實驗數(shù)據(jù)[15],考核了k-ε、k-ω和SST湍流模型的計算精度,圖4給出了3種湍流模型的計算結(jié)果。與實驗對比可以看出:k-ε湍流模型對肩壁和尾緣的傳熱系數(shù)預(yù)測偏高(見圖4中區(qū)域A、B);SST湍流模型對凹槽底部高傳熱區(qū)范圍預(yù)測偏大(見圖4中區(qū)域C)、對壓力側(cè)肩壁的傳熱系數(shù)預(yù)測偏低;k-ω湍流模型的計算結(jié)果與實驗結(jié)果吻合良好。

    (a)實驗[15] (b)k-ε湍流模型

    (c)k-ω湍流模型 (d)SST湍流模型

    圖5給出了采用葉頂+壓力面?zhèn)葰饽た讜r的實驗結(jié)果與k-ω湍流模型的計算結(jié)果。由圖可知,利用k-ω湍流模型能夠正確預(yù)測存在氣膜冷卻時葉頂?shù)膫鳠嵝阅?。因?本文最終采用k-ω湍流模型進行計算。

    (a)傳熱系數(shù)分布

    1.4 網(wǎng)格無關(guān)性驗證

    采用k-ω湍流模型對帶單排葉頂氣膜孔的透平級凹槽葉頂傳熱性能進行網(wǎng)格無關(guān)性驗證。為了確定動、靜葉網(wǎng)格數(shù)之間的關(guān)系,首先進行靜葉網(wǎng)格無關(guān)性驗證,然后在此基礎(chǔ)上進行級網(wǎng)格無關(guān)性驗證。

    進行靜葉網(wǎng)格無關(guān)性驗證時,動葉網(wǎng)格數(shù)為2 000萬,靜葉網(wǎng)格數(shù)分別為動葉網(wǎng)格數(shù)的1/5、2/5、3/5和4/5。計算得到的葉頂平均傳熱系數(shù)分別為1 688.04、1 665.53、1 665.33、1 663.28 W/(m2·K)。當(dāng)靜葉網(wǎng)格數(shù)達到800萬(即動葉網(wǎng)格數(shù)的2/5)時,計算結(jié)果與Richardson外推值[17]間的誤差小于0.18%,繼續(xù)增加網(wǎng)格數(shù)對計算結(jié)果的影響較小。因此,進行級網(wǎng)格無關(guān)性驗證時,靜葉網(wǎng)格數(shù)最終取動葉網(wǎng)格數(shù)的2/5。

    進行級網(wǎng)格無關(guān)性驗證時,級網(wǎng)格數(shù)分別為810萬、1 540萬、2 800萬和5 700萬。計算得到的葉頂平均傳熱系數(shù)分別為1 772.74、1 813.22、1 665.53、1 702.12 W/(m2·K)。當(dāng)級網(wǎng)格數(shù)達到2 800萬時,計算結(jié)果與Richardson外推值[17]間的誤差小于2.85%,繼續(xù)增加網(wǎng)格數(shù)對計算結(jié)果的影響較小。因此,最終確定帶單排葉頂氣膜孔的透平級網(wǎng)格數(shù)為2 800萬,帶葉頂+壓力面?zhèn)葰饽た椎耐钙郊壘W(wǎng)格數(shù)為3 500萬。

    2 結(jié)果與討論

    2.1 帶單排葉頂氣膜孔時凹槽葉頂?shù)睦鋮s傳熱性能

    2.1.1 吹風(fēng)比對當(dāng)?shù)卮碉L(fēng)比的影響 由于采用冷氣腔為多個氣膜孔供氣,所以各氣膜孔的當(dāng)?shù)卮碉L(fēng)比有所不同。為了研究不同氣膜孔的出流狀態(tài)對葉頂冷卻傳熱性能的影響,首先計算了不同吹風(fēng)比條件下各氣膜孔的當(dāng)?shù)卮碉L(fēng)比。

    圖6給出了吹風(fēng)比為0.5、1.0、2.0時各氣膜孔的當(dāng)?shù)卮碉L(fēng)比。由圖可知,當(dāng)?shù)卮碉L(fēng)比隨吹風(fēng)比的增大而增大,MT1~MT8基本呈現(xiàn)先增大后減小的趨勢,其中MT4最大。當(dāng)M較小時,各氣膜孔的Mi相近;當(dāng)M較大時,各氣膜孔的Mi差異增加。這是因為當(dāng)?shù)卮碉L(fēng)比受氣膜孔位置、氣膜孔尺寸和主流當(dāng)?shù)貕毫Φ挠绊?當(dāng)M較小時,雖然各氣膜孔的位置、孔徑大小不同,但由于冷卻流的流量較小,不同位置孔徑的變化對冷卻流流動狀態(tài)和凹槽內(nèi)部流場結(jié)構(gòu)的影響較小,氣膜孔的進出口壓差不會發(fā)生顯著變化,因此各氣膜孔的Mi相近。當(dāng)M較大時,冷卻流的流量增大,冷卻流的注入對凹槽內(nèi)部流場結(jié)構(gòu)的影響增加,因此氣膜孔位置、孔徑大小等因素對氣膜孔進出口壓差的影響更為顯著,各氣膜孔的Mi差異增加。

    圖6 帶單排葉頂氣膜孔時各氣膜孔的當(dāng)?shù)卮碉L(fēng)比

    2.1.2 凹槽葉頂?shù)闹饕鲌鼋Y(jié)構(gòu) 圖7給出了M=2.0時透平級凹槽葉頂?shù)娜S流線。由圖7a可知,葉頂流場結(jié)構(gòu)主要由5部分組成:

    (a)葉頂主要流場結(jié)構(gòu)

    (1)直接掠過肩壁的前緣泄漏流所形成的吸力側(cè)前緣渦;

    (2)進入凹槽的前緣泄漏流所形成的凹槽腔室渦;

    (3)由于端壁與葉頂之間存在相對運動,端壁附近的流體受到剪切力作用所形成的刮削渦;

    (4)由于刮削渦靠近端壁,阻礙了泄漏流的流動,一部分流體從刮削渦與凹槽腔室渦之間進入凹槽,沖擊凹槽底部,形成沖擊流動;

    (5)刮削渦掠過吸力側(cè)肩壁后,凹槽腔室內(nèi)的流動對泄漏流的阻礙作用消失,泄漏流直接流過葉頂形成穿越流。

    由圖7b可知,上述流體流出葉頂后形成或匯入上通道渦和泄漏渦。根據(jù)圖7b中的流動特性,可將流出區(qū)域分為3部分:Ⅰ區(qū)域的流體形成上通道渦;Ⅰ、Ⅱ交界處附近的流體主要形成泄漏渦;Ⅱ區(qū)域的流體掠過泄漏渦后匯入上通道渦;Ⅲ區(qū)域的流體主要匯入泄漏渦。

    圖8給出了M=2.0時透平級單排葉頂氣膜孔的冷卻流流線。結(jié)合圖7a和圖8可以看出,冷卻流的流動軌跡主要受凹槽腔室渦的影響:T1~T3的冷卻流首先匯集于凹槽前緣,然后向下游流動,最終掠過吸力側(cè)肩壁離開葉頂;T4~T8的冷卻流先向壓力側(cè)流動,然后掠過吸力側(cè)肩壁最終離開葉頂。由此可見,在氣膜孔布置方案的設(shè)計中,為了減小冷卻流與主流間的干涉作用,增加冷卻流的貼壁性,氣膜孔的射流方向應(yīng)盡可能與凹槽腔室渦的流動方向一致。在凹槽底面附近,凹槽腔室渦由吸力側(cè)流向壓力側(cè),因此葉頂氣膜孔的方向也應(yīng)指向壓力側(cè)。

    圖8 M=2.0時葉頂冷卻流的三維流線

    2.1.3 帶單排葉頂氣膜孔的凹槽葉頂傳熱性能 圖9給出了3種吹風(fēng)比條件下透平級凹槽葉頂?shù)膫鳠嵯禂?shù)分布,該分布具有以下4個特點。

    (1)凹槽底部前緣存在低傳熱區(qū)A。結(jié)合圖8可知,這是因為T1~T3的冷卻流匯集于凹槽前緣,形成了覆蓋整個前緣區(qū)域的冷卻氣膜。

    (2)凹槽底部壓力側(cè)的傳熱系數(shù)明顯低于吸力側(cè)。結(jié)合圖8可知,T4~T8的冷卻流從氣膜孔射出后向壓力側(cè)流動,然后脫離壁面從吸力側(cè)肩壁離開葉頂,因此冷卻流對凹槽底部壓力側(cè)的冷卻效果較好,對吸力側(cè)的冷卻效果較差。

    (3)凹槽底部吸力側(cè)存在高傳熱區(qū)B。結(jié)合圖7可知,這是由于高溫沖擊流從刮削渦和凹槽腔室渦之間進入凹槽腔室,沖擊凹槽底部。

    (4)吸力側(cè)肩壁C1、C2區(qū)域傳熱系數(shù)較小。這是由于一部分T1~T3的冷卻流從C1區(qū)域離開葉頂,其余大部分冷卻流從C2區(qū)域離開葉頂(見圖8)。

    (a)M=0.5

    圖10 凹槽底部節(jié)距平均傳熱系數(shù)沿軸向的分布

    圖11 凹槽肩壁節(jié)距平均傳熱系數(shù)沿軸向的分布

    2.1.4 帶單排葉頂氣膜孔的凹槽葉頂冷卻性能 圖12給出了3種吹風(fēng)比條件下透平級凹槽葉頂?shù)臍饽だ鋮s效率分布??梢钥闯?高氣膜冷卻效率區(qū)集中在凹槽底部前緣(D區(qū)域)、壓力側(cè)(E區(qū)域)及尾緣(F區(qū)域)附近。此外,由于T6~T8的冷卻流較為集中,擴散性較弱,氣膜孔下游存在氣膜冷卻效率極高的帶狀區(qū)域(G區(qū)域)。M=2.0時,由于大量的冷卻流與主流混合后降低了主流的溫度,凹槽底部吸力側(cè)的氣膜冷卻效率明顯增大。

    (a)M=0.5

    圖13 凹槽底部節(jié)距平均氣膜冷卻效率沿軸向的分布

    由以上分析可知,凹槽前緣(T1~T3附近)存在冷卻流匯集區(qū),采用較小孔徑的中弧線氣膜孔就可以獲得較好的冷卻效果。但是,由于前緣冷卻流能夠降低主流溫度,提高整體氣膜冷卻效率,因此冷卻流的流量不宜過小。綜合考慮以上兩點,應(yīng)該在凹槽前緣采用中等孔徑的中弧線氣膜孔;凹槽尾緣(T7~T8附近)易形成冷卻流匯集區(qū),采用較小孔徑的中弧線氣膜孔;凹槽中間區(qū)域(T4~T6附近)無法形成冷卻流匯集區(qū),不易于冷卻,且存在吸力側(cè)高傳熱系數(shù)區(qū)B,采用靠近吸力側(cè)的大孔徑氣膜孔。

    圖14 凹槽肩壁節(jié)距平均氣膜冷卻效率沿軸向的分布

    2.2 帶葉頂+壓力面?zhèn)葰饽た椎陌疾廴~頂冷卻傳熱性能

    2.2.1 吹風(fēng)比對當(dāng)?shù)卮碉L(fēng)比的影響 圖15給出了吹風(fēng)比為0.5、1.0、2.0時透平級凹槽葉頂各氣膜孔的當(dāng)?shù)卮碉L(fēng)比。由圖可知,MT1~MT8基本呈現(xiàn)先增大后減小的趨勢,其中MT4最大。吹風(fēng)比為0.5、1.0時,部分壓力面?zhèn)葰饽た椎漠?dāng)?shù)卮碉L(fēng)比為負(M=0.5時,MPS1~MPS7為負;M=1.0時,MPS1~MPS6為負),出現(xiàn)了主流吸入現(xiàn)象。M=2.0時,主流吸入現(xiàn)象消失,從壓力面?zhèn)葰饽た咨涑龅睦鋮s流流量顯著增加,導(dǎo)致部分葉頂氣膜孔的當(dāng)?shù)卮碉L(fēng)比減小。

    (a)各葉頂氣膜孔的當(dāng)?shù)卮碉L(fēng)比

    2.2.2 帶葉頂+壓力面?zhèn)葰饽た讜r冷卻流的流場結(jié)構(gòu) 圖16給出了M=2.0時冷卻流的三維流線。結(jié)合圖8與圖16可知,葉頂冷卻流的流動軌跡與帶單排葉頂氣膜孔時基本相同。壓力面?zhèn)壤鋮s流的流動軌跡受到泄漏流的影響,主要分為4種:①PS1的冷卻流進入凹槽內(nèi)部;②PS2~PS5的冷卻流從凹槽上方流過葉頂;③PS6~PS8的冷卻流從尾緣上方流過葉頂;④PS9的大部分冷卻流直接匯入尾跡流。

    圖16 M=2.0時葉頂及壓力面?zhèn)壤鋮s流的流線

    此外,不同流動軌跡的壓力面?zhèn)壤鋮s流對葉頂流場的影響不同:①PS1的冷卻流進入凹槽后匯入凹槽腔室渦;②PS2~PS5的冷卻流有利于阻礙凹槽內(nèi)部流體的泄漏;③PS6~PS8的冷卻流可以形成葉頂尾緣冷卻氣膜;④PS9的冷卻流對葉頂流場結(jié)構(gòu)的影響較小。由此可見,為了最大程度地發(fā)揮壓力面?zhèn)壤鋮s流的作用,壓力面?zhèn)葰饽た讘?yīng)該布置在PS9所在位置(軸向弦長的76%)的上游。

    2.2.3 帶葉頂+壓力面?zhèn)葰饽た椎陌疾廴~頂傳熱性能 圖17給出了3種吹風(fēng)比條件下透平級凹槽葉頂?shù)膫鳠嵯禂?shù)分布。與圖9相比,凹槽底部高傳熱區(qū)B的范圍及傳熱系數(shù)減小,這是因為MT1~MT3增大,匯入凹槽腔室渦的冷卻流增加,對凹槽底部吸力側(cè)的冷卻效果增強。M=2.0時,由于壓力面?zhèn)壤鋮s流的覆蓋,壓力側(cè)肩壁與葉頂尾緣的傳熱系數(shù)顯著減小。

    (a)M=0.5

    圖18給出了存在葉頂+壓力面?zhèn)葰饽た讜r透平級凹槽底部節(jié)距平均傳熱系數(shù)沿軸向的分布。與圖10相似,圖18中的曲線呈現(xiàn)先增大后減小的趨勢,但最大值明顯降低,這是因為高傳熱區(qū)B的范圍及傳熱系數(shù)減小(見圖17)。

    圖18 凹槽底部節(jié)距平均傳熱系數(shù)沿軸向的分布

    圖19 凹槽肩壁節(jié)距平均傳熱系數(shù)沿軸向的分布

    圖20給出了3種吹風(fēng)比條件下凹槽葉頂壓力面?zhèn)葌鳠嵯禂?shù)分布。當(dāng)氣膜孔中幾乎無冷卻流流出時,氣膜孔下游存在高傳熱區(qū)J。Haydt等的實驗研究[18]也觀察到了類似的現(xiàn)象。這是由于氣膜孔使主流發(fā)生了流動分離,主流與壁面間的對流換熱增加。M=2.0時,壓力面?zhèn)葰饽た紫掠未嬖趲畹蛡鳠釁^(qū),且越靠近尾緣,低傳熱區(qū)向尾緣方向的傾斜角越大。這是因為越靠近尾緣,壓力側(cè)與吸力側(cè)間的驅(qū)動壓差越小,冷卻流向吸力側(cè)方向流動的趨勢越小,向尾緣方向流動的趨勢越大。

    (a)M=0.5

    值得注意的是,有冷卻流流出的壓力面?zhèn)葰饽た赘浇嬖诟邆鳠釁^(qū)K(見圖20c)。由圖21可知,這是因為氣膜孔下端流出的冷卻流會形成出口漩渦,將主流卷吸至壁面附近,增加了主流與壁面間的對流換熱。越靠近尾緣,高傳熱區(qū)K的范圍越大,這是因為越靠近尾緣,主流向尾緣方向的流動趨勢越大,與冷卻流間的夾角越大。Aga和Abhari的研究表明,隨著冷卻流與主流間的夾角增大,出口漩渦的強度增強,對主流的卷吸作用增強,因此越靠近尾緣,高傳熱區(qū)K的范圍越大[19]。為了減小主流與冷卻流間的夾角,從而減小K的范圍,可以將靠近尾緣的壓力面?zhèn)葰饽た紫蛭簿壏较蛐D(zhuǎn)適當(dāng)?shù)慕嵌取?/p>

    圖21 M=2.0時PS8、PS9的冷卻流流線及附近的主流流線

    2.2.4 帶葉頂+壓力面?zhèn)葰饽た椎陌疾廴~頂冷卻性能 圖22給出了3種吹風(fēng)比條件下透平級凹槽葉頂氣膜冷卻效率分布。與圖12對比可發(fā)現(xiàn),M=0.5,1.0時,采用葉頂+壓力面?zhèn)葰饽た撞贾每梢杂行г龃蟀疾鄣撞课?cè)的氣膜冷卻效率。這是因為MT1~MT3增大,匯入凹槽腔室渦并到達凹槽底部吸力側(cè)的冷卻流增加。M=2.0時,壓力側(cè)肩壁與凹槽葉頂尾緣的氣膜冷卻效率顯著增大。由圖16可知,這是因為PS1~PS5和PS6~PS9的冷卻流分別形成了壓力側(cè)肩壁和葉頂尾緣冷卻氣膜。但是,由于PS1位于17%軸向弦長處,冷卻氣膜對前緣的覆蓋不足,因此應(yīng)該增加前緣壓力面?zhèn)葰饽た讈碓鰪娗熬壚鋮s效果。

    (a)M=0.5

    圖23給出了葉頂+壓力面?zhèn)葰饽た撞贾脮r透平級凹槽底部節(jié)距平均氣膜冷卻效率沿軸向的分布。圖中的曲線變化趨勢與圖13相似,但變化速度趨緩,說明氣膜冷卻效率的分布更加均勻。

    圖23 凹槽底部節(jié)距平均氣膜冷卻效率沿軸向的分布

    圖24給出了葉頂+壓力面?zhèn)葰饽た撞贾脮r凹槽肩壁節(jié)距平均氣膜冷卻效率沿軸向的分布。M=0.5,1.0時,圖中的曲線與帶單排葉頂氣膜孔時相似(見圖14);M=2.0時,由于壓力面?zhèn)壤鋮s流對壓力側(cè)肩壁及葉頂尾緣的冷卻,圖中的曲線顯著上移,但由于壓力側(cè)肩壁的冷卻氣膜是分散分布的(見圖22c),曲線呈波浪狀。

    圖24 凹槽肩壁節(jié)距平均氣膜冷卻效率沿軸向的分布

    圖25給出了3種吹風(fēng)比條件下凹槽葉頂壓力面?zhèn)鹊臍饽だ鋮s效率分布。M=0.5,1.0時,由于存在主流吸入現(xiàn)象,只有靠近尾緣的氣膜孔下游會形成冷卻氣膜;M=2.0時,主流吸入現(xiàn)象消失。因此,應(yīng)在采用壓力面?zhèn)葰饽た椎耐瑫r采用較大的吹風(fēng)比。M=2.0時,PS1~PS3的冷卻流向尾緣方向的偏轉(zhuǎn)角度較小,冷卻氣膜是分散分布的;PS4~PS9的冷卻流向尾緣方向的偏轉(zhuǎn)角度較大,冷卻氣膜是連續(xù)的。因此,為了增加冷卻氣膜的連續(xù)性,可以適當(dāng)減小前緣及中間區(qū)域氣膜孔間的距離。

    (a)M=0.5

    2.3 2種氣膜孔布置時凹槽葉頂冷卻傳熱性能對比

    圖26給出了2種氣膜孔布置條件下凹槽葉頂?shù)拿娣e平均傳熱系數(shù)和氣膜冷卻效率。與帶單排葉頂氣膜孔相比,葉頂+壓力面?zhèn)葰饽た撞贾脤ν钙郊壈疾鄣撞坷鋮s傳熱性能的影響不大,但M=2.0時,可以使凹槽肩壁的面積平均傳熱系數(shù)減小13.89%,平均氣膜冷卻效率增大61.45%。因此,雖然在吹風(fēng)比相同的條件下,葉頂+壓力面?zhèn)葰饽た撞贾玫睦錃庥昧扛?但壓力面?zhèn)壤鋮s流顯著增強了凹槽肩壁區(qū)域的氣膜冷卻效果。

    (a)凹槽底部

    3 結(jié) 論

    本文用數(shù)值方法研究了發(fā)動機運行工況下2種氣膜孔布置方案(單排葉頂氣膜孔、葉頂+壓力面?zhèn)葰饽た?時透平級凹槽葉頂?shù)睦鋮s傳熱性能,獲得了3種吹風(fēng)比(0.5、1.0和2.0)條件下凹槽葉頂?shù)膫鳠嵯禂?shù)和氣膜冷卻效率,主要結(jié)論如下。

    (1)凹槽葉頂?shù)闹饕鲌鼋Y(jié)構(gòu)包括吸力側(cè)前緣渦、凹槽腔室渦、刮削渦、沖擊流和穿越流5個部分。凹槽腔室渦使凹槽底部前緣和壓力側(cè)得到了有效的冷卻,沖擊流導(dǎo)致凹槽底部吸力側(cè)形成了高傳熱區(qū)。

    (2)采用單排葉頂氣膜孔時,冷卻流對凹槽底部前緣、壓力側(cè)及葉頂尾緣區(qū)域的冷卻效果較好。為了提高冷卻流的利用率,葉頂前緣區(qū)域應(yīng)該采用中等孔徑的中弧線氣膜孔,中間區(qū)域應(yīng)該采用靠近吸力側(cè)的較大孔徑氣膜孔,尾緣區(qū)域應(yīng)該采用較小孔徑的中弧線氣膜孔,所有氣膜孔的射流方向都應(yīng)指向壓力側(cè)。

    (3)吹風(fēng)比為2.0時,采用葉頂+壓力面?zhèn)葰饽た撞贾媚軌蝻@著增大凹槽葉頂壓力側(cè)肩壁和葉頂尾緣的冷卻效果。為了提高壓力面?zhèn)葰饽た椎睦鋮s性能,應(yīng)該增加前緣壓力面?zhèn)葰饽た?減小前緣及中間區(qū)域氣膜孔間的距離,將靠近尾緣的氣膜孔向尾緣方向旋轉(zhuǎn)適當(dāng)?shù)慕嵌取?/p>

    (4)壓力面?zhèn)壤鋮s流主要對壓力側(cè)肩壁和葉頂尾緣形成冷卻,葉頂冷卻流主要對腔室底部和吸力側(cè)肩壁形成冷卻。吹風(fēng)比為2.0時,與單排葉頂氣膜孔布置方案相比,葉頂+壓力面?zhèn)葰饽た椎牟贾梅桨改軌蚴拱疾奂绫诘拿娣e平均傳熱系數(shù)降低13.89%,氣膜冷卻效率增大61.45%。

    猜你喜歡
    葉頂尾緣氣膜
    平面葉柵多凹槽葉頂傾斜圓柱孔氣膜冷卻與氣動特性研究
    分離渦模擬在葉頂間隙流模擬中的可行性分析
    T 型槽柱面氣膜密封穩(wěn)態(tài)性能數(shù)值計算研究
    透平動葉多凹槽葉頂氣膜冷卻特性的研究
    基于強化換熱的偏斜尾緣設(shè)計
    能源工程(2021年1期)2021-04-13 02:05:46
    氣膜孔堵塞對葉片吸力面氣膜冷卻的影響
    靜葉柵上游端壁雙射流氣膜冷卻特性實驗
    火箭推進(2020年2期)2020-05-06 02:53:56
    翼型湍流尾緣噪聲半經(jīng)驗預(yù)測公式改進
    具有尾緣襟翼的風(fēng)力機動力學(xué)建模與恒功率控制
    躲避霧霾天氣的氣膜館
    亚洲一级一片aⅴ在线观看| 少妇熟女欧美另类| 一边亲一边摸免费视频| 波多野结衣高清无吗| 欧美最黄视频在线播放免费| 男人狂女人下面高潮的视频| 国产精品久久久久久av不卡| 国产精品,欧美在线| 国产高潮美女av| 日韩视频在线欧美| a级毛色黄片| 91在线精品国自产拍蜜月| 精品一区二区免费观看| 日本-黄色视频高清免费观看| 国产精品99久久久久久久久| 亚洲人成网站在线播| 最近2019中文字幕mv第一页| 女人十人毛片免费观看3o分钟| 久久精品国产清高在天天线| 亚洲电影在线观看av| 九九在线视频观看精品| 中文资源天堂在线| 国产久久久一区二区三区| 欧美极品一区二区三区四区| 99久久人妻综合| 成人一区二区视频在线观看| 精品人妻一区二区三区麻豆| 亚洲丝袜综合中文字幕| 国产欧美日韩精品一区二区| 日韩大尺度精品在线看网址| 午夜久久久久精精品| 好男人在线观看高清免费视频| 国产精品久久视频播放| 成人无遮挡网站| 午夜老司机福利剧场| 亚洲aⅴ乱码一区二区在线播放| а√天堂www在线а√下载| 国产精品三级大全| 我的老师免费观看完整版| 国产探花在线观看一区二区| 在线免费观看不下载黄p国产| 毛片女人毛片| 中出人妻视频一区二区| 天堂网av新在线| 亚洲经典国产精华液单| 国产精品久久久久久精品电影| 国产精品久久久久久精品电影| 国产精品,欧美在线| 噜噜噜噜噜久久久久久91| 女人十人毛片免费观看3o分钟| 给我免费播放毛片高清在线观看| 最近2019中文字幕mv第一页| 亚洲久久久久久中文字幕| 99热这里只有是精品50| 中文资源天堂在线| 女同久久另类99精品国产91| 午夜福利高清视频| 看非洲黑人一级黄片| 波多野结衣高清作品| 18禁在线播放成人免费| 一本久久精品| 青春草国产在线视频 | 99久久精品国产国产毛片| 免费大片18禁| 久久精品国产亚洲av香蕉五月| 丰满的人妻完整版| 国产精品免费一区二区三区在线| 国内揄拍国产精品人妻在线| 男女做爰动态图高潮gif福利片| 国产精品无大码| 日韩中字成人| 大香蕉久久网| 一级毛片aaaaaa免费看小| 久久久久性生活片| 欧美3d第一页| 日韩欧美国产在线观看| 国产在视频线在精品| 高清毛片免费看| 老司机福利观看| 欧美色欧美亚洲另类二区| 国产精品一区www在线观看| 岛国在线免费视频观看| 毛片一级片免费看久久久久| 日本爱情动作片www.在线观看| 国产精品国产三级国产av玫瑰| 高清日韩中文字幕在线| 亚州av有码| 久久亚洲精品不卡| 久久久午夜欧美精品| 免费观看在线日韩| 亚洲精品成人久久久久久| 少妇人妻一区二区三区视频| 日韩精品青青久久久久久| 亚洲精品日韩在线中文字幕 | 亚洲四区av| 99久国产av精品国产电影| 亚洲人与动物交配视频| 久久久国产成人精品二区| 国内久久婷婷六月综合欲色啪| 观看免费一级毛片| 国产乱人偷精品视频| 在线播放国产精品三级| 久久久a久久爽久久v久久| 国产精品久久久久久久久免| 99热6这里只有精品| 国产一级毛片在线| 午夜福利成人在线免费观看| 人妻久久中文字幕网| 国产久久久一区二区三区| 日韩大尺度精品在线看网址| 可以在线观看毛片的网站| 小蜜桃在线观看免费完整版高清| 精品国内亚洲2022精品成人| 国产伦在线观看视频一区| 久久鲁丝午夜福利片| 亚洲av二区三区四区| 亚洲人成网站在线播放欧美日韩| 一边亲一边摸免费视频| 亚洲美女视频黄频| 一级黄片播放器| 免费大片18禁| 国产精品女同一区二区软件| 欧美一区二区国产精品久久精品| 亚洲中文字幕日韩| 夫妻性生交免费视频一级片| 全区人妻精品视频| 国产精品一区二区性色av| 22中文网久久字幕| 日韩人妻高清精品专区| 日韩欧美 国产精品| 国产精品一区二区性色av| 日本-黄色视频高清免费观看| 观看免费一级毛片| 热99re8久久精品国产| 最近中文字幕高清免费大全6| 成年女人永久免费观看视频| av又黄又爽大尺度在线免费看 | 久久精品91蜜桃| 久久精品91蜜桃| 国产黄片美女视频| 高清午夜精品一区二区三区 | 我要搜黄色片| 国产一区二区在线av高清观看| 久久精品国产清高在天天线| 亚洲18禁久久av| 免费不卡的大黄色大毛片视频在线观看 | 变态另类成人亚洲欧美熟女| 男的添女的下面高潮视频| 免费av毛片视频| a级毛片a级免费在线| 舔av片在线| 免费观看人在逋| 亚洲欧美日韩卡通动漫| 性插视频无遮挡在线免费观看| 观看美女的网站| 一级毛片久久久久久久久女| 青青草视频在线视频观看| 国产一区二区三区在线臀色熟女| 免费人成视频x8x8入口观看| 亚洲高清免费不卡视频| 午夜精品在线福利| 毛片女人毛片| 欧美色欧美亚洲另类二区| 久久精品国产鲁丝片午夜精品| 久久人妻av系列| 国产亚洲av片在线观看秒播厂 | 国产探花极品一区二区| 在线观看免费视频日本深夜| 我的老师免费观看完整版| 熟女电影av网| 亚洲真实伦在线观看| 日韩人妻高清精品专区| 国内精品美女久久久久久| 国产真实伦视频高清在线观看| 国产精品国产三级国产av玫瑰| 男人舔女人下体高潮全视频| 亚洲美女视频黄频| 国产伦精品一区二区三区视频9| 美女国产视频在线观看| 美女xxoo啪啪120秒动态图| 噜噜噜噜噜久久久久久91| 欧美性猛交黑人性爽| 久久精品影院6| 亚洲内射少妇av| 床上黄色一级片| 欧美zozozo另类| 最后的刺客免费高清国语| 好男人视频免费观看在线| 毛片女人毛片| 春色校园在线视频观看| 精品欧美国产一区二区三| 亚洲国产精品成人久久小说 | 看十八女毛片水多多多| 色哟哟哟哟哟哟| 在线观看一区二区三区| 久久99热这里只有精品18| 免费不卡的大黄色大毛片视频在线观看 | 国产精品精品国产色婷婷| 麻豆国产av国片精品| 亚洲人成网站高清观看| 午夜福利高清视频| 久久国产乱子免费精品| 久久6这里有精品| 一区福利在线观看| 男人舔女人下体高潮全视频| 午夜a级毛片| 麻豆精品久久久久久蜜桃| 国产中年淑女户外野战色| 免费观看的影片在线观看| 国产片特级美女逼逼视频| 男人和女人高潮做爰伦理| 一进一出抽搐动态| 午夜爱爱视频在线播放| 国产精品久久电影中文字幕| 99热这里只有是精品50| 欧美最黄视频在线播放免费| 少妇的逼水好多| 国产伦在线观看视频一区| 国产精品一区二区在线观看99 | 日产精品乱码卡一卡2卡三| 亚洲成人中文字幕在线播放| 黄片wwwwww| 99热6这里只有精品| 午夜精品国产一区二区电影 | 午夜福利视频1000在线观看| 成人av在线播放网站| 精品久久久久久成人av| 免费人成在线观看视频色| 久久亚洲国产成人精品v| 国产午夜福利久久久久久| 久久欧美精品欧美久久欧美| 日韩制服骚丝袜av| 丰满人妻一区二区三区视频av| 久久精品综合一区二区三区| 国产精品久久久久久av不卡| 国产亚洲精品久久久com| 久久久久网色| 啦啦啦啦在线视频资源| 国产激情偷乱视频一区二区| 久久久久久久久久久免费av| 亚洲国产高清在线一区二区三| 国产极品天堂在线| 亚洲欧洲日产国产| 日本一二三区视频观看| 欧美成人免费av一区二区三区| 啦啦啦观看免费观看视频高清| www.av在线官网国产| 久久欧美精品欧美久久欧美| 麻豆国产97在线/欧美| 国产黄片视频在线免费观看| 欧美色视频一区免费| 亚洲三级黄色毛片| 久久热精品热| 99九九线精品视频在线观看视频| 青春草国产在线视频 | 99视频精品全部免费 在线| 久久精品国产清高在天天线| 午夜亚洲福利在线播放| 99热6这里只有精品| 国产精品麻豆人妻色哟哟久久 | 性欧美人与动物交配| 不卡视频在线观看欧美| 亚洲精品成人久久久久久| 毛片一级片免费看久久久久| 精品99又大又爽又粗少妇毛片| 色综合色国产| 日韩制服骚丝袜av| 一本久久中文字幕| 国产又黄又爽又无遮挡在线| 人体艺术视频欧美日本| 九九久久精品国产亚洲av麻豆| 禁无遮挡网站| 99久久人妻综合| 一进一出抽搐gif免费好疼| 欧美性感艳星| 人妻久久中文字幕网| 小说图片视频综合网站| 亚洲av男天堂| 国产中年淑女户外野战色| 国产成人freesex在线| 久久婷婷人人爽人人干人人爱| 成人性生交大片免费视频hd| 少妇的逼水好多| 欧美日韩精品成人综合77777| 国产精品日韩av在线免费观看| 免费看av在线观看网站| 成人亚洲欧美一区二区av| 国产成年人精品一区二区| 韩国av在线不卡| 看片在线看免费视频| 国产女主播在线喷水免费视频网站 | 给我免费播放毛片高清在线观看| 国产午夜精品论理片| 丰满人妻一区二区三区视频av| 国产精品一二三区在线看| 国产亚洲精品av在线| 精品国产三级普通话版| 最近中文字幕高清免费大全6| 国内久久婷婷六月综合欲色啪| 国产高潮美女av| 免费大片18禁| 日本三级黄在线观看| 女同久久另类99精品国产91| 精品久久国产蜜桃| 老司机福利观看| 亚洲成人久久性| av视频在线观看入口| 亚洲在线观看片| 亚洲乱码一区二区免费版| 一边摸一边抽搐一进一小说| 国产精品人妻久久久久久| 国产精品福利在线免费观看| 能在线免费看毛片的网站| 久久久久免费精品人妻一区二区| 国产成人91sexporn| 亚洲成a人片在线一区二区| 国产白丝娇喘喷水9色精品| 午夜a级毛片| av女优亚洲男人天堂| 亚洲性久久影院| 国产蜜桃级精品一区二区三区| 国产亚洲精品久久久com| 亚洲无线观看免费| 97人妻精品一区二区三区麻豆| 精品久久久久久久末码| 国产色爽女视频免费观看| 麻豆久久精品国产亚洲av| 91久久精品国产一区二区三区| 女人被狂操c到高潮| 在线a可以看的网站| 久久人人爽人人爽人人片va| 午夜激情福利司机影院| 国产精品不卡视频一区二区| av免费观看日本| 美女内射精品一级片tv| 国产亚洲av片在线观看秒播厂 | 精品免费久久久久久久清纯| 久久久a久久爽久久v久久| 亚洲欧美日韩无卡精品| 午夜精品在线福利| 国产精品乱码一区二三区的特点| 亚洲国产欧洲综合997久久,| 精品99又大又爽又粗少妇毛片| 91午夜精品亚洲一区二区三区| 卡戴珊不雅视频在线播放| 大又大粗又爽又黄少妇毛片口| 午夜福利在线在线| 国产精品综合久久久久久久免费| 亚洲综合色惰| 国产麻豆成人av免费视频| 网址你懂的国产日韩在线| 亚洲国产欧洲综合997久久,| 欧美3d第一页| 国产亚洲精品久久久久久毛片| 亚洲国产欧洲综合997久久,| 国产视频内射| 国产探花极品一区二区| 久久婷婷人人爽人人干人人爱| 乱系列少妇在线播放| 久久久久国产网址| 国产精品人妻久久久影院| 最近2019中文字幕mv第一页| 精品少妇黑人巨大在线播放 | 午夜a级毛片| 最近2019中文字幕mv第一页| 亚洲经典国产精华液单| 91久久精品国产一区二区三区| 欧美最黄视频在线播放免费| 丰满乱子伦码专区| 成人午夜精彩视频在线观看| 国国产精品蜜臀av免费| 18+在线观看网站| 久久久精品欧美日韩精品| 18+在线观看网站| 国产精品日韩av在线免费观看| 亚洲av中文字字幕乱码综合| av在线天堂中文字幕| kizo精华| 日本爱情动作片www.在线观看| 亚洲欧美日韩东京热| 午夜亚洲福利在线播放| 亚洲色图av天堂| 久久韩国三级中文字幕| 九九在线视频观看精品| 好男人在线观看高清免费视频| 亚洲中文字幕一区二区三区有码在线看| 久久九九热精品免费| 色哟哟哟哟哟哟| 国产精品久久久久久精品电影小说 | 三级国产精品欧美在线观看| 乱系列少妇在线播放| 日本免费一区二区三区高清不卡| 国产中年淑女户外野战色| 蜜臀久久99精品久久宅男| 色视频www国产| 国产一区二区三区av在线 | 亚洲av.av天堂| 国产 一区 欧美 日韩| 美女黄网站色视频| 亚洲丝袜综合中文字幕| 欧美性猛交╳xxx乱大交人| 伦精品一区二区三区| 天美传媒精品一区二区| 国产精品.久久久| videossex国产| 全区人妻精品视频| 能在线免费观看的黄片| 2021天堂中文幕一二区在线观| 干丝袜人妻中文字幕| av免费观看日本| 成年av动漫网址| 少妇被粗大猛烈的视频| 乱码一卡2卡4卡精品| av又黄又爽大尺度在线免费看 | 特级一级黄色大片| 91在线精品国自产拍蜜月| 国产黄片美女视频| 99久国产av精品国产电影| 一区二区三区高清视频在线| 久久精品国产亚洲网站| 3wmmmm亚洲av在线观看| 精品人妻一区二区三区麻豆| 亚洲熟妇中文字幕五十中出| 国产成人午夜福利电影在线观看| 69人妻影院| 看十八女毛片水多多多| 给我免费播放毛片高清在线观看| 又粗又爽又猛毛片免费看| 国产日本99.免费观看| 三级国产精品欧美在线观看| 亚洲五月天丁香| 国产伦一二天堂av在线观看| 久久久久久久午夜电影| av在线老鸭窝| 婷婷色综合大香蕉| 99热全是精品| 国产午夜福利久久久久久| 亚洲欧美日韩卡通动漫| 97超碰精品成人国产| 国产精品福利在线免费观看| 亚洲久久久久久中文字幕| 九九久久精品国产亚洲av麻豆| 国产色爽女视频免费观看| 九色成人免费人妻av| 小蜜桃在线观看免费完整版高清| 日本av手机在线免费观看| 国语自产精品视频在线第100页| 高清毛片免费观看视频网站| 1000部很黄的大片| 日本黄色片子视频| 热99在线观看视频| 欧美一级a爱片免费观看看| 九九热线精品视视频播放| 麻豆成人av视频| 丝袜美腿在线中文| 免费观看在线日韩| 久久精品国产亚洲网站| 日韩精品有码人妻一区| 国产成人精品一,二区 | 久久鲁丝午夜福利片| 18禁在线无遮挡免费观看视频| 亚洲中文字幕一区二区三区有码在线看| 亚洲人与动物交配视频| 日韩中字成人| 最近最新中文字幕大全电影3| 长腿黑丝高跟| 国产精品99久久久久久久久| 日韩欧美精品免费久久| 国内少妇人妻偷人精品xxx网站| 国产日本99.免费观看| 色噜噜av男人的天堂激情| 国产 一区 欧美 日韩| 成年女人永久免费观看视频| 亚洲最大成人中文| 日韩欧美精品免费久久| 亚洲激情五月婷婷啪啪| 午夜福利成人在线免费观看| 精品午夜福利在线看| 在线播放无遮挡| 亚洲aⅴ乱码一区二区在线播放| 在线观看免费视频日本深夜| 午夜激情欧美在线| 亚洲精品日韩av片在线观看| 午夜久久久久精精品| 少妇的逼好多水| 熟妇人妻久久中文字幕3abv| 国产成人影院久久av| 久久热精品热| 校园春色视频在线观看| 国产亚洲91精品色在线| 一个人看视频在线观看www免费| 午夜老司机福利剧场| 欧美一区二区精品小视频在线| 婷婷色av中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 国产精品嫩草影院av在线观看| av天堂在线播放| 亚洲国产精品久久男人天堂| 色哟哟·www| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品永久免费网站| АⅤ资源中文在线天堂| 搞女人的毛片| 毛片女人毛片| 91aial.com中文字幕在线观看| 又爽又黄a免费视频| 97在线视频观看| 国产毛片a区久久久久| 在现免费观看毛片| 亚洲激情五月婷婷啪啪| 国产久久久一区二区三区| 国产成人福利小说| 国产毛片a区久久久久| 精品欧美国产一区二区三| 91aial.com中文字幕在线观看| av黄色大香蕉| 青春草亚洲视频在线观看| 变态另类成人亚洲欧美熟女| 三级毛片av免费| 久久这里只有精品中国| 欧美潮喷喷水| 亚洲国产精品成人综合色| 国产高潮美女av| 欧美成人免费av一区二区三区| 亚洲熟妇中文字幕五十中出| 免费看光身美女| 在现免费观看毛片| 18禁黄网站禁片免费观看直播| 亚洲天堂国产精品一区在线| 成人漫画全彩无遮挡| 人体艺术视频欧美日本| 日本撒尿小便嘘嘘汇集6| 99热网站在线观看| 亚洲国产日韩欧美精品在线观看| 最后的刺客免费高清国语| 午夜福利在线在线| 中文在线观看免费www的网站| 久久久精品欧美日韩精品| 在线观看免费视频日本深夜| 国产视频内射| 97超视频在线观看视频| 精品少妇黑人巨大在线播放 | 大又大粗又爽又黄少妇毛片口| 亚洲av第一区精品v没综合| 亚洲中文字幕一区二区三区有码在线看| 午夜激情欧美在线| 丰满的人妻完整版| 国产伦理片在线播放av一区 | 成年版毛片免费区| 狂野欧美白嫩少妇大欣赏| 国产蜜桃级精品一区二区三区| 国产精品人妻久久久久久| 最后的刺客免费高清国语| 人人妻人人澡欧美一区二区| 久久人妻av系列| 一个人免费在线观看电影| 国产成人精品久久久久久| 国产精品蜜桃在线观看 | 99久久成人亚洲精品观看| 少妇的逼好多水| 国产白丝娇喘喷水9色精品| 国产麻豆成人av免费视频| www.色视频.com| 色5月婷婷丁香| 午夜精品在线福利| 看非洲黑人一级黄片| a级一级毛片免费在线观看| 国产精品久久视频播放| 自拍偷自拍亚洲精品老妇| 亚洲av成人精品一区久久| 日日干狠狠操夜夜爽| 午夜福利高清视频| 国产av一区在线观看免费| 18禁在线播放成人免费| 亚洲一区二区三区色噜噜| 亚洲av电影不卡..在线观看| 色哟哟·www| 国产高清不卡午夜福利| 亚洲一区高清亚洲精品| 中文精品一卡2卡3卡4更新| 99久久九九国产精品国产免费| 日韩 亚洲 欧美在线| 精品一区二区免费观看| 成人永久免费在线观看视频| 国产片特级美女逼逼视频| 亚洲欧美成人综合另类久久久 | 成人毛片60女人毛片免费| 午夜精品国产一区二区电影 | 亚洲欧美日韩高清在线视频| 免费人成视频x8x8入口观看| 2022亚洲国产成人精品| 99久久人妻综合| 精品午夜福利在线看| 国产在线男女| 欧美zozozo另类| 国产老妇伦熟女老妇高清| 永久网站在线| 韩国av在线不卡| 国产 一区精品| 国产av一区在线观看免费| 国产精品三级大全| 国产成人一区二区在线| 亚洲精品影视一区二区三区av| 欧美精品一区二区大全| av女优亚洲男人天堂| 久久精品国产自在天天线| 久久精品国产鲁丝片午夜精品| 一边亲一边摸免费视频| 欧美+日韩+精品| 能在线免费看毛片的网站| ponron亚洲| 中文字幕免费在线视频6| 国产精品免费一区二区三区在线| 麻豆一二三区av精品| а√天堂www在线а√下载|