王 浩，徐晶磊

(1.北京航空航天大學(xué) 能源與動(dòng)力工程學(xué)院，北京 100191；2.上海交通大學(xué) 航空航天學(xué)院，上海 200240；3.北京航空航天大學(xué)杭州創(chuàng)新研究院（余杭）航空學(xué)院，杭州 310023)

0 引言

高超聲速轉(zhuǎn)捩對(duì)高超聲速飛行器的飛行狀態(tài)有顯著影響，主要體現(xiàn)在氣動(dòng)熱/力與動(dòng)力方面。一方面，轉(zhuǎn)捩后湍流的壁面摩擦力和傳熱速率遠(yuǎn)高于層流，大約是3～5 倍[1]，使得高超聲速飛行器熱防護(hù)系統(tǒng)的設(shè)計(jì)非常具有挑戰(zhàn)性；轉(zhuǎn)捩后的氣動(dòng)載荷分布也會(huì)發(fā)生顯著變化，影響氣動(dòng)配平、帶來(lái)飛行安全問(wèn)題。另一方面，湍流邊界層的抗分離能力更強(qiáng)，可以有效地克服超燃沖壓發(fā)動(dòng)機(jī)入口處的邊界層分離，并且湍流邊界層的摻混能力也更強(qiáng)，可以增強(qiáng)發(fā)動(dòng)機(jī)中燃料與空氣的摻混，使燃料可以充分燃燒。如果能夠準(zhǔn)確預(yù)測(cè)到轉(zhuǎn)捩的位置和區(qū)域，就更有把握通過(guò)優(yōu)化外形、改變壁面條件等多種手段來(lái)延遲或促進(jìn)轉(zhuǎn)捩的發(fā)生，以滿足種種工程需求。在真實(shí)飛行器的三維邊界層中，受壓力梯度和后掠角的共同作用，邊界層內(nèi)近壁區(qū)往往會(huì)出現(xiàn)與邊界層外勢(shì)流平面內(nèi)流線方向垂直的流動(dòng)分量，稱為橫流流動(dòng)[2]。當(dāng)橫流速度出現(xiàn)拐點(diǎn)時(shí)，將會(huì)產(chǎn)生流動(dòng)失穩(wěn)并容易發(fā)展為三維邊界層轉(zhuǎn)捩。因此，高超聲橫流轉(zhuǎn)捩的研究有著迫切的實(shí)際需求，發(fā)展相應(yīng)的預(yù)測(cè)或計(jì)算方法具有重要的工程意義。

高超聲速橫流轉(zhuǎn)捩是集來(lái)流湍流-激波相互作用、可壓縮湍流、橫流、轉(zhuǎn)捩等復(fù)雜流動(dòng)因素于一身的、長(zhǎng)期困擾科研人員的難題。轉(zhuǎn)捩本身的難點(diǎn)在于，其起始位置對(duì)擾動(dòng)的演化十分敏感，只有精確計(jì)算出擾動(dòng)的演化過(guò)程才能得到準(zhǔn)確的起始位置。而演化過(guò)程涉及到上述4 個(gè)問(wèn)題的前3 個(gè)，即使這3 個(gè)問(wèn)題得到準(zhǔn)確建模，建模公式的同時(shí)使用也可能產(chǎn)生耦合誤差。目前在工程應(yīng)用層面，基于雷諾平均N-S方程（RANS）的湍流/轉(zhuǎn)捩模型是預(yù)測(cè)轉(zhuǎn)捩較為實(shí)際和有效的選擇。2000 年Suzen 和Huang[3]構(gòu)造出間歇因子輸運(yùn)方程，它與湍流輸運(yùn)方程合并成為轉(zhuǎn)捩模型。2006 年Menter 等[4-5]構(gòu)造出γ-Reθ-SST 四方程轉(zhuǎn)捩模型，只采用當(dāng)?shù)亓鲃?dòng)變量，從而與現(xiàn)代CFD 程序兼容。隨后該模型集成在ANSYS 商業(yè)軟件，達(dá)到工業(yè)級(jí)成熟度，成為轉(zhuǎn)捩模型的標(biāo)桿。該模型沒(méi)有考慮轉(zhuǎn)捩機(jī)理，而是基于經(jīng)驗(yàn)關(guān)系式構(gòu)造轉(zhuǎn)捩判據(jù)，而經(jīng)驗(yàn)關(guān)系式是通過(guò)低速湍流邊界層、壓氣機(jī)等標(biāo)定的。2009 年王亮和符松[6]提出k-ω-γ轉(zhuǎn)捩模型，對(duì)超聲速平板邊界層和超/高超聲速圓錐邊界層的預(yù)測(cè)結(jié)果與試驗(yàn)結(jié)果定量符合。直到2016 年，周玲在高超聲速轉(zhuǎn)捩模型中引入了橫流模態(tài)時(shí)間尺度和橫流轉(zhuǎn)捩準(zhǔn)則，對(duì)原模型預(yù)測(cè)的HIFiRE-5 橫流轉(zhuǎn)捩圖像有所改善，是高超聲速橫流轉(zhuǎn)捩的先驅(qū)建模工作。2016 年張毅鋒等[7]對(duì)γ-Reθ-SST 轉(zhuǎn)捩模型的低速經(jīng)驗(yàn)關(guān)系式進(jìn)行間歇因子的馬赫數(shù)修正和湍流普朗特?cái)?shù)修正，顯著改進(jìn)了原模型的性能，對(duì)多個(gè)高超聲速尖錐的計(jì)算結(jié)果與試驗(yàn)符合較好。2017 年張毅鋒等[8]進(jìn)一步引入以橫流雷諾數(shù)為基礎(chǔ)的轉(zhuǎn)捩準(zhǔn)則，對(duì)HIFiRE-5 橫流轉(zhuǎn)捩標(biāo)模，在不同雷諾數(shù)下都取得了與噪聲風(fēng)洞實(shí)驗(yàn)非常吻合的計(jì)算結(jié)果。2021 年向星皓等[9]使用eN方法對(duì)高超聲速橫流轉(zhuǎn)捩數(shù)據(jù)進(jìn)行拓展，結(jié)合橫流強(qiáng)度與表面粗糙度構(gòu)造當(dāng)?shù)鼗母叱曀贆M流轉(zhuǎn)捩判據(jù)，構(gòu)造的新模型對(duì)多狀態(tài)下的高超聲速尖錐進(jìn)行橫流轉(zhuǎn)捩預(yù)測(cè)，取得了與實(shí)驗(yàn)結(jié)果符合較好的預(yù)測(cè)效果。

從上述間歇因子轉(zhuǎn)捩模型的進(jìn)化路線來(lái)看，為了捕捉高超聲速橫流轉(zhuǎn)捩，在保證常規(guī)轉(zhuǎn)捩的預(yù)測(cè)精度的同時(shí)，還需引入高超聲速（或可壓縮）效應(yīng)和橫流判據(jù)等，建模難度和工作量很大?；赗ANS 的“轉(zhuǎn)捩模型”還有一條古老的路線，那就是直接用湍流方程（k、ω等）計(jì)算轉(zhuǎn)捩。2019 年徐晶磊等將KDO 模型[10-11]推廣至轉(zhuǎn)捩的計(jì)算[12-13]，具體做法是，以輸運(yùn)變量r=μt/μ重新標(biāo)定模型參數(shù)，使模型整體達(dá)到流動(dòng)結(jié)構(gòu)自適應(yīng)。而模型捕捉轉(zhuǎn)捩的機(jī)制是，通過(guò)流動(dòng)結(jié)構(gòu)的自適應(yīng)和湍流輸運(yùn)特性的保存，自發(fā)演化出轉(zhuǎn)捩。這樣，不用引入特定的轉(zhuǎn)捩機(jī)理，卻有可能捕捉到特定類型的轉(zhuǎn)捩，具有潛在的多機(jī)理普適性。目前，該模型已成功預(yù)測(cè)了T3A 與T3B 平板邊界層旁路轉(zhuǎn)捩、T3A-平板邊界層自然轉(zhuǎn)捩、Aero-A 翼型分離泡轉(zhuǎn)捩、DLR-F5 橫流轉(zhuǎn)捩、6 ∶ 1 橢球橫流轉(zhuǎn)捩、超聲速平板邊界層轉(zhuǎn)捩以及超聲速尖錐轉(zhuǎn)捩。本文將該模型推廣至高超聲速橫流轉(zhuǎn)捩的預(yù)測(cè)，采用經(jīng)典標(biāo)模HIFiRE-5，評(píng)估其捕捉高超聲速橫流轉(zhuǎn)捩的潛力與不足，為高超聲速橫流轉(zhuǎn)捩模型的發(fā)展提供參考。

1 CKDO 轉(zhuǎn)捩模型

1.1 CKDO 湍流模型

精確的可壓縮湍動(dòng)能控制方程[14]為：

對(duì)于可壓縮湍流，可壓縮效應(yīng)的建模是十分必要的。對(duì)于可壓縮修正的KDO 模型，即CKDO 模型，方程右端各項(xiàng)依次分別?；癁椋?/p>

為了維持計(jì)算的穩(wěn)定性，ε2的建模被分為兩層建模，分別黏性子層和外層:

若Rek＜10，

若Rek＞10，

模型中的許多經(jīng)驗(yàn)數(shù)值，只是為了精準(zhǔn)刻畫(huà)從DNS 數(shù)據(jù)得到的r～Rb與Rek～Aε分布曲線而采用的多層結(jié)構(gòu)標(biāo)定方法，導(dǎo)致的多層系數(shù)[12]。本質(zhì)上，模型的經(jīng)驗(yàn)參數(shù)只有2 個(gè)，即雷諾應(yīng)力本構(gòu)關(guān)系的Bradshaw 參數(shù)Rb及ε方程毀滅項(xiàng)系數(shù)Aε，這兩個(gè)經(jīng)驗(yàn)參數(shù)經(jīng)平板邊界層的Reθ=4 060 剖面標(biāo)定后，就再未改變。

1.2 CKDO 轉(zhuǎn)捩模型

注意到，雷諾應(yīng)力本構(gòu)關(guān)系中的Bradshaw 系數(shù)，在公式（9）被標(biāo)定為關(guān)于Rek的函數(shù)，Rek包含壁面距離，而無(wú)論流動(dòng)形態(tài)如何，壁面距離是固定不變的，不反映流動(dòng)結(jié)構(gòu)，這就是為什么原始 KDO 在T3B 平板算例中展現(xiàn)轉(zhuǎn)捩行為，但轉(zhuǎn)捩過(guò)早。KDO 以及CKDO 轉(zhuǎn)捩模型引入湍流黏性與層流黏性之比r，它傳達(dá)湍流的強(qiáng)度，是一個(gè)輸運(yùn)變量。再次使用平板邊界層的Reθ=4 060 截面把Rb標(biāo)定為：

CKDO 和KDO 模型的經(jīng)驗(yàn)參數(shù)標(biāo)定于充分發(fā)展平板湍流邊界層的Reθ=4 060 剖面，但又能用于充分轉(zhuǎn)捩區(qū)域甚至橫流區(qū)域，這一點(diǎn)是令人疑惑的。解釋如下：首先，理論上允許這樣的CFD 計(jì)算系統(tǒng)存在，正如N-S 方程一樣，并不考慮任何的轉(zhuǎn)捩機(jī)理，卻能捕捉轉(zhuǎn)捩，CKDO-tran 計(jì)算轉(zhuǎn)捩的機(jī)制與N-S 方程是一致的。其次，可以注意到，CKDO 和KDO 捕捉轉(zhuǎn)捩的關(guān)鍵在黏性比r=μt/μ＜ 1 的涵蓋層流區(qū)域和層流-湍流轉(zhuǎn)捩區(qū)域（“pre-transition”區(qū)域）。只要r＜1 時(shí)，模型采用的相似性規(guī)律是普適的，那么模型就能普適預(yù)測(cè)轉(zhuǎn)捩。實(shí)際情況是，當(dāng)r＜ 1 時(shí)，公式（14）在Reθ=150～900 的“pre-transition”區(qū)域仍然普適成立[12]。為什么普適?r＜ 1 的區(qū)域常常對(duì)應(yīng)于y+＜20 的非常近壁區(qū)，即黏性子層和部分過(guò)渡層，而該區(qū)域的流動(dòng)屬于壁面主導(dǎo)的流動(dòng)，該區(qū)域的流體微團(tuán)只感受到壁面的影響而基本無(wú)視外界的影響，形成一個(gè)普適的“局部平板邊界層”。將采用公式（14）的CKDO 模型記作CKDO-tran，將采用公式（14）的KDO 模型記作KDO-tran。

2 超聲速轉(zhuǎn)捩算例評(píng)估

為了展示對(duì)于超聲速轉(zhuǎn)捩流動(dòng)，考慮可壓縮效應(yīng)的CKDO-tran 模型相對(duì)于KDO-tran 模型的優(yōu)勢(shì)，本節(jié)使用超聲速平板邊界層的DNS 數(shù)據(jù)作為基準(zhǔn)驗(yàn)證算例。德州大學(xué)的Jiang 等[15]進(jìn)行了來(lái)流馬赫4.5 的平板直接數(shù)值模擬（U∞=705.093 m/s，T∞=61.111 K），壁溫恒定等于層流流態(tài)下的絕熱溫度。文獻(xiàn)用oblique-mode 擾動(dòng)激發(fā)湍流，這意味著來(lái)流的湍流度未知，因此不同的研究者采用不同的來(lái)流湍流度，來(lái)使他們的模型獲得最佳表現(xiàn)。對(duì)于KDO-tran 和CKDO-tran，來(lái)流湍流度分別取1%和0.73%，計(jì)算結(jié)果見(jiàn)圖1。CKDO-tran 預(yù)測(cè)的湍流區(qū)摩阻大小與DNS 結(jié)果更接近，且捕捉到了轉(zhuǎn)捩后的摩阻波動(dòng)；而KDO-tran 預(yù)測(cè)的摩阻偏低且過(guò)于平滑，原因在于未考慮可壓縮效應(yīng)。低估超聲速流動(dòng)的摩阻，這一現(xiàn)象在文獻(xiàn)[16]也有體現(xiàn)，就是因?yàn)楦鬓D(zhuǎn)捩模型搭載的湍流方程忽略了可壓縮效應(yīng)的緣故。然而KDOtran 計(jì)算的轉(zhuǎn)捩過(guò)渡段好于CKDO-tran，這與KDOtran 的一貫表現(xiàn)不一致（KDO-tran 的轉(zhuǎn)捩過(guò)渡段總是很短），推測(cè)是KDO-tran 的湍流黏性過(guò)小，對(duì)過(guò)渡段造成了緩沖，歪打正著。

圖1 平板表面摩阻分布Fig.1 Skin friction distribution on the flat plate

3 高超聲速橫流轉(zhuǎn)捩算例評(píng)估

HIFiRE（Hypersonic International Flight Research Experimentation）是由美國(guó)空軍研究實(shí)驗(yàn)室和澳大利亞防衛(wèi)科研組織聯(lián)合舉辦的高超聲速飛行器測(cè)試項(xiàng)目，其目的是發(fā)展和驗(yàn)證下一代高超聲速系統(tǒng)的關(guān)鍵技術(shù)。其中HIFiRE-5 外形設(shè)計(jì)為長(zhǎng)短軸之比為2 ∶ 1 的橢圓錐，使得三維效應(yīng)和橫向流動(dòng)更加明顯，被專門(mén)用于研究高超聲速三維邊界層轉(zhuǎn)捩。本文計(jì)算采用Juliano 等[17-20]風(fēng)洞試驗(yàn)中的HIFiRE-5 縮比模型，總長(zhǎng)為328 mm，橫截面長(zhǎng)軸半徑為82 mm，短軸半錐角為7°，短軸所在截面頭部半徑為0.95 mm，模型三維示意圖如圖2 所示。HIFiRE-5 模型的計(jì)算網(wǎng)格采用從短軸中心線平面劃分的二分之一模型，計(jì)算采用三重網(wǎng)格并行計(jì)算技術(shù)，共分8 個(gè)網(wǎng)格塊，每個(gè)網(wǎng)格塊的各邊網(wǎng)格點(diǎn)數(shù)均控制為4n+1，經(jīng)過(guò)網(wǎng)格無(wú)關(guān)性驗(yàn)證后，最終模型計(jì)算網(wǎng)格流向共181 個(gè)網(wǎng)格點(diǎn)，法向共113 個(gè)網(wǎng)格點(diǎn)，周向共129 個(gè)網(wǎng)格點(diǎn)，總網(wǎng)格量約為252 萬(wàn)。法向第一層網(wǎng)格高度約為5×10-4mm，并在激波位置附近加密。HIFiRE-5 結(jié)構(gòu)網(wǎng)格三維示意圖如圖3 所示。

圖2 HIFiRE-5 模型三維示意圖Fig.2 3D model of HIFiRE-5

圖3 HIFiRE-5 結(jié)構(gòu)網(wǎng)格示意圖Fig.3 Structured mesh of HIFiRE-5

本節(jié)主要模擬HIFiRE-5 模型在不同雷諾數(shù)、湍流度和迎角的工況下的轉(zhuǎn)捩結(jié)果，并與Juliano[17-20]等的風(fēng)洞實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比分析。實(shí)驗(yàn)來(lái)流馬赫數(shù)Ma=5.8，來(lái)流總溫T0=410 K，等溫壁面溫度Twall=300 K。變化的單位雷諾數(shù)Re、迎角條件α、以及來(lái)流湍流度Tu，其工況在表1 中給出。

表1 HIFiRE-5 模型實(shí)驗(yàn)工況Table 1 Experimental setup of HIFiRE-5

邊界層發(fā)生轉(zhuǎn)捩，最顯著的特征以及影響是壁面摩擦阻力和熱流系數(shù)的急劇上升，這也是目前判斷邊界層轉(zhuǎn)捩與否的常見(jiàn)方法之一。Juliano 等使用溫敏漆（temperature-sensitive paint，TSP）技術(shù)進(jìn)行測(cè)量分析，溫敏漆顯示的強(qiáng)度與溫度變化成正比，所以可以用來(lái)判斷壁面熱流系數(shù)的相對(duì)大小，如果觀察到熱流急劇升高，即為發(fā)生邊界層轉(zhuǎn)捩現(xiàn)象。圖4（a）顯示了HIFiRE-5 模型在單位雷諾為10.2 × 106/m 工況的壁面溫度變化ΔT分布，對(duì)應(yīng)計(jì)算算例Case3。Case3正是各文獻(xiàn)研究較多的工況，本文亦圍繞這一工況開(kāi)展多種效應(yīng)研究。

圖4 Case3 轉(zhuǎn)捩圖像Fig.4 Transition pattern of Case3

3.1 可壓縮效應(yīng)

高超聲速狀態(tài)下，湍流的可壓縮效應(yīng)十分顯著，需要對(duì)湍流模型進(jìn)行可壓縮修正。本節(jié)主要評(píng)估可壓和不可壓兩個(gè)版本的KDO 轉(zhuǎn)捩模型，即CKDO-tran和KDO-tran 的預(yù)測(cè)能力。圖4（b）和圖4（c）顯示分別的是CKDO-tran 和KDO-tran 預(yù)測(cè)的Stanton（St）分布，可通過(guò)St計(jì)算出實(shí)際熱流。可見(jiàn)，CKDO-tran 預(yù)測(cè)的雙肺葉轉(zhuǎn)捩圖像與實(shí)驗(yàn)符合較好，而在此湍流度下KDO-tran 也預(yù)測(cè)出雙肺葉轉(zhuǎn)捩圖像，但轉(zhuǎn)捩位置顯著推遲了，造成轉(zhuǎn)捩形狀、大小與實(shí)驗(yàn)差距很大。此外，對(duì)KDO-tran 還進(jìn)行了變來(lái)流湍流度計(jì)算，無(wú)論怎么調(diào)整來(lái)流湍流度，都無(wú)法取得與實(shí)驗(yàn)符合的轉(zhuǎn)捩圖像。CKDO-tran 和KDO-tran 的計(jì)算結(jié)果表明，對(duì)于用湍流模型直接計(jì)算轉(zhuǎn)捩的方法，可壓縮效應(yīng)的引入對(duì)高超聲速橫流轉(zhuǎn)捩的預(yù)測(cè)至關(guān)重要。

3.2 多雷諾數(shù)工況評(píng)估

Juliano 的風(fēng)洞實(shí)驗(yàn)在保證多數(shù)來(lái)流參數(shù)的條件下，進(jìn)行了多個(gè)變雷諾數(shù)工況研究。文獻(xiàn)多是對(duì)Case3 的雷諾數(shù)工況進(jìn)行研究，本文將轉(zhuǎn)捩模型拓展至更多雷諾數(shù)的研究，一方面通過(guò)更多工況驗(yàn)證模型的適用性、另一方面為他人的對(duì)照研究提供參考。圖5、圖6 給出Case1、Case2 工況下，實(shí)驗(yàn)溫敏漆顯示的轉(zhuǎn)捩圖像和CKDO-tran 計(jì)算的St云圖。聯(lián)合Case1、Case2、Case3 的實(shí)驗(yàn)轉(zhuǎn)捩圖像來(lái)看，隨著雷諾數(shù)的增加，轉(zhuǎn)捩形成的高溫雙肺葉區(qū)域面積增加，并且形狀變得更加瘦長(zhǎng)。CKDO-tran 的計(jì)算結(jié)果符合這種趨勢(shì)，且預(yù)測(cè)的轉(zhuǎn)捩圖像與實(shí)驗(yàn)總體符合較好。

圖5 Case1 轉(zhuǎn)捩圖像Fig.5 Transition pattern of Case1

圖6 Case2 轉(zhuǎn)捩圖像Fig.6 Transition pattern of Case2

3.3 迎角效應(yīng)

盡管實(shí)驗(yàn)給出了帶迎角的HIFiRE-5 高超聲速轉(zhuǎn)捩工況，文獻(xiàn)中卻缺乏相應(yīng)的轉(zhuǎn)捩模型研究，可見(jiàn)帶迎角的高超聲速轉(zhuǎn)捩具有一定挑戰(zhàn)性。本節(jié)在此進(jìn)行Case4 計(jì)算評(píng)估，用以展示迎角效應(yīng)、評(píng)估模型。

圖7 給出Case4 工況下，實(shí)驗(yàn)溫敏漆顯示的轉(zhuǎn)捩圖像和CKDO-tran 計(jì)算的St云圖。轉(zhuǎn)捩圖像來(lái)自于HIFiRE-5 模型的迎風(fēng)面，此時(shí)轉(zhuǎn)捩圖像從原來(lái)的雙肺葉形狀變成頂部帶尖角的鐘形。CKDO-tran 預(yù)測(cè)的轉(zhuǎn)捩圖像與實(shí)驗(yàn)差距較大，主要是因?yàn)檗D(zhuǎn)捩位置整體提前。然而，CKDO-tran 捕捉到的中軸線上的藍(lán)色細(xì)長(zhǎng)泡、藍(lán)色細(xì)長(zhǎng)泡上下方的細(xì)長(zhǎng)泡（低溫區(qū)域）、以及隨后的鐘形及尖角，都與實(shí)驗(yàn)圖像一一對(duì)應(yīng)，表明CKDO-tran 預(yù)測(cè)的轉(zhuǎn)捩圖像與實(shí)驗(yàn)具有相似性。向星皓等[9]對(duì)高超聲速尖錐橫流轉(zhuǎn)捩的計(jì)算，同樣發(fā)現(xiàn)了隨著迎角增大，計(jì)算的轉(zhuǎn)捩位置靠前、轉(zhuǎn)捩圖像相似但有所變形的現(xiàn)象?？梢?jiàn)，真實(shí)情況的迎角效應(yīng)是轉(zhuǎn)捩圖像變化；而湍流/轉(zhuǎn)捩模型的迎角效應(yīng)，是隨著迎角增加，預(yù)測(cè)精度下降。

圖7 Case4 轉(zhuǎn)捩圖像Fig.7 Transition pattern of Case4

鑒于CKDO-tran 計(jì)算轉(zhuǎn)捩的機(jī)制是流動(dòng)結(jié)構(gòu)自適應(yīng)自發(fā)捕捉，并已考慮可壓縮效應(yīng)，變迎角導(dǎo)致的迎風(fēng)面壓縮性變化不應(yīng)導(dǎo)致如此大的計(jì)算誤差。鑒于各向同性湍流-正激波這種最簡(jiǎn)單的湍流過(guò)激波問(wèn)題，各種標(biāo)準(zhǔn)湍流模型的計(jì)算誤差都很大[21]，本文推測(cè)，誤差應(yīng)來(lái)自來(lái)流湍流與頭部脫體激波的相互作用。兩者間的作用在迎角0°時(shí)導(dǎo)致的誤差尚可接受，一旦有迎角，湍流-激波作用的各向異性造成的波后湍流度不均勻性就會(huì)增強(qiáng)，反映在迎風(fēng)面上的轉(zhuǎn)捩面的劇烈變化。因此，為了精準(zhǔn)預(yù)測(cè)高超聲速橫流轉(zhuǎn)捩，需加強(qiáng)湍動(dòng)能或間歇因子過(guò)激波過(guò)程的建模。通過(guò)改進(jìn)來(lái)流湍流過(guò)激波的建模，也能提高轉(zhuǎn)捩模型對(duì)迎角0°工況的計(jì)算精度。

3.4 湍流度的影響

圖8 給出了湍流度對(duì)轉(zhuǎn)捩的影響。在前文，分析了湍流-激波作用的影響，認(rèn)為湍流-激波作用建模的誤差，導(dǎo)致飛行器頭部感受到的來(lái)流湍流度失真,失真的湍流度造成失真的轉(zhuǎn)捩圖像。這一觀點(diǎn)成立的前提是，來(lái)流湍流度是影響轉(zhuǎn)捩的主要因素之一。事實(shí)上，轉(zhuǎn)捩圖像對(duì)湍流度的變化十分敏感，本節(jié)在此展示。仍然采用帶迎角的Case4 算例，僅變化其中的湍流度，形成工況Case5 和Case6，計(jì)算的轉(zhuǎn)捩圖像分別見(jiàn)圖8（a）和圖8（b）?？梢?jiàn)，隨著湍流度的降低，CKDO-tran 計(jì)算的轉(zhuǎn)捩位置逐漸后移；轉(zhuǎn)捩圖像逐漸偏離實(shí)驗(yàn)結(jié)果，與實(shí)驗(yàn)圖像的相似性最終消失。

圖8 湍流度對(duì)轉(zhuǎn)捩的影響Fig.8 Effect of turbulence intensity on transition

4 結(jié)論

本文采用HIFiRE-5 經(jīng)典標(biāo)模，評(píng)估了CKDO-tran模型預(yù)測(cè)高超聲速橫流轉(zhuǎn)捩的性能，并研究了可壓縮因素、雷諾數(shù)因素和迎角因素對(duì)轉(zhuǎn)捩的影響。結(jié)論如下:

1）CKDO-tran 的高超聲速橫流轉(zhuǎn)捩預(yù)測(cè)精度顯著高于KDO-tran 模型，表明可壓縮效應(yīng)的捕捉，對(duì)基于流動(dòng)結(jié)構(gòu)自適應(yīng)的轉(zhuǎn)捩模型至關(guān)重要；

2）采用CKDO-tran 計(jì)算了迎角0°、多雷諾數(shù)的HIFiRE-5 標(biāo)模，豐富了雷諾數(shù)工況，為他人研究提供了參考。計(jì)算得出的橫流轉(zhuǎn)捩與實(shí)驗(yàn)符合較好，表明該模型對(duì)0°迎角的橫流轉(zhuǎn)捩具有較好的適用性；

3）4°迎角下，CKDO-tran 計(jì)算的橫流轉(zhuǎn)捩圖樣結(jié)構(gòu)與實(shí)驗(yàn)相似，但轉(zhuǎn)捩線過(guò)于提前，表明CKDO-tran對(duì)非0°迎角的高超聲速橫流轉(zhuǎn)捩的預(yù)測(cè)具有一定潛力，但仍需進(jìn)一步完善該模型。

4）基于間歇因子γ的轉(zhuǎn)捩模型，用湍流模型提供全湍流解，用γ指示層流區(qū)域和轉(zhuǎn)捩區(qū)域。然而，γ方程里顯含k與μt，這正是來(lái)流k與μt影響轉(zhuǎn)捩位置的原理。基于γ的轉(zhuǎn)捩模型是有可能在k方程過(guò)激波計(jì)算失誤的情況下，通過(guò)調(diào)整γ方程的其他參數(shù)而獲得準(zhǔn)確的預(yù)測(cè)結(jié)果的。因此，間歇因子轉(zhuǎn)捩模型的轉(zhuǎn)捩位置同時(shí)取決于湍流方程和γ方程。

本文完全通過(guò)湍流方程的演化捕捉轉(zhuǎn)捩，排除了γ方程的因素，所得的結(jié)論完全針對(duì)于湍流方程，可為各模型的湍流方程的改進(jìn)提供可靠的信息，如需加強(qiáng)來(lái)流湍流過(guò)激波的研究等。