夏陳超，姜婷婷，郭中州，陳偉芳

(浙江大學(xué)航空航天學(xué)院，浙江杭州 310027)

壓縮性修正對(duì)γ-Reθ轉(zhuǎn)捩模型的影響研究

夏陳超，姜婷婷，郭中州，陳偉芳*

(浙江大學(xué)航空航天學(xué)院，浙江杭州 310027)

將基于流場(chǎng)當(dāng)?shù)刈兞康摩?Reθ轉(zhuǎn)捩模型加入可壓縮RANS方程計(jì)算程序，對(duì)平均流場(chǎng)控制方程、湍流和轉(zhuǎn)捩控制方程進(jìn)行基于LU-SGS的隱式緊耦合求解。結(jié)合三種壓縮性修正方法對(duì)高超聲速平板、雙楔和圓錐流動(dòng)進(jìn)行了數(shù)值模擬，給出了壁面斯坦頓數(shù)、熱流值和實(shí)驗(yàn)值的比較，以及湍動(dòng)能和間歇因子分布云圖。數(shù)值計(jì)算結(jié)果表明，相同來(lái)流湍流度下不同壓縮性修正方法得到的轉(zhuǎn)捩位置差別較大。相比于未修正的模型，基于當(dāng)?shù)伛R赫數(shù)壓縮性修正后的轉(zhuǎn)捩位置最靠后，其對(duì)帶有分離的雙楔流動(dòng)預(yù)測(cè)較為準(zhǔn)確;基于來(lái)流馬赫數(shù)的壓縮性修正僅使得轉(zhuǎn)捩位置稍有延遲;對(duì)湍動(dòng)能源項(xiàng)的壓縮性修正可提高模型在轉(zhuǎn)捩后湍流段的熱流預(yù)測(cè)精度。

高超聲速;RANS方程;γ-Reθ轉(zhuǎn)捩模型;緊耦合;壓縮性修正

0 引言

邊界層轉(zhuǎn)捩是目前流體力學(xué)領(lǐng)域的研究熱點(diǎn)和難點(diǎn)，也是航空航天工程中的關(guān)鍵問(wèn)題。雖然國(guó)內(nèi)外研究者已經(jīng)做了大量卓有成效的工作，但對(duì)轉(zhuǎn)捩相關(guān)機(jī)理的認(rèn)識(shí)和理論分析仍有待提升。對(duì)高超聲速飛行器而言，流動(dòng)形態(tài)的變化關(guān)系到飛行器的升阻特性、熱防護(hù)以及發(fā)動(dòng)機(jī)中的流動(dòng)品質(zhì)等。準(zhǔn)確預(yù)測(cè)轉(zhuǎn)捩位置、轉(zhuǎn)捩區(qū)長(zhǎng)度以及氣動(dòng)力/熱環(huán)境對(duì)飛行器設(shè)計(jì)有重要意義。

相比于直接數(shù)值模擬(DNS)和大渦模擬(LES)，基于RANS的數(shù)值模擬是當(dāng)前工程上預(yù)測(cè)湍流和轉(zhuǎn)捩較為經(jīng)濟(jì)和高效的方法，其中Huang［1］、Walters［2］和Lantry［3］等的研究具有代表性。Lantry等提出的基于經(jīng)驗(yàn)關(guān)系式的γ-Reθ轉(zhuǎn)捩模型在k-ω SST兩方程湍流模型的基礎(chǔ)上添加了兩個(gè)輸運(yùn)方程，一個(gè)是間歇因子γ的輸運(yùn)方程，一個(gè)是動(dòng)量厚度雷諾數(shù)Reθ的輸運(yùn)方程。γ-Reθ模型是基于流場(chǎng)當(dāng)?shù)刈兞康霓D(zhuǎn)捩模型，可以和現(xiàn)代CFD程序相兼容，相比于目前常用的eN方法［4］，具有更強(qiáng)的適用性。針對(duì)低速平板、翼型/機(jī)翼和渦輪機(jī)械的眾多研究［5-7］表明γ-Reθ模型對(duì)自然轉(zhuǎn)捩和分離轉(zhuǎn)捩具有較好的預(yù)測(cè)能力。

一方面，雖然原始的γ-Reθ模型基于不可壓實(shí)驗(yàn)發(fā)展而來(lái)，但 Krause［8］、You［9］、孔維萱［10］以及鄭赟［11］等采用γ-Reθ模型對(duì)高超聲速雙楔和尖錐等流動(dòng)的數(shù)值計(jì)算表明該模型也能較準(zhǔn)確地預(yù)測(cè)高超聲速流動(dòng)中的轉(zhuǎn)捩問(wèn)題。另一方面，為增強(qiáng)該模型對(duì)高速流動(dòng)問(wèn)題的適用性，有必要進(jìn)行適當(dāng)?shù)膲嚎s性修正。Kaynak［12］給出了一種基于來(lái)流馬赫數(shù)的壓縮性修正關(guān)系式，對(duì)超聲速平板的邊界層轉(zhuǎn)捩進(jìn)行了預(yù)測(cè);張曉東［13］則對(duì)經(jīng)驗(yàn)關(guān)系式進(jìn)行了基于當(dāng)?shù)伛R赫數(shù)的修正，通過(guò)對(duì)雙楔模型的數(shù)值計(jì)算，得到的壁面壓力系數(shù)與實(shí)驗(yàn)值吻合較好?？傮w而言，該模型目前主要應(yīng)用于低速流動(dòng)問(wèn)題，其在高速流動(dòng)情況下的壓縮性修正研究還很少。本文將γ-Reθ模型編入課題組發(fā)展的基于SST湍流模型的可壓縮RANS方程耦合求解程序，結(jié)合高超聲速平板和圓錐等典型流動(dòng)問(wèn)題，通過(guò)數(shù)值模擬和對(duì)比，初步探索了三種壓縮性修正方法對(duì)γ-Reθ模型預(yù)測(cè)轉(zhuǎn)捩的影響，為該模型的進(jìn)一步改進(jìn)提供參考。

1 轉(zhuǎn)捩模型

1.1 γ-Reθ轉(zhuǎn)捩模型

γ-Reθ轉(zhuǎn)捩模型包含間歇因子和轉(zhuǎn)捩動(dòng)量厚度雷諾數(shù)兩個(gè)輸運(yùn)方程，其中間歇因子的輸運(yùn)方程為:

其中，F(xiàn)length為轉(zhuǎn)捩區(qū)長(zhǎng)度控制函數(shù)，F(xiàn)onset為轉(zhuǎn)捩發(fā)生函數(shù)，F(xiàn)turb用于抑制再層流化，相關(guān)參數(shù)詳見文獻(xiàn)［3］。

1.2 壓縮性修正

由于原始的γ-Reθ模型中轉(zhuǎn)捩動(dòng)量厚度雷諾數(shù)Reθt關(guān)系式是通過(guò)不可壓平板實(shí)驗(yàn)獲得，因此將該模型應(yīng)用于高速可壓縮流問(wèn)題時(shí)有必要進(jìn)行適當(dāng)?shù)目蓧嚎s性修正(Compressibility Correction，文中用CC表示)。本文實(shí)現(xiàn)并對(duì)比了三種壓縮性修正方法。

(1)CC=1

張曉東［13］基于高超聲速風(fēng)洞實(shí)驗(yàn)數(shù)據(jù)，給出了一種基于當(dāng)?shù)伛R赫數(shù)的壓縮性修正關(guān)系式:

(2)CC=2

Kaynak［12］給出了一種基于來(lái)流馬赫數(shù)的修正關(guān)系式:

雖然文獻(xiàn)［12］給出了該關(guān)系式的適用范圍(來(lái)流馬赫數(shù)小于2.4)，本文仍將其視為一種潛在的修正方法，評(píng)估其壓縮性修正效果。

(3)CC=3

Papp等人［14］在對(duì)其發(fā)展的基于SSGZ k-ε模型的工程轉(zhuǎn)捩模型進(jìn)行壓縮性修正時(shí)在湍動(dòng)能方程中添加了一個(gè)源項(xiàng)，形式如下:

本文參考該方法，將上述源項(xiàng)加入SST湍流模型的湍動(dòng)能方程中進(jìn)行修正。其中，，湍流馬赫數(shù)，M't為考慮轉(zhuǎn)捩效應(yīng)的轉(zhuǎn)捩湍流馬赫數(shù)，系數(shù)λ的作用是當(dāng)轉(zhuǎn)捩湍流馬赫數(shù)小于λ時(shí)不起作用，取值0.2，α1=2.5，α2=2.0，ε=βωk。

壓縮性修正1和2的實(shí)施是將F(Ma)乘到轉(zhuǎn)捩動(dòng)量厚度雷諾數(shù)，壓縮性修正3則是將Pe直接加到湍動(dòng)能k方程的源項(xiàng)中。

2 數(shù)值方法

將 γ-Reθ轉(zhuǎn)捩模型加入課題組發(fā)展的可壓縮RANS程序，流動(dòng)控制方程為雷諾平均的三維NS方程、SST兩方程湍流模型方程以及γ-Reθ轉(zhuǎn)捩輸運(yùn)方程。計(jì)算程序基于結(jié)構(gòu)網(wǎng)格的有限體積方法，其中無(wú)黏通量采用計(jì)算效率和分辨率均較高的AUSMPW+格式［15］并采用MUSCL進(jìn)行高階插值，黏性通量采用中心差分進(jìn)行離散。

為提高計(jì)算效率，時(shí)間推進(jìn)采用Yoon和Jameson［16］提出的LU-SGS隱式方法，同時(shí)采用當(dāng)?shù)貢r(shí)間步長(zhǎng)加速收斂。為克服源項(xiàng)剛性，對(duì)湍流和轉(zhuǎn)捩輸運(yùn)方程中的負(fù)值耗散項(xiàng)進(jìn)行隱式處理以增強(qiáng)矩陣對(duì)角占優(yōu)［12］:

3 計(jì)算結(jié)果與分析

3.1 平板流動(dòng)

高超聲速平板流動(dòng)算例取自Mee［17］的激波風(fēng)洞實(shí)驗(yàn)結(jié)果。平板長(zhǎng)1.5m，來(lái)流馬赫數(shù)6.2，來(lái)流靜溫690K，來(lái)流靜壓5.4kPa，單位雷諾數(shù)為2.6×106。計(jì)算網(wǎng)格為224×150(流向和法向)，在壁面和前緣進(jìn)行了加密，壁面第一層網(wǎng)格y+小于1(約0.5)。

圖1和圖2為兩種來(lái)流湍流度情況下的平板壁面斯坦頓數(shù)分布，斯坦頓數(shù)定義為:

其中，qw為熱流，ρ∞和u∞分別為來(lái)流密度和速度，h0∞為來(lái)流總焓，hw壁面焓。圖中對(duì)比了無(wú)壓縮性修正的γ-Reθ模型、三種壓縮性修正的γ-Reθ模型、完全層流狀態(tài)和完全湍流狀態(tài)(SST模型)的結(jié)果。可以看出，層流和湍流的斯坦頓數(shù)與實(shí)驗(yàn)值吻合較好，轉(zhuǎn)捩模型實(shí)現(xiàn)了層流到湍流的過(guò)渡，隨著來(lái)流湍流度的增大，轉(zhuǎn)捩位置前移，但相比于實(shí)驗(yàn)值，轉(zhuǎn)捩區(qū)長(zhǎng)度過(guò)長(zhǎng)。相同湍流度條件下不同壓縮性修正方法的轉(zhuǎn)捩起始位置差異較大，方法1轉(zhuǎn)捩位置最靠后，方法2居中，相比于無(wú)壓縮性模型均延遲了轉(zhuǎn)捩位置，方法3由于沒有改變轉(zhuǎn)捩模型經(jīng)驗(yàn)關(guān)系式，其轉(zhuǎn)捩位置與無(wú)壓縮性模型基本一致，但其在一定程度上降低了湍流段的熱流。

圖3和圖4分別為湍流度3%時(shí)無(wú)壓縮性修正和采用壓縮性方法2的流場(chǎng)湍動(dòng)能云圖，湍動(dòng)能的急劇增大意味著轉(zhuǎn)捩的發(fā)生，可以看出，進(jìn)行壓縮性修正后轉(zhuǎn)捩位置有了明顯的后移，即該壓縮性修正方法限制了湍動(dòng)能的發(fā)展。數(shù)值計(jì)算也表明，無(wú)壓縮性修正時(shí)最大湍動(dòng)能為61 611m2/s2，壓縮性修正方法2得到的最大湍動(dòng)能減少到了60 257 m2/s2。

圖1 湍流度2%的平板壁面斯坦頓數(shù)Fig.1 Stanton number on flat plate(Tu=2%)

圖2 湍流度3%的平板壁面斯坦頓數(shù)Fig.2 Stanton number on flat plate(Tu=3%)

圖3 無(wú)壓縮性修正的平板湍動(dòng)能云圖Fig.3 Contour of turbulence kinetic energy without compressibility correction

圖4 壓縮性方法2的平板湍動(dòng)能云圖Fig.4 Contour of turbulence kinetic energy with compressibility correction method 2

3.2 雙楔流動(dòng)

雙楔幾何外形［18］如圖5所示，其第一個(gè)和第二個(gè)坡的角度分別為 9°和 20.5°，前緣鈍化半徑0.5mm。計(jì)算馬赫數(shù)為8.1，來(lái)流單位雷諾數(shù)為3.8 ×106，來(lái)流靜壓520Pa，來(lái)流靜溫106K，壁面溫度為300K，來(lái)流湍流度0.9%。鈍頭雙楔計(jì)算網(wǎng)格如圖6所示，網(wǎng)格在壁面、頭部和拐角附近進(jìn)行了加密，并保證第一層y+值小于1(約為0.3)。

圖5 雙楔流動(dòng)示意圖Fig.5 Schematic of blunted double wedge flow

圖6 雙楔模型計(jì)算網(wǎng)格Fig.6 Grid of blunted double wedge

圖7 雙楔模型壁面壓力分布Fig.7 Comparison of pressure coefficient of double wedge with blunt leading edge

圖7 和圖8分別為雙楔模型壁面壓力系數(shù)和斯坦頓數(shù)分布。可以看出，完全層流狀態(tài)在拐角的分離區(qū)過(guò)大，斯坦頓數(shù)過(guò)低，完全湍流狀態(tài)則未發(fā)生分離，斯坦頓數(shù)過(guò)高。轉(zhuǎn)捩模型的結(jié)果與實(shí)驗(yàn)值更接近，其中，壓縮性修正方法1的分離區(qū)大小、壓力和斯坦頓數(shù)與實(shí)驗(yàn)值最為吻合，壓縮性修正方法3與未修正的轉(zhuǎn)捩模型相比，拐角后的湍流段壁面斯坦頓數(shù)較低。

圖8 雙楔模型壁面斯坦頓數(shù)分布Fig.8 Comparison of Stanton number of double wedge with blunt leading edge

圖9 至圖11給出了不同壓縮性方法在拐角處的局部放大流線圖和間歇因子云圖，可以看出，間歇因子分布和流動(dòng)的分離是相對(duì)應(yīng)的，分離起始點(diǎn)之前的間歇因子保持很小的值，表明流動(dòng)還是層流;分離之后，間歇因子迅速增長(zhǎng)。從壓縮性修正方法1到3，分離區(qū)逐漸減小，壓縮性減弱，這與壁面壓力和熱流的分布趨勢(shì)是一致的。

圖9 壓縮性方法1拐角處流線和間歇因子Fig.9 Streamline and intermittency factor around corner with compressibility correction method 1

圖10 壓縮性方法2拐角處流線和間歇因子Fig.10 Streamline and intermittency factor around corner with compressibility correction method 2

圖11 壓縮性方法3拐角處流線和間歇因子Fig.11 Streamline and intermittency factor around corner with compressibility correction method 3

3.3 圓錐流動(dòng)

針對(duì)高超聲速圓錐流動(dòng)［19］進(jìn)行數(shù)值計(jì)算，圓錐錐角為7°，頭部鈍化半徑5 mm，如圖12所示。來(lái)流馬赫數(shù)7.15，來(lái)流靜溫214 K，壁溫300 K，來(lái)流靜壓7 722 Pa，單位雷諾數(shù)為9.64×106。采用半模計(jì)算，網(wǎng)格量為106×78×29(流向、法向和周向)，壁面第一層網(wǎng)格間距為1×10-6m。

圖12 圓錐幾何外形Fig.12 Geometry of cone

由于不易實(shí)現(xiàn)計(jì)算與實(shí)驗(yàn)條件的完全一致，本文計(jì)算了多個(gè)來(lái)流湍流度的情況，圖13至圖15分別為來(lái)流湍流度0.2%、0.35%和0.5%條件下的圓錐表面熱流分布對(duì)比?？梢钥闯?，完全層流的計(jì)算結(jié)果在層流段與實(shí)驗(yàn)值吻合，完全湍流的計(jì)算結(jié)果在湍流段熱流比實(shí)驗(yàn)值偏高，轉(zhuǎn)捩模型的結(jié)果與實(shí)驗(yàn)值更為接近，且隨著來(lái)流湍流度的提高，轉(zhuǎn)捩位置前移。相同湍流度條件下三種壓縮性修正方法的轉(zhuǎn)捩位置相差較大，方法1轉(zhuǎn)捩位置最靠后(來(lái)流湍流度為0.2%時(shí)未發(fā)生轉(zhuǎn)捩)，方法2居中，方法3轉(zhuǎn)捩位置與無(wú)壓縮性修正基本相同，趨勢(shì)與3.1中的平板一致。壓縮性修正方法3在轉(zhuǎn)捩后的湍流段熱量值較低，與實(shí)驗(yàn)值更接近。

圖13 來(lái)流湍流度0.2%壁面熱流分布Fig.13 Heat flux on wall(Tu=0.2%)

圖14 來(lái)流湍流度0.35%壁面熱流分布Fig.14 Heat flux on wall(Tu=0.35%)

圖15 來(lái)流湍流度0.5%壁面熱流分布Fig.15 Heat flux on wall(Tu=0.5%)

4 結(jié)論

基于可壓縮RANS方程計(jì)算程序，引入γ-Reθ轉(zhuǎn)捩模型，采用三種壓縮性修正方法對(duì)高超聲速平板、雙楔和圓錐進(jìn)行了數(shù)值模擬并與實(shí)驗(yàn)結(jié)果進(jìn)行了比較。研究表明，三種壓縮性修正后的轉(zhuǎn)捩模型在相同湍流度條件下的轉(zhuǎn)捩位置相差較大?；诋?dāng)?shù)伛R赫數(shù)的壓縮性修正有效地預(yù)測(cè)了帶有分離的雙楔流動(dòng)，同時(shí)，相比于未修正的轉(zhuǎn)捩模型，該修正會(huì)較大程度地延遲轉(zhuǎn)捩位置?；趤?lái)流馬赫數(shù)的壓縮性修正僅使得轉(zhuǎn)捩位置稍有延遲。湍動(dòng)能方程源項(xiàng)的壓縮性修正能降低轉(zhuǎn)捩后湍流段壁面熱流，使結(jié)果與實(shí)驗(yàn)值更加接近，對(duì)轉(zhuǎn)捩位置沒有影響。此外還可看出，對(duì)于不同的流動(dòng)特性，相同的壓縮性修正方式效果是不一樣的，即很難得到對(duì)任何流動(dòng)情況都適用的模型，實(shí)際情況中，應(yīng)根據(jù)不同的流動(dòng)特性進(jìn)行湍流和轉(zhuǎn)捩模型的選擇。

總體而言，要將γ-Reθ轉(zhuǎn)捩模型應(yīng)用于高速流動(dòng)問(wèn)題還有許多值得改進(jìn)的工作。相比于壁面壓力系數(shù)，高超聲速條件下壁面熱流值預(yù)測(cè)精度的提高仍是一個(gè)挑戰(zhàn)。本文中研究的壓縮性修正方法均較為簡(jiǎn)單，將當(dāng)?shù)伛R赫數(shù)和湍動(dòng)能源項(xiàng)修正相結(jié)合可能是一個(gè)更好的嘗試。考慮更多轉(zhuǎn)捩影響因素、適用范圍更廣的壓縮性修正方法是今后的一個(gè)探索方向。

參考文獻(xiàn):

［1］Suzen Y B，Huang P G.An intermittency transport equation for modeling flow transition［R］.AIAA 2000-0287.2000

［2］Walters D K，Leylek J H.A new model for boundary-layer transition using a single-point RANS approach［J］.Journal of Turbomachinery，2002，126(1):193-202.

［3］Langtry R B，Menter F R.Transition modeling for general CFD applications in aeronautics［C］.Reno:43rd AIAA Aerospace Sciences Meeting and Exhibit，2005.

［4］Malik M，Balakumar P.Instability and transition in three-dimensional supersonic boundary layers［R］.AIAA-92-5049.1992

［5］Malan P，Suluksna K，Juntasaro E.Calibrating the γ-Reθtransition model for commercial CFD［R］.AIAA 2009-1142.2009

［6］Langtry R B，Menter F R.Correlation-based transition modeling for unstructured parallelized computational fluid dynamics codes［J］.AIAA Journal，2009，47(12):2894-2906.

［7］Wang Gang，Liu Yi，Wang Guangqiu，et al.Transitional flow simulation based on γ-Reθttransition model［J］.Acta Aeronautica et Astronautica Sinica，2014，35(1):70-79.(in Chinese)王剛，劉毅，王光秋，等.采用γ-Reθt模型的轉(zhuǎn)捩流動(dòng)計(jì)算分析［J］.航空學(xué)報(bào).2014，35(1):70-79.

［8］Krause M，Behr M，Ballmann J.Modeling of transition effects in hypersonic intake flows using a correlation-based intermittency model［C］.Dayton:15th AIAA International Space Planes and Hypersonic Systems and Technologies Conference，2008.

［9］You Y，Luedeke H，Eggers T，Hannemann K.Application of the γ-Reθttransition model in high speed flows［C］.Tours:18th AIAA/3AF International Space Planes and Hypersonic Systems and Technologies Conference，2012.

［10］Kong Weixuan，Yan Chao，Zhao Rui.γ-Reθmodel research for high-speed boundary layer transition［J］.Acta Aerodynamica Sinica，2013，31(1):120-126.(in Chinese)孔維萱，閻超，趙瑞.γ-Reθ模式應(yīng)用于高速邊界層轉(zhuǎn)捩的研究［J］.空氣動(dòng)力學(xué)學(xué)報(bào).2013，31(1):120-126.

［11］Zheng Yun，Li Hongyang，Liu Daxiang.Application and analysis of γ-Reθtransition model in hypersonic flow［J］.Journal of Propulsion Technology，2014，35(3):296-304.(in Chinese)鄭赟，李虹楊，劉大響.γ-Reθ轉(zhuǎn)捩模型在高超聲速下的應(yīng)用及分析［J］.推進(jìn)技術(shù)，2014，35(3):296-304.

［12］Kaynak U.Supersonic boundary-layer transition prediction under the effect of compressibility using a correlation-based model［J］.Proceedings of the Institution of Mechanical Engineers，Part G:Journal of Aerospace Engineering，2012，226(7):722-739.

［13］Zhang X D，Gao Z H.Numerical discuss to complete empirical correlation in Langtry＇s transition model［J］.Applied Mathematics and Mechanics，2010，31(5):544-552.(in Chinese)張曉東，高正紅.關(guān)于補(bǔ)充Langtry的轉(zhuǎn)捩模型經(jīng)驗(yàn)修正式的數(shù)值探討［J］.應(yīng)用數(shù)學(xué)和力學(xué).2010，31(5):544-552.

［14］Papp J L，Kenzakowski D C，Dash S M.Extensions of a rapid engineering approach to modeling hypersonic laminar to turbulent transitional flows［R］.AIAA 2005-0892.

［15］Kim K H，Kim C，Rho O.Methods for the accurate computations of hypersonic flows:I.AUSMPW+scheme［J］.Journal of Computational Physics，2001，174(1):38-80.

［16］Yoon S，Jameson A.Lower-upper symmetric-Gauss-Seidel method for the Euler and Navier-Stokes equations［J］.AIAA Journal，1988，26(9):1025-1026.

［17］Mee D J.Boundary-layer transition measurements in hypervelocity flows in a shock tunnel［J］.AIAA Journal，2002，40(8):1542-1548.

［18］Neuenhahn T，Olivier H.Influence of the wall temperature and the entropy layer effects on double wedge shock boundary layer interactions［R］.AIAA 2006-8136.

［19］Maclean M，Mundy E，Wadhams T，Holden M，Johnson H，Candler G.Comparisons of transition prediction using PSE-chem to measurements for a shock tunnel environment［R］.AIAA 2007-4490.

Effects of compressibility correction on γ-Reθtransition model

Xia Chenchao，Jiang Tingting，Guo Zhongzhou，Chen Weifang*
(School of Aeronautics and Astronautics，Zhejiang University，Hangzhou 310027，China)

Boundary layer transition plays a significant role in the prediction of aerodynamic and aerothermodynamics characteristics of hypersonic vehicle.The commonly used correlation based on γ-Reθtransition model is implemented into a Reynolds-Averaged Navier-Stokes solver in order to evaluate its capability in hypersonic flows.The mean flow and turbulent flow equations are solved simultaneously based on the LU-SGS method.Three compressibility correction methods are adopted to simulate the hypersonic flows around a flat plate，a double wedge and a cone body.Results of Stanton number，heat flux on walls，contours of turbulence kinetic energy and intermittency factor are presented.The results of γ-Reθtransition model show better flow essential than that of full laminar or full turbulent results.The transition onset positions obtained by different compressibility methods under the same free stream turbulence intensity varied significantly.Compressibility correction based on local Mach number delays the transition onset greatly，while the correction of source term of turbulent kinetic energy equation improves the accuracy of heat flux on turbulent section of the wall.Rational selection of compressibility correction method should rely on the type of flow.Elaborate and reliable compressibility correction method still need to be further investigated.

hypersonic;RANS equations;γ-Reθtransition model;coupled method;compressibility correction

V211.3

:Adoi:10.7638/kqdlxxb-2014.0068

0258-1825(2015)05-0603-07

2014-07-14;

:2014-10-10

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(2014CB340201)

夏陳超(1989-)，男，福建寧德人，博士生，流體力學(xué)專業(yè).E-mail:aeroxia@zju.edu.cn

陳偉芳*(1970-)，男，博士，教授，研究方向?yàn)楦叱曀倏諝鈩?dòng)力學(xué).E-mail:chenwfnudt@163.com

夏陳超，姜婷婷，郭中州，等.壓縮性修正對(duì)γ-Reθ轉(zhuǎn)捩模型的影響研究［J］.空氣動(dòng)力學(xué)學(xué)報(bào)，2015，33(5):603-609.

10.7638/kqdlxxb-2014.0068 Xia C C，Jiang T T，Guo Z Z，et al.Effects of compressibility correction on γ-Reθtransition model［J］.Acta Aerodynamica Sinica，2015，33(5):603-609.