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

    高超聲速飛行器氣動(dòng)熱預(yù)測(cè)技術(shù)研究進(jìn)展

    2021-04-08 03:21:11喻成璋劉衛(wèi)華
    航空科學(xué)技術(shù) 2021年2期
    關(guān)鍵詞:數(shù)值計(jì)算試驗(yàn)研究

    喻成璋 劉衛(wèi)華

    摘要:氣動(dòng)熱預(yù)測(cè)技術(shù)是制約高超聲速飛行器發(fā)展的關(guān)鍵技術(shù)之一。飛行器在高速飛行過程中,氣動(dòng)加熱對(duì)其結(jié)構(gòu)強(qiáng)度影響顯著,嚴(yán)重時(shí)甚至可能導(dǎo)致結(jié)構(gòu)損傷,因此,為保障飛行器飛行安全,必須采取有效的熱防護(hù)措施,而掌握氣動(dòng)熱變化規(guī)律是合理設(shè)計(jì)高超聲速飛行器熱防護(hù)措施的基礎(chǔ),它對(duì)于飛行器結(jié)構(gòu)設(shè)計(jì)、材料選擇均有重要的指導(dǎo)意義。本文從試驗(yàn)、工程計(jì)算與數(shù)值仿真三個(gè)方面系統(tǒng)地歸納、總結(jié)國(guó)內(nèi)外學(xué)者在氣動(dòng)熱預(yù)測(cè)方面的研究成果,并展望其未來的發(fā)展,以期為國(guó)內(nèi)高超聲速飛行器的研制工作提供有益參考與借鑒。

    關(guān)鍵詞:氣動(dòng)熱;高超聲速飛行器;試驗(yàn)研究;工程算法;數(shù)值計(jì)算

    中圖分類號(hào):V434+.11文獻(xiàn)標(biāo)識(shí)碼:ADOI:10.19452/j.issn1007-5453.2021.02.002

    基金項(xiàng)目:國(guó)家自然科學(xué)基金民航聯(lián)合基金(U1933121);中央高?;究蒲袠I(yè)務(wù)基金(NC2020001);江蘇高校優(yōu)勢(shì)學(xué)科建設(shè)工程

    因氣動(dòng)熱所引發(fā)的熱障問題已成為制約高超聲速飛行器發(fā)展的瓶頸。當(dāng)飛行器作超聲速飛行時(shí),機(jī)身表面與氣流發(fā)生劇烈的摩擦,空氣受到阻滯和壓縮導(dǎo)致溫度急劇上升,對(duì)飛行器形成強(qiáng)烈的熱障[1-2];當(dāng)飛行速度進(jìn)一步提高到高超聲速時(shí),空氣流動(dòng)特性還將發(fā)生本質(zhì)的改變,其物理現(xiàn)象主要體現(xiàn)為薄激波層、熵層、黏性干擾、低密度性及高溫激波層內(nèi)的真實(shí)氣體效應(yīng)[3]。以馬赫數(shù)Ma>20的再入飛行器為例,其前緣駐點(diǎn)溫度可高達(dá)1×104K以上,在此高溫環(huán)境下,不僅其周圍空氣會(huì)發(fā)生電離,而且飛行器表面材料也將發(fā)生化學(xué)反應(yīng),對(duì)飛行器造成不可恢復(fù)的損傷。因此,為保證飛行安全,必須采取相應(yīng)的熱防護(hù)措施以保證飛行器結(jié)構(gòu)和內(nèi)部設(shè)置正常工作。

    掌握氣動(dòng)熱變化規(guī)律是高超聲速飛行器熱防護(hù)的關(guān)鍵技術(shù),但氣動(dòng)熱的影響因素眾多,要想實(shí)現(xiàn)氣動(dòng)熱精確的預(yù)測(cè)十分困難。從20世紀(jì)50年代開始,人們分別采用試驗(yàn)研究、工程計(jì)算,以及數(shù)值仿真等技術(shù)手段針對(duì)氣動(dòng)熱問題開展了大量的研究工作,研究?jī)?nèi)容也隨著高超聲速飛行器的發(fā)展不斷深化,本文旨在從這三個(gè)方面系統(tǒng)地總結(jié)當(dāng)前的研究進(jìn)展,并明確后續(xù)工作的重點(diǎn)。

    1試驗(yàn)研究

    熱流測(cè)試技術(shù)的進(jìn)步為氣動(dòng)熱試驗(yàn)研究奠定了良好的基礎(chǔ),目前,氣動(dòng)熱測(cè)量的地面試驗(yàn)都是以各類激波風(fēng)洞和常規(guī)高超聲速風(fēng)洞為手段,其試驗(yàn)研究?jī)?nèi)容主要涉及迎角、鈍比率、表面熱流分布、激波形狀、邊界層轉(zhuǎn)捩和非平衡效應(yīng)等因素。

    20世紀(jì)60年代,美國(guó)國(guó)家航空航天局(NASA)對(duì)15°鈍錐體在Ma10.6下的氣動(dòng)熱問題進(jìn)行了試驗(yàn)研究,結(jié)果表明,在迎角為0°時(shí),增大鈍比率可降低熱流密度并阻止邊界層轉(zhuǎn)捩[4]。

    1992年,King[5]對(duì)5°角錐在Ma3.5下的邊界層轉(zhuǎn)捩位置進(jìn)行了測(cè)量,結(jié)果表明,隨著橫向流動(dòng)的增加,轉(zhuǎn)捩點(diǎn)的位置發(fā)生了明顯改變。

    1998年,S.A. Berry[6]等就X-34模型進(jìn)行了氣動(dòng)熱試驗(yàn)及氣動(dòng)特性分析,借助于磷光測(cè)量技術(shù)實(shí)測(cè)了0°~35°迎角范圍內(nèi)模型表面的熱流分布、流線形式與激波形狀,結(jié)果顯示,當(dāng)迎角從8°增至23°時(shí),迎風(fēng)面熱流略有增加,弓形激波與機(jī)翼的相互作用逐漸向機(jī)翼內(nèi)側(cè)移動(dòng);且隨著馬赫數(shù)的不斷增大,迎風(fēng)面上的流動(dòng)狀態(tài)從層流轉(zhuǎn)向湍流,中心線上的轉(zhuǎn)捩位置不斷向前緣移動(dòng)。

    2000年,Nonaka[7]對(duì)速度2.44 ~3.85km/s、壓力0.56~ 20kPa范圍的球體繞流進(jìn)行了試驗(yàn)研究,采用陰影法和紋影法對(duì)球體表面的激波形狀及流場(chǎng)進(jìn)行了可視化處理,試驗(yàn)結(jié)果清楚地展示了凍結(jié)流、平衡流、非平衡流之間激波脫體距離的差距,表明了非平衡效應(yīng)對(duì)流場(chǎng)的影響。

    2018年,Borg[8]等在Ma6的靜風(fēng)洞中對(duì)飛行姿態(tài)下的HIFiRE-5縮比模型進(jìn)行了測(cè)試,結(jié)果表明,偏航角對(duì)橫向流失穩(wěn)的影響要大于迎角的影響,在低噪聲水平下,改變1°偏航角可使得最小橫向流轉(zhuǎn)捩雷諾數(shù)減小20%。

    事實(shí)上,邊界層轉(zhuǎn)捩和非平衡效應(yīng)都是飛行器熱防護(hù)系統(tǒng)設(shè)計(jì)需要面對(duì)的難題,其中邊界層轉(zhuǎn)捩將產(chǎn)生復(fù)雜的、不確定的縱橫向氣流擾動(dòng)及氣動(dòng)加熱,但由于風(fēng)洞試驗(yàn)中各種擾動(dòng)的存在,針對(duì)邊界層轉(zhuǎn)捩的試驗(yàn)研究結(jié)果與實(shí)際情況存在較大差距;同時(shí),非平衡效應(yīng)試驗(yàn)條件的不易滿足也制約了氣動(dòng)熱試驗(yàn)研究的發(fā)展,雖然有學(xué)者進(jìn)行了一些嘗試,但目前這方面的研究無論是質(zhì)量上還是數(shù)量上都無法滿足發(fā)展的需求,未來還需對(duì)試驗(yàn)設(shè)備及試驗(yàn)?zāi)芰ψ鲞M(jìn)一步的完善。

    由于受限于風(fēng)洞試驗(yàn)設(shè)備的模擬能力,地面試驗(yàn)在模型尺寸、雷諾數(shù)和真實(shí)氣體效應(yīng)方面均有一定的限制,無法完整地對(duì)各類飛行器真實(shí)氣動(dòng)熱環(huán)境進(jìn)行模擬,為此,歐美等國(guó)家也開展了大量的飛行試驗(yàn)研究,如1968年NASA進(jìn)行的Reentry F項(xiàng)目[9]、1994年日本HOPE-X計(jì)劃中的軌道再入飛行試驗(yàn)[10]等。雖然飛行試驗(yàn)?zāi)軌颢@得較為準(zhǔn)確的氣動(dòng)熱數(shù)據(jù),但其周期長(zhǎng)、難度大、耗資巨的特點(diǎn)也決定了它不可能成為氣動(dòng)熱預(yù)測(cè)研究的主要手段,目前飛行試驗(yàn)數(shù)據(jù)大多用于數(shù)值仿真結(jié)果驗(yàn)證、模型完善及飛行器設(shè)計(jì)的最終校核。

    2工程計(jì)算方法

    雖然氣動(dòng)熱的復(fù)雜性使得其尚無普適的理論預(yù)測(cè)方法,但工程計(jì)算以其獨(dú)特的優(yōu)勢(shì),在高超聲速飛行器概念設(shè)計(jì)階段扮演著重要角色。工程算法將擬設(shè)計(jì)的飛行器劃分為無迎角和有迎角兩種類型。

    對(duì)于無迎角的氣動(dòng)熱預(yù)測(cè)問題,工程算法將其劃分為駐點(diǎn)、非駐點(diǎn)層流區(qū)、非駐點(diǎn)湍流區(qū)分別進(jìn)行處理。其中,駐點(diǎn)熱流密度采用邊界層相似解的近似分析方法或風(fēng)洞試驗(yàn)數(shù)據(jù)關(guān)聯(lián)分析獲得的經(jīng)驗(yàn)公式來進(jìn)行計(jì)算;非駐點(diǎn)層流流動(dòng)則假定各熱力學(xué)參數(shù)沿物面切向的變化率遠(yuǎn)小于沿法向的變化率,再通過數(shù)學(xué)變換得到局部相似解,從而推導(dǎo)出層流的熱流密度分布公式;非駐點(diǎn)湍流區(qū)主要通過動(dòng)量積分方程求解出局部摩擦因數(shù),再根據(jù)雷諾比獲得傳熱系數(shù)進(jìn)而確定湍流邊界層加熱情況[11]。

    由于不同流動(dòng)區(qū)域的氣動(dòng)熱計(jì)算方法不同,邊界層轉(zhuǎn)捩位置的預(yù)測(cè)將直接決定工程計(jì)算的準(zhǔn)確性,對(duì)此,人們開發(fā)了多種方法來對(duì)轉(zhuǎn)捩點(diǎn)位置進(jìn)行預(yù)測(cè),如經(jīng)驗(yàn)方法、eN方法、轉(zhuǎn)捩模式等。值得注意的是,雖然目前大多數(shù)計(jì)算都采用了地面和飛行試驗(yàn)獲得的轉(zhuǎn)捩準(zhǔn)則來對(duì)轉(zhuǎn)捩點(diǎn)的位置進(jìn)行判斷,然而,不同轉(zhuǎn)捩準(zhǔn)則下氣動(dòng)熱的結(jié)果存在一定偏差。

    對(duì)于有迎角情況,則必須首先根據(jù)飛行器的氣動(dòng)外形分為小寬鈍比與大寬鈍比兩種預(yù)測(cè)類型,然后再變換至零迎角條件下進(jìn)行氣動(dòng)熱的計(jì)算。

    小寬鈍比氣動(dòng)熱工程預(yù)測(cè)方法主要有等價(jià)錐法和軸對(duì)稱比擬法。等價(jià)錐法是將有迎角的錐體流動(dòng)用零迎角的等價(jià)錐替代,適合于小迎角計(jì)算,并對(duì)于迎風(fēng)面母線和背風(fēng)面母線處的熱流均有較為準(zhǔn)確的預(yù)測(cè)能力,但在處理非軸對(duì)稱體時(shí)則存在較大偏差[12];軸對(duì)稱比擬法基于小橫向流假設(shè),通過Manger變化將三維邊界層方程簡(jiǎn)化為準(zhǔn)軸對(duì)稱邊界層方程,再通過一定條件將各條流線上的熱流密度與零迎角下軸對(duì)稱物體的熱流密度聯(lián)系起來,雖然軸對(duì)稱比擬法被認(rèn)為是當(dāng)前工程上計(jì)算三維邊界層傳熱問題最有效的方法之一,但該方法中間量迭代計(jì)算復(fù)雜,通用性較差[13]。

    大寬鈍比氣動(dòng)熱預(yù)測(cè)主要采用了“片條理論”,它把飛行器切成二維的片條,并假定這些片條之間沒有干擾,然后將這些二維片條作為二維鈍頭體來計(jì)算其表面熱流,再綜合三維效應(yīng)來對(duì)結(jié)果進(jìn)行修正。

    對(duì)于稀薄流動(dòng)氣動(dòng)熱的預(yù)測(cè),橋函數(shù)法是目前應(yīng)用最為廣泛的工程計(jì)算方法,它將連續(xù)流方法預(yù)測(cè)值和自由分子流方法預(yù)測(cè)值通過函數(shù)橋連接獲得近似熱流,其中Mattinng[14]和Nomura[15]建立的橋函數(shù)最為常用。

    橋函數(shù)法多用于駐點(diǎn)熱流密度以及鈍錐體表面熱流密度計(jì)算[16-18],具有較高的計(jì)算精度,但也有研究表明,橋函數(shù)法對(duì)尖長(zhǎng)體及細(xì)長(zhǎng)體表面熱流密度的預(yù)測(cè)并不理想[19]。

    事實(shí)上,工程計(jì)算方法的不斷完善也催生了工程預(yù)測(cè)軟件的開發(fā)。早期著名氣動(dòng)熱計(jì)算軟件MINIVER[20]采用了常見的熱流公式來對(duì)駐點(diǎn)、層流區(qū)、湍流區(qū)分別進(jìn)行氣動(dòng)熱計(jì)算,它適用于完全氣體和平衡氣體兩種情況,但三維計(jì)算有一定的局限性,且無法計(jì)算鈍錐模型的下游效應(yīng)[21]。

    AEROHEAT是人們?cè)缙诟鶕?jù)軸對(duì)稱比擬法開發(fā)的一款經(jīng)典的氣動(dòng)熱計(jì)算程序;Zoby等采用適體坐標(biāo)系,開發(fā)了LATCH計(jì)算程序,與AEROHEAT相比,它增強(qiáng)了處理復(fù)雜外形氣動(dòng)熱的能力;2006年,Hamilton等對(duì)LATCH算法進(jìn)行改進(jìn),發(fā)展了三角形結(jié)構(gòu)化網(wǎng)格的UNLATCH2算法,它改善了LATCH在結(jié)構(gòu)化網(wǎng)格計(jì)算無黏流場(chǎng)方面的局限性[22]。

    3數(shù)值模擬技術(shù)

    由于工程算法對(duì)飛行器外形的適用性較差,對(duì)轉(zhuǎn)捩、激波以及黏性邊界層等物理現(xiàn)象難以準(zhǔn)確描述,因此,早期人們只能通過試驗(yàn)獲得復(fù)雜流動(dòng)的流動(dòng)細(xì)節(jié),周期長(zhǎng)且耗費(fèi)巨大。20世紀(jì)80年代后,數(shù)值模擬技術(shù)的長(zhǎng)足進(jìn)步改變了這種不利狀況,通過求解控制方程繼而獲得表面熱流分布情況,不僅花費(fèi)更少、效率更高,且對(duì)流動(dòng)細(xì)節(jié)的刻畫也更為細(xì)致,目前數(shù)值模擬方法已經(jīng)成為氣動(dòng)熱預(yù)測(cè)的重要研究手段。

    不同流動(dòng)區(qū)域的氣動(dòng)特性不同,因此當(dāng)飛行器處于連續(xù)流區(qū)和稀薄流區(qū)等不同流動(dòng)區(qū)域時(shí),其數(shù)值計(jì)算方法也有所差異[23]。

    3.1連續(xù)流區(qū)氣動(dòng)熱數(shù)值模擬

    在連續(xù)流區(qū)對(duì)N-S方程及其各種近似方程的求解是獲取熱流分布的主要途徑,隨著計(jì)算機(jī)性能的不斷提高,氣動(dòng)熱數(shù)值模擬經(jīng)歷了邊界層方程、黏性激波層方程、拋物化N-S方程、全N-S方程等不同的發(fā)展歷程。

    1982年,Hamilton[24]發(fā)展了一種求解二階邊界層方程的計(jì)算方法,并以此獲得了飛行器熱流分布情況,計(jì)算結(jié)果表明邊界層方程能夠較好地模擬迎角為25°~ 40°的流動(dòng)狀況。

    邊界層理論為氣動(dòng)熱研究提供了有效的計(jì)算方法,但由于邊界層方程本身的局限性,它無法反映出無黏流和黏流之間的相互干擾機(jī)理,因此邊界層方程僅適于處理一些馬赫數(shù)不太高的簡(jiǎn)單工況;為了精確地考慮黏性作用,有學(xué)者將激波和物面之間的整個(gè)流場(chǎng)作為全黏性處理,建立了統(tǒng)一的控制方程,形成了黏性激波層理論。

    1992年,歐陽(yáng)水吾[25]等通過對(duì)多組元黏性激波層方程的求解,就包含鈍頭體氣動(dòng)加熱率在內(nèi)等多個(gè)問題進(jìn)行分析計(jì)算,所得出的各物理量變化規(guī)律表明了該計(jì)算方法的合理性。但值得注意的是,由于高雷諾數(shù)流動(dòng)的物面附近物理量變化非常劇烈,黏性激波層方程的求解需要采用許多特殊的數(shù)值方法,計(jì)算過程復(fù)雜,通用性較差。

    由定常N-S方程舍去含有流向?qū)?shù)黏性項(xiàng)的拋物化N-S方程(PN-S方程),在各類高超聲速飛行器繞流問題中獲得了較好的應(yīng)用,如20世紀(jì)90年代,Lawrance[26]將TVD格式應(yīng)用于PN-S方程的空間推進(jìn)算法中,對(duì)半角為10°的圓錐流動(dòng)進(jìn)行求解,給出了三個(gè)不同迎角下圓錐的熱流密度,計(jì)算結(jié)果與試驗(yàn)結(jié)果較好吻合;2002年,NASA通過求解PN-S方程獲得了X-43前機(jī)身,以及整流罩前端的氣動(dòng)熱分布情況[27],進(jìn)一步豐富了PN-S方程數(shù)值解法在高超聲速飛行器領(lǐng)域的應(yīng)用。

    雖然通過簡(jiǎn)化N-S方程獲取氣動(dòng)熱分布的PN-S方法在早期取得了豐富的研究成果,但與全N-S方程相比,兩者還是存在一定的誤差,隨著計(jì)算機(jī)性能的提升,對(duì)全N-S方程進(jìn)行求解成為現(xiàn)實(shí),該方法越來越多地被應(yīng)用到氣動(dòng)熱數(shù)值計(jì)算中。在21世紀(jì)初期,賀國(guó)宏[28]通過求解N-S方程對(duì)鈍錐和鈍雙錐流動(dòng)進(jìn)行了數(shù)值模擬,計(jì)算結(jié)果不僅精確地模擬了迎風(fēng)區(qū)的熱流分布,還準(zhǔn)確地刻畫了鈍錐背風(fēng)區(qū)氣流分離、雙錐交接區(qū)域熱流變化等一系列現(xiàn)象。

    事實(shí)上,除控制方程外,氣動(dòng)熱數(shù)值模擬計(jì)算的結(jié)果還受到差分格式、網(wǎng)格分布、限制器和湍流模型等多個(gè)因素的影響,為此,國(guó)內(nèi)外學(xué)者針對(duì)這些影響因素開展了較為系統(tǒng)的研究工作。

    當(dāng)流動(dòng)中存在間斷時(shí),中心型格式通常需要添加人工黏性來抑制振蕩,它極易造成污染、影響計(jì)算精度,因此從20世紀(jì)80年代開始,各類迎風(fēng)格式成為數(shù)值計(jì)算的主流格式。

    2003年,李君哲[29]等分別選用FDS格式、FVS格式、AUSM+三種迎風(fēng)格式以及一種中心差分格式對(duì)二維圓柱繞流和鈍雙錐繞流進(jìn)行計(jì)算,結(jié)果表明,三種迎風(fēng)格式的數(shù)值計(jì)算精度明顯優(yōu)于中心差分格式,其中AUSM+格式和FDS格式更為逼近試驗(yàn)數(shù)據(jù)。

    氣動(dòng)熱計(jì)算網(wǎng)格的選擇十分苛刻,有時(shí)需要根據(jù)數(shù)值計(jì)算結(jié)果對(duì)網(wǎng)格進(jìn)行多次調(diào)整,這也是數(shù)值計(jì)算過程中較為耗時(shí)的部分。

    物面第一層網(wǎng)格高度對(duì)數(shù)值計(jì)算結(jié)果影響較大,許多文獻(xiàn)將其作為重要參數(shù)加以分析。1998年,Lee[30]等在鈍頭體氣動(dòng)熱計(jì)算中探討了網(wǎng)格相關(guān)性問題,結(jié)果表明,在其他條件均不改變的情況下,壁面網(wǎng)格雷諾數(shù)擴(kuò)大一倍,計(jì)算結(jié)果偏差值能達(dá)到20%以上。由此可見,氣動(dòng)熱能否精確求解,與近壁面的網(wǎng)格尺寸有著直接關(guān)系。

    2016年,張翔[31]以二維圓柱繞流為例分析了網(wǎng)格雷諾數(shù)對(duì)熱流計(jì)算結(jié)果的影響。結(jié)果表明,氣動(dòng)熱計(jì)算中通常需要保持網(wǎng)格雷諾數(shù)小于10這一最低要求,而網(wǎng)格雷諾數(shù)小于8即可獲得收斂結(jié)果。目前人們廣泛根據(jù)網(wǎng)格雷諾數(shù)來確定第一層計(jì)算網(wǎng)格高度,但也有研究表明這種方法并不總是可靠的,例如,Gao[32]在圓柱體繞流計(jì)算中發(fā)現(xiàn),根據(jù)網(wǎng)格雷諾數(shù)確定的第一層網(wǎng)格高度偏大,若要達(dá)到計(jì)算要求至少還需要減小一個(gè)數(shù)量級(jí)。

    為了避免結(jié)果振蕩,計(jì)算過程中通常需引入限制器對(duì)插值梯度進(jìn)行限制[33],各種限制器的精度和耗散性不同,曾有學(xué)者認(rèn)為限制器對(duì)計(jì)算精度的影響甚至要大于計(jì)算格式。

    楊建龍[34]分別采用minmod、Van Leer和Osher-C三種限制器對(duì)雙錐模型外部流場(chǎng)進(jìn)行了數(shù)值模擬,結(jié)果表明,采用minmod限制器獲得的最大熱流值與熱流分布特性明顯優(yōu)于其余兩種限制器。Zhang[35]將傳統(tǒng)的MUSCL型限制器與多維限制器進(jìn)行對(duì)比分析,結(jié)果表明,Van Leer限制器是傳統(tǒng)MUSCL型限制器中考慮魯棒性和準(zhǔn)確性的最佳選擇,而多維限制器則在這兩方面比其他所有傳統(tǒng)限制器表現(xiàn)得更為優(yōu)秀。趙雅甜[36]發(fā)展了一種新型三階TVD限制器并對(duì)其進(jìn)行了性能分析,結(jié)果表明,該限制器在復(fù)雜外形氣動(dòng)熱算例中均表現(xiàn)出更好的適應(yīng)性以及氣動(dòng)熱預(yù)測(cè)能力,與double minmod限制器相比,計(jì)算精度更高;與minmod限制器相比,具有良好的間斷分辨率且避免了過多的數(shù)值耗散。

    自流體力學(xué)問世以來,湍流一直是各類流動(dòng)問題的難點(diǎn),但到目前為止,人們對(duì)于湍流氣動(dòng)熱的預(yù)測(cè)結(jié)果并不理想[37-38],如GASP程序?qū)-33模型湍流中心線上的計(jì)算結(jié)果就比試驗(yàn)值偏小30% ~ 40%。目前湍流的數(shù)值模擬方法包括雷諾平均N-S方程求解、大渦模擬、分離渦流模擬以及直接數(shù)值模擬等,但考慮到魯棒性以及計(jì)算效率的限制,求解雷諾平均N-S方程(RAN-S)仍是目前采用的主要研究手段。

    隨著數(shù)值計(jì)算方法日趨成熟,對(duì)于復(fù)雜流動(dòng)的工程算法逐步讓位于數(shù)值計(jì)算,因此,開發(fā)專業(yè)計(jì)算軟件成為了熱門,對(duì)此,NASA Langley研究中心、Ames研究中心及波音公司做了大量工作,其中比較著名的有LAURA[39]程序和GASP[40]程序,其有效性已被多個(gè)飛行試驗(yàn)以及地面試驗(yàn)所證實(shí)[41]。目前國(guó)外這些氣動(dòng)熱數(shù)值計(jì)算軟件已較為成熟,對(duì)于特定構(gòu)型有較高的計(jì)算精度,相比之下,我國(guó)仍缺乏高效、可靠的氣動(dòng)熱模擬工具。

    3.2稀薄流區(qū)氣動(dòng)熱數(shù)值模擬

    從6×104m高空開始,氣體性質(zhì)逐漸偏離連續(xù)介質(zhì),飛行器的氣動(dòng)特性亦發(fā)生了顯著變化,在各類稀薄流動(dòng)中,自由分子流由于可以忽略氣體速度分布函數(shù)的變化,通??梢愿鶕?jù)氣體分子動(dòng)理論獲得解析解[42],而滑移流區(qū)以及過渡流區(qū)的氣動(dòng)熱預(yù)測(cè)則較為困難。

    目前針對(duì)滑移流區(qū)和過渡流區(qū)的數(shù)值模擬方法主要有兩種:采用滑移邊界條件計(jì)算流體力學(xué)(CFD)方法和基于分子動(dòng)力學(xué)的直接模擬蒙特卡羅(DSMC)法。

    在滑移流域,其主流場(chǎng)中的流動(dòng)依然適用于連續(xù)介質(zhì)假設(shè),因此相較于復(fù)雜的DSMC算法,人們更傾向于應(yīng)用滑移條件來對(duì)N-S方程的應(yīng)用范圍進(jìn)行拓展。目前廣泛采用的有Maxwell滑移條件、Gokcen滑移條件以及Lockerby滑移條件。Lofthouse[43]分析了該方法在滑移流域的適用性并對(duì)三種滑移條件進(jìn)行比較,結(jié)果表明,當(dāng)克努森數(shù)Kn∞≤0.25時(shí),帶滑移邊界的CFD方法與DSMC方法預(yù)測(cè)值誤差<5%,而且極大地縮短了計(jì)算時(shí)間。在三種滑移條件中,Gocken條件所得到的計(jì)算結(jié)果最好,但也最為耗時(shí)。

    DSMC方法是由Bird[44]首先提出的。該方法采用了大量模擬分子代替真實(shí)氣體分子來進(jìn)行分子運(yùn)動(dòng)和碰撞過程的解耦運(yùn)算,待流場(chǎng)中分子數(shù)量趨于穩(wěn)定后再通過統(tǒng)計(jì)采樣的方式獲得宏觀計(jì)算結(jié)果。由于DSMC方法是從微觀氣體分子入手,本質(zhì)上更接近于真實(shí)的氣體運(yùn)動(dòng)狀態(tài)和能量交換過程,因此在處理稀薄流氣動(dòng)熱問題上具有天然的優(yōu)勢(shì)。

    2013年,Prasanth[45]對(duì)DSMC方法的限制條件進(jìn)行了系統(tǒng)的總結(jié),提出了計(jì)算結(jié)果的準(zhǔn)確性必須滿足的三個(gè)條件:網(wǎng)格小于分子的平均自由程、碰撞網(wǎng)格大小(每個(gè)網(wǎng)格中模擬的分子數(shù))滿足計(jì)算需求以及時(shí)間步長(zhǎng)小于分子的平均碰撞時(shí)間。依據(jù)Prasanth研究結(jié)果,小Kn流動(dòng)需要采用高精度網(wǎng)格并消耗大量的計(jì)算資源,所以DSMC方法多用于稀薄程度較大的過渡流區(qū)。

    DSMC方法的有效性已被多次證明并取得了一系列成果,例如,2007年,Liechty[46]運(yùn)用DSMC方法計(jì)算了火星偵察軌道飛行器(MRO)在不同迎角、側(cè)滑角下的傳熱系數(shù),計(jì)算結(jié)果與飛行試驗(yàn)結(jié)果較為吻合,證實(shí)了該方法在過渡流區(qū)的有效性;2015年,Sengil[47]采用二維DSMC求解器分析了不同幾何外形對(duì)氣動(dòng)特性的影響;2019年,NASA采用DSMC并行計(jì)算軟件DAC對(duì)空心圓柱和空心方形棱柱的氣動(dòng)熱環(huán)境進(jìn)行模擬,并對(duì)迎角、壁厚以及Kn等重要參數(shù)對(duì)氣動(dòng)加熱的影響進(jìn)行了分析[48]。

    目前,DSMC方法已被公認(rèn)為處理稀薄流動(dòng)最成功的辦法,但龐大的計(jì)算量限制了其發(fā)展,為此,許多學(xué)者對(duì)其進(jìn)行了優(yōu)化,在諸多優(yōu)化算法中以DSMC/N-S耦合算法最為常見。2005年,Lian[49]提出一種基于非結(jié)構(gòu)網(wǎng)格的DSMC/N-S耦合算法,它不僅能夠細(xì)致刻畫連續(xù)介質(zhì)失效帶來的一系列物理現(xiàn)象,而且能夠?qū)πㄐ伪砻媲熬壍姆瞧胶庑?yīng)進(jìn)行準(zhǔn)確的預(yù)測(cè);2018年,Xu[50]提出了一種基于自適應(yīng)結(jié)構(gòu)化/非結(jié)構(gòu)化重疊網(wǎng)格的耦合算法,該方法綜合了結(jié)構(gòu)化網(wǎng)格和非結(jié)構(gòu)化網(wǎng)格的優(yōu)點(diǎn),不僅可以高精度地研究流動(dòng)中的各種特性,在計(jì)算速度方面比全DSMC方法提高了近三倍左右。

    除了對(duì)耦合算法進(jìn)行改進(jìn),完善連續(xù)介質(zhì)失效邊界的準(zhǔn)確判定和分析連續(xù)假設(shè)失效所帶來的影響也是研究熱點(diǎn)。由于連續(xù)介質(zhì)失效準(zhǔn)則與耦合算法的計(jì)算效率息息相關(guān),為此,研究者提出了一系列準(zhǔn)則數(shù)來對(duì)失效邊界進(jìn)行準(zhǔn)確劃分,如P參數(shù)、當(dāng)?shù)乜伺瓟?shù)Knl、局部克努森數(shù)Knp等。多數(shù)情況下采用的是Knp作為判定準(zhǔn)則,當(dāng)Knp≥0.05時(shí)則認(rèn)為連續(xù)方程失效[51];但也有研究表明,該結(jié)論并不一定適用。2019年,午辛暄[52]對(duì)不同物理問題所體現(xiàn)的失效準(zhǔn)則數(shù)進(jìn)行了區(qū)分,其結(jié)論是,對(duì)于平板流動(dòng)結(jié)構(gòu)采用全局Kn作為失效準(zhǔn)則更為合適,其失效大小為Kn≥0.01;而采用Knp則不夠準(zhǔn)確。

    連續(xù)假設(shè)失效的影響通過量化CFD方法計(jì)算誤差來體現(xiàn)。Lofthouse[53]分別采用CFD方法和DSMC方法對(duì)不同稀薄程度的氬氣圓柱繞流進(jìn)行模擬。結(jié)果表明,隨著稀薄程度的不斷增大,兩種方法對(duì)氣動(dòng)熱峰值的預(yù)測(cè)差異從Kn = 0.002時(shí)不足1%逐漸增加到Kn = 0.25時(shí)的30%以上。連續(xù)介質(zhì)失效的研究對(duì)耦合算法有重要意義,隨著該部分理論的逐漸完善,算法的計(jì)算效率將進(jìn)一步提高。

    目前帶滑移邊界條件的CFD方法以及DSMC類算法在稀薄氣體氣動(dòng)熱中廣泛使用,但也有學(xué)者嘗試通過高階Boltzmann方程來獲得稀薄流動(dòng)氣動(dòng)熱的分布情況,如Burnett方程[54]及Eu方程[55]等。雖然有相關(guān)研究表明這些方程在計(jì)算中有一定的優(yōu)勢(shì),但一些難以解決的問題仍限制著其在氣動(dòng)熱領(lǐng)域的應(yīng)用。隨著該部分理論的進(jìn)一步完善,拓展的流體動(dòng)力學(xué)方程也許將成為預(yù)測(cè)稀薄流區(qū)氣動(dòng)熱的又一重要工具。

    4結(jié)束語

    目前的風(fēng)洞試驗(yàn)設(shè)備無法對(duì)真實(shí)飛行環(huán)境進(jìn)行高品質(zhì)的復(fù)現(xiàn),尤其是對(duì)稀薄流動(dòng)以及非平衡效應(yīng)難以模擬,因此,對(duì)現(xiàn)有風(fēng)洞進(jìn)行改造升級(jí),拓展其功能,研制開發(fā)一些新概念風(fēng)洞來適用于未來高超聲速飛行器發(fā)展將勢(shì)在必行。

    (1)發(fā)展可適用于不同幾何體稀薄流動(dòng)的橋函數(shù),完善與拓展工程計(jì)算方法的使用范圍是有現(xiàn)實(shí)意義的。

    (2)采用高分辨率、高精度的差分形式求解全N-S方程是連續(xù)流區(qū)氣動(dòng)熱預(yù)測(cè)技術(shù)的發(fā)展方向,同時(shí)計(jì)算格式、網(wǎng)格效應(yīng)、限制器以及湍流模型等數(shù)值計(jì)算理論的完善將進(jìn)一步增強(qiáng)對(duì)復(fù)雜工況的氣動(dòng)熱預(yù)測(cè)的準(zhǔn)確性。

    (3)稀薄流動(dòng)的連續(xù)假設(shè)失效對(duì)氣動(dòng)熱預(yù)測(cè)技術(shù)提出了新的挑戰(zhàn),其中滑移流區(qū)和過渡流區(qū)的氣動(dòng)熱預(yù)測(cè)較為困難,因此繼續(xù)完善滑移邊界條件以及DSMC方法將會(huì)是未來稀薄流動(dòng)氣動(dòng)熱發(fā)展的重點(diǎn)。

    參考文獻(xiàn)

    [1]趙旭升,楊志斌,張新霞.一種一體化熱防護(hù)系統(tǒng)綜合效能量化評(píng)價(jià)方法[J].航空科學(xué)技術(shù),2019,30(5):65-69. Zhao Xusheng, Yang Zhibin, Zhang Xinxia. A quantitative evaluation method for comprehensive efficiency of integrated thermalprotectionsystem[J].AeronauticalScience& Technology, 2019, 30(5): 65-69.(in Chinese)

    [2]羅楚養(yǎng),張朋,李偉東,等.高溫復(fù)合材料在空空導(dǎo)彈上的應(yīng)用研究[J].航空科學(xué)技術(shù),2017,28(1):19-24. Luo Chuyang, Zhang Peng, Li Weidong, et al. Application research of high temperature composite on airborne missile[J]. Aeronautical Science & Technology. 2017, 28(1): 19-24.(in Chinese)

    [3]Anderson J D. Hypersonic and high-temperature gas dynamics[R].Blackburg,Virgina,US:AmericanInstituteof Aeronautics andAstronautics,2006.

    [4]Cleary J W. Effects of angle of attack and bluntness on laminar heating-rate distributions of a 15 cone at a Mach number of 10.6[R]. NASATND-5450,1969:1-4.

    [5]King R A. Three-dimensional boundary-layer transition on a cone at Mach 3.5[J]. Experiments in Fluids,1992,13(5):305-314.

    [6]Berry S A,Horvath T J,DiFulvio M,et al. X-34 experimental aeroheating at Mach 6 and 10[J]. Journal of Spacecraft and Rockets,1999,36(2):171-178.

    [7]Nonaka S,Mizuno H,Takayama K,et al. Measurement of shock standoff distance for sphere in ballistic range[J]. Journal of Thermophysics and Heat Transfer,2000,14(2):225-229.

    [8]Borg M P,Kimmel R L. Ground test of transition for HIFiRE-5b at flight-relevant attitudes[J]. Journal of Spacecraft and Rockets,2018,55(6):1329-1340.

    [9]Sinha K,Barnhardt M,Candler G. Detached eddy simulation of hypersonic base flows with application to fire ii experiments[C]//34th AIAA Fluid Dynamics Conference and Exhibit,2004:2633.

    [10]Yanagihara M,Miyazawa Y,Alimoto T,et al. HOPE-X high speed flight demonstration program phase II[C]//10th AIAA/NALNASDA-ISASInternationalSpacePlanesandHypersonic Systems and Technologies Conference,2001:1805.

    [11]彭治雨,石義雷,龔紅明,等.高超聲速氣動(dòng)熱預(yù)測(cè)技術(shù)及發(fā)展趨勢(shì)[J].航空學(xué)報(bào), 2015, 36(1): 325-345. Peng Zhiyu, Shi Yilei, Gong Hongming, et al. Hypersonic aeroheating prediction technique and its trend of development[J]. Acta Aeronautica et Astronautica Sinica,2015, 36(1):325-345.(in Chinese)

    [12]高翔.迎角下高超聲速?gòu)楊^氣動(dòng)熱和溫度場(chǎng)的計(jì)算與研究[D].南京:南京理工大學(xué),2016. GaoXiang.Calculationofaerodynamicheatingand temperature field for hypersonic warhead under attack angle[D]. Nanjing:Nanjing University of Science&Technology, 2016.(in Chinese)

    [13]Parhizkar H,Karimian S M H. Application of axisymmetric analog to unstructured grid for aeroheating prediction of hypersonic vehicles[J]. International Journal of Numerical Methods for Heat & Fluid Flow,2009,19(3/4):501-520.

    [14]Matting F W. Approximate bridging relations in the transitional regime between continuum and free-molecule flows[J]. Journal of Spacecraft and Rockets,1971,8(1):35-40.

    [15]Nomura S. Correlation of hypersonic stagnation point heat transfer at low Reynolds numbers[J]. AIAA Journal,1983,21(11):1598-1600.

    [16]戎宜生,陳偉芳,石于中,等.過渡流區(qū)圓柱體駐點(diǎn)熱流的工程計(jì)算[J].國(guó)防科技大學(xué)學(xué)報(bào),2008(5):6-9. Rong Yisheng, Chen Weifang, Shi Yuzhong, et al. Engineer calculcation of heat flux at the stagnation point of cylinder in transitional region[J]. Journal of National University of Defense Technology, 2008(5):6-9.(in Chinese)

    [17]戎宜生,劉偉強(qiáng).過渡流區(qū)鈍錐體Linear橋函數(shù)調(diào)節(jié)參數(shù)研究[J].物理學(xué)報(bào), 2012, 61(4): 40508. Rong Yisheng, Liu Weiqiang. Reseach on accommodation parameter of linear bridging relation for blunt cone in transitional regime[J]. Acta Physica Sinica, 2012, 61(4): 40508.(in Chinese)

    [18]Lutz S. Heating correlations for bluff cylinder hypersonic rarefied flows[C]//36th AIAA Thermophysics Conference,2003:4060.

    [19]Swaminathan P K,Taylor J C,Rault D F G,et al. Transition regimeaerodynamicheatingofmissiles[J].Journalof Spacecraft and Rockets,1996,33(5):607-613.

    [20]Engel C D. Miniver upgrade for the avid system,volume 1:lanmin users manual[R]. NASACR-172212,1983.

    [21]Hamilton H H,Greene F A,DeJarnette F R. Approximate method for calculating heating rates on three-dimensional vehicles[J]. Journal of Spacecraft and Rockets,1994,31(3): 345-354.

    [22]陳雄昕,劉衛(wèi)華,羅智勝,等.高超聲速飛行器氣動(dòng)熱研究進(jìn)展[J].航空兵器, 2014(6):8-13. Chen Xiongxin, Liu Weihua, Luo Zhisheng, et al. Research advances of aerodynamic heating for hypersonic aircraft[J]. Aero Weaponry, 2014(6):8-13.(in Chinese)

    [23]Tsien H S. Superaerodynamics,mechanics of rarefied gases[J]. Journal of the Aeronautical Sciences,1946,13(12):653-664.

    [24]Hamilton II H. Approximate method of predicting heating on the windward side of space shuttle orbiter and comparisons with flight data[C]//3rd Joint Thermophysics,F(xiàn)luids,Plasma and Heat Transfer Conference,1982:823.

    [25]歐陽(yáng)水吾,蘇玉宏.高超聲速有迎角鈍頭體三維化學(xué)非平衡黏性激波層流動(dòng)數(shù)值計(jì)算[J].宇航學(xué)報(bào),1992(3):1-12. Ouyang Shuiwu, Su Yuhong. Viscous shock-layer numerical calculations of three dimensional nonequilibrium flows over hypersonic blunt bodies at high angle of attack[J]. Journal of Astronautics, 1992(3):1-12.(in Chinese)

    [26]Lawrence S,Chaussee D,Tannehill J. Application of an upwind algorithm to the three-dimensional parabolized Navier-Stokes equations[C]//8th Computational Fluid Dynamics Conference,1987:1112.

    [27]Cockre C,Auslender A,White J,et al. Aeroheating predictions for the X-43 hyper-X cowl-closed configuration at Mach 7 and 10[C]//40th AIAA Aerospace Sciences Meeting & Exhibit,2002:218.

    [28]賀國(guó)宏,高曉成,龐勇.高超聲速再入體氣動(dòng)熱環(huán)境數(shù)值模擬研究[C]//第十屆全國(guó)計(jì)算流體力學(xué)會(huì)議, 2000. He Guohong, Gao Xiaocheng, Pang Yong. Numerical simulation research on aerodynamic thermal environment of hypersonic reentry body[C]//The 10th National Conference on Fluid Mechanics, 2000.(in Chinese)

    [29]李君哲,閻超,柯倫,等.氣動(dòng)熱CFD計(jì)算的格式效應(yīng)研究[J].北京航空航天大學(xué)學(xué)報(bào),2003(11):1022-1025. Li Junzhe, Yan Chao, Ke Lun, et al. Research on scheme effect of computational fluid dynamics in aerothermal[J]. Journal of Beijing University of Aeronautics and Astronautics, 2003(11): 1022-1025.(in Chinese)

    [30]Lee J,Rho O H. Accuracy of AUSM+ scheme in hypersonic blunt body flow calculations[C]//8th AIAA International Space Planes and Hypersonic Systems and Technologies Conference,1998:1538.

    [31]張翔,閻超,楊威,等.高超聲速飛行器氣動(dòng)熱網(wǎng)格依賴性研究[J].戰(zhàn)術(shù)導(dǎo)彈技術(shù), 2016(3): 21-27. Zhang Xiang, Yan Chao, Yang Wei, et al. Investigation of the grid-dependency in heat transfer simulation for hypersonic vehicle[J]. Tactical Missile Technology, 2016(3): 21-27.(in Chinese)

    [32]Gao Z,Jiang C,Lee C. Aeroheating study of hypersonic chemical nonequilibrium flows around a reentry blunt body[C]//AIAA Space 2014 Conference and Exposition,2014:4415.

    [33]Scott J,Niu Y Y. Comparison of limiters in flux-split algorithms for Euler equations[C]//31st Aerospace Sciences Meeting,1993:68.

    [34]楊建龍,劉猛.限制器對(duì)高超聲速體表面熱流數(shù)值模擬的影響[J].北京航空航天大學(xué)學(xué)報(bào),2014,40(3):417-421. Yang Jianlong, Liu Meng. Influence of limiters on numerical simulation of heating distributions for hypersonic bodies[J]. Journal of Beijing University of Aeronautics and Astronautics, 2014, 40(3):417-421.(in Chinese)

    [35]Zhang S,Chen F,Liu H. Assessment of limiting processes of numerical schemes on hypersonic aeroheating predictions[J]. Journal of Thermophysics and Heat Transfer,2016,30(4):754-769.

    [36]趙雅甜,閻超,孫迪,等.新型三階TVD限制器性能分析[J].北京航空航天大學(xué)學(xué)報(bào), 2017, 43(4): 800. Zhao Yatian, Yan Chao, Sun Di, et al. Performance analysis of a new-type third-order TVD limiter[J]. Journal of Beijing University of Aeronautics and Astronautics, 2017, 43(4): 800.(in Chinese)

    [37]Hollis B,Nowak R,Thompson R,et al. X-33 aerodynamic and aeroheating computations for wind tunnel and flight conditions[C]//24th Atmospheric Flight Mechanics Conference,1999:4163.

    [38]Berry S A,Horvath T J,Weilmuenster K J,et al. X-38 experimental aeroheating at Mach 10[J]. Journal of Spacecraft and Rockets,2004,41(2):293-301.

    [39]Gnoffo P A. An upwind-biased,point-implicit relaxation algorithm for viscous,compressible perfect-gas flows[R]. NASA-TP-2953,1990.

    [40]Aerosoft.GASP version 3The general aerodynamic simulation computational flow analysis software for the scientist and engineer. Users Manual[S].Aerosoft,1996.

    [41]William L K,Peter A G,Stephen J A. Computational aeroheating predictions for X-34[R].AIAA98-0879,1998.

    [42]謝福壽,雷剛,徐元元,等.自由分子流區(qū)氣體外掠平板/圓柱體流動(dòng)與傳熱預(yù)測(cè)[J].西安交通大學(xué)學(xué)報(bào), 2019, 53(3): 6-12. Xie Fushou, Lei Gang, Xu Yuanyuan, et al. Prediction on the flow and heat transfer performance of rarefied gas over objects in free molecular flow regime[J]. Journal of Xian Jiaotong University, 2019, 53(3): 6-12.(in Chinese)

    [43]Lofthouse A J,Scalabrin L C,Boyd I D. Velocity slip and temperaturejumpinhypersonicaerothermodynamics[J]. Journal of thermophysics and heat transfer,2008,22(1):38-49.

    [44]Bird G A. Molecular gas dynamics[R]. NASA STI/Recon Technical ReportA,1976.

    [45]Prasanth P S,Kakkassery J K. Direct simulation Monte Carlo(DSMC):a numerical method for transition-regime flows-a review[J]. Journal of the Indian Institute of Science,2013,86(3):169.

    [46]Liechty D S. Aeroheating analysis for the Mars reconnaissance orbiter with comparison to flight data[J]. Journal of Spacecraft and Rockets,2007,44(6):1226-1231.

    [47]Sengil N,Sengil U. Aerothermal load and drag force analysis of the electromagnetically launched projectiles under rarefied gas conditions[J]. IEEE Transactions on Plasma Science,2015,43(5):1131-1135.

    [48]Marichalar J J,Ostrom C L. Estimating drag and heating coefficients for hollow reentry objects in transitional flow using DSMC[C]// International Orbital Debris Conference,2019:6019.

    [49]Lian Y Y,Wu J S,Cheng G,et al. Development of a parallel hybrid method for the DSMC and NS solver[C]//43rd AIAA Aerospace Sciences Meeting and Exhibit,2005:435.

    [50]Xu X,Wang X,Zhang M,et al. A parallelized hybrid NS/ DSMC-IP approach based on adaptive structured/unstructured overlapping grids for hypersonic transitional flows[J]. Journal of Computational Physics,2018,371:409-433.

    [51]唐振宇,蔡國(guó)飆.兩種用于NS-DSMC耦合方法的連續(xù)失效參數(shù)對(duì)比[J].北京航空航天大學(xué)學(xué)報(bào),2014, 40(3):389-393. Tang Zhenyu, Cai Guobiao. Comparison analysis between two continuum breakdown parameters applied in coupled NS-DSMC method[J]. Journal of Beijing University of Aeronautics and Astronautics, 2014, 40(3):389-393.(in Chinese)

    [52]午辛暄,陳浩,張斌,等.基于DSMC和NS的稀薄高超聲速平板流場(chǎng)分析[J].計(jì)算機(jī)仿真, 2019 (2): 10. Wu Xinxuan, Chen Hao, Zhang Bin, et al. Analysis of rarefied hypersonic flow around flat plate based on the DSMC and NS method[J]. Computer Simulation, 2019 (2): 10.(in Chinese)

    [53]Lofthouse A J,Boyd I D,Wright M J. Effects of continuum breakdown on hypersonic aerothermodynamics[J]. Physics of Fluids,2007,19(2):027105.

    [54]Zhong X. On numerical solutions of Burnett equations for hypersonic flow past 2D circular blunt leading edges in continuumtransitionregime[C]//23rdFluidDynamics,Plasmadynamics,and Lasers Conference,1993:3092.

    [55]Xiao H,He Q J. Aero-heating in hypersonic continuum and rarefied gas flows[J]. Aerospace Science and Technology,2018,82:566-574.

    (責(zé)任編輯王昕)

    作者簡(jiǎn)介

    喻成璋(1997-)男,碩士研究生。主要研究方向:飛行器熱防護(hù)技術(shù)。

    Tel:15950538812

    E-mail:yucz-1248@nuaa.edu.cn

    劉衛(wèi)華(1965-)男,教授,博士生導(dǎo)師。主要研究方向:飛行器環(huán)境與生命保障工程、飛行器燃油系統(tǒng)研究。

    Research Status of Aeroheating Prediction Technology For Hypersionic Aircraft

    Yu Chengzhang*,Liu Weihua

    Key Laboratory of Aircraft Environmental Control and Life Support Industry and Information Technology,Nanjing University of Aeronautics and Astronautics,Nanjing 210016,China

    Abstract: Aeroheating prediction technology is one of the key problems which restrict the development of hypersonic aircraft. When the aircraft is flying at high speed, aerodynamic heating affects the structure intensely, and even cause structural damage in some severe cases. Therefore, in order to ensure the flight safety of aircraft, the thermal protection system must be adopted, and to master the changing rules of aeroheating prediction is the basis for it, which has important guiding significance for aircraft structure design and material selection. From the aspects of experiments, engineering computation methods and numerical simulations, this paper aims to systematically summarize the research results, and explore future development in order to provide useful references for domestic hypersonic aircraft development.

    Key Words: aeroheating; hypersonic aircraft; experimental research; engineering computation method; numeral simulation

    猜你喜歡
    數(shù)值計(jì)算試驗(yàn)研究
    淺談MATLAB在數(shù)學(xué)建模中的應(yīng)用
    卷宗(2016年10期)2017-01-21 16:30:04
    單電機(jī)調(diào)節(jié)無級(jí)變速器試驗(yàn)研究
    加筋土技術(shù)的研究現(xiàn)狀與展望
    彎槽段冰塞形成臨界條件的試驗(yàn)研究
    應(yīng)力加載路徑對(duì)風(fēng)積沙的臨界動(dòng)應(yīng)力影響研究
    矩形迷宮泵內(nèi)部流場(chǎng)數(shù)值模擬及試驗(yàn)研究
    控制棒驅(qū)動(dòng)機(jī)構(gòu)傳動(dòng)件耐磨性研究
    科技視界(2016年12期)2016-05-25 00:40:52
    調(diào)速器導(dǎo)葉開啟時(shí)間對(duì)水電站過渡過程的影響
    南水北調(diào)中線工程輝縣段水泥改性膨脹土工程特性的試驗(yàn)研究
    科技視界(2016年1期)2016-03-30 11:55:29
    平衡流量計(jì)流動(dòng)特性數(shù)值計(jì)算分析
    科技視界(2015年25期)2015-09-01 17:51:38
    人人妻人人澡欧美一区二区| 国产精品一区二区三区四区久久| 国内揄拍国产精品人妻在线| 精品免费久久久久久久清纯| 热99re8久久精品国产| 国产 一区 欧美 日韩| 国产精品久久久久久久久免 | 国产主播在线观看一区二区| 中文在线观看免费www的网站| 真人一进一出gif抽搐免费| 国产精品三级大全| 欧美av亚洲av综合av国产av| 午夜福利高清视频| 中文字幕av成人在线电影| 身体一侧抽搐| 国产精品99久久久久久久久| 亚洲午夜理论影院| 久久精品亚洲精品国产色婷小说| 激情在线观看视频在线高清| 给我免费播放毛片高清在线观看| 欧美黑人巨大hd| 听说在线观看完整版免费高清| 少妇丰满av| 国产99白浆流出| 亚洲中文字幕日韩| 黄色成人免费大全| 好看av亚洲va欧美ⅴa在| 亚洲国产日韩欧美精品在线观看 | 最新中文字幕久久久久| av欧美777| 美女免费视频网站| 久久99热这里只有精品18| 少妇裸体淫交视频免费看高清| 亚洲人与动物交配视频| 国产av麻豆久久久久久久| 成年女人看的毛片在线观看| 国产主播在线观看一区二区| 性色avwww在线观看| 国产成人啪精品午夜网站| 又黄又爽又免费观看的视频| 一区福利在线观看| 色精品久久人妻99蜜桃| 丰满人妻熟妇乱又伦精品不卡| 熟妇人妻久久中文字幕3abv| 亚洲精品一区av在线观看| 久久久国产精品麻豆| 最近在线观看免费完整版| 国产色爽女视频免费观看| 婷婷精品国产亚洲av| 欧美在线黄色| 欧美激情在线99| 特大巨黑吊av在线直播| 特级一级黄色大片| 午夜福利欧美成人| 亚洲精华国产精华精| 色av中文字幕| av天堂中文字幕网| 美女高潮喷水抽搐中文字幕| 国产真实乱freesex| 两个人看的免费小视频| 99久久久亚洲精品蜜臀av| 亚洲18禁久久av| 一本一本综合久久| 国产激情偷乱视频一区二区| 久久欧美精品欧美久久欧美| or卡值多少钱| 免费在线观看亚洲国产| 国产成+人综合+亚洲专区| 亚洲美女黄片视频| а√天堂www在线а√下载| 首页视频小说图片口味搜索| 两性午夜刺激爽爽歪歪视频在线观看| 欧美+日韩+精品| 91字幕亚洲| 亚洲精品国产精品久久久不卡| tocl精华| 人人妻人人看人人澡| 国内少妇人妻偷人精品xxx网站| 亚洲av成人av| x7x7x7水蜜桃| 一本综合久久免费| 国产精品av视频在线免费观看| 狠狠狠狠99中文字幕| 国产成人欧美在线观看| 免费看光身美女| 国产精品久久久久久亚洲av鲁大| 国内毛片毛片毛片毛片毛片| 精品午夜福利视频在线观看一区| 网址你懂的国产日韩在线| 久久亚洲精品不卡| 禁无遮挡网站| 色综合婷婷激情| ponron亚洲| 成人特级av手机在线观看| 18禁黄网站禁片午夜丰满| 亚洲中文日韩欧美视频| 色在线成人网| 亚洲精华国产精华精| 尤物成人国产欧美一区二区三区| 国产午夜福利久久久久久| 黄片小视频在线播放| 免费观看的影片在线观看| 天堂网av新在线| 国产免费一级a男人的天堂| 少妇的逼好多水| 欧美日韩一级在线毛片| 一区二区三区国产精品乱码| 久久久久亚洲av毛片大全| 亚洲狠狠婷婷综合久久图片| 国产免费男女视频| 久久久久久久久大av| 亚洲av免费在线观看| 91久久精品国产一区二区成人 | 亚洲av成人不卡在线观看播放网| 51午夜福利影视在线观看| 午夜激情福利司机影院| 欧美激情久久久久久爽电影| 国产欧美日韩一区二区三| 国产精品久久久久久人妻精品电影| 一卡2卡三卡四卡精品乱码亚洲| 熟妇人妻久久中文字幕3abv| 九九久久精品国产亚洲av麻豆| 久久亚洲真实| 黑人欧美特级aaaaaa片| 中文字幕av成人在线电影| 他把我摸到了高潮在线观看| 国产精华一区二区三区| 国产在线精品亚洲第一网站| 大型黄色视频在线免费观看| 免费观看人在逋| av片东京热男人的天堂| 岛国在线观看网站| 黄片小视频在线播放| 亚洲五月天丁香| 美女高潮的动态| 精品人妻1区二区| 村上凉子中文字幕在线| 18禁黄网站禁片免费观看直播| 国产一区在线观看成人免费| 国产真实伦视频高清在线观看 | 中出人妻视频一区二区| 日韩欧美在线二视频| 国内久久婷婷六月综合欲色啪| 国产欧美日韩精品一区二区| 亚洲无线观看免费| 国产精品1区2区在线观看.| 亚洲中文字幕日韩| 亚洲成人久久性| 男女下面进入的视频免费午夜| 乱人视频在线观看| 天天添夜夜摸| 狠狠狠狠99中文字幕| 夜夜爽天天搞| 91久久精品国产一区二区成人 | 亚洲avbb在线观看| 观看免费一级毛片| a级一级毛片免费在线观看| 一个人看视频在线观看www免费 | 91久久精品国产一区二区成人 | 国产黄片美女视频| 午夜免费观看网址| 午夜精品久久久久久毛片777| 成人精品一区二区免费| 精品国产三级普通话版| 在线观看免费午夜福利视频| 无限看片的www在线观看| 久久99热这里只有精品18| 日韩亚洲欧美综合| 一个人看视频在线观看www免费 | 在线播放无遮挡| 精华霜和精华液先用哪个| 欧美成人a在线观看| 观看免费一级毛片| 日本五十路高清| 欧美zozozo另类| 观看免费一级毛片| 日本与韩国留学比较| 男女视频在线观看网站免费| 国产探花在线观看一区二区| 搡老岳熟女国产| 18禁裸乳无遮挡免费网站照片| 国产精品电影一区二区三区| 人人妻人人澡欧美一区二区| 99riav亚洲国产免费| 日韩大尺度精品在线看网址| 欧美在线黄色| 午夜精品久久久久久毛片777| 9191精品国产免费久久| 午夜a级毛片| a在线观看视频网站| 99精品在免费线老司机午夜| 亚洲精品成人久久久久久| 久久天躁狠狠躁夜夜2o2o| 夜夜躁狠狠躁天天躁| 91字幕亚洲| 国产精品自产拍在线观看55亚洲| 性欧美人与动物交配| 亚洲av二区三区四区| 久久精品影院6| av片东京热男人的天堂| 动漫黄色视频在线观看| 精品久久久久久久毛片微露脸| 国产精品亚洲av一区麻豆| 国产精品99久久99久久久不卡| 免费电影在线观看免费观看| 夜夜夜夜夜久久久久| 黄色日韩在线| 久久久久久国产a免费观看| 久久草成人影院| 变态另类成人亚洲欧美熟女| 偷拍熟女少妇极品色| 欧美国产日韩亚洲一区| 国产精品久久久久久亚洲av鲁大| 嫩草影院入口| 香蕉久久夜色| 国产黄a三级三级三级人| 欧美日韩亚洲国产一区二区在线观看| 亚洲av成人不卡在线观看播放网| 男人和女人高潮做爰伦理| 欧美日韩国产亚洲二区| АⅤ资源中文在线天堂| www.熟女人妻精品国产| 国产真实伦视频高清在线观看 | 51午夜福利影视在线观看| 国产高清三级在线| 男女下面进入的视频免费午夜| 女警被强在线播放| 岛国在线观看网站| 中文字幕熟女人妻在线| 中文字幕久久专区| 国产91精品成人一区二区三区| 午夜激情欧美在线| 一区二区三区激情视频| 一夜夜www| 亚洲精品乱码久久久v下载方式 | а√天堂www在线а√下载| 国产亚洲av嫩草精品影院| 99久久精品一区二区三区| 国产真人三级小视频在线观看| 一本一本综合久久| 男女那种视频在线观看| 国产精品影院久久| 最近在线观看免费完整版| 激情在线观看视频在线高清| 露出奶头的视频| 亚洲精品乱码久久久v下载方式 | 色尼玛亚洲综合影院| 久久久久久九九精品二区国产| 欧美在线一区亚洲| 免费一级毛片在线播放高清视频| 色吧在线观看| 亚洲人成网站在线播| 国产爱豆传媒在线观看| 高清毛片免费观看视频网站| 婷婷精品国产亚洲av| 亚洲av电影在线进入| 亚洲第一电影网av| av中文乱码字幕在线| 午夜免费激情av| 两性午夜刺激爽爽歪歪视频在线观看| 琪琪午夜伦伦电影理论片6080| 亚洲av电影不卡..在线观看| 无遮挡黄片免费观看| 一本一本综合久久| 麻豆成人av在线观看| 日韩欧美免费精品| a级一级毛片免费在线观看| 久久婷婷人人爽人人干人人爱| 51国产日韩欧美| 欧美最新免费一区二区三区 | 色播亚洲综合网| 国产伦精品一区二区三区四那| 热99在线观看视频| 成人鲁丝片一二三区免费| 精品一区二区三区视频在线观看免费| 精品免费久久久久久久清纯| 一本综合久久免费| 成人特级av手机在线观看| 一本精品99久久精品77| 99久久久亚洲精品蜜臀av| 精品乱码久久久久久99久播| 九色国产91popny在线| 欧美在线一区亚洲| 特大巨黑吊av在线直播| 757午夜福利合集在线观看| 免费在线观看影片大全网站| 精品人妻1区二区| 国产伦精品一区二区三区视频9 | 精品免费久久久久久久清纯| e午夜精品久久久久久久| 国产激情欧美一区二区| 精品国产超薄肉色丝袜足j| 午夜福利欧美成人| 男女视频在线观看网站免费| 亚洲成a人片在线一区二区| 一级毛片高清免费大全| 99视频精品全部免费 在线| 日韩 欧美 亚洲 中文字幕| 国产亚洲欧美在线一区二区| 操出白浆在线播放| 国产精品一区二区三区四区免费观看 | 看片在线看免费视频| 少妇丰满av| 一进一出抽搐动态| 欧美精品啪啪一区二区三区| 久久久久久国产a免费观看| 在线播放国产精品三级| avwww免费| 国产97色在线日韩免费| 悠悠久久av| 午夜福利在线观看吧| 天天一区二区日本电影三级| 国产精品女同一区二区软件 | 亚洲av不卡在线观看| 国产三级中文精品| 一进一出抽搐动态| 国产乱人视频| 日韩亚洲欧美综合| 精品国产超薄肉色丝袜足j| 亚洲精品粉嫩美女一区| 尤物成人国产欧美一区二区三区| 在线观看一区二区三区| 99久久综合精品五月天人人| 美女大奶头视频| 亚洲欧美日韩高清专用| 香蕉久久夜色| 色噜噜av男人的天堂激情| 欧美一区二区亚洲| 给我免费播放毛片高清在线观看| 亚洲第一电影网av| 精品久久久久久,| 99热6这里只有精品| 九九在线视频观看精品| 国产精品久久视频播放| 免费无遮挡裸体视频| 欧美日韩乱码在线| 特级一级黄色大片| 国产伦人伦偷精品视频| 午夜免费男女啪啪视频观看 | 香蕉久久夜色| 色综合站精品国产| 久久九九热精品免费| 波野结衣二区三区在线 | 日韩大尺度精品在线看网址| 99热只有精品国产| 一本精品99久久精品77| 午夜免费观看网址| 国产综合懂色| 天堂动漫精品| 成人欧美大片| 亚洲av电影在线进入| 看片在线看免费视频| 一级a爱片免费观看的视频| aaaaa片日本免费| 国内精品美女久久久久久| 午夜老司机福利剧场| 激情在线观看视频在线高清| 色播亚洲综合网| 日韩国内少妇激情av| 久久久久久久午夜电影| 看片在线看免费视频| 久久精品影院6| 亚洲人成网站高清观看| 岛国视频午夜一区免费看| 制服人妻中文乱码| 久久婷婷人人爽人人干人人爱| 欧美极品一区二区三区四区| 丁香欧美五月| 色av中文字幕| 国产一区在线观看成人免费| 日本三级黄在线观看| 亚洲精品成人久久久久久| 婷婷丁香在线五月| 亚洲成人中文字幕在线播放| 亚洲av不卡在线观看| 久久精品国产自在天天线| 免费看日本二区| 午夜福利18| 又爽又黄无遮挡网站| 村上凉子中文字幕在线| 亚洲欧美精品综合久久99| 18禁黄网站禁片免费观看直播| 中出人妻视频一区二区| eeuss影院久久| 婷婷丁香在线五月| 五月伊人婷婷丁香| 成熟少妇高潮喷水视频| 男女下面进入的视频免费午夜| 老司机福利观看| 亚洲成人中文字幕在线播放| 香蕉丝袜av| 看黄色毛片网站| 90打野战视频偷拍视频| 国产免费男女视频| 久久亚洲真实| 动漫黄色视频在线观看| 亚洲成人久久爱视频| 日韩欧美国产一区二区入口| 国产97色在线日韩免费| 18禁黄网站禁片免费观看直播| 亚洲avbb在线观看| 精品国产美女av久久久久小说| 久久伊人香网站| 欧美午夜高清在线| 久久国产精品人妻蜜桃| 亚洲国产欧洲综合997久久,| 看片在线看免费视频| 久久久精品欧美日韩精品| netflix在线观看网站| 中文字幕熟女人妻在线| 老司机深夜福利视频在线观看| 成人午夜高清在线视频| 国产黄片美女视频| 禁无遮挡网站| 国产爱豆传媒在线观看| 日本撒尿小便嘘嘘汇集6| 伊人久久大香线蕉亚洲五| 热99re8久久精品国产| 日韩欧美 国产精品| 非洲黑人性xxxx精品又粗又长| 97碰自拍视频| 亚洲欧美激情综合另类| 一级a爱片免费观看的视频| 老司机福利观看| 午夜福利免费观看在线| 人人妻人人看人人澡| 黄色视频,在线免费观看| 搞女人的毛片| 嫩草影院入口| 亚洲国产高清在线一区二区三| 久久九九热精品免费| 我的老师免费观看完整版| 欧美日韩一级在线毛片| 午夜影院日韩av| 中出人妻视频一区二区| www.色视频.com| 制服人妻中文乱码| 国产黄a三级三级三级人| 国产亚洲精品综合一区在线观看| 观看美女的网站| 亚洲av五月六月丁香网| 淫秽高清视频在线观看| 人妻丰满熟妇av一区二区三区| 婷婷六月久久综合丁香| 国产视频内射| 在线十欧美十亚洲十日本专区| 波多野结衣高清作品| 国产亚洲欧美98| 日本黄色视频三级网站网址| 淫秽高清视频在线观看| 最近最新中文字幕大全电影3| 97人妻精品一区二区三区麻豆| 91久久精品国产一区二区成人 | 又粗又爽又猛毛片免费看| 成人无遮挡网站| 一级毛片女人18水好多| 精品乱码久久久久久99久播| 蜜桃久久精品国产亚洲av| 亚洲乱码一区二区免费版| 国产美女午夜福利| 午夜福利视频1000在线观看| 亚洲最大成人中文| 嫩草影院入口| 国产免费一级a男人的天堂| 无人区码免费观看不卡| 亚洲成人中文字幕在线播放| 人妻久久中文字幕网| 国产主播在线观看一区二区| 国产色婷婷99| 淫秽高清视频在线观看| 色老头精品视频在线观看| 中亚洲国语对白在线视频| 久久国产精品影院| 日韩有码中文字幕| 美女大奶头视频| 国产在视频线在精品| 国产精品,欧美在线| 日韩大尺度精品在线看网址| 偷拍熟女少妇极品色| 白带黄色成豆腐渣| 国产亚洲av嫩草精品影院| 精品电影一区二区在线| 成人鲁丝片一二三区免费| 精品国产三级普通话版| 日韩大尺度精品在线看网址| 叶爱在线成人免费视频播放| 小说图片视频综合网站| 久久久国产成人免费| 亚洲成人精品中文字幕电影| 欧洲精品卡2卡3卡4卡5卡区| 久久精品影院6| 亚洲无线观看免费| 久久久久国内视频| 老熟妇仑乱视频hdxx| 欧美3d第一页| 久久婷婷人人爽人人干人人爱| 夜夜躁狠狠躁天天躁| 国产精品自产拍在线观看55亚洲| 久久久久免费精品人妻一区二区| 男女视频在线观看网站免费| 国产爱豆传媒在线观看| 欧美一级毛片孕妇| 88av欧美| 亚洲国产欧美网| 性色av乱码一区二区三区2| 天天一区二区日本电影三级| 舔av片在线| 亚洲成a人片在线一区二区| 欧美日韩亚洲国产一区二区在线观看| 国产乱人伦免费视频| 国产精品日韩av在线免费观看| 国产av一区在线观看免费| 三级国产精品欧美在线观看| 欧美黑人巨大hd| 国产毛片a区久久久久| 51午夜福利影视在线观看| 一个人看视频在线观看www免费 | av黄色大香蕉| 国产欧美日韩一区二区三| 精品人妻一区二区三区麻豆 | 九九久久精品国产亚洲av麻豆| 高清在线国产一区| 草草在线视频免费看| 国产精品久久电影中文字幕| 国产主播在线观看一区二区| 亚洲国产中文字幕在线视频| 国产一区二区在线观看日韩 | 国产精品一区二区三区四区免费观看 | 欧美黑人欧美精品刺激| 国产色婷婷99| 美女黄网站色视频| 一级毛片女人18水好多| 午夜福利在线在线| 美女高潮的动态| 麻豆久久精品国产亚洲av| 99riav亚洲国产免费| 午夜福利在线在线| 99热这里只有精品一区| 久久精品综合一区二区三区| 香蕉久久夜色| 国产成人影院久久av| 久久久久亚洲av毛片大全| 国产97色在线日韩免费| 国产成年人精品一区二区| 91在线观看av| 91久久精品国产一区二区成人 | 久久性视频一级片| 国产野战对白在线观看| xxx96com| 亚洲熟妇中文字幕五十中出| 在线十欧美十亚洲十日本专区| 欧美成狂野欧美在线观看| 日韩成人在线观看一区二区三区| 日本三级黄在线观看| 亚洲av日韩精品久久久久久密| 国产精品乱码一区二三区的特点| 亚洲av第一区精品v没综合| 国内精品美女久久久久久| 免费看日本二区| 国产国拍精品亚洲av在线观看 | 内射极品少妇av片p| 亚洲最大成人手机在线| 亚洲精品成人久久久久久| 国内精品久久久久精免费| av视频在线观看入口| 动漫黄色视频在线观看| 免费看日本二区| 波多野结衣高清作品| 成人国产综合亚洲| 欧美色视频一区免费| 国产真人三级小视频在线观看| 一级毛片女人18水好多| 日韩av在线大香蕉| 一进一出抽搐gif免费好疼| 老司机福利观看| 久久精品国产清高在天天线| 国产伦一二天堂av在线观看| 国产精华一区二区三区| 五月伊人婷婷丁香| 国产精品野战在线观看| 91久久精品国产一区二区成人 | 久久午夜亚洲精品久久| 欧美在线黄色| 国产激情偷乱视频一区二区| 欧美zozozo另类| 精品人妻偷拍中文字幕| 日韩欧美在线二视频| 亚洲av电影在线进入| 九九在线视频观看精品| www国产在线视频色| 国产av不卡久久| h日本视频在线播放| 51国产日韩欧美| 国产私拍福利视频在线观看| 中文资源天堂在线| 国产欧美日韩一区二区精品| 久久精品夜夜夜夜夜久久蜜豆| 亚洲成av人片在线播放无| 黄片大片在线免费观看| 啦啦啦免费观看视频1| 国产伦精品一区二区三区四那| 最近最新免费中文字幕在线| 亚洲美女黄片视频| 黑人欧美特级aaaaaa片| 色播亚洲综合网| 日本免费一区二区三区高清不卡| svipshipincom国产片| 亚洲av电影在线进入| 亚洲 欧美 日韩 在线 免费| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利免费观看在线| 国产精品久久久久久久电影 | 深夜精品福利| 久久久久国产精品人妻aⅴ院| 久久午夜亚洲精品久久|