局部擾動對平板邊界層流動穩(wěn)定性影響的研究

陸昌根，趙玲慧，沈露予

（南京信息工程大學(xué) 大氣科學(xué)學(xué)院，江蘇 南京210044）

0 引 言

平板邊界層穩(wěn)定性問題對飛行器、船舶等動力特性有著重要的影響，人們已經(jīng)對平板邊界層流動穩(wěn)定性及轉(zhuǎn)捩問題作了大量的理論工作，但局部擾動對平板邊界層流動穩(wěn)定性影響和控制的研究較少；目前，仍有許多問題需要我們?nèi)フJ(rèn)識和解決［1－2］。Herbert等人對平板邊界層流中穩(wěn)定性問題進(jìn)行了深入的研究［3］，其某些結(jié)論已被實驗證實。影響平板邊界層流動穩(wěn)定性的因素很多，包括溫度、對流、壓力梯度、壁面條件等等，且壁面局部擾動是普遍存在的，它對流動中的流體運動必然會產(chǎn)生影響。該局部擾動必將影響到邊界層流動的改變，這必將改變擾動波在平板邊界層運動中的演化規(guī)律，誘導(dǎo)出 展向速度以及產(chǎn)生更復(fù)雜的運動，最終將導(dǎo)致橫流不穩(wěn)定等［4－6］。實驗表明，改變局部粗糙的分布及幾何形狀，可改變低速流動的層流分離點。Floryan利用連續(xù)分布的抽吸結(jié)構(gòu)來模擬壁面粗糙的特征，并對Poiseuille流、Couette流及Blasius邊界層流中的擾動波的穩(wěn)定特性進(jìn)行了數(shù)值研究，獲得了一些有意義的成果［7］。本文將在考慮前人研究的基礎(chǔ)上，在平板上設(shè)計不同形式的局部擾動結(jié)構(gòu)，數(shù)值研究該結(jié)構(gòu)產(chǎn)生的擾動對平板邊界層流動穩(wěn)定性影響和控制機理的問題。對該問題的進(jìn)一步數(shù)值研究，將會讓我們能做到有效地抑制平板邊界層轉(zhuǎn)捩過程中湍流發(fā)生的物理機制。

1 理論方程和數(shù)值方法

采用的理論方程為不可壓、無量綱的N－S方程和連續(xù)方程：

其中▽為梯度算子，▽2為拉普拉斯算子，Re為雷諾數(shù)。u、p分別為速度和壓力。

對理論方程（1）進(jìn)行直接數(shù)值模擬，其時間導(dǎo)數(shù)為三階混合顯－隱分裂格式，空間導(dǎo)數(shù)在x和y向分別采用非等間距網(wǎng)格的緊致有限差分離散和z向空間導(dǎo)數(shù)用Fourier譜展開相雜交的方法逼近，具體數(shù)值方法詳見文獻(xiàn)［8］。其中無量綱特征長度和特征速度分別為進(jìn)口邊界層厚度δ和無窮遠(yuǎn)來流速度U∞。雷諾數(shù)為：Re＝U∞·δ／υ＝3250，υ為運動粘性系數(shù)。

在x、y、z方向上的計算區(qū)域和網(wǎng)格點數(shù)分別為120×120×14.6和250×250×32。

基本流為Blasius解。入口速度邊界條件，x＝0：u＝a0［（y）］ei（αx＋βz－ωt）＋c.c；其中c.c為共扼復(fù)數(shù)；初始幅值：a0＝0.006；流向波數(shù)：α＝（0.374696，0.001682）；展向波數(shù)：β＝0.43；

頻率：ω＝0.13；［（y）］＝｛u（y），v（y），w（y）｝是由Orr－Sommerfeld方程特征值問題獲得的特征速度，其歸一化條件 max｜u（y）｜＝1；?p／?x＝0；

出口邊界條件：x＝120，速度為無反射出流邊界條件，?p／?x＝0；

邊界層外緣邊界條件為：y＝6，?u／?y＝?v／?y＝?w／?y＝0，p＝0；

壁面邊界條件為：y＝0，無滑移條件，?p／?y＝0。

光滑壁面（1）：無滑移條件；

局部抽吸（2）：v1＝0.004（網(wǎng)格x＝1～1.8，z＝－1.37～1.37），其余網(wǎng)格：無滑移條件；

局部噴出（3）：v1＝－0.004（網(wǎng)格x＝1～1.8，z＝－1.37～1.37），其余網(wǎng)格：無滑移條件；

局部噴出和局部抽吸組合（4）：v1＝－0.004（網(wǎng)格x＝1～1.8，z＝－1.37～1.37），v1＝0.004（網(wǎng)格x＝3.4～4.2，z＝－1.37～1.37），其余網(wǎng)格：無滑移條件；

局部抽吸和局部噴出組合（5）：v1＝0.004（網(wǎng)格x＝1～1.8，z＝－1.37～1.37），v1＝－0.004（網(wǎng)格x＝3.4～4.2，z＝－1.371.37），其余網(wǎng)格：無滑移條件；

局部抽吸（6）：v1＝0.004（網(wǎng)格x＝1～1.8，z＝－7.3～7.3），其余網(wǎng)格：無滑移條件；

局部噴出（7）：v1＝－0.004（網(wǎng)格x＝1～1.8，z＝－7.3～7.3），其余網(wǎng)格：無滑移條件；

在展向上設(shè)計兩個局部抽吸組合（8）：v＝0.004（網(wǎng)格x＝1～1.8，z＝－4.10～－1.37），v＝0.004（網(wǎng)格x＝1～1.8，z＝1.37～4.10），其余網(wǎng)格：無滑移條件。

2 數(shù)值結(jié)果

圖1給出了在不同局部擾動結(jié)構(gòu)情況下，三維擾動波的幅值隨雷諾數(shù)的變化。幅值定義為：

其中｛u，v，w｝T表示擾動速度，且1、2、3、4、5、6、7、8分別與上述壁面結(jié)構(gòu)情況相對應(yīng)，其在對應(yīng)上述結(jié)構(gòu)邊界層流中三維擾動波的幅值隨雷諾數(shù)演化的規(guī)律。

圖1 幅值的隨雷諾數(shù)的演化曲線Fig.1 Amplitudes with Reynolds number

與光滑壁面相比較，2、3、4、5和8壁面的局部抽吸、局部噴出以及其它組合結(jié)構(gòu)均能促使擾動波幅值的增大，其原因是上述壁面產(chǎn)生的擾動使原二維基本流變成三維基本流，并且不同的類形、強度及結(jié)構(gòu)分布各自對三維擾動波幅值演化的作用能力是不同的。展向構(gòu)建的兩個局部抽吸的壁面8促使三維擾動波的幅值增長最快，其次是流向局部噴出和局部抽吸組合的壁面4和流向局部抽吸和噴出組合的壁面5對擾動波的幅值貢獻(xiàn)率大于單個局部抽吸的壁面2。但是局部噴出的壁面3也能使擾動波幅值增長，其激勵能力相對弱一些。局部抽吸的壁面6在局部抽吸附近對幅值貢獻(xiàn)率較小，隨著擾動波向下游的演化，擾動幅值的增長迅速增強，其原因是該壁面產(chǎn)生的擾動使二維基本流比原基本流瘦小或虧損；在局部噴出的壁面7中，三維擾動波的幅值增長小于光滑壁面的三維擾動波幅值增長的演化，它對流動的演化有較好的穩(wěn)定和控制作用，其原因是該壁面產(chǎn)生的擾動使二維基本流更加飽滿。

由于局部抽吸、局部噴出以及其它組合結(jié)構(gòu)誘導(dǎo)出的擾動速度改變了原邊界層的二維基本流特性，從這些壁面結(jié)構(gòu)類型可知，除壁面6、7外（這兩種壁面設(shè)計的是二維局部擾動，未改變原基本流的二維特性），二維基本流特性巳發(fā)生改變，即基本流已具有三維特性，這可能是使得沿空間演化的三維擾動波增長率得以增長的一個重要原因。對于不同局部類型、強度大小及分布結(jié)構(gòu)特征的不同所呈現(xiàn)的增長機制不同而言。圖2所示，壁面2、3、4、5、8及壁面6的平板邊界層流中，三維擾動波幅值增長率均大于光滑壁面，壁面7的邊界層流中擾動波幅值增長率小于光滑壁面；壁面8的三維擾動波幅值增長率貢獻(xiàn)最大；其次是壁面4和壁面5，壁面3邊界層中的增長率雖然大于光滑壁1，但貢獻(xiàn)率較小。不管是什么性質(zhì)的局部壁面，只要該結(jié)構(gòu)能誘導(dǎo)產(chǎn)生三維擾動，總時加速流動的失穩(wěn)，且在局部擾動區(qū)域附近，幅值增長率加速增長；隨著局部擾動向下游傳播時，其局部擾動產(chǎn)生的三維擾動逐漸減弱，這時幅值增長率的增長趨于平緩，這是由于擾動波向下游發(fā)展過程中展向速度的減少導(dǎo)致的。壁面局部擾動（7）能控制流動的失穩(wěn)?？傊?，除局部擾動壁面（7）外，上述局部擾動與光滑壁面比較，對流動失穩(wěn)都起激勵作用，并促使轉(zhuǎn)捩的快速發(fā)生，其轉(zhuǎn)捩發(fā)生的速度與貢獻(xiàn)率的大小有關(guān)，貢獻(xiàn)率大的，失穩(wěn)越易發(fā)生；貢獻(xiàn)率小的，失穩(wěn)越難發(fā)生。

圖2 不同局部粗糙壁面情況下幅值增長率隨雷諾數(shù)的變化Fig.2 Amplitudes'growth－rates with Reynolds numbers in boundary layers of the various localized roughness

壁面2、4、8誘導(dǎo)產(chǎn)生的擾動速度使原來二維基本流發(fā)生改變，即具有三維性（展向速度的存在），其速度剖面存有拐點，詳見圖3給出的壁面2、4、8中絕對值最大處的展向速度分布，其三維性強度遠(yuǎn)遠(yuǎn)大于其它壁面的強度。由于展向速度的存在，并很快形成復(fù)雜的流向渦、法向渦和展向渦，使得三維擾動波快速增長，產(chǎn)生橫流不穩(wěn)定，并加速轉(zhuǎn)捩，加劇流動的失穩(wěn)，有利于湍流的發(fā)生［9－10］。從圖1中也可知，壁面2、3、4、5、8的局部分布的區(qū)域（Re＝3500～3750），原因是該局部范圍內(nèi)基本速度場的三維性較強，且流動易于失穩(wěn)，在離局部位置較遠(yuǎn)處，三維特性逐漸減弱，流動的不穩(wěn)定性影響較弱。該結(jié)果與Croke等人的結(jié)果是一致的。

圖3 壁面2、4、8誘導(dǎo)出展向速度的分布曲線Fig.3 Spanwise velocity distribution in boundary layers induced by the 2、4、8wall

以三維擾動波的演化為例（a0＝0.006），利用線性理論分析Re＝4300，z＝0.91處的光滑壁面、壁面4的邊界層以及非線性作用后平均流得以修正后的中性曲線，分別如圖（4）a、b、c所示，由于壁面特性的不同使得中性曲線向低頻方向移動，而擾動后的平均流得以修正后進(jìn)一步使得中性曲線的下移，而且不穩(wěn)定區(qū)域也逐漸減小，這說明在基本流得以修正后，外界噪聲更容易促進(jìn)低頻不穩(wěn)定的形成，使得流動更易趨于失穩(wěn)，這與周恒、曹偉等人的研究成果一致［11］。

圖4 中性曲線Fig.4 Neutral curves

圖5中的曲線分別表示光滑壁面和壁面4中的流向擾動速度的空間演化（y＝0.1，z＝－1.83）。任意選取曲線4的流向擾動速度的傳播波形毗鄰正峰值位置x分別為19.4、36.6；與此對應(yīng)的曲線1流向擾動速度的傳播波形的毗鄰正峰值位置x分別為20.7、37.2，由于局部噴出和局部抽吸組合結(jié)構(gòu)誘導(dǎo)出三維擾動速度的影響，使上游三維擾動波的傳播產(chǎn)生了滯后?？紤]到擾動波的平均流向波數(shù)計算公式：＝2·π／（x2－x1），其數(shù)值結(jié)果分別為0.382和0.365，再定義三維擾動波與基本流動主流的夾角計算公式為：θ＝tan－1（β／），其數(shù)值結(jié)果分別為48.38°和49.67°。計算結(jié)果顯示，由于局部噴出和局部抽吸組合結(jié)構(gòu)誘導(dǎo)出三維擾動速度的影響，使三維擾動波的傳播方向發(fā)生了改變。而在邊界層的其他位置，擾動波的傳播方向也受到不同程度的變化，這可能是局部噴出和局部抽吸組合結(jié)構(gòu)誘導(dǎo)出三維擾動速度加速了邊界層流中流向渦、法向渦的形成及分布結(jié)構(gòu)等因素對三維擾動波傳播方向的影響相關(guān)，該結(jié)論與Jean－Francois Pinton等結(jié)果基本吻合［12］。但壁面6、7誘導(dǎo)產(chǎn)生的擾動流場是二維的，即三維擾動波在向下游演化時，傳播方向變化不明顯，這里不再一一表述。

圖5 流向擾動速度的空間演化Fig.5 Spatial evolution of stream wise disturbance velocity

3 結(jié) 論

本文采用直接數(shù)值模擬研究三維擾動波在光滑壁面和局部壁面特性改變的邊界層流中的非線性演化規(guī)律，分析了不同局部壁面特性改變對擾動波失穩(wěn)機制的影響。在局部壁面2、3、4、5、8特性改變的邊界層流中，三維擾動波幅值增長和幅值增長率均大于光滑壁面，這說明原二維基本流特性已改變，即具有三維性，且展向速度的存在改變原光滑壁面邊界層流中穩(wěn)定性的特征，使流動更容易失穩(wěn)，加速了湍流的形成和發(fā)展。另外，局部壁面6同樣也促使三維擾動波速度幅值增長，可是其增長率變化與光滑壁面相比較變化很小，這說明原二維基本流特性沒改變，基本保持光滑壁面邊界層流中穩(wěn)定性的特征。但局部壁面7對三維擾動波的增長起抑制作用，有易于流動得到控制和穩(wěn)定。數(shù)值結(jié)果顯示，局部壁面特性變化誘導(dǎo)產(chǎn)生的展向速度而形成橫流不穩(wěn)定特性是改變擾動波失穩(wěn)的另一重要因素。局部壁面特性的改變使三維擾動波向下游演化時的傳播方向的發(fā)生偏移現(xiàn)象?？傊?，局部壁面的類型、分布結(jié)構(gòu)（如壁面2、3、4、5、6、8）對擾動波增長都起到一定的激勵作用，激勵速率的快慢取決于局部壁面的類型、分布結(jié)構(gòu)及強度大小等因素；壁面7對擾動波的演化起抑制作用，有易于流動得到控制或?qū)α鲃悠鸬椒€(wěn)定作用的目的。

［1］SARIC W S，KOZLOV V V，LEVCHENKO V Y.Forced and.unforced subharmonic resonance in boundary layer transition［R］.AIAA paper 84－0007.

［2］SARIC W S，REED H L，WHITE E B.Stability and transition of three－dimensional boundary layers［J］.AnnualReviewofFluidMechanics，2003，（35）：413－440.

［3］HERBERT T.Parabolized stability equations［J］.AnnualReviewofFluidMechanics，1997，（29）：244－283.

［4］DIETMAR REMPFER.Low－dimensional modeling and numerical simulation of transition in simple shear flows［J］.AnnualReviewofFluidMechanics，2003，（35）：229－265.

［5］HUEBSCH W W，ROTHMAYER A P.Small－scale roughness effects on laminar separation［A］.Theoretical and Computational Fluid Dynamics［M］.Published online 17September 2003.

［6］BLOCH D R，MUELLER T J.Effects of distributed grit roughness on separation and transition on an airfoil at low Reynolds number［R］.AIAA paper 86－1788.

［7］FLORYAN J M.Stability of wall－bounded shear layers in the presence of simulated distributed surface roughness［J］.JournalofFluidMechanics，1997，335：29－55.

［8］LU Changgen，CAO Weidong，QIAN jianhua.A study on numerical method of Navier－Stokes equations and nonlinear－evolution of the coherent structure in a laminar boundary layer［J］.JournalofHydrodynamics，Ser.B，2006，18（3）：372－377.

［9］CHOUDLHAIR M.Roughness－induced generation of crossflow boundary vortices in three－dimensional boundary layers［J］.TheoreticalandComputational FluidDynamics，1994，6：1－30.

［10］COLISS S S，LELE S K.Receptivity to surface roughness near a swept leading edge［J］.JournalofFluid Mechanics，19998，380：285－315.

［11］CAO Wei，HUANG Zhang feng，ZHOU Heng.Study of the mechanism of breakdown in laminar－turbulent transition of a supersonic boundary layer on a flat plate［J］.AppliedMathematicsandMechanics，2006，27（4）：379－386.

［12］JEAN－FRANCOIS PINTON，GUILLAUME BRILLANT.Sound and vorticity：transmission and scatter［J］.TheoreticalandComputationalFluidDynamics，2005，18：413－433.