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

    V肋對尾緣劈縫氣膜冷卻特性的影響

    2021-07-07 10:18:06葉林劉存良楊寓全黃蓉朱安冬
    航空學(xué)報(bào) 2021年6期
    關(guān)鍵詞:尾緣氣膜對流

    葉林,劉存良, 2,*,楊寓全,黃蓉,朱安冬

    1. 西北工業(yè)大學(xué) 動力與能源學(xué)院,西安 710072

    2. 西北工業(yè)大學(xué) 陜西省航空動力系統(tǒng)熱科學(xué)重點(diǎn)實(shí)驗(yàn)室,西安 710129

    不加力超聲速巡航是第4代戰(zhàn)斗機(jī)最重要的技術(shù)特征之一,要實(shí)現(xiàn)這一技術(shù),在發(fā)動機(jī)方面需要提高渦輪前溫度[1]。追蹤先進(jìn)飛行器的發(fā)展歷程,第4代戰(zhàn)斗機(jī)的熱負(fù)荷相比之前會有較為明顯的增加,因此高效的熱管理系統(tǒng)設(shè)計(jì)受到越來越多的重視[2]。尾緣區(qū)域的有效冷卻是渦輪葉片冷卻中最困難的挑戰(zhàn)之一,這是由于狹小的通道中難以布置冷卻結(jié)構(gòu),且尾緣的壓力面?zhèn)群臀γ鎮(zhèn)鹊臒嶝?fù)荷都很高。因此為了保證尾緣的結(jié)構(gòu)完整性在發(fā)動機(jī)服役期間內(nèi)不被高溫燃?xì)鉄g破壞,必須對其進(jìn)行高效冷卻。

    Cunha等[3-4]最先對比了全縫、離散孔、劈縫3種冷卻結(jié)構(gòu)在尾緣有限空間內(nèi)的冷卻效果,基于一維方法建立了不同尾緣冷卻設(shè)計(jì)方案的溫度分布解析關(guān)系式,結(jié)果表明相同工況下所計(jì)算出的一維溫度分布中劈縫結(jié)構(gòu)的壁溫最低。劈縫冷卻結(jié)構(gòu)是通過切除葉片尾緣壓力面?zhèn)鹊牟糠直诿?,且保留吸力面?zhèn)鹊谋诿娑纬傻?,壁面間的若干分隔肋所形成的切向出流縫使得冷氣從縫口噴射出,在劈縫表面上覆蓋形成冷氣膜。相對于其余2種冷卻結(jié)構(gòu),劈縫結(jié)構(gòu)在確保結(jié)構(gòu)完整性和空氣動力學(xué)的需求的同時(shí),可有效增強(qiáng)其冷卻性能。

    目前有關(guān)結(jié)構(gòu)參數(shù)對尾緣劈縫氣膜冷卻特性和機(jī)理的研究主要集中在兩方面,一是劈縫上游強(qiáng)化換熱結(jié)構(gòu)的影響。Martini等[5]研究了上游擾流柱對劈縫氣膜冷卻效率的影響,他們發(fā)現(xiàn)劈縫出口的氣膜核心區(qū)的長度與上游結(jié)構(gòu)有關(guān),增加擾流柱排并未提高氣膜冷卻效率,反而提高了氣膜的衰弱速率,沒有擾流柱排的延伸段起到了均化氣流的作用,可以減小湍流的摻混;Wu等[6]使用瞬態(tài)熱色液晶測量技術(shù)研究了尾緣劈縫上游帶有直肋和擾流柱的內(nèi)部通道的傳熱特性和壓降系數(shù),Shen等[7]將凹窩和突脊結(jié)構(gòu)作為內(nèi)部冷卻通道的強(qiáng)化換熱裝置,采用瞬態(tài)紅外熱成像和穩(wěn)態(tài)數(shù)值模擬考察了尾緣的狹窄梯形通道的傳熱特性及機(jī)理,但文獻(xiàn)[6-7]中并未關(guān)注劈縫表面的冷卻性能。

    影響尾緣劈縫氣膜冷卻特性的另一結(jié)構(gòu)方面是外部結(jié)構(gòu)參數(shù),比如唇板造型、劈縫間隔、分隔肋形狀等。Martini等[5]考察了唇板厚度對劈縫氣膜冷卻效率的影響,與Mukherjee[8]的試驗(yàn)數(shù)據(jù)對比表明了氣膜的衰減過程主要受控于唇板厚度,薄唇板主導(dǎo)氣膜冷卻效率會顯著增大,而厚唇板的非定常脫落渦是導(dǎo)致劈縫表面氣膜冷卻效率下降的一個(gè)重要因素;Horbach等[9-10]采用紅外熱成像技術(shù)研究了唇板造型及厚度對劈縫氣膜冷卻特性的影響,結(jié)果表明未經(jīng)倒圓的唇板導(dǎo)致劈縫出口下游的氣流分離會產(chǎn)生較大的渦旋,倒圓造型為氣流提供了導(dǎo)引,降低脫落渦的強(qiáng)度,使得內(nèi)側(cè)倒圓唇板結(jié)構(gòu)的氣膜冷卻效率較大;魏建生等[11]采用瞬態(tài)熱色液晶測量技術(shù)研究了不同加工造型的收縮型分隔肋對劈縫氣膜冷卻特性的影響;王瑞琴等[12]數(shù)值模擬對比了有無分隔肋下的劈縫區(qū)域的流動傳熱性能。

    以上研究有助于說明和揭示尾緣劈縫的流動特性及機(jī)理,但目前仍缺乏提高尾緣冷卻性能的具體建議。與葉片其他區(qū)域氣膜冷卻相比,尾緣劈縫結(jié)構(gòu)的壓力面?zhèn)群臀γ鎮(zhèn)韧瑫r(shí)受到燃?xì)鉄崃鞯挠绊懀p出口噴射形成的氣膜不僅覆蓋在劈縫表面從而保護(hù)壓力面,還需充當(dāng)對流散熱器以將吸力面?zhèn)葌飨驂毫γ鎮(zhèn)鹊臒崃繋ё?。因此尾緣區(qū)域的冷卻設(shè)計(jì)追求的目標(biāo)是:高氣膜有效性和強(qiáng)對流換熱性。迄今為止,通過改變唇板與分隔肋結(jié)構(gòu)使得尾緣劈縫的氣膜貼壁性極好[13-14],氣膜冷卻效率的提升空間十分有限。因此,研究具有更強(qiáng)換熱性以達(dá)到更高的綜合冷卻性能的尾緣劈縫冷卻結(jié)構(gòu)是極其必需的。

    從已公開發(fā)表的文獻(xiàn)來看,突脊擾流結(jié)構(gòu)的研究通常著眼于熱端部件的內(nèi)冷通道[15]或是氣膜孔附近[16-17],V肋作為擾流強(qiáng)化換熱的一種基礎(chǔ)肋型,其表現(xiàn)出了較直肋和斜肋更佳的傳熱性能,目前并未應(yīng)用于尾緣劈縫表面。因此,本文突破單純增強(qiáng)劈縫表面氣膜冷卻效率的傳統(tǒng)做法,提出了在劈縫表面布置V肋型擾流結(jié)構(gòu)以增強(qiáng)劈縫表面對流換熱系數(shù)來提升尾緣綜合冷卻效果的新思路,設(shè)計(jì)了一種帶有V肋的尾緣擾流劈縫冷卻結(jié)構(gòu)。

    本文報(bào)道了一項(xiàng)試驗(yàn)研究以探究V肋對尾緣劈縫表面氣膜冷卻特性的影響,采用壓力敏感漆(Pressure-Sensitive Paint, PSP)技術(shù)和瞬態(tài)熱色液晶(Transient thermochromic Liquid-Crystal, TLC)技術(shù)分別獲得了劈縫表面的氣膜冷卻效率和對流換熱系數(shù)分布,詳細(xì)對比分析了吹風(fēng)比及V肋寬度對2種不同尾緣劈縫表面氣膜冷卻特性的影響,并采用凈熱流密度值評估對比了帶有V肋的劈縫結(jié)構(gòu)的綜合冷卻性能。

    1 試驗(yàn)系統(tǒng)及測量方法

    1.1 試驗(yàn)臺系統(tǒng)

    尾緣劈縫冷卻平板試驗(yàn)在低速風(fēng)洞中進(jìn)行,低速環(huán)境下的尾緣冷卻試驗(yàn)數(shù)據(jù)已由Krueckels等[18]驗(yàn)證過,與高速葉柵試驗(yàn)數(shù)據(jù)近乎一致,可應(yīng)用于真實(shí)渦輪葉片。試驗(yàn)系統(tǒng)如圖1所示,其主要由主流、空氣二次流和氮?dú)舛瘟鹘M成。主流由離心風(fēng)機(jī)提供,冷卻穩(wěn)定后進(jìn)入試驗(yàn)臺的主流段,主流穩(wěn)壓段與試驗(yàn)測量段之前由收縮腔過渡連接,兩段收縮腔中夾有一個(gè)絲網(wǎng)快速加熱器,測量段的截面尺寸為220 mm×80 mm。空氣二次流和氮?dú)舛瘟鞣謩e由螺桿式壓縮機(jī)和高壓氮?dú)夤藿?jīng)減壓閥降壓后提供,二次流氣體流過質(zhì)量流量計(jì)和加熱罐后進(jìn)入劈縫冷氣(文中的“冷氣”代指二次流)通道,再從尾緣劈縫縫口射出。

    圖1 試驗(yàn)系統(tǒng)示意圖

    試驗(yàn)采集系統(tǒng)中,主流速度Ug和溫度Tg分別由布置在試驗(yàn)段中的皮托管和4根熱電偶獲得;尾緣劈縫測量板在展向設(shè)置了5個(gè)劈縫表面以形成周期性,展向兩側(cè)的縫出口處均布有熱電偶以測量劈縫出口的二次流溫度Tc;采用圖像采集系統(tǒng)獲取中間3個(gè)劈縫表面的數(shù)據(jù)以檢驗(yàn)周期性。試驗(yàn)過程中,將發(fā)射光波長范圍為390~410 nm 的藍(lán)紫光激發(fā)光源和F2000 Scientific相機(jī)應(yīng)用于PSP技術(shù),將泛光燈和型號為HDR-CX350的索尼相機(jī)應(yīng)用于TLC技術(shù)。

    1.2 試驗(yàn)件

    由于渦輪葉片尾緣實(shí)物尺寸很小,劈縫出口高度約0.5 mm,本文按照相似原理[19]將實(shí)物模型化,在原始結(jié)構(gòu)的基礎(chǔ)上簡化并放大14倍的劈縫模型由唇板(頂部)、多條分隔肋(中間)和測量板(底部)組成,如圖2所示。其中,唇板的上表面(其代表了尾緣的壓力面)和分隔肋上表面在同一水平面以接近真實(shí)渦輪葉片的尾緣結(jié)構(gòu)。沿展向中心線方向截出的A-A視圖展示了劈縫模型的主要參數(shù)尺寸,主流從左側(cè)流入,二次流由唇板與測量板的中間流道進(jìn)入高度H為7 mm的劈縫冷氣通道,經(jīng)延伸長度為14H、寬為2H的矩形肋將進(jìn)入冷卻通道的流體分開,以形成截面尺寸為4H×H的獨(dú)立二次流通道,二次流的噴射角為14°,唇板形狀為內(nèi)側(cè)倒圓結(jié)構(gòu),唇板厚度等于通道高度H,其結(jié)構(gòu)選取原則根據(jù)文獻(xiàn)[9-10]。為滿足瞬態(tài)換熱試驗(yàn)的一維半無限大假設(shè)理論,設(shè)計(jì)測量板厚度為35 mm。

    圖2 尾緣劈縫模型示意圖

    劈縫表面為主流與二次流摻混區(qū)域,是本研究的主要對象。縫出口外分隔肋的展向?qū)挾扔?H變化至不同數(shù)值以形成中心軸對稱的不同形狀的劈縫表面。如圖3所示,半劈縫表面的流向長度為8H,兩邊的展向?qū)挾确謩e為4H和S。本文研究的結(jié)構(gòu)包含2種劈縫表面形狀:一種為未擴(kuò)張型劈縫表面,即S=4H;一種為分隔肋的展向?qū)挾扔?H收縮至0所形成的擴(kuò)張型劈縫表面,即S=6H。 所研究的帶有V肋的擾流劈縫結(jié)構(gòu)是在劈縫表面上布置2條V肋,靠近縫出口的V肋起始端距劈縫出口為0.6H,肋高為0.2H,肋間距p=1.6H,2種肋寬e=0.4H, 0.2H。本文共研究了6種結(jié)構(gòu),如表1所示,編號為CASE#1~CASE#6,尾緣半劈縫試驗(yàn)件如圖4所示。

    圖3 原始及帶有V肋的劈縫表面示意圖

    表1 結(jié)構(gòu)參數(shù)

    圖4 尾緣劈縫試驗(yàn)件實(shí)物圖片

    1.3 測量方法及原理

    PSP技術(shù)最初被應(yīng)用于壁面靜壓的測量,根據(jù)傳熱傳質(zhì)類比原理又被推廣為氣膜冷卻研究中氣膜冷卻效率的測量[20],PSP發(fā)射光的強(qiáng)度是隨涂料表面的氧氣分壓遞減的函數(shù),可將修正背景噪聲后的光強(qiáng)比和氧氣分壓比數(shù)據(jù)進(jìn)行擬合,從而獲得兩者之間的定量關(guān)系式為

    (1)

    式中:IR和 (PO2)R分別為參考條件下(有光無風(fēng))的發(fā)射光強(qiáng)度和氧氣分壓;IB為背景條件下(無光無風(fēng))的發(fā)射光強(qiáng)度;I和PO2分別為試驗(yàn)條件下(有光有風(fēng))的發(fā)射光強(qiáng)度和氧氣分壓;系數(shù)A、B、C、D可通過標(biāo)定試驗(yàn)[20]測得,本試驗(yàn)溫度303 K條件下標(biāo)定表達(dá)式的系數(shù)分別為0.220 8、0.703 3、0.220 5和-0.137 8。

    試驗(yàn)過程中,通過調(diào)節(jié)主流水冷與二次流加熱罐將Tg和Tc的溫差均控制在0.3 K以內(nèi)。待流動狀態(tài)穩(wěn)定后,分別在主流為空氣、二次流為空氣的試驗(yàn)條件和主流為空氣、二次流為氮?dú)獾脑囼?yàn)條件獲得I;在無風(fēng)有光的參考條件和無風(fēng)無光的黑暗條件下分別獲得IR和IB,結(jié)合標(biāo)定表達(dá)式,可獲得待測表面的氣膜冷卻效率。PSP技術(shù)測得的局部氣膜冷卻效率η越大,不確定度越小。本試驗(yàn)中當(dāng)η=0.3時(shí)不確定度為2.38%,當(dāng)η=0.7時(shí)不確定度為1.36%,當(dāng)η=0.9時(shí)不確定度為1.12%。

    本文的對流換熱系數(shù)h由基于瞬態(tài)導(dǎo)熱理論的窄帶熱色液晶技術(shù)測量獲得,關(guān)于雙參數(shù)求解的瞬態(tài)傳熱測量理論的具體介紹可參考文獻(xiàn)[21-22],待測表面的壁溫Tw隨時(shí)間t的變化可根據(jù)一維半無限大理論的非穩(wěn)態(tài)導(dǎo)熱方程推導(dǎo)得到:

    [ηTc(t)+(1-η)Tg(t)-T0]

    (2)

    式中:T0為試驗(yàn)件表面的初始壁溫;a為熱擴(kuò)散系數(shù);λ為試驗(yàn)件材料的導(dǎo)熱系數(shù);erfc(·)為余誤差函數(shù),其表達(dá)式可參考文獻(xiàn)[19]。式(2)中,η和h為2個(gè)未知量。

    由于劈縫表面的氣膜冷卻效率已通過PSP技術(shù)測得,故可通過設(shè)置氣流溫度將式(2)近似為常規(guī)單參數(shù)的求解式。試驗(yàn)過程中,通過調(diào)整加熱功率,使得主流和二次流獲得相同的溫度階躍,即Tg(t)=Tc(t),此時(shí)僅需再明確某一時(shí)刻下的壁溫Tw(t)和變色時(shí)間t即可通過一次試驗(yàn)求解出對流換熱系數(shù),本研究中已通過多次試驗(yàn)驗(yàn)證了試驗(yàn)重復(fù)性。TLC技術(shù)的不確定度[23]主要源于時(shí)間、溫度及材料物性,通過隨機(jī)誤差傳遞進(jìn)行分析可得本試驗(yàn)中換熱系數(shù)的不確定度為10%~15%。

    2 結(jié)果分析與討論

    2.1 試驗(yàn)工況及數(shù)據(jù)處理

    本試驗(yàn)中設(shè)置入口速度為23.16 m/s以確保主流雷諾數(shù)Reg為10 000,通過調(diào)節(jié)二次流的質(zhì)量流量mc使得吹風(fēng)比BR在0.25~2.00范圍內(nèi)變化,Reg和BR的定義分別為

    (3)

    (4)

    式中:Uc為縫出口射流的速度;ρg和ρc分別為主流和二次流的密度;μg為主流的動力黏度;Ac為劈縫出口的面積。

    由于分隔肋影響下所產(chǎn)生的二次流渦流以及主流和二次流流經(jīng)劈縫出口(唇板)造成的渦脫落等物理過程使得劈縫出口下游的氣流具有很高的非定常性質(zhì),這些渦流使得其表面的對流換熱特性有別于一般平板結(jié)構(gòu)。因此本文用無量綱參數(shù)Nu/Nu∞表達(dá)劈縫表面的對流換熱強(qiáng)化性能,其中Nu為劈縫表面的努塞爾數(shù),Nu∞參考Kays和Crawford[24]的充分發(fā)展段管內(nèi)湍流換熱經(jīng)驗(yàn)表達(dá)式描述:

    (5)

    式中:k為流體的導(dǎo)熱系數(shù);Prc為二次流通道的冷氣普朗特?cái)?shù);Rec為二次流通道的冷氣雷諾數(shù),本文選取BR=1.00時(shí)的Rec。

    采用基于恒定綜合冷卻效率φ=0.6計(jì)算的凈熱流密度值NHFR來評估V肋對尾緣劈縫綜合冷卻性能的影響,數(shù)值越高表明整體冷卻性能越好。由于尾緣劈縫壁面的熱流量是由劈縫表面?zhèn)飨蛏淞?,因此NHFR值大于1.0。

    (6)

    式中:q為有冷氣射流時(shí)劈縫表面的熱流密度;q∞是基于式(5)的Nu∞計(jì)算得到的。

    2.2 原始尾緣劈縫表面的氣膜冷卻特性

    圖5展示了小吹風(fēng)比時(shí)的原始尾緣劈縫表面(CASE#1、CASE#2)的氣膜冷卻效率η分布,同一吹風(fēng)比時(shí),擴(kuò)張型的劈縫表面(CASE#2)的氣膜覆蓋性弱于未擴(kuò)張型結(jié)構(gòu)(CASE#1),劈縫出口的核心區(qū)的η接近于1,且其長度受劈縫表面擴(kuò)張程度的影響,這是由于擴(kuò)張較大的結(jié)構(gòu)使得冷氣射流的流向動量降低,同樣的原因顯著影響了下游區(qū)域的氣膜覆蓋性。當(dāng)BR=0.25時(shí),CASE#1下游X/H=6.0~8.0區(qū)域的低η的分布較為均勻,數(shù)值約為0.7,而CASE#2下游呈現(xiàn)出明顯的“雙駝峰”式低η分布,數(shù)值約為0.4~0.5, 且范圍較大,該現(xiàn)象在文獻(xiàn)[14]中也出現(xiàn)。隨著吹風(fēng)比的增大,冷氣流量隨之增大,尾緣劈縫下游的氣膜覆蓋范圍增大,且更加穩(wěn)定,高η核心區(qū)逐漸擴(kuò)大。在BR=0.50時(shí),CASE#1的氣膜已幾乎完全覆蓋劈縫表面,更大的吹風(fēng)比工況下的分布此處不再贅述。而對于CASE#2,遠(yuǎn)下游區(qū)域的中心位置處的氣膜覆蓋性相對較好,而兩側(cè)的低η值在0.7~0.8范圍內(nèi),隨著吹風(fēng)比的繼續(xù)增大,低η值雖有所改善,但始終未能完全被核心氣膜覆蓋,BR=1.00與BR=0.75時(shí)相比,低η區(qū)域增大,這是由于在該種結(jié)構(gòu)下縫出口渦流的脫落隨著吹風(fēng)比的提高而有所加強(qiáng),使得遠(yuǎn)下游區(qū)域的主流與冷氣的摻混加劇,類似的現(xiàn)象在文獻(xiàn)[25]有所報(bào)道。

    圖5 原始尾緣劈縫表面的氣膜冷卻效率分布云圖

    由于原始劈縫表面對流換熱系數(shù)隨著吹風(fēng)比的增大單調(diào)增大[23],因此本文選取了 CASE#1和CASE#2在3種典型吹風(fēng)比工況下的對流換熱強(qiáng)化性能Nu/Nu∞分布云圖進(jìn)行展示對比,如圖6所示。在BR=0.50時(shí),X/H<3范圍的對流換熱強(qiáng)度很弱,且在中心線附近存在一個(gè)極其低的換熱區(qū),這是由于劈縫出口處的壁面附近的渦旋引起壁面附近流體的速度梯度降低,從而使得該處換熱降低。下游區(qū)域展向兩側(cè)的低換熱區(qū)的面積逐漸減小,中心區(qū)域的換熱較為均勻,CASE#2的遠(yuǎn)下游區(qū)域隱約出現(xiàn)了“雙駝峰”狀的高換熱區(qū),這是由于該處的氣膜層受到主流脫落渦的擾動加劇,使得近壁面流體的速度梯度提升,進(jìn)而換熱強(qiáng)度增大;而該工況下的CASE#1的遠(yuǎn)下游氣膜覆蓋性很好,主流未能擾動到該處的流體,CASE#1并未出現(xiàn)高換熱區(qū)。隨著吹風(fēng)比的增大,核心低換熱區(qū)和下游展向兩側(cè)的低換熱區(qū)的面積均逐漸減小,換熱分布更為均勻。CASE#1和CASE#2的換熱分布在規(guī)律上無明顯差異,這說明劈縫表面的形狀不是影響換熱特性的主要因素。

    圖6 原始尾緣劈縫表面的對流換熱強(qiáng)化性能分布云圖

    圖7展示了吹風(fēng)比為0.50和1.00時(shí)劈縫表面的展向換熱強(qiáng)度,BR=0.50時(shí)X/H<2區(qū)域換熱相對較低,隨后急劇上升,CASE#2的換熱較強(qiáng),且與CASE#1的差距逐漸拉大,這是由于在小吹風(fēng)比工況下,CASE#2遠(yuǎn)下游區(qū)域的主流與冷氣摻混較為劇烈,增強(qiáng)了壁面附近的擾動;在BR=1.00時(shí),換熱系數(shù)的變化較為平緩,兩者相差較為穩(wěn)定,此時(shí)主流對壁面的擾動影響并不強(qiáng)烈。與Murata等[26-27]的試驗(yàn)數(shù)據(jù)對比,整體變化趨勢與換熱分布吻合很好,驗(yàn)證了本文換熱試驗(yàn)的合理性。

    圖7 本文實(shí)驗(yàn)得到的展向平均努塞爾數(shù)比值Nu/Nu∞與文獻(xiàn)數(shù)據(jù)的對比

    2.3 帶有V肋的尾緣劈縫表面的氣膜冷卻效率

    圖8展示了小吹風(fēng)比工況下帶有V肋的尾緣劈縫表面的氣膜冷卻效率η分布,在BR=0.25 時(shí), 2條V肋對下游X/H<3范圍內(nèi)的氣膜覆蓋基本沒有影響,由于V肋肋條的導(dǎo)流作用,使得在壁面附近的部分流體產(chǎn)生了沿肋條方向的分速度,其對主流與V肋間及第2條V肋的末端附近射流有著微弱的擾動,導(dǎo)致該處的η相對較低;沿著中心線Y/H=0的射流能保持較好的速度向下游發(fā)展,從而導(dǎo)致V肋的夾角附近有較好的冷流覆蓋;由于凸起型擾流結(jié)構(gòu)顯著影響了遠(yuǎn)下游區(qū)域的主流與冷氣的摻混,X/H=6.0~8.0區(qū)域出現(xiàn)了明顯的“雙駝峰”狀的低η分布,但其面積和幅度受劈縫表面的擴(kuò)張程度和V肋寬度影響。對比不同V肋寬度的結(jié)構(gòu),窄肋結(jié)構(gòu)比寬肋結(jié)構(gòu)的低η范圍明顯較小。

    圖8 小吹風(fēng)比工況下帶有V肋的劈縫表面的氣膜冷卻效率分布云圖

    冷流流量增大至BR=0.50時(shí),氣膜已能覆蓋大部分區(qū)域,且更加穩(wěn)定,僅有遠(yuǎn)下游區(qū)域未能完全覆蓋,對于未擴(kuò)張的劈縫表面(CASE#3和CASE#4),此工況下V肋寬度的影響較小,而對比CASE#5和CASE#6可以看出,帶有寬V肋的擴(kuò)張型劈縫表面的低η范圍明顯大于窄V肋結(jié)構(gòu);隨著吹風(fēng)比繼續(xù)增大到0.75,低η分布依然集中在遠(yuǎn)下游區(qū)域,寬V肋結(jié)構(gòu)的低η范圍略大于窄V肋結(jié)構(gòu),但數(shù)值均在0.8附近,與BR=0.50相比,未擴(kuò)張的劈縫表面結(jié)構(gòu)的下游氣膜覆蓋有著降低的趨勢,尤其是CASE#3,這是由于脫落渦在V肋的影響下改變了流動結(jié)構(gòu),使得隨著吹風(fēng)比的升高對壁面的擾動增強(qiáng)。

    圖9展示了BR=0.25和0.75時(shí)劈縫表面的展向平均氣膜冷卻效率η曲線,BR=0.25時(shí),帶有V肋的劈縫表面展向平均η數(shù)值在上游X/H<3與原始結(jié)構(gòu)相同,CASE#1的氣膜覆蓋性最好,V肋的加入使得下游展向平均η數(shù)值降低了0.1,CASE#3和CASE#4的差異并不明顯。帶有窄V肋的擴(kuò)張型表面(CASE#6)可達(dá)到與帶V肋的未擴(kuò)張型表面相近的水平,且與CASE#5在遠(yuǎn)下游區(qū)域的差異逐漸拉大至0.05,但均略高于CASE#2,這是由于小吹風(fēng)比工況下V肋的擾動作用使得主流對擴(kuò)張型劈縫表面附近射流的擾動減弱,使得氣膜有效性略有增強(qiáng)。BR=0.75時(shí),帶有V肋的劈縫表面在下游的展向平均η始終低于原始結(jié)構(gòu),帶有窄V肋的η下降得較為平緩,略占有優(yōu)勢。

    圖9 小吹風(fēng)比工況下展向平均氣膜冷卻效率曲線

    圖10為劈縫表面的面積平均氣膜冷卻效率η曲線,整體上隨著吹風(fēng)比的增大而增大,對于未擴(kuò)張型結(jié)構(gòu),V肋的加入對劈縫表面的氣膜覆蓋產(chǎn)生了不利影響,在小吹風(fēng)比工況下,V肋寬度對面積平均η無明顯影響,當(dāng)BR>0.75時(shí),寬V肋結(jié)構(gòu)的面積平均η較窄V肋結(jié)構(gòu)略低,對于CASE#3,存在較為明顯的隨吹風(fēng)比增大,面積平均η降低的變化規(guī)律。相同V肋寬度結(jié)構(gòu)下,未擴(kuò)張型結(jié)構(gòu)的面積平均η始終高于擴(kuò)張型結(jié)構(gòu),小吹風(fēng)比BR=0.25時(shí),帶窄V肋的CASE#6可達(dá)到與帶有V肋的CASE#3和CASE#4相近的面積平均η,V肋的加入對劈縫表面的η有較小的促進(jìn)作用,而隨著吹風(fēng)比的增大,這種積極影響不再,而CASE#5和CASE#6消除了因脫落渦影響導(dǎo)致的隨著吹風(fēng)比增大而η降低的現(xiàn)象,不同V肋寬度的擴(kuò)張型表面的面積平均η差距逐漸減小??偟膩碚f,大吹風(fēng)比工況下V肋對尾緣劈縫表面的氣膜冷卻效率影響不大。

    圖10 劈縫表面面積平均氣膜冷卻效率曲線

    2.4 帶有V肋的尾緣劈縫表面的對流換熱系數(shù)

    圖11選取展示了3個(gè)吹風(fēng)比(0.50,1.00和1.50)時(shí)帶有V肋的尾緣劈縫表面的對流換熱強(qiáng)化性能分布。在BR=0.50時(shí),第1條V肋上游的換熱很弱,甚至比與無擾流結(jié)構(gòu)的原始劈縫表面還要低;每個(gè)V肋肋條后均存在一條低換熱系數(shù)帶,這是由于肋后回流造成的低湍流度區(qū)域,對壁面邊界層的擾動較小,隨著劈縫表面擴(kuò)張程度的增大,V肋后的低換熱系數(shù)帶寬度無明顯變化;V肋表面為高換熱區(qū),且第2條V肋較第1條V肋的換熱系數(shù)有所提升;第2條V肋下游隱約出現(xiàn)了高換熱核心區(qū),此處是由射流流動結(jié)構(gòu)主導(dǎo),遠(yuǎn)下游區(qū)域的換熱較為均勻,且努塞爾數(shù)比值均大于1.1。隨著吹風(fēng)比增大至1.00,V肋上游的低換熱僅表現(xiàn)在劈縫出口的核心區(qū),肋后的低換熱區(qū)有了明顯減小,肋間逐漸出現(xiàn)了核心高換熱區(qū),第2條V肋肋條表面基本都為核心高換熱區(qū);其下游的高換熱區(qū)呈現(xiàn)明顯的“心”狀分布,由于擴(kuò)張劈縫表面的展向延伸性,使得CASE#5的高換熱區(qū)較CASE#3略大。在大吹風(fēng)比BR=1.50時(shí),核心高換熱區(qū)有了明顯的增大,且CASE#3與CASE#5在遠(yuǎn)下游區(qū)域的差異增大,但整體換熱規(guī)律無明顯變化??偟膩碚f,與小吹風(fēng)比工況相比,大吹風(fēng)比時(shí)的整體換熱強(qiáng)度更加均勻。

    圖12展示了吹風(fēng)比為0.50和1.50時(shí)的展向平均對流換熱強(qiáng)化性能曲線,帶有V肋的劈縫表面沿流向共有3處換熱峰值,分別位于第1條V肋尖端、第2條V肋附近及下游區(qū)域。在BR=0.50 時(shí),劈縫出口X/H<1范圍的換熱略低于原始劈縫結(jié)構(gòu),不同劈縫形狀結(jié)構(gòu)僅在遠(yuǎn)下游X/H>6出現(xiàn)微小差異,窄V肋結(jié)構(gòu)的展向平均換熱較寬V肋略低,但無明顯差異;在BR=1.50時(shí),第一峰值不再顯著,V肋后的核心高換熱區(qū)的范圍和幅度有了明顯增大,不同劈縫形狀結(jié)構(gòu)在遠(yuǎn)下游區(qū)域的差異有所增大,但V肋寬度依然對劈縫表面換熱無明顯影響。

    圖13為劈縫表面的面積平均對流換熱強(qiáng)化性能曲線,其中Nus為無擾流下的原始劈縫表面的努塞爾數(shù),寬V肋結(jié)構(gòu)略占有優(yōu)勢,但兩者差異極其不明顯。隨著吹風(fēng)比的升高,CASE#3和CASE#4的換熱系數(shù)比呈先上升后下降的趨勢,均在BR=1.25時(shí)達(dá)到了峰值,在BR<0.75時(shí),由于主流仍能擾動影響劈縫表面,此時(shí)V肋的加入有效地改變了摻混流動結(jié)構(gòu),BR>1.00時(shí),氣膜已能夠完全覆蓋劈縫表面,此時(shí)僅受射流流動結(jié)構(gòu)主導(dǎo),而當(dāng)吹風(fēng)比繼續(xù)增大至1.50和2.00,由于劈縫出口渦流結(jié)構(gòu)和V肋的共同作用,使得換熱強(qiáng)度不及BR=1.25。對于擴(kuò)張型劈縫表面結(jié)構(gòu),分別在BR=0.75和BR=1.50時(shí)達(dá)到峰值,原因是在小吹風(fēng)比工況下,遠(yuǎn)下游區(qū)域的核心換熱區(qū)受流體摻混擾動和射流共同影響,而在大吹風(fēng)比工況下僅受射流流動結(jié)構(gòu)主導(dǎo)。對比不同劈縫表面形狀可以看出,V肋在未擴(kuò)張劈縫表面結(jié)構(gòu)上展現(xiàn)出的強(qiáng)換熱性優(yōu)于擴(kuò)張型劈縫表面結(jié)構(gòu)。

    圖13 劈縫表面面積平均對流換熱強(qiáng)化性能曲線

    2.5 凈熱流密度值

    圖14展示了本文研究的6種結(jié)構(gòu)尾緣劈縫表面的面積平均凈熱流密度值NHFR曲線,帶有V肋的劈縫表面結(jié)構(gòu)的NHFR值始終遠(yuǎn)遠(yuǎn)大于原始尾緣劈縫結(jié)構(gòu),且隨著吹風(fēng)比的增大優(yōu)勢越發(fā)明顯,由于尾緣區(qū)域的熱流量由劈縫表面流向射流,熱流量的增大有效說明了帶有V肋的尾緣劈縫結(jié)構(gòu)的冷卻優(yōu)勢。

    圖14 劈縫表面面積平均凈熱流密度值曲線

    原始劈縫表面結(jié)構(gòu)(CASE#1和CASE#2)的NHFR值無明顯的差異,這說明僅改變劈縫表面形狀對尾緣冷卻的綜合冷卻性能的提升有限。當(dāng)BR≤1.00時(shí),帶有V肋的4種結(jié)構(gòu)的NHFR值相近,但需考慮到在小吹風(fēng)比工況下,劈縫表面遠(yuǎn)下游區(qū)域呈現(xiàn)“低冷效高換熱”分布,盡管NHFR值相對較高,但在工程設(shè)計(jì)中應(yīng)當(dāng)充分考慮這些因素的影響,以防止尾緣遠(yuǎn)下游區(qū)域的燒蝕,因此可在小吹風(fēng)比工況下選用帶有V肋的未擴(kuò)張的劈縫表面結(jié)構(gòu);在大吹風(fēng)比工況(BR>1.00)下V肋對尾緣劈縫表面的氣膜冷卻效率無明顯影響,此時(shí)表面對流換熱系數(shù)的大小主導(dǎo)壁面熱流密度的強(qiáng)弱,由于劈縫渦流結(jié)構(gòu)和V肋擾動的共同作用,不同劈縫表面形狀結(jié)構(gòu)的優(yōu)勢相當(dāng),因此在大吹風(fēng)比工況下應(yīng)用V肋時(shí)可無需考慮劈縫表面形狀。對比同一劈縫表面形狀結(jié)構(gòu),V肋寬度對劈縫表面的NHFR值無明顯影響。

    3 結(jié) 論

    本文使用PSP技術(shù)和TLC技術(shù)試驗(yàn)研究了V肋對尾緣劈縫表面的氣膜冷卻特性的影響,得到如下結(jié)論:

    1) 對于未擴(kuò)張型劈縫表面,V肋的加入對劈縫表面的氣膜覆蓋產(chǎn)生了不利影響,在小吹風(fēng)比工況下,V肋寬度對面積平均氣膜冷卻效率無明顯影響,而對于擴(kuò)張型劈縫表面,V肋的加入對劈縫表面的氣膜覆蓋有較小的促進(jìn)作用;相同V肋寬度結(jié)構(gòu)下,未擴(kuò)張型劈縫表面的氣膜冷卻效率始終高于擴(kuò)張型劈縫表面。

    2) V肋寬度對劈縫表面的換熱強(qiáng)度的影響不明顯,V肋在未擴(kuò)張劈縫表面結(jié)構(gòu)上展現(xiàn)出的強(qiáng)換熱性優(yōu)于擴(kuò)張型劈縫表面結(jié)構(gòu)。

    3) 相比原始尾緣劈縫結(jié)構(gòu),帶有2條V肋的尾緣劈縫冷卻結(jié)構(gòu)可有效增大6.9%~26.6%的凈熱流密度值。小吹風(fēng)比工況下,宜將V肋應(yīng)用于未擴(kuò)張的劈縫表面結(jié)構(gòu),大吹風(fēng)比工況下無需考慮劈縫表面形狀;V肋寬度對劈縫表面的凈熱流密度值無明顯影響。

    猜你喜歡
    尾緣氣膜對流
    T 型槽柱面氣膜密封穩(wěn)態(tài)性能數(shù)值計(jì)算研究
    齊口裂腹魚集群行為對流態(tài)的響應(yīng)
    基于強(qiáng)化換熱的偏斜尾緣設(shè)計(jì)
    能源工程(2021年1期)2021-04-13 02:05:46
    氣膜孔堵塞對葉片吸力面氣膜冷卻的影響
    靜葉柵上游端壁雙射流氣膜冷卻特性實(shí)驗(yàn)
    翼型湍流尾緣噪聲半經(jīng)驗(yàn)預(yù)測公式改進(jìn)
    具有尾緣襟翼的風(fēng)力機(jī)動力學(xué)建模與恒功率控制
    躲避霧霾天氣的氣膜館
    基于ANSYS的自然對流換熱系數(shù)計(jì)算方法研究
    二元驅(qū)油水界面Marangoni對流啟動殘余油機(jī)理
    99热国产这里只有精品6| 岛国视频午夜一区免费看| 久久伊人香网站| 精品久久蜜臀av无| 大香蕉久久成人网| 久久人人精品亚洲av| 久久久久国产精品人妻aⅴ院| 久久久久久久久中文| 国产人伦9x9x在线观看| 午夜成年电影在线免费观看| 最新在线观看一区二区三区| 性少妇av在线| bbb黄色大片| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美精品一区二区免费开放| 日韩av在线大香蕉| 热re99久久国产66热| 级片在线观看| 亚洲精品国产一区二区精华液| 黄片大片在线免费观看| 亚洲av熟女| 免费高清在线观看日韩| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲国产欧美日韩在线播放| 麻豆久久精品国产亚洲av | 这个男人来自地球电影免费观看| 交换朋友夫妻互换小说| 日本免费a在线| 丝袜美腿诱惑在线| 国产精品美女特级片免费视频播放器 | 久久 成人 亚洲| 老汉色∧v一级毛片| 欧美黄色片欧美黄色片| 欧美日韩视频精品一区| 亚洲情色 制服丝袜| 亚洲五月色婷婷综合| avwww免费| 欧美成人性av电影在线观看| 成人亚洲精品av一区二区 | 日韩大码丰满熟妇| 天堂影院成人在线观看| 国产亚洲av高清不卡| 久久亚洲精品不卡| 黄片大片在线免费观看| 免费观看精品视频网站| ponron亚洲| 亚洲精品美女久久久久99蜜臀| 中亚洲国语对白在线视频| 欧美丝袜亚洲另类 | 国产乱人伦免费视频| 日本欧美视频一区| 午夜亚洲福利在线播放| 免费久久久久久久精品成人欧美视频| 国产精品免费视频内射| 免费在线观看完整版高清| 视频在线观看一区二区三区| 久久久久久久久免费视频了| 午夜精品国产一区二区电影| 亚洲性夜色夜夜综合| 99久久国产精品久久久| 无限看片的www在线观看| 国产三级在线视频| 国产成人精品无人区| 欧美成人免费av一区二区三区| 韩国av一区二区三区四区| 亚洲伊人色综图| 午夜福利免费观看在线| 少妇的丰满在线观看| 超碰成人久久| 亚洲国产精品999在线| 亚洲九九香蕉| 欧美日韩亚洲国产一区二区在线观看| 一边摸一边做爽爽视频免费| www.999成人在线观看| 麻豆国产av国片精品| 人妻丰满熟妇av一区二区三区| 色哟哟哟哟哟哟| 亚洲av成人av| 国产熟女午夜一区二区三区| 亚洲中文日韩欧美视频| 亚洲第一av免费看| av网站免费在线观看视频| 亚洲国产中文字幕在线视频| 老司机在亚洲福利影院| 国产精品二区激情视频| 欧美日韩一级在线毛片| 国产免费男女视频| 亚洲全国av大片| 亚洲 欧美 日韩 在线 免费| 高清黄色对白视频在线免费看| 久久99一区二区三区| 日本a在线网址| 交换朋友夫妻互换小说| 一级作爱视频免费观看| 桃色一区二区三区在线观看| 女人精品久久久久毛片| 精品久久久久久成人av| 看免费av毛片| 丝袜美足系列| 久久久久国产精品人妻aⅴ院| 精品少妇一区二区三区视频日本电影| 精品一区二区三区视频在线观看免费 | 欧美日韩精品网址| 99久久久亚洲精品蜜臀av| 久久久久久久久久久久大奶| videosex国产| 91精品三级在线观看| 国产成年人精品一区二区 | 国产成人影院久久av| 国产黄色免费在线视频| 久久人妻熟女aⅴ| 亚洲av日韩精品久久久久久密| 午夜福利一区二区在线看| 黄色女人牲交| 亚洲熟妇熟女久久| 免费人成视频x8x8入口观看| 久久九九热精品免费| 国产亚洲欧美98| 人妻丰满熟妇av一区二区三区| 色婷婷久久久亚洲欧美| 亚洲一区二区三区欧美精品| 欧美色视频一区免费| 欧美黄色淫秽网站| 最近最新中文字幕大全电影3 | 久久人人爽av亚洲精品天堂| 欧美成人免费av一区二区三区| 欧美久久黑人一区二区| av国产精品久久久久影院| 日韩大尺度精品在线看网址 | 国产精品亚洲一级av第二区| 中文字幕精品免费在线观看视频| 9191精品国产免费久久| av电影中文网址| 看黄色毛片网站| 91九色精品人成在线观看| 757午夜福利合集在线观看| 女性被躁到高潮视频| 国产精品乱码一区二三区的特点 | 黄色毛片三级朝国网站| 又黄又粗又硬又大视频| 色婷婷久久久亚洲欧美| 国产99白浆流出| 新久久久久国产一级毛片| 丝袜美足系列| 这个男人来自地球电影免费观看| 人人妻人人澡人人看| 色婷婷av一区二区三区视频| 欧美激情高清一区二区三区| 国产精品 国内视频| 亚洲性夜色夜夜综合| 亚洲熟妇中文字幕五十中出 | 国产一区二区三区在线臀色熟女 | 欧美丝袜亚洲另类 | 国产精品永久免费网站| 美国免费a级毛片| 在线观看www视频免费| 国产在线观看jvid| 露出奶头的视频| 亚洲人成电影观看| 国产片内射在线| 美国免费a级毛片| 国产一区二区三区视频了| 亚洲自拍偷在线| 免费在线观看完整版高清| 51午夜福利影视在线观看| 男人操女人黄网站| 男女高潮啪啪啪动态图| 黄色丝袜av网址大全| 国产精品 国内视频| 一个人免费在线观看的高清视频| 岛国视频午夜一区免费看| 一个人免费在线观看的高清视频| 一级片免费观看大全| 国产精品成人在线| 日本撒尿小便嘘嘘汇集6| 欧美成人免费av一区二区三区| 丝袜美足系列| 一区二区三区激情视频| 成人国产一区最新在线观看| 日本 av在线| 国产精品av久久久久免费| 1024视频免费在线观看| 国产成人av教育| 国产成人啪精品午夜网站| 又黄又粗又硬又大视频| 亚洲第一av免费看| 51午夜福利影视在线观看| 亚洲五月婷婷丁香| 色在线成人网| 亚洲第一青青草原| 淫秽高清视频在线观看| 中国美女看黄片| 国产单亲对白刺激| 久久影院123| 咕卡用的链子| 丝袜在线中文字幕| 精品一区二区三卡| 99精国产麻豆久久婷婷| 波多野结衣一区麻豆| 成人18禁在线播放| 国产成+人综合+亚洲专区| 国产99白浆流出| 午夜福利一区二区在线看| 久久久久久久久久久久大奶| 怎么达到女性高潮| 美女午夜性视频免费| 国产成人影院久久av| 看片在线看免费视频| 成人三级黄色视频| 午夜免费激情av| 色播在线永久视频| 亚洲成人国产一区在线观看| 首页视频小说图片口味搜索| 欧美一区二区精品小视频在线| 狠狠狠狠99中文字幕| 久久久久久久久免费视频了| 一级a爱片免费观看的视频| 超色免费av| netflix在线观看网站| 欧美黑人精品巨大| 女人精品久久久久毛片| 成人特级黄色片久久久久久久| 久久久国产精品麻豆| 欧美另类亚洲清纯唯美| 女人高潮潮喷娇喘18禁视频| 一区在线观看完整版| 一本大道久久a久久精品| 亚洲第一青青草原| 亚洲国产中文字幕在线视频| 国产三级在线视频| 操出白浆在线播放| 国产视频一区二区在线看| 久久久精品欧美日韩精品| 日韩欧美在线二视频| 嫩草影视91久久| 在线观看午夜福利视频| 欧美激情久久久久久爽电影 | 成在线人永久免费视频| 我的亚洲天堂| 日韩 欧美 亚洲 中文字幕| 男人操女人黄网站| 亚洲欧美日韩另类电影网站| 亚洲欧美一区二区三区黑人| 脱女人内裤的视频| 国产真人三级小视频在线观看| 久久中文字幕人妻熟女| 国产亚洲欧美精品永久| 一边摸一边做爽爽视频免费| 男女床上黄色一级片免费看| 国产精品美女特级片免费视频播放器 | 亚洲精品一二三| 一级片免费观看大全| www.精华液| 国产精品国产高清国产av| 成人三级黄色视频| 1024视频免费在线观看| 欧美最黄视频在线播放免费 | 国产精品自产拍在线观看55亚洲| www国产在线视频色| 两性午夜刺激爽爽歪歪视频在线观看 | 国产熟女午夜一区二区三区| 老司机在亚洲福利影院| а√天堂www在线а√下载| av天堂久久9| 亚洲欧美日韩无卡精品| 日本wwww免费看| 久久久久久免费高清国产稀缺| 51午夜福利影视在线观看| av国产精品久久久久影院| 黄色女人牲交| 欧洲精品卡2卡3卡4卡5卡区| 一a级毛片在线观看| 亚洲欧美激情综合另类| 欧美精品亚洲一区二区| 亚洲成av片中文字幕在线观看| 一级a爱视频在线免费观看| www.自偷自拍.com| 免费观看人在逋| 午夜老司机福利片| 欧美精品啪啪一区二区三区| 男女做爰动态图高潮gif福利片 | 在线观看一区二区三区| 日韩欧美免费精品| 精品日产1卡2卡| 两个人看的免费小视频| 97超级碰碰碰精品色视频在线观看| av天堂久久9| 水蜜桃什么品种好| 少妇被粗大的猛进出69影院| 久久99一区二区三区| 正在播放国产对白刺激| 亚洲精品在线美女| 色综合站精品国产| 成人av一区二区三区在线看| 一区二区日韩欧美中文字幕| 18禁裸乳无遮挡免费网站照片 | 久久久久久人人人人人| 19禁男女啪啪无遮挡网站| 国产欧美日韩一区二区精品| 韩国av一区二区三区四区| 免费不卡黄色视频| 成人影院久久| 一二三四在线观看免费中文在| 国产免费av片在线观看野外av| 69精品国产乱码久久久| 久久精品国产99精品国产亚洲性色 | 亚洲自偷自拍图片 自拍| 国产乱人伦免费视频| 日韩 欧美 亚洲 中文字幕| 免费高清在线观看日韩| 在线av久久热| 日韩国内少妇激情av| avwww免费| 国产熟女午夜一区二区三区| 99国产精品免费福利视频| 国产黄色免费在线视频| 国产精品电影一区二区三区| 黄色 视频免费看| 琪琪午夜伦伦电影理论片6080| 一进一出抽搐动态| 欧美日韩国产mv在线观看视频| 大型黄色视频在线免费观看| 国产亚洲精品久久久久5区| 久久久国产一区二区| 国产亚洲精品综合一区在线观看 | 午夜91福利影院| 日韩精品青青久久久久久| 窝窝影院91人妻| 在线观看www视频免费| 亚洲成av片中文字幕在线观看| 一区二区三区精品91| 级片在线观看| 亚洲三区欧美一区| 亚洲精品美女久久av网站| 三上悠亚av全集在线观看| 男女下面进入的视频免费午夜 | 女警被强在线播放| 欧美日韩国产mv在线观看视频| 亚洲第一av免费看| 麻豆av在线久日| 日日摸夜夜添夜夜添小说| 免费看a级黄色片| 欧美日韩乱码在线| 色综合欧美亚洲国产小说| 欧美不卡视频在线免费观看 | 日韩 欧美 亚洲 中文字幕| 丝袜人妻中文字幕| xxx96com| 亚洲精品在线观看二区| 一边摸一边抽搐一进一小说| 最近最新中文字幕大全电影3 | 午夜福利欧美成人| 搡老乐熟女国产| 黑人巨大精品欧美一区二区mp4| 天堂动漫精品| av福利片在线| 亚洲欧美一区二区三区久久| 美女午夜性视频免费| 一边摸一边抽搐一进一小说| 人妻丰满熟妇av一区二区三区| 国产深夜福利视频在线观看| 日日干狠狠操夜夜爽| 久久草成人影院| 精品国产亚洲在线| 91成年电影在线观看| 性欧美人与动物交配| 成年版毛片免费区| 亚洲精品一二三| 18禁黄网站禁片午夜丰满| 不卡av一区二区三区| 他把我摸到了高潮在线观看| 亚洲人成伊人成综合网2020| 9热在线视频观看99| 19禁男女啪啪无遮挡网站| 美女扒开内裤让男人捅视频| 女人被躁到高潮嗷嗷叫费观| 国产伦人伦偷精品视频| 精品午夜福利视频在线观看一区| 在线观看免费视频日本深夜| 日本精品一区二区三区蜜桃| 一进一出抽搐动态| 精品国产国语对白av| 少妇被粗大的猛进出69影院| 99热国产这里只有精品6| 午夜福利影视在线免费观看| 亚洲男人天堂网一区| 精品一区二区三区四区五区乱码| av免费在线观看网站| 美女高潮到喷水免费观看| 欧美黄色淫秽网站| 热re99久久精品国产66热6| 欧美激情高清一区二区三区| 中文字幕最新亚洲高清| 女人被躁到高潮嗷嗷叫费观| 精品国产乱码久久久久久男人| 黄色丝袜av网址大全| 涩涩av久久男人的天堂| 久久亚洲真实| 91大片在线观看| 国产精品影院久久| 国产色视频综合| 欧美人与性动交α欧美软件| 国产精品一区二区三区四区久久 | 性色av乱码一区二区三区2| 极品人妻少妇av视频| 国产成+人综合+亚洲专区| 18禁国产床啪视频网站| 最新在线观看一区二区三区| 欧美 亚洲 国产 日韩一| 五月开心婷婷网| 国产欧美日韩综合在线一区二区| 首页视频小说图片口味搜索| 久久香蕉精品热| 国产高清激情床上av| 亚洲av第一区精品v没综合| 国产aⅴ精品一区二区三区波| 欧美性长视频在线观看| 国产精品久久久人人做人人爽| 国产成人av教育| 国产亚洲av高清不卡| 99精国产麻豆久久婷婷| 国产高清国产精品国产三级| 国产一区在线观看成人免费| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美精品综合一区二区三区| 叶爱在线成人免费视频播放| 精品久久久久久久久久免费视频 | 麻豆国产av国片精品| 亚洲专区中文字幕在线| www国产在线视频色| 男女下面进入的视频免费午夜 | 国产精华一区二区三区| 日韩免费av在线播放| 精品国产一区二区久久| 亚洲中文av在线| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美日韩无卡精品| 成人18禁在线播放| 日本黄色日本黄色录像| 亚洲国产欧美网| 91av网站免费观看| 黑人操中国人逼视频| 国产一区二区激情短视频| 久久久久久久久免费视频了| 欧美在线黄色| x7x7x7水蜜桃| 亚洲五月色婷婷综合| 亚洲性夜色夜夜综合| 丁香六月欧美| 午夜a级毛片| bbb黄色大片| 亚洲欧美日韩另类电影网站| 久久中文字幕人妻熟女| 久久国产亚洲av麻豆专区| 亚洲精品久久成人aⅴ小说| 亚洲成a人片在线一区二区| 一边摸一边做爽爽视频免费| 亚洲精品中文字幕一二三四区| 丰满迷人的少妇在线观看| 国产在线观看jvid| 美女 人体艺术 gogo| 每晚都被弄得嗷嗷叫到高潮| 一边摸一边抽搐一进一小说| 国产精品久久久久成人av| 亚洲伊人色综图| 欧美在线黄色| 日本撒尿小便嘘嘘汇集6| 欧美日韩精品网址| 一个人免费在线观看的高清视频| 国产在线精品亚洲第一网站| 深夜精品福利| 久久国产亚洲av麻豆专区| 欧美日韩乱码在线| 欧美精品亚洲一区二区| 国产精品偷伦视频观看了| 午夜精品国产一区二区电影| www国产在线视频色| 他把我摸到了高潮在线观看| 99国产精品99久久久久| 亚洲aⅴ乱码一区二区在线播放 | 麻豆国产av国片精品| 男女午夜视频在线观看| 精品日产1卡2卡| 日韩欧美三级三区| 久久 成人 亚洲| 亚洲av第一区精品v没综合| 久热爱精品视频在线9| av福利片在线| 精品一区二区三卡| 久久久久久久久免费视频了| 久久中文字幕人妻熟女| 欧美乱码精品一区二区三区| 国产精品日韩av在线免费观看 | 亚洲久久久国产精品| 亚洲国产中文字幕在线视频| 亚洲美女黄片视频| 国产麻豆69| 三上悠亚av全集在线观看| 别揉我奶头~嗯~啊~动态视频| 国产视频一区二区在线看| 亚洲国产欧美一区二区综合| 欧美激情久久久久久爽电影 | bbb黄色大片| 欧美+亚洲+日韩+国产| 热99re8久久精品国产| 国产欧美日韩一区二区三区在线| 99久久国产精品久久久| 女人被躁到高潮嗷嗷叫费观| 国产国语露脸激情在线看| 欧美日本亚洲视频在线播放| 两人在一起打扑克的视频| 天天影视国产精品| 午夜成年电影在线免费观看| 日韩欧美一区二区三区在线观看| 国产97色在线日韩免费| 一进一出抽搐gif免费好疼 | 妹子高潮喷水视频| 亚洲一卡2卡3卡4卡5卡精品中文| 久久香蕉精品热| 国产精品偷伦视频观看了| 精品国产亚洲在线| 久久精品亚洲av国产电影网| 国产无遮挡羞羞视频在线观看| 精品熟女少妇八av免费久了| 亚洲全国av大片| 99久久国产精品久久久| 三上悠亚av全集在线观看| 女警被强在线播放| 精品人妻1区二区| 人人妻人人澡人人看| 999久久久国产精品视频| 国产精品美女特级片免费视频播放器 | 成年人黄色毛片网站| 最好的美女福利视频网| 最新美女视频免费是黄的| 久久久久久大精品| 久久国产精品影院| 老司机午夜十八禁免费视频| 午夜免费观看网址| 国产亚洲精品一区二区www| 桃红色精品国产亚洲av| 久久国产乱子伦精品免费另类| av天堂久久9| 久久中文字幕一级| 亚洲精品在线美女| www国产在线视频色| 91精品三级在线观看| 91九色精品人成在线观看| 嫁个100分男人电影在线观看| 大型黄色视频在线免费观看| 亚洲欧美激情综合另类| 成人三级黄色视频| 国产精品永久免费网站| a级片在线免费高清观看视频| 日韩三级视频一区二区三区| 最新美女视频免费是黄的| 国产黄a三级三级三级人| 看片在线看免费视频| 黑丝袜美女国产一区| 美女福利国产在线| 久久久久久大精品| 在线观看舔阴道视频| 一个人免费在线观看的高清视频| 人人澡人人妻人| 成人亚洲精品一区在线观看| www日本在线高清视频| 电影成人av| 在线观看www视频免费| 男人操女人黄网站| 欧美日韩亚洲高清精品| 桃色一区二区三区在线观看| 亚洲第一av免费看| 国产精品 欧美亚洲| 国产精品久久久人人做人人爽| 夜夜躁狠狠躁天天躁| 欧美激情极品国产一区二区三区| 色综合婷婷激情| 中文字幕人妻丝袜制服| 国产亚洲av高清不卡| 如日韩欧美国产精品一区二区三区| 成人国产一区最新在线观看| 亚洲成国产人片在线观看| 两个人免费观看高清视频| 50天的宝宝边吃奶边哭怎么回事| 黑丝袜美女国产一区| 777久久人妻少妇嫩草av网站| 色精品久久人妻99蜜桃| 999精品在线视频| 久99久视频精品免费| 精品人妻1区二区| 成熟少妇高潮喷水视频| 国产蜜桃级精品一区二区三区| 亚洲在线自拍视频| 淫秽高清视频在线观看| 超色免费av| 变态另类成人亚洲欧美熟女 | 国产av精品麻豆| 国产精品久久久久成人av| 免费人成视频x8x8入口观看| 水蜜桃什么品种好| 在线观看免费日韩欧美大片| 好男人电影高清在线观看| www国产在线视频色| 美女午夜性视频免费| 亚洲国产中文字幕在线视频| 天堂√8在线中文| 十八禁网站免费在线| 国产精品电影一区二区三区| 日韩欧美免费精品| 亚洲欧美精品综合一区二区三区| 波多野结衣一区麻豆| 神马国产精品三级电影在线观看 | 一本综合久久免费| 久久精品91无色码中文字幕| 老熟妇仑乱视频hdxx|