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

    進(jìn)入火星大氣的高溫真實(shí)氣體效應(yīng)與氣動(dòng)加熱研究

    2018-10-13 03:05:52楊肖峰國(guó)義軍桂業(yè)偉杜雁霞
    宇航學(xué)報(bào) 2018年9期
    關(guān)鍵詞:激波超聲速氣動(dòng)

    楊肖峰,國(guó)義軍,唐 偉,桂業(yè)偉,杜雁霞

    (中國(guó)空氣動(dòng)力研究與發(fā)展中心空氣動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,綿陽(yáng) 621000)

    0 引 言

    火星探測(cè)是當(dāng)前國(guó)際深空探測(cè)活動(dòng)的熱點(diǎn),得到各航天大國(guó)的高度關(guān)注[1],以美國(guó)為首的航天大國(guó)曾進(jìn)行數(shù)次著陸和巡視任務(wù),并積極實(shí)施樣本返回、載人登火等新的火星探測(cè)計(jì)劃[2-3]。火星任務(wù)在我國(guó)同樣備受關(guān)注[4],《2016年中國(guó)的航天》白皮書(shū)[5]指出,實(shí)施中國(guó)首次火星探測(cè)任務(wù),突破火星環(huán)繞、著陸、巡視探測(cè)等關(guān)鍵技術(shù);2020年發(fā)射首顆火星探測(cè)器,實(shí)施環(huán)繞和巡視聯(lián)合探測(cè)。

    然而,火星探測(cè)器在進(jìn)入火星大氣層的關(guān)鍵階段存在巨大的技術(shù)挑戰(zhàn)。據(jù)統(tǒng)計(jì),截至目前共有16次登陸火星嘗試,僅成功7次,成功率極低[6]。特殊且缺乏足夠數(shù)據(jù)庫(kù)支撐的火星大氣環(huán)境及由此引起的進(jìn)入飛行不確定性是造成探測(cè)器進(jìn)入飛行安全隱患的主要原因之一[7]。返回/進(jìn)入型的航天飛行器穿過(guò)大氣層過(guò)程要經(jīng)歷自由分子流區(qū)、過(guò)渡流區(qū)、滑移流區(qū)和連續(xù)流區(qū),在各流區(qū)的飛行過(guò)程中,氣體分子會(huì)發(fā)生振動(dòng)、離解、電離和復(fù)合等化學(xué)物理反應(yīng)的真實(shí)氣體效應(yīng)?;鹦沁M(jìn)入高超聲速飛行的真實(shí)氣體效應(yīng)[8-9]及其對(duì)氣動(dòng)熱環(huán)境的影響是火星探測(cè)器熱防護(hù)系統(tǒng)設(shè)計(jì)和結(jié)構(gòu)安全可靠性設(shè)計(jì)亟需解決的關(guān)鍵問(wèn)題之一。

    與常規(guī)的空氣飛行環(huán)境相比,火星大氣相對(duì)稀薄且主要含有CO2等氣體成分[6],此類氣體分子構(gòu)型及高溫條件下分子振動(dòng)和離解特性與空氣分子有顯著區(qū)別[10]。為此進(jìn)入火星大氣層的流場(chǎng)需考慮CO2氣體主導(dǎo)的高超聲速化學(xué)反應(yīng)流動(dòng)。考慮火星大氣分子的振動(dòng)激發(fā)和化學(xué)反應(yīng)的時(shí)間尺度,火星進(jìn)入流場(chǎng)的高超聲速熱力學(xué)和化學(xué)非平衡特性及與地球環(huán)境的區(qū)別值得研究[8]?;鹦沁M(jìn)入問(wèn)題的非平衡等真實(shí)氣體效應(yīng)對(duì)探測(cè)器氣動(dòng)加熱的精確預(yù)測(cè)有重要影響[11-12],因此,需要深入開(kāi)展進(jìn)入火星大氣的真實(shí)氣體效應(yīng)分析,獲得氣動(dòng)加熱的主要影響因素,進(jìn)而揭示火星進(jìn)入氣動(dòng)加熱機(jī)理。

    本文針對(duì)火星大氣環(huán)境的分子振動(dòng)和離解行為特征,采用理論分析和數(shù)值模擬兩種手段,開(kāi)展探測(cè)器高超聲速進(jìn)入飛行的高溫真實(shí)氣體效應(yīng)研究與高超聲速氣動(dòng)加熱分析,同時(shí)獲得火星大氣環(huán)境與空氣環(huán)境的相似之處與差異性,以期為火星進(jìn)入數(shù)值模擬的求解器選擇提供技術(shù)參考,同時(shí)為探測(cè)器熱防護(hù)系統(tǒng)設(shè)計(jì)提供數(shù)據(jù)支持。

    1 研究方法

    1.1 理論分析方法

    考察常規(guī)尺度飛行器駐點(diǎn)線上通用氣體介質(zhì)的局部平衡物理化學(xué)過(guò)程,建立高超聲速化學(xué)反應(yīng)流動(dòng)和氣動(dòng)加熱的理論分析方法。激波前后化學(xué)反應(yīng)氣體混合物滿足化學(xué)平衡的Rankine-Hugoniot激波關(guān)系式,即滿足連續(xù)性方程、動(dòng)量方程和能量方程等三大守恒方程

    ρ∞V∞=ρ2V2

    (1)

    (2)

    (3)

    式中:Δh0為生成焓,h為顯焓,根據(jù)氣體模型由溫度求得,下標(biāo)∞為自由來(lái)流參數(shù),下標(biāo)2為激波后參數(shù)。激波后參數(shù)由上述守恒方程通過(guò)迭代求解來(lái)獲得,其中激波前后氣體混合物各組分的質(zhì)量分?jǐn)?shù)由化學(xué)反應(yīng)平衡常數(shù)和當(dāng)?shù)販囟韧ㄟ^(guò)內(nèi)部迭代求解獲得。已知激波后的參數(shù),激波層內(nèi)從激波后到邊界層外緣的滯止過(guò)程為等熵流動(dòng)且滿足

    (4)

    通過(guò)離散求解獲得激波層內(nèi)滯止過(guò)程的顯焓。因比熱為溫度的函數(shù),等熵關(guān)系式無(wú)直觀表達(dá)式,需要離散求解。

    針對(duì)高焓氣體高超聲速化學(xué)反應(yīng)流動(dòng)問(wèn)題,定義化學(xué)反應(yīng)和振動(dòng)激發(fā)特征時(shí)間尺度與流場(chǎng)特征時(shí)間尺度之比為Damokhler數(shù)[13-14](Da)和振動(dòng)數(shù)(Vi)

    (5)

    (6)

    分別用以表征高超聲速流動(dòng)的化學(xué)和熱力學(xué)非平衡程度。當(dāng)激波層內(nèi)化學(xué)反應(yīng)很快時(shí),特征時(shí)間很短,Da→ ∞,流動(dòng)處于化學(xué)平衡狀態(tài);當(dāng)Da為較小的有限值,流動(dòng)處于化學(xué)非平衡狀態(tài);當(dāng)化學(xué)反應(yīng)很慢時(shí),Da→ 0,流動(dòng)處于化學(xué)凍結(jié)狀態(tài)。Vi數(shù)同理。高超聲速流動(dòng)、化學(xué)反應(yīng)和分子振動(dòng)激發(fā)的特征時(shí)間由激波后平衡流動(dòng)參數(shù)來(lái)確定[8, 15-16]

    (7)

    (8)

    (9)

    式中:火星大氣和空氣介質(zhì)的化學(xué)反應(yīng)速率kf由化學(xué)反應(yīng)動(dòng)力學(xué)模型[15, 17]求得,振動(dòng)特征尺度相關(guān)系數(shù)C1和C2由實(shí)驗(yàn)測(cè)得[15, 18-19]。

    對(duì)火星大氣和空氣介質(zhì)的高超聲速鈍頭體繞流氣動(dòng)加熱問(wèn)題,駐點(diǎn)熱流均可推導(dǎo)出其自相似解,即Fay-Riddell公式[18]

    (10)

    (11)

    (12)

    以上理論分析方法是針對(duì)通用氣體化學(xué)反應(yīng)流動(dòng)推導(dǎo)出來(lái)的,對(duì)常規(guī)尺度飛行器的地球再入和火星進(jìn)入問(wèn)題具有普適性,但涉及的模型參數(shù)存在差異。上述方法同樣適用于完全氣體情況。

    1.2 數(shù)值模擬方法

    數(shù)值模擬控制方程為考慮化學(xué)反應(yīng)的三維層流可壓縮Navier-Stokes方程組。對(duì)完全氣體而言,控制方程僅包括連續(xù)性方程、動(dòng)量和能量方程,而對(duì)化學(xué)反應(yīng)氣體,控制方程還包括Ns-Ne個(gè)組分連續(xù)性方程(Ns和Ne為組分和元素個(gè)數(shù)),界面或壁面上熱流通量除含熱傳導(dǎo)項(xiàng)外,還包括因組分?jǐn)U散所致的化學(xué)熱通量項(xiàng)

    (13)

    使用TVD型的有限體積法對(duì)控制方程作數(shù)值離散??臻g無(wú)黏通量使用二階Van-Leer方法作矢通量分裂,黏性通量采用中心格式作離散,界面變量通量采用帶有van Albada限制器的MUSCL方法插值求得。時(shí)間推進(jìn)上采用LUSGS隱式方法。由于所計(jì)算的火星大氣和空氣介質(zhì)條件下飛行器進(jìn)入雷諾數(shù)較低,數(shù)值計(jì)算僅考慮層流流動(dòng)。求解完全氣體和化學(xué)反應(yīng)氣體的數(shù)值程序已通過(guò)大量火星大氣和空氣環(huán)境的算例驗(yàn)證[6,20-21]。

    1.3 氣體模型

    針對(duì)火星大氣環(huán)境,僅考慮其主要成分CO2氣體;而針對(duì)空氣環(huán)境,考慮由質(zhì)量分?jǐn)?shù)分別為77%和23%的N2和O2混合物。CO2為直線型三原子分子,具有3個(gè)平動(dòng)和2個(gè)轉(zhuǎn)動(dòng)自由度,在高溫條件下振動(dòng)能會(huì)部分或完全激發(fā),其三個(gè)振動(dòng)模態(tài)分別為彎曲、對(duì)稱拉伸和反對(duì)稱拉伸模態(tài)。N2和O2均為雙原子分子,同樣具有3個(gè)平動(dòng)和2個(gè)轉(zhuǎn)動(dòng)自由度,在高溫條件下僅具有1個(gè)拉伸模態(tài)。高溫火星大氣和空氣及其離解物的物理化學(xué)數(shù)據(jù)見(jiàn)表1。

    表1 高溫火星大氣和空氣物理化學(xué)數(shù)據(jù)Table 1 Physical and chemical data for the high-temperature Martian gas and air

    根據(jù)統(tǒng)計(jì)熱力學(xué)理論,氣體的熱力學(xué)特性由氣體分子的平動(dòng)能、轉(zhuǎn)動(dòng)能、振動(dòng)能等決定,氣體混合物各組分的內(nèi)能關(guān)系式[18]為

    (14)

    量熱完全氣體的熱容和比熱比是恒定的,考慮氣體分子振動(dòng)未激發(fā),火星大氣和空氣的比熱比均為1.4,高溫條件氣體分子振動(dòng)完全激發(fā),比熱比分別降至1.154和1.286。熱完全氣體或化學(xué)反應(yīng)氣體混合物的熱容隨溫度而改變,比熱比不再恒定。基于地球和火星大氣分子的高溫物理化學(xué)行為,建立火星大氣和空氣的量熱完全、熱完全和多種化學(xué)反應(yīng)氣體模型,詳見(jiàn)表2和表3。

    在數(shù)值模擬方面,采用基于Park化學(xué)動(dòng)力學(xué)的火星大氣5組分6化學(xué)反應(yīng)(M-5S6R)模型[15, 17]和基于Dunn-Kang化學(xué)動(dòng)力學(xué)的空氣5組分5化學(xué)反應(yīng)(A-5S5R)模型[22-23],涉及的反應(yīng)速率系數(shù)見(jiàn)文獻(xiàn)[17]。在理論分析方面,考慮描述火星大氣環(huán)境的3組分1化學(xué)反應(yīng)(M-3S1R)、4組分2化學(xué)反應(yīng)(M- 4S2R)和5組分3化學(xué)反應(yīng)(M-5S3R)模型,并考慮描述空氣環(huán)境的5組分3化學(xué)反應(yīng)(A-5S3R)模型。

    表2 高溫火星大氣和空氣流動(dòng)涉及的化學(xué)反應(yīng)Table 2 Chemical reactions in the Martian gas and the air for high-temperature flows

    表3 火星大氣和空氣氣體模型Table 3 Models of the Martian gas and air

    關(guān)于火星大氣和空氣介質(zhì)及其離解產(chǎn)物的輸運(yùn)特性,黏性系數(shù)由Blottner擬合關(guān)系式獲得,熱傳導(dǎo)系數(shù)由Prandtl數(shù)求得,取Pr=0.71,混合氣體參數(shù)由Wilke公式計(jì)算。組分?jǐn)U散系數(shù)由Schmitt數(shù)獲得,取Sc=0.5。針對(duì)理論分析,假設(shè)黏性系數(shù)為恒定值。

    2 火星進(jìn)入流場(chǎng)的真實(shí)氣體效應(yīng)

    2.1 激波層氣體分子的振動(dòng)激發(fā)與化學(xué)反應(yīng)

    探測(cè)器以極高速度進(jìn)入以CO2為主要成分的火星大氣層并產(chǎn)生高超聲速流動(dòng)。隨著大氣層由稀薄到稠密,探測(cè)器歷經(jīng)自由分子流、過(guò)渡流、滑移流、連續(xù)介質(zhì)流等不同流區(qū)。圖1給出火星探測(cè)器標(biāo)稱進(jìn)入軌道及等Knudsen數(shù)隨飛行速度和飛行高度的變化曲線。按連續(xù)介質(zhì)流區(qū)Kn<0.03的定義,火星大氣在40 km以下才進(jìn)入了連續(xù)介質(zhì)流區(qū),流區(qū)高度界限低于地球大氣情況(60 km),高空稀薄效應(yīng)明顯??紤]到火星大氣比較稀薄、氣溫較低,火星探測(cè)器流動(dòng)呈現(xiàn)高馬赫數(shù)、低雷諾數(shù)、流動(dòng)稀薄等特點(diǎn)。

    圖1 Kn和離解率隨進(jìn)入速度和高度的變化曲線Fig.1 Kn and dissociation percentage curves

    進(jìn)入火星大氣層的高超聲速流動(dòng)與空氣流動(dòng)相似,在極高速度條件下激波層內(nèi)伴有雙/多原子分子的離解和電離等化學(xué)反應(yīng)過(guò)程。激波層內(nèi)高焓火星大氣以CO2離解為主,同時(shí)存在離解物O2和CO的離解反應(yīng)過(guò)程。圖2為M-5S3R機(jī)理作用下的沿進(jìn)入軌道各組分質(zhì)量分?jǐn)?shù)分布曲線。從圖2可以看出,CO2在整個(gè)高超聲速進(jìn)入過(guò)程均發(fā)生離解,尤其在進(jìn)入速度4 km/s以上的早期進(jìn)入段幾乎完全離解。混合物中O2在4 km/s以上的進(jìn)入段發(fā)生較大規(guī)模的離解,在6 km/s以上的早期進(jìn)入段幾乎完全離解?;旌衔镏蠧O只在6 km/s以上的早期進(jìn)入段才發(fā)生離解。在整個(gè)進(jìn)入段,CO成分因CO2的離解而占較大比重;而在進(jìn)入早期O質(zhì)量分?jǐn)?shù)占優(yōu),在進(jìn)入后期O2質(zhì)量分?jǐn)?shù)占優(yōu)。

    圖2 沿進(jìn)入軌道的波后質(zhì)量分?jǐn)?shù)Fig.2 Post-shock mass fraction along the trajectory

    圖3 典型反應(yīng)的極限離解率隨進(jìn)入速度的變化曲線Fig.3 Dissociation limitation for typical dissociation reactions along the trajectory

    圖3為典型化學(xué)反應(yīng)過(guò)程的極限離解率隨進(jìn)入速度的變化曲線,因單位質(zhì)量CO2離解焓低于O2和N2,故CO2離解為CO和O2的極限離解率高于其它離解過(guò)程。

    多原子分子振動(dòng)能的影響比雙原子/單原子分子復(fù)雜。盡管常溫下CO2的平動(dòng)和轉(zhuǎn)動(dòng)自由度與N2,O2相同,但高溫條件下振動(dòng)能激發(fā)模式有差異。由于高溫CO2振動(dòng)模態(tài)較多且振動(dòng)特征溫度從中低溫度到高溫均有散布,故火星大氣分子振動(dòng)能激發(fā)特性明顯強(qiáng)于空氣。

    圖4為使用不同的氣體模型獲得的沿進(jìn)入軌道的激波后溫度變化曲線。使用量熱完全氣體(CPG)模型獲得的激波層內(nèi)溫度最高,且在進(jìn)入早期出現(xiàn)數(shù)萬(wàn)度的“虛假”溫度;考慮高溫條件下CO2分子振動(dòng)激發(fā),使用熱完全氣體(TPG)模型獲得的激波層內(nèi)溫度相對(duì)較低。使用化學(xué)反應(yīng)模型獲得的激波后溫度最低,且不同的化學(xué)反應(yīng)模型在高速段的激波后溫度有差別,說(shuō)明除CO2發(fā)生離解外,在4.2 km/s以上O2發(fā)生離解,6.1 km/s以上CO發(fā)生離解。對(duì)比各個(gè)氣體模型可知,沿進(jìn)入軌道的高超聲速流動(dòng)必須考慮激波層內(nèi)的分子振動(dòng)激發(fā)和化學(xué)反應(yīng)過(guò)程,尤其在高空高速流動(dòng)區(qū)域需要考慮更多組分的化學(xué)反應(yīng)過(guò)程。

    圖4 不同的氣體模型獲得的沿進(jìn)入軌道的波后溫度Fig.4 Post-shock temperature with different gas models along the trajectory

    運(yùn)用數(shù)值模擬手段,針對(duì)直徑76 mm的球頭繞流(V∞= 2772 m/s,p∞= 2132 Pa,T∞= 712 K),研究火星大氣與空氣環(huán)境下的氣體模型對(duì)高超聲速流場(chǎng)的影響。圖5給出了火星大氣與空氣環(huán)境下不同氣體模型獲得的脫體激波形狀和比熱比。從圖5可以看出,比熱比降低(熱容增大),激波脫體距離減小,熱完全氣體模型獲得的激波脫體距離與分子振動(dòng)完全激發(fā)條件下的量熱完全氣體模型結(jié)果相近,同時(shí)也與化學(xué)反應(yīng)氣體模型獲得的結(jié)果相近?;鹦谴髿夥肿诱駝?dòng)激發(fā)強(qiáng)于空氣,故火星大氣流場(chǎng)比熱比較低,激波脫體距離較薄。

    圖5 火星大氣與空氣環(huán)境的脫體激波形狀與比熱比Fig.5 Shock shapes and specific heat ratios for the Martian atmosphere and the air

    2.2 化學(xué)和熱力學(xué)非平衡特性

    火星探測(cè)器高超聲速進(jìn)入大氣層的過(guò)程產(chǎn)生的弓形激波使波后氣體溫度升高,促使激波后氣體發(fā)生分子振動(dòng)、離解等物理化學(xué)反應(yīng),分子振動(dòng)和離解存在一定的弛豫過(guò)程,其特征時(shí)間尺度隨進(jìn)入速度和高度而改變,激波層內(nèi)氣體呈現(xiàn)化學(xué)/熱力學(xué)非平衡或平衡狀態(tài)。

    圖6給出了駐點(diǎn)附近高超聲速火星大氣流動(dòng)在飛行速度-高度域內(nèi)的流動(dòng)狀態(tài)。從圖6可以看出,火星探測(cè)器僅在在高速低空域內(nèi)處于化學(xué)平衡狀態(tài),而在較低速度或較高高度域內(nèi)飛行均處于化學(xué)非平衡狀態(tài);火星探測(cè)器在高速低空大部分域內(nèi)均處于熱力學(xué)平衡狀態(tài),而僅在低速高空域內(nèi)飛行處于熱力學(xué)非平衡狀態(tài)。

    圖7進(jìn)一步對(duì)比了火星大氣介質(zhì)和空氣介質(zhì)的熱力學(xué)與化學(xué)非平衡特性。從圖7可以看出,保持兩種介質(zhì)的大氣壓力和溫度條件相同(密度因介質(zhì)差異而不同),火星大氣介質(zhì)較空氣介質(zhì)具有更強(qiáng)的化學(xué)非平衡特性,但具有相對(duì)較弱的熱力學(xué)非平衡特性。圖7還給出了火星探路者號(hào)探測(cè)器的標(biāo)稱進(jìn)入軌道。從圖7可以看出,在進(jìn)入初期,探測(cè)器處于化學(xué)和熱力學(xué)非平衡狀態(tài);在進(jìn)入大部分時(shí)段,探測(cè)器保持化學(xué)非平衡狀態(tài),分子振動(dòng)激發(fā)逐漸趨于平衡狀態(tài)。

    圖7 駐點(diǎn)附近高超聲速火星氣體和空氣流動(dòng)的熱力學(xué)與化學(xué)非平衡特性對(duì)比Fig.7 Comparison of thermal and chemical non-equilibrium characteristics between hypersonic Martian gas and air flow in the stagnation region

    3 火星探測(cè)器高超聲速氣動(dòng)加熱

    進(jìn)入火星大氣層的真實(shí)氣體效應(yīng)影響火星探測(cè)器物面的高超聲速氣動(dòng)加熱量。針對(duì)火星探路者號(hào)進(jìn)入器70°球錐型防熱大底及其標(biāo)稱彈道[6],開(kāi)展理論計(jì)算和數(shù)值模擬分析。

    考慮完全氣體和化學(xué)反應(yīng)氣體模型,圖8為基于火星大氣介質(zhì)的沿火星進(jìn)入軌道的駐點(diǎn)氣動(dòng)加熱熱流變化曲線。從圖8可以看出,沿進(jìn)入軌道的氣動(dòng)加熱量先升高再降低,在兼顧較高的進(jìn)入速度和較大的大氣密度的6~7 km/s段出現(xiàn)氣動(dòng)加熱峰值。對(duì)比各個(gè)氣體模型可知,使用化學(xué)反應(yīng)模型獲得的駐點(diǎn)熱流因表面催化特性的不同而產(chǎn)生較大的差異,完全催化(FCW)和完全非催化(NCW)壁熱流相差3~4倍,完全氣體模型的結(jié)果介于完全催化和完全非催化壁結(jié)果的中間,其中熱完全氣體(TPG)模型結(jié)果略高于量熱完全氣體(CPG)模型。

    圖8 駐點(diǎn)熱流(火星大氣介質(zhì)與火星進(jìn)入軌道)Fig.8 Trajectory based stagnation heatfluxes (Martian gas and Mars entry trajectory)

    CO2的離解作用造成激波層內(nèi)溫升低于完全氣體情況,進(jìn)而降低溫度梯度所致的熱傳導(dǎo)熱流量,并造成完全非催化壁條件下的氣動(dòng)加熱量低于完全氣體情況。完全催化壁條件強(qiáng)制壁面混合氣體化學(xué)焓最低,促使能量全部通過(guò)壁面催化反應(yīng)以化學(xué)生成焓的形式在壁面釋放,反應(yīng)所致組分?jǐn)U散熱流成為氣動(dòng)加熱的重要組成部分。與量熱完全氣體模型相比,熱完全氣體模型使激波層減薄,駐點(diǎn)邊界層外緣密度升高,進(jìn)而造成壁面熱流高于量熱完全氣體情況。在中后期大部分進(jìn)入段,使用M-5S3R和M- 4S 2R化學(xué)反應(yīng)模型獲得的氣動(dòng)加熱量相同,僅在6 km以上的早期進(jìn)入段存在差別,說(shuō)明在此期間CO才發(fā)生離解反應(yīng)。

    保持沿進(jìn)入軌道的大氣壓力和溫度不變,將氣體介質(zhì)置換為空氣介質(zhì),圖9為基于空氣介質(zhì)的沿火星進(jìn)入軌道的駐點(diǎn)氣動(dòng)加熱熱流變化曲線。由圖9可知,在相同來(lái)流條件下,使用不同的氣體模型獲得的火星大氣環(huán)境的物面熱流值均高于空氣環(huán)境。由于空氣介質(zhì)的分子振動(dòng)激發(fā)特性弱于CO2氣體,量熱與熱完全氣體模型獲得的駐點(diǎn)熱流差值低于火星大氣介質(zhì)情況。在進(jìn)入段后期,探測(cè)器飛行高度和速度均較低,完全催化與完全非催化壁結(jié)果的差值低于火星大氣介質(zhì)情況,說(shuō)明在該階段的N2和O2分子離解程度較CO2分子弱??傊?,CO2為主的火星大氣環(huán)境下的高超聲速氣動(dòng)加熱強(qiáng)于N2和O2為主的地球大氣。

    圖9 駐點(diǎn)熱流(空氣介質(zhì)與火星進(jìn)入軌道)Fig.9 Trajectory based stagnation heatfluxes (the air and Mars entry trajectory)

    盡管如此,但火星大氣密度較低,真實(shí)進(jìn)入過(guò)程的氣動(dòng)熱載荷未必強(qiáng)于地球再入熱環(huán)境。在地球大氣環(huán)境設(shè)計(jì)相同的進(jìn)入軌道,即進(jìn)入速度和高度保持不變,相同高度下的大氣壓力和密度高于火星環(huán)境。圖10為基于空氣介質(zhì)的沿地球進(jìn)入軌道的駐點(diǎn)氣動(dòng)加熱熱流變化曲線。由圖10可知,使用各種氣體模型獲得的地球再入熱環(huán)境均強(qiáng)于火星進(jìn)入熱環(huán)境。盡管CO2氣體的參與加劇了氣動(dòng)加熱量,但稀薄的火星大氣造成真實(shí)進(jìn)入過(guò)程的氣動(dòng)熱載荷弱于地球再入過(guò)程。

    圖10 駐點(diǎn)熱流(空氣介質(zhì)與地球再入軌道)Fig.10 Trajectory based stagnation heatfluxes (the air and Earth reentry trajectory)

    圖11給出了高度為29.4 km、速度為4241.6 m/s的彈道點(diǎn)上(處于連續(xù)流區(qū)內(nèi))在300 K壁溫條件下的氣動(dòng)加熱數(shù)值模擬結(jié)果,氣體模型包括量熱完全(比熱比1.29和1.15)、熱完全(TPG)和M-5S6R化學(xué)反應(yīng)(完全非催化NCW和完全催化FCW)模型。從圖11可以看出,完全氣體模型條件下分子振動(dòng)激發(fā)和化學(xué)反應(yīng)氣體模型的壁面催化復(fù)合作用均能增加氣動(dòng)加熱量,且完全氣體模型結(jié)果處于化學(xué)反應(yīng)模型的完全非催化與完全催化壁結(jié)果之間。由于理論和數(shù)值計(jì)算參數(shù)差別甚大,兩種手段計(jì)算結(jié)果有一定的差別,但數(shù)值模擬獲得的氣動(dòng)加熱規(guī)律與理論分析結(jié)果相一致。

    圖11 不同氣體模型下壁面熱流CFD結(jié)果對(duì)比Fig.11 Comparison of CFD results of surface heatfluxes with different gas models

    從能量守恒和分配的角度分析,來(lái)流氣體經(jīng)過(guò)激波后分子振動(dòng)激發(fā),氣體熱容增加,相同溫升條件下單位質(zhì)量氣體承載的能量增加,氣體分子儲(chǔ)能能力的提升會(huì)造成較薄的脫體激波層以抵抗高超聲速來(lái)流的動(dòng)能沖擊,激波后溫度顯著降低且激波層內(nèi)氣體被壓縮,密度升高,平衡狀態(tài)下的較高密度的儲(chǔ)能氣體在壁面上需要釋放能量以回歸壁面較低溫度的熱狀態(tài),勢(shì)必產(chǎn)生較大的壁面熱流。因此,熱完全氣體模型的熱流結(jié)果介于分子振動(dòng)完全未激發(fā)與完全激發(fā)的量熱完全氣體模型之間,其與振動(dòng)完全未激發(fā)熱流結(jié)果的差值反映分子振動(dòng)激發(fā)儲(chǔ)能能力。

    激波層內(nèi)的化學(xué)反應(yīng)是另外一種儲(chǔ)能方式,通過(guò)化學(xué)生成熱將高速來(lái)流的動(dòng)能以化學(xué)能的形式儲(chǔ)存于離解混合物中,該部分能量需要離解混合物的物面催化復(fù)合行為來(lái)予以釋放。完全催化壁面使得離解混合氣體的化學(xué)生成焓在壁面完全釋放,故氣動(dòng)加熱量最大;完全非催化壁面使得離解混合氣體的化學(xué)生成焓在壁面完全不釋放,而隨繞流離開(kāi)壁面,故氣動(dòng)加熱量最??;實(shí)際情況應(yīng)為有限速率的部分催化壁,熱流理應(yīng)介于二者之間。

    4 結(jié) 論

    基于高焓大氣介質(zhì)的物理化學(xué)特性,本文理論分析了飛行器進(jìn)入火星大氣所遇到的真實(shí)氣體效應(yīng),并獲得了火星進(jìn)入飛行器高超聲速氣動(dòng)加熱的理論和數(shù)值預(yù)測(cè)結(jié)果。分析認(rèn)為:探測(cè)器進(jìn)入火星大氣層的稀薄效應(yīng)明顯;因CO2分子振動(dòng)激發(fā)模態(tài)多且特征溫度較低,再者發(fā)生離解反應(yīng)的溫度閾值較低,火星進(jìn)入流動(dòng)為多種振動(dòng)模態(tài)激發(fā)的高超聲速高焓化學(xué)反應(yīng)流動(dòng);激波層內(nèi)發(fā)生CO2氣體為主的大規(guī)模離解,隨著激波層內(nèi)溫度的升高,O2和CO也將發(fā)生離解反應(yīng);探測(cè)器在大部分進(jìn)入軌道時(shí)段處于化學(xué)非平衡但熱力學(xué)平衡狀態(tài)。

    采用多種氣體模型對(duì)比分析來(lái)流氣體介質(zhì)和真實(shí)氣體效應(yīng)對(duì)氣動(dòng)加熱的影響,結(jié)果表明:激波層能量?jī)?chǔ)存和分配模式因氣體分子振動(dòng)激發(fā)和化學(xué)反應(yīng)而改變,氣體較強(qiáng)的儲(chǔ)能能力可有效降低激波層內(nèi)的溫度,但高溫儲(chǔ)能介質(zhì)在物面釋放能量會(huì)增大氣動(dòng)加熱量;分子振動(dòng)激發(fā)會(huì)增強(qiáng)氣動(dòng)加熱量,但均介于化學(xué)反應(yīng)模型的完全非催化和完全催化壁結(jié)果之間;CO2為主的火星大氣環(huán)境下的高超聲速氣動(dòng)加熱強(qiáng)于N2和O2為主的地球大氣,但稀薄的火星大氣造成真實(shí)進(jìn)入過(guò)程的氣動(dòng)熱載荷弱于地球再入過(guò)程。最后,表4系統(tǒng)總結(jié)了火星大氣和空氣介質(zhì)相關(guān)的高溫真實(shí)氣體效應(yīng)的差異性。

    表4 火星大氣和空氣相關(guān)高溫真實(shí)氣體效應(yīng)的差異性Table 4 Difference of high-temperature real-gas effects of the Martian gas and air

    猜你喜歡
    激波超聲速氣動(dòng)
    中寰氣動(dòng)執(zhí)行機(jī)構(gòu)
    高超聲速出版工程
    高超聲速飛行器
    基于NACA0030的波紋狀翼型氣動(dòng)特性探索
    一種基于聚類分析的二維激波模式識(shí)別算法
    基于HIFiRE-2超燃發(fā)動(dòng)機(jī)內(nèi)流道的激波邊界層干擾分析
    基于反饋線性化的RLV氣動(dòng)控制一體化設(shè)計(jì)
    斜激波入射V形鈍前緣溢流口激波干擾研究
    超聲速旅行
    適于可壓縮多尺度流動(dòng)的緊致型激波捕捉格式
    亚洲无线观看免费| 天美传媒精品一区二区| 国产午夜福利久久久久久| 乱码一卡2卡4卡精品| 国产人妻一区二区三区在| av在线天堂中文字幕| 午夜福利网站1000一区二区三区| 国产又黄又爽又无遮挡在线| 亚洲欧美中文字幕日韩二区| 听说在线观看完整版免费高清| 久久99热这里只有精品18| 22中文网久久字幕| 亚洲欧美一区二区三区国产| 成年女人永久免费观看视频| 亚洲欧美成人综合另类久久久 | 一区二区三区高清视频在线| 不卡视频在线观看欧美| 91精品国产九色| 午夜福利在线观看吧| 久久久久久久久久黄片| 午夜福利高清视频| 尾随美女入室| 蜜桃亚洲精品一区二区三区| 国产真实乱freesex| 免费观看性生交大片5| 九草在线视频观看| 国产精品熟女久久久久浪| 色噜噜av男人的天堂激情| 亚洲成人av在线免费| 大又大粗又爽又黄少妇毛片口| 成人性生交大片免费视频hd| 久久综合国产亚洲精品| 午夜a级毛片| 久久精品91蜜桃| kizo精华| 亚洲精品一区蜜桃| 午夜激情欧美在线| 国产精品综合久久久久久久免费| 免费观看精品视频网站| 国产69精品久久久久777片| 18禁在线无遮挡免费观看视频| 精品久久国产蜜桃| 亚洲国产精品sss在线观看| h日本视频在线播放| 国产淫语在线视频| 老女人水多毛片| 国产淫片久久久久久久久| 国产一区二区亚洲精品在线观看| 久久久成人免费电影| 久久久久久九九精品二区国产| 色网站视频免费| 两性午夜刺激爽爽歪歪视频在线观看| 老司机福利观看| 午夜精品在线福利| 青青草视频在线视频观看| 国产亚洲91精品色在线| 91av网一区二区| 22中文网久久字幕| 国内精品美女久久久久久| 人妻制服诱惑在线中文字幕| 啦啦啦观看免费观看视频高清| 建设人人有责人人尽责人人享有的 | 亚洲国产精品成人久久小说| 国产精品一区二区三区四区久久| 国产亚洲5aaaaa淫片| 精品不卡国产一区二区三区| 日韩欧美在线乱码| 中文字幕亚洲精品专区| 久久久久网色| 毛片一级片免费看久久久久| 色哟哟·www| 亚洲怡红院男人天堂| 午夜日本视频在线| 小说图片视频综合网站| 91狼人影院| 午夜精品一区二区三区免费看| 国产精品久久久久久久久免| 中文乱码字字幕精品一区二区三区 | 欧美激情久久久久久爽电影| 日韩视频在线欧美| 亚洲精品乱码久久久v下载方式| 日韩 亚洲 欧美在线| 丰满少妇做爰视频| 最近视频中文字幕2019在线8| 少妇丰满av| 边亲边吃奶的免费视频| 黄色日韩在线| 国产黄片美女视频| 欧美高清成人免费视频www| 国产乱人偷精品视频| 久热久热在线精品观看| av女优亚洲男人天堂| 国产亚洲av嫩草精品影院| 一级毛片我不卡| videossex国产| 国产精品人妻久久久影院| 日韩精品有码人妻一区| 精华霜和精华液先用哪个| 神马国产精品三级电影在线观看| 边亲边吃奶的免费视频| 日本与韩国留学比较| 91aial.com中文字幕在线观看| 亚洲av日韩在线播放| 三级男女做爰猛烈吃奶摸视频| 99久久九九国产精品国产免费| 欧美日韩国产亚洲二区| 久久久精品大字幕| 久久久a久久爽久久v久久| 国产在线一区二区三区精 | 国产一区有黄有色的免费视频 | 亚洲一级一片aⅴ在线观看| 亚洲精品日韩av片在线观看| 中文天堂在线官网| 日本免费在线观看一区| 色综合色国产| 久久鲁丝午夜福利片| 在线播放无遮挡| 亚洲av成人精品一区久久| 精品人妻一区二区三区麻豆| 十八禁国产超污无遮挡网站| av线在线观看网站| 久久精品久久久久久久性| 97超碰精品成人国产| 成年av动漫网址| 中文字幕免费在线视频6| 国产精品久久电影中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 最近手机中文字幕大全| 性插视频无遮挡在线免费观看| 日本av手机在线免费观看| 日本免费a在线| 欧美一区二区亚洲| 91久久精品国产一区二区成人| 美女内射精品一级片tv| 亚洲av熟女| 亚洲av一区综合| av视频在线观看入口| 欧美激情在线99| 性色avwww在线观看| 久久久色成人| 国产精品国产三级国产专区5o | 久久久久久久久中文| 99热精品在线国产| 国产免费又黄又爽又色| 国产在视频线在精品| 蜜桃久久精品国产亚洲av| 久热久热在线精品观看| 免费看av在线观看网站| 青春草国产在线视频| 少妇熟女欧美另类| 中文字幕精品亚洲无线码一区| 精品人妻熟女av久视频| 国产真实乱freesex| 久久久久网色| 国产老妇伦熟女老妇高清| 国产亚洲精品久久久com| 男插女下体视频免费在线播放| 久久午夜福利片| 欧美成人a在线观看| 丝袜美腿在线中文| 欧美激情在线99| 免费搜索国产男女视频| av在线播放精品| 一区二区三区四区激情视频| 一边摸一边抽搐一进一小说| 亚洲五月天丁香| 日本欧美国产在线视频| 中文乱码字字幕精品一区二区三区 | 成人毛片a级毛片在线播放| 精品不卡国产一区二区三区| 国产成人精品久久久久久| 欧美三级亚洲精品| av在线观看视频网站免费| 久久这里只有精品中国| 日韩欧美三级三区| 我的女老师完整版在线观看| 国产私拍福利视频在线观看| 亚洲自偷自拍三级| 青青草视频在线视频观看| 一个人观看的视频www高清免费观看| 久久午夜福利片| 成人高潮视频无遮挡免费网站| 亚洲国产精品国产精品| 国产又黄又爽又无遮挡在线| 免费一级毛片在线播放高清视频| 天天一区二区日本电影三级| 久久精品91蜜桃| 国产成人精品久久久久久| 两个人的视频大全免费| 91久久精品电影网| 久久久亚洲精品成人影院| 2022亚洲国产成人精品| 哪个播放器可以免费观看大片| 床上黄色一级片| 日韩欧美精品免费久久| 99在线视频只有这里精品首页| 人妻少妇偷人精品九色| 久久久久久久国产电影| 国产综合懂色| 国产激情偷乱视频一区二区| 在线播放无遮挡| 一个人看视频在线观看www免费| 精品少妇黑人巨大在线播放 | 亚洲怡红院男人天堂| 黑人高潮一二区| 免费大片18禁| av国产久精品久网站免费入址| 天天躁日日操中文字幕| 欧美色视频一区免费| 国产精品美女特级片免费视频播放器| 一卡2卡三卡四卡精品乱码亚洲| 真实男女啪啪啪动态图| 在线播放无遮挡| 在线a可以看的网站| 国产精品福利在线免费观看| 日韩制服骚丝袜av| 少妇的逼好多水| 69av精品久久久久久| 97热精品久久久久久| 久99久视频精品免费| 嘟嘟电影网在线观看| 一级毛片我不卡| 青春草国产在线视频| 国产精品久久电影中文字幕| 尤物成人国产欧美一区二区三区| 久久久a久久爽久久v久久| 亚洲av中文字字幕乱码综合| 欧美日韩国产亚洲二区| 亚洲成色77777| 国产女主播在线喷水免费视频网站 | 国产单亲对白刺激| 亚洲精品一区蜜桃| 国产欧美日韩精品一区二区| 国产成人一区二区在线| 亚洲经典国产精华液单| 国产免费福利视频在线观看| 只有这里有精品99| 婷婷六月久久综合丁香| 国产国拍精品亚洲av在线观看| 欧美性猛交╳xxx乱大交人| 一个人观看的视频www高清免费观看| 大香蕉97超碰在线| 韩国av在线不卡| 久久久久久久久久久丰满| 欧美潮喷喷水| 国语自产精品视频在线第100页| 只有这里有精品99| 亚洲精品乱久久久久久| 丝袜美腿在线中文| 国产精品一区www在线观看| 久久久久久久久久久丰满| 91精品一卡2卡3卡4卡| 久久久久久久久久久免费av| 两个人的视频大全免费| 久久亚洲精品不卡| 免费观看人在逋| 精品欧美国产一区二区三| 99久久成人亚洲精品观看| 免费av观看视频| 欧美三级亚洲精品| av播播在线观看一区| 国产一区亚洲一区在线观看| 精品免费久久久久久久清纯| 男人和女人高潮做爰伦理| 国产精品伦人一区二区| 又爽又黄a免费视频| 搞女人的毛片| 亚洲精品久久久久久婷婷小说 | 成年女人永久免费观看视频| 亚洲精品影视一区二区三区av| 99在线视频只有这里精品首页| 午夜精品国产一区二区电影 | 国内少妇人妻偷人精品xxx网站| 亚洲怡红院男人天堂| 国国产精品蜜臀av免费| 99国产精品一区二区蜜桃av| 亚洲欧美中文字幕日韩二区| 欧美另类亚洲清纯唯美| 国产伦精品一区二区三区四那| 国产高清不卡午夜福利| 亚洲精品,欧美精品| 久久婷婷人人爽人人干人人爱| 少妇裸体淫交视频免费看高清| 麻豆成人午夜福利视频| 国产乱人偷精品视频| 小蜜桃在线观看免费完整版高清| 欧美激情久久久久久爽电影| 久久精品人妻少妇| 亚洲国产欧洲综合997久久,| 91精品伊人久久大香线蕉| 午夜精品在线福利| 亚洲va在线va天堂va国产| 韩国av在线不卡| 看黄色毛片网站| 精品人妻偷拍中文字幕| 天堂√8在线中文| 丰满少妇做爰视频| 天堂影院成人在线观看| 日本黄大片高清| 男女边吃奶边做爰视频| 国产黄色视频一区二区在线观看 | 久久欧美精品欧美久久欧美| 亚洲国产欧洲综合997久久,| 插阴视频在线观看视频| 18禁动态无遮挡网站| 天天一区二区日本电影三级| 欧美日本亚洲视频在线播放| 91久久精品国产一区二区三区| 亚洲国产色片| 国产乱人偷精品视频| 日韩国内少妇激情av| 国产高潮美女av| 九九久久精品国产亚洲av麻豆| 狂野欧美激情性xxxx在线观看| 亚洲综合精品二区| 亚洲av不卡在线观看| 午夜视频国产福利| 九九在线视频观看精品| 看非洲黑人一级黄片| 激情 狠狠 欧美| 国产片特级美女逼逼视频| 直男gayav资源| 熟妇人妻久久中文字幕3abv| 国产女主播在线喷水免费视频网站 | 欧美成人午夜免费资源| 欧美+日韩+精品| 日韩精品青青久久久久久| 欧美三级亚洲精品| 蜜臀久久99精品久久宅男| 国产av不卡久久| 最后的刺客免费高清国语| 国产91av在线免费观看| 欧美激情国产日韩精品一区| 亚洲精品自拍成人| 国产av不卡久久| 99久久精品国产国产毛片| 国产在线一区二区三区精 | 成人毛片a级毛片在线播放| 国产私拍福利视频在线观看| 一二三四中文在线观看免费高清| 免费看av在线观看网站| 亚洲最大成人中文| 日韩中字成人| 狠狠狠狠99中文字幕| 成人亚洲精品av一区二区| 成人一区二区视频在线观看| 亚洲精品aⅴ在线观看| 精品久久国产蜜桃| 免费大片18禁| 能在线免费观看的黄片| 国产色婷婷99| 人人妻人人澡欧美一区二区| 91在线精品国自产拍蜜月| 婷婷色av中文字幕| 大香蕉久久网| 亚洲av成人精品一二三区| 午夜福利在线观看吧| 色尼玛亚洲综合影院| 波野结衣二区三区在线| 国产不卡一卡二| av卡一久久| 成人鲁丝片一二三区免费| 久久久久免费精品人妻一区二区| 日本黄色视频三级网站网址| 亚洲精品色激情综合| 亚洲国产精品成人久久小说| 日韩视频在线欧美| 国产熟女欧美一区二区| 国产高清视频在线观看网站| 91精品伊人久久大香线蕉| 久久精品国产亚洲av天美| 国产黄片美女视频| 亚洲一级一片aⅴ在线观看| 日韩亚洲欧美综合| 日韩视频在线欧美| 亚洲av中文av极速乱| 激情 狠狠 欧美| 亚洲成人久久爱视频| 国产欧美日韩精品一区二区| 永久网站在线| 最近最新中文字幕免费大全7| 国模一区二区三区四区视频| 亚洲国产精品久久男人天堂| 免费av不卡在线播放| 午夜日本视频在线| 偷拍熟女少妇极品色| 91精品国产九色| 国产精品.久久久| 亚洲一级一片aⅴ在线观看| 免费观看的影片在线观看| 亚洲在线观看片| 国产免费福利视频在线观看| 成年女人永久免费观看视频| 国产单亲对白刺激| 日韩国内少妇激情av| 99热网站在线观看| 国产精华一区二区三区| 午夜视频国产福利| 国产女主播在线喷水免费视频网站 | 国产日韩欧美在线精品| 欧美最新免费一区二区三区| 亚洲国产欧洲综合997久久,| 三级毛片av免费| 91狼人影院| 九九久久精品国产亚洲av麻豆| 久久久色成人| 中国美白少妇内射xxxbb| 久久国产乱子免费精品| 国模一区二区三区四区视频| 国产成人a区在线观看| 午夜福利在线观看免费完整高清在| 免费不卡的大黄色大毛片视频在线观看 | 亚洲av日韩在线播放| 中国国产av一级| 中文字幕亚洲精品专区| 色综合色国产| 精品久久久久久电影网 | www.av在线官网国产| 久久婷婷人人爽人人干人人爱| av卡一久久| 国产高清三级在线| 成人三级黄色视频| 午夜福利视频1000在线观看| 欧美一级a爱片免费观看看| 日韩欧美 国产精品| 乱人视频在线观看| 中文字幕制服av| 欧美日韩综合久久久久久| 亚洲国产欧洲综合997久久,| 麻豆成人午夜福利视频| 日韩一区二区视频免费看| 一区二区三区高清视频在线| 乱人视频在线观看| 午夜精品国产一区二区电影 | av国产免费在线观看| videossex国产| 99热精品在线国产| av播播在线观看一区| 91精品伊人久久大香线蕉| av在线老鸭窝| 久久欧美精品欧美久久欧美| 一区二区三区四区激情视频| 寂寞人妻少妇视频99o| 九九热线精品视视频播放| 国产91av在线免费观看| 日韩欧美三级三区| 国产精品一区www在线观看| 国产精品99久久久久久久久| 亚洲综合色惰| 一级黄片播放器| 亚洲电影在线观看av| 99九九线精品视频在线观看视频| 国产av一区在线观看免费| 日韩av在线免费看完整版不卡| 91精品伊人久久大香线蕉| 精品国内亚洲2022精品成人| 97人妻精品一区二区三区麻豆| 成人无遮挡网站| 麻豆一二三区av精品| 久久亚洲精品不卡| 乱系列少妇在线播放| 99视频精品全部免费 在线| 久久综合国产亚洲精品| 亚洲av成人精品一区久久| 亚洲欧美日韩无卡精品| 99久国产av精品国产电影| 国产精品乱码一区二三区的特点| 美女黄网站色视频| 国产成人aa在线观看| 精品少妇黑人巨大在线播放 | 1024手机看黄色片| 色5月婷婷丁香| 国产精品,欧美在线| 国产69精品久久久久777片| 久久精品91蜜桃| 欧美zozozo另类| 国产在线男女| 菩萨蛮人人尽说江南好唐韦庄 | videos熟女内射| 亚洲精品久久久久久婷婷小说 | 亚洲一区高清亚洲精品| 在线天堂最新版资源| 国产精品女同一区二区软件| 国产亚洲av片在线观看秒播厂 | 秋霞伦理黄片| 午夜a级毛片| 国产麻豆成人av免费视频| 亚洲经典国产精华液单| 男女国产视频网站| 亚洲国产精品国产精品| 欧美又色又爽又黄视频| 能在线免费看毛片的网站| 欧美日本视频| 美女黄网站色视频| 又爽又黄无遮挡网站| 我的老师免费观看完整版| 亚洲av熟女| 国产欧美日韩精品一区二区| 亚洲欧美成人综合另类久久久 | 精品一区二区三区人妻视频| 亚洲av熟女| 18+在线观看网站| 亚洲精品乱久久久久久| 女的被弄到高潮叫床怎么办| 国产爱豆传媒在线观看| 能在线免费看毛片的网站| 亚洲最大成人中文| 少妇的逼好多水| 亚洲不卡免费看| 国产三级在线视频| 在现免费观看毛片| 国产黄片美女视频| 观看美女的网站| 可以在线观看毛片的网站| 3wmmmm亚洲av在线观看| 久久久久性生活片| 亚洲精品aⅴ在线观看| 一边摸一边抽搐一进一小说| 两性午夜刺激爽爽歪歪视频在线观看| 一边亲一边摸免费视频| 国产乱人偷精品视频| 国产一区二区三区av在线| 麻豆成人午夜福利视频| 插逼视频在线观看| 99久久无色码亚洲精品果冻| 麻豆成人av视频| 久久久久久大精品| 中文字幕久久专区| 久久久成人免费电影| 国产v大片淫在线免费观看| 神马国产精品三级电影在线观看| 亚洲欧美日韩高清专用| 精品不卡国产一区二区三区| 一个人免费在线观看电影| 亚洲人成网站在线观看播放| 免费电影在线观看免费观看| 中文字幕亚洲精品专区| 日本五十路高清| 男人和女人高潮做爰伦理| 白带黄色成豆腐渣| 男女国产视频网站| 免费人成在线观看视频色| av免费在线看不卡| 国产精品永久免费网站| av国产免费在线观看| 麻豆av噜噜一区二区三区| 欧美又色又爽又黄视频| 国产淫片久久久久久久久| .国产精品久久| 国产高清有码在线观看视频| 亚洲欧美成人精品一区二区| 长腿黑丝高跟| 国产亚洲av片在线观看秒播厂 | 国产精品麻豆人妻色哟哟久久 | 亚洲欧美精品专区久久| 日韩 亚洲 欧美在线| 国产免费福利视频在线观看| 最近手机中文字幕大全| 日本黄大片高清| 精品人妻熟女av久视频| 又粗又爽又猛毛片免费看| 国产乱来视频区| 亚洲乱码一区二区免费版| 99在线视频只有这里精品首页| av天堂中文字幕网| av国产久精品久网站免费入址| 黄色一级大片看看| 狠狠狠狠99中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 天美传媒精品一区二区| 欧美激情国产日韩精品一区| 国产高清国产精品国产三级 | 亚洲精品乱码久久久v下载方式| 久久婷婷人人爽人人干人人爱| 男女那种视频在线观看| 97超视频在线观看视频| 真实男女啪啪啪动态图| 亚洲精品乱久久久久久| 精品人妻一区二区三区麻豆| 亚洲美女视频黄频| 欧美成人午夜免费资源| 色噜噜av男人的天堂激情| 18禁动态无遮挡网站| 国产一区二区在线观看日韩| 欧美成人免费av一区二区三区| 美女大奶头视频| 美女高潮的动态| av播播在线观看一区| 一级av片app| 国内揄拍国产精品人妻在线| 99热这里只有精品一区| 日韩精品有码人妻一区| 精品久久久久久成人av| 精品国产三级普通话版| av线在线观看网站| 国产成人免费观看mmmm| 久久久a久久爽久久v久久| 日韩三级伦理在线观看| 亚洲av成人精品一区久久| 长腿黑丝高跟| 成人午夜精彩视频在线观看| 91久久精品国产一区二区三区| 日韩在线高清观看一区二区三区| 国产精品精品国产色婷婷| 日韩成人伦理影院| 午夜亚洲福利在线播放| 亚洲国产精品成人综合色| 搞女人的毛片| 69av精品久久久久久| av天堂中文字幕网| 国产精华一区二区三区| 永久网站在线| 能在线免费观看的黄片| 日日啪夜夜撸| 特大巨黑吊av在线直播|