黃章峰，逯學(xué)志，于高通

（1．天津大學(xué)力學(xué)系，天津 300072；2．空氣動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，四川 綿陽 621000）

機(jī)翼邊界層的橫流穩(wěn)定性分析和轉(zhuǎn)捩預(yù)測(cè)

黃章峰1，2，逯學(xué)志1，于高通1

（1．天津大學(xué)力學(xué)系，天津 300072；2．空氣動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，四川 綿陽 621000）

通過求解經(jīng)典O-S方程（LST）、擴(kuò)展的O-S方程（EOS）和線性拋物化穩(wěn)定性方程（LPSE），對(duì)展向無限長(zhǎng)、后掠角25°、迎角0°、來流Mach數(shù)0．8、單位Reynolds數(shù)6．79×106／m的機(jī)翼邊界層進(jìn)行了穩(wěn)定性分析，結(jié)合eN方法進(jìn)行了轉(zhuǎn)捩預(yù)測(cè)。研究發(fā)現(xiàn)無限長(zhǎng)后掠機(jī)翼在（x，ω）平面上的中性曲線沒有下支，在（x，β）平面上的中性曲線呈反拇指的形狀，橫流不穩(wěn)定性在機(jī)翼前緣占主導(dǎo)作用。當(dāng)外界擾動(dòng)進(jìn)入邊界層后，幅值將被直接放大，對(duì)于頻率相同的擾動(dòng)，首先是展向波數(shù)β大的增長(zhǎng)起來，演化到一定位置開始衰減，然后是展向波數(shù)β小的逐漸增長(zhǎng)起來，并且增長(zhǎng)的指數(shù)N逐步超過波數(shù)大的擾動(dòng)。轉(zhuǎn)捩在機(jī)翼前緣完成，引起轉(zhuǎn)捩的擾動(dòng)波的展向波長(zhǎng)約為2mm。

流動(dòng)穩(wěn)定性分析；轉(zhuǎn)捩預(yù)測(cè)；橫流不穩(wěn)定性；后掠機(jī)翼；拋物化穩(wěn)定性方程

0 引 言

機(jī)翼是飛機(jī)的重要部件之一，其最主要作用是產(chǎn)生升力，同時(shí)也產(chǎn)生一部分阻力。氣流在機(jī)翼表面的流動(dòng)狀態(tài)是層流還是湍流直接決定了摩擦阻力的大小。研究結(jié)果表明層流翼型的阻力比湍流翼型的阻力可以減少一半以上。若大型運(yùn)輸機(jī)翼能大部分保持層流，則可節(jié)省25%的燃料，具有很高的經(jīng)濟(jì)價(jià)值［1］。深入理解后掠機(jī)翼的轉(zhuǎn)捩機(jī)理、預(yù)測(cè)轉(zhuǎn)捩位置對(duì)氣動(dòng)力計(jì)算和翼型優(yōu)化設(shè)計(jì)是非常必要的，對(duì)機(jī)翼的氣動(dòng)減阻有著重要的經(jīng)濟(jì)價(jià)值。

后掠機(jī)翼的邊界層是典型的三維邊界層，其主要特點(diǎn)是在邊界層外緣勢(shì)流流速的橫向有壓力梯度。因此，邊界層內(nèi)的橫向速度不為零，稱為橫流。橫流的速度剖面存在拐點(diǎn)，因此流動(dòng)是不穩(wěn)定的，即所謂橫流不穩(wěn)定性。橫流不穩(wěn)定性的本質(zhì)和T-S波不同，T-S波的產(chǎn)生依靠粘性。除了不穩(wěn)定行進(jìn)波外，橫流不穩(wěn)定性還可能產(chǎn)生不穩(wěn)定的定常渦。通常在背景擾動(dòng)較小時(shí)，定常渦起主導(dǎo)作用，在背景擾動(dòng)較大時(shí)，行進(jìn)波在起主導(dǎo)作用。

早在十幾年前，機(jī)翼橫流不穩(wěn)定性和轉(zhuǎn)捩問題就是國(guó)外研究的熱點(diǎn)，如Malik［2-3］、Saric［4］、Haynes［5］等人的工作。近五年，除了橫流不穩(wěn)定性機(jī)理的研究外，添加粗糙度或吹吸的人工方法對(duì)橫流和轉(zhuǎn)捩位置進(jìn)行控制是熱點(diǎn)，如Mack［6］、Guha［7］等研究了橫流不穩(wěn)定性與壓縮性、Couette不穩(wěn)定性的關(guān)系，Chernoray［8］采用實(shí)驗(yàn)的方法、Li［9］采用拋物化穩(wěn)定性方程（PSE）、Nishino［10］和Rizzetta［11］采用直接數(shù)值模擬（DNS）方法研究了粗糙度對(duì)橫流不穩(wěn)定性的影響和轉(zhuǎn)捩位置的控制。

近五年國(guó)內(nèi)大部分的研究工作是橫流不穩(wěn)定性機(jī)理和轉(zhuǎn)捩預(yù)測(cè)，如楊永等人［12-13］在采用升華法測(cè)量轉(zhuǎn)捩位置的基礎(chǔ)上，結(jié)合線性穩(wěn)定性理論（LST）做了一系列的轉(zhuǎn)捩控制實(shí)驗(yàn)，王運(yùn)濤［14］、錢煒祺［15］、陳奕［16］采用不同的湍流模式研究了翼型轉(zhuǎn)捩的特點(diǎn)和轉(zhuǎn)捩位置的預(yù)測(cè)，張坤［17］、左歲寒［18］通過線性穩(wěn)定性分析（LST）、線性拋物化穩(wěn)定性方程（LPSE）研究了橫流不穩(wěn)定性，徐國(guó)亮［19］采用非線性拋物化穩(wěn)定性方程（NPSE）研究了曲率對(duì)機(jī)翼邊界層二次失穩(wěn)的影響。

轉(zhuǎn)捩預(yù)測(cè)eN方法的理論基礎(chǔ)是線性穩(wěn)定性理論。線性穩(wěn)定性分析方法可以分為兩類，一類是求解經(jīng)典的O-S方程，即線性穩(wěn)定性分析（LST），歸結(jié)為特征值問題，另一類是求解線性拋物化穩(wěn)定性方程（LPSE），歸結(jié)為初值問題。

借鑒LPSE方法，本文對(duì)LST進(jìn)行了改進(jìn)，提出用擴(kuò)展的O-S方程（EOS）來進(jìn)行穩(wěn)定性分析的方法。該方法仍然是求解特征值問題，部分考慮了邊界層內(nèi)非平行流等因素對(duì)擾動(dòng)幅值演化的影響。然后結(jié)合LST、EOS、LPSE和eN方法研究展向無限長(zhǎng)的后掠機(jī)翼穩(wěn)定性特性、擾動(dòng)演化規(guī)律及轉(zhuǎn)捩機(jī)理，對(duì)轉(zhuǎn)捩位置進(jìn)行預(yù)測(cè)。

1 數(shù)值方法

1.1 基本流場(chǎng)

采用CFD軟件計(jì)算后掠翼的三維流場(chǎng)，其中x方向?yàn)闄C(jī)翼弦長(zhǎng)方向，z方向?yàn)闄C(jī)翼展向方向，遠(yuǎn)離機(jī)身為正。收斂后，取展向均勻區(qū)域的流場(chǎng)截面，重新劃分并加密網(wǎng)格，其中流向方向沿機(jī)翼剖面的弧線方向，采用變間距網(wǎng)格，機(jī)翼兩頭密中間稀；法向方向的網(wǎng)格全部垂直于機(jī)翼壁面，y的大小定義為坐標(biāo)點(diǎn)到壁面的距離。采用直接數(shù)值模擬（DNS）的方法計(jì)算收斂后得到后掠翼的基本流場(chǎng)。

1.2 經(jīng)典的O-S方程（LST）

在正交曲線坐標(biāo)系下，考慮流向曲率drx和展向曲率drz，將N-S方程寫成擾動(dòng)形式，并進(jìn)行線化就得到了經(jīng)典的O-S方程。采用2點(diǎn)4階的Malik差分法，將O-S方程離散后得到特征關(guān)系式：

1.3 擴(kuò)展的O-S方程（EOS）

考慮，且基本流是x的函數(shù)將NS方程寫成擾動(dòng)形式，并進(jìn)行線化就得到了擴(kuò)展的O-S方程，離散后得到特征關(guān)系式：

1.4 線性拋物化穩(wěn)定性方程（LPSE）

無論是求解經(jīng)典的O-S方程還是擴(kuò)展的O-S方程，都是采用某個(gè)x位置的基本流求解當(dāng)?shù)氐奶卣髦郸?，β，ω和特征函?shù)φ，并對(duì)特征函數(shù)φ進(jìn)行歸一化，并沒有考慮擾動(dòng)幅值A(chǔ)本身的演化。將N-S方程拋物化和線化就得到線性拋物化穩(wěn)定性方程：

1.5 穩(wěn)定性分析及轉(zhuǎn)捩預(yù)測(cè)

在得到基本流的基礎(chǔ)上，通過求解經(jīng)典的O-S方程（LST）、擴(kuò)展的O-S方程（EOS）或線性拋物化方程（LPSE），就可以得到某個(gè)擾動(dòng)波的特征值α，β，ω和特征函數(shù)φ，進(jìn)行穩(wěn)定性分析，進(jìn)而進(jìn)行轉(zhuǎn)捩預(yù)測(cè)。轉(zhuǎn)捩預(yù)測(cè)方法常用的是eN方法，其基本概念是：若在邊界層內(nèi)存在各種頻率的小擾動(dòng)，這些擾動(dòng)向下游傳播過程中，當(dāng)其進(jìn)入各個(gè)擾動(dòng)相應(yīng)的不穩(wěn)定區(qū)域時(shí)，擾動(dòng)的幅值會(huì)被逐步放大，將各個(gè)頻率的擾動(dòng)從其幅值開始放大的位置起，沿傳播路徑累計(jì)其幅值的放大倍數(shù)en，取對(duì)數(shù)后為n值。在所有擾動(dòng)中，若有某一擾動(dòng)的幅值放大倍數(shù)的對(duì)數(shù)最先達(dá)到預(yù)設(shè)值N時(shí)，則可判定轉(zhuǎn)捩發(fā)生。該方法的理論基礎(chǔ)是線性穩(wěn)定性理論，而作為轉(zhuǎn)捩判定的域值N，則是由大量的實(shí)驗(yàn)和計(jì)算數(shù)據(jù)相結(jié)合給出。

在本問題中，擾動(dòng)演化是空間模式，因此ω為實(shí)數(shù)。展向?yàn)闊o限長(zhǎng)，流場(chǎng)在展向均勻，僅在流向變化，因此擾動(dòng)僅在流向有變化，故β為也實(shí)數(shù)，而α＝αr+iαi為復(fù)數(shù)。給定β和ω，通過LST、EOS或LPSE方法就可以得到不同流向位置下的αi，積分得到eN方法中的n值：

其中x0為擾動(dòng)開始增長(zhǎng)的位置或參考位置，積分到該擾動(dòng)開始衰減，即-αi＜0。改變?chǔ)潞挺?，求出n（β，ω，x）的包絡(luò)值，得到轉(zhuǎn)捩位置的N值：

1.6 物理模型

選用NACA0012翼型，弦長(zhǎng)1m，厚度0．12m，后掠角25°，迎角0°。取高空10000m處無量綱參數(shù)，溫度223．3K，密度0．4134kg／m3，粘性系數(shù)1．4579× 10-5Pa·s，來流Mach數(shù)0．8，單位Reynolds數(shù)6．79 ×106／m，無量綱長(zhǎng)度1mm，無量綱Reynolds數(shù)6790。下文除非特殊說明，均為無量綱量。

2 結(jié)果分析

2.1 網(wǎng)格及基本流場(chǎng)

用有限長(zhǎng)度機(jī)翼的流場(chǎng)代替無限長(zhǎng)機(jī)翼的流場(chǎng)。計(jì)算發(fā)現(xiàn)展向前緣的影響長(zhǎng)度大約是弦長(zhǎng)的2倍，后緣影響的長(zhǎng)度大約是弦長(zhǎng)的0．5倍，本文選取展向長(zhǎng)度為弦長(zhǎng)的4倍。圖1（a）給出了CFD整體網(wǎng)格示意圖，在32位系統(tǒng)中網(wǎng)格達(dá)到最大值70萬，在64位系統(tǒng)中網(wǎng)格為250萬。取出機(jī)翼展向均勻區(qū)域的流場(chǎng)截面，重新劃分并加密網(wǎng)格。圖1（b）給出了加密后局部網(wǎng)格示意圖。在流向，靠近機(jī)翼頭部的間距Δx＜0．1，在機(jī)翼中部的間距Δx＜1；在法向，1個(gè)邊界層內(nèi)有100個(gè)網(wǎng)格點(diǎn)，壁面處網(wǎng)格間距Δy＜0．008，確保網(wǎng)格滿足穩(wěn)定性分析的需要。

圖1 計(jì)算網(wǎng)格Fig.1 Computational grid

圖2（a-c）給出了機(jī)翼頭部不同x位置基本流的速度剖面u、v、w沿y方向的分布，其中u和w是平行于當(dāng)?shù)乇诿嫜亓飨蚝驼瓜虻乃俣确至?，v是垂直于當(dāng)?shù)乇诿嫜胤ㄏ虻乃俣确至俊？梢钥闯鲈趛＞1以上的速度u和w的剖面接近于常數(shù)。在機(jī)翼頭部出現(xiàn)了加速效應(yīng)，邊界層外的速度u隨著x的增大而增大，而速度w在邊界層外接近于常數(shù)，即來流速度在展向的分量。法向速度分量v要比流向速度u小2個(gè)量級(jí)，滿足局部平行流假設(shè)的要求。圖2（d）給出了邊界層名義、排移和動(dòng)量厚度沿流向的分布曲線。在x＝100處的邊界層名義厚度約為0．4mm，排移厚度約為0．2mm。機(jī)翼前部分的邊界層名義厚度是毫米量級(jí)。

圖2 基本流的速度剖面和邊界層厚度Fig.2 Pr of ile of mean flow and boundary layer thickness

2.2 特征值和特征函數(shù)

表1給出了采用LST、EOS、LPSE三種方法計(jì)算出來的不同流向位置的特征值?？梢钥闯鯡OS計(jì)算出的幅值增長(zhǎng)率-αi略大于LST的計(jì)算結(jié)果，而LPSE的結(jié)果比LST的結(jié)果要大30%左右，表明對(duì)無限長(zhǎng)后掠翼而言，除了擾動(dòng)波特征值-αi引起的幅值增長(zhǎng)外，還需要考慮幅值本身的增長(zhǎng)。

表1 擾動(dòng)波（β＝2.0，ω＝0.05）的增長(zhǎng)率Table 1 Growth rate of the wave（β＝2.0，ω＝0.05）

圖3（a）給出了擾動(dòng)波（β＝2．0，ω＝0．05）在x＝150．7處的特征函數(shù)沿法向的分布。速度｜u｜的最大值遠(yuǎn)離壁面，在邊界層排移厚度外，屬于第一模態(tài)的不穩(wěn)定波。通過LST和EOS計(jì)算得到的特征函數(shù)符合的很好。圖3（b）給出了該擾動(dòng)波的幅值增長(zhǎng)曲線，同樣可以看出LPSE的結(jié)果要比EOS和LST的結(jié)果大很多，LST和EOS的幅值最大增長(zhǎng)值n在5左右，而LPSE的結(jié)果達(dá)到了8左右。說明非平行性效應(yīng)在機(jī)翼頭部必須考慮。

圖3 擾動(dòng)波（β＝2.0，ω＝0.05）的特征函數(shù)和幅值演化曲線Fig.3 Eigen-function and growth curve of the wave（β＝2.0，ω＝0.05）

2.3 中性曲線和最大增長(zhǎng)率

圖4（a）給出了LST得到的不同β的（x，ω）平面上的中性曲線，即在該線上αi＝0。本文中取x0＝25作為穩(wěn)定性分析的開始位置。對(duì)于平板邊界層，中性曲線有上下兩支，并且存在一最小的位置xcr，在該位置上下兩支重合，形成一個(gè)類似于大拇指的外形，也稱為拇指曲線。在中性曲線內(nèi)，αi＜0，擾動(dòng)波是增長(zhǎng)的，在中性曲線外，αi＞0，擾動(dòng)波是衰減的。不同于平板邊界層，有后掠角的機(jī)翼邊界層不一定存在下支中性曲線。

從圖4（a）可以看出，當(dāng)β＝4．0時(shí)，中性曲線呈半開放式，即在中性曲線的左邊，αi＜0，擾動(dòng)波都是增長(zhǎng)的，而在中性曲線的右邊，αi＞0，擾動(dòng)是衰減的。隨著β的減小，中性曲線出現(xiàn)左半支，在左半支和右半支之間擾動(dòng)波是增長(zhǎng)的，之外是衰減的，但是仍然沒有下支。不同的β，在一定的x范圍內(nèi)，出現(xiàn)了ω＝0的擾動(dòng)波是增長(zhǎng)的，即出現(xiàn)不穩(wěn)定的駐波。從圖4（a）可以看出，當(dāng)β＝4，x＜150，而當(dāng)β＝2．5、2．0和1時(shí)，在x＜300的位置都出現(xiàn)了增長(zhǎng)的駐波，是典型的橫流不穩(wěn)定性問題。

圖4（b）給出了不同ω的（x，β）平面上的中性曲線。不同的ω，存在上、下兩支中性曲線，但是兩支曲線在右邊重合封閉，而向左邊是開口的，形成一個(gè)反拇指的形狀。而且在頻率ω從0到0．10的范圍內(nèi)，在x＜200位置的中性曲線非常接近，表明在這個(gè)范圍內(nèi)的擾動(dòng)都是增長(zhǎng)的。外界擾動(dòng)進(jìn)入邊界層后將直接被放大。

圖4 中性曲線Fig.4 Neutral curve

圖5（a）和（b）分別給出了不同β的最大增長(zhǎng)率-αi和相應(yīng)的ω。從圖5（a）可看出，當(dāng)β＝4時(shí)，越靠近前緣最大增長(zhǎng)率-αi越大，最大能到0．085，相應(yīng)的ω也越大，超過了0．20，但是增長(zhǎng)的范圍較小，在x＜150以內(nèi)。當(dāng)β＝2．5、2和1時(shí)，最大增長(zhǎng)率小于0．04，但是增長(zhǎng)的范圍很寬，從x＝50到x＝400，而且當(dāng)x＞150后最大增長(zhǎng)率曲線重合的很好，不同的展向波數(shù)β在該范圍內(nèi)的最大增長(zhǎng)率相同。從圖5（b）可以看出最不穩(wěn)定的擾動(dòng)波的頻率在ω＝0．05附近。

圖5 最大增長(zhǎng)率和相應(yīng)的頻率Fig.5 Maxinum growth rate and its corresponding frequence

2.4 擾動(dòng)波的幅值演化曲線

圖6（a）給出了LST計(jì)算得到的不同擾動(dòng)的幅值演化曲線，可以看出，當(dāng)擾動(dòng)頻率相同時(shí)，首先是展向波數(shù)β大的擾動(dòng)波先增長(zhǎng)起來，隨著往下游的演化，到一定位置又開始衰減，展向波數(shù)β小的波也逐漸增長(zhǎng)起來，并且增長(zhǎng)的指數(shù)n逐步超過波數(shù)大的擾動(dòng)。由于中性曲線向左是開口的，這些擾動(dòng)在開始積分的起點(diǎn)位置x0＝25處都是不穩(wěn)定的，因此所有擾動(dòng)從積分起點(diǎn)就開始增長(zhǎng)。EOS和LPSE的計(jì)算結(jié)果與LST的一致，僅僅是擾動(dòng)曲線幅值n的大小不同。圖6（b）給出了LPSE的計(jì)算結(jié)果，趨勢(shì)與LST相同，不同的是增長(zhǎng)幅度n值要大。

2.5 轉(zhuǎn)捩預(yù)測(cè)

給定擾動(dòng)的展向波數(shù)β和頻率ω就能得到一條幅值演化曲線，這些幅值演化曲線的包絡(luò)線就是轉(zhuǎn)捩預(yù)測(cè)eN方法中N值曲線。從圖4可以看出幅值演化曲線能達(dá)到最大值的擾動(dòng)波的展向波數(shù)范圍為0．5到8，頻率范圍為0到0．5。本文中取波數(shù)間隔為0．1，頻率間隔為0．05，計(jì)算了全部擾動(dòng)波的幅值演化曲線，在每個(gè)位置上取全部擾動(dòng)波中最大的n值，作為該處的包絡(luò)值N。

圖6 幅值增長(zhǎng)曲線Fig.6 Growth curve

圖7給出了基于LST、EOS、LPSE三種方法得到了包絡(luò)曲線，其中線條表示所有擾動(dòng)波的包絡(luò)線，而符號(hào)表示駐波的包絡(luò)線?？梢钥闯霾捎肔ST方法預(yù)測(cè)的N值相對(duì)EOS和LPSE的值偏小，最大值小于5，EOS的結(jié)果略大于LST的結(jié)果，最大值接近5，而LPSE的結(jié)果比LST和EOS均要大，最大值接近9。

圖7 轉(zhuǎn)捩預(yù)測(cè)方法eN中的N值Fig.7 N-factor in eNmethod

在后掠機(jī)翼中存在橫流不穩(wěn)定性，其擾動(dòng)波的波矢方向跟主速度矢量方向垂直，而且擾動(dòng)頻率很小，接近于0，常以駐波的形式出現(xiàn)。對(duì)比所有擾動(dòng)波的包絡(luò)線和駐波的包絡(luò)線，可以看出，橫流不穩(wěn)定波的預(yù)測(cè)N值略小，表明橫流不穩(wěn)定性占主導(dǎo)作用。

采用LPSE的方法預(yù)測(cè)的轉(zhuǎn)捩N值達(dá)到了9，表明轉(zhuǎn)捩位置在機(jī)翼的前緣。實(shí)際的轉(zhuǎn)捩位置跟感受性問題有關(guān)，取決于外界擾動(dòng)進(jìn)行邊界層后的擾動(dòng)展向波數(shù)、頻率和幅值大小。不同的x位置，對(duì)應(yīng)的N值不同，相應(yīng)的最不穩(wěn)定波的展向波數(shù)和頻率也不同。研究結(jié)果表明，如果擾動(dòng)波的無量綱的幅值達(dá)到20%就發(fā)生轉(zhuǎn)捩了［20］。表2給出了x＝50、100和200位置處發(fā)生轉(zhuǎn)捩時(shí)的N值及其可能引發(fā)轉(zhuǎn)捩的擾動(dòng)波的展向波數(shù)β、頻率ω和在x0＝25處的擾動(dòng)幅值u0。

表2 在x處發(fā)生轉(zhuǎn)捩時(shí)（u＝0.2U∞）的N值，以及引起轉(zhuǎn)捩發(fā)生的擾動(dòng)波的波數(shù)β、頻率ω和在x0＝25處的幅值u0（m／s）Table 2 N-factor at the transition location x where u＝0.2U∞，and the wave numberβ，frequencyωand the amplitude u0（m／s）at x0＝25 of the inducing wave

以LST的計(jì)算結(jié)果為例，如果展向波數(shù)為4，頻率為0．15的擾動(dòng)波進(jìn)入了x0＝25處的幅值為7．93m／s，那么該擾動(dòng)傳播到達(dá)x＝50位置時(shí)的幅值將被放大e1．8＝6倍，即48．8m／s，達(dá)到了來流速度U∞＝240m／s的20%，可以認(rèn)為轉(zhuǎn)捩發(fā)生了。而EOS和LPSE的計(jì)算結(jié)果表明在x0＝25處的幅值分別為5．88m／s和2．16m／s就能在x＝50處引發(fā)轉(zhuǎn)捩。而若轉(zhuǎn)捩在x＝200處發(fā)生，LST、EOS、LPSE三種方法計(jì)算出擾動(dòng)在x0＝25的幅值僅為0．65m／s、0．36m／s和0．04m／s。

從表2還可以看出，雖然三種穩(wěn)定性分析方法獲得的N值不完全相同，但預(yù)測(cè)的最可能引發(fā)轉(zhuǎn)捩的擾動(dòng)波的展向波數(shù)在3．0到3．5之間，頻率在0．035到0．060之間，相應(yīng)的展向波長(zhǎng)為約為2mm。

3 結(jié) 論

本文在采用CFD和DNS相結(jié)合的方法計(jì)算得到無限長(zhǎng)后掠翼基本流的基礎(chǔ)上，通過求解經(jīng)典的O-S方程（LST）、擴(kuò)展的O-S方程（EOS）和線性拋物化穩(wěn)定性方程（LPSE）進(jìn)行穩(wěn)定分析和轉(zhuǎn)捩預(yù)測(cè)，得到以下結(jié)論：

（1）后掠翼在（x，ω）平面上的中性曲線沒有下支，在（x，β）平面上的中性曲線呈反拇指的形狀，橫流不穩(wěn)定性在機(jī)翼前緣占主導(dǎo)作用。

（2）當(dāng)外界擾動(dòng)進(jìn)入邊界層后，幅值將被直接放大，首先是展向波數(shù)β大的擾動(dòng)波增長(zhǎng)起來，然后是波數(shù)β小的波逐漸增長(zhǎng)起來，并且增長(zhǎng)的指數(shù)N逐步超過波數(shù)大的擾動(dòng)。

（3）轉(zhuǎn)捩在機(jī)翼前緣完成，引起轉(zhuǎn)捩的擾動(dòng)波的展向波長(zhǎng)約為2mm。

［1］周恒，趙耕夫．流動(dòng)穩(wěn)定性［M］．北京：國(guó)防工業(yè)出版社，2004．

［2］MALIK M R，LI F，CHANG C L．Crossflow disturbances in three-dimensional boundary layers：nonlinear development，wave interaction and secondary instability［J］．J．Flutd Mech．，1994，268：1-36．

［3］MALIK M R，LI F，CHOUDHARI M M，et al．Secondary instability of cross-flow vortices and swept-wing boundary-layer transition［J］．J．Flutd Mech．，1999，399：85-115．

［4］SARIC W S，CARRILLO J R，REIBERT M S．Nonlinear stability and transition in 3-D boundary layers［J］．Meccantca，1998，33：469-487．

［5］HAYNES T S，REED H L．Simulation of swept-wing vortices using nonlinear parabolized stability equations［J］．J．Flutd Mech．，2000，405：325-349．

［6］MACK C J，SCHMID P J，SESTERHENN J L．Global stability of swept flow around a parabolic body：connecting attachment-line and crossflow modes［J］．J．Flutd Mech．，2008，611：205-214．

［7］GUHA A，F(xiàn)RIGAARD I A．On the stability of plane Couette-Poiseuille flow with uniform crossflow［J］．J．Flutd Mech．，2010，656：417-447．

［8］CHERNORAY V G，DOVGAL A V，KOZLOV V V，et al．Secondary instability of a swept-wing boundary layer disturbed by controlled roughness elements［J］．Journal of Vtsualtzatton，2010，13（3）：251-256．

［9］LI F，CHOUDHARI M M，CHANG C L，et al．Computational modeling of roughness-based laminar flow control on a subsonic swept wing［J］．AIAA Journal，2011，49（3）：520-529．

［10］NISHINO T，SHARIFF K．Direct numerical simulation ofa swept-wing boundary layer with an array of discrete roughness elements［A］．Seventh Iutam Symposium on Laminar-Turbulent Transition［C］，2010，18：289-294．

［11］RIZZETTA D P，VISBAL M R，REED H L，et al．Direct numerical simulation of discrete roughness on a sweptwing leading edge［J］．AIAA Journal，2011，48（11）：2660-2673．

［12］楊永，左歲寒，李喜樂，等．基于升華法實(shí)驗(yàn)研究后掠翼三維邊界層的轉(zhuǎn)捩［J］．實(shí)驗(yàn)流體力學(xué)，2009，23（3）：40-43．

［13］左歲寒，楊永，李棟，等．后掠翼邊界層橫流不穩(wěn)定性的實(shí)驗(yàn)研究（英文）［J］．空氣動(dòng)力學(xué)學(xué)報(bào)，2010，28（5）：495-502．

［14］王運(yùn)濤，王光學(xué)，張玉倫．30P-30N多段翼型復(fù)雜流場(chǎng)數(shù)值模擬技術(shù)研究［J］．空氣動(dòng)力學(xué)學(xué)報(bào)，2010，28（1）：99-103．

［15］錢煒祺，RANDOLPH C，LEUNG K．考慮轉(zhuǎn)捩影響的翼型動(dòng)態(tài)失速數(shù)值模擬［J］．空氣動(dòng)力學(xué)學(xué)報(bào)，2008，26（1）：50-55

［16］陳奕，高正紅．Gamma-Theta轉(zhuǎn)捩模型在繞翼型流動(dòng)問題中的應(yīng)用［J］．空氣動(dòng)力學(xué)學(xué)報(bào)，2009，27（4）：411-418．

［17］張坤，宋文萍．基于線性穩(wěn)定性分析的e～N方法在準(zhǔn)確預(yù)測(cè)翼型氣動(dòng)特性中的應(yīng)用［J］．西北工業(yè)大學(xué)學(xué)報(bào)，2009，27（3）：294-299．

［18］左歲寒，楊永，李棟．基于線性拋物化穩(wěn)定性方程的后掠翼邊界層內(nèi)橫流穩(wěn)定性研究［J］．計(jì)算物理，2010，27（5）：665-670．

［19］徐國(guó)亮，符松．曲率對(duì)機(jī)翼邊界層二次失穩(wěn)影響［J］．力學(xué)學(xué)報(bào)，2010，42（6）：995-1005．

［20］黃章峰，曹偉，周恒．超音速平板邊界層轉(zhuǎn)捩中層流突變?yōu)橥牧鞯臋C(jī)理-時(shí)間模式［J］．中國(guó)科學(xué)，G輯，2005，35（5）：537-547．

Cross-flow instability analysis and transition prediction of airfoil boundary layer

HUANG Zhangfeng1，2，LU Xuezhi1，YU Gaotong1
（1．Department of Mechanics，Tianjin University，Tianjin 300072，China；2．State Key Laboratory of Aerodynamics，Mianyang Sichuan 621000，China）

Cross-flow instability analysis and transition prediction have been done on a 25°swept，0°attacked airfoil with Mach number 0．8，unit Reynolds number 6．79×106／m by solving the classic O-S equation（LST），extended O-S equation（EOS）and linear parabolic stable equation（LPSE）together with theeNmethod．Results show that the lower branch of the neutral curve in（x，ω）plane does not exist and the shape of the neutral curve in（x，β）plane appears as an anti-thumb．The cross-flow instability plays a key role in the leading edge．As disturbances enter the boundary layer，their amplitude will be enlarged directly．For the same frequence，the disturbances with big wave numberβwill grow firstly and propogate to a certain place where they begin to decay，then the disturbances with small wave numberβincreased gradually and their increasing factorNwill exceed that with big wave number finally．Transition occurs in the leading edge of the airfoil and the wave lengh of the inducing wave is about 2mm．

instability analysis；transition prediction；cross-flow instability；swept airfoil；parabolic stable equation

V211.1+9；O357.4+1

A doi:10.7638／kqdlxxb-2012.0062

0258-1825（2014）01-0014-07

2012-08-20；

2012-11-15

空氣動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室開放課題資助（SKLA201101，SKLA201201）

黃章峰（1977-），男，湖南宜章人，副教授，研究方向：流動(dòng)穩(wěn)定性、轉(zhuǎn)捩及湍流、直接數(shù)值模擬．E-mail：hzf＠tju．edu．cn

黃章峰，逯學(xué)志，于高通．機(jī)翼邊界層的橫流穩(wěn)定性分析和轉(zhuǎn)捩預(yù)測(cè)［J］．空氣動(dòng)力學(xué)學(xué)報(bào)，2014，32（1）：14-20．

10．7638／kqdlxxb-2012．0062．HUANG Z F，LU X Z，YU G T．Cross-flow instability analysis and transition prediction of airfoil boundary layer［J］．ACTA Aerodynamica Sinica，2014，32（1）：14-20．