李 強(qiáng)，張扣立，莊 宇，趙金山

(中國空氣動力研究與發(fā)展中心超高速所，綿陽 621000)

?

激波風(fēng)洞邊界層強(qiáng)制轉(zhuǎn)捩試驗研究

李 強(qiáng)，張扣立，莊 宇，趙金山

(中國空氣動力研究與發(fā)展中心超高速所，綿陽 621000)

針對升力體模型設(shè)計了渦流發(fā)生器，在中國空氣動力研究與發(fā)展中心(CARDC)Φ2 m激波風(fēng)洞上開展風(fēng)洞試驗，研究了高超聲速邊界層強(qiáng)制轉(zhuǎn)捩問題。試驗來流名義馬赫數(shù)分別為10、12，單位雷諾數(shù)分別為2.4×106/m、2.1×106/m，模型攻角10°。試驗中應(yīng)用鉑薄膜熱流傳感器技術(shù)和溫敏熱圖(TSP)技術(shù)測量了模型表面熱流，證明渦流發(fā)生器實(shí)現(xiàn)了模型邊界層強(qiáng)制轉(zhuǎn)捩，使Φ2 m激波風(fēng)洞擁有了模擬高馬赫數(shù)低雷諾數(shù)湍流邊界層的能力。試驗結(jié)果表明,不同形狀不同高度渦流發(fā)生器對邊界層完全轉(zhuǎn)捩成湍流后的熱流影響不明顯，由此可提出一種新的激波風(fēng)洞試驗方法，即利用渦流發(fā)生器開展相同來流條件下不同邊界層流態(tài)對模型表面熱流等邊界層參數(shù)分布影響的試驗研究。

激波風(fēng)洞；渦流發(fā)生器；強(qiáng)制轉(zhuǎn)捩；湍流模擬

0 引 言

激波風(fēng)洞等大型地面氣動熱環(huán)境模擬試驗設(shè)備，由于受到風(fēng)洞運(yùn)行參數(shù)的制約，對飛行條件的模擬受到較大限制。而風(fēng)洞模型長度相對較短，某些條件下模型表面邊界層無法轉(zhuǎn)捩形成湍流，或者無法在需要的位置轉(zhuǎn)捩形成湍流，導(dǎo)致風(fēng)洞試驗無法模擬飛行器表面邊界層流動狀態(tài)。為了能準(zhǔn)確地將高超聲速飛行器風(fēng)洞試驗結(jié)果外推到全尺寸飛行器，就需要利用邊界層轉(zhuǎn)捩裝置強(qiáng)迫形成湍流流動，使高超聲速邊界層按可控和可預(yù)測的方式轉(zhuǎn)捩。前期的一些CFD計算工作發(fā)現(xiàn)，層流流態(tài)條件下，升力體方向舵舵軸迎風(fēng)面熱流比湍流流態(tài)條件高，這與迎風(fēng)大面積區(qū)域熱流分布規(guī)律相反。分析認(rèn)為導(dǎo)致這種現(xiàn)象的原因可能是相同來流條件下，方向舵位置處的層流邊界層比湍流邊界層薄，因而進(jìn)入方向舵縫隙的氣流速度更高。為了在風(fēng)洞試驗中復(fù)現(xiàn)這一現(xiàn)象，首先需要解決相同來流條件下，實(shí)現(xiàn)模型表面層流和湍流邊界層模擬。這需要在較低單位雷諾數(shù)流場條件下，通過強(qiáng)制轉(zhuǎn)捩控制裝置，將模型表面層流邊界層強(qiáng)制轉(zhuǎn)捩為湍流邊界層。

目前風(fēng)洞試驗中采用較多的強(qiáng)制轉(zhuǎn)捩裝置[1-4]是被動控制裝置，主要通過采用增加模型壁面粗糙度和安裝渦流發(fā)生器的方式，使有限尺寸模型表面邊界層強(qiáng)制轉(zhuǎn)捩，以實(shí)現(xiàn)地面風(fēng)洞設(shè)備按照需求對湍流流動狀態(tài)的模擬。近幾年在中國空氣動力研究與發(fā)展中心超高速所Φ2 m激波風(fēng)洞(FD-14A)中，開展了一系列的強(qiáng)制轉(zhuǎn)捩試驗，采用增加模型壁面粗糙度和安裝渦流發(fā)生器的方式，在風(fēng)洞試驗來流單位雷諾數(shù)1×107/m條件下，使模型表面邊界層轉(zhuǎn)捩位置提前?？紤]到很多文獻(xiàn)[5-9]都提到粗糙元高度是影響強(qiáng)制轉(zhuǎn)捩的關(guān)鍵因素，為了在激波風(fēng)洞更低單位雷諾數(shù)流場條件下，使模型表面層流邊界層強(qiáng)制轉(zhuǎn)捩為湍流邊界層，本文采用渦流發(fā)生器作為強(qiáng)制轉(zhuǎn)捩裝置。

近些年來，吸氣式高超聲速飛行器發(fā)展非常迅速，由于其采用的吸氣式超燃發(fā)動機(jī)要求氣流在進(jìn)入進(jìn)氣道前形成湍流，因此高超聲速邊界層強(qiáng)制轉(zhuǎn)捩試驗技術(shù)得到了深入的研究。美國的Berry等[2]、Borg等[10]、Schneider等[11]，在20inM6風(fēng)洞、31inM10風(fēng)洞、M6靜音風(fēng)洞等風(fēng)洞設(shè)備上對Hyper-X、X- 43和X-51A飛行器模型前體進(jìn)行了高超聲速邊界層控制方法研究，采用了鉆石型(包括圓柱型)和斜坡型(包括不同參數(shù)的變形)渦流發(fā)生器作為強(qiáng)制轉(zhuǎn)捩控制裝置。國內(nèi)的戰(zhàn)培國等[1]、戰(zhàn)培國[12]和趙慧勇等[13]對吸氣式高超聲速飛行器進(jìn)氣道強(qiáng)制轉(zhuǎn)捩問題及風(fēng)洞試驗情況作了比較系統(tǒng)的研究，發(fā)表了綜述性的文章。趙慧勇等[14]在中國空氣動力研究與發(fā)展中心高速所Φ0.5 m高超聲速風(fēng)洞(FL-31)、趙俊波等[15]在中國航天空氣動力技術(shù)研究院Φ0.5 m高超聲速風(fēng)洞(FD-07)分別對超燃進(jìn)氣道邊界層強(qiáng)制轉(zhuǎn)捩進(jìn)行了研究，都獲得了比較好的試驗效果。

上述試驗大多是在常規(guī)高超聲速風(fēng)洞較低馬赫數(shù)流場條件下，針對進(jìn)氣道壓縮面平面區(qū)域開展強(qiáng)制轉(zhuǎn)捩試驗研究工作。本文在脈沖型設(shè)備激波風(fēng)洞高馬赫數(shù)低雷諾數(shù)流場條件下，針對升力體弧面區(qū)域開展強(qiáng)制轉(zhuǎn)捩試驗研究，研究渦流發(fā)生器的強(qiáng)制轉(zhuǎn)捩特性。

1 渦流發(fā)生器強(qiáng)制轉(zhuǎn)捩原理

渦流發(fā)生器依據(jù)粗糙元誘導(dǎo)邊界層轉(zhuǎn)捩，影響粗糙元強(qiáng)制轉(zhuǎn)捩效果的因素非常多，包括粗糙元的高度、外形和間距，另外馬赫數(shù)、雷諾數(shù)、模型攻角、來流噪聲、壁面溫度、逆壓梯度、壁面結(jié)構(gòu)特征等因素也會對粗糙元強(qiáng)制轉(zhuǎn)捩效果帶來影響。粗糙元對高超聲速邊界層轉(zhuǎn)捩影響的研究工作持續(xù)了超過50年，目前仍然是國際上的研究熱點(diǎn)。并且以渦流發(fā)生器為代表的粗糙元研究結(jié)果，已經(jīng)被成功應(yīng)用到吸氣式高超聲速飛行器上。但是目前對粗糙元誘導(dǎo)轉(zhuǎn)捩的物理機(jī)理仍未完全研究清楚，沒有統(tǒng)一完整的解釋。國內(nèi)外研究粗糙元誘導(dǎo)強(qiáng)制轉(zhuǎn)捩機(jī)理的文獻(xiàn)很多，下面挑選部分典型的文獻(xiàn)來說明粗糙元強(qiáng)制轉(zhuǎn)捩機(jī)理問題。

Schneider[16]認(rèn)為粗糙元在強(qiáng)制轉(zhuǎn)捩過程中的作用為粗糙元會在尾跡區(qū)產(chǎn)生流向渦結(jié)構(gòu)和不穩(wěn)定剪切層，流向渦可能通過一些失穩(wěn)機(jī)制發(fā)展，也可能通過感受性機(jī)理，與聲波或者來流的其他擾動相互作用進(jìn)而發(fā)展出失穩(wěn)波。文獻(xiàn)[2,17-20]認(rèn)為在粗糙元上游會形成反向旋轉(zhuǎn)的渦對，下游會形成不穩(wěn)定尾跡渦和側(cè)邊的馬蹄渦結(jié)構(gòu)，尾跡中的渦結(jié)構(gòu)和剪切層相互作用而導(dǎo)致邊界層轉(zhuǎn)捩。流向渦對邊界層的失穩(wěn)和轉(zhuǎn)捩過程十分重要，它通常對應(yīng)于瞬態(tài)增長不穩(wěn)定機(jī)制，因此很多研究[21-22]認(rèn)為粗糙元引發(fā)的邊界層轉(zhuǎn)捩是由瞬態(tài)增長引起的旁路轉(zhuǎn)捩。而Merkle等[23]認(rèn)為，粗糙元對基本流的改變以及其自身提供的擾動在轉(zhuǎn)捩中占據(jù)主導(dǎo)作用。另外，文獻(xiàn)[24]提到，粗糙元的外形對轉(zhuǎn)捩效果影響不大，這也是本文通過風(fēng)洞試驗所驗證的結(jié)果。

國內(nèi)近幾年對粗糙元誘導(dǎo)強(qiáng)制轉(zhuǎn)捩機(jī)理也開展了很多研究工作[25-29]，基本都認(rèn)為粗糙元尾跡流向渦與剪切層的相互作用是強(qiáng)制轉(zhuǎn)捩的主要機(jī)理，并且還分析了鉆石型和斜坡型粗糙元誘導(dǎo)轉(zhuǎn)捩機(jī)理的異同點(diǎn)。天津大學(xué)周恒院士等[30]認(rèn)為，渦流發(fā)生器這種大粗糙元引發(fā)的旁路(bypass)轉(zhuǎn)捩，是由于大粗糙元下游產(chǎn)生的流向渦旋，其作用可能與流向條帶類似，會導(dǎo)致某種二次不穩(wěn)定，由此引發(fā)某些能快速增長的不穩(wěn)定波，這些波會修正平均流剖面，使得其線性穩(wěn)定性特征發(fā)生有利于更多T-S波產(chǎn)生并快速增長的變化，并形成一種正反饋的作用而導(dǎo)致轉(zhuǎn)捩。

圖1 斜坡型粗糙元尾跡的紅外、油流、升華法顯示圖[18]Fig.1 Infra-red, oil and sublimation visualization of ramp roughness wake[18]

在試驗研究方面，除了引言中所提到的試驗，Tirtey等[17-18]對鉆石型、斜坡型、圓柱型及圓球形粗糙元繞流流場，采用數(shù)值計算、紅外熱圖顯示、油流顯示、升華法顯示等方法，進(jìn)行了對比研究。如圖 1所示，試驗結(jié)果顯示了粗糙元上游的反向旋轉(zhuǎn)渦對及下游尾跡流向渦情況，以及尾跡中反向旋轉(zhuǎn)渦對相互吸引與交互作用、邊界層分離與再附，并完成轉(zhuǎn)捩的情況，其中紅外熱圖結(jié)果證實(shí)了分離區(qū)的低熱流以及再附后完成轉(zhuǎn)捩。Danehy等[19]和Danehy等[20]采用一氧化氮激光誘導(dǎo)熒光(Planar laser induced fluorescence,PLIF)技術(shù)測量了馬赫數(shù)10高超聲速條件下半球形和圓柱形兩種粗糙元的尾跡流動結(jié)構(gòu)。Wheaton等[31]在普渡大學(xué)靜音風(fēng)洞馬赫數(shù)6試驗條件下測量了圓柱型粗糙元尾跡不穩(wěn)定波的發(fā)展情況。

2 試驗設(shè)備及流場條件

2.1 試驗設(shè)備及測量手段

試驗在中國空氣動力研究與發(fā)展中心超高速所Φ2 m激波風(fēng)洞(FD-14A)上開展，其由內(nèi)徑為150 mm，高壓段、低壓段長度分別為9 m和18 m的激波管和相應(yīng)的噴管、試驗段、真空箱組成，其型面噴管出口直徑為1.2 m。風(fēng)洞試驗氣體為氮?dú)?，采用氫氣或氫氣和氮?dú)饣旌蠚怏w驅(qū)動，驅(qū)動壓力目前可達(dá)50 MPa。通過更換喉道或噴管可獲得不同的來流馬赫數(shù)，通過調(diào)節(jié)高低壓段的壓力可獲得不同的雷諾數(shù)，實(shí)現(xiàn)不同的模擬環(huán)境。目前該風(fēng)洞所能模擬的馬赫數(shù)范圍是6～16，雷諾數(shù)范圍是2.1×105～6.7×107m-1，試驗段的橫截面積是2.6 m×2.6 m，試驗的有效時間為4～18 ms。

試驗所采用的測量手段有點(diǎn)式鉑薄膜熱流傳感器技術(shù)、紋影技術(shù)、溫敏熱圖技術(shù)。對邊界層流態(tài)的判斷主要依靠點(diǎn)式鉑薄膜熱流傳感器技術(shù)[32]，其依據(jù)是邊界層流態(tài)由層流轉(zhuǎn)捩為湍流后，氣體向壁面?zhèn)鳠嵩鰪?qiáng)，壁面熱流升高，而這一變化可根據(jù)單點(diǎn)熱流時變特性及沿流向熱流分布特性識別出來。采用紋影技術(shù)顯示渦流發(fā)生器尾跡流場結(jié)構(gòu)，以觀察其對邊界層的擾動情況。采用溫敏熱圖技術(shù)[32]可定性顯示渦流發(fā)生器下游流向渦形成的條帶結(jié)構(gòu)。

2.2 試驗?zāi)Ｐ图傲鲌鰲l件

采用的升力體示意圖如圖 2所示，渦流發(fā)生器安裝在模型迎風(fēng)面距頭部x/L=0.32位置。取鉆石型粗糙元中心點(diǎn)連線與模型迎風(fēng)中心線交點(diǎn)為基準(zhǔn)點(diǎn)，以表1中設(shè)計流場的參數(shù)計算了該基準(zhǔn)點(diǎn)處法向速度分布(見圖 3)，選取速度分布圖中第一個拐點(diǎn)處對應(yīng)的Y值為邊界層厚度值(4.5 mm)。以該數(shù)值為基準(zhǔn)，設(shè)計了三種渦流發(fā)生器構(gòu)型，分別為2.0 mm高度、4.5 mm高度鉆石型以及4.5 mm高度斜坡型渦流發(fā)生器(見圖 4)。

圖2 升力體示意圖Fig.2 The lifting body model

為了檢驗渦流發(fā)生器對激波風(fēng)洞中模型邊界層強(qiáng)制轉(zhuǎn)捩的效果，選取激波風(fēng)洞高馬赫數(shù)低雷諾數(shù)流場條件，以確保在無渦流發(fā)生器時模型表面邊界層為層流流態(tài)，選取的試驗流場參數(shù)如表1所示。模型攻角為10°，側(cè)滑角0°。需要說明的是，這兩個試驗流場的單位雷諾數(shù)，逼近Φ2 m激波風(fēng)洞兩個馬赫數(shù)條件下的下限單位雷諾數(shù)。當(dāng)流場單位雷諾數(shù)在1.0×107/m以上時，Φ2 m激波風(fēng)洞中模型邊界層比較容易轉(zhuǎn)捩。因此在設(shè)計試驗方案時，預(yù)計表1的設(shè)計流場能保證模型表面邊界層為層流流態(tài)，但風(fēng)洞試驗結(jié)果表明，迎風(fēng)面中心線后部區(qū)域測點(diǎn)在該流場條件下開始轉(zhuǎn)捩，邊界層處于轉(zhuǎn)捩過渡狀態(tài)。為了保證模型表面邊界層為層流流態(tài)，為了驗證渦流發(fā)生器對低雷諾數(shù)流場的強(qiáng)制轉(zhuǎn)捩效果，選取的試驗流場如表1所示的兩個高馬赫數(shù)低雷諾數(shù)流場，圖 3表明了兩個試驗流場中基準(zhǔn)點(diǎn)處法向速度分布與原設(shè)計流場相差不大。

流場參數(shù)設(shè)計流場試驗流場Ma10Ma10Ma12總壓P0/MPa11．555．6610．09溫度T0/K98213481471馬赫數(shù)Ma∞9．959．7011．82單位雷諾數(shù)Re∞/m-18．5×1062．4×1062．1×106邊界層厚度δ/mm1．012．843．12k/δ(k=4．5mm)-1．581．44k/δ(k=2．0mm)-0．700．64

圖4 三種渦流發(fā)生器Fig.4 Three kinds of vortex generators

文獻(xiàn)[2,14,33]中在Ma6～Ma10流場條件下，在進(jìn)氣道壓縮面或平板模型表面安裝了渦流發(fā)生器，其誘導(dǎo)轉(zhuǎn)捩的有效高度所對應(yīng)的k/δ值大致為0.67～1.73。對于該δ值，是以99.5%總焓來定義邊界層厚度，如表1所示基準(zhǔn)點(diǎn)處邊界層厚度值，Ma10試驗流場δ=2.84 mm，Ma12試驗流場δ=3.12 mm，4.5 mm高度渦流發(fā)生器對應(yīng)的k/δ值分別為1.58和1.44，2 mm高度渦流發(fā)生器對應(yīng)的k/δ值分別為0.70和0.64。對比上述風(fēng)洞試驗所采用的渦流發(fā)生器，考慮到本文試驗中模型外形、來流條件、模型攻角等因素的差異，本文選取的渦流發(fā)生器高度比較合適。

2.3 計算方法

為了計算中心線邊界層參數(shù)，預(yù)測試驗條件下的流場結(jié)構(gòu)和氣動熱，本文采用計算流體力學(xué)(Computational fluid dynamics,CFD)數(shù)值方法進(jìn)行了計算分析。采用的計算軟件是中國空氣動力研究與發(fā)展中心超高速所自主研發(fā)的AHENS并行計算軟件。該軟件采用有限體積方法，求解三維坐標(biāo)系下的N-S方程組[34]，氣體模型為完全氣體。在采用有限體積方法對控制方程組進(jìn)行空間離散時，通過對網(wǎng)格控制體單元內(nèi)的無黏通量與黏性通量進(jìn)行積分來實(shí)現(xiàn)空間離散。對控制方程進(jìn)行積分，結(jié)合Gauss定理，可得：

(1)

3 試驗結(jié)果與分析

試驗獲得了比較理想的結(jié)果，驗證了渦流發(fā)生器的強(qiáng)制轉(zhuǎn)捩效果。下面針對鉑薄膜熱流傳感器技術(shù)、紋影技術(shù)和溫敏熱圖技術(shù)獲得的試驗結(jié)果，分別進(jìn)行討論分析。

圖 5給出了模型迎風(fēng)面中心線測點(diǎn)在無渦流發(fā)生器、2.0 mm鉆石型、4.5 mm鉆石型、4.5 mm斜坡型渦流發(fā)生器作用下的熱流測量結(jié)果，以及無渦流發(fā)生器時層流和湍流條件下迎風(fēng)面中心線熱流計算結(jié)果，可見計算結(jié)果與試驗結(jié)果吻合較好。邊界層流態(tài)由層流轉(zhuǎn)捩到湍流，流體的動量和能量輸運(yùn)能力增強(qiáng)，體現(xiàn)在氣動熱上即是使得物面熱流明顯升高。由圖 5可知，兩個流場條件下無渦流發(fā)生器時，迎風(fēng)中心線測點(diǎn)熱流呈現(xiàn)單調(diào)降低的趨勢，邊界層流動為層流，計算結(jié)果與試驗結(jié)果相符。對于存在渦流發(fā)生器的情況，圖 5(a)中，來流馬赫數(shù)為10，在x/L=0.32之前，中心線測點(diǎn)邊界層為層流；x/L=0.32～0.55之間，測點(diǎn)受渦流發(fā)生器的影響，熱流逐漸增加，邊界層流動處在轉(zhuǎn)捩過渡區(qū)；在x/L=0.55之后，熱流趨于平穩(wěn)并逐漸降低，邊界層完全轉(zhuǎn)捩為湍流，熱流峰值為轉(zhuǎn)捩開始之前熱流值的2.9倍；在轉(zhuǎn)捩過渡區(qū)，2.0 mm鉆石型和4.5 mm鉆石型渦流發(fā)生器引起的測點(diǎn)熱流爬升速度不一致，但在邊界層完全轉(zhuǎn)捩成為湍流之后，兩種情況下的熱流值相當(dāng)。圖 5(b)中，來流馬赫數(shù)為12，在x/L=0.32之前，中心線測點(diǎn)邊界層為層流；4.5 mm高度渦流發(fā)生器引起的轉(zhuǎn)捩過渡區(qū)較長，在x/L=0.32～0.55之間；2.0 mm渦流發(fā)生器引起的轉(zhuǎn)捩過渡區(qū)較短，在x/L=0.32～0.36之間很快完成轉(zhuǎn)捩。4.5 mm高度渦流發(fā)生器在x/L=0.55之后，邊界層完全轉(zhuǎn)捩為湍流，熱流峰值為轉(zhuǎn)捩開始之前熱流值的2.9倍；2.0 mm高度渦流發(fā)生器在x/L=0.36之后，邊界層為湍流，熱流峰值為轉(zhuǎn)捩開始之前熱流值的1.67倍。相同高度的4.5 mm鉆石型和4.5 mm斜坡型渦流發(fā)生器對轉(zhuǎn)捩區(qū)熱流的影響效果基本一致，2.0 mm高度渦流發(fā)生器影響的邊界層很快完成轉(zhuǎn)捩，完成轉(zhuǎn)捩之后的熱流在三種渦流發(fā)生器情況下基本一致，說明不同高度渦流發(fā)生器雖然會對轉(zhuǎn)捩過程產(chǎn)生不同程度的影響，但轉(zhuǎn)捩完成之后的熱流基本一致。這種情況可用于在相同的風(fēng)洞自由來流參數(shù)條件下，利用渦流發(fā)生器實(shí)現(xiàn)邊界層層流流態(tài)和湍流流態(tài)的模擬，研究不同邊界層流態(tài)對模型表面熱流分布的影響。

圖5 中心線熱流分布Fig.5 Heat flux distribution of the center line

圖6給出了Ma12流場條件下，三種渦流發(fā)生器下游(x/L=0.32～0.54)尾跡影響下邊界層紋影圖像，該圖像由5000幀的高速相機(jī)拍攝獲取。圖像中可以清晰分辨尾跡對邊界層的擾動情況，該擾動就是流向渦的發(fā)生發(fā)展，并導(dǎo)致轉(zhuǎn)捩的發(fā)生及完成轉(zhuǎn)捩。紋影照片表征邊界層密度梯度變化，可認(rèn)為粗糙元上游模型表面的黑帶即是邊界層邊緣，圖中可見2.0 mm高度渦流發(fā)生器浸沒在邊界層內(nèi)，4.5 mm高度渦流發(fā)生器則略高于邊界層。從圖6可以看出，渦流發(fā)生器上下游的激波結(jié)構(gòu)，流動在渦流發(fā)生器上游分離形成分離激波，并在前端經(jīng)過壓縮形成前緣激波，經(jīng)過渦流發(fā)生器之后，流動出現(xiàn)膨脹和再壓縮過程，形成相應(yīng)的波系結(jié)構(gòu)。

圖6 渦流發(fā)生器尾跡紋影圖像Fig.6 Schlieren visualization of the wake

圖7為Ma12流場條件下，鉆石型渦流發(fā)生器尾跡的溫敏熱圖，從中可以清晰分辨渦流發(fā)生器下游流向渦生成的條帶結(jié)構(gòu)，隨后條帶結(jié)構(gòu)消失、流向渦破碎，最終完成轉(zhuǎn)捩。該溫敏熱圖顯示試驗是在另一較小尺寸模型上開展的，計算獲得的當(dāng)?shù)厮俣冗吔鐚雍穸葹?.0 mm，渦流發(fā)生器基準(zhǔn)尺度也設(shè)計為3.0 mm。雖然邊界層相關(guān)結(jié)構(gòu)并不是按模型縮比的相同比例縮小，但該圖同樣能夠定性說明強(qiáng)制轉(zhuǎn)捩的過程。圖中渦流發(fā)生器中間區(qū)域引起的條帶結(jié)構(gòu)和兩側(cè)的不一樣，兩側(cè)條帶結(jié)構(gòu)引起的熱流升高幅度更大。這是由于該模型迎風(fēng)面展向是弧形的，其邊界層厚度由中心線向外側(cè)降低，但渦流發(fā)生器高度是以中心點(diǎn)位置的邊界層厚度為基準(zhǔn)設(shè)計的，這就導(dǎo)致渦流發(fā)生器外側(cè)粗糙元高度與當(dāng)?shù)剡吔鐚雍穸炔黄ヅ?，引起了更?qiáng)的渦流尾跡。

圖7 尾跡溫敏熱圖顯示結(jié)果Fig.7 Hear flux distribution of the wake with TSP technique

4 結(jié) 論

1)在Φ2 m激波風(fēng)洞的Ma10和Ma12、單位雷諾數(shù)2×106/m量級的流場條件下，渦流發(fā)生器有效地促使模型表面邊界層由層流轉(zhuǎn)捩成湍流。由于激波風(fēng)洞模擬能力的限制，現(xiàn)有高馬赫數(shù)流場基本只能在模型表面模擬層流邊界層，也曾經(jīng)嘗試過通過增加模型壁面粗糙度的方式來實(shí)現(xiàn)強(qiáng)制轉(zhuǎn)捩，但效果不理想，在粗糙帶下游一定距離，粗糙帶擾動耗散之后，邊界層恢復(fù)到層流流態(tài)。本文渦流發(fā)生器的強(qiáng)制轉(zhuǎn)捩能力，可使激波風(fēng)洞具備高馬赫數(shù)條件下模擬湍流邊界層的試驗?zāi)芰Α?/p>

2)本文試驗結(jié)果驗證了文獻(xiàn)[24]中的觀點(diǎn)，渦流發(fā)生器的外形對轉(zhuǎn)捩效果影響不大。本文試驗結(jié)果更進(jìn)一步表明不同形狀不同高度渦流發(fā)生器對邊界層完全轉(zhuǎn)捩成湍流后的熱流影響不明顯，說明強(qiáng)制轉(zhuǎn)捩完成之后，相關(guān)參數(shù)主要受邊界層流動特征影響，渦流發(fā)生器粗糙元擾動的影響不明顯。

3)在激波風(fēng)洞等地面模擬設(shè)備中，模擬層流和湍流邊界層的風(fēng)洞自由來流參數(shù)不一樣，較難開展不同流態(tài)對模型表面參數(shù)影響的試驗研究工作。本文研究結(jié)果表明，可提出一種新的風(fēng)洞模擬試驗方法，利用渦流發(fā)生器在相同風(fēng)洞來流參數(shù)條件下，實(shí)現(xiàn)邊界層層流和湍流流態(tài)的模擬，以利于開展不同邊界層流態(tài)對模型表面流動特性影響的試驗研究。

[1] 戰(zhàn)培國, 胥家常. 吸氣式高超聲速飛行器邊界層控制研究概述[J]. 試驗流體力學(xué), 2008, 22(1): 95-100. [Zhan Pei-guo, Xu Jia-chang. The review of boundary layer control for hypersonic air-breathing vehicle[J]. Journal of Experiments in Fluid Mechanics, 2008, 22(1): 95-100.]

[2] Berry S A, Auslender A H, Dilley A D, et al. Hypersonic boundary-layer trip development for Hyper-X[J]. Journal of Spacecraft and Rockets, 2001, 38(6): 853-864.

[3] Zuo S H, Yang Y, Li D, et al. Experimental investigation of cross-flow instability in swept-wing boundary layers[J]. Acta Automatica Sinica, 2010, 28(5): 495-502.

[4] Keisuke F J. An experiment of two dimensional roughness effect on hypersonic boundary-layer transition[C]. The 43rd AIAA Aerospace Sciences Meeting and Exhibit, Reno, USA, January 10-13, 2005.

[5] Stetson K F. Comments on hypersonic boundary layer transition[R]. AD/A227 242, 1990.

[6] Berry S A, Bouslog S A, Brauckmann G J, et al. Boundary layer transition due to isolated roughness-shuttle results from the LaRC 20-inch Mach 6 tunnel[C]. The 35th AIAA Aerospace Sciences Meeting and Exhibit, Reno, USA, January 6-9, 1997.

[7] Balakumar P. Transition in a supersonic boundary layer due to roughness and acoustic disturbances[C]. The 33rd AIAA Fluid Dynamics Conference and Exhibit, Orlando, USA, June 23-26, 2003.

[8] Balakumar P. Boundary layer receptivity due to roughness and freestream sound for supersonic flows over axisymmetric cones[C]. The 38th AIAA Fluid Dynamics Conference and Exhibit, Seattle, USA, June 23-26, 2008.

[9] Estruch-S D, Vanstone L, Ganapathisubramani B, et al. Effect of roughness-induced disturbances on an axisymmetric hypersonic laminar boundary layer[C]. The 51st AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Grapevine, USA, January 7-10, 2013.

[10] Borg M P, Schneider S P, Juliano T J. Effect of freestream noise on roughness-induced transition for the X-51A forebody[C]. The 46th AIAA Aerospace Sciences Meeting and Exhibit, Reno, USA, January 7-10, 2008.

[11] Schneider S P. Effects of roughness on hypersonic boundary-layer transition[C]. The 45th AIAA Aerospace Sciences Meeting and Exhibit, Reno, USA, January 8-11, 2007.

[12] 戰(zhàn)培國. 超燃沖壓發(fā)動機(jī)前體邊界層轉(zhuǎn)捩風(fēng)洞試驗方法[J]. 航空科學(xué)技術(shù), 2012(6): 22-26. [Zhan Pei-guo. Wind tunnel test methods for boundary layer transition on scramjet engine forebody[J]. Aeronautical Science & Technology, 2012(6): 22-26.]

[13] 趙慧勇, 易淼榮. 高超聲速進(jìn)氣道強(qiáng)制轉(zhuǎn)捩裝置設(shè)計綜述[J]. 空氣動力學(xué)學(xué)報, 2014, 32(5): 623-627. [Zhao Hui-yong, Yi Miao-rong. Review of design for forced-transition trip of hypersonic inlet[J]. Acta Aerodynamica Sinica, 2014, 32(5): 623-627. ]

[14] 趙慧勇, 周瑜, 倪鴻禮, 等. 高超聲速進(jìn)氣道邊界層強(qiáng)制轉(zhuǎn)捩試驗[J]. 實(shí)驗流體力學(xué), 2012, 26(1):1-6. [Zhao Hui-yong, Zhou Yu, Ni Hong-li, et al. Test of forced boundary-layer transition on hypersonic inlet[J]. Journal of Experiments in Fluid Mechanics, 2012, 26(1): 1-6.]

[15] 趙俊波, 沈清, 張紅軍, 等. 基于T-S波諧頻共振的超燃進(jìn)氣道邊界層轉(zhuǎn)捩[J]. 航空動力學(xué)報, 2010, 25(11): 2420-2424. [Zhao Jun-bo, Shen Qing, Zhang Hong-jun, et al. Boundary transition research of scramjet inlet based on the Tollimien-Schiliting(T-S) wave syntony[J].Journal of Aerospace Power, 2010, 25(11): 2420-2424.]

[16] Schneider S P. Effects of roughness on hypersonic boundary-layer transition[J]. Journal of Spacecraft and Rockets, 2008, 45(2): 193-209.

[17] Tirtey S C, Chazot O, Walpot L. Characterization of hypersonic roughness-induced boundary-layer transition[J]. Experiments in Fluids, 2011, 50(2): 407-418.

[18] Tirtey S C, Chazot O, Walpot L. Characterization of hypersonic roughness induced transition for the EXPERT Flight Experiment[C]. The 16th AIAA/DLR/ DGLR International Space Planes and Hypersonic Systems and Technologies Conference, Bremen, Germany, October 19-22, 2009.

[19] Danehy P M, Ivey C B, Inman J A, et al. High-speed PLIF imaging of hypersonic transition over discrete cylindrical roughness[C]. The 48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Orlando, USA, January 4-7, 2010.

[20] Danehy P M, Bathel B F, Ivey C B, et al. NO PLIF study of hypersonic transition over a discrete hemispherical roughness element[C]. The 47th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Orlando, USA, January 5-8, 2009.

[21] Redford J A, Sandham N D, Roberts G T. Compressibility effects on boundary-layer transition induced by an isolated roughness element[J]. AIAA Journal, 2010, 48: 2818-2830.

[22] Wang X W, Zhong X L. Receptivity of a hypersonic flat-plate boundary layer to three-dimensional surface roughness[J]. Journal of Spacecraft and Rockets, 2008, 45(6): 1165-1175.

[23] Merkle C L, Kubota T. An analytical study of the effects of surface roughness on boundary layer transition[R]. AD/A004 786, 1974.

[24] Van den Dynde J P, Sandham N D. Numerical simulation of rouhgness induced instability growth and transition at Mach 6[C]. The 7th AIAA Theoretical Fluid Mechanics Conference, Atlanta, USA, June 16-20, 2014.

[25] 趙云飛, 劉 偉, 岡敦殿, 等. 粗糙物面引起的超聲速邊界層轉(zhuǎn)捩現(xiàn)象研究[J]. 宇航學(xué)報, 2015, 36(6): 739-746. [Zhao Yun-fei, Liu Wei, Gang Dun-dian, et al. Study of surface roughness induced supersonic boundary layer transition[J]. Journal of Astronautics, 2015, 36(6): 739-746.]

[26] 段志偉, 肖志祥. 粗糙元誘導(dǎo)的高超聲速邊界層轉(zhuǎn)捩[J]. 航空學(xué)報, 2016, 37(8): 2454-2463. [Duan Zhi-wei, Xiao Zhi-xiang. Roughness element induced hypersonic boundary layer transition[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(8): 2454- 2463.]

[27] 段志偉. 直接數(shù)值模擬粗糙單元誘導(dǎo)的高超聲速邊界層轉(zhuǎn)捩[D]. 北京: 清華大學(xué), 2014. [Duan Zhi-wei. Direct numerical simulation of rougness-induced transition in hypersonic boundary layer[D]. Beijing: Tsinghua University, 2014.]

[28] 朱德華, 袁湘江, 沈清, 等. 高超聲速粗糙元誘導(dǎo)轉(zhuǎn)捩的數(shù)值模擬及機(jī)理分析[J]. 力學(xué)學(xué)報, 2015, 47(3): 381-388. [ Zhu De-hua, Yuan Xiang-jiang, Shen Qing, et al. Numerical simulation and mechanism analysis of hypersonic roughness induced transition[J]. Chinese Journal of Theoretical and Applied Mechanics, 2015, 47(3): 381-388.]

[29] 趙曉慧, 鄧小兵, 毛枚良, 等. 高超聲速進(jìn)氣道強(qiáng)制轉(zhuǎn)捩流動的大渦模擬[J]. 航空學(xué)報, 2016, 37(8): 2445-2453. [Zhao Xiao-hui, Deng Xiao-bing, Mao Mei-liang, et al. Large eddy simulation for forced transition flow at hypersonic inlet[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(8): 2445-2453.]

[30] 周恒, 蘇彩虹, 張永明. 超聲速/高超聲速邊界層的轉(zhuǎn)捩機(jī)理及預(yù)測[M]. 北京: 科學(xué)出版社, 2015.

[31] Wheaton B M, Schneider S P. Hypersonic boundary-layer instabilities due to near-critical roughness[J]. Journal of Spacecraft and Rockets, 2014, 51: 327-342．

[32] 張扣立, 常雨, 孔榮宗, 等. 溫敏漆技術(shù)及其在邊界層轉(zhuǎn)捩測量中的應(yīng)用[J]. 宇航學(xué)報, 2013, 34(6): 860-865. [Zhang Kou-li, Chang Yu, Kong Rong-zong, et al. Temperature sensitive paint technique and its application in measurement of boundary layer transition[J]. Journal of Astronautics, 2013, 34(6): 860- 865.]

[33] Wise D J, Smart M K. Forced transition of hypervelocity boundary layers[C]. The 18th AIAA/3AF International Space Planes and Hypersonic Systems and Technologies Conference, Tours, France, September 24-28, 2012.

[34] 傅德薰, 馬延文. 計算流體力學(xué)[M].北京: 高等教育出版社, 2002.

通信地址：四川綿陽高新區(qū)二環(huán)路南段6號15信箱2分箱(621000)

電話：(0816)2465269

E-mail:409195707@qq.com

Experimental Investigation on Forced Boundary-Layer Transition in Shock Tunnel

LI Qiang, ZHANG Kou-li, ZHUANG Yu, ZHAO Jin-shan

(Hypervelocity Aerodynamics Institute, China Aerodynamics Research and Development Center, Mianyang 621000, China)

The experiment of the hypersonic boundary-layer forced-transition is carried out in aΦ2 m shock tunnel which is affiliated to China Aerodynamics Research and Development Center (CARDC). The forced-transition trip is a vortex generator which is designed and mounted on a lifting body model. The test nominal Mach numbers are 10 and 12, the unit Reynolds numbers are 2.4×106/m and 2.1×106/m respectively, and the angle of attack is 10°. The model surface heat transfer is measured by the thin film thermal gauges and temperature sensitive paint (TSP) technique. Tests results show that the initial motivation of the artificial transition is triggered. Thus, it is able to carry out the measurement campaigns which require turbulent boundary layer in high Mach number and low Reynolds number flows. A series of vortex generators, differing from shape and height, are investigated in this project. By comparing heat transfer downstream of these vortex generators, no significant difference is found in regions where artificial transition is totally completed. And also, the results give us confidence to carry out the experimental study aiming to find out some different flow characteristics such as heat transfer distribution under laminar and turbulent boundary layer in the identical inflow condition.

Shock tunnel; Vortex generator; Forced transition; Turbulent flow simulation

2017-01-22;

2017-05-11

國家重點(diǎn)研發(fā)計劃“大科學(xué)裝置前沿研究”重點(diǎn)專項(2016YFA0401201)

V211

A

1000-1328(2017)07-0758-08

10.3873/j.issn.1000-1328.2017.07.000

李 強(qiáng)(1982-)，男，助理研究員，主要從事高超聲速氣動熱環(huán)境試驗、高超聲速邊界層轉(zhuǎn)捩試驗研究。