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

    不同射流狀態(tài)下射流氣體與高超聲速主流相互作用影響

    2020-07-10 07:42:46王麗燕檀妹靜蔣云淞王振峰
    關(guān)鍵詞:馬赫數(shù)超聲速熱流

    王麗燕,檀妹靜,聶 亮,蔣云淞,袁 野,王振峰

    (中國運(yùn)載火箭技術(shù)研究院 空間物理重點(diǎn)實(shí)驗(yàn)室, 北京 100076)

    嚴(yán)重的氣動(dòng)加熱是高超聲速飛行器無法避免而又必須妥善解決的一個(gè)重大技術(shù)難題,稱為“熱障”[1],目前通常采用被動(dòng)式防熱設(shè)計(jì)來解決高超聲速飛行器的這一問題。被動(dòng)式熱防護(hù)系統(tǒng)提高防熱或隔熱能力的主要方法為提高材料防熱或隔熱性能和增加材料厚度等,但是新型材料的研制周期較長,材料厚度的增加又無疑增加了熱防護(hù)系統(tǒng)的重量,這使得單一的被動(dòng)式熱防護(hù)系統(tǒng)難以完全應(yīng)對(duì)未來高超聲速飛行器長時(shí)間、高熱流和可重復(fù)使用的服役環(huán)境。主動(dòng)式熱防護(hù)系統(tǒng)具有冷卻能力強(qiáng)、可抵抗高熱流密度的長時(shí)間加熱、可重復(fù)使用等優(yōu)點(diǎn),可以滿足未來高超聲速飛行器的熱防護(hù)需求。主動(dòng)式熱防護(hù)系統(tǒng)防熱方法主要有發(fā)汗冷卻、膜冷卻、對(duì)流冷卻和噴霧冷卻。對(duì)于發(fā)汗冷卻和膜冷卻來說,飛行器表面引射氣體會(huì)與高超聲速主流來流發(fā)生相互作用,產(chǎn)生復(fù)雜的激波和漩渦剪切流動(dòng),這使得高超聲速飛行器流場(chǎng)內(nèi)流動(dòng)更為復(fù)雜。引射效應(yīng)一直以來都是國內(nèi)外學(xué)者研究的熱點(diǎn)[2-5]。

    Heufer和Olivier等[6-8]通過實(shí)驗(yàn)和數(shù)值模擬對(duì)高超聲速主流中的楔形體狹縫式氣膜冷卻進(jìn)行了研究,表明氣膜冷卻在高超聲速流動(dòng)中能有效減少機(jī)體的熱負(fù)荷;隨著吹風(fēng)比、主流雷諾及馬赫數(shù)的增大,氣膜冷卻效率提高。Wimberly等[9]針對(duì)細(xì)圓錐體模型在來流馬赫數(shù)12~17條件下進(jìn)行了風(fēng)洞實(shí)驗(yàn),結(jié)果表明壁面換熱系數(shù)及流阻與吹風(fēng)比呈負(fù)相關(guān),吹風(fēng)比增大,二者隨之減小。熊宴斌[10]在馬赫數(shù)3的風(fēng)洞中研究分析了青銅、不銹鋼粉末燒結(jié)多孔平板、燒結(jié)金屬絲網(wǎng)多孔平板以及不銹鋼燒結(jié)粉末多孔曲面結(jié)構(gòu)的引射規(guī)律,研究表明,二次流注入能夠有效減小壁面速度梯度,降低壁面溫度。數(shù)值模擬研究方面,Yang等[11]針對(duì)槽縫式平板進(jìn)行了不同流態(tài)下不同影響因素(射流率、縫高、吹風(fēng)比)狀態(tài)下的模擬研究,發(fā)現(xiàn)氣膜冷卻是壁面進(jìn)行主動(dòng)熱防護(hù)的一種有效途徑,在高超聲速流動(dòng)中降溫減阻效果更加明顯。Konopka等[12-13]采用大渦模擬方法對(duì)氣膜冷卻問題進(jìn)行了研究,發(fā)現(xiàn)順壓梯度有利于降低剪切層湍流度并推遲轉(zhuǎn)捩,從而提高冷卻效率;湍流流態(tài)下氣膜冷卻的整體效率比層流流態(tài)下的稍低,冷卻長度稍短[13]。王建[14]、余磊[15]等的研究表明氣體引射進(jìn)入邊界層使邊界層增厚,有效降低了壁面摩擦系數(shù);隨著引射氣體流量增加,壁面溫度和對(duì)流換熱系數(shù)顯著減小。另外,還有學(xué)者使用直接數(shù)值模擬[16-18]對(duì)壁面存在二次流注入的流動(dòng)傳熱進(jìn)行計(jì)算,獲得了射流與主流相互作用的精細(xì)流動(dòng)圖畫[19]。上述研究豐富了人們對(duì)氣體引射效應(yīng)的認(rèn)識(shí)。工程上高超聲速飛行器常出現(xiàn)局部熱流過高的情況,氣體引射效應(yīng)能有效解決這一狀況,因此系統(tǒng)研究氣體引射與高超聲速主流相互作用有重要意義。

    本文針對(duì)方孔和圓孔橫向引射模型、高超聲速來流條件下冷卻氣體射流問題,采用數(shù)值模擬的方法,開展氣體射流條件下的高超聲速飛行器表面氣動(dòng)加熱機(jī)理研究,分析射流氣體壓強(qiáng)、射流速度以及射流方向等因素對(duì)壁面溫度和附面層內(nèi)流動(dòng)結(jié)構(gòu)的影響,研究氣體射流表面降熱特性,并根據(jù)模擬結(jié)果對(duì)高溫邊界層流場(chǎng)進(jìn)行深入分析。

    1 計(jì)算方法與模型

    1.1 計(jì)算模型

    圖1 物理模型Fig.1 Physical model

    圖1給出了本文研究模型,其中高超聲速主流從左側(cè)進(jìn)入,射流氣體與壁面成一定角度從射流孔射出,與主流相互作用,向下游發(fā)展。x、y、z分別為主流的流向、法向和展向。射流孔(包括圓孔和方孔)中心取為坐標(biāo)原點(diǎn),主流進(jìn)口位于x=-20 mm處,出口位于x=80 mm處,射流孔出口位于y=-10 mm處??臻g流場(chǎng)高度Ly=30 mm,流場(chǎng)寬度Lz=6 mm。為了便于分析引射表面與光滑表面流動(dòng)特性、降低計(jì)算誤差同時(shí)保證數(shù)值計(jì)算結(jié)果是在同一流場(chǎng)狀態(tài)下得到,將引射表面對(duì)應(yīng)的壁面作為光滑壁面[19]。考慮主流和射流流動(dòng)特性,流動(dòng)為層流,控制方程由經(jīng)典的Navier-Stokes方程組來描述,見式(1)~(3):

    (1)

    (2)

    (3)

    1.2 網(wǎng)格劃分及邊界條件設(shè)置

    對(duì)于高超聲速主流來流條件下的氣體引射流動(dòng)的計(jì)算流體動(dòng)力學(xué)(Computational Fluid Dynamics, CFD)模擬,網(wǎng)格是至關(guān)重要的。本文開展了包含基準(zhǔn)網(wǎng)格在內(nèi)的5個(gè)網(wǎng)格狀態(tài)的計(jì)算,以考察第一層網(wǎng)格高度、網(wǎng)格數(shù)量對(duì)熱環(huán)境計(jì)算結(jié)果的影響。壁面第一層網(wǎng)格法向距離計(jì)算公式為:

    (4)

    式中:y+為第一層網(wǎng)格到壁面的無量綱距離,取y+=1;uτ為壁面摩擦系數(shù)。

    為排除網(wǎng)格數(shù)量的影響,分別對(duì)240萬、140萬和410萬網(wǎng)格進(jìn)行計(jì)算。240萬與410萬網(wǎng)格數(shù)量下壁面中心線溫度相對(duì)接近而與140萬網(wǎng)格數(shù)量下壁面溫度存在一定差異,即在240萬基準(zhǔn)網(wǎng)格基礎(chǔ)上增加網(wǎng)格數(shù)量不會(huì)引起壁面溫度的較大變化。為提高計(jì)算效率,網(wǎng)格設(shè)置為第一層網(wǎng)格間距0.001 mm、總網(wǎng)格量為240萬的分塊對(duì)接網(wǎng)格,如圖2所示。求解控制方程采用有限體積法離散,通量離散采用總變差減小格式以及MinMod限制器,保證空間為二階精度。時(shí)間離散采用點(diǎn)隱式方法[19]。

    圖2 計(jì)算區(qū)域網(wǎng)格Fig.2 Grid of computational domain

    邊界條件設(shè)置:光滑壁面及射流區(qū)域壁面為無滑移絕熱壁面條件;主流區(qū)兩側(cè)為對(duì)稱邊界條件;主流入口為高超聲速入口,給定壓力、速度和溫度條件;射流氣體入口為質(zhì)量流率入口,給定質(zhì)量流率和溫度條件,氣體溫度為221 K;主流和引射氣體工質(zhì)均為空氣;出口超聲速,所有邊界條件外推;初始條件采用壓力值為2550 Pa。

    1.3 方法驗(yàn)證

    為了驗(yàn)證計(jì)算軟件氣動(dòng)熱模擬的準(zhǔn)確性,采用二維鈍頭圓柱作為驗(yàn)證算例,該算例包含Holden的風(fēng)洞試驗(yàn)數(shù)據(jù)[20]。其中圓柱半徑為0.038 1 m,來流馬赫數(shù)為8,靜壓為855 Pa,靜溫為125.07 K,壁面溫度為294 K。計(jì)算網(wǎng)格如圖3所示,網(wǎng)格量100×80(周向×徑向),第一層網(wǎng)格間距3.81×10-6m。

    圖3 二維圓柱計(jì)算網(wǎng)格Fig.3 Two-dimensional grid of computational domain

    圖4給出了圓柱表面斯坦頓數(shù)分布與風(fēng)洞實(shí)驗(yàn)值比較,從計(jì)算結(jié)果來看熱流誤差均在10%以內(nèi),證明本方法計(jì)算熱流的精度是可靠的。

    圖4 圓柱表面斯坦頓數(shù)分布曲線Fig.4 Stanton number distribution at cylinder surface

    2 結(jié)果與分析

    考慮氣體引射與主流摻混后流場(chǎng)結(jié)構(gòu)與引射氣體壓強(qiáng)、引射速度和射流方向等的變化密切相關(guān),本文以主流0°攻角,主流來流靜溫221 K、靜壓2550 Pa、馬赫數(shù)6,射流來流方向90°、靜溫221 K、靜壓19 125 Pa、馬赫數(shù)1.2為基準(zhǔn)狀態(tài),通過改變射流壓強(qiáng)、射流速度和射流方向,得到了摻混后流場(chǎng)結(jié)構(gòu)、壁面溫度分布和壁面中心線溫度等參數(shù),并對(duì)其變化規(guī)律進(jìn)行分析。

    2.1 流場(chǎng)特征

    圖5展示了基準(zhǔn)狀態(tài)下氣體引射的基本流動(dòng)結(jié)構(gòu)。由圖可知:一方面,自入口處邊界層開始發(fā)展,主流首先產(chǎn)生一道斜激波,而到了射流孔附近,受射流垂直進(jìn)入主流影響,受到阻礙的主流產(chǎn)生一道弓形激波;另一方面,射流從方孔進(jìn)入主流,在射流出口附近形成明顯的馬赫盤。斜激波和弓形激波下面存在明顯的速度剪切層,射流垂直進(jìn)入主流,受到該剪切層的阻礙,在靠近剪切層位置附近滯止,向四周流動(dòng)。流向上游的流體產(chǎn)生回流渦,沿展向左右流動(dòng)的流體,在x截面形成一對(duì)“腎型”渦,其余大部分流體則流向下游。

    圖5 氣體引射基本流動(dòng)結(jié)構(gòu)Fig.5 Basic flow structure in the gas ejection

    圖6給出了主流壁面溫度分布及射流孔附近表面極限流線。由圖可知,入口附近流動(dòng)滯止,該位置附近的壁溫極高。射流孔前方某位置處壁溫迅速降低,因?yàn)樵撐恢么嬖谝粋€(gè)較大的回流區(qū),一直延伸到高溫的滯止區(qū)域,這與空間流線顯示出的射流孔上游回流渦相對(duì)應(yīng)。圍繞射流孔一周直到遠(yuǎn)下游,存在明顯的“腎型”渦流動(dòng)結(jié)構(gòu),而“腎型”渦流動(dòng)結(jié)構(gòu)與壁面的低溫區(qū)相對(duì)應(yīng)。

    圖6 主流壁面溫度分布及射流孔附近表面極限流線Fig.6 Wall temperature distribution and the limit streamline near the surface of ejecting hole

    2.2 射流壓強(qiáng)影響

    為了研究射流壓強(qiáng)對(duì)流動(dòng)結(jié)構(gòu)和壁面熱流的影響,開展了包含基準(zhǔn)狀態(tài)在內(nèi)的5個(gè)不同射流壓強(qiáng)狀態(tài)的模擬。其中主流來流靜溫221 K,馬赫數(shù)6,靜壓2550 Pa。射流靜溫221 K,馬赫數(shù)1.2,靜壓分別為:2550 Pa、5100 Pa、10 200 Pa、19 125 Pa和40 800 Pa。

    圖7為不同射流壓強(qiáng)下z=0 mm截面和x=0 mm截面流線及馬赫數(shù)分布。由圖可知,射流對(duì)橫向來流產(chǎn)生阻擋,形成一道弓形激波,弓形激波使迎風(fēng)面近壁面發(fā)生流動(dòng)分離1,隨著壓強(qiáng)增大,弓形激波強(qiáng)度、分離程度均增強(qiáng)。此外,射流與來流之間形成明顯邊界。射流下游背風(fēng)區(qū)的低壓特性會(huì)造成此處近壁面產(chǎn)生流動(dòng)分離2,隨著壓強(qiáng)增大,此處流動(dòng)分離加劇。同時(shí)射流與來流剪切摻混形成的正反向旋轉(zhuǎn)流向“腎型”渦,隨著壓強(qiáng)增大,靠近射流孔的“腎型”渦增大,較遠(yuǎn)的一對(duì)“腎型”渦逐漸消失。因?yàn)椤澳I型”渦的產(chǎn)生主要依靠射流與主流在邊界處的剪切摻混,而射流壓強(qiáng)增大使邊界范圍變大,速度剪切層向遠(yuǎn)離壁面方向移動(dòng),導(dǎo)致渦的攪動(dòng)范圍增大。

    圖7 不同射流壓強(qiáng)下流線及馬赫數(shù)分布Fig.7 Distribution of streamline and Mach number under different ejecting pressure

    圖8為不同射流壓強(qiáng)下壁面中心線溫度變化。由圖可知,射流孔附近低溫區(qū)降溫效果較強(qiáng)。因?yàn)閴簭?qiáng)較小時(shí),引射氣體對(duì)主流的阻礙作用較小,對(duì)上下游影響范圍較小,引射氣體被主流壓向壁面,在壁面很近的區(qū)域內(nèi)向下游流動(dòng)。隨著射流壓強(qiáng)增大,射流孔上游回流渦增大,上游壁面低溫區(qū)范圍增大,但低溫區(qū)內(nèi)降溫效果減弱;射流孔下游低溫區(qū)范圍也增大,靠近射流孔處的壁面溫度升高,遠(yuǎn)離射流孔處的壁面溫度降低,即增大射流壓強(qiáng),能增大下游低溫區(qū)范圍,增強(qiáng)下游整體降溫效果,即單位質(zhì)量的冷卻介質(zhì)冷卻效果最好。若射流孔數(shù)目增多,則低壓射流對(duì)壁面降溫效果更優(yōu)。

    圖8 不同射流壓強(qiáng)下壁面中心線溫度對(duì)比Fig.8 Comparison of wall temperature of center line under different ejecting pressure

    2.3 射流方向影響

    為了研究射流方向?qū)α鲃?dòng)結(jié)構(gòu)和壁面熱流的影響,開展了不同射流孔與主流壁面夾角研究。其中射流方向與主流方向的夾角分別為30°、45°、60°、90°、120°、135°和150°。

    圖9給出了不同射流孔方向下流線及馬赫數(shù)分布。由圖9(a)可知,在截面z=0 mm處,射流角度越大,射流與主流作用產(chǎn)生的弓形激波越大,射流孔上下游區(qū)域的回流區(qū)越大。射流角度為銳角時(shí),射流對(duì)上游的影響范圍隨角度增大,下游的影響范圍也隨角度增大,但強(qiáng)度減弱。引射角度為鈍角時(shí),射流對(duì)上游的影響范圍隨角度增大,強(qiáng)度也增大,對(duì)下游影響則相反。由圖9(b)可知,隨著射流角度增大,射流孔兩側(cè)的“腎型”渦增大。射流角度為銳角時(shí),角度越大,“腎型”渦越貼近壁面。射流角度為鈍角時(shí),“腎型”渦產(chǎn)生位置較銳角時(shí)遠(yuǎn),且角度越大,產(chǎn)生位置距壁面越遠(yuǎn)。由于射流角度的存在,引射氣體沿主流方向的速度分量導(dǎo)致主流流場(chǎng)較基準(zhǔn)流場(chǎng)有所變化,引射角度為銳角時(shí)沿主流方向的速度分量增強(qiáng),鈍角時(shí)則相反。

    圖9 不同射流孔方向下流線及馬赫數(shù)分布Fig.9 Distribution of streamline and Mach number under different ejecting direction

    圖10為壁面中心線溫度對(duì)比。由圖可知,射流方向與主流方向夾角為銳角時(shí),射流孔上游降溫效果劣于基準(zhǔn)態(tài);下游降溫效果差別較小,都優(yōu)于基準(zhǔn)態(tài)。隨著射流方向與主流方向夾角增大,上游低溫區(qū)增大,降溫效果增強(qiáng)。射流方向與主流方向夾角為鈍角時(shí),射流上游降溫效果優(yōu)于基準(zhǔn)態(tài),下游則劣于基準(zhǔn)態(tài)。隨著角度增大,射流孔下游降溫效果減弱;射流方向與主流方向夾角為銳角時(shí)降溫范圍大于鈍角,效果也優(yōu)于鈍角時(shí)狀態(tài)。角度越小,射流越貼近壁面,射流孔附近降溫效果越好,射流角為30°時(shí)較射流角為90°時(shí)壁面溫度降低約30%。

    圖10 不同射流方向下壁面中心線溫度對(duì)比Fig.10 Comparison of the wall temperature of the center line under different ejecting direction

    2.4 射流孔形狀影響

    為了研究射流孔形狀對(duì)流動(dòng)結(jié)構(gòu)和壁面熱流的影響,研究了方孔和圓孔兩種射流孔。

    圖11給出了不同射流孔形狀下流線及馬赫數(shù)分布。由圖可知,相同條件下,射流孔形狀的改變對(duì)流場(chǎng)結(jié)構(gòu)影響不大。相同工況不同引射孔情況下,弓形激波大小位置、分流區(qū)大小、強(qiáng)度基本相同。可見,孔結(jié)構(gòu)的改變,對(duì)射流與主流摻混作用影響較小。由圖12不同射流孔形狀下壁面中心線溫度變化可知,相同工況下,圓孔的降溫效果略優(yōu)于方孔。因?yàn)榭捉Y(jié)構(gòu)的改變,對(duì)射流與主流摻混作用影響較小,所以對(duì)流場(chǎng)的降溫效果影響也較小。

    圖11 不同射流孔形狀下z=0 mm截面流線及馬赫數(shù)分布Fig.11 Distribution of streamline and Mach number at z=0 mm under different ejecting holes

    圖12 不同射流孔形狀下壁面中心線溫度對(duì)比Fig.12 Comparison of the wall temperature of the center line under different ejecting holes

    2.5 射流速度影響

    2.5.1 高速射流

    以方孔模型為研究對(duì)象,研究射流來流馬赫數(shù)為2.4、1.2、0.6時(shí)的流場(chǎng)狀態(tài)。通過對(duì)比分析,研究高速射流速度情況對(duì)流動(dòng)及冷卻效果的影響。

    圖13給出了不同射流速度下z=0 mm截面和x=0 mm截面流線及馬赫數(shù)分布。隨著射流速度增大,射流對(duì)主流的阻礙作用增大,產(chǎn)生的弓形激波也較強(qiáng),使主流遠(yuǎn)離壁面。隨著射流速度增大,射流孔上游和下游流動(dòng)分離加劇,射流孔兩側(cè)的“腎型”渦增大。

    圖14為不同射流速度下壁面中心線溫度變化。由圖可知,隨著射流速度增大,射流孔上游回流渦增大,上游壁面低溫區(qū)范圍增大,但低溫區(qū)內(nèi)降溫效果減弱;引射下游靠近射流孔處的壁面溫度升高,遠(yuǎn)離射流孔處的壁面溫度先降低后增大,推測(cè)可能存在一個(gè)最優(yōu)射流速度,使下游壁面降溫效果最佳。

    圖13 不同射流速度下流線及馬赫數(shù)分布Fig.13 Distribution of streamline and Mach number under different ejecting velocity

    圖14 不同射流速度下壁面中心線溫度對(duì)比Fig.14 Comparison of the wall temperature of the center line under different ejecting velocity

    2.5.2 中低速射流

    以圓孔模型為研究對(duì)象,研究射流來流馬赫數(shù)為1.2、0.6、0.3、0.1、0.05時(shí)的流場(chǎng)狀態(tài),壁面取300 K等溫壁面邊界條件。通過對(duì)比分析,研究中低速射流速度情況對(duì)流動(dòng)及冷卻效果的影響。

    圖15給出了不同射流速度情況下射流孔附近流動(dòng)對(duì)比,可以看出,速度越低,其對(duì)主流的改變?cè)叫 I淞黢R赫數(shù)為0.05(速度約15 m·s-1)時(shí),流體從射流孔流出之后幾乎是貼壁繼續(xù)向后流動(dòng),對(duì)主流阻礙較小,幾乎不產(chǎn)生弓形激波。而射流馬赫數(shù)為1.2時(shí),在射流孔附近區(qū)域形成較強(qiáng)的弓形激波,在其上游形成一個(gè)較大的回流渦,受回流渦的影響,在上游形成高熱流區(qū)域,如圖16所示。

    圖15 不同射流速度下射流孔附近流動(dòng)對(duì)比Fig.15 Comparison of flow structure near the ejecting hole under different ejecting velocity

    圖16給出了中心線上的熱流分布,可以看出,在射流的下游區(qū)域熱流均有所降低。射流馬赫數(shù)小于0.3時(shí),射流對(duì)主流影響較小,基本集中在射流孔附近,冷卻效果較弱。射流馬赫數(shù)大于等于0.3時(shí),下游區(qū)域熱流顯著降低。射流馬赫數(shù)為0.6時(shí)在射流孔附近的降溫效果優(yōu)于馬赫數(shù)為1.2時(shí)的,而下游稍遠(yuǎn)區(qū)域冷卻效果基本相當(dāng)。但射流馬赫數(shù)為1.2時(shí),射流孔上游區(qū)域形成不利干擾,使得熱流不降反增,可以推測(cè),存在最佳的射流速度,使冷卻效果最好。就本文研究工況下,單孔射流馬赫數(shù)為0.6左右,冷卻效果最好。射流的影響呈現(xiàn)明顯的三維效應(yīng),其影響范圍明顯受射流速度的影響,局部呈現(xiàn)條帶狀的熱流分布,因此,若要大范圍地冷卻,則需要更多的射流孔出流進(jìn)行冷卻。

    圖16 不同射流速度下壁面中心線熱流對(duì)比Fig.16 Comparison of the wall heat flux of the center line under different ejecting velocity

    2.5.3 多孔射流影響

    進(jìn)一步對(duì)多孔(壁面整體)射流情況進(jìn)行研究,其中,主流來流0°攻角、靜溫221 K、靜壓2550 Pa、馬赫數(shù)6,射流氣體靜溫221 K、靜壓19 125 Pa,射流速度為1 m·s-1、2 m·s-1、15 m·s-1、30 m·s-1。壁面為等溫壁300 K。

    圖17給出了不同射流速度的流動(dòng)對(duì)比,可以看出,增加射流速度后,會(huì)引起邊界層流動(dòng)的明顯改變。低速射流情況下,射流氣體直接進(jìn)入邊界層貼壁流動(dòng),使得熱流明顯降低。高速射流情況下,射流氣體將原平板邊界層“抬離”壁面,使得壁面附近充滿射流氣體,從而使得壁面熱流接近于0。

    圖17 不同射流速度下流動(dòng)對(duì)比Fig.17 Comparison of flow structure under different ejecting velocity

    圖18給出了不同射流速度下壁面熱流對(duì)比,可以看出,考慮壁面整體射流之后,熱流大幅降低。1 m·s-1的壁面射流速度即使得壁面熱流降低接近三分之二,2 m·s-1的壁面射流速度即使得壁面熱流降低接近90%,而15 m·s-1的射流速度即使得壁面熱流降為接近0的量值,由此可見,壁面射流即使射流速度很低,降溫冷卻效果也非常顯著。

    圖18 不同射流速度下壁面中心線熱流對(duì)比Fig.18 Comparison of the wall heat flux of thecenter line under different ejecting velocity

    2.5.4 小結(jié)

    1)高速射流的冷卻機(jī)理:高速射流一般用于小孔噴流,當(dāng)冷卻工質(zhì)以高速狀態(tài)(特別是Ma>1)噴出時(shí),噴流氣體形成類似“氣柱”的一團(tuán)流體,來流和氣柱相互作用,在氣柱前方形成渦系流動(dòng),使得局部熱流降低或升高,而在氣柱后方,來流將氣柱壓在物體表面,使得后方壁面得到冷卻,如圖19(a)所示;當(dāng)氣流速度進(jìn)一步增大,氣流對(duì)主流的沖破作用加強(qiáng),對(duì)流動(dòng)結(jié)構(gòu)的改變也更加明顯,如圖19(b)所示。實(shí)際上,噴流和主流相互作用的流動(dòng)結(jié)構(gòu)很大程度上與兩種流動(dòng)的靜壓、動(dòng)壓兩種關(guān)系的匹配相關(guān)。這種高速噴流對(duì)流動(dòng)結(jié)構(gòu)的改變大,需要的冷卻工質(zhì)也多,對(duì)于可以短時(shí)改變流動(dòng)結(jié)構(gòu)從而降低熱流的情況較為適用。

    圖19 高速射流與主流相互作用的流動(dòng)結(jié)構(gòu)Fig.19 Flow structure of high velocity ejecting interaction with mainstream

    2)低速射流的冷卻機(jī)理:低速射流可用于小孔噴流和多孔介質(zhì)引射,冷卻工質(zhì)以低速噴出時(shí),噴流在主流的作用下被壓在物體表面,在邊界層內(nèi)貼壁流動(dòng),進(jìn)而冷卻飛行器外表面。由于具有需要冷卻工質(zhì)相對(duì)較少和對(duì)主流流動(dòng)改變較少的特點(diǎn),這種冷卻方式的應(yīng)用范圍較廣,并且可較長時(shí)間持續(xù)工作。綜合考慮實(shí)際情況,可對(duì)孔的大小、孔的分布數(shù)量、流速大小等因素進(jìn)行綜合優(yōu)化設(shè)計(jì),選擇整體效果最好的方案。

    3 結(jié)論

    1)高超聲速來流條件下的射流流場(chǎng)中存在主流斜激波、弓形激波、馬赫盤、速度剪切層、流向回流渦、展向“腎型”渦等復(fù)雜流動(dòng)結(jié)構(gòu),流動(dòng)結(jié)構(gòu)的改變和射流的壓強(qiáng)、速度、方向相關(guān),其中射流速度對(duì)流動(dòng)結(jié)構(gòu)的改變尤為明顯。

    2)射流壓強(qiáng)增大,射流孔下游低溫區(qū)范圍增大,靠近射流孔處的壁面溫度升高,遠(yuǎn)離射流孔處的壁面溫度降低,增強(qiáng)了下游整體降溫效果,若射流孔數(shù)目增多,則低壓射流對(duì)壁面降溫效果更優(yōu)。射流方向與主流方向夾角為銳角時(shí),利于射流孔下游降溫,角度越小,射流越貼近壁面,射流孔附近降溫效果越好,射流角為30°時(shí)較射流角為90°時(shí)壁面溫度降低約30%。夾角為鈍角時(shí),利于射流孔上游降溫。圓孔的降溫效果略優(yōu)于方孔。

    3)在高速(Ma>1)射流情況下,一定范圍內(nèi),壓強(qiáng)和速度適當(dāng)增大,均可使射流下游的冷卻效果加強(qiáng)。當(dāng)射流速度達(dá)到一定量值后,主流流場(chǎng)結(jié)構(gòu)明顯改變,流動(dòng)結(jié)構(gòu)改變大,需要的冷卻工質(zhì)多,較為適用于可以短時(shí)改變流動(dòng)結(jié)構(gòu)從而降低熱流的情況。

    4)在中低速(Ma<0.6)射流情況下,射流對(duì)流場(chǎng)的改變相對(duì)較小,低速時(shí),射流基本上不改變主流流場(chǎng)而在邊界層內(nèi)流動(dòng),達(dá)到冷卻壁面的效果;在這個(gè)速度區(qū)間范圍內(nèi),流速越大,冷卻范圍越大,冷卻效果也相對(duì)較好。壁面整體引射時(shí), 1 m·s-1的壁面引射速度即使得壁面熱流降低接近三分之二,2 m·s-1的壁面引射速度即使得壁面熱流降低接近90%,而15 m·s-1的引射速度即使得壁面熱流降為接近0的量值。

    猜你喜歡
    馬赫數(shù)超聲速熱流
    高馬赫數(shù)激波作用下單模界面的Richtmyer-Meshkov不穩(wěn)定性數(shù)值模擬
    爆炸與沖擊(2024年7期)2024-11-01 00:00:00
    高超聲速出版工程
    高超聲速飛行器
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    載荷分布對(duì)可控?cái)U(kuò)散葉型性能的影響
    超聲速旅行
    內(nèi)傾斜護(hù)幫結(jié)構(gòu)控釋注水漏斗熱流道注塑模具
    空調(diào)溫控器上蓋熱流道注塑模具設(shè)計(jì)
    聚合物微型零件的熱流固耦合變形特性
    中國塑料(2017年2期)2017-05-17 06:13:24
    透明殼蓋側(cè)抽模熱流道系統(tǒng)的設(shè)計(jì)
    中國塑料(2014年5期)2014-10-17 03:02:17
    国产高清视频在线观看网站| 国产精品亚洲美女久久久| 插逼视频在线观看| 精品一区二区免费观看| 亚洲精品久久国产高清桃花| 又爽又黄a免费视频| 伊人久久精品亚洲午夜| 搡老岳熟女国产| 国产精品一区二区三区四区久久| 美女cb高潮喷水在线观看| 高清毛片免费看| 熟妇人妻久久中文字幕3abv| av国产免费在线观看| 乱人视频在线观看| 少妇熟女欧美另类| 欧美激情国产日韩精品一区| 一a级毛片在线观看| 听说在线观看完整版免费高清| 国产真实乱freesex| 久久久久久久亚洲中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 在线播放无遮挡| 老熟妇仑乱视频hdxx| 精品午夜福利在线看| 久久久久国内视频| 成人美女网站在线观看视频| av.在线天堂| 我的女老师完整版在线观看| 亚洲va在线va天堂va国产| 亚洲av电影不卡..在线观看| 欧美日韩综合久久久久久| 午夜a级毛片| 国产高清有码在线观看视频| 九九爱精品视频在线观看| 亚洲国产精品合色在线| 国产真实伦视频高清在线观看| 高清日韩中文字幕在线| 日日撸夜夜添| 国产高清有码在线观看视频| 搡老妇女老女人老熟妇| 亚洲高清免费不卡视频| 俺也久久电影网| 午夜爱爱视频在线播放| aaaaa片日本免费| 99久国产av精品| 久久精品国产亚洲av天美| 免费观看精品视频网站| 国产精品电影一区二区三区| 久久久久久久久中文| 免费看光身美女| 国产视频一区二区在线看| 两个人的视频大全免费| 免费看a级黄色片| 久久99热这里只有精品18| 女生性感内裤真人,穿戴方法视频| 99国产精品一区二区蜜桃av| 精品午夜福利视频在线观看一区| 老熟妇乱子伦视频在线观看| 日韩欧美三级三区| 日韩高清综合在线| 真实男女啪啪啪动态图| 亚洲国产高清在线一区二区三| 一a级毛片在线观看| 夜夜爽天天搞| 日韩欧美国产在线观看| 免费在线观看影片大全网站| 露出奶头的视频| 美女 人体艺术 gogo| 天堂av国产一区二区熟女人妻| 成人毛片a级毛片在线播放| 国产亚洲91精品色在线| 国产欧美日韩一区二区精品| 看黄色毛片网站| 此物有八面人人有两片| 一夜夜www| 午夜福利成人在线免费观看| 内射极品少妇av片p| 国产一级毛片七仙女欲春2| 欧美日韩在线观看h| 91麻豆精品激情在线观看国产| 国产一区二区亚洲精品在线观看| 欧美一区二区亚洲| 六月丁香七月| 亚洲欧美中文字幕日韩二区| 成年女人毛片免费观看观看9| 在线播放国产精品三级| 美女被艹到高潮喷水动态| 久久精品国产亚洲av涩爱 | 看免费成人av毛片| 夜夜看夜夜爽夜夜摸| 神马国产精品三级电影在线观看| 欧美bdsm另类| 成人特级av手机在线观看| 六月丁香七月| 欧美色欧美亚洲另类二区| 国产精品一区二区免费欧美| 亚洲国产精品合色在线| 无遮挡黄片免费观看| 久久久精品94久久精品| 午夜精品在线福利| 尤物成人国产欧美一区二区三区| 亚洲国产精品合色在线| 亚洲乱码一区二区免费版| 有码 亚洲区| 亚洲欧美日韩高清专用| 尤物成人国产欧美一区二区三区| 亚洲欧美日韩东京热| 嫩草影视91久久| 亚洲国产欧洲综合997久久,| videossex国产| 成年av动漫网址| 亚洲国产精品合色在线| 亚洲国产精品sss在线观看| 国产精品三级大全| 午夜精品在线福利| 久久久久久久久久久丰满| 久久鲁丝午夜福利片| 神马国产精品三级电影在线观看| 少妇的逼好多水| 少妇猛男粗大的猛烈进出视频 | 久久草成人影院| 成人高潮视频无遮挡免费网站| 国产不卡一卡二| 插阴视频在线观看视频| 免费av观看视频| 国产 一区精品| 亚洲国产精品成人综合色| 久久久久久久久久久丰满| 日本-黄色视频高清免费观看| 九色成人免费人妻av| 中文字幕精品亚洲无线码一区| 国产精品99久久久久久久久| 菩萨蛮人人尽说江南好唐韦庄 | 精品久久久久久久久av| 成年女人毛片免费观看观看9| 人人妻人人澡人人爽人人夜夜 | 亚洲图色成人| 日本黄大片高清| 久久久精品94久久精品| 国产精品永久免费网站| 日韩欧美 国产精品| 欧美性感艳星| 男女那种视频在线观看| 又粗又爽又猛毛片免费看| 少妇裸体淫交视频免费看高清| 亚洲国产精品久久男人天堂| 亚洲一区二区三区色噜噜| 国产私拍福利视频在线观看| 69av精品久久久久久| 色视频www国产| 亚洲在线观看片| 欧美激情久久久久久爽电影| 日韩一本色道免费dvd| 九九久久精品国产亚洲av麻豆| 男女之事视频高清在线观看| ponron亚洲| 日韩欧美一区二区三区在线观看| 欧美一区二区亚洲| 亚洲成人久久爱视频| 免费观看精品视频网站| 午夜福利在线观看免费完整高清在 | 在线天堂最新版资源| eeuss影院久久| 长腿黑丝高跟| 听说在线观看完整版免费高清| videossex国产| 久久久国产成人精品二区| 精品久久久久久久人妻蜜臀av| 国产午夜福利久久久久久| avwww免费| 狠狠狠狠99中文字幕| 老司机午夜福利在线观看视频| 国产精品野战在线观看| 免费在线观看影片大全网站| 亚洲精品色激情综合| 久久久久九九精品影院| 一本精品99久久精品77| 永久网站在线| 丰满的人妻完整版| 欧美xxxx黑人xx丫x性爽| 国产91av在线免费观看| 97超碰精品成人国产| 国产免费一级a男人的天堂| 亚洲av.av天堂| 精品久久久久久久末码| 成熟少妇高潮喷水视频| 大香蕉久久网| 露出奶头的视频| 蜜桃久久精品国产亚洲av| 我的老师免费观看完整版| 少妇裸体淫交视频免费看高清| 午夜久久久久精精品| 久久精品人妻少妇| 亚洲精品国产av成人精品 | 亚洲av一区综合| 一本精品99久久精品77| 3wmmmm亚洲av在线观看| 在线免费十八禁| 亚洲最大成人中文| 丰满乱子伦码专区| 国产在线精品亚洲第一网站| 最好的美女福利视频网| 1024手机看黄色片| 变态另类成人亚洲欧美熟女| 国产黄色小视频在线观看| 观看免费一级毛片| 成人国产麻豆网| 99在线视频只有这里精品首页| 久久精品国产亚洲av涩爱 | 在线免费十八禁| 看免费成人av毛片| 少妇的逼好多水| 校园春色视频在线观看| 两个人的视频大全免费| 日韩高清综合在线| 免费看光身美女| 久久精品影院6| 欧美极品一区二区三区四区| 悠悠久久av| 深爱激情五月婷婷| 色视频www国产| 黄色一级大片看看| 久久人人精品亚洲av| a级毛片免费高清观看在线播放| 亚洲欧美中文字幕日韩二区| 亚洲丝袜综合中文字幕| 99久久久亚洲精品蜜臀av| 成人特级av手机在线观看| 99热这里只有是精品在线观看| 久久久久国产网址| 男女下面进入的视频免费午夜| 搡老妇女老女人老熟妇| 极品教师在线视频| 看免费成人av毛片| 久久亚洲精品不卡| 欧美成人a在线观看| 国产综合懂色| 日韩一区二区视频免费看| 亚洲久久久久久中文字幕| 久久久久国产网址| 日产精品乱码卡一卡2卡三| 国产一区二区三区在线臀色熟女| 麻豆av噜噜一区二区三区| 天堂网av新在线| 精品一区二区三区人妻视频| 日日撸夜夜添| 男人和女人高潮做爰伦理| 亚洲在线观看片| 中文字幕av成人在线电影| 欧美性猛交╳xxx乱大交人| 性欧美人与动物交配| 精品少妇黑人巨大在线播放 | 老女人水多毛片| 午夜福利在线在线| 日韩 亚洲 欧美在线| 女同久久另类99精品国产91| 我要看日韩黄色一级片| 欧美3d第一页| 全区人妻精品视频| 99九九线精品视频在线观看视频| 老师上课跳d突然被开到最大视频| 国产爱豆传媒在线观看| 十八禁网站免费在线| 在现免费观看毛片| 天堂av国产一区二区熟女人妻| 99国产极品粉嫩在线观看| 日本a在线网址| 亚洲精品乱码久久久v下载方式| 黄色一级大片看看| 91麻豆精品激情在线观看国产| 男人和女人高潮做爰伦理| 看片在线看免费视频| 日本五十路高清| 搡老妇女老女人老熟妇| 国产aⅴ精品一区二区三区波| 在线观看av片永久免费下载| 综合色av麻豆| 成人鲁丝片一二三区免费| 日韩国内少妇激情av| 搡女人真爽免费视频火全软件 | 国产精品一区二区免费欧美| ponron亚洲| 日韩欧美免费精品| 成人欧美大片| 国产精品人妻久久久影院| 老司机午夜福利在线观看视频| 国产精品不卡视频一区二区| 欧美3d第一页| 亚洲欧美成人精品一区二区| 欧美激情国产日韩精品一区| 黄片wwwwww| 给我免费播放毛片高清在线观看| 99精品在免费线老司机午夜| 国产精品久久久久久亚洲av鲁大| 亚洲七黄色美女视频| 亚洲在线自拍视频| 日本精品一区二区三区蜜桃| 亚洲第一电影网av| 久久这里只有精品中国| 国产91av在线免费观看| 黄片wwwwww| 国产69精品久久久久777片| 久久久久久九九精品二区国产| 最近手机中文字幕大全| 在线a可以看的网站| 欧美激情在线99| 伦精品一区二区三区| 久久久久久国产a免费观看| 久久久久久久亚洲中文字幕| 亚洲av免费在线观看| 直男gayav资源| 美女免费视频网站| 老司机福利观看| 2021天堂中文幕一二区在线观| 最近2019中文字幕mv第一页| 精品一区二区三区视频在线观看免费| av视频在线观看入口| 久久久精品大字幕| 精品人妻视频免费看| 亚洲精品国产av成人精品 | 午夜福利在线在线| 在现免费观看毛片| 免费在线观看影片大全网站| 国产淫片久久久久久久久| 亚洲国产高清在线一区二区三| 亚洲五月天丁香| 国产毛片a区久久久久| 免费看美女性在线毛片视频| 中国美白少妇内射xxxbb| 成人特级av手机在线观看| 99精品在免费线老司机午夜| 国产一区二区激情短视频| 黄色欧美视频在线观看| 国产亚洲91精品色在线| 久久久久久久亚洲中文字幕| 亚洲熟妇中文字幕五十中出| 成人一区二区视频在线观看| 精品久久久久久久久av| 亚洲av中文字字幕乱码综合| 国产精品日韩av在线免费观看| 免费观看精品视频网站| 男女边吃奶边做爰视频| 午夜免费男女啪啪视频观看 | 天天一区二区日本电影三级| av视频在线观看入口| 日韩国内少妇激情av| 我要搜黄色片| 内射极品少妇av片p| 黄色视频,在线免费观看| 日韩在线高清观看一区二区三区| 亚洲性夜色夜夜综合| 午夜视频国产福利| or卡值多少钱| 中文资源天堂在线| 日本免费a在线| 国产午夜精品论理片| av在线天堂中文字幕| 18+在线观看网站| 国国产精品蜜臀av免费| 午夜日韩欧美国产| 国产精品精品国产色婷婷| 在线播放无遮挡| 露出奶头的视频| 免费av观看视频| 激情 狠狠 欧美| 97超碰精品成人国产| 久久久精品欧美日韩精品| 国产探花在线观看一区二区| 91久久精品国产一区二区三区| 秋霞伦理黄片| 中国国产av一级| 亚洲国产精品999| 久久综合国产亚洲精品| www.av在线官网国产| 国产成人免费无遮挡视频| av有码第一页| 99九九在线精品视频 | 亚洲人与动物交配视频| 亚洲av国产av综合av卡| 国产69精品久久久久777片| 国产午夜精品久久久久久一区二区三区| 丰满饥渴人妻一区二区三| 亚洲美女黄色视频免费看| 亚洲av综合色区一区| 日本av免费视频播放| 嘟嘟电影网在线观看| 3wmmmm亚洲av在线观看| 高清在线视频一区二区三区| 精品午夜福利在线看| 日韩电影二区| 精品一区在线观看国产| 高清视频免费观看一区二区| 91久久精品国产一区二区成人| 色94色欧美一区二区| 亚洲国产最新在线播放| 全区人妻精品视频| 中文乱码字字幕精品一区二区三区| 边亲边吃奶的免费视频| 亚洲欧美日韩另类电影网站| av播播在线观看一区| 如何舔出高潮| 99热网站在线观看| 青春草视频在线免费观看| 精品久久国产蜜桃| 国产成人一区二区在线| av在线播放精品| 亚洲国产精品999| 啦啦啦在线观看免费高清www| 一个人免费看片子| 亚洲欧洲日产国产| 成人综合一区亚洲| 最近最新中文字幕免费大全7| 精品亚洲成国产av| 多毛熟女@视频| 日韩成人伦理影院| 久久国产亚洲av麻豆专区| 久久国产精品大桥未久av | freevideosex欧美| 特大巨黑吊av在线直播| 国产成人免费观看mmmm| 日韩电影二区| 9色porny在线观看| 亚洲成色77777| 婷婷色综合www| 交换朋友夫妻互换小说| 丰满饥渴人妻一区二区三| 日韩中字成人| 丝袜喷水一区| 久久久精品94久久精品| 伊人久久精品亚洲午夜| 精品久久久精品久久久| 伊人亚洲综合成人网| 三级国产精品片| 欧美日韩精品成人综合77777| 丝瓜视频免费看黄片| 一本—道久久a久久精品蜜桃钙片| 成人美女网站在线观看视频| 久久久久人妻精品一区果冻| 看免费成人av毛片| 高清欧美精品videossex| av网站免费在线观看视频| 99热全是精品| 国内揄拍国产精品人妻在线| 观看免费一级毛片| 91久久精品国产一区二区三区| 亚洲在久久综合| 国产免费视频播放在线视频| 亚洲欧美一区二区三区黑人 | 亚洲欧洲国产日韩| tube8黄色片| 免费久久久久久久精品成人欧美视频 | 在线观看国产h片| 免费大片黄手机在线观看| 最后的刺客免费高清国语| 两个人的视频大全免费| 91aial.com中文字幕在线观看| 妹子高潮喷水视频| 日韩制服骚丝袜av| 亚洲丝袜综合中文字幕| 国产一区亚洲一区在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲熟女精品中文字幕| 亚洲欧美一区二区三区国产| av福利片在线观看| 国产日韩欧美视频二区| 寂寞人妻少妇视频99o| a级毛片免费高清观看在线播放| 99热这里只有是精品50| 精品一区二区免费观看| 免费看日本二区| 国产一区二区在线观看av| 中文字幕久久专区| 欧美日韩视频高清一区二区三区二| 国产精品一区二区三区四区免费观看| 日本欧美国产在线视频| 只有这里有精品99| 黑人巨大精品欧美一区二区蜜桃 | 精品国产乱码久久久久久小说| 一个人免费看片子| 国产日韩一区二区三区精品不卡 | 男的添女的下面高潮视频| 欧美丝袜亚洲另类| 99久久中文字幕三级久久日本| 秋霞伦理黄片| 啦啦啦视频在线资源免费观看| 日本-黄色视频高清免费观看| 少妇熟女欧美另类| 精品久久久久久电影网| 亚洲人成网站在线观看播放| 久久韩国三级中文字幕| 中文天堂在线官网| 嘟嘟电影网在线观看| 色吧在线观看| tube8黄色片| 亚洲国产欧美在线一区| 男人添女人高潮全过程视频| 亚洲内射少妇av| 成人美女网站在线观看视频| 精品久久久久久电影网| 成年av动漫网址| 国产欧美日韩综合在线一区二区 | 黄色日韩在线| 久久97久久精品| 久久鲁丝午夜福利片| 大香蕉97超碰在线| 如日韩欧美国产精品一区二区三区 | 国模一区二区三区四区视频| 最近中文字幕2019免费版| 国产免费一区二区三区四区乱码| 黄色日韩在线| 精品久久久精品久久久| 人妻夜夜爽99麻豆av| 三级国产精品片| www.色视频.com| 少妇被粗大的猛进出69影院 | 菩萨蛮人人尽说江南好唐韦庄| av专区在线播放| 久久6这里有精品| 成人毛片60女人毛片免费| 国产精品久久久久久久久免| 国产欧美另类精品又又久久亚洲欧美| 久久久久久久久久成人| 亚洲欧美清纯卡通| 只有这里有精品99| 丝袜喷水一区| 人人妻人人添人人爽欧美一区卜| 制服丝袜香蕉在线| 欧美成人精品欧美一级黄| 亚洲精品乱久久久久久| 人妻一区二区av| 一本大道久久a久久精品| 国产黄频视频在线观看| 国产高清三级在线| 久久人人爽av亚洲精品天堂| 午夜福利在线观看免费完整高清在| 欧美丝袜亚洲另类| xxx大片免费视频| 国产一区二区三区综合在线观看 | 一区二区三区乱码不卡18| 亚洲精品亚洲一区二区| 又爽又黄a免费视频| 精品一区在线观看国产| 成人亚洲欧美一区二区av| 国产日韩欧美视频二区| 中国三级夫妇交换| 亚洲在久久综合| 女性被躁到高潮视频| 亚洲天堂av无毛| 国产淫语在线视频| 中文字幕av电影在线播放| 欧美精品人与动牲交sv欧美| 亚洲图色成人| 成人漫画全彩无遮挡| 免费观看av网站的网址| 极品教师在线视频| 精品久久国产蜜桃| 夫妻性生交免费视频一级片| 国产中年淑女户外野战色| 国产一区二区在线观看av| 蜜桃在线观看..| 春色校园在线视频观看| 日本午夜av视频| av专区在线播放| 好男人视频免费观看在线| 国产极品天堂在线| 国产高清三级在线| 一区二区三区四区激情视频| 久久久亚洲精品成人影院| 日日啪夜夜爽| 人人澡人人妻人| 午夜福利网站1000一区二区三区| 九九爱精品视频在线观看| 亚洲精品国产色婷婷电影| 一级毛片黄色毛片免费观看视频| 国产一级毛片在线| 多毛熟女@视频| 欧美日韩综合久久久久久| 国产69精品久久久久777片| 国产亚洲91精品色在线| 大香蕉久久网| 少妇熟女欧美另类| 国产一区亚洲一区在线观看| 大香蕉久久网| 国产日韩欧美视频二区| 五月玫瑰六月丁香| 日韩熟女老妇一区二区性免费视频| 亚洲精品第二区| 免费不卡的大黄色大毛片视频在线观看| 久久久久久久久久成人| 亚洲欧美中文字幕日韩二区| 日本vs欧美在线观看视频 | 人妻人人澡人人爽人人| 亚洲天堂av无毛| 80岁老熟妇乱子伦牲交| 插阴视频在线观看视频| 哪个播放器可以免费观看大片| av专区在线播放| 性色avwww在线观看| 国产一区二区三区av在线| 亚洲丝袜综合中文字幕| 一级毛片我不卡| 日本av免费视频播放| 久久这里有精品视频免费| 国产免费一区二区三区四区乱码| 国产精品无大码| 80岁老熟妇乱子伦牲交| 高清欧美精品videossex| 又大又黄又爽视频免费| 尾随美女入室| 丝袜脚勾引网站| 国产成人免费观看mmmm| 免费大片黄手机在线观看| 99热全是精品| 国产国拍精品亚洲av在线观看| 色网站视频免费|