張冬, 胡孟權(quán), 王旭, 吳章沅

(空軍工程大學(xué) 航空航天工程學(xué)院, 陜西 西安 710038)

?

前后掠翼鴨式布局中鴨翼渦的流動(dòng)機(jī)理

張冬, 胡孟權(quán), 王旭, 吳章沅

(空軍工程大學(xué) 航空航天工程學(xué)院, 陜西 西安 710038)

摘要:基于前后掠鴨式布局的簡化模型,通過求解雷諾平均N-S方程,模擬了前后掠鴨式布局的繞流結(jié)構(gòu),得到了不同布局下鴨翼的升力系數(shù)曲線。通過空間流線圖,分析了單獨(dú)鴨翼漩渦的發(fā)展特點(diǎn),以及不同布局中鴨翼渦與機(jī)翼前緣渦的干擾機(jī)理。結(jié)果表明:在后掠翼鴨式布局中,鴨翼渦在大迎角時(shí)受到機(jī)翼前緣渦的有利干擾,增大了鴨翼的升力系數(shù),提高了失速迎角;在前掠翼鴨式布局中,鴨翼的最大升力系數(shù)有所提高,失速迎角基本保持不變。

關(guān)鍵詞:鴨式布局; 鴨翼渦; 流動(dòng)機(jī)理

0引言

近距耦合鴨式布局渦系之間可以產(chǎn)生有利干擾,提高布局的升力系數(shù)和失速迎角[1-4]?，F(xiàn)代戰(zhàn)斗機(jī)追求機(jī)動(dòng)性和敏捷性,尤其是大迎角和過失速機(jī)動(dòng)能力,近距耦合鴨式布局是現(xiàn)代戰(zhàn)斗機(jī)經(jīng)常采用的先進(jìn)布局之一。

瑞典的Behrbohm[5]在20世紀(jì)60年代中期研究發(fā)現(xiàn), 近距耦合鴨式布局的鴨翼渦可對(duì)機(jī)翼前緣渦產(chǎn)生有利干擾, 從而增加布局的升力系數(shù), 并提高布局的失速迎角,SAAB-37飛機(jī)的設(shè)計(jì)就成功運(yùn)用了該研究成果。此后,許多戰(zhàn)斗機(jī)采用近距耦合鴨式布局,如瑞典的JAS-39、法國的幻影和陣風(fēng)等,都屬于后掠翼鴨式布局;美國的X-29和俄羅斯的Su-47則是前掠翼鴨式布局中的典型代表。

在近距耦合鴨式布局中,鴨翼兼有氣動(dòng)增升部件及操縱面雙重功能[6],鴨翼的氣動(dòng)特性對(duì)主機(jī)翼以及整機(jī)的性能有著十分重要的影響。對(duì)氣動(dòng)布局中的鴨翼氣動(dòng)特性進(jìn)行深入分析,有助于充分利用鴨翼帶來的氣動(dòng)收益。在前后掠翼鴨式布局中,由于機(jī)翼的流動(dòng)存在不同的特點(diǎn),鴨翼渦受到機(jī)翼前緣渦的干擾也有所不同,對(duì)這兩種布局下鴨翼渦流動(dòng)機(jī)理的研究是十分必要的。對(duì)于復(fù)雜渦系的干擾分析,實(shí)驗(yàn)研究有很大的難度,而數(shù)值仿真及其后處理則有著獨(dú)特的優(yōu)勢,本文通過數(shù)值計(jì)算,研究了前后掠鴨式布局中鴨翼渦的發(fā)展過程,分析了機(jī)翼前緣渦對(duì)鴨翼渦的干擾作用。

1計(jì)算模型和網(wǎng)格

由于本文主要研究鴨翼渦的流動(dòng)機(jī)理以及機(jī)翼對(duì)鴨翼流動(dòng)的干擾,為方便分析,計(jì)算采用平板模型(C代表單獨(dú)鴨翼,BC代表后掠翼鴨式布局, BSC代表帶邊條的后掠翼鴨式布局, FC代表前掠翼鴨式布局,FSC代表帶邊條的前掠翼鴨式布局),模型厚度4 mm,迎風(fēng)面邊緣皆倒角45°,從而在小迎角下就能產(chǎn)生穩(wěn)定的漩渦。鴨翼與機(jī)翼共面,鴨翼前緣后掠角56°,后緣后掠角25°,鴨翼翼根弦長4 mm,翼尖弦長24 mm,展長100 mm。保持前后掠機(jī)翼的翼面積、展弦比、根梢比和1/4弦線后掠角相同,邊條后掠70°,圖1為不同鴨翼布局的平面示意圖。

圖1　不同鴨翼布局平面示意圖Fig.1　Diagram of different canard configurations

由于模型的幾何對(duì)稱性,本文在計(jì)算時(shí)采用半模,并采用結(jié)構(gòu)網(wǎng)格,圖2為BSC和FSC布局的網(wǎng)格分布。網(wǎng)格劃分時(shí),邊界層內(nèi)第1層網(wǎng)格高度控制為機(jī)翼平均氣動(dòng)弦長的10-6,y+值在0～1之間,以滿足飛機(jī)表面粘性邊界層的計(jì)算要求[7],從而保證能夠模擬邊界層內(nèi)的流動(dòng)特征。各模型網(wǎng)格數(shù)量控制在1.5×106～3.0×106。

圖2　BSC和FSC布局網(wǎng)格分布Fig.2　Surface mesh of BSC and FSC configuration

2計(jì)算方法

數(shù)值模擬的控制方程采用三維Navier-Stokes方程。在直角坐標(biāo)系中,無熱源的三維Navier-Stokes方程守恒形式為[8]:

式中:w為狀態(tài)矢量;f為無粘(對(duì)流)通矢量項(xiàng);fv為粘性(耗散)通矢量項(xiàng)。

湍流模型選用SST (Shear-Stress Transport)k-ω模型,考慮了低雷諾數(shù)和剪切流,加入了橫向耗散導(dǎo)數(shù)項(xiàng),并在定義湍流粘度時(shí)考慮了湍流剪切應(yīng)力的輸運(yùn)過程,適合有逆壓梯度的流動(dòng)計(jì)算[9]。

采用有限體積法對(duì)控制方程進(jìn)行離散,對(duì)流項(xiàng)選用二階迎風(fēng)差分格式[10]。遠(yuǎn)場條件為壓力遠(yuǎn)場,對(duì)稱面為對(duì)稱邊界條件,物面為無滑移壁面條件,計(jì)算殘差收斂精度為10-5。

3計(jì)算結(jié)果及分析

3.1不同鴨翼的升力系數(shù)

圖3給出了不同布局下鴨翼氣動(dòng)升力系數(shù)的計(jì)算結(jié)果?？梢钥闯?后掠翼鴨式布局能大幅提高鴨翼最大升力系數(shù),并延遲失速,尤其是帶邊條后掠翼鴨式布局,不僅能大幅提高最大升力系數(shù),還能使失速后鴨翼的升力系數(shù)下降變得平緩,極大地改善了鴨翼的失速特性,這是因?yàn)榇笥窍馒喴頊u與機(jī)翼前緣渦之間存在著有利干擾;對(duì)于前掠翼鴨式布局,鴨翼的最大升力系數(shù)有少量提高,失速迎角基本不變,失速后升力系數(shù)下降緩慢;在帶邊條的后掠翼鴨式布局中,鴨翼的最大升力系數(shù)和失速迎角都有明顯提高,這說明在大迎角下,鴨翼渦、邊條渦和機(jī)翼前緣渦之間的有利干擾更加明顯。

圖3　升力系數(shù)曲線Fig.3　Lift coefficient curves

3.2單獨(dú)鴨翼布局

對(duì)于在中小迎角下的單獨(dú)鴨翼,鴨翼渦隨著迎角的增加而增強(qiáng),對(duì)翼面的控制范圍也逐漸擴(kuò)大,升力系數(shù)隨迎角增加呈線性增長;當(dāng)迎角繼續(xù)增加,鴨翼渦渦核開始從后向前破裂,渦核破裂后的漩渦能量下降,但漩渦在不斷擴(kuò)大,對(duì)翼面的控制區(qū)域也相應(yīng)擴(kuò)大,所以升力系數(shù)隨迎角增加有緩慢的增長;當(dāng)迎角進(jìn)一步增加,鴨翼渦渦核破裂至前緣,鴨翼渦對(duì)翼面的控制作用達(dá)到最大,升力系數(shù)達(dá)到最高點(diǎn);當(dāng)迎角繼續(xù)增加,鴨翼渦外側(cè)氣流能量下降,對(duì)翼面的控制能力減弱,升力系數(shù)下降,鴨翼進(jìn)入失速狀態(tài);進(jìn)一步增加迎角情況下,鴨翼渦外側(cè)氣流破散,氣流出現(xiàn)大規(guī)模倒流,升力系數(shù)明顯下降,失速進(jìn)一步惡化。

圖4給出了不同迎角下鴨翼渦的發(fā)展變化。為方便對(duì)比,本文所有布局都采用相同的壓力云圖,圖例分布如圖4左側(cè)所示,后面的壓力云圖都采用相同的圖例分布。

圖4　鴨翼空間流線圖Fig.4　Streamlines of canard

可以看出:迎角10°時(shí),鴨翼渦渦核開始從后緣破裂;當(dāng)迎角15°時(shí),渦核破裂至鴨翼中部;當(dāng)迎角達(dá)到21°時(shí),渦核破裂至前緣,升力系數(shù)達(dá)到最大;繼續(xù)增大迎角至23°,鴨翼渦破裂,鴨翼進(jìn)入失速狀態(tài)。

3.3后掠翼鴨式布局

圖5　BC布局的空間流線圖Fig.5　Streamlines of BC configuration

圖5為BC布局的空間流線圖。圖5(a)給出了迎角23°時(shí)的情況,可以看到BC布局中的鴨翼渦和機(jī)翼前緣渦都保持穩(wěn)定狀態(tài),而此時(shí)單獨(dú)鴨翼和單獨(dú)后掠機(jī)翼的漩渦早已破裂,如圖5(b)所示。當(dāng)迎角繼續(xù)增加,兩渦能量減小,漩渦破裂,鴨翼出現(xiàn)明顯失速,如圖5(c)所示,在25°迎角時(shí),鴨翼渦和機(jī)翼前緣渦都已破裂。

圖6為機(jī)翼70%翼根弦長位置上的截面流線圖?？梢钥闯觯壶喴頊u與機(jī)翼前緣渦相互切洗,且有卷繞的趨勢,所以兩渦能量得到加強(qiáng),漩渦變得穩(wěn)定。

圖6　截面流線圖Fig.6　Streamlines on cross section

3.4帶邊條后掠翼鴨式布局

圖7為不同迎角下BSC布局的空間流線圖?？梢钥闯觯涸?5°迎角時(shí),邊條渦和機(jī)翼前緣渦相互卷并,卷并后的漩渦能量得到增強(qiáng),變得更加穩(wěn)定,與鴨翼渦的有利誘導(dǎo)作用更加明顯。一方面,鴨翼渦控制著卷并后的漩渦向機(jī)翼前緣破裂的速度減慢,使之變得更加穩(wěn)定;另一方面,卷并后的漩渦對(duì)鴨翼渦形成有利干擾,使鴨翼渦變得穩(wěn)定不易破裂,并在很大迎角范圍仍可保持較大的升力系數(shù),且無明顯失速。在28°迎角時(shí),鴨翼渦依然比較穩(wěn)定,還未破裂,只是能量有所減小,升力系數(shù)稍有下降,但沒有出現(xiàn)明顯失速,直到30°迎角時(shí),鴨翼渦才開始破裂。

圖7　BSC布局的空間流線圖Fig.7　Streamlines of BSC configuration

3.5前掠翼鴨式布局

對(duì)于FC布局,機(jī)翼前緣渦對(duì)鴨翼渦存在兩方面的作用,圖8給出了15°迎角時(shí)FC布局和單獨(dú)鴨翼及機(jī)翼的空間流線圖。可以看出:FC布局中的機(jī)翼前緣渦和鴨翼渦呈一定角度發(fā)展,鴨翼渦對(duì)機(jī)翼前緣渦產(chǎn)生外推和擠壓作用,機(jī)翼前緣渦不能正常發(fā)展,對(duì)翼根處的控制作用降低;反過來,機(jī)翼前緣渦的上洗作用也阻礙著鴨翼渦的向后流動(dòng),形成兩渦之間的不利干擾。

圖8　FC布局和單獨(dú)部件的空間流線圖Fig.8　Streamlines of FC and single canard and wing configuration

圖9為15°迎角時(shí)FC布局中機(jī)翼20% 翼根弦長位置上的截面流線圖?？梢钥闯?此時(shí)鴨翼渦繞至機(jī)翼前緣渦上方,由于兩渦旋轉(zhuǎn)方向相反,相互切洗,能量有所增強(qiáng),形成兩渦之間的有利干擾。

圖9　截面流線圖Fig.9　Streamlines on cross section

圖10為不同迎角下FC布局的空間流線圖。可以看出:21°迎角時(shí),鴨翼升力系數(shù)達(dá)到最大,鴨翼渦渦核已經(jīng)破裂至前緣且渦核出現(xiàn)倒流,而單獨(dú)鴨翼渦依然保持穩(wěn)定,渦核沒有出現(xiàn)倒流,這是因?yàn)轼喴頊u在發(fā)展的過程中遇到機(jī)翼前緣渦的不利干擾,但當(dāng)鴨翼渦繞到機(jī)翼前緣渦上方時(shí),兩渦相互切洗,使得鴨翼渦外側(cè)氣流能量有所提高,所以整個(gè)鴨翼渦并沒有迅速破裂,還在一定程度上保持穩(wěn)定;迎角增大到23°時(shí),鴨翼渦外側(cè)氣流能量降低,整個(gè)鴨翼渦破裂,升力系數(shù)下降。

圖10　FC布局的空間流線圖Fig.10　Streamlines of FC configuration

3.6帶邊條后掠翼鴨式布局

對(duì)于FSC布局,鴨翼渦也受到兩方面的影響:一方面鴨翼渦外側(cè)氣流受到邊條渦的卷繞,能量有所提高;另一方面從鴨翼渦渦核發(fā)展來的氣流受到機(jī)翼前緣渦及繞機(jī)翼前緣渦流動(dòng)的氣流的上洗作用,阻礙鴨翼渦的發(fā)展,導(dǎo)致鴨翼渦能量下降。兩方面的共同作用使得鴨翼渦最大升力系數(shù)有所提高,而失速迎角基本保持不變。

圖11為不同迎角下FSC布局的空間流線圖?？梢钥闯?21°迎角時(shí),鴨翼升力系數(shù)最大,邊條卷繞著鴨翼渦外側(cè)氣流,增大了鴨翼渦能量,而鴨翼渦中心的氣流則受到后面氣流上洗的不利影響;迎角增大到23°時(shí),鴨翼渦渦核出現(xiàn)倒流,升力系數(shù)開始下降。

4結(jié)束語

本文通過數(shù)值計(jì)算,對(duì)比研究了前后掠鴨式布局中鴨翼渦的發(fā)展過程,分析了機(jī)翼前緣渦對(duì)鴨翼渦的干擾作用,得到以下結(jié)論:

(1)在后掠翼鴨式布局中,鴨翼渦在大迎角時(shí)受到機(jī)翼前緣渦的有利干擾,增大了鴨翼的升力系數(shù),提高了失速迎角,尤其是帶邊條后掠翼鴨式布局的效果非常明顯。

(2)在前掠翼鴨式布局中,鴨翼渦在大迎角時(shí)受到有利和不利兩方面的干擾,鴨翼的最大升力系數(shù)有所提高,失速迎角基本不變。

參考文獻(xiàn)：

[1]Er-El J,Seginer A.Vortex trajectories and break down on wing-canard configurations [J].Journal of Aircraft,1985,22(8):641-648.

[2]Er-El J.Effect of wing/canard interference on the loading of a delta wing [J].Journal of Aircraft,1988,25(1):18-24.

[3]Richard M H,John F O.Flow field study of a closed-coupled canard configuration [J].Journal of Aircraft,1994,31(4):908-914.

[4]劉沛清,魏園.在近距耦合鴨式布局中的渦系結(jié)構(gòu)[J].實(shí)驗(yàn)流體力學(xué),2005,19(3):85-89.

[5]Behrbohm H.Basic low speed aerodynamic of short-coupled canard configuration of small aspect ratio[R].SAAB,TN-60,Linko-ping,Sweden,1965.

[6]馬寶峰,鄧學(xué)鎣,劉沛清.近距耦合鴨式布局氣動(dòng)研究進(jìn)展[J].空氣動(dòng)力學(xué)學(xué)報(bào),2012,38(7):320-329.

[7]于沖,王旭,董福安.y+值對(duì)翼型氣動(dòng)參數(shù)計(jì)算精度的影響研究[J].空軍工程大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,13(3):25-29.

[8]王福軍.計(jì)算流體動(dòng)力學(xué)分析[M].北京:清華大學(xué)出版社,2004:31-35.

[9]劉沛清,王亞平,劉杰,等.近距耦合鴨式布局復(fù)雜渦系的干擾機(jī)理[J].北京航空航天大學(xué)學(xué)報(bào),2003,21(3):873-881.

[10]李人憲.有限體積法基礎(chǔ)[M].北京:國防工業(yè)出版社,2008:67-73.

(編輯：崔立峰)

Flow mechanism of canard vortex over forward and backward-swept wing configurations

ZHANG Dong, HU Meng-quan, WANG Xu, WU Zhang-yuan

(Aeronautics and Astronautics Engineering College, AFEU, Xi’an 710038, China)

Abstract:Based on simplified models of forward and backward-swept wing configuration, numerical investigation on the vortical structure was carried out by solving Reynolds averaged Navier-Stokes (N-S) equations, and obtained the lift coefficient curves of canard on different configurations. The characteristics of the vortex over single canard and the vortex interaction mechanisms of canard vortex and main wing vortex were analyzed by streamlines. The analysis indicates that the main wing vortex would perform a favorable impact on the canard vortex at high AOA region for backward-swept wing configuration,which enhanced the lift coefficient and increased stall angle of attack. The lift coefficient is also enhanced for forward-swept wing configuration, while the stall angle of attack remains unchanged.

Key words:canard configuration; canard vortex; flow mechanism

中圖分類號(hào)：V211.3

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1002-0853(2016)01-0036-04

作者簡介:張冬(1990-)，男，四川遂寧人，碩士，研究方向?yàn)轱w行器設(shè)計(jì)及氣動(dòng)仿真。

基金項(xiàng)目:國家自然科學(xué)基金資助(11402301)；航空科學(xué)基金資助(20110596007)

收稿日期:2015-05-25;

修訂日期:2015-09-21; 網(wǎng)絡(luò)出版時(shí)間:2015-09-29 12:56