劉 燕，岳巍澎

（1.華北電力大學(xué)，能源動(dòng)力與機(jī)械工程學(xué)院，北京 102206；2.國(guó)網(wǎng)新源張家口風(fēng)光儲(chǔ)示范電站有限公司，河北張家口 075000）

海上濕氣對(duì)翼型氣動(dòng)特性的影響初探

劉 燕1，＊，岳巍澎1，2

（1.華北電力大學(xué)，能源動(dòng)力與機(jī)械工程學(xué)院，北京 102206；2.國(guó)網(wǎng)新源張家口風(fēng)光儲(chǔ)示范電站有限公司，河北張家口 075000）

研究高濕度空氣對(duì)風(fēng)力機(jī)翼型氣動(dòng)性能影響對(duì)海上風(fēng)力機(jī)翼型優(yōu)化、提高發(fā)電效率具有重要意義。首先采用組分輸運(yùn)模型數(shù)值計(jì)算濕空氣條件下翼型升阻力特性，其次在組分輸運(yùn)模型基礎(chǔ)上，考慮濕空氣中水蒸氣的凝結(jié)，建立用戶自定義函數(shù)定義凝結(jié)體積源項(xiàng)，求解輸運(yùn)方程。對(duì)濕空氣流過(guò)翼型表面氣動(dòng)性能的初步研究發(fā)現(xiàn)：升阻力系數(shù)隨濕度變化而改變，不考慮凝結(jié)時(shí)，升阻力系數(shù)隨濕度增大而減小，當(dāng)發(fā)生凝結(jié)時(shí)，升力系數(shù)增加。翼型表面前緣點(diǎn)附近凝結(jié)流量最大，翼型表面凝結(jié)流量隨濕度增大而增大，溫度越高，凝結(jié)流量越大。

翼型；海上風(fēng)力機(jī)；濕空氣；氣動(dòng)特性；凝結(jié)

0 引 言

由于全球氣候變暖和能源危機(jī)，傳統(tǒng)能源模式已難以為繼，加快新能源開(kāi)發(fā)利用勢(shì)在必行。海上風(fēng)電容量大、效率高，已成為全球風(fēng)電發(fā)展的最前沿方向。與此同時(shí)，惡劣的海洋環(huán)境制約著海上風(fēng)力機(jī)的發(fā)展，強(qiáng)烈的日光照射、高濃度的鹽霧環(huán)境、高濕度、水分侵蝕等因素都會(huì)對(duì)葉片表面造成損害，影響機(jī)組安全，增加運(yùn)行維護(hù)成本［1－3］。

海洋環(huán)境濕度較大，我國(guó)渤海灣長(zhǎng)年濕度在50%－95%之間，翼型及葉片周?chē)鷫毫?、溫度的變化在一定條件下會(huì)引起水蒸汽發(fā)生凝結(jié)現(xiàn)象。凝結(jié)不僅會(huì)顯著改變翼型表面附近的壓力、溫度等流動(dòng)參數(shù)，甚至還會(huì)引起穩(wěn)定或振蕩的激波，導(dǎo)致翼型升力、阻力特性發(fā)生改變，影響其氣動(dòng)性能［4－5］。此外，葉片表面結(jié)冰問(wèn)題也與濕空氣及凝結(jié)現(xiàn)象緊密相關(guān)，濕空氣凝結(jié)是形成冰晶的重要機(jī)制之一［6］。

國(guó)內(nèi)外許多研究機(jī)構(gòu)和學(xué)者對(duì)濕空氣凝結(jié)流動(dòng)問(wèn)題開(kāi)展研究。NASA實(shí)驗(yàn)研究了濕空氣凝結(jié)對(duì)機(jī)翼氣動(dòng)性能的影響［7］。Schnerr研究了均勻凝結(jié)流動(dòng)中湍流邊界層和沖擊波相互作用［8］；Rusak等［9］學(xué)者理論分析了凝結(jié)起始馬赫數(shù)和凝結(jié)放熱量等方面與凝結(jié)流動(dòng)之間的相似性；Karabelas和Markato［10］數(shù)值研究了高對(duì)流條件下翼型的水蒸氣凝結(jié)造成的影響；孫秀玲、李亮等［11－12］數(shù)值模擬了跨聲速條件下RAE 2822翼型及ONERA M6機(jī)翼濕空氣非平衡凝結(jié)流動(dòng)對(duì)氣動(dòng)性能的影響。

研究表明：流動(dòng)與凝結(jié)過(guò)程相互耦合，速度、迎角、空氣濕度等參數(shù)與凝結(jié)過(guò)程相互影響，濕空氣凝結(jié)對(duì)翼型升力、阻力特性有顯著影響。盡管如此，在濕空氣凝結(jié)流動(dòng)研究方面，仍有一些重要的問(wèn)題沒(méi)有解決，溫度、壓力、相對(duì)濕度、速度及幾何條件等對(duì)凝結(jié)發(fā)生及氣動(dòng)性能的影響，目前尚沒(méi)有一般性的結(jié)論，凝結(jié)影響機(jī)翼氣動(dòng)性能的一般規(guī)律也并不清楚；且之前的研究多是針對(duì)航空翼型或機(jī)翼，針對(duì)風(fēng)力機(jī)領(lǐng)域的研究還沒(méi)有，特別是隨著海上風(fēng)電的大力發(fā)展，高惡劣海上氣候條件對(duì)風(fēng)力機(jī)葉片設(shè)計(jì)的要求更加苛刻，需要進(jìn)一步了解濕空氣對(duì)海上風(fēng)力機(jī)葉片氣動(dòng)性能、噪聲的產(chǎn)生傳播乃至結(jié)構(gòu)的影響。

本文以海上風(fēng)力機(jī)高濕度環(huán)境為背景，考慮濕空氣的物理特性及凝結(jié)過(guò)程，根據(jù)濕度不同，將濕空氣考慮為水蒸氣和空氣的混合氣體，采用輸運(yùn)模型研究翼型的升阻力，并與實(shí)驗(yàn)對(duì)比；在輸運(yùn)模型基礎(chǔ)上，考慮凝結(jié)發(fā)生，通過(guò)FLUENT編寫(xiě)表面凝結(jié)自定義函數(shù)（UDF）實(shí)現(xiàn)翼型表面凝結(jié)，研究凝結(jié)發(fā)生情況下的翼型氣動(dòng)性能。

1 模 型

常規(guī)翼型及風(fēng)力機(jī)CFD計(jì)算，以干空氣為流動(dòng)介質(zhì)?？紤]實(shí)際海上濕度環(huán)境，將濕空氣（干空氣和水蒸氣的混合氣體）作為流動(dòng)介質(zhì)，采用組分輸運(yùn)模型，更接近實(shí)際情況，計(jì)算結(jié)果較干空氣更加可靠。在輸運(yùn)模型基礎(chǔ)上，考慮濕空氣中水蒸氣的凝結(jié)，通過(guò)編寫(xiě)用戶自定義凝結(jié)函數(shù)，實(shí)現(xiàn)水蒸氣在壁面凝結(jié)。

凝結(jié)原理（圖1）：當(dāng)濕空氣流過(guò)物體表面，壁面溫度低于該流體的飽和溫度時(shí)，表面發(fā)生冷凝，形成液膜，且在流體重力和剪切力作用下流動(dòng)。凝結(jié)過(guò)程中，非凝氣體（空氣）在氣液交界面逐漸累積，阻礙蒸汽擴(kuò)散到壁面，降低了水汽的凝結(jié)速率。盡管凝結(jié)是涉及流動(dòng)、傳質(zhì)換熱的多相耦合過(guò)程，但在某些條件下，凝結(jié)速率主要取決于蒸汽相的擴(kuò)散過(guò)程。如果凝結(jié)液膜的熱阻非常小，可以認(rèn)為氣液界面的溫度與壁面溫度相同，凝結(jié)速率可以通過(guò)求解蒸汽流動(dòng)和非凝氣體的分布預(yù)測(cè)。壁面冷凝邊界條件允許蒸汽相質(zhì)量傳遞，非凝氣體為零凈傳質(zhì)速率。

圖1 蒸汽壁面凝結(jié)原理Fig.1 Water vapor condensation on the wall

基于以上原理，作以下假設(shè)：

① 冷凝發(fā)生由于水蒸汽通過(guò)組分邊界層擴(kuò)散。

② 液膜和壁面存在熱平衡。

③ 氣液界面是飽和狀態(tài)。即交界面水蒸氣分壓力等于壁面溫度下的飽和蒸汽壓力。

④ 不考慮液膜運(yùn)動(dòng)帶來(lái)的影響。由于流體剪切力和重力，液膜可能運(yùn)動(dòng)，但是運(yùn)動(dòng)不影響空氣流動(dòng)和凝結(jié)。

⑤ 分析僅限于膜狀凝結(jié)，不考慮滴狀冷凝。

⑥ 液膜的熱阻忽略不計(jì)。

空氣和水蒸氣各組分的質(zhì)量流量方程為：

其中，ω為質(zhì)量分?jǐn)?shù)，m″為質(zhì)量流量（kg／m2·s），ρ為密度（kg／m3），v為界面速度（m／s），D為質(zhì)量擴(kuò)散率（m2／s）。

由二元混合物質(zhì)量分?jǐn)?shù)定義可得：

氣液交界面混合物質(zhì)量為：

根據(jù)假設(shè)②由于液相只包含水，因此：

將公式（3）、（4）帶入（1），聯(lián)合公式（5）可得：

通過(guò)在近壁面使用質(zhì)量（連續(xù)性）源項(xiàng)來(lái)代表冷凝率，即：

UDF中，假定壁面溫度等于或小于表面水蒸氣分壓力對(duì)應(yīng)的飽和溫度時(shí)發(fā)生凝結(jié)。如果溫度高于飽和溫度，水蒸氣質(zhì)量分?jǐn)?shù)等于壁面鄰近單元的值。為了滿足氣液交界面局部熱力學(xué)平衡的假設(shè)，當(dāng)溫度低于或等于飽和溫度時(shí)，代表氣液交界面的水蒸氣質(zhì)量分?jǐn)?shù)被分配一個(gè)值，使得水蒸汽的分壓等于水在局部壁面溫度下的飽和壓力。為此，近壁面鄰近單元的組分方程中的源項(xiàng)必須包含擴(kuò)散通量。

由公式（8）結(jié)合公式（2）、（5）、（6）得：

組分方程中移除壁面鄰近單元的水汽量，等于連續(xù)方程中移除壁面鄰近單元的水汽量，即：

公式（10）減去公式（9）重排后得：

公式（11）中混合物的體積質(zhì)量源項(xiàng)為：

水蒸氣的體積質(zhì)量源項(xiàng)為：

2 結(jié)果驗(yàn)證

2.1 翼型CFD數(shù)值計(jì)算模型驗(yàn)證

為驗(yàn)證翼型數(shù)值模擬方法的可行性，以NACA63418翼型為測(cè)試翼型，計(jì)算翼型的升力系數(shù)、阻力系數(shù)、表面壓力系數(shù)等氣動(dòng)參數(shù)，并與實(shí)驗(yàn)數(shù)據(jù)對(duì)比，以驗(yàn)證本次計(jì)算的數(shù)值方法。為了便于與實(shí)驗(yàn)數(shù)據(jù)對(duì)比，網(wǎng)格劃分方法及網(wǎng)格數(shù)參考文獻(xiàn)［13］，外流場(chǎng)計(jì)算邊界取30倍弦長(zhǎng)。測(cè)試翼型為標(biāo)準(zhǔn)翼型，弦長(zhǎng)為1m，雷諾數(shù)Re＝3 000 000，常溫下空氣密度ρ＝1.225 kg／m3，粘性系數(shù)μ＝1.7894×10－5kg／m·s。

采用有限體積法求解偏微分方程，離散方程采用基于壓力的simple算法求解，對(duì)流通量選取二階迎風(fēng)格式。為了解不同湍流模型對(duì)氣動(dòng)性能的影響，分別采用k－ωSST模型、S－A流模型、transition SST（4eqn）模型計(jì)算翼型的氣動(dòng)性能，并與實(shí)驗(yàn)數(shù)據(jù)對(duì)比。

圖2為不同模型對(duì)應(yīng)的升力系數(shù)和阻力系數(shù)。對(duì)比實(shí)驗(yàn)數(shù)據(jù)可以看出，三種模型整體上模擬結(jié)果與實(shí)驗(yàn)結(jié)果吻合較好，特別是在迎角小于10°時(shí)，各模型計(jì)算結(jié)果與實(shí)驗(yàn)值一致。升力系數(shù)方面，k－ωSST模型和S－A模型，其計(jì)算值與實(shí)驗(yàn)結(jié)果吻合度最高；而阻力系數(shù)Transition SST模型與實(shí)驗(yàn)結(jié)果最接近。迎角大于10°時(shí)，各模型計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果相差較大。綜合考慮計(jì)算精度和計(jì)算時(shí)間，本次計(jì)算選取k－ωSST模型。

圖2 升力系數(shù)和阻力系數(shù)Fig.2 Lift coefficient and drag coefficient

2.2 凝結(jié)模型驗(yàn)證

濕空氣中水蒸氣的凝結(jié)在許多工程應(yīng)用中發(fā)揮重要的作用，例如空調(diào)設(shè)備、電廠換熱器、濕煙氣余熱回收以及窗戶結(jié)霧都涉及到水汽凝結(jié)。通過(guò)計(jì)算濕空氣流過(guò)水平板后凝結(jié)與Sparrow［14］實(shí)驗(yàn)結(jié)果對(duì)比，來(lái)驗(yàn)證凝結(jié)UDF函數(shù)的可用性。圖3為計(jì)算幾何模型。指定入口流體速度、溫度、蒸汽質(zhì)量分?jǐn)?shù)及壁面溫度，蒸汽質(zhì)量分?jǐn)?shù)根據(jù)相對(duì)濕度確定。通過(guò)表面凝結(jié)UDF計(jì)算氣液交界面處凝結(jié)后的水蒸氣質(zhì)量分?jǐn)?shù)。

圖3 平板幾何及邊界條件Fig.3 Geometry and boundary conditions for condensation on a horizontal flat plate

圖4是計(jì)算結(jié)果與強(qiáng)制對(duì)流凝結(jié)解析解［15］的對(duì)比。x軸為沿平板方向的距離，y軸為凝結(jié)速率（kg／m2·s）。CFD預(yù)測(cè)凝結(jié)速率除了平板前緣區(qū)之外，與解析解能夠很好的吻合。結(jié)果表明，自定義函數(shù)能夠滿足凝結(jié)邊界的求解。

圖4 CFD計(jì)算結(jié)果與解析解Fig.4 Analytical and CFD results for condensation rates

3 結(jié)果分析

為全面地理解濕空氣對(duì)翼型氣動(dòng)性能的影響，翼型的迎角計(jì)算2°、6°、10°、16°四種工況，翼型弦長(zhǎng)為1m?？紤]實(shí)際海上濕度范圍一般為55%～95%，溫度－20°～45°，本次計(jì)算選取相對(duì)濕度RH為25%、50%、70%、90%四種濕度工況，分別對(duì)應(yīng)溫度為30°和40°。濕空氣根據(jù)溫度和濕度考慮為干空氣和水蒸氣按相應(yīng)比例的混合氣體。

圖5為兩種模型計(jì)算的對(duì)應(yīng)三個(gè)迎角和兩個(gè)溫度下濕度對(duì)阻力系數(shù)的影響。由圖看出，三種迎角阻力系數(shù)隨濕度變化趨勢(shì)不完全相同。2°迎角時(shí)，阻力系數(shù)先降低后增大。10°迎角和16°迎角時(shí)，除了溫度為40°不考慮凝結(jié)的工況，其余工況阻力系數(shù)基本上隨濕度增大而減小。阻力系數(shù)較升力系數(shù)，考慮凝結(jié)和不考慮凝結(jié)，兩種模型變化趨勢(shì)基本一致。同時(shí)可以看出，溫度對(duì)阻力系數(shù)的影響，相同濕度條件時(shí)，溫度高阻力系數(shù)高。

圖6為不同濕度下翼型升力系數(shù)曲線。可以看出，不考慮凝結(jié)時(shí)，組分輸運(yùn)模型計(jì)算升力系數(shù)變化趨勢(shì)一致，隨濕度增大升力系數(shù)減小。考慮凝結(jié)時(shí)，升力系數(shù)變化均先降低，再增大。主要是兩種溫度下，相對(duì)濕度為25%時(shí)，沒(méi)有發(fā)生凝結(jié)，流動(dòng)認(rèn)為單相流動(dòng)，氣動(dòng)性能改變主要是由于濕空氣密度的變化所引起的，計(jì)算結(jié)果趨勢(shì)與輸運(yùn)模型結(jié)果趨勢(shì)一致；濕度增大到50%時(shí)，溫度30°的工況沒(méi)有發(fā)生凝結(jié)，而溫度40°的工況已發(fā)生凝結(jié)，當(dāng)濕度增大到70%以后，計(jì)算達(dá)到穩(wěn)定后，兩種溫度下均發(fā)生凝結(jié)，產(chǎn)生液體，空氣中的水蒸氣在翼型表面凝結(jié)成水膜，氣動(dòng)性能發(fā)生改變。

圖5 不同迎角阻力系數(shù)隨濕度變化曲線Fig.5 Variation of the drag coefficient as a function of humidity level

圖6 不同迎角升力系數(shù)隨濕度變化曲線Fig.6 Variation of the lift coefficient as a function of humidity level

圖7是相同迎角不同濕度沿翼型表面水蒸汽質(zhì)量分?jǐn)?shù)分布。可以看出，兩種迎角下，30°溫度時(shí)，相對(duì)濕度為50%的翼型表面水蒸汽質(zhì)量分?jǐn)?shù)為定值，說(shuō)明沒(méi)有發(fā)生凝結(jié)；迎角相同，翼型表面水質(zhì)量分布一致；溫度相同時(shí)，相同迎角翼型的表面水蒸汽質(zhì)量分?jǐn)?shù)隨濕度增大，凈質(zhì)量分?jǐn)?shù)增大。濕度相同時(shí)，溫度高凈質(zhì)量分?jǐn)?shù)越大。比較凈質(zhì)量分?jǐn)?shù)數(shù)值，發(fā)現(xiàn)溫度對(duì)凝結(jié)的影響高于濕度對(duì)凝結(jié)的影響。

圖7 不同濕度翼型表面水蒸汽凈質(zhì)量分?jǐn)?shù)分布Fig.7 Water vapor mass fraction distribution across the airfoil surface for different Relative Humidity

圖8是相同濕度和相同溫度下不同迎角翼型表面水蒸氣質(zhì)量分?jǐn)?shù)分布和凈質(zhì)量分?jǐn)?shù)分布圖。凈質(zhì)量分?jǐn)?shù)是凝結(jié)前質(zhì)量分?jǐn)?shù)減去凝結(jié)后的質(zhì)量分?jǐn)?shù)。發(fā)現(xiàn)，兩種溫度下迎角越大，翼型表面水蒸氣的質(zhì)量分?jǐn)?shù)越高。從第三張凈質(zhì)量分?jǐn)?shù)分布圖也看出，迎角越大，凈水蒸氣質(zhì)量分?jǐn)?shù)越大，說(shuō)明凝結(jié)掉的水蒸氣越多。

圖8 不同迎角翼型表面水蒸汽質(zhì)量分?jǐn)?shù)分布Fig.8 Water vapor mass fraction distribution across the airfoil surface for different angles of attack

圖9為不同工況下翼型表面水蒸氣的凝結(jié)質(zhì)量流率。可以看出，溫度相同時(shí)，相對(duì)濕度越高，凝結(jié)質(zhì)量流率越大；相同濕度時(shí)，溫度越高，凝結(jié)質(zhì)量流率越大?？拷熬壧?，凝結(jié)質(zhì)量流率越大。

圖9 不同濕度翼型表面凝結(jié)質(zhì)量流率Fig.9 Condensation mass flux across the airfoil surface for different Relative Humidity

翼型上的諸多力來(lái)源于壓力分布及翼型與流體的表面摩擦，表面摩擦主要形成阻力，對(duì)壓力積分求得的力包括升力和阻力兩個(gè)分量。因此通過(guò)分析流動(dòng)參數(shù)（壓力、速度等）的變化可以了解氣動(dòng)性能發(fā)生變化的機(jī)理。圖10為16°迎角翼型凝結(jié)前后參數(shù)差云圖，此參數(shù)差為90%相對(duì)濕度下，凝結(jié)后各參數(shù)值與凝結(jié)前對(duì)應(yīng)參數(shù)的差值，參數(shù)差能夠反應(yīng)凝結(jié)前后參數(shù)的變化情況。由圖10（a）、（b）、（c）可以看出，開(kāi)啟凝結(jié)模型后，整體的壓力、速度和溫度值增大，這些改變是凝結(jié)導(dǎo)致的，壓力變化主要發(fā)生在翼型上表面前半部分對(duì)應(yīng)的上部區(qū)域，該區(qū)域壓力減小，其他大部分區(qū)域壓力增大；開(kāi)啟凝結(jié)模型后，翼型周?chē)^大部區(qū)域速度值增大，增速最大區(qū)域發(fā)生在流動(dòng)分離區(qū)，漩渦的外側(cè)，速度增大的原因可能是凝結(jié)發(fā)生后，空氣中的水蒸氣析出，空氣的粘性阻力減小，速度增大；相對(duì)應(yīng)的，凝結(jié)放出熱量，使翼型周?chē)黧w溫度升高，且最高溫度差區(qū)對(duì)應(yīng)最大速度差區(qū)，同時(shí)發(fā)現(xiàn)，漩渦中心區(qū)溫度降低。翼型截面相對(duì)濕度差分布與溫度分布差相對(duì)應(yīng)，最高溫度差區(qū)對(duì)應(yīng)于最低相對(duì)濕度差區(qū)，說(shuō)明凝結(jié)發(fā)生后，水蒸氣析出，放出熱量，溫度升高，相對(duì)濕度減小。圖10（f）、（g）、（h）為翼型表面摩擦系數(shù)差、水質(zhì)量分?jǐn)?shù)濕度差及相對(duì)濕度差云圖，可以發(fā)現(xiàn)，翼型表面摩擦系數(shù)差在翼型前緣到上表面約二分之一位置處為正值，說(shuō)明該表面區(qū)域凝結(jié)后，表面摩擦系數(shù)增大，其他翼型表面區(qū)域摩擦系數(shù)變?。灰硇捅砻嫠|(zhì)量分?jǐn)?shù)差分布與相對(duì)濕度差分布近似相同，最大水質(zhì)量分?jǐn)?shù)差位于翼型下表面前半部分，與圖9翼型表面凝結(jié)質(zhì)量流率曲線相對(duì)應(yīng)；水質(zhì)量分?jǐn)?shù)差越大，凝結(jié)質(zhì)量流率越大，凝結(jié)出的水蒸汽越多。

圖10 16°迎角翼型（T＝313K，RH＝90%）凝結(jié)前后參數(shù)差云圖Fig.10 Contours of parameter difference before and after condensation at 16degree angle of attack airfoil（T＝313K，RH＝90%）

圖11為16°迎角翼型90%相對(duì)濕度與50%相對(duì)濕度的參數(shù)變量差，此變量差為開(kāi)啟凝結(jié)模型后，90%相對(duì)濕度下各變量與50%相對(duì)濕度下對(duì)應(yīng)變量的差值，通過(guò)變量差可以了解翼型氣動(dòng)參數(shù)隨濕度變化的趨勢(shì)，有助于流動(dòng)機(jī)理的了解。由圖11可以看出，隨著濕度增大，流過(guò)翼型周?chē)鷫毫υ龃?，?duì)應(yīng)速度增加，相對(duì)濕度增大40%，速度最大增加約0.4m／s，增速較大的區(qū)域在漩渦分離外側(cè)區(qū)；濕度增大，溫度差增大，最大溫度差達(dá)到10.4℃，對(duì)應(yīng)兩種濕度下凝結(jié)露點(diǎn)溫度的差值，且發(fā)生溫度變化的區(qū)域?yàn)橐硇蜕舷卤砻鎸恿鬟吔鐚雍屯牧鬟吔鐚?，觀察圖11（g）相對(duì)濕度差分布圖，與溫度差分布變化趨勢(shì)相同，越靠近翼型表面，濕度差越小，說(shuō)明凝結(jié)量越大，釋放的熱量越多，溫度變化越大。由圖11（e）翼型表面壓力差發(fā)現(xiàn)，翼型表面壓力隨濕度變化的區(qū)域主要在翼型上表面前半部分。由圖11（f）翼型表面摩擦系數(shù)差分布看出，隨濕度增大，表面摩擦系數(shù)增大，且在前緣點(diǎn)附近增幅最大。由圖11（g）翼型表面水質(zhì)量分?jǐn)?shù)差分布可以看出，90%相對(duì)濕度（質(zhì)量分?jǐn)?shù)0.042）較50%相對(duì)濕度（質(zhì)量分?jǐn)?shù)0.023）上表面相差約0.01887，上表面質(zhì)量分?jǐn)?shù)差值較下表面大，說(shuō)明上表面凝結(jié)量少，下表面凝結(jié)量多，這與圖9翼型表面凝結(jié)質(zhì)量流率相對(duì)應(yīng)。

圖11 16°迎角翼型（T＝313K，RH＝90%－50%）不同濕度參數(shù)差云圖Fig.11 Contours of parameter difference of two humidity condition at 16degree angle of attack airfoil（T＝313K，RH＝90%－50%）

4 結(jié) 論

濕空氣中的水汽凝結(jié)影響翼型的氣動(dòng)性能。采用組分輸運(yùn)模型和壁面凝結(jié)模型研究濕空氣流過(guò)翼型表面的氣動(dòng)性能。通過(guò)求解壁面邊界條件凝結(jié)組分的體積源項(xiàng)，實(shí)現(xiàn)凝結(jié)，通過(guò)預(yù)測(cè)平板表面蒸汽凝結(jié)與理論結(jié)果對(duì)比分析，驗(yàn)證了凝結(jié)模型的正確性。對(duì)濕空氣流過(guò)翼型表面氣動(dòng)性能的初步研究發(fā)現(xiàn)，升阻力系數(shù)隨濕度變化而改變，不考慮凝結(jié)時(shí)，升阻力系數(shù)隨濕度增大而減小，當(dāng)發(fā)生凝結(jié)時(shí)，升力系數(shù)增加。迎角不同，翼型表面凝結(jié)分布不同，翼型前緣點(diǎn)附近凝結(jié)量最大；相同迎角時(shí)，凝結(jié)質(zhì)量流率隨濕度升高增大；溫度越高，凝結(jié)量越大。

濕空氣凝結(jié)流動(dòng)涉及復(fù)雜的相變過(guò)程，如果考慮大氣中塵埃、離子等各種微小懸浮顆粒對(duì)凝結(jié)過(guò)程的影響，模型將更為復(fù)雜。因此，有必要在計(jì)算和實(shí)驗(yàn)兩方面進(jìn)一步研究凝結(jié)物理過(guò)程對(duì)翼型乃至葉片氣動(dòng)性能和流動(dòng)的影響。

［1］James F C，Joseph R C，Christopher L R.Oberservation of air plane flow fields by natural condensation effects［J］.Journal of Aircraft，1989，26（7）：593－604.

［2］Schenrr G H，Mundinger G.Similarity，drag and lift in transonic flow with given internal heat addition［J］.European Journal of Mechanics－B／Fluids，1993，12（5）：597－612.

［3］Schenrr G H.Transonic aerodynamics including strong effects from heat addition［J］.Computers Fluids，1993，22（3）：103－116.

［4］Li Liang，Sun Xiuling，F(xiàn)eng Zhenping，et al.Transonic flow of moist air around an NACA 0012airfoil with non－equilibrium condensation［J］.Progress in Natural Science，2005，15（9）：838－842.

［5］Sun Xiuling，Li Liang，Li Guojun.Transonic flow of moist air around the ONERA M6wing with non－equilibrium and homogeneous condensation［J］.Research Journal of Applied Sciences，Engineering and Technology，2013，6（10）：1825－1833.

［6］William L W，Glenn G，Thomas J H.Aircraft－Produced Ice Particles（APIPs）：additional results and further insights［J］.Journal of Applied Meteorology and Climatology，2003，42（5）：640－651.

［7］Jordan F L.Investigation at near－sonic speed of some effects of humidity on the longitudinal aero dynamic characteristics of a NASA supercritical wing research airplane model［R］.Hampton：NASA Langley Research Center，1972.

［8］Schenrr G H，Li P.Shock wave／boundary layer interaction with heat addition by non－equilibrium phase transition［J］.International Journal of Multiphase Flow，1993，19（5）：737－749.

［9］Rusak Z，Lee J C.Transonic flow of moist air around a thin airfoil with nonequilibrium and homogeneous condensation［J］.Journal of Fluid Mechanics，2000，43（1）：173－199.

［10］Karableas S J，Markatos N C.Water vapor condensation in forced convection flow over an airfoil［J］.Aerospace Science and Technology，2008，12（2）：150－158.

［11］Sun Xiuling，Li Liang，Li Guojun.Numerical study of the effects of condensation of moist air on the aerodynamic characteristics of the ONERA M6wing［J］.Acta Aerodynamica Sinica，2010，28（2）：143－148.（in Chinese）孫秀玲，李亮，李國(guó)君.濕空氣凝結(jié)對(duì)ONERA M6機(jī)翼氣動(dòng)特性影響的數(shù)值研究［J］.空氣動(dòng)力學(xué)學(xué)報(bào)，2010，28（2）：143－148.

［12］Sun Xiuling，Li Liang，Li Guojun.Numerical investigation of moist air flow with non－equilibrium condensation around airfoils［J］.Acta Aeronautica et Astronautica Sinica，2010，31（8）：1546－1551.（in Chinese）孫秀玲，李亮，李國(guó)君.翼型濕空氣非平衡凝結(jié)流動(dòng)的數(shù)值研究［J］.航空學(xué)報(bào)，2010，31（8）：1546－1551.

［13］Franck B，Niels S，Jeppe J.Wind turbine airfoil catalogue［R］.Roskilde：RisNational Laboratory，2001.

［14］Sparrow E M，Minkowycz W J，Saddy M.Forced convection condensation in the presence of noncondensables and interfacial resistance［J］.International Journal of Heat and Mass Transfer，1967，10（12）：1829－1845.

［15］Bell B.Application Brief：Film Condensation of Water Vapor［R］.Lebanon，2003.

Study of the influence of moist air on airfoil aerodynamic characteristics

Liu Yan1，＊，Yue Weipeng1，2
（1.North China Electric Power University，Beijing 102206，China；2.Zhangjiakou Wind and Solar Power Energy Demonstration Station CO.LTD State Grid Xin Yuan Company，Zhangjiakou 075000，China）

The purpose of this paper is to study the effect of the moist air over an airfoil.The influence of moist air on the airfoil aerodynamic characteristics is very important for the optimization of offshore wind turbine airfoil and improving the efficiency of wind turbine.Firstly，the lift and drag of airfoil is calculated by using the component transport model.Then，the condensation of water vapor in the moist air is considered，and the condensation volume source is defined by the user defined function in fluent.The preliminary study on the aerodynamic performance of the moist air flow over airfoil shows that the lift and drag coefficient will change with the humidity.The condensation mass flux near the leading edge of the airfoil is the largest，and condensation mass flux of the airfoil increases with the increase of humidity.The condensation mass flux is higher at high temperature.

airfoil；offshore wind turbine；moist airs；aerodynamic characteristics；condensation

V211.3

Adoi：10.7638／kqdlxxb－2015.0205

0258－1825（2016）04－0541－07

2015－12－07；

2016－03－15

劉燕＊（1986－），男，博士研究生，研究方向：海上風(fēng)力機(jī)及風(fēng)場(chǎng)氣動(dòng)性能優(yōu)化.E－mail：sdailiuyan＠163.com

劉燕，岳巍澎.海上濕氣對(duì)翼型氣動(dòng)特性的影響初探［J］.空氣動(dòng)力學(xué)學(xué)報(bào)，2016，34（4）：541－548.

10.7638／kqdlxxb－2015.0205 Liu Y，Yue W P.Study of the influence of moist air on airfoil aerodynamic characteristics［J］.Acta Aerodynamica Sinica，2016，34（4）：541－548.