帶襟翼的機(jī)翼尾渦合并數(shù)值計(jì)算

王志博，孫 剛

(復(fù)旦大學(xué) 航空航天系，上海 200433)

帶襟翼的機(jī)翼尾渦合并數(shù)值計(jì)算

王志博，孫 剛

(復(fù)旦大學(xué) 航空航天系，上海 200433)

為獲得帶有襟翼的機(jī)翼尾渦的合并動(dòng)力學(xué)過程，在驗(yàn)證數(shù)值方法的基礎(chǔ)上，數(shù)值模擬帶有襟翼的機(jī)翼繞流尾流場，根據(jù)渦合并特征將合并過程劃分為誘導(dǎo)共轉(zhuǎn)階段、合并階段和軸對稱化階段.采用渦間距量綱—化尾流區(qū)域描述二渦誘導(dǎo)合并.變換弦長雷諾數(shù)、襟翼翼梢與機(jī)翼翼梢的間距、襟翼角度，改變襟翼翼梢渦與機(jī)翼翼梢渦的強(qiáng)度比，得到尾渦合并差異的特征.計(jì)算結(jié)果表明：隨著雷諾數(shù)的增加渦的強(qiáng)度增加，渦量的擴(kuò)散程度減低，渦合并過程被推遲，空間誘導(dǎo)運(yùn)動(dòng)過程得到延長，渦系空間誘導(dǎo)運(yùn)動(dòng)增強(qiáng)，渦合并的雷諾數(shù)效應(yīng)隨著雷諾數(shù)增加而減弱, 縮短渦間距，加速渦合并過程，但合并后的遠(yuǎn)場渦尺寸和形態(tài)沒有顯著改變; 襟翼渦隨著襟翼角的減小而減弱，強(qiáng)度逐漸減弱的襟翼渦逐步被翼梢渦拉伸卷吸，微弱的襟翼渦系在近場中的合并過程完全改變，翼梢渦的運(yùn)動(dòng)軌跡并未受到誘導(dǎo)運(yùn)動(dòng).翼梢渦合并的雷諾數(shù)效應(yīng)表現(xiàn)為誘導(dǎo)運(yùn)動(dòng)過程的增強(qiáng)，尾流中強(qiáng)度小的渦系起不到明顯的誘導(dǎo)作用.

襟翼渦; 翼梢渦; 雷諾數(shù); 渦合并; 環(huán)量

帶有翼面控制的航行體模型風(fēng)洞試驗(yàn)顯示，航行過程中機(jī)翼的翼梢渦和因襟翼打開導(dǎo)致翼面不連續(xù)而造成的渦系在近場的尾流中相互作用，這兩個(gè)渦系是尾流的主要成分.此外還包括從主機(jī)翼短艙結(jié)合部脫落的渦、尾翼翼梢渦、翼身組合體渦等各渦系，這些強(qiáng)度和尺寸各異的渦系在近場中最終合并為兩個(gè)對稱的集中渦.由于風(fēng)洞中均勻來流的整流作用，將這些低動(dòng)量高渦量的流體通過對流作用快速合并在一起.在近場尾流中，由于相互的誘導(dǎo)運(yùn)動(dòng)，脫體渦系發(fā)生了短波失穩(wěn)包括已知的橢圓失穩(wěn)、合作失穩(wěn)和其他形式的失穩(wěn)，渦在近場的互誘導(dǎo)作用使得渦核心區(qū)域渦量略有減低，渦發(fā)生了翻轉(zhuǎn)和擴(kuò)散，渦系逐步趨向于合并.

相互誘導(dǎo)形成的渦系運(yùn)動(dòng)對尾流中最終形成的渦的尺寸和強(qiáng)度以及持續(xù)時(shí)間等均起著決定性作用，研究這些渦系的空間演化，可揭示尾跡的演化規(guī)律.通過文獻(xiàn)歸納，按照渦雷諾數(shù)的大小將渦的合并可分為準(zhǔn)二維層流合并和三維短波橢圓合并兩種合并的形式.一般認(rèn)為200

針對高雷諾數(shù)情況下的飛行器三維空間繞流，為了抽象出合理的尾渦結(jié)構(gòu)的理論和實(shí)驗(yàn)?zāi)Ｐ?，F(xiàn)abre等[3]推導(dǎo)兩個(gè)平行的有渦核結(jié)構(gòu)的無限長渦絲的誘導(dǎo)失穩(wěn),并解釋了飛行器遠(yuǎn)場尾跡中經(jīng)常觀察的失穩(wěn)形式.Crouch[4]建立了二渦對的共轉(zhuǎn)失穩(wěn)模型，而Spalart[5]應(yīng)用對轉(zhuǎn)二渦對描述飛機(jī)尾流的優(yōu)化擾動(dòng)，用最優(yōu)擾動(dòng)分析來描述實(shí)際觀察到的飛機(jī)尾渦構(gòu)型是內(nèi)襟翼的脫落渦和平尾翼梢渦的失穩(wěn)可能性.Leweke等[6]研究了對轉(zhuǎn)渦對的三維短波失穩(wěn)，失穩(wěn)波長約為渦核的尺寸，通過細(xì)致的流場可視化和測速揭示了失穩(wěn)渦對的空間結(jié)構(gòu).文獻(xiàn)[6]發(fā)現(xiàn)了各渦的擾動(dòng)演化遵守一個(gè)運(yùn)動(dòng)匹配條件，每一個(gè)渦的耦合失穩(wěn)被稱為合作失穩(wěn).Crow[7]構(gòu)造了兩對對轉(zhuǎn)渦對空間演化，利用文獻(xiàn)[4]的對轉(zhuǎn)渦模型描述飛機(jī)遠(yuǎn)場尾渦的失穩(wěn).

近場生成的渦系不滿足文獻(xiàn)[7]關(guān)于柱渦間距離遠(yuǎn)大于渦的半徑的假設(shè)，也不滿足長波合作失穩(wěn)的基本假設(shè)，也不是初始小擾動(dòng)造成渦系失穩(wěn)，近場渦系存在于機(jī)翼脫體的自由剪切流中，連通集中渦的剪切層起到了加速渦系合并的作用.航行器的尾流中存在大量的共轉(zhuǎn)渦合并的現(xiàn)象，尾渦合并不僅是遠(yuǎn)場尾跡的形成的主要內(nèi)容，而且與近場渦致阻力緊密相關(guān)，航行器布局設(shè)計(jì)應(yīng)與渦致阻力之間通過渦合并動(dòng)力學(xué)過程建立橋梁；遠(yuǎn)場中存在的合作失穩(wěn)則對飛機(jī)的尾流尺度效應(yīng)產(chǎn)生影響，決定了機(jī)場的起降容積率.深入研究渦在近場的合并，將有助于揭示尾跡的形成機(jī)制和規(guī)模.Breitsamter[8]對飛機(jī)模型尾流風(fēng)洞試驗(yàn)進(jìn)行了系統(tǒng)總結(jié)，指出飛機(jī)機(jī)翼的翼面不連續(xù)處的脫落尾渦在近場中以共轉(zhuǎn)和對轉(zhuǎn)合并是尾流結(jié)構(gòu)形成的主要機(jī)制.Delbende等[9]擴(kuò)展到三維螺旋形渦的合并，指出螺旋形渦的自誘導(dǎo)作用對渦的合并起到了延遲作用.Laurent等[10]指出共轉(zhuǎn)合并后的渦系不穩(wěn)定性特征，Churchfield等[11]和Devenport等[12]對翼梢渦的合并進(jìn)行了系統(tǒng)實(shí)驗(yàn)和數(shù)值分析給出了翼尖渦初生的細(xì)節(jié).Bristol等[13]對平行渦系誘導(dǎo)合并過程進(jìn)行了模擬.

不同于經(jīng)典的二維情況下的共轉(zhuǎn)渦合并，渦合并的三維效應(yīng)顯示了獨(dú)特的三維渦拉伸，誘導(dǎo)共轉(zhuǎn)和剪切變形失穩(wěn)等獨(dú)有的物理現(xiàn)象.調(diào)整舵面布局，例如調(diào)整內(nèi)外襟翼的構(gòu)型和布局，合理配比設(shè)計(jì)翼面積分布與后掠角、扭轉(zhuǎn)角等構(gòu)型參數(shù)，制造合理的尾流渦系近場分布，延緩各類渦系共轉(zhuǎn)和對轉(zhuǎn)合并，對制造有益于航行性能的誘導(dǎo)下洗流場起到至關(guān)重要的作用，本文針對翼梢和襟翼脫落的渦系在尾流中的誘導(dǎo)共轉(zhuǎn)與合并運(yùn)動(dòng)特征，系統(tǒng)的描述尾流結(jié)構(gòu)演化.

1 數(shù)值模型與驗(yàn)證

1.1 計(jì)算設(shè)置

為研究繞流場中相鄰的兩個(gè)渦的相互誘導(dǎo)運(yùn)動(dòng)與合并，Bruin等[14-15]設(shè)計(jì)了如圖1所示的模型，并進(jìn)行了系統(tǒng)的風(fēng)洞試驗(yàn)，測試了機(jī)翼的翼梢脫體渦與襟翼翼梢發(fā)放的渦在尾流中的相互作用的尾流速度場.試驗(yàn)的弦長雷諾數(shù)為2.98×106，測量距機(jī)翼尾緣10倍的機(jī)翼展長至30倍的機(jī)翼展長范圍內(nèi)若干橫截面的速度分布.可供開展數(shù)值模擬的驗(yàn)證對比，試驗(yàn)?zāi)Ｐ偷脑敿?xì)尺寸見表1.

圖1 實(shí)驗(yàn)?zāi)Ｐ统叽?/p>

表1 SWIM模型參數(shù)

計(jì)算采用有限體積法求解雷諾平均方程結(jié)合SSTk-ω湍流模型與近壁面處理模型，近壁面區(qū)域網(wǎng)格加密滿足給定來流速度條件下y+≈1，速度與壓力的耦合采用PISO算法，采用二階迎風(fēng)對流格式離散動(dòng)量方程，求解格式細(xì)節(jié)與離散數(shù)值處理方法在此不做詳述.比照風(fēng)洞試驗(yàn)，選取來流速度為60 m/s的常溫空氣作為計(jì)算介質(zhì)，來流攻角為0°，遠(yuǎn)場邊界采用速度入口，對應(yīng)于風(fēng)洞試驗(yàn)中的來流湍流度為2.5%，遠(yuǎn)場壓力出口和對稱面邊界，壁面采用了無滑移邊界條件.為獲得更為詳細(xì)的尾流結(jié)構(gòu)，對尾流場的遠(yuǎn)場采用200倍的機(jī)翼展長作為擴(kuò)展計(jì)算域，并對尾流區(qū)域做了網(wǎng)格加密.

除了用于驗(yàn)證風(fēng)洞試驗(yàn)的算例之外，計(jì)算中省略了風(fēng)洞試驗(yàn)中存在的支柱和導(dǎo)流罩，這些支撐結(jié)構(gòu)對翼梢和襟翼后方的尾流流場造成的影響小，本文省略這些次要因素，網(wǎng)格無關(guān)性是通過渦核位置捕捉精度的比較來確定，最終選取網(wǎng)格為軸向×周向×展向 =450×100×50，網(wǎng)格劃分和邊界設(shè)置如圖2所示.

圖2網(wǎng)格劃分與邊界條件設(shè)置

1.2 數(shù)值計(jì)算結(jié)果驗(yàn)證

本研究選取機(jī)翼展向距離導(dǎo)流罩0.15 m處的截面壓力測試數(shù)據(jù)與數(shù)值計(jì)算進(jìn)行比較，壓力曲線分布比較如圖3所示.帶有襟翼的機(jī)翼增升效果顯著，機(jī)翼尾緣沒有顯著的壓力下降.

圖3 y/b1=0.5截面處壓力系數(shù)試驗(yàn)值與計(jì)算值比較

Fig.3 Pressure coefficients comparison from numerical simulation and experiment aty/b1=0.5

比較圖4、5可知，采用PIV試驗(yàn)獲得的尾渦位置與形態(tài)與數(shù)值計(jì)算結(jié)果基本一致，集中渦量的核心區(qū)位置略有差別，數(shù)值模擬襟翼渦與翼梢渦的形態(tài)和互誘導(dǎo)運(yùn)動(dòng)與實(shí)驗(yàn)基本一致.本文建立的數(shù)值方法能夠準(zhǔn)確模擬襟翼與翼梢形成的尾渦流動(dòng)結(jié)構(gòu)的演化.數(shù)值方法可獲得二渦在尾流中出現(xiàn)的合并現(xiàn)象，并描述近場渦系共轉(zhuǎn)渦對的合并造成渦的尺寸在尾流中迅速的擴(kuò)大等三維誘導(dǎo)運(yùn)動(dòng)現(xiàn)象.

圖4 風(fēng)洞試驗(yàn)測得的x/b1=1.25處尾渦系位置與數(shù)值結(jié)果比較

圖5 風(fēng)洞試驗(yàn)測得x/b1=5.0處尾渦系位置與數(shù)值結(jié)果比較

2 尾流場結(jié)構(gòu)特征

2.1翼梢渦與襟翼渦的誘導(dǎo)共轉(zhuǎn)階段x/Δ<20

在該階段中翼梢渦與襟翼渦各自具有獨(dú)立的渦系結(jié)構(gòu)，可利用經(jīng)典的蘭金渦理論模型描述渦核結(jié)構(gòu)，隨著渦脫體后，受橫流效應(yīng)和機(jī)翼脫體的自由剪切層影響，渦呈現(xiàn)軸對稱結(jié)構(gòu)，渦的尺寸增長但渦量衰減.隨后兩個(gè)柱狀渦受到互相的誘導(dǎo)而翻轉(zhuǎn)，渦核軌跡形成了螺旋線，柱狀集中渦由于互誘導(dǎo)而逐漸失去軸對稱的特征，兩個(gè)渦逐漸靠近.

2.2 共轉(zhuǎn)渦的合并20

由于二渦間存在機(jī)翼尾流的自由剪切層，該剪切層起到了連接二渦的作用，剪切層卷起后，二渦的渦核靠近，渦首先出現(xiàn)橢圓失穩(wěn)變形，渦相互靠近，直至不能區(qū)分兩個(gè)渦的邊界，無法辨識出渦核的各自位置.

2.3 結(jié)果渦的軸對稱化階段50

合并的渦在尾流中進(jìn)一步軸對稱化，在軸對稱化的過程中，結(jié)果渦的尺寸保持一定.渦的結(jié)構(gòu)逐步由橢圓形過渡到具有近似蘭金渦速度剖面的形狀.隨著尾渦在流場中的擴(kuò)散，尾渦的尺寸達(dá)到極限尺寸，渦量的峰值緩慢衰減下降.

可見上述合并的動(dòng)力學(xué)過程具有自相似的特征，即具有相似結(jié)構(gòu)的襟翼渦和翼梢渦合并后的結(jié)果渦也具有與子渦相似的結(jié)構(gòu)特征.尾流結(jié)構(gòu)劃分如圖6所示.

圖6 帶有襟翼的機(jī)翼尾渦流場劃分

在第1個(gè)階段由于二渦的獨(dú)立性強(qiáng)，渦間相互作用以互誘導(dǎo)的翻轉(zhuǎn)運(yùn)動(dòng)為主，那么渦系可利用下述渦辨識方法得到渦核心運(yùn)動(dòng)軌跡.對于柱狀渦而言，沿著下游發(fā)展的渦系軸向渦量在任意橫截面的分布為

對應(yīng)的環(huán)量為

單個(gè)渦的渦心位置為:

式中，積分上、下限以渦量的軸對稱分布區(qū)域而定.那么渦系的核心軌跡通過對截取有限個(gè)尾流中近場截面得到如圖7所示的互誘導(dǎo)翻轉(zhuǎn)運(yùn)動(dòng).

圖7 尾流場中翼梢渦核與襟翼渦核的翻轉(zhuǎn)運(yùn)動(dòng)軌跡

采用求解雷諾平均方程的方法求解流場.雖然實(shí)驗(yàn)中觀察到的渦在x/b1>10之后仍能存在高渦量的集中渦，但是由于過分?jǐn)?shù)值耗散，本文中給出的量綱—化的渦量耗散較試驗(yàn)顯著，圖8顯示了數(shù)值模擬給出的渦量在尾流中的耗散和擴(kuò)散.在合并之前的階段，設(shè)定襟翼渦渦強(qiáng)度Γ1和機(jī)翼的渦強(qiáng)度Γ2，根據(jù)數(shù)值計(jì)算得到渦在近場運(yùn)動(dòng)的軌跡，忽略黏性耗散作用，根據(jù)渦管強(qiáng)度保持定理，在近場認(rèn)為渦系的總的渦強(qiáng)度保持不變，那么有

圖8 襟翼渦與翼梢渦渦量衰減(U∞=60 m/s)

不考慮渦系的三維誘導(dǎo)作用造成的環(huán)量衰減，以平面渦系運(yùn)動(dòng)來近似模擬二渦的互誘導(dǎo)翻轉(zhuǎn)運(yùn)動(dòng).渦系中的第m個(gè)渦受其他渦誘導(dǎo)產(chǎn)生的復(fù)速度為

兩邊同乘以Γm，兩邊對m求和得到：

積分得到：

那么復(fù)數(shù)的實(shí)部和虛部分別為常數(shù)，可寫出點(diǎn)渦系的重心位置為:

實(shí)際在尾流中不同的位置，處于機(jī)翼脫體的剪切層中的集中渦量迅速的耗散.故每一個(gè)流場剖面內(nèi)的重心位置均不同，重心位置分布如圖9所示.重心位置的變化反映出二渦共轉(zhuǎn)誘導(dǎo)運(yùn)動(dòng)使渦處于不穩(wěn)定狀態(tài)，相互靠近逐漸誘導(dǎo)擴(kuò)散的特征.

3 渦間距對二渦合并影響

同一個(gè)雷諾數(shù)下，采用渦間距量綱—化尾流場截面位置，可完整的描述尾渦合并的過程，改變襟翼翼梢與機(jī)翼翼梢的間距Δ對二渦合并的影響，通過間距調(diào)整使二渦合并的過程得到改變，擴(kuò)展計(jì)算了襟翼翼梢與機(jī)翼翼梢不同間距情況下的二渦合并的情況， 圖10顯示了采用x/Δ作為量綱—化參數(shù)描述不同渦間距的二渦合并的過程，上述劃分適用于描述尾渦合并的3個(gè)階段.如圖10所示，縮短襟翼翼梢與機(jī)翼翼梢的間距，加速翼梢渦與襟翼渦的合并過程，使得合并渦的尺寸得到快速的增長并迅速達(dá)到極限尺寸.

圖9 襟翼渦與翼梢渦共轉(zhuǎn)誘導(dǎo)運(yùn)動(dòng)的重心軌跡

如圖11在x/Δ=60處的尾流場截面處，渦已經(jīng)處于合并后的軸對稱化階段，對于柱狀軸對稱渦而言，結(jié)果渦的尺寸并不隨著渦間距減小而顯著的增加，但渦量的峰值具有較大差異.表明量綱一的參數(shù)x/Δ無法描述軸對稱化后的渦量擴(kuò)散的過程.

圖10 不同翼梢間距的近場尾渦合并

圖11 x/Δ=60處合并后渦的渦量剖面

4 渦合并的雷諾數(shù)效應(yīng)

隨著來流速度的增加，尾流中集中渦的強(qiáng)度急劇的增加，合并過程必然隨著渦強(qiáng)度增加而產(chǎn)生差別.對于全尺度的實(shí)際流動(dòng)而言，渦的耗散過程是極為緩慢的，但是對于模型尺度的渦而言，渦的耗散相對較快.本文改變翼的展長對比發(fā)現(xiàn)渦的合并依賴于渦間距離和渦的初始尺寸與形態(tài).對于較低的雷諾數(shù)而言，集中渦總是尺寸快速增長，因此渦總會在近場首先合并.隨著雷諾數(shù)增加，渦的合并過程被推遲到更大的x/Δ值，因此近場的渦互誘導(dǎo)運(yùn)動(dòng)現(xiàn)象與低雷諾數(shù)略有不同，即使弦長雷諾數(shù)大于105后，尾流中渦的合并過程繼續(xù)隨著雷諾數(shù)的增長而被推遲.本文選取對應(yīng)于風(fēng)洞試驗(yàn)的標(biāo)準(zhǔn)算例對應(yīng)風(fēng)速的0.5倍和2.0倍進(jìn)行模擬，進(jìn)而分析渦系運(yùn)動(dòng)軌跡與合并雷諾數(shù)效應(yīng).圖12、13 顯示在同一個(gè)截面位置不同量級的雷諾數(shù)情況下兩個(gè)渦的合并狀態(tài).對比圖12、13可知隨著渦的強(qiáng)度增加渦量的擴(kuò)散程度減低，合并被推遲，獨(dú)立渦量的空間互誘導(dǎo)運(yùn)動(dòng)得到延長.

圖12 Δ=75 mm，x/Δ=20 處軸向渦量等值線云圖

圖13 Δ=25 mm，x/Δ=20 處軸向渦量等值線云圖

在與風(fēng)洞實(shí)驗(yàn)對應(yīng)的數(shù)值算例中發(fā)現(xiàn)，由于襟翼渦的渦量集中，初始強(qiáng)度約是翼梢渦的2倍，那么翼梢渦受到襟翼渦誘導(dǎo)而發(fā)生翻轉(zhuǎn)卷曲運(yùn)動(dòng)，如圖14所示隨著雷諾數(shù)的增加而逐漸增強(qiáng)，渦的卷曲效應(yīng)增加，由圖14可知翼梢渦翻轉(zhuǎn)軌跡曲率隨雷諾數(shù)而增大.在圖14中顯示二渦強(qiáng)度比在流向趨于接近的同時(shí)，二渦的誘導(dǎo)作用形成的軌跡趨于一致.

圖14 襟翼渦與翼梢渦近場軌跡(x/Δ<20)

5 襟翼渦強(qiáng)度減弱時(shí)的渦合并

通過改變襟翼角度，使得二渦的強(qiáng)度比在初始生成時(shí)減小，分析二渦強(qiáng)度變化對合并過程的影響.渦間誘導(dǎo)效應(yīng)決定于渦強(qiáng)度的大小和二者強(qiáng)度之比Γ1/Γ2.在上述驗(yàn)證模型算例的基礎(chǔ)上調(diào)整襟翼角度分別為20°、10°、5°，隨著襟翼角的增加，襟翼渦的強(qiáng)度不斷的增大，襟翼渦對翼梢渦的誘導(dǎo)作用變得越來越顯著，結(jié)果渦的尺寸有所增長.襟翼角為5°時(shí)，處于自由剪切尾流中的襟翼渦受到機(jī)翼的翼梢渦的拉伸而變形，形成快速的橢圓失穩(wěn)，逐漸拉伸成為剪切層被翼梢渦卷起，最終形成自由剪切層與機(jī)翼翼梢渦合并，翼梢渦并不發(fā)生翻轉(zhuǎn)運(yùn)動(dòng).隨著襟翼角度的增加，至10°時(shí)，襟翼渦的強(qiáng)度增大，襟翼渦與機(jī)翼翼梢渦的互誘導(dǎo)運(yùn)動(dòng)增強(qiáng)，翼梢渦出現(xiàn)了翻轉(zhuǎn)運(yùn)動(dòng)，襟翼渦與翼梢渦合并后的渦尺寸顯著的增大.在如圖15所示的尾流云圖可觀察到襟翼20°時(shí)二渦呈現(xiàn)了互誘導(dǎo)翻轉(zhuǎn)運(yùn)動(dòng).

圖15 尾流中渦量分布云圖對比(Δ=100 mm，截面位置x/Δ=2.5至15間隔2.5)

6 結(jié) 論

1)采用渦間距值進(jìn)行量綱—化，可劃分含襟翼渦的翼梢尾流結(jié)構(gòu)，根據(jù)運(yùn)動(dòng)特征可劃分為誘導(dǎo)共轉(zhuǎn)階段、共轉(zhuǎn)渦的合并、結(jié)果渦的軸對稱化等3個(gè)主要階段.但在對襟翼渦與翼梢渦的合并過程中未觀察到二渦互誘導(dǎo)翻轉(zhuǎn)運(yùn)動(dòng)軌跡超過180°的情況，說明襟翼渦系與翼梢渦系的三維空間拉伸失穩(wěn)存在極限.并且在同一個(gè)襟翼角和雷諾數(shù)下，雖然襟翼渦和翼梢渦的間距不同，但合并之后的結(jié)果渦具有相近的尺寸和形態(tài).

2)襟翼渦的強(qiáng)度遠(yuǎn)小于翼梢渦的強(qiáng)度時(shí)，在近場襟翼渦被拉伸變形并被翼梢渦卷起，而翼梢渦的運(yùn)動(dòng)軌跡并未受到襟翼渦的誘導(dǎo)運(yùn)動(dòng)，隨著襟翼角度的增加，襟翼渦強(qiáng)度增大，二渦形成了相互誘導(dǎo)運(yùn)動(dòng)與共轉(zhuǎn)合并.

3)當(dāng)Re<106時(shí)，隨著Re的增加，可觀察到尾流中一定程度的Re效應(yīng).隨著Re的增大，渦的強(qiáng)度增加，渦量的擴(kuò)散程度減低，合并過程被推遲，空間誘導(dǎo)運(yùn)動(dòng)得到延長，Re效應(yīng)減低.

[1] PAULO J, FERREIRA S, PEREIRA J. Reynolds number dependence of two-dimensional laminar co-rotating vortex merging[J]. Theoretical and Computational Fluid Dynamics, 2005, 19(1):65-75. DOI: 10.1007/s00162-004-0154-0.

[2] CERRETELLI C, WILLIAMSON C. The physical mechanism for vortex merging [J]. Journal of Fluid Mechanics, 2003, 475(1)：41-77. DOI: 10.1017/S0022112002002847.

[3] FABRE D, JACQUIN L, LOOF A. Optimal perturbations in a four-vortex aircraft wake in counter-rotating configuration [J]. Journal of Fluid Mechanics, 2002, 451(1):319-328.

[4] CROUCH J D. Instability and transient growth for two trailing-vortex pairs [J]. Journal of Fluid Mechanics, 1997, 350(1):311-330.

[5] SPALART R. Airplane trailing vortices[J]. Annual Review of Fluid Mechanics,1998,30:107-138.

[6] LEWEKE T, WILLIAMSON C H K. Cooperative elliptic instability of a vortex pair [J]. Journal of Fluid Mechanics, 1998, 360(1): 85-119.

[7] CROW S. Stability theory for a pair of trailing vortices[J]. AIAA Journal, 1970, 8(12):2172-2179. DOI: 10.2514/3.6083.

[8] BREITSAMTER C. Wake vortex characteristics of transport aircraft [J]. Progress in Aerospace Sciences, 2011, 47(2):89-134.DOI: 10.1016/j.paerosci.2010.09.002.

[9] DELBENDE I, PITON B, ROSSI M. Merging of two helical vortices [J]. European Journal of Mechanics-B/Fluids, 2015, 49(2):363-372.DOI: 10.1016/j.euromechflu.2014.04.005.

[10]LAURENT J, DAVID F, DENIS S, et al. Unsteadiness, instability and turbulence in trailing vortices[J]. Comptes Rendus Physique, 2005, 6(5):399-414. DOI:10.1016/j.crhy.2005.05.007.

[11]CHURCHFIELD J, BLAISDELL A. Numerical simulations of a wingtip vortex in the near field[J]. Journal of Aircraft, 2012, 46(1):230-243.DOI:10.2514/1.38086.

[12]DEVENPORT J, RIFE C, LIAPIS I. The structure and development of a wing-tip vortex[J]. Journal of Fluid Mechanics, 1996, 312(4):67-106. DOI: 10.1017/S0022112096001929.

[13]BRISTOL R L, ORTEGA J M, MARCUS P S. On cooperative instabilities of parallel vortex pairs[J]. Journal of Fluid Mechanics, 2004, 517(517):331-358. DOI:10.1017/S0022112004001016.

[14]BRUIN A. Test report for wake survey tests behind SWIM model geometry in DNWLST and DNW-LLF wind tunnels: NLR-TR-2001-183 [R]. Netherlands: DNW-LLF/ LST wind tunnels, 1999.

[15]BRUIN A,Ganzevles F. Data analysis of wake survey tests behind SWIM model geometry in DNW-LST and DNW-LLF wind tunnels: NLR-TR-2001-201 [R]. Netherlands: DNW-LLF/ LST wind tunnels, 2001.

(編輯 張 紅)

Numerical simulation of trailing vortex merge in a flapped wing wake

WANG Zhibo, SUN Gang

(Department of Aeronautics and Astronautics, Fudan University, Shanghai 200433, China)

A validated numerical simulation modeling of flapped wing wake is carried out by using RANS discretion. The flapped wing wake region is divided into three regions with significant different dynamic features. Three regions including co-rotation field, merging field and axial symmetry far field are induced. A non dimensional parameter of vortex interval is used to describe above regions. The Reynolds number, vortex interval, flap angle are used as control variables to change ratio of circulation of flap vortex and tip vortex. Reynolds effect is found as prolonged co-rotation field and enhanced rotation motion in wake. Decreased vortices interval leads an accelerated merging in near field but a limited dimension of resulted vortex. Small flap angle leads a weak flap vortex. The dominated tip vortex stretches and filament flap vortex in free share layer. The merging process is totally changed. The tip vortex trajectory does not expand. Reynolds number effect is the enhancement of induction in tip vortices merge. Induction of vortices with low intensity is not obvious in the wake.

flap vortex; wingtip vortex; Reynolds number; vortex merge; circulation

10.11918/j.issn.0367-6234.201510028

2015-10-12

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

王志博(1983—)，男，博士研究生; 孫 剛(1966—)，男，教授，博士生導(dǎo)師

孫 剛，Gang_sun@fudan.edu.cn

O355

A

0367-6234(2017)04-0073-07