雙粗糙元對高超音速邊界層感受性的影響

史明方，徐立丹，唐小軍，王振清

（1.內(nèi)蒙古科技大學(xué) 土木工程學(xué)院，內(nèi)蒙古 包頭 014000；2.哈爾濱工程大學(xué) 航天與建筑工程學(xué)院，黑龍江 哈爾濱 150001；3.中國空間技術(shù)研究院 北京衛(wèi)星制造廠有限公司，北京 100094）

高超音速武器以其優(yōu)越的突防性能和損傷能力引起各國的重點關(guān)注［1］。在高超音速武器的設(shè)計過程中，高超音速邊界層流動狀態(tài)會直接影響高超音速武器的性能［2］，因為邊界層層流轉(zhuǎn)捩湍流會導(dǎo)致邊界層內(nèi)摩擦阻力和壁面加熱顯著增加。鑒于此，高超音速邊界層的轉(zhuǎn)捩控制是高超音速武器設(shè)計的關(guān)鍵技術(shù)之一。感受性是轉(zhuǎn)捩的第1 個階段，定義了邊界層內(nèi)擾動波的初始模態(tài)和幅值等參數(shù)。感受性階段受到很多因素的影響，簡單的可以分為來流擾動和壁面擾動［3］。其中來流擾動對感受性的影響研究較早，并已經(jīng)取得了一系列的成果。

壁面粗糙元作為壁面擾動的重要組成部分，成為近些年的研究熱點。在高超音速飛行器表面不可避免地存在壁面粗糙元，而壁面粗糙元又會對高超音速流場及邊界層造成顯著的影響，因此開展壁面粗糙元對邊界層感受性作用機理的研究是十分必要的。董昊［4］總結(jié)了壁面粗糙元對高超音速邊界層轉(zhuǎn)捩的影響，指出壁面粗糙元對高超音速邊界層感受性的作用機理十分復(fù)雜，包括壁面粗糙元和邊界層內(nèi)擾動波的相互作用，涉及瞬態(tài)增長和尾跡等現(xiàn)象。Saikishan 等［5］研究了壁面粗糙元造成的壓力梯度對邊界層轉(zhuǎn)捩的影響。Clark 等［6］利用開源軟件SU2 優(yōu)化了渦流發(fā)生器的形狀，以達(dá)到抑制邊界層內(nèi)Mack 模態(tài)演化的目的。壁面粗糙元可以依據(jù)相對光滑壁面的位置進(jìn)行分類，在本文中，將高于光滑壁面的粗糙元稱為K型粗糙元，相應(yīng)的，將低于光滑壁面的粗糙元稱為D 型粗糙元，K 型粗糙元的研究相對較早。段志偉等［7］分析了壁面粗糙元形狀對轉(zhuǎn)捩模式的影響。Ryota 等［8］結(jié)合fhf 直接數(shù)值模擬（direct numerical simulation，DNS）和非線性拋物化穩(wěn)定性方程（nolinear parabolizd stability equations，NPSE）方法分析了三維圓柱狀粗糙元高度對邊界層轉(zhuǎn)捩的影響，并討論了邊界層轉(zhuǎn)捩時壁面粗糙元的臨界高度。Vignes 等［9］分析了渦流發(fā)生器對圓柱形粗糙元周圍高超音速流場氣動特性的影響，發(fā)現(xiàn)渦流發(fā)生器有助于降低圓柱形粗糙元造成的氣動阻力。D型粗糙元（主要集中在矩形腔結(jié)構(gòu)）對高超音速邊界層感受性的影響成為研究熱點。李闖等［10］采用直接數(shù)值模擬方法討論了空腔結(jié)構(gòu)對6.0 Ma 高超音速邊界層內(nèi)擾動波演化的作用機理，發(fā)現(xiàn)空腔深度的增加會改變邊界層內(nèi)第2 模態(tài)的演化機理。劉俊等［11］系統(tǒng)地研究了空腔結(jié)構(gòu)內(nèi)波系的傳播規(guī)律，討論了壓力振蕩反饋模型的異同。李強等［12］利用激波風(fēng)洞分析了展向凹槽對高超音速平板邊界層轉(zhuǎn)捩的影響，發(fā)現(xiàn)凹槽會促進(jìn)轉(zhuǎn)捩，并且隨著凹槽深度和寬度的增加而增強。上述學(xué)者所分析的壁面粗糙元均為獨立粗糙元，與之相對應(yīng)的分布式粗糙元的研究在獨立粗糙元的基礎(chǔ)上也在逐步開展。董昊等［13］利用數(shù)值模擬和油膜干涉技術(shù)分析了壁面粗糙元高度和間距對轉(zhuǎn)捩位置的影響，發(fā)現(xiàn)高度的影響要大于間距的影響。Irimpan 等［14］利用分布式粗糙元來表征燒蝕后的熱防護(hù)系統(tǒng)，發(fā)現(xiàn)分布式粗糙元可以降低壁面加熱率17%～42%。Zhou 等［15］指出波浪形壁面可以有效地抑制邊界層內(nèi)第2模態(tài)不穩(wěn)定性，并且抑制效應(yīng)隨波浪形壁面的延長、加深和向下游移動而增強。Tong等［16-17］研究了正弦狀壁面粗糙元對高超音速邊界層轉(zhuǎn)捩的影響。

目前，獨立粗糙元和分布式粗糙元的研究都已取得了一定成果，但雙粗糙元對高超音速邊界層感受性影響的研究相對較少，并且大部分研究集中在分析雙粗糙元間距的影響，缺乏雙粗糙元組合類型影響的研究。本文對高超音速鈍錐繞流流場進(jìn)行直接數(shù)值模擬，討論了雙粗糙元組合類型對高超音速流場及邊界層的影響；向流場中引入有效振幅脈沖慢聲波，借助快速傅里葉變換方法分析了雙粗糙元組合類型對脈沖慢聲波作用下邊界層內(nèi)擾動波傳播過程的影響，以揭示雙粗糙元組合類型對高超音速邊界層感受性的作用機理。

1 控制方程及數(shù)值方法

為了分析雙粗糙元對脈沖有限振幅自由來流擾動作用下高超音速流場及邊界層內(nèi)流動及邊界層的影響，本文采用直接數(shù)值模擬方法對高超音速流場進(jìn)行直接數(shù)值模擬，揭示雙粗糙元的組合類型對高超音速邊界層感受性的影響。

對高超音速流場的數(shù)值模擬，控制方程為二維N-S方程：

式中：U為解向量；F和G為對流項；Eυ與Fυ為黏性項。解向量的具體表達(dá)式為：

對流項的表達(dá)式為：

黏性項的表達(dá)式為：

式中：ρ、u、v、T和p分別為密度、x方向速度、y方向速度、溫度和壓強；導(dǎo)熱系數(shù)k可以通過假設(shè)普朗特數(shù)為常數(shù)獲得；E為單位流體的總能，可以表示為單位流體的內(nèi)能e與動能的和，具體為：

對于量熱氣體，氣體內(nèi)能e可計算為：

式中：比熱比γ=1.4；τ為切應(yīng)力，一般將空氣看作牛頓流體：

式中μ為黏性系數(shù)，可以通過Sutherland公式獲得：

式中：T0和μ0分別為參考溫度和黏性系數(shù)；Ts=110.4 K。

本文考慮空氣為完全氣體，補充狀態(tài)方程為：

式中：R為通用氣體常數(shù)，R=286.94 J/(kg·K)。此外，本文將笛卡爾坐標(biāo)系下控制方程進(jìn)行坐標(biāo)變換，轉(zhuǎn)換到曲線坐標(biāo)系下，以便于對復(fù)雜外形下的流場進(jìn)行數(shù)值求解。

為了分析雙壁面粗糙元組合形式對高超音速邊界層感受性的影響，本文以自由來流擾動作用下高超音速鈍錐繞流流場為研究對象，鈍錐鼻部半徑r=1 mm，半錐角θ=5°。本節(jié)以K型粗糙元（上游）和D型粗糙元（下游）的組合為例，給出了本文的計算模型及網(wǎng)格，計算網(wǎng)格分布為801×251，計算模型如圖1 所示。K 型粗糙元（上游）和D 型粗糙元（下游）的雙壁面粗糙元組合形式表示為Kf Db，下標(biāo)f和b分別代表上游粗糙元和下游粗糙元，其余3 種組合形式以此類推。壁面粗糙元形狀方程（以弧線ABCDE為例）為：

圖1 計算模型及網(wǎng)格Fig.1 Computing model and grid

點A和點E分別為壁面粗糙元的起點和終點；點B和點D是壁面粗糙元形狀的拐點(xⅠf或xⅠb和xⅡf或xⅡb)，分別稱為第1拐點(xⅠf或xⅠb)和第2拐點(xⅡf或xⅡb)；點C為壁面粗糙元中點(xcf或xcb)。需要說明的是，K型粗糙元和D型粗糙元在同一位置處（壁面粗糙元中點相同）關(guān)于光滑壁面對稱。計算模型的幾何參數(shù)在表1中給出。本文使用的變量將分別利用頭部半徑r、來流參數(shù)ρ∞、u∞、T∞、p∞，u∞/r，r/u∞、k∞和μ∞對長度、密度、速度、溫度、壓強、頻率、時間、導(dǎo)熱系數(shù)和黏性系數(shù)進(jìn)行無量綱化處理，自由來流條件在表2中給出。

表1 壁面粗糙元幾何參數(shù)Table 1 Geometric parameters of roughness element mm

表2 自由來流條件Table 2 Conditions of freestream

自由來流參數(shù)，如擾動類型、持續(xù)時間等，會顯著地改變高超音速邊界層感受機理［18-19］。較大幅值的自由來流擾動與邊界層的作用更加強烈，導(dǎo)致感受性階段更復(fù)雜，甚至改變轉(zhuǎn)捩路徑。擾動波的作用時間也會顯著地影響邊界層感受性機理，短時間高強度的自由來流擾動又稱為有限振幅脈沖擾動，對高超音速邊界層感受性的影響機理與長時間低強度的自由來流擾動存在很大的差異。本文在分析壁面粗糙元對高超音速邊界層感受性的影響時，向流場中引入脈沖慢聲波擾動，其表達(dá)式為：

式中：ε、k和f分別為擾動波的幅值、波數(shù)和頻率，ε=6×10-2,k=3.144 6×10-4和f=50π，因此本文自由來流擾動的周期為4。本文中擾動量是指非穩(wěn)態(tài)高超音速流場的瞬時值與穩(wěn)態(tài)流場值的差值，擾動量利用上標(biāo)“'”表示，如溫度擾動值利用“T'”表示。在分析壁面粗糙元對高超音速邊界層感受性的作用機理時向高超音速流場引入脈沖擾動波，即向流場中引入半個周期(Δt=2)的脈沖慢聲波擾動。

高階精度有限差分方法可以穩(wěn)定準(zhǔn)確地對高超音速流場進(jìn)行數(shù)值模擬。本文基于N-S方程采用S-W（steger-warming）通量矢量分裂方法將對流項分裂為正、負(fù)2 個流通量矢量。利用五階加權(quán)本質(zhì)無振蕩（weighted essentially non-scillatory，WENO），格式對上述正、負(fù)流通量矢量進(jìn)行離散；對于粘性項采用六階精度中心差分格式進(jìn)行離散；時間項利用三步三階TVD 型Runge-Kutta 格式進(jìn)行推進(jìn)。本文采用的數(shù)值方法已經(jīng)廣泛使用在高超音速流場及邊界層的數(shù)值模擬中［20-21］，并已對數(shù)值方法和網(wǎng)格獨立性進(jìn)行過驗證，結(jié)果表明計算過程具有良好的穩(wěn)定性，并可以保證計算結(jié)果的準(zhǔn)確性。

2 雙粗糙元對穩(wěn)態(tài)高超音速流場的影響

壁面粗糙元會改變高超音速邊界層感受性機理，為揭示雙粗糙元對高超音速邊界層感受性的影響，首先需對雙粗糙元作用下穩(wěn)態(tài)高超音速流場進(jìn)行分析。圖2 為壁面粗糙元作用下高超音速流場流向速度分布，可以看出壁面粗糙元顯著地改變了流場結(jié)構(gòu)，在流場中形成壓縮波和膨脹波。獨立K型粗糙元導(dǎo)致流場中形成壓縮波-膨脹波-壓縮波波系；而獨立D 型粗糙元導(dǎo)致流場中形成膨脹波-壓縮波-膨脹波波系。雙粗糙元會導(dǎo)致流場中形成更復(fù)雜的波系，上游粗糙元造成的壓縮波或膨脹波會與下游粗糙元造成的壓縮波或膨脹波相“融合”，例如KfKb型粗糙元對導(dǎo)致流場中形成壓縮波-膨脹波-壓縮波-膨脹波-壓縮波波系，Kf粗糙元造成的下游壓縮波和Kb粗糙元造成的上游壓縮波發(fā)生了融合；相似的，KfDb型粗糙元則在流場中形成了壓縮波-膨脹波-壓縮波-膨脹波波系。需要指出的是K型粗糙元形成的波系結(jié)構(gòu)與Xue等［22］分析在前臺階對高超音速流場的影響是相似的，同時也指出前臺階對高超音速流動的影響局限在一定范圍內(nèi)。

圖2 穩(wěn)定高超音速流場流向速度云圖Fig.2 Contours of velocity along x-axis in syeady hypersonic flowfield

壁面粗糙元對高超音速流場結(jié)構(gòu)產(chǎn)生了一定影響，壁面粗糙元也會對高超音速邊界層內(nèi)流動狀態(tài)的影響，圖3 給出了壁面雙粗糙元對穩(wěn)態(tài)高超音速邊界層內(nèi)壁面壓力的影響?？梢钥闯霰诿娲植谠獙Υ植谠浇诿鎵毫Ψ植籍a(chǎn)生了顯著的影響。以獨立K 型粗糙元為例，壁面壓力在K 型粗糙元前緣區(qū)域顯著增加，在第1拐點處達(dá)到最大值，之后壁面壓力沿流向方向逐漸減小至第2 拐點，在粗糙元中心恢復(fù)到光滑壁面水平。在第2拐點至壁面粗糙元后緣區(qū)域，壁面壓力逐漸增加至光滑壁面狀態(tài)。壁面壓力的變化情況是由壁面曲率變化導(dǎo)致的，本文所采用的壁面粗糙元由3 次多項式控制，并且壁面粗糙元中點處與光滑壁面平行，因此在壁面曲率增大的區(qū)域壁面壓力增加，并在壁面粗糙元中心處恢復(fù)到光滑壁面狀態(tài)。D型粗糙元對壁面壓力的影響情況相似。此外，壁面粗糙元的影響范圍有限，僅局限在壁面粗糙元附近，并且K 型粗糙元的影響范圍要大于D 型粗糙元。K 型粗糙元的出現(xiàn)導(dǎo)致其上游和下游區(qū)域的壓強分布產(chǎn)生了變化，該現(xiàn)象在D 型粗糙元中并不明顯。對于雙粗糙元，雖然壁面粗糙元的影響范圍有限，但本文中下游粗糙元均位于上游粗糙元的影響范圍之內(nèi)，因此在K 型粗糙元下游的壁面粗糙元區(qū)域的壁面壓強均發(fā)生了一定程度的減小，而在D 型粗糙元下游的壁面粗糙元區(qū)域的壁面壓強均發(fā)生了一定程度的增長。DfDb型粗糙元在除兩粗糙元之間的區(qū)域內(nèi)，壁面壓強分布與獨立粗糙元相同。

圖3 高超音速邊界層內(nèi)壁面壓力比較Fig.3 Comparisons of wall pressure in hypersonic bound‐ary layer

壁面摩擦系數(shù)表征了邊界層內(nèi)的剪切狀態(tài)，并且壁面摩擦力會顯著地影響高超音速飛行器的性能。鑒于此，對邊界層內(nèi)壁面摩擦系數(shù)進(jìn)行分析是有必要的。圖4給出了邊界層內(nèi)壁面摩擦系數(shù)沿流向的分布。計算結(jié)果表明壁面粗糙元顯著地改變了邊界層內(nèi)的流動狀態(tài)，在K型粗糙元的中心處，壁面摩擦系數(shù)到達(dá)極大值。這是因為對于穩(wěn)態(tài)流場來說，邊界層外邊界的速度是相同的，因此壁面摩擦系數(shù)與邊界層厚度呈反比。在邊界層內(nèi)粘性的作用下，K型粗糙元中心處邊界層的厚度最薄，壁面摩擦系數(shù)達(dá)到極大值。同理，壁面摩擦系數(shù)在D型粗糙元的前緣和后遠(yuǎn)處達(dá)到極大值，在D型粗糙元中心處達(dá)到了極小值。當(dāng)壁面摩擦系數(shù)小于零時，說明在區(qū)域內(nèi)發(fā)生了回流，產(chǎn)生了流動分離。在本文中，獨立K型粗糙元并未造成流動分離，而獨立D型粗糙元均導(dǎo)致邊界層流動分離。這說明D型粗糙元更容易導(dǎo)致流動的分離。在雙粗糙元中包含D型粗糙元時，在D型粗糙元區(qū)域內(nèi)邊界層都發(fā)生了流動分離。

圖4 高超音速邊界層內(nèi)壁面摩擦系數(shù)比較Fig.4 Comparisons of friction coefficient in hypersonic boundary layer

此外，KfKb型雙粗糙元也造成了流動分離，流動分離發(fā)生在雙粗糙元之間。因為在Kf粗糙元后緣下游至Kb粗糙元上游區(qū)域內(nèi)，壓力顯著增加，逆壓梯度很大，導(dǎo)致流動在該區(qū)域內(nèi)發(fā)生分離。流動分離現(xiàn)象會顯著地改變邊界層內(nèi)的流動結(jié)構(gòu)，進(jìn)而影響到邊界層內(nèi)的壁面摩擦力和熱流密度。獨立粗糙元導(dǎo)致壁面摩擦系數(shù)沿流向的減小值要超過增加值，也就是說壁面粗糙元可以起到減小阻力的作用。也有學(xué)者得出相似的結(jié)論，Brian等［23］指出分布式粗糙元可以在保持邊界層內(nèi)摩擦阻力基本不變的情況下穩(wěn)定邊界層內(nèi)擾動波。因此合理地布置壁面粗糙元可以在一定程度上達(dá)到減阻，甚至降熱的目的。結(jié)合Zhou 等［15］指出壁面粗糙元可以穩(wěn)定邊界層的作用，壁面粗糙元的合理布置對高超音速飛行器的設(shè)計具有一定的指導(dǎo)意義。因此，對雙粗糙元作用下高超音速邊界層感受性機理的研究是有意義的。

3 雙粗糙元對非穩(wěn)態(tài)高超音速流場的影響

為了分析雙粗糙元對高超音速邊界層感受性的影響，向高超音速流場中引入脈沖慢聲波擾動。圖5 以KfKb型粗糙元為例給出了脈沖擾動波在高超音速流場中的傳播過程。

圖5 脈沖擾動作用下高超音速溫度擾動流場Fig.5 Contours of transient temperature fluctuation at different time

高超音速流場可以分為4 個區(qū)域：R1區(qū)域-激波外區(qū)域；R2區(qū)域-激波區(qū)域；R3區(qū)域-邊界層與激波之間的區(qū)域；R4區(qū)域-邊界層區(qū)域，如圖5(h)所示。脈沖擾動波與弓形激波的作用更加強烈，在脈沖慢聲波經(jīng)過弓形激波后會在高超音速流場中形成“拖拽波”，拖拽波的形成是因為脈沖慢聲波在經(jīng)過激波后的速度不同。脈沖慢聲波在與斜激波相互作用時，并未在流場中形成明顯的拖拽波，而是在流場中依次形成低溫-高溫結(jié)構(gòu)。壁面粗糙元影響了脈沖慢聲波在高超音速流場中R3和R4區(qū)域的傳播。

第2 節(jié)已經(jīng)指出壁面粗糙元改變了穩(wěn)態(tài)高超音速流場結(jié)構(gòu)，在流場中形成一組新的波系，并導(dǎo)致邊界層發(fā)生流動分離?？梢钥闯鲭p粗糙元造成的壓縮波和膨脹波波系改變了R3區(qū)域內(nèi)擾動波的傳播，在流場中明顯的看出3道壓縮波對流場內(nèi)擾動波的放大作用，如圖5(k)所示。自由來流擾動在經(jīng)過KfKb型粗糙元造成的壓縮波時被放大，并在離開這3道壓縮波時恢復(fù)到光滑壁面狀態(tài)。脈沖慢聲波在經(jīng)過激波后進(jìn)入到邊界層內(nèi)，激發(fā)出邊界層內(nèi)擾動波，導(dǎo)致邊界層內(nèi)依次形成了低溫斑和高溫斑結(jié)構(gòu)，如圖5(g)中矩形區(qū)域所示。邊界層外邊界處形成“繩索”狀高溫斑-低溫斑結(jié)構(gòu)，如圖5(g)中橢圓形區(qū)域所示，Ma等［24-26］指出這是邊界層內(nèi)慢聲波結(jié)構(gòu)的典型結(jié)構(gòu)。邊界層內(nèi)擾動波在向下游傳播的過程中逐漸耗散，但在邊界層內(nèi)形成殘余反射波，反射波在邊界層內(nèi)進(jìn)行多次反射，最終在邊界層內(nèi)粘性作用下耗散掉，如圖5（l）所示。

為了更清晰地揭示壁面粗糙元對高超音速邊界層感受性的影響，圖6給出了邊界層內(nèi)不同位置處壁面壓強擾動值隨時間的變化情況。邊界層內(nèi)越下游區(qū)域越晚感受到邊界層內(nèi)擾動波，擾動波在向下游傳播過程中，持續(xù)時間越長，在x=0.90 處，擾動波持續(xù)時間Δt=3.8；而在x=18.89 處，擾動波持續(xù)時間Δt=17.9。邊界層內(nèi)擾動波持續(xù)時間都超過了自由來流擾動波的時間(Δt=0.5T=2)。說明邊界層內(nèi)反射波在邊界層內(nèi)多次反射后才逐漸耗散消失。此外，擾動波向下游傳播過程中幅值逐漸減小，擾動波波形越復(fù)雜。邊界層內(nèi)擾動波在進(jìn)入到高超音速流場后呈現(xiàn)“正弦波”分布，在向下游傳播的過程中，在粘性的作用下，擾動波的幅值逐漸減小，波形逐漸變得復(fù)雜。邊界層內(nèi)擾動波在通過第1粗糙元(Kf或Df)時，邊界層內(nèi)擾動波在通過Kf粗糙元前緣(x=0.90)處被放大，而Df粗糙元前緣處壁面壓強擾動并未受到影響，因為D 型粗糙元的影響范圍較小。在Kf粗糙元中點(x=1.39)處，邊界層內(nèi)擾動波稍被放大。邊界層擾動波在通過第1粗糙元后緣(x=1.79)處，Kf、Kf Kb和Kf Db粗糙元降低了該處邊界層擾動波的幅值，而Df、Df Db和Df Kb粗糙元放大了該處擾動波的幅值。Kb粗糙元也放大了擾動波的幅值，但Db粗糙元并未對該處產(chǎn)生影響。在兩粗糙元中點(x=1.91)處，Kf和Db粗糙元抑制了邊界層內(nèi)擾動波的傳播，因此，Kf Db粗糙元對邊界層內(nèi)擾動波傳播的抑制作用更明顯。同理，Df和Kb粗糙元促進(jìn)了邊界層內(nèi)擾動波的傳播，因此，Df Kb粗糙元對邊界層內(nèi)擾動波傳播的促進(jìn)作用更明顯。Kf Kb和Df Db雙粗糙元對該點擾動波的傳播基本沒影響，這可以說明同類型的上游粗糙元和下游粗糙元對兩粗糙元中點處的影響是相同的，二者的影響相互抵消。下游粗糙元對邊界層內(nèi)擾動波的作用機理與上游粗糙元相同，這里不再贅述。壁面粗糙元影響的局限性導(dǎo)致在兩粗糙元下游區(qū)域的擾動波的傳播并未受到影響。

圖6 各流向位置處壁面壓強擾動幅值-時間對比Fig.6 Time history trace of wall pressure fluctuation at various streamwise locations

為了更準(zhǔn)確地分析高超音速邊界層內(nèi)擾動波的傳播過程，本文采用快速傅里葉變換方法將邊界層內(nèi)的時域信號轉(zhuǎn)換到頻域內(nèi)，以分析邊界層內(nèi)不同模態(tài)擾動波的演化過程?？焖俑道锶~變換公式為：

式中：f為自由來流擾動波的頻率；|p'|為任意瞬時壓力擾動；|p'n|和φn分別為局部擾動幅值和相角，下標(biāo)整數(shù)n代表了擾動場的波模態(tài)，當(dāng)n=0代表平均流場變形，n=1代表基頻擾動，n=2代表二階諧頻擾動，以此類推。第2節(jié)已經(jīng)指出壁面粗糙元對高超音速流場的影響具有局限性，因此本節(jié)僅給出x=0.90～10.01范圍內(nèi)邊界層內(nèi)擾動波的演化情況，如圖7所示。可以看出邊界層內(nèi)擾動波在向下游傳播過程中擾動波的波形變得越來越復(fù)雜，擾動波的幅值在逐漸減小，在x=10.01處稍有增加。此外，邊界層內(nèi)擾動波的主導(dǎo)模態(tài)在逐漸增大，如在x=0.90處，主導(dǎo)模態(tài)f=0.033；在x=1.79處，主導(dǎo)模態(tài)f=0.028；在x=2.40處，主導(dǎo)模態(tài)f=0.025，并且主導(dǎo)模態(tài)在邊界層內(nèi)擾動波所占比例在逐漸減小。也就是說高超音速邊界層內(nèi)擾動波存在模態(tài)競爭現(xiàn)象，主導(dǎo)模態(tài)進(jìn)行了模態(tài)遷移。

圖7 邊界層內(nèi)擾動波頻率振幅Fig.7 Frequency-amplitude diagram of perturbations in boundary layer

壁面粗糙元對邊界層內(nèi)不同模態(tài)擾動波傳播過程的影響相同，即在K 型粗糙元前緣區(qū)域和D 型粗糙元的后緣區(qū)域內(nèi)，邊界層內(nèi)擾動波的幅值被放大；相應(yīng)的在D型粗糙元前緣區(qū)域和K型粗糙元的后緣區(qū)域內(nèi)，邊界層內(nèi)擾動波的傳播被抑制。值得注意的是Kf獨立粗糙元抑制了其下游區(qū)域內(nèi)擾動的演化。雙粗糙元的組合類型對邊界層內(nèi)不同模態(tài)擾動波的傳播具有一定的影響，Kf Kb型雙粗糙元導(dǎo)致壁面粗糙元下游區(qū)域內(nèi)擾動波的主導(dǎo)模態(tài)顯著降低，而且這種抑制效應(yīng)還要強于Kf獨立粗糙元。下游的Kb粗糙元增強了上游的Kf粗糙元對Kf Kb雙粗糙元下游邊界層內(nèi)擾動波的抑制效應(yīng)。這就證明了K型粗糙元對下游的邊界層內(nèi)擾動波演化的抑制具有協(xié)同效應(yīng)。

4 結(jié)論

1）雙壁面粗糙元會顯著地改變穩(wěn)態(tài)高超音速流場及邊界層內(nèi)的流動狀態(tài)，導(dǎo)致流場中形成一組新的波系，并導(dǎo)致邊界層內(nèi)出現(xiàn)流動分離。如Kf Kb型粗糙元會導(dǎo)致流場中形成壓縮波-膨脹波-壓縮波-膨脹波-壓縮波波系，并在兩粗糙元之間形成流動分離。

2）雙壁面粗糙元導(dǎo)致流場中形成新的波系顯著地影響了高超音速流場中自由來流擾動的傳播過程，使擾動波在壓縮波作用下發(fā)生了一定的變形；當(dāng)自由來流擾動離開壁面粗糙元造成的波系后，新的波系影響也隨之消失。

3）雙壁面粗糙元也改變了邊界層內(nèi)擾動波的演化過程，在K 型粗糙元前緣區(qū)域和D 型粗糙元的后緣區(qū)域內(nèi)，邊界層內(nèi)擾動波的幅值被放大；在D型粗糙元前緣區(qū)域和K 型粗糙元的后緣區(qū)域內(nèi)，邊界層內(nèi)擾動波的傳播被抑制。壁面粗糙元對邊界層內(nèi)擾動波的作用也局限在一定范圍內(nèi)。

4）邊界層內(nèi)擾動波在向下游傳播過程中主導(dǎo)模態(tài)在逐漸減小，并存在模態(tài)競爭現(xiàn)象。雙粗糙元的組合類型對邊界層內(nèi)不同模態(tài)擾動波的傳播具有一定的影響。其中，KfKb型雙粗糙元顯著抑制了其下游區(qū)域內(nèi)主導(dǎo)模態(tài)擾動波的演化，并且2個K型粗糙元對其下游邊界層內(nèi)擾動波演化的抑制具有協(xié)同效應(yīng)。