跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究

錢麗麗　鄧長喜

【摘 要】飛機(jī)在大氣中跨音速飛行時，機(jī)翼附近的空氣壓力會發(fā)生劇烈的變化，引起局部飽和水汽壓的波動，使其中含有的水蒸氣可能高于飽和線而發(fā)生凝結(jié)，這會對機(jī)翼的空氣動力學(xué)特性造成影響。本文以NACA0012翼型為研究，建立在濕空氣中飛行時的氣流及水汽凝結(jié)的數(shù)學(xué)模型，分析了在跨音速飛行時，攻角α為0°、2°、4°，空氣相對濕度φ為20%、40%、60%、80%的情況下水汽凝結(jié)情況進(jìn)行了模擬分析，得出在上述幾種狀態(tài)下的水霧分布情況，為針對水霧對機(jī)翼空氣動力學(xué)特性影響的研究提供了一定的依據(jù)。

【關(guān)鍵詞】機(jī)翼；跨音速；凝結(jié)；數(shù)值模擬

Numerical Simulation Research on Condensation of Water Vapor Around Airfoil in Transonic Flying

QIAN Li-li DENG Chang-xi

（Avic Xian Aircraft Industry（Group）Company， Yanliang Shaanxi 710089， China）

【Abstract】The air pressure around airfoils would be changing severely during transonic flying in atmosphere， that would vary the partial saturation vapor pressure lower to make the vapor carried by moist air condense for it may expand across the saturation line， which could affect the aerodynamic characteristics of airfoils. The numerical model of vapor condensation and airflow during flying in the moist air about the NACA0012 airfoil was established， the situations of vapor condensation in transonic flying under the angle of attack were 0°， 2° and 4°， and the relative humidity were 20%， 40%， 60% and 80% were analysed， the distributions of vapor under situations above were got， which provide the certain basis to the studies about water mist effect on airfoil aerodynamic characteristics.

【Key words】Airfoil； Transonic； Condensation； Numerical simulation

空氣作為一種最普遍的氣體介質(zhì)，除了在少數(shù)實驗室條件下，其中通常會含有一定量的水蒸氣。飛機(jī)在大氣中飛行時，空氣中的水蒸氣在未發(fā)生相變即始終為氣體狀態(tài)的情況下，其影響通?？梢院雎?。但在跨音速飛行時，由于機(jī)身周圍局部氣流膨脹，氣壓降低，飽和水汽壓減小，使得水蒸氣高于飽和線發(fā)生析出、凝結(jié)，在機(jī)身周圍形成水霧，俗稱“音爆云”。圖1為殲10戰(zhàn)斗機(jī)機(jī)身周圍出現(xiàn)的水霧。

圖1 殲10機(jī)身周圍水汽凝結(jié)現(xiàn)象

這種水汽凝結(jié)現(xiàn)象不僅會改變機(jī)身附近氣流的溫度、壓力等特征，還會對機(jī)翼的升阻特性等空氣動力學(xué)性質(zhì)有較大的影響。另一方面，跨音速飛行中水汽的凝結(jié)現(xiàn)象不僅伴隨相變、傳熱等情況，同時凝結(jié)潛熱產(chǎn)生的凝結(jié)沖波還會與激波、邊界層之間相互影響，流動過程非常復(fù)雜，相關(guān)實驗研究難度較大[1-2]。因此，對該現(xiàn)象在不同飛行條件下的產(chǎn)生、分布情況進(jìn)行數(shù)值模擬研究，有著重要的實際意義。

本文以NACA0012型機(jī)翼為研究對象，對其在不同攻角、不同空氣相對濕度條件下跨音速飛行時，機(jī)身周圍產(chǎn)生的水霧分布情況進(jìn)行了分析。

1 水蒸氣凝結(jié)的數(shù)學(xué)模型

水蒸氣與干空氣組成的混合氣體稱為濕空氣，在流動中可將其視為一個兩相流場。在二維無黏可壓縮流動情況下，氣相流動的控制方程為[3-4]：

其中，ρ為濕空氣與霧?；旌衔锏拿芏龋裧為氣相密度，β為霧粒在混合物中的質(zhì)量分?jǐn)?shù)，p為濕空氣的靜壓，ρ為其密度，e為其總能，h0為總焓，h為混合物焓，hg為水蒸氣焓，hl為霧粒焓。

當(dāng)?shù)乜諝庵泻瑵窳縟與空氣溫度t間的關(guān)系如下述經(jīng)驗公式[5]：

式中，Pb為飽和水汽壓，P為當(dāng)?shù)卮髿鈮?，φ為空氣相對濕度?/p>

2 機(jī)翼水汽凝結(jié)模擬分析

根據(jù)前述數(shù)學(xué)模型，本文使用Fluent軟件對NACA0012型機(jī)翼在大氣跨音速飛行時的水汽凝結(jié)分布情況進(jìn)行了數(shù)值模擬。

建立NACA0012翼型橫截面二維幾何模型，采用非結(jié)構(gòu)網(wǎng)格劃分計算區(qū)域，節(jié)點總數(shù)為19660。計算中的流場條件為：機(jī)翼弦長c=1.0m，當(dāng)?shù)販囟萾=298k，當(dāng)?shù)卮髿鈮簽闃?biāo)準(zhǔn)大氣壓P=101325Pa，馬赫數(shù)Ma=0.8，攻角α分別取0°、2°、4°，相對濕度φ分別取20%、40%、60%、80%，遠(yuǎn)場邊界與機(jī)翼的距離為15倍弦長。

跨音速飛行時，在不同攻角情況下，機(jī)翼附近氣流馬赫數(shù)分布情況如圖2所示。

由圖2可看出，在零升力即攻角α=0°時，機(jī)翼兩側(cè)馬赫數(shù)為對稱分布，隨著攻角的增大，上翼面馬赫數(shù)變化情況明顯較下翼面更劇烈。

在濕空氣中跨音速飛行過程中，機(jī)翼附近局部空氣飽和水汽壓降低，水蒸氣越過飽和線形成水霧。

當(dāng)飛行攻角α及空氣相對濕度φ發(fā)生變化時，機(jī)翼附近水霧分布情況相應(yīng)發(fā)生變化，如圖3所示。

（a）攻角α=0°

（b）攻角α=2°

（c）攻角α=4°

圖3 機(jī)翼附近水霧分布情況

由圖3可以看出，在α=0°時水霧在機(jī)翼附近以翼弦為中心呈對稱分布，主要集中在機(jī)翼前半部分，且隨著空氣相對濕度的升高，水霧含量相應(yīng)增加。

當(dāng)α增加至2°和4°時，機(jī)翼上下側(cè)水霧出分呈明顯的不對稱。隨著α的逐漸增大，水霧含量高的區(qū)域越集中與翼面前半部的上側(cè)空間內(nèi)。

3 機(jī)翼表面含濕量分布分析

考慮到水霧越靠近機(jī)翼表面，對機(jī)翼氣動特性及飛行狀態(tài)的影響越直接，因此本文以相對濕度φ=80%的情況為例，對機(jī)翼表面的水霧分布情況進(jìn)行了分析，如圖4、圖5所示。

圖4 機(jī)翼上表面水霧含量分布

圖5 機(jī)翼下表面水霧含量分布

在圖4、圖5中，縱軸0點位置可認(rèn)為是水蒸氣飽和線，高于此點則水蒸氣凝結(jié)為水霧，低于此點則水蒸氣仍為氣態(tài)。

由圖可知，在機(jī)翼前緣附近空氣隨沿弦向快速逼近飽和線，在距前緣0.1m處附近開始，上翼面逐漸開始發(fā)生凝結(jié)現(xiàn)象，出現(xiàn)水霧，隨后水霧含量繼續(xù)增加，在機(jī)翼中部即距前緣0.5m附近開始快速下降，將至飽和線以下時，水霧變回氣態(tài)。

在上翼面，不同攻角情況下水霧分布趨勢基本一致，但是攻角越大，水霧最前端距機(jī)翼前緣越近；α=2°時，水霧在上翼面覆蓋率最大，α=4°時，水霧在上翼面前半部最集中，即覆蓋率最小；α=0°時水霧的含量峰值相對最小。

在α=0°時，由于水霧分布的對稱性，下翼面的水霧分布與上翼面幾乎無差別；α=2°及4°時，機(jī)翼下表面空氣含濕量均未越過飽和線，水蒸氣未出現(xiàn)凝結(jié)而形成水霧。

4 結(jié)論

本文以NACA0012翼型為研究對象，建立了其在濕空氣中飛行時氣流運動及水蒸氣凝結(jié)的數(shù)學(xué)模型。對機(jī)翼在Ma=0.8時，在不同攻角及空氣相對濕度情況下的跨音速飛行中水汽凝結(jié)情況進(jìn)行了模擬，得出了機(jī)翼附近產(chǎn)生水霧的分布情況，為針對水霧對機(jī)翼空氣動力學(xué)特性影響的研究提供了一定的依據(jù)。

跨音速飛行時濕空氣中的水汽凝結(jié)情況非常復(fù)雜，包含相變、傳熱以及霧粒蒸發(fā)等物理過程，同時會對激波、邊界層等產(chǎn)生影響，若考慮到懸浮顆粒等影響，將使模型更加復(fù)雜，因此很有必要針對此問題進(jìn)行進(jìn)一步的實驗研究及模型完善等工作。

【參考文獻(xiàn)】

[1]Hill P G. Condensation of water vapor during supersonic expansion in nozzles[J]. J Fluid Mech， 1966， 25（3）： 593-620.

[2]YOUNG J B.Two-dimensional， nonequilibrium， wet-steam calculations for nozzles and turbine cascades[J]. ASMEJ Turbomach， 1992， 114：569-579.

[3]孫秀玲，李亮，李國君.濕空氣凝結(jié)對ONERA M6機(jī)翼氣動特性影響的數(shù)值研究[J].空氣動力學(xué)學(xué)報，2010，02：143-148.

[4]楊勇，沈勝強，孔泰佑.水蒸氣跨音速流動中非平衡現(xiàn)象的數(shù)值模擬[J].熱科學(xué)與技術(shù)，2008，03：197-202.

[5]于新峰.基于雙流體模型的濕蒸汽兩相流動數(shù)值模擬[D].哈爾濱工業(yè)大學(xué)，2011.

[6]蔡頤年，王乃寧.濕蒸汽兩相流[M].西安交通大學(xué)出版社，1985.

[責(zé)任編輯：楊玉潔]