陳 蘭,王晉軍,左林玄,馮立好

(1.北京航空航天大學(xué)流體力學(xué)教育部重點(diǎn)實(shí)驗(yàn)室,北京 100191;2.中國空氣動力研究與發(fā)展中心,四川 綿陽 621000)

0 引 言

近年來,隨著微型飛行器、無人飛行器等新型氣動布局設(shè)計(jì)發(fā)展的需要,人們的研究興趣開始集中到對非細(xì)長三角翼氣動特性和流動結(jié)構(gòu)的研究。有關(guān)實(shí)驗(yàn)和數(shù)值模擬研究表明,中小后掠三角翼在低雷諾數(shù)下呈現(xiàn)出復(fù)雜的流動現(xiàn)象,與傳統(tǒng)的細(xì)長三角翼的繞流不完全相同。

在低雷諾數(shù)(Re=4.3×103～3.47×104)條件下,Taylor等[1]在水洞中對50°后掠三角翼進(jìn)行了染色液流動顯示和PIV測量研究,分析了雷諾數(shù)、攻角和滾轉(zhuǎn)角對前緣渦的形成和破裂的影響。他們的研究表明,非細(xì)長三角翼的旋渦形成與雷諾數(shù)密切相關(guān)。實(shí)驗(yàn)和相關(guān)計(jì)算結(jié)果[2]發(fā)現(xiàn),在Re=1.3×104、α=7.5°時,前緣主渦非?？拷砻?致使邊界層與前緣分離的剪切層之間產(chǎn)生相互作用。當(dāng)二次分離流從物面離開,撞擊到主剪切層時,將主渦分割為兩個獨(dú)立的同向集中渦,形成了獨(dú)特的雙渦結(jié)構(gòu)。

魯素芬等[3-4]針對50°后掠三角翼進(jìn)行了實(shí)驗(yàn)研究,包括低速風(fēng)洞測力、測壓和流動顯示以及前緣剖面形狀、模型厚度等對三角翼氣動特性和前緣渦流動結(jié)構(gòu)的影響和規(guī)律。氫氣泡流動顯示實(shí)驗(yàn)雷諾數(shù)為1.2×104,在較小攻角下,觀察到了雙渦現(xiàn)象。Wang和Zhang[5]實(shí)驗(yàn)研究了不同后掠角(Λ=45°～65°)三角翼的前緣渦空間流動結(jié)構(gòu),發(fā)現(xiàn)在低雷諾數(shù)、小攻角下所有實(shí)驗(yàn)?zāi)Ｐ投即嬖陔p渦結(jié)構(gòu),并且出現(xiàn)雙渦結(jié)構(gòu)時模型后掠角比Gursul[6]的結(jié)果要大10°。此外,Wang和Zhang[5]指出,后掠角和雷諾數(shù)影響出現(xiàn)雙渦流動結(jié)構(gòu)的攻角及渦結(jié)構(gòu)的發(fā)展和演化。

Ol和 Gharib[7]采用體視 PIV,在水洞中對 Λ=50°的三角翼進(jìn)行了速度場的測量和流動顯示實(shí)驗(yàn),Re=6.0×103～ 1.5×104,攻角范圍 α=2.5°～20°。通過與 Λ=65°三角翼的對比,揭示了中小后掠三角翼的流動特征,如較早發(fā)生前緣渦破裂;靠近頂端的流動呈錐形流特點(diǎn);伴隨著穩(wěn)定的前緣主渦,其渦核的軸向速度剖面均為尾跡型分布,而不是通常在細(xì)長翼研究中觀察到的射流型速度分布。

目前的研究顯示,相對于細(xì)長三角翼繞流,低雷諾數(shù)下非細(xì)長三角翼繞流具有一些特有的結(jié)構(gòu),人們對此的認(rèn)識還相當(dāng)有限。本文通過流動顯示實(shí)驗(yàn)和數(shù)值模擬的方法對50°后掠三角翼繞流進(jìn)行了研究,重點(diǎn)探討了在低雷諾數(shù)下,穩(wěn)定的前緣雙渦結(jié)構(gòu)流動現(xiàn)象以及渦破裂特征。

1 模型與方法

實(shí)驗(yàn)?zāi)Ｐ蜑榧馇熬壍钠桨迦且?后掠角為50°,根弦長為150mm,厚度為2mm,前緣與后緣的迎風(fēng)面均為45°倒角。模型材料采用硬鋁,三角翼模型表面噴涂亞光漆,以增強(qiáng)流動顯示效果。模型示意見圖1(a)。

圖1 模型示意圖Fig.1 Schematic of the delta wing

氫氣泡流動顯示實(shí)驗(yàn)是在北京航空航天大學(xué)自循環(huán)水槽中進(jìn)行的,水槽觀察窗為有機(jī)玻璃。實(shí)驗(yàn)段長4.8m,截面0.6m×0.6m,水流速度在20cm/s內(nèi)可調(diào),來流湍流度小于 1%。本實(shí)驗(yàn)的來流速度為8.4cm/s,基于三角翼根弦長的雷諾數(shù)Re=1.2×104。

本文的數(shù)值模擬采用完全氣體的NS方程,有限體積法離散,隱式方法求解。方程中的無粘通量項(xiàng)采用Roe的矢通量分裂格式離散,粘性部分采用中心格式[8]。假設(shè)流動對稱,各算例都采用半模、層流計(jì)算,如圖1(b)。

2 結(jié)果及分析

2.1 流動現(xiàn)象的觀察

氫氣泡流動顯示實(shí)驗(yàn)的攻角范圍為α=0°～18°,用數(shù)字?jǐn)z像機(jī)記錄流動現(xiàn)象隨攻角的變化。從錄像中可以觀察到,當(dāng)α=2°時,三角翼上沒有出現(xiàn)任何集中渦(圖2a)。α=5°時,一對穩(wěn)定的前緣渦從前緣拖出,上翼面無渦破裂(圖2b)。α=7°時,在前緣渦的外側(cè)出現(xiàn)新的集中渦,旋轉(zhuǎn)方向與內(nèi)側(cè)的主渦一致,本文稱為外渦,由此可見形成了雙渦結(jié)構(gòu)[1-3,5-6,9]。攻角繼續(xù)增大到α=8°時,外渦明顯增強(qiáng),同時主渦在上翼面破裂。α=11°時,雙渦結(jié)構(gòu)依然存在,且靠近內(nèi)側(cè)的主渦軸向速度明顯快于外渦。另外,在上翼面中下游,外渦的流動被卷入主渦(圖2c)。攻角增大到α=15°時,外渦完全消失,只剩下一對前緣主渦,并出現(xiàn)螺旋型渦破裂(圖2d)。從各個攻角下的流動顯示來看,左右前緣渦基本對稱,主渦破裂位置隨攻角增大而不斷前移。

2.2 流動結(jié)構(gòu)的分析

本文對圖1(b)所示的三角翼繞流進(jìn)行了數(shù)值模擬研究。自由來流馬赫數(shù)M∞=0.3,雷諾數(shù)Re=1.3×104,攻角 α=2°～15°。隨著攻角增大,非定常流動區(qū)域逐漸增大。因此,對于α＞5°的算例,采用了時間推進(jìn)的二階精度非定常計(jì)算模型。

圖2 氫氣泡流動顯示圖案Fig.2 Hydrogen bubble visualization at different α

2.2.1 旋渦結(jié)構(gòu)

圖3給出了部分攻角下的空間流線圖。從圖3(a)中可以看出,α=2°時流動從前緣分離,但沒有形成明顯的集中渦。圖3(b)則表明在 α=5°時,從前緣分離的流動在離開頂點(diǎn)不遠(yuǎn)的下游形成了清晰的前緣渦,一直延伸到后緣。圖3(c)展現(xiàn)了α=11°時的雙渦結(jié)構(gòu)。仔細(xì)觀察后發(fā)現(xiàn),部分從前緣分離的流線繞過主渦后,在主渦外側(cè)又相互纏繞而形成外渦。圖3(d)顯示 α=15°時,前緣渦破裂,外渦消失。本文的實(shí)驗(yàn)和數(shù)值模擬結(jié)果都清楚地表明,隨著攻角增大,低雷諾數(shù)中等后掠三角翼繞流經(jīng)歷了從無集中渦到穩(wěn)定前緣渦,再到雙渦結(jié)構(gòu),最后蛻變?yōu)閱我磺熬墱u結(jié)構(gòu)的演化過程。

圖3 空間流線圖Fig.3 Streamlines visualization of vortex structure

表1給出了 x/c=0.5橫截面處,用當(dāng)?shù)匕胝归L無量綱化的前緣渦渦核的位置,ye是當(dāng)?shù)匕胝归L?？梢婋S著攻角的增大,渦核遠(yuǎn)離三角翼上表面,同時向三角翼對稱面移動,其變化規(guī)律與文獻(xiàn)[5]一致。

表1 渦核位置Table1 Vortex core location at x/c=0.5

α=5°時的表面極流線如圖4所示。流動從前緣分離,進(jìn)而再附,然后發(fā)生二次分離和二次再附,圖中的 S1、A1、S2、A2 分別表示主分離線、主再附線 、二次分離線和二次再附線。圖中顯示,大約從15%根弦長的位置才開始逐步形成二次渦。計(jì)算表明,隨著攻角增大,二次渦形成的起始點(diǎn)隨攻角增加不斷向頂點(diǎn)前移。α=10°時,二次渦起始點(diǎn)前移到6%根弦長的位置。當(dāng)α=15°時,該點(diǎn)已前移至4%根弦長處。表2列出了x/c=0.3、0.5截面處的主再附點(diǎn)(PA)、二次分離點(diǎn)(SS)、二次再附點(diǎn)(SA)的相對展向位置。由表可見,隨著攻角的增加,二次分離點(diǎn)向?qū)ΨQ面移動,而二次再附點(diǎn)向前緣移動,二次分離區(qū)逐步擴(kuò)大。隨攻角增大,主再附點(diǎn)不斷向內(nèi)側(cè)中心線移動。

圖 4 表面極流線(α=5°)Fig.4 Surface streamline pattern(α=5°)

表2 分離點(diǎn)、再附點(diǎn)位置Table2 Separation and attachment locations

2.2.2 雙渦現(xiàn)象

本文以α=10°為例,分析雙渦結(jié)構(gòu)的流動特征。以前緣渦渦軸為法向,如圖1(b)所示分別在 x/c=0.05、0.10和0.30處截取三個平面。圖5繪制出渦量在渦軸方向分量的等值線圖,清晰地給出了主渦、二次渦、外渦沿流向的生成發(fā)展以及相對位置。由圖可見,在 x/c=0.05截面處主渦已經(jīng)生成(圖5a)。在x/c=0.10處(圖5b),微弱的二次渦出現(xiàn)在主渦外側(cè)下方,與主渦旋轉(zhuǎn)方向相反。在 x/c=0.30處(圖5c),可以清晰地看出在主渦外側(cè)、二次渦上方,又形成一個與主渦同向旋轉(zhuǎn)的外渦,即出現(xiàn)了雙渦結(jié)構(gòu),外渦具有與主渦相同的渦量符號,且比主渦弱。

圖6給出的是靠近頂點(diǎn)處的空間流線局部圖,可見主渦幾乎從三角翼頂端發(fā)出,而外渦卻是在主渦形成后的下游再形成,圍繞外渦渦核的流線相互纏繞且較主渦松散,氫氣泡實(shí)驗(yàn)也觀察到相同的現(xiàn)象。為跟蹤旋渦的生成,繪制空間流線時在前緣剪切層內(nèi)布點(diǎn),發(fā)現(xiàn)流動從前緣分離,靠近剪切層底層的大部分流線被卷入主渦,部分中間層流線卷入外渦,少量最外層流線進(jìn)入二次渦。

圖5 三角翼的旋渦結(jié)構(gòu)(α=10°)Fig.5 Vortex structure over delta wing at α=10°

圖6 空間流線局部圖Fig.6 Detail of streamlines

圖7、圖8為數(shù)值模擬得到的渦核在翼面的投影和距離翼面的高度。與主渦不同的是外渦不是在三角翼的頂點(diǎn)形成,而是在其下游外側(cè)形成。另外,主渦渦核在渦破裂之前近似為直線,而外渦渦核略顯彎曲。外渦形成后,離物面高度稍低于主渦。隨著向下游發(fā)展,兩個渦都在不斷遠(yuǎn)離模型上表面。在弦長x/c≈0.50～0.70范圍內(nèi),從渦核高度上看,外渦出現(xiàn)了波動。因此,由圖7和圖8可知,外渦沒有主渦穩(wěn)定,在渦核的空間分布上呈現(xiàn)波動狀。

圖7 渦核在翼面的投影Fig.7 Projection of vortex core

圖8 渦核距離翼面Fig.8 Height of vortex core

圖9比較了雙渦的渦軸速度。主渦破裂前,主渦渦核的軸向速度遠(yuǎn)遠(yuǎn)高于外渦,與氫氣泡流動顯示實(shí)驗(yàn)觀察到的現(xiàn)象一致,主渦軸向運(yùn)動速度大于外渦。一般而言,若在三角翼背風(fēng)面上形成前緣渦,其對物面壓力分布的誘導(dǎo)作用,將產(chǎn)生吸力峰值,由此產(chǎn)生了渦升力,進(jìn)而可導(dǎo)致升力的增加。圖10給出了x/c=0.3截面處壓力系數(shù)分布,圖中虛線表示該截面的物面壓力分布,實(shí)線代表了該截面離物面距離為渦核高度z/c=0.0257處的壓力分布。在主渦渦核y/ye=0.536處,出現(xiàn)負(fù)壓峰值。外渦渦核在y/ye=0.816附近也有負(fù)壓的抬升,但相比于主渦要小得多。主渦在物面上引起了明顯的吸力峰值而外渦并沒有引起物面的壓力系數(shù)出現(xiàn)第二個峰值,因此我們得出結(jié)論:外渦對升力的貢獻(xiàn)不顯著。

圖9 雙渦的渦核軸向速度分布Fig.9 Axial velocity of primary and outer vortex

圖10 壓力分布圖Fig.10 Pressure distributions

2.2.3 渦破裂的特征

一般認(rèn)為,渦破裂有兩種形式,一是泡型破裂,一是螺旋型破裂[10]。計(jì)算中發(fā)現(xiàn)了這兩種破裂方式,圖11分別展示了8°攻角時的泡型破裂和10°攻角時的螺旋型破裂。

圖11 前緣主渦的破裂Fig.11 Streamline visualization of vortex breakdown

圖12比較了各攻角下前緣主渦的渦軸速度。根據(jù)圖中曲線的變化趨勢,渦核軸向流動可分為三個階段:起始階段的加速區(qū),達(dá)到峰值后的緩慢減小區(qū),進(jìn)入急劇減速區(qū)。攻角越大,從前緣分離的剪切層卷起形成集中渦的起始點(diǎn)越接近三角翼頂端,渦核最大軸向速度也越大。前緣渦沒有破裂的5°攻角,旋渦加、減速的梯度都很平緩,一直到后緣都保持為近似平臺。其它發(fā)生渦破裂的攻角,都有沿軸向流動的急劇減速區(qū),攻角越大,減速區(qū)開始得越早,渦破裂也越提前。例如,α=7°時,渦破裂點(diǎn)位于上翼面后緣。α=10°時,渦軸速度從 x/c=0.4開始迅速下降,在 x/c=0.65附近破裂。實(shí)驗(yàn)和計(jì)算表明,破裂方式的不同,在破裂點(diǎn)附近,渦軸速度沿軸向的變化也有所不同。如圖12所示,α=9°和10°時為螺旋型破裂,其破裂點(diǎn)的軸向速度大于零,而α=8°時則呈現(xiàn)為泡型破裂,軸向速度一直降低到過零。數(shù)值計(jì)算結(jié)果與文獻(xiàn)[10]中關(guān)于旋渦沿其軸線的分叉演化理論相符。

圖12 沿渦核軸向的速度分布Fig.12 Axial velocities along core

4 結(jié) 論

本文通過水槽流動顯示和數(shù)值模擬,對低雷諾數(shù)下50°后掠三角翼的流動特點(diǎn)進(jìn)行了初步研究,結(jié)論如下:

(1)50°后掠三角翼在一定攻角范圍內(nèi)存在明顯的雙渦結(jié)構(gòu),外渦的軸向渦量比主渦小、渦核軸向速度比主渦慢,其對翼面壓力系數(shù)的影響有限,對升力的貢獻(xiàn)不明顯。

(2)隨著攻角的增大,前緣主渦的渦核位置不斷升高,主再附線向三角翼中心線移動;低雷諾數(shù)下,流動的二次分離和再附是在三角翼頂點(diǎn)下游逐漸形成的,并非從頂點(diǎn)產(chǎn)生,隨著攻角增加,其形成起始點(diǎn)不斷向頂點(diǎn)前移,二次分離區(qū)也逐漸擴(kuò)大。

(3)該三角翼在低雷諾數(shù)下發(fā)生前緣主渦破裂的起始攻角較小。隨攻角增大,破裂點(diǎn)提前。另外,泡型和螺旋型破裂形式都可能出現(xiàn)。

[1]TAYLOR G,SCHNORBUS T,GURSUL I.An investigation of vortex flows over low sweep delta wings[R].AIAA 2003-4021.

[2]GORDNIER R,VISBAL M.Higher-order compact difference scheme applied to the simulation of a low sweep delta wing flow[R].AIAA 2003-620.

[3]魯素芬.50°后掠角三角翼流動結(jié)構(gòu)及氣動特性實(shí)驗(yàn)研究[D].[碩士學(xué)位論文].北京航空航天大學(xué).2005.(LU S F.Experimental investigation on the vortex structures and aerodynamics for flow over 50°delta wings[D].[Master's Thesis].Beijing University of Aeronautics and Astronautics,2005.)

[4]WANG J J,LU S F.Effects of leading-edge bevel angle on the aerodynamic forces of a non-slender 50°delta wing[J].Aeronautical Journal,2005,109(1098):403-407.

[5]WANG J J,ZHANG W.Experimental investigations on leading-edge vortex structures for flow over non-slender delta wings[J].Chinese Physical Letters,2008,25(7):2550-2553.

[6]GURSUL I.Recent developments in delta wing aerodynamics[J].The Aeronautical J ournal,2004,108:437-452.

[7]MICHAEL V Ol,MORTEZA G.Leading-edge vortex structure of non-slender delta wings at low Reynolds number[R].AIAA 1930-656.

[8]周乃春.非圓截面機(jī)身氣動特性研究[D].[博士學(xué)位論文].中國空氣動力研究與發(fā)展中心.2003.

[9]WANG J J,ZHAO X,LIU W C,T U J Q.Experimental investigation on flow structures over non-slender delta wings at low Reynolds numbers[J].Journal of Experiments in Fluid Mechanics,2007,21(2):1-7.

[10]張涵信.分離流與旋渦運(yùn)動的結(jié)構(gòu)分析[M].國防工業(yè)出版社,2005.(ZHANG H X.Structural analysis of separated flows and vortex motion[M].Nation Defence Industry Press,2005.)