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

    延伸沖擊-扇形氣膜孔復合結(jié)構耦合換熱實驗與數(shù)值研究

    2024-09-29 00:00:00吳航楊星蔡海揚劉釗豐鎮(zhèn)平
    西安交通大學學報 2024年10期

    摘要:針對沖擊-氣膜復合冷卻結(jié)構,采用實驗與數(shù)值方法,對比了基于圓形氣膜孔的傳統(tǒng)沖擊、延伸沖擊、小距離沖擊、基于扇形氣膜孔的傳統(tǒng)沖擊、延伸沖擊5種復合冷卻系統(tǒng)的耦合換熱和流動特性。采用紅外熱成像技術,獲得了5種復合冷卻結(jié)構在吹風比分別為0.6、1.0、1.5時外壁面的綜合冷卻有效度,并通過數(shù)值計算進一步揭示了內(nèi)部冷卻的流動和換熱特征。研究結(jié)果表明:在不產(chǎn)生額外氣動損失的前提下,延伸沖擊孔結(jié)構可提升內(nèi)部冷卻的換熱系數(shù),進而小幅度提升壁面的綜合冷卻有效度,幅度為1.2% ~ 4.6%,但隨著冷氣量的增大綜合冷卻有效度提升幅度有所減小;減小沖擊距離能夠提升內(nèi)部沖擊換熱效果,但不會對綜合冷卻有效度產(chǎn)生明顯影響;采用扇形氣膜孔可大幅度提升外部氣膜冷卻性能,且提升幅度大于采用延伸沖擊內(nèi)部改進結(jié)構的。相較于傳統(tǒng)的圓形氣膜孔-沖擊復合冷卻結(jié)構,在相同冷氣量條件下,基于扇形孔的延伸沖擊改進方案可將壁面的面平均綜合冷卻有效度提高7.6% ~ 8.5%,并將系統(tǒng)的流量系數(shù)提升30%以上。

    關鍵詞:耦合換熱;延伸沖擊;氣膜冷卻;綜合冷卻有效度

    中圖分類號:V231.1"文獻標志碼:A

    DOI:10.7652/xjtuxb202410014"文章編號:0253-987X(2024)10-0156-12

    Conjugate Heat Transfer of Extended Jet Impingement and Fan-Shaped

    Film Hole Composite Structure Experimental and Computational Approaches

    WU Hang1,2, YANG Xing1,2, CAI Haiyang1,2, LIU Zhao1,2, FENG Zhenping1,2

    (1. School of Energy and Power Engineering, Xi’an Jiaotong University, Xi’an 710049, China; 2. Shaanxi Engineering

    Laboratory of Turbomachinery and Power Equipment, Xi’an 710049, China)

    Abstract:For the impingement-film composite cooling structure, the conjugate heat transfer and flow characteristics of five composite cooling systems, including traditional jet impingement, extended jet impingement, and small-distance jet impingement, in combination with cylindrical film holes and fan-shaped film holes, are compared by using experimental and computational approaches. Infrared thermography is used in experiments to obtain overall cooling effectiveness of the outer wall surface under three blowing ratios: 0.6, 1.0 and 1.5. The computational approach is used to further reveal the flow and heat exchange characteristics of internal cooling details. Results show that the extended jet impingement hole structure can enhance the internal heat transfer coefficient without generating additional aerodynamic losses, which in turn improves overall cooling effectiveness by a small amount, ranging from 1.2% to 4.6%, but the improvement decreased with the increase of the coolant mass flow rate, that reducing the impingement distance can enhance the internal impingement heat transfer performance, which had no significant effect on overall cooling effectiveness and that the use of fan-shaped film holes can significantly improve the external film cooling performance, and the enhancement was greater than that through the internal optimization structure using extended jet holes. Compared with the traditional composite cooling structure of cylindrical film holes combined with impingement, the extended jet impingement scheme based on fan-shaped holes can improve the area-averaged overall cooling effectiveness by 7.6% to 8.5% and improve discharge coefficient by more than 30% under the same amount of coolant.

    Keywords:conjugate heat transfer; extended jet impingement; film cooling; overall cooling effectiveness

    沖擊-氣膜復合冷卻結(jié)構是一種常見的葉片冷卻強化技術,其結(jié)合內(nèi)、外多種冷卻方式,對冷氣的冷卻潛能利用更充分,因此廣泛應用于渦輪發(fā)動機燃燒室襯板和渦輪葉片前緣、中弦和端壁等區(qū)域的熱防護中。復合冷卻結(jié)構涉及沖擊冷卻、擾流柱、孔內(nèi)對流換熱和氣膜冷卻,各冷卻環(huán)節(jié)相互影響,產(chǎn)生了復雜的耦合換熱過程。掌握影響其冷卻性能的關鍵因素,有助于研發(fā)人員有針對性地優(yōu)化各冷卻環(huán)節(jié)。

    目前,考慮固體導熱的耦合換熱是研究復合冷卻結(jié)構綜合冷卻性能的主要方法。文獻[1-2]對比了有、無沖擊冷卻時復合層板結(jié)構的冷卻特性,發(fā)現(xiàn)沖擊-氣膜共同作用時冷卻效果更好,外表面溫度更低。Jung等[3]指出畢渥數(shù),即沖擊靶面的熱導率也是影響復合冷卻結(jié)構金屬表面溫度分布的重要因素。Tan等[4]通過調(diào)整沖擊孔與氣膜孔之間的排列方式與節(jié)距,形成了叉排和直列順排兩種方式,研究結(jié)果表明:隨著吹氣比增加,叉排和順排的冷效均升高,但叉排的冷效值始終比順排更高且更均勻,這種現(xiàn)象在孔節(jié)距較小時比較明顯。文獻[5-6]強調(diào),當沖擊孔出口正對氣膜孔進口時,冷氣會直接從氣膜孔中排出,從而減小了沖擊高換熱區(qū)的面積,降低了冷氣的利用率,在復合冷卻設計中應盡量避免。Rao等[7]在沖擊-氣膜復合冷卻中進一步增加擾流柱結(jié)構,通過增加內(nèi)部換熱面積略微提升了冷卻水平,但壓力損失大幅度增加,得不償失。文獻[8-9]拓展了沖擊-氣膜復合冷卻在葉片前緣的應用前景。Liu等[10]系統(tǒng)研究了調(diào)整氣膜孔直徑、復合角、吹風比等外部冷卻參數(shù)對雙層壁冷效的影響趨勢??梢钥闯觯鲜鰧秃侠鋮s結(jié)構的研究往往著眼于外部冷卻設計的調(diào)整,鮮有針對內(nèi)部冷卻設計的優(yōu)化。事實上,內(nèi)部沖擊冷卻對于復合結(jié)構冷卻性能的貢獻同樣重要,因此有必要同時關注內(nèi)部冷卻的強化設計。

    陣列沖擊射流的主要影響因素有沖擊雷諾數(shù)Re、沖擊孔間距xn/d和yn/d、沖擊距離zn/d以及橫流等[11]。為強化內(nèi)部沖擊冷卻性能,并降低上游橫流對下游射流的不利影響,學者們提出了多種新穎的強化沖擊結(jié)構,包括傾斜沖擊[12]、變節(jié)距陣列沖擊[13]、波紋沖擊板[14-15]、多級沖擊冷卻[16]與延伸沖擊孔[14,17-20]。其中,延伸沖擊孔使射流出口貼近沖擊靶面,能夠削弱沖擊射流動量的耗散,提升沖擊換熱水平,且不改變原有的沖擊腔高度,降低了橫流的不利影響,在內(nèi)部冷卻應用中效果顯著。Esposito等[14]通過實驗對比了傳統(tǒng)沖擊冷卻、波紋沖擊板與延伸沖擊的靶面換熱分布,結(jié)果表明,延伸沖擊孔可緩解橫流對下游沖擊的偏轉(zhuǎn)作用,相對于傳統(tǒng)沖擊冷卻與波紋沖擊板,換熱系數(shù)分別提升了40%、20%以上。Tepe等[17]從流動阻力和綜合換熱性能的角度全面評價了延伸沖擊孔,當沖擊腔高度為6倍孔直徑時,延伸孔的最佳長度為3~4倍孔徑。吳航等[18]則在更接近發(fā)動機幾何結(jié)構的小沖擊高度下探索了最佳的延伸孔長度。Yang等[19-20]設計了一種沿流向逐漸伸長的漸進延伸孔,在強橫流條件下開展了換熱實驗驗證,結(jié)果表明,與傳統(tǒng)沖擊冷卻相比,漸進延伸孔的局部換熱峰值明顯更高,盡管延伸孔需要更高的泵功來驅(qū)動,但強化換熱帶來的增益遠超流動損失。但是,上述公開文獻主要針對延伸沖擊孔自身的單獨作用開展研究,而沒有考慮其與外部氣膜抽吸和固體導熱共同存在時的耦合換熱特性,這與燃燒室和渦輪葉片等熱端部件的工程設計需求明顯不符。

    由于內(nèi)部沖擊冷卻、孔內(nèi)對流、外部氣膜抽吸和固體導熱等組成的耦合效應非常復雜,會顯著改變流場結(jié)構,進而影響延伸沖擊的局部換熱水平與分布,本文通過實驗與數(shù)值方法對比驗證了傳統(tǒng)沖擊、延伸沖擊、小距離沖擊、扇形氣膜孔-延伸沖擊和扇形氣膜孔-傳統(tǒng)沖擊5種復合冷卻系統(tǒng)的耦合換熱和流動特性,揭示了內(nèi)、外冷卻性能的耦合作用規(guī)律,最終得到了一種冷卻性能更佳、流動損失更小的延伸沖擊-扇形氣膜孔復合冷卻結(jié)構,以期為下一代航空發(fā)動機燃燒室和高壓渦輪葉片的高效冷卻設計提供借鑒。

    1"實驗與數(shù)值方法

    1.1"沖擊-氣膜復合冷卻實驗臺

    沖擊-氣膜復合冷卻實驗系統(tǒng)包括主流、冷氣與數(shù)據(jù)采集系統(tǒng),如圖1所示[21-22]。主流由4臺壓氣機提供,壓縮空氣依次經(jīng)過擴壓段和整流蜂窩網(wǎng)以確保來流均勻,并流經(jīng)湍流柵格后進入寬為180mm、高為120mm的實驗段。依據(jù)文獻[23]歸納的湍流柵格設計經(jīng)驗公式,選擇合適的湍流柵格直徑與節(jié)距,確保復合冷卻實驗件第一排氣膜孔處平均湍流強度為5.0%。主流進口的溫度、壓力、速度分別采用T型熱電偶、總壓探針和皮托管測量。主流通道上方開設紅外玻璃窗口用于紅外熱像儀拍攝。主流入口平均溫度為288K,平均速度為17.6m/s,根據(jù)入口速度和氣膜孔直徑,得到的進口雷諾數(shù)為3500,主流出口為1個大氣壓。

    實驗采用CO2作為二次流,依次經(jīng)過質(zhì)量流量計與換熱器后溫度達到323K。二次流與主流的密度比rDR為1.5,接近發(fā)動機實際條件。在供氣腔側(cè)壁布置兩個引壓管與T型熱電偶,用于測量冷氣入口的壓力和溫度。二次流管道和腔室均被保溫棉包裹,以減少熱損失并提供更均勻的溫度。

    本研究復合冷卻結(jié)構圖2(a)所示,氣膜孔與沖擊孔均叉排布置,流向上布置7排孔,且沖擊孔與氣膜孔之間交錯排列。氣膜孔直徑Df=3mm,射流角θ=30°,沖擊孔直徑Di =3.5mm,流向間距Px=10Df,展向間距Py=5Df,沖擊板和氣膜板的厚度均為10mm。為避免邊緣效應的影響,以中間2排孔所在位置作為測量區(qū)域,并作為后續(xù)線平均與面平均數(shù)據(jù)的統(tǒng)計區(qū)域。將第一排氣膜孔出口圓心為無量綱軸向位置x/Df=0,則沖擊孔實際布置區(qū)域大致為x/Df=-5~25。

    本文重點關注耦合換熱條件下延伸沖擊孔的強化冷卻性能,考慮了5組冷卻結(jié)構,分別為基于圓形氣膜孔的延伸沖擊(CEX)、傳統(tǒng)沖擊(C3D)、小距離沖擊(C1D),以及基于扇形氣膜孔的傳統(tǒng)沖擊(F3D)、延伸沖擊(FEX),分別如圖2(b)所示。在延伸沖擊孔結(jié)構中,沖擊高度H=3Di,延伸孔長度L=2Di,將傳統(tǒng)沖擊的沖擊高度H=3Di與小距離沖擊H=Di作為內(nèi)部冷卻對照組。與C3D相比,F(xiàn)3D僅將外部氣膜孔型變?yōu)樯刃?,具體采用7-7-7型孔[24]。金屬外壁面為本文關注的測量面,內(nèi)表面為沖擊靶面。

    為使實驗獲得的綜合冷卻有效度更好地代表航空發(fā)動機渦輪端壁的實際耦合換熱過程,將實驗模型的外、內(nèi)換熱系數(shù)之比h∞/hi、畢渥數(shù)Bi與實際情況進行了匹配,如表1所示。在幾何參數(shù)、流動工況合理縮放和匹配后,h∞/hi將自動匹配,氣膜板選用316L不銹鋼材料。為了盡可能地降低熱損失,冷氣供氣腔和沖擊板采用導熱系數(shù)為0.2W·m-1·K-1的樹脂材料3D打印而成。

    1.2"測量方法與不確定度分析

    采用紅外熱像儀FLIR T865測量金屬的外壁面溫度,儀器的分辨率為640×480像素,實際測量時的物理分辨率為4像素/mm,工作波長為7.5~14μm。選用鍍膜鍺片作為紅外窗口,可保證在工作波長范圍內(nèi)有90%以上的透射率。采用紅外熱像儀測溫前,使用T型熱電偶先對其進行標定。實驗前,使用啞光黑漆均勻地噴涂于金屬外壁面,以保證全表面發(fā)射率均勻。標定時,利用加熱膜使表面溫度在實驗測量范圍內(nèi)變化,并在壁面相同位置處同時讀取熱電偶測點溫度T和紅外熱像儀的溫度TIR,得到多組一一對應的溫度數(shù)據(jù),處理得到如圖3所示的標定曲線。圖中給出了擬合得到標定曲線的具體表達式以及決定系數(shù)R2=0.9999,表明擬合曲線可精確完成紅外溫度數(shù)據(jù)的標定。

    綜合冷卻有效度(以下簡稱綜合冷效)是一個無量綱溫度參數(shù),用來評估耦合換熱條件下復合冷卻結(jié)構的冷卻效果,定義如下

    Φ=T∞-TwT∞-Tc, in(1)

    式中:Tw為金屬壁面溫度;T∞、Tc,in分別為T型熱電偶測量的主流溫度和冷氣溫度。

    本實驗采用I級精度的T型熱電偶,在-40~125℃測量范圍內(nèi),測量誤差為0.5℃,由標定公式可知,紅外熱像儀測得壁溫的不確定度絕對值不大于0.5℃。根據(jù)文獻[25]的實驗誤差分析方法,綜合冷卻有效度為 0.15、0.75兩個典型值時的相對不確定度分別為14.0%、2.6%。

    1.3"參數(shù)定義

    外部氣膜冷卻吹風比定義為

    M=c/Ac∞/A∞(2)

    式中:c、∞分別為冷氣和主流的質(zhì)量流量;Ac、A∞分別為氣膜孔和主流通道的通流面積。

    內(nèi)部沖擊雷諾數(shù)是與外部氣膜冷卻吹風比對應的參數(shù),定義為

    Rej=4cπDiμN(3)

    式中:μ為射流的動力黏性系數(shù);N為沖擊孔個數(shù)。本文選取了3組吹風比M=0.6,1.0,1.5,對應的沖擊雷諾數(shù)Rej=2050,3400,5100。

    采用靶面努塞爾數(shù)衡量內(nèi)部沖擊冷卻的換熱系數(shù),可寫為

    Nu=qwDi(Tw,in-Tc)λ(4)

    式中:qw為壁面熱流密度;Tc為射流入口總溫;Tw,in為靶面溫度;λ為冷卻工質(zhì)的導熱系數(shù)。

    采用無量綱溫度θ來衡量主流與冷氣的局部混合程度,該數(shù)值越接近1,表明流體域局部冷氣越多,固體域冷卻效果越好,定義為

    θ=T∞-TlT∞-Tc, in(5)

    式中:Tl為局部溫度;Tc,in為冷氣進口溫度。

    采用流量系數(shù)Cd衡量整個沖擊-氣膜復合冷卻系統(tǒng)的流動損失,定義為

    Cd=mi(6)

    i=AcPt,c(Ps/Pt,c)(κ+1)/2κ·

    2κ(κ-1)RTc,in[(Pt,c/Ps)(κ-1)/κ-1] (7)

    式中:m、i分別為實際測量和理想的冷氣質(zhì)量流量;Ac為冷氣出流面積,即氣膜孔出口面積之和;R、κ分別為冷氣的氣體常數(shù)和絕熱系數(shù);Pt,c、Ps分別為沖擊入口總壓與氣膜孔出口靜壓,兩者之差即為復合冷卻系統(tǒng)的壓降。

    1.4"數(shù)值計算模型與網(wǎng)格劃分

    本文結(jié)合數(shù)值計算來補充實驗中難以獲得的內(nèi)部強化換熱特征與流動機理。沖擊-氣膜復合冷卻數(shù)值計算模型的計算域與邊界條件如圖4所示,模型只繪制中間兩排孔,兩側(cè)為周期性邊界條件,沖擊孔、氣膜孔與主流流道等幾何尺寸均與實驗一致。主流工質(zhì)為空氣,入口溫度和速度邊界條件與實驗相同,主流出口為大氣壓力。冷氣工質(zhì)設置為CO2,在冷氣入口處給定與實驗相同的質(zhì)量流量和總溫條件,壁面材料選用316L不銹鋼,導熱系數(shù)設置為溫度的函數(shù)。固體-流體接觸面,如外壁面、內(nèi)部沖擊靶面與氣膜孔面等設置為熱耦合交界面,其余面均設置為絕熱面。利用商用軟件 ANSYS CFX求解穩(wěn)態(tài)條件下的三維可壓縮雷諾時均方程。

    數(shù)值模型采用ANSYS Mesh軟件生成非結(jié)構化網(wǎng)格,主流區(qū)域為四面體,在氣膜孔、沖擊孔與壁面局部加密,壁面邊界層用棱柱體填充,第一層網(wǎng)格高度為0.005mm,保證y+lt;1.5。針對延伸沖擊結(jié)構,分別采用630萬、1180萬、1950萬依次細化的網(wǎng)格數(shù)開展網(wǎng)格無關性驗證。驗證工況的冷氣吹風比為1.0,將綜合冷效沿著橫向(y軸)進行平均并沿著流向(x軸)繪制,得到3套網(wǎng)格的橫向平均綜合冷卻有效度,如圖5所示,可知當網(wǎng)格數(shù)達到1180萬時,繼續(xù)提升網(wǎng)格數(shù),壁面綜合冷效基本不再發(fā)生變化。因此,選用網(wǎng)格數(shù)為1180萬開展后續(xù)計算,其中流體域1080萬,固體域100萬,其他模型的網(wǎng)格參數(shù)設置均與此相同。

    1.5"湍流模型驗證

    前人的研究結(jié)果表明,SST k-ω湍流模型對沖擊-氣膜復合冷卻結(jié)構的耦合換熱特性具有良好的預測精度[1,5,7,10]。為此,本文采用SST k-ω湍流模型,以延伸沖擊M=1.5工況下的實驗結(jié)果來驗證數(shù)值計算方法的準確性。圖6對比了實驗測量和數(shù)值模擬的壁面綜合冷效分布,數(shù)值選取的湍流模型對壁溫的計算結(jié)果與實驗具有高度相似性,正確捕捉到了固體橫向?qū)岬姆秶c強度,可以準確預測壁面的綜合冷效水平和流場細節(jié)。

    圖6(b)進一步定量比較了中間一排氣膜孔中心線上的綜合冷效,可知數(shù)值模擬與實驗測量之間的綜合冷效局部值非常接近,且沿流向的變化趨勢相似。因此,本文所有后續(xù)計算均采用 SST k-ω湍流模型。

    2"結(jié)果與討論

    本節(jié)將結(jié)合實驗測量與數(shù)值模擬結(jié)果,充分對比基于外部扇形氣膜孔改進與內(nèi)部沖擊結(jié)構改進下的壁面綜合冷效、內(nèi)部強化換熱分布、壓力損失系數(shù)、氣膜孔出口流場與溫度場,以全面衡量扇形氣膜孔與延伸沖擊組合應用的優(yōu)勢。

    2.1"耦合換熱特性

    圖7展示了實驗測量的5種冷卻結(jié)構在3個吹風比工況下的壁面綜合冷效分布云圖。相應地,為了量化冷氣質(zhì)量流量的影響,綜合冷效的橫向平均值沿流向x軸的變化曲線如圖8所示,圖中豎直方向的7條灰色條帶對應7排氣膜孔的出口位置。由圖7、8可知,在不同內(nèi)部沖擊冷卻結(jié)構下,隨著吹風比的增加,整個壁面的綜合冷卻效果都有所提高,且延伸沖擊和高低沖擊距下沿流向的橫向平均綜合冷效的變化趨勢均很相似,即綜合冷效先增加,然后再向下游逐漸降低,該值的增加主要得益于內(nèi)部沖擊換熱、孔內(nèi)對流換熱的提升及固體的導熱作用。在絕熱條件下,圓形氣膜冷效的最佳吹風比約在0.6附近,扇形氣膜冷效的最佳吹風比約在1.5附近,繼續(xù)提高冷氣量會降低冷氣射流的附壁性,反而惡化氣膜冷卻效果[26]。在本研究的耦合換熱條件下,由于內(nèi)部沖擊冷卻與孔內(nèi)換熱的存在,綜合冷效會隨吹風比的提高而持續(xù)提升,這也表明內(nèi)部沖擊冷卻特性是影響沖擊/氣膜復合結(jié)構冷卻效果的主導因素。另一方面,在固體導熱作用下,外壁面也不會呈現(xiàn)出和絕熱氣膜冷效一樣明顯的高冷效射流軌跡,綜合冷效更加均勻。以上兩點也是氣膜冷卻流熱耦合研究與絕熱研究的主要區(qū)別,總體來看,在3個吹風比工況下,結(jié)合扇形孔與延伸沖擊兩種內(nèi)外改進的FEX模型均擁有最高的綜合冷效值。

    由圖8可以看出:所有工況下綜合冷效高值區(qū)域更多地集中在壁面中心區(qū)域(x/Df=15~30,y/Df= -2~2)且向四周遞減,邊緣處最低,這主要是因為中心區(qū)域是沖擊孔作用的核心位置,冷卻最好;壁面四周位于沖擊邊緣,內(nèi)部換熱系數(shù)降低,且固體會向更遠處未被冷卻措施覆蓋的區(qū)域?qū)?,進而中心溫度高,四周溫度低,這是耦合換熱的機理所致。

    在相同吹風比條件下,圖9給出了5種冷卻結(jié)構的橫向平均綜合冷效。通過對比C3D與C1D的實驗結(jié)果可知,在H=1~3Di范圍內(nèi),增大或減小沖擊距離對綜合冷效的影響有限。而將CEX與上述兩者進行對比可以發(fā)現(xiàn),延伸沖擊綜合冷效的局部提升僅在壁面中心區(qū)域幅度較大,最高可達5%,而在x/Dflt;5和x/Dfgt;30的上下游沖擊邊緣位置,改進沖擊結(jié)構的提升并不明顯,這說明僅在沖擊孔作用的核心區(qū)域x/Df=15~30對內(nèi)部冷卻進行改進,才能對冷卻性能的提升發(fā)揮作用。

    進一步對比兩種扇形氣膜孔與其他3種圓形氣膜孔結(jié)構可以發(fā)現(xiàn),扇形氣膜孔-傳統(tǒng)沖擊結(jié)構與圓形氣膜孔-延伸沖擊在x/Df=15~20的位置有一個冷卻性能的分水嶺,在x/Dflt;15~20的上游處CEX結(jié)構的綜合冷效更高,而x/Dfgt;15~20下游處F3D結(jié)構的綜合冷效更高。這主要由于上游前兩排孔處冷卻氣膜并沒有充分形成,扇形氣膜孔的外部冷卻優(yōu)勢沒有完全展現(xiàn),因而上游處擁有更強內(nèi)部冷卻性能的CEX結(jié)構的綜合冷效會更高。在下游第5、6排氣膜孔的位置,冷卻氣膜充分發(fā)展,扇形氣膜孔的冷卻優(yōu)勢得以展現(xiàn),并且x/Dfgt;15~20的下游已位于內(nèi)部沖擊冷卻的邊緣區(qū),因而下游處F3D結(jié)構的綜合冷效更高。這也表明沖擊-氣膜復合冷卻結(jié)構中,有效的沖擊與氣膜冷卻改進設計均可以發(fā)揮重要作用。得益于附壁性更好的外部氣膜冷卻特性與內(nèi)部強化換熱特性,扇形孔-延伸沖擊結(jié)構在流向10lt;x/Dflt;30范圍內(nèi)的冷卻性能明顯高于其他4種結(jié)構,且隨著吹風比的增大,冷效提升愈加明顯;相反地,圓形氣膜孔-延伸沖擊結(jié)構與傳統(tǒng)沖擊結(jié)構之間的冷效差值隨著吹風比的增大而減小,在高吹風比時幾乎可忽略不計。

    圖10給出了圖7中云圖展示區(qū)域的面積平均綜合冷效對比。由圖10可知:FEX結(jié)構的面平均綜合冷效比F3D的高約2%,比CEX的高3.5%~7.2%,比傳統(tǒng)的C3D結(jié)構高7.6%~8.5%;CEX結(jié)構較傳統(tǒng)沖擊結(jié)構C3D有小幅度提升,提升幅度約為1.2%~4.6%。

    2.2"內(nèi)部強化換熱特性

    為了更好地理解延伸沖擊和沖擊間距對綜合冷效的影響機理,本節(jié)選取圓形氣膜孔對應的3種冷卻結(jié)構,討論數(shù)值計算得到的內(nèi)部換熱分布和流場、溫度場。

    3種內(nèi)部沖擊冷卻結(jié)構的靶面努塞爾數(shù)分布云圖如圖11所示,可以看到叉排陣列射流的典型Nu數(shù)的分布特征為:沖擊駐點區(qū)湍流度高且邊界層被破壞,發(fā)生強烈的換熱,此時存在換熱峰值,在駐點區(qū)外的射流形成壁面射流,因而速度逐漸減低,換熱系數(shù)逐漸降低。圖11(b)、11(c)兩種沖擊距離的對比表明,隨著沖擊距離增大,滯止區(qū)附近的換熱顯著減小,射流從沖擊孔流出后,由于剪切流動與周圍流體產(chǎn)生動量交換,隨著射流距離的增加,射流核心區(qū)速度逐漸下降,導致?lián)Q熱降低,這種現(xiàn)象在高射流雷諾數(shù)下的對比更為明顯。圖11(a)所示的延伸沖擊在沖擊駐點區(qū)也形成明顯的換熱峰值,換熱分布情況與小距離沖擊的結(jié)果更加接近,這體現(xiàn)了延伸沖擊的關鍵優(yōu)勢,即在不改變沖擊腔高度的前提下,提升射流抵達靶面時的峰值速度。此外,射流孔之間的邊緣區(qū)域存在換熱低谷值,且x/Dfgt;25后的區(qū)域由于沒有布置沖擊孔,下個節(jié)距不再存在強烈的沖擊換熱,導致外部綜合冷效在x/Dfgt;25后快速降低。

    圖12展示了不同射流雷諾數(shù)下3種冷卻結(jié)構靶面的橫向平均努塞爾數(shù)分布,無量綱軸向坐標與研究外部綜合冷效時選取一致,圖中豎直方向的7條藍色條帶對應7排沖擊孔的位置。射流雷諾數(shù)作為最基本的氣動參數(shù),主要通過影響射流速度來影響換熱特性。由圖12可見,3種冷卻結(jié)構下,靶面努塞爾數(shù)均隨雷諾數(shù)的增大而增大,射流駐點區(qū)域為換熱峰值,而谷值位于兩排沖擊孔之間,射流駐點區(qū)的換熱系數(shù)隨雷諾數(shù)增加尤為明顯。

    3種冷卻結(jié)構在不同雷諾數(shù)下靶面的橫向平均努塞爾數(shù)分布如圖13所示,可知大沖擊距離時靶面換熱系數(shù)最低,而延伸沖擊與小距離沖擊主要提升了靶面滯止區(qū)的換熱峰值。Rej=2050時,相較于大沖擊距離的傳統(tǒng)沖擊,延伸與小距離沖擊的靶面面積平均努塞爾數(shù)分別提升了12.5%、11.1%;Rej=5100時,提升幅度縮減為9.6%、7.2%。總體而言,隨著雷諾數(shù)的增大,采用延伸沖擊或減小距離沖擊,均降低了內(nèi)部冷卻的提升幅度。

    圖14、圖15分別展示了Rej=3400時3種冷卻結(jié)構孔中心截面上的速度與無量綱溫度分布云圖,此截面位于中間排氣膜孔與沖擊孔之間的中心位置。對比圖14(b)、(c)可知:在大沖擊距離下,射流抵達壁面后速度有所降低,最終導致如圖11所反映出的冷卻性能的衰退;而在延伸沖擊與小距離沖擊下,抵達靶面時射流核心區(qū)速度較高,換熱較強。由圖15可知,在氣膜射流與主流摻混后,在近外壁面區(qū)域形成一層冷卻氣膜,并且沿著流向氣膜逐漸增厚且溫度逐漸降低,但壁面外表面溫度卻沒有隨之下降,這是因為下游位置已經(jīng)超出內(nèi)部沖擊的作用范圍,這也反映出沖擊冷卻是影響壁面綜合冷效的重要因素。此外,由于固體的橫向?qū)嶙饔?,壁面外?cè)的溫度梯度略小于內(nèi)側(cè)靶面。

    2.3"流動損失特性

    本節(jié)采用流量系數(shù)Cd衡量整個沖擊-氣膜復合冷卻系統(tǒng)的流動損失。圖16給出了不同冷卻結(jié)構下系統(tǒng)流量系數(shù)的實驗結(jié)果。沖擊冷卻孔的典型Cd為0.84[19],傾斜角為30°的圓形氣膜孔的典型Cd為0.74[26]。由于本研究考慮了沖擊射流和氣膜射流共同作用下的壓力損失,因此流量系數(shù)更小。

    由圖16可知,延伸沖擊與減小距離沖擊等內(nèi)部改進措施不會產(chǎn)生明顯的額外壓力損失,而氣膜孔型的改變對Cd有著顯著影響。一方面,扇形孔有更大的通流面積,大大降低了氣膜孔內(nèi)的流動損失;另一方面,沖擊系統(tǒng)的壓力損失主要發(fā)生在沖擊孔入口收縮流動和出口擴張流動過程,而沖擊距離的影響相對較小。因此,CEX、C3D、C1D結(jié)構的流量系數(shù)幾乎一致,而FEX、F3D結(jié)構的流動損失較前三者明顯更低,降低幅度達30%以上。

    3"結(jié)"論

    采用實驗測量和數(shù)值模擬方法,研究了改進內(nèi)部沖擊及外部氣膜孔型對沖擊-氣膜復合冷卻結(jié)構流動與換熱特性的影響,詳細分析了外壁面綜合冷效、靶面內(nèi)部換熱系數(shù)、流場和溫度場以及流動損失的變化規(guī)律,并與傳統(tǒng)復合冷卻結(jié)構進行對比,表明結(jié)合內(nèi)外兩種改進的延伸沖擊-扇形氣膜孔復合冷卻結(jié)構具有最佳的冷卻性能。得到的主要結(jié)論如下。

    (1) 與采用傳統(tǒng)內(nèi)部沖擊冷卻的復合冷卻結(jié)構相比,延伸沖擊可以提升內(nèi)部冷卻的換熱系數(shù),進而小幅度提升壁面的綜合冷效。然而,隨著冷氣量的增大,冷效的提高幅度變小,減小沖擊距離會提升內(nèi)部沖擊換熱效果,但不會對綜合冷效產(chǎn)生明顯的影響。

    (2) 采用扇形氣膜孔可以有效提高復合冷卻系統(tǒng)的冷卻性能。結(jié)合扇形孔和延伸沖擊兩種內(nèi)、外改進的FEX結(jié)構方案具有最高的綜合冷效效果,相較于傳統(tǒng)的圓形氣膜孔-沖擊復合冷卻結(jié)構,效果提高了7.6%~8.5%,較圓形孔-延伸沖擊結(jié)構的提升幅度也可達3.5%~7.2%。

    (3) 在沖擊距離H/Di=1~3時,減小沖擊距離或者采用延伸沖擊等內(nèi)部冷卻改進措施不會對復合冷卻系統(tǒng)產(chǎn)生明顯的額外壓力損失。采用扇形氣膜孔外部改進方案的氣動性能有明顯的提升,其中流量系數(shù)的提升幅度達30%以上。

    參考文獻:

    [1]PANDA R K, PRASAD B. Conjugate heat transfer from a flat plate with combined impingement and film cooling [C]//ASME Turbo Expo 2012: Turbine Technical Conference and Exposition. New York, USA: ASME, 2012: 347-356.

    [2]費昕陽, 王新軍, 陸??? 平板沖擊發(fā)散冷卻流動與換熱特性的數(shù)值模擬 [J]. 西安交通大學學報, 2017, 51(7): 57-61.

    FEI Xinyang, WANG Xinjun, LU Haikong. Numerical simulation for flow and heat transfer characteristics of impingement/effusion cooling on a flat plane [J]. Journal of Xi’an Jiaotong University, 2017, 51(7): 57-61.

    [3]JUNG E Y, CHUNG H, CHOI S M, et al. Conjugate heat transfer on full-coverage film cooling with array jet impingements with various Biot numbers [J]. Experimental Thermal and Fluid Science, 2017, 83: 1-8.

    [4]TAN Xiaoming, ZHANG Jingzhou, XU Huasheng. Experimental investigation on impingement/effusion cooling with short normal injection holes [J]. International Communications in Heat and Mass Transfer, 2015, 69: 1-10.

    [5]XIE Gang, LIU Cunliang, YE Lin, et al. Effects of impingement gap and hole arrangement on overall cooling effectiveness for impingement/effusion cooling [J]. International Journal of Heat and Mass Transfer, 2020, 152: 119449.

    [6]CHO H H, RHEE D H, GOLDSTEIN R J. Effects of hole arrangements on local heat/mass transfer for impingement/effusion cooling with small hole spacing [J]. Journal of Turbomachinery, 2008, 130(4): 041003.

    [7]RAO Yu, LIU Yuyang, WAN Chaoyi. Multiple-jet impingement heat transfer in double-wall cooling structures with pin fins and effusion holes [J]. International Journal of Thermal Sciences, 2018, 133: 106-119.

    [8]劉釗, 賈哲, 張志欣, 等. 透平動葉前緣沖擊-氣膜復合冷卻與旋流-氣膜復合冷卻的熱流耦合對比研究 [J]. 西安交通大學學報, 2021, 55(4): 116-125.

    LIU Zhao, JIA Zhe, ZHANG Zhixin, et al. A comparative study on conjugate heat transfer of impingement-film composite cooling and swirl-film composite cooling on leading edge of a turbine blade [J]. Journal of Xi’an Jiaotong University, 2021, 55(4): 116-125.

    [9]鄧清華, 王輝輝, 何偉, 等. 氣膜抽吸與葉片前緣內(nèi)部壁面射流的作用機制 [J]. 西安交通大學學報, 2022, 56(10): 160-169.

    DENG Qinghua, WANG Huihui, HE Wei, et al. Mechanism of action between film suction and wall jet in leading edge of blades [J]. Journal of Xi’an Jiaotong University, 2022, 56(10): 160-169.

    [10]LIU Xuebin, ZHANG Chao, SONG Liming, et al. Influence of Biot number and geometric parameters on the overall cooling effectiveness of double wall structure with pins [J]. Applied Thermal Engineering, 2021, 198: 117439.

    [11]WEIGAND B, SPRING S. Multiple jet impingement: a review [J]. Heat Transfer Research, 2011, 42(2): 101-142.

    [12]LI Weihong, XU Minghe, REN Jing, et al. Experimental investigation of local and average heat transfer coefficients under an inline impinging jet array, including jets with low impingement distance and inclined angle [J]. Journal of Heat Transfer,2017, 139(1): 012201.

    [13]GAO Lujia, EKKAD S V, BUNKER R S. Impingement heat transfer part I: linearly stretched arrays of holes [J]. Journal of Thermophysics and Heat Transfer, 2005, 19(1): 57-65.

    [14]ESPOSITO E I, EKKAD S V, KIM Y, et al. Novel jet impingement cooling geometry for combustor liner backside cooling [J]. Journal of Thermal Science and Engineering Applications, 2009, 1(2): 021001.

    [15]KIM T, JUNG E Y, BANG M, et al. Heat transfer measurements for array jet impingement with castellated wall [J]. Journal of Turbomachinery,2022, 144(3): 031009.

    [16]LIU Kexin, ZHANG Qiang. A novel multi-stage impingement cooling scheme: part Ⅰ"concept study [J]. Journal of Turbomachinery, 2020, 142(12): 121008.

    [17]TEPE A , YETIKEN Y, UYSAL , et al. Experimental and numerical investigation of jet impingement cooling using extended jet holes [J]. International Journal of Heat and Mass Transfer, 2020, 158: 119945.

    [18]吳航, 楊星, 趙強, 等. 延伸沖擊孔沖擊冷卻流動與換熱特性的數(shù)值研究 [J]. 推進技術, 2022, 43(10): 291-298.

    WU Hang, YANG Xing, ZHAO Qiang, et al. Numerical study on flow and heat transfer characteristics of impingement cooling with extended jet holes [J]. Journal of Propulsion Technology, 2022, 43(10): 291-298.

    [19]YANG Xing, WU Hang, FENG Zhenping. Jet impingement heat transfer characteristics with variable extended jet holes under strong crossflow conditions [J]. Aerospace, 2022, 9(1): 44.

    [20]YANG Xing, WU Hang, FENG Zhenping. Experimental and numerical investigations of extended jet effects on multiple-jet impingement heat transfer characteristics [J]. Journal of Thermal Science and Engineering Applications, 2022, 14(10): 101015.

    [21]YANG Xing, WU Hang, HAO Zihan, et al. Assessment of a virtual boundary concept for efficient modeling of turbine film cooling [J/OL]. Numerical Heat Transfer: Part A"Applications. [2023-06-23]. https://doi.org/10.1080/10407782.2023.2299286.

    [22]WU Hang, YANG Xing, WU Yongqiang, et al. Experimental decoupled-analysis of overall cooling effectiveness for a turbine endwall with internal and external cooling configurations [J]. Applied Thermal Engineering, 2023, 228: 120435.

    [23]ROACH P E. The generation of nearly isotropic turbulence by means of grids [J]. International Journal of Heat and Fluid Flow, 1987, 8(2): 82-92.

    [24]THOLE K, GRITSCH M, SCHULZ A, et al. Flowfield measurements for film-cooling holes with expanded exits [C]//ASME 1996 International Gas Turbine and Aeroengine Congress and Exhibition. New York, USA: ASME, 1996: V004T09A010.

    [25]MOFFAT R J. Describing the uncertainties in experimental results [J]. Experimental Thermal and Fluid Science, 1988, 1(1): 3-17.

    [26]SAUMWEBER C, SCHULZ A. Effect of geometry variations on the cooling performance of fan-shaped cooling holes [J]. Journal of Turbomachinery, 2012, 134(6): 061008.

    (編輯"趙煒"李慧敏)

    大香蕉久久网| 毛片女人毛片| 国产成人a∨麻豆精品| 久久久色成人| 99国产精品一区二区蜜桃av| 韩国高清视频一区二区三区| 欧美日韩国产亚洲二区| 国产免费福利视频在线观看| 我的女老师完整版在线观看| 亚洲精品乱码久久久v下载方式| 女人久久www免费人成看片 | 久久精品91蜜桃| 村上凉子中文字幕在线| 国产免费福利视频在线观看| 久久精品国产亚洲网站| 国产在视频线精品| 午夜爱爱视频在线播放| 久久精品夜色国产| 三级毛片av免费| 99热这里只有精品一区| 国产精品伦人一区二区| 亚洲国产精品国产精品| 观看免费一级毛片| 视频中文字幕在线观看| 成人毛片a级毛片在线播放| 欧美高清成人免费视频www| 久久久a久久爽久久v久久| av专区在线播放| 级片在线观看| 国产高清三级在线| 国产一级毛片在线| 中文字幕亚洲精品专区| 如何舔出高潮| 最近的中文字幕免费完整| 日韩三级伦理在线观看| 国产精品日韩av在线免费观看| av免费在线看不卡| 亚洲精品日韩在线中文字幕| 国产伦在线观看视频一区| 国产综合懂色| 色尼玛亚洲综合影院| 欧美极品一区二区三区四区| 国产91av在线免费观看| 日本黄色片子视频| 国产v大片淫在线免费观看| 成年av动漫网址| 国产成人精品一,二区| 中国美白少妇内射xxxbb| 国产精品一区二区三区四区免费观看| kizo精华| 亚洲综合精品二区| 联通29元200g的流量卡| 天天躁夜夜躁狠狠久久av| 亚洲av成人av| 狂野欧美白嫩少妇大欣赏| 成人综合一区亚洲| 男女啪啪激烈高潮av片| 青春草亚洲视频在线观看| 亚洲欧美中文字幕日韩二区| 久久亚洲国产成人精品v| 国产精品综合久久久久久久免费| 久久亚洲精品不卡| 免费看美女性在线毛片视频| 国产精品国产三级国产专区5o | 国产熟女欧美一区二区| 亚洲怡红院男人天堂| 综合色丁香网| 久热久热在线精品观看| 国产黄片美女视频| 色吧在线观看| 亚洲国产精品成人久久小说| 久久这里有精品视频免费| 永久免费av网站大全| 精品午夜福利在线看| 最近中文字幕高清免费大全6| 国产69精品久久久久777片| 国产免费又黄又爽又色| 91午夜精品亚洲一区二区三区| 国产成人福利小说| 亚洲一区高清亚洲精品| 色综合色国产| 99久久中文字幕三级久久日本| 我要搜黄色片| 久久精品国产亚洲av天美| 久久久久久久久久久免费av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精品自拍成人| 日本免费a在线| 啦啦啦啦在线视频资源| 国产毛片a区久久久久| 亚洲精品日韩av片在线观看| 99久久无色码亚洲精品果冻| 久久久a久久爽久久v久久| 日本-黄色视频高清免费观看| 99久久精品热视频| 欧美色视频一区免费| 免费观看精品视频网站| av在线观看视频网站免费| 青春草亚洲视频在线观看| 国产单亲对白刺激| 日本熟妇午夜| 婷婷六月久久综合丁香| 在现免费观看毛片| 真实男女啪啪啪动态图| 中文字幕av成人在线电影| 综合色av麻豆| 一卡2卡三卡四卡精品乱码亚洲| 尤物成人国产欧美一区二区三区| 欧美三级亚洲精品| 淫秽高清视频在线观看| 2021少妇久久久久久久久久久| av线在线观看网站| 边亲边吃奶的免费视频| 欧美一区二区亚洲| 秋霞在线观看毛片| 欧美不卡视频在线免费观看| 国产一区亚洲一区在线观看| 久久精品夜色国产| 亚洲色图av天堂| 亚洲怡红院男人天堂| 亚洲精品国产av成人精品| 99热这里只有是精品50| 日韩av在线免费看完整版不卡| 男女边吃奶边做爰视频| 如何舔出高潮| 精品一区二区三区人妻视频| 一级毛片久久久久久久久女| 国产又黄又爽又无遮挡在线| 亚洲精品日韩在线中文字幕| 亚洲五月天丁香| 国产大屁股一区二区在线视频| 国产精品国产高清国产av| 九九久久精品国产亚洲av麻豆| 日韩欧美精品免费久久| 赤兔流量卡办理| 国产三级在线视频| 欧美xxxx黑人xx丫x性爽| 欧美高清性xxxxhd video| 久久久欧美国产精品| 91久久精品电影网| 九色成人免费人妻av| 日本av手机在线免费观看| 成人三级黄色视频| 国产高清不卡午夜福利| 国产极品天堂在线| 精品国产三级普通话版| 一个人观看的视频www高清免费观看| 欧美3d第一页| 久久久久久久久中文| 人体艺术视频欧美日本| 亚洲伊人久久精品综合 | 免费观看的影片在线观看| 18禁动态无遮挡网站| 国产欧美日韩精品一区二区| 亚洲性久久影院| 永久网站在线| 99久久中文字幕三级久久日本| 男人的好看免费观看在线视频| 最近最新中文字幕大全电影3| 少妇熟女aⅴ在线视频| 久久国产乱子免费精品| 国产黄色小视频在线观看| 久久久久性生活片| 国产 一区 欧美 日韩| 久久精品国产亚洲av天美| 日韩国内少妇激情av| 男女啪啪激烈高潮av片| 亚洲国产精品专区欧美| 久久精品久久精品一区二区三区| 三级毛片av免费| 亚洲精华国产精华液的使用体验| 22中文网久久字幕| 最近视频中文字幕2019在线8| 在线免费观看的www视频| 18禁在线播放成人免费| a级毛片免费高清观看在线播放| 在现免费观看毛片| 亚洲av中文字字幕乱码综合| 亚洲第一区二区三区不卡| av又黄又爽大尺度在线免费看 | 99久国产av精品| 非洲黑人性xxxx精品又粗又长| 欧美成人午夜免费资源| 色综合色国产| 久久欧美精品欧美久久欧美| 国产一区亚洲一区在线观看| 免费观看人在逋| 白带黄色成豆腐渣| 内射极品少妇av片p| 日本一本二区三区精品| 亚洲电影在线观看av| 欧美精品国产亚洲| 久久久国产成人免费| 麻豆一二三区av精品| 三级国产精品欧美在线观看| 久久精品影院6| 精品久久久久久久末码| 晚上一个人看的免费电影| 黄色日韩在线| 哪个播放器可以免费观看大片| 久久久久久久午夜电影| 少妇人妻精品综合一区二区| 亚洲在线自拍视频| 免费观看在线日韩| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | av卡一久久| 亚洲高清免费不卡视频| 国产69精品久久久久777片| 尾随美女入室| av播播在线观看一区| 少妇裸体淫交视频免费看高清| 免费看av在线观看网站| 国产成人freesex在线| 精品久久久久久久久亚洲| 亚洲欧美日韩卡通动漫| 亚洲欧美精品综合久久99| 岛国毛片在线播放| 最新中文字幕久久久久| 欧美成人一区二区免费高清观看| 又爽又黄a免费视频| 九九在线视频观看精品| 中文欧美无线码| 丝袜美腿在线中文| 精品久久久久久成人av| 国产白丝娇喘喷水9色精品| 精品国产三级普通话版| 亚洲婷婷狠狠爱综合网| 乱码一卡2卡4卡精品| 国内精品美女久久久久久| 国产成人精品婷婷| 最近2019中文字幕mv第一页| 小蜜桃在线观看免费完整版高清| 全区人妻精品视频| 床上黄色一级片| 少妇熟女欧美另类| 国产精品99久久久久久久久| 黄片wwwwww| 天天躁日日操中文字幕| 色尼玛亚洲综合影院| 欧美日韩综合久久久久久| 精品久久久久久电影网 | 少妇丰满av| 老司机影院成人| 国内精品一区二区在线观看| 欧美性猛交╳xxx乱大交人| 日本三级黄在线观看| 国产精品国产三级专区第一集| 丝袜美腿在线中文| 国产又色又爽无遮挡免| 国产黄a三级三级三级人| 亚洲图色成人| eeuss影院久久| 能在线免费观看的黄片| 欧美区成人在线视频| 免费观看人在逋| 99热精品在线国产| 亚洲精品色激情综合| 如何舔出高潮| 国产午夜精品久久久久久一区二区三区| 国产亚洲av嫩草精品影院| 亚洲国产精品合色在线| 美女xxoo啪啪120秒动态图| 国产免费又黄又爽又色| 丰满人妻一区二区三区视频av| 欧美3d第一页| 蜜臀久久99精品久久宅男| 男插女下体视频免费在线播放| 国产成人freesex在线| 国国产精品蜜臀av免费| 久久久久久久亚洲中文字幕| 七月丁香在线播放| 国产久久久一区二区三区| 欧美又色又爽又黄视频| 国产真实乱freesex| 欧美日韩综合久久久久久| 日本免费a在线| a级一级毛片免费在线观看| 久久这里只有精品中国| 少妇熟女欧美另类| 亚洲精品自拍成人| 亚洲国产精品成人久久小说| 国内精品美女久久久久久| 国产成人a∨麻豆精品| 色5月婷婷丁香| 国产成人精品婷婷| 天堂网av新在线| 国内揄拍国产精品人妻在线| 超碰97精品在线观看| 国产高潮美女av| 18禁裸乳无遮挡免费网站照片| 九草在线视频观看| 精品免费久久久久久久清纯| 久久亚洲国产成人精品v| АⅤ资源中文在线天堂| 99久久人妻综合| 亚洲自拍偷在线| 欧美xxxx性猛交bbbb| 精品不卡国产一区二区三区| 日本爱情动作片www.在线观看| 久久久久久久久中文| 亚洲av成人av| 中文字幕久久专区| 高清日韩中文字幕在线| 99久久成人亚洲精品观看| av在线观看视频网站免费| 成人性生交大片免费视频hd| 91久久精品电影网| 亚洲国产精品成人久久小说| 国产高潮美女av| 美女cb高潮喷水在线观看| 亚洲欧美日韩东京热| 少妇丰满av| 男人狂女人下面高潮的视频| av播播在线观看一区| 欧美人与善性xxx| 精品久久久久久久久av| 国产日韩欧美在线精品| 亚洲一级一片aⅴ在线观看| 欧美另类亚洲清纯唯美| 亚洲成人精品中文字幕电影| 国产视频内射| 亚洲精品乱码久久久久久按摩| 老司机影院成人| 亚洲中文字幕日韩| 男人的好看免费观看在线视频| 国语自产精品视频在线第100页| 青青草视频在线视频观看| 赤兔流量卡办理| 国产成人福利小说| 国产伦精品一区二区三区四那| 亚洲18禁久久av| 亚洲天堂国产精品一区在线| 色吧在线观看| 又黄又爽又刺激的免费视频.| 免费播放大片免费观看视频在线观看 | 国产片特级美女逼逼视频| 国产精品永久免费网站| 成人亚洲欧美一区二区av| 国产亚洲午夜精品一区二区久久 | 亚洲精品乱码久久久久久按摩| 亚洲aⅴ乱码一区二区在线播放| 国产免费又黄又爽又色| 国产一区二区在线观看日韩| 中文乱码字字幕精品一区二区三区 | 亚洲最大成人手机在线| 日韩成人伦理影院| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲国产日韩欧美精品在线观看| 色噜噜av男人的天堂激情| 国产精品美女特级片免费视频播放器| 少妇高潮的动态图| 亚洲内射少妇av| 高清毛片免费看| av在线观看视频网站免费| 最近手机中文字幕大全| 日韩成人伦理影院| 免费不卡的大黄色大毛片视频在线观看 | 欧美一区二区国产精品久久精品| 老司机福利观看| 国产精品综合久久久久久久免费| 69av精品久久久久久| 联通29元200g的流量卡| av免费观看日本| 水蜜桃什么品种好| 三级经典国产精品| 久久精品综合一区二区三区| 高清在线视频一区二区三区 | av在线播放精品| 日韩成人伦理影院| 床上黄色一级片| 久久精品国产亚洲av天美| 国产 一区精品| 欧美成人一区二区免费高清观看| 国产亚洲精品av在线| 久久婷婷人人爽人人干人人爱| 亚洲国产色片| 日韩,欧美,国产一区二区三区 | 又爽又黄a免费视频| 亚洲av成人精品一二三区| 精品人妻熟女av久视频| 欧美三级亚洲精品| 久久午夜福利片| 亚洲国产精品合色在线| 观看美女的网站| 国产精品一区www在线观看| 特级一级黄色大片| 我的女老师完整版在线观看| 成人av在线播放网站| 伦理电影大哥的女人| 97在线视频观看| 又爽又黄a免费视频| 亚洲欧美日韩无卡精品| 日韩成人av中文字幕在线观看| 联通29元200g的流量卡| 精品国产三级普通话版| 桃色一区二区三区在线观看| 久久精品91蜜桃| 嫩草影院新地址| 亚洲欧美成人精品一区二区| 亚洲av免费在线观看| 亚洲人成网站在线播| 人人妻人人澡欧美一区二区| 简卡轻食公司| 国产av码专区亚洲av| 最近视频中文字幕2019在线8| 国产一区有黄有色的免费视频 | 纵有疾风起免费观看全集完整版 | 国产成人福利小说| 亚洲av不卡在线观看| 午夜日本视频在线| 日本欧美国产在线视频| 黄色日韩在线| 欧美激情在线99| 欧美日本视频| 国产亚洲91精品色在线| 久久久久久久久久久丰满| 只有这里有精品99| 国产精品久久久久久精品电影小说 | 国产黄a三级三级三级人| 蜜桃久久精品国产亚洲av| 日韩av不卡免费在线播放| 国产人妻一区二区三区在| 欧美潮喷喷水| 国产成人免费观看mmmm| 特大巨黑吊av在线直播| 91久久精品电影网| 久久国内精品自在自线图片| 国产精品一区二区三区四区免费观看| 欧美一区二区国产精品久久精品| 国产大屁股一区二区在线视频| 中文欧美无线码| 欧美色视频一区免费| 久久久久久久久久久免费av| 校园人妻丝袜中文字幕| 看片在线看免费视频| 久久精品国产自在天天线| 日韩一区二区三区影片| ponron亚洲| 1024手机看黄色片| 亚洲国产高清在线一区二区三| 丝袜美腿在线中文| 亚洲最大成人中文| 国产精品人妻久久久影院| 非洲黑人性xxxx精品又粗又长| 国产女主播在线喷水免费视频网站 | 日本色播在线视频| 建设人人有责人人尽责人人享有的 | 97人妻精品一区二区三区麻豆| 在线观看美女被高潮喷水网站| 国产一区有黄有色的免费视频 | 99热6这里只有精品| 精品少妇黑人巨大在线播放 | 国产亚洲av片在线观看秒播厂 | 两个人的视频大全免费| 丰满少妇做爰视频| 99久久九九国产精品国产免费| 中国国产av一级| 亚洲人成网站在线观看播放| 欧美一区二区国产精品久久精品| 啦啦啦观看免费观看视频高清| 丰满人妻一区二区三区视频av| 亚洲国产精品久久男人天堂| 自拍偷自拍亚洲精品老妇| 欧美xxxx性猛交bbbb| 国产免费一级a男人的天堂| 97热精品久久久久久| 免费无遮挡裸体视频| 日日干狠狠操夜夜爽| 搞女人的毛片| 亚洲av成人av| 亚洲av成人精品一区久久| 99热网站在线观看| 亚洲图色成人| 国产精品av视频在线免费观看| 午夜精品国产一区二区电影 | 最后的刺客免费高清国语| www.av在线官网国产| 天美传媒精品一区二区| 日日摸夜夜添夜夜添av毛片| a级毛片免费高清观看在线播放| 中文字幕免费在线视频6| 午夜久久久久精精品| 国产av码专区亚洲av| 亚洲综合精品二区| 精品不卡国产一区二区三区| 午夜亚洲福利在线播放| av在线亚洲专区| 亚洲欧洲日产国产| 国产色婷婷99| 国产精品三级大全| 欧美高清成人免费视频www| 欧美一区二区国产精品久久精品| 小蜜桃在线观看免费完整版高清| 国产视频首页在线观看| 日韩欧美在线乱码| 偷拍熟女少妇极品色| 久久精品国产亚洲av涩爱| 国产一区二区亚洲精品在线观看| 精品久久久噜噜| 哪个播放器可以免费观看大片| 免费看av在线观看网站| 国产午夜精品久久久久久一区二区三区| 色视频www国产| 97人妻精品一区二区三区麻豆| 久热久热在线精品观看| 久久这里只有精品中国| 国产乱人偷精品视频| 精品国产三级普通话版| 久99久视频精品免费| 男人舔奶头视频| 精品无人区乱码1区二区| 亚洲精品久久久久久婷婷小说 | 少妇熟女欧美另类| 国产 一区精品| 免费看日本二区| 久久精品久久久久久噜噜老黄 | 亚洲欧洲国产日韩| 色尼玛亚洲综合影院| 久久99热这里只有精品18| 高清在线视频一区二区三区 | 三级毛片av免费| 麻豆久久精品国产亚洲av| 老司机影院成人| 午夜日本视频在线| 性插视频无遮挡在线免费观看| 人体艺术视频欧美日本| 国产老妇伦熟女老妇高清| 成人毛片60女人毛片免费| 1024手机看黄色片| 中文天堂在线官网| 亚洲人与动物交配视频| 久久久久国产网址| 身体一侧抽搐| 日日摸夜夜添夜夜爱| 不卡视频在线观看欧美| 91久久精品国产一区二区成人| 亚洲色图av天堂| 在线观看美女被高潮喷水网站| 欧美日本视频| 最近中文字幕高清免费大全6| 色综合色国产| 男人的好看免费观看在线视频| 国产精品久久电影中文字幕| 久久精品人妻少妇| 亚洲国产精品成人久久小说| 午夜日本视频在线| 国产精品久久久久久久久免| 久久久久精品久久久久真实原创| 九九爱精品视频在线观看| 亚洲怡红院男人天堂| videos熟女内射| 免费电影在线观看免费观看| 国产精品永久免费网站| 最近最新中文字幕免费大全7| 国产成人a区在线观看| 亚洲18禁久久av| av专区在线播放| 噜噜噜噜噜久久久久久91| 亚洲国产日韩欧美精品在线观看| 天美传媒精品一区二区| 国产精品乱码一区二三区的特点| 亚洲aⅴ乱码一区二区在线播放| 三级男女做爰猛烈吃奶摸视频| 日韩亚洲欧美综合| 亚洲精品自拍成人| 一级av片app| 少妇丰满av| 欧美极品一区二区三区四区| 国产午夜精品久久久久久一区二区三区| 特大巨黑吊av在线直播| 日日撸夜夜添| 美女被艹到高潮喷水动态| 天堂√8在线中文| 18+在线观看网站| 日本av手机在线免费观看| 三级经典国产精品| 2021少妇久久久久久久久久久| 日韩人妻高清精品专区| 婷婷色av中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 国产高清三级在线| 成人特级av手机在线观看| 搞女人的毛片| 欧美xxxx黑人xx丫x性爽| 国产免费一级a男人的天堂| 少妇熟女欧美另类| eeuss影院久久| 国产又黄又爽又无遮挡在线| 国产一区二区在线av高清观看| 日韩欧美在线乱码| 亚洲欧美日韩卡通动漫| 九九在线视频观看精品| 乱系列少妇在线播放| 只有这里有精品99| 精品久久久久久成人av| 97超视频在线观看视频| 国产色爽女视频免费观看| 少妇的逼好多水| 久久亚洲国产成人精品v| 91精品伊人久久大香线蕉| av在线蜜桃| 久久精品夜色国产| 亚洲真实伦在线观看| 麻豆av噜噜一区二区三区| 热99在线观看视频| 真实男女啪啪啪动态图| 国产国拍精品亚洲av在线观看| 国产免费男女视频| 最后的刺客免费高清国语| 亚洲精品aⅴ在线观看| a级毛色黄片| 少妇丰满av| 国产乱人视频| 人妻系列 视频| 日本与韩国留学比较|