自然風速下低雷諾數風力機翼型氣動特性隨機數值分析

唐新姿 袁可人 王效禹 彭銳濤

摘? ?要：為量化隨機自然風速條件下風力機翼型氣動特性不確定程度，以S809風力機翼型為研究對象，基于非嵌入式概率配置點法和Transition SST轉捩方程，建立了低雷諾數風力機翼型氣動特性隨機數值分析模型，獲得了自然風速條件下風力機翼型氣動力確定性和不確定性成分比例，并揭示了風速大小和方向隨機耦合作用對翼型流場結構、壓力系數和摩阻分布及湍動能的影響及不確定傳播機制.結果表明，隨機風速風向對翼型升阻氣動因子不確定度影響顯著，在計算攻角范圍內S809翼型升阻比3σ置信區間相對不確定度最大為±35.13%;隨機風速風向耦合作用下翼型升阻比不確定度分別是單隨機因素下的4.76倍和1.08倍;翼型對來流不確定性敏感區域為前緣，可以考慮在翼型前緣部分進行氣動穩健性優化設計.

關鍵詞：風力機翼型;風速;不確定性;隨機數值分析;非嵌入式概率配置點法

Abstract：To quantify the uncertainty of wind turbine airfoil aerodynamic performance under random fluctuation of natural wind， a stochastic numerical analysis model for the aerodynamic characteristics of wind airfoil based on the Transition SST model combining the non-intrusive probabilistic collocation method was established for the S809 airfoil. The determined and uncertain components of aerodynamic characteristics of S809 airfoil under natural wind conditions were obtained. The effects of wind speed and direction uncertainty on the flow field structure， pressure coefficient， friction distribution， and turbulence kinetic energy were quantitated and the uncertainty propagation mechanism was revealed. The results show that the influence of stochastic wind speed and direction on the uncertainty of airfoil lift drag aerodynamic factor is significant. The maximum relative uncertainty of the 3σ confidence interval of lift drag ratio for S809 airfoil within the range of calculated attack angle is ± 35.13%. The uncertainty of lift drag ratio under the coupling effect of stochastic wind speed and direction is 4.76 times and 1.08 times that with single stochastic factor， respectively. The uncertainty sensitive region is observed at the leading edge. Aerodynamic stability optimization design is recommended on the leading edge of the airfoil.

Key words：wind turbine airfoil;wind speed;uncertainty;stochastic numerical analysis;non-intrusive probabilistic collocation method

翼型氣動特性是風力機設計決定性因素，其特性參數通常在確定性條件下獲得.低風型風力機工作雷諾數較低，常年工作在地表層自然來流環境中，氣流極不穩定，風速大小和風向均存在隨機不確定特點，其翼型氣動特性受到不確定風速影響，給風力機實際運行性能與載荷預測帶來困難.研究隨機自然來流對低雷諾數風力機翼型氣動性能的影響，對于風電穩定可靠發展具有重要意義.

對于風力機翼型氣動特性研究，目前國內外主要基于確定性工況條件進行風洞試驗或數值模擬.黃宸武等[1]開展了S809翼型低雷諾數下的氣動特性風洞試驗研究. Wang W C等[2]通過風洞測試研究了有無湍流條件下風力機靜動態性能. 李仁年等[3]采用數值計算方法分析了翼型后緣厚度對翼型升阻力系數的影響. 姜鑫等[4]采用數值模擬計算了多種工況不同曲率葉片翼型升力系數.為考慮非定常流動，朱志斌等[5]采用大渦模擬方法計算低雷諾數下翼型升阻力系數.如果直接在數值模擬中加入不確定邊界條件計算耗時耗力，傳統確定性氣動數值分析已經不能滿足實際工程需求.

為考慮不確定因素，國內外學者將傳統確定性氣動數值分析和不確定理論相結合，發展了不確定數值分析方法解決氣動參數隨機問題.如Abdallah I等[6]將數值模擬與蒙特卡洛法相結合量化翼型升阻力系數不確定性. 該方法依賴于抽樣數量，對于復雜流動分析計算效率不理想，計算時間較長，難以實際應用. 為解決計算效率問題，董世充等[7]、琚亞平等[8]和劉智益等[9]基于譜展開法、嵌入式和非嵌入式混沌多項式方法，分別研究隨機表面粗糙度對風力機翼型氣動特性影響和隨機失諧葉片安裝角對風力機流場及氣動特性影響. 趙軻等[10]基于混沌多項式展開并結合數值模擬方法分析了馬赫數對翼型氣動性能影響. 鄔曉敬等[11]采用非嵌入式混沌多項式方法分析了飛行狀態不確定性對翼型氣動力影響. Zhong X P等[12]與丁繼鋒等[13]采用多項式響應面模型進行不確定馬赫數條件下翼型穩健設計.

風力機隨機氣動問題正受到越來越多研究者的重視，研究發現自然來流風速不確定會導致風力機理論設計與實際運行產生偏差[14]，但由于自然隨機風速很難通過實驗方式產生，風速大小和方向不確定雙因素耦合作用對風力機翼型氣動特性影響及其不確定性在流場中的傳播機制亟待明確.

本文基于多項式混沌法和隨機配點法以及 Transition SST湍流方程，構建隨機翼型氣動數值分析模型，以S809風力機翼型為研究對象，從不確定風速大小和方向均值、方差與氣動力等方面討論隨機自然來流對翼型升力和阻力不確定性影響，量化隨機風速大小和方向工況下風力機翼型氣動特性的不確定程度，為極限載荷和疲勞載荷分布譜分析和可靠性評估提供依據;獲得隨機風速條件下風力機翼型氣動力確定性和不確定性成分比例，通過壓力系數摩阻系數和湍動能分布揭示不確定性在流場中的傳播機制，為風電葉片氣動穩健性研究提供參考.

1? ?研究方法

1.1? ?低雷諾數數值分析方法及驗證

在低雷諾數流動下，準確判斷層流分離與轉捩位置至關重要. 經過湍流模型對比分析，本文采用四方程轉捩模型Transition SST[15]預測層流分離.計算雷諾數為3 × 105，計算模型S809翼型弦長c為1 m，對應自然風速為4.38 m/s，采用二維“C”型結構網格，翼型上下及來流上游方向取15倍翼型弦長，翼型下游方向取20倍弦長作為計算域邊界.采用ICEM生成結構化網格，近壁面進行局部加密，物面法向第一層網格高度為10-5倍弦長，保證壁面處y+<1，網格增長率為1.1，近翼面網格如圖1所示.流域左邊界設為速度入口，右邊界設為壓力出口，壓力設置為絕對壓力101 325 Pa，翼面為無滑移條件.采用二階迎風差分格式進行離散，壓力速度耦合采用SIMPLE算法，收斂精度標準均為殘差小于1 × 10-5.

為驗證網格無關性，計算了風速4.38 m/s、風向6.11°時4種不同網格數量的升阻力系數，計算結果對比如表1所示. 由表1可知，網格數量達到89 700時，增加網格數量計算結果不再明顯變化，綜合考慮計算效率，最終確定網格總數為89 700.

為進一步驗證計算方法的正確性，分別計算不同攻角下S809翼型升阻力系數并與實驗值[16]進行對比，計算結果如圖2所示.結果表明采用上述方法和網格計算所得的升阻力系數計算值與實驗值在低攻角下基本一致，驗證了該方法的可行性.

1.2? ?隨機分析方法與驗證

本文將多項式混沌法和隨機配點法結合，基于Transition SST湍流方程，形成隨機翼型氣動分析方法.多項式混沌法和隨機配點法聯合又稱為非嵌入式概率配置點法（Non-intrusive Probabilistic Collocation Method，NIPRCM），參照文獻[17]，得出任意樣本空間上的變量關系如式（1）所示，其中yk（x，t）是由空間自變量x，時間t和參數θ組成的y在第k個配置點的值，Np為配置點的個數，hk（ξ）為對應的拉格朗日插值多項式混沌.

根據隨機變量分布函數，將選取的高斯積分點及對應的求積系數映射到隨機空間P上，即得到配置點及其權重.得到權重系數后，隨機變量y的均值 y和方差σ2可以通過式（2）和（3）計算，其中yk為物理參數的確定性計算結果，wk為權重系數.

根據文獻[18]可知，用公式（2）和（3）分別采用二階與三階方法計算翼型氣動參數所得計算結果基本一致，說明二階方法計算精度已滿足要求，因此后續采用二階NIPRCM方法進行隨機分析.

為進一步驗證所采用的隨機分析方法，分別計算在風速服從高斯分布的不確定條件下S809翼型升阻力系數并與實驗值[16]進行對比，計算雷諾數為3 ×105，攻角分別為1.99°、4.08°、6.11°、8.14°，標準差為10%均值，計算結果如圖3所示.采用上述隨機方法計算所得的升阻力系數均值與確定性計算值基本一致，且兩者與實驗值吻合，驗證了該方法的可行性.

2? ?隨機風速對翼型氣動特性的影響

為量化風速大小隨機性對風力機翼型氣動性能的影響，計算雷諾數為3 × 105，典型攻角1.99°、6.11°、8.14°時風速大小標準差分別為均值的5%、10%、15%，即湍流強度為5%、10%、15%三種工況下S809翼型升阻氣動特性.

2.1? ?升阻特性

圖4為三種隨機風速條件下攻角為1.99°、6.11°和8.14°時S809翼型的升阻力正態分布，經無關性分析確定隨機計算總頻次為5 000. 標準差越大，即風速隨機性越大，升阻比不確定程度增加，但升阻比均值基本保持不變，且均值對應頻次減少. 攻角8.14°時，隨機性對升阻比影響最大. 由圖4（c）可知，

風速大小標準差σv為15%時，翼型升阻比標準差分別是標準差為5%和10%時升阻比標準差的4.15倍和1.75倍，升阻比3σ不確定度為39.576 7±11.947 5.

2.2? ?壓力分布

圖5給出了攻角1.99°時隨機風速大小條件下S809翼型上翼面壓力系數Cp的標準差分布，圖中橫坐標x/c表示翼型表面點的相對弦長位置，即翼型表面點在弦長方向的投影坐標x與弦長c的比值. 風速標準差越大，壓力系數標準差峰值越大;風速大小標準差由5%增加到15%，壓力系數標準差最大值由0.010 2增加到0.047，且各標準差峰值位置始終出現在翼型0.6倍弦長附近位置，這是由于在此位置翼型邊界層出現層流分離泡[15]，流動分離轉捩，對風速大小波動更為敏感.

2.3? ?摩阻分布

圖6是攻角為1.99°時隨機風速條件下S809翼型上翼面摩阻系數Cf的標準差分布.風速大小標準差由5%增加到15%，摩阻系數標準差最大值由5.58 × 10-4增加到1.66 × 10-3. 翼型前緣出現極大值，且在0.6倍弦長附近位置再次出現峰值，規律與前述壓力分布一致.

2.4? ?湍動能分布

圖7是攻角為1.99°時隨機風速條件下S809翼型上翼面湍動能k的標準差分布.由圖7可知，風速大小標準差由5%增加到15%，湍動能最大標準差由1.21 × 10-4增加到2.91 × 10-4. 湍動能標準差峰值始終位于翼型0.6倍弦長附近位置，與壓力和摩阻系數分布規律一致.

3? ?隨機風向對翼型氣動特性的影響

為分析風速方向隨機性對風力機翼型氣動性能的影響，在雷諾數為3 × 105，攻角1.99°、6.11°、8.14°工況，分別取標準差為均值的5%、10%、15%這三種風向分布，計算風速方向不確定S809翼型氣動性能.

3.1? ?升阻特性

圖8為三種隨機風向條件下攻角為1.99°、6.11°和8.14°時S809翼型升阻比正態分布，計算總頻數為5 000. 風向標準差越大，即風向隨機性越大，升阻比不確定程度增加. 攻角8.14°時，風向隨機性對升阻比影響最大，風向標準差為15%時，翼型升阻比標準差σ 分別是5%和10%時的1.95倍和1.18倍，其升阻比3σ不確定度為39.891 9 ± 24.742 5.

3.2? ?壓力分布

圖9給出了攻角1.99°時隨機風向條件下S809翼型上翼面壓力系數標準差分布. 由圖9可知，風向標準差由5%增加到15%，壓力系數最大標準差由0.032 3增加到0.100 1.壓力系數最大標準差位于前緣，翼型0.6倍弦長附近位置出現小峰值.

3.3? ?摩阻分布

圖10是攻角為1.99°時隨機風向條件下S809翼型上翼面摩阻系數標準差分布. 風向標準差由5%增加到15%，摩阻系數最大標準差由9.29 × 10-4增加到2.92 × 10-3，最大標準差位置位于前緣，且在翼型0.6倍弦長附近位置再次出現峰值.

圖11是攻角為1.99°時隨機風向條件下S809翼型上翼面湍動能標準差分布. 由圖11可知，風向標準差由5%增加到15%，湍動能最大標準差由7.8 × 10-6增加到5.33 × 10-5，其峰值位置位于0.6倍弦長附近位置，敏感位置與壓力摩阻系數分析結論一致.

4? ?隨機風速風向耦合對翼型氣動特性的影響

為分析風速和風向同時存在隨機性對風力機翼型氣動性能的影響，在雷諾數為3 × 105，攻角1.99°、6.11°、8.14°工況，各計算標準差為5%均值，計算S809翼型氣動性能.

4.1? ?升阻特性

圖12給出了風速、風向單隨機因素以及風速和風向同時存在隨機性，攻角分別為1.99°、6.11°和8.14°時S809翼型升阻比正態分布，計算總頻數為5 000. 相同攻角下，隨機耦合作用下的升阻比標準差大于單隨機因素下的升阻比標準差. 攻角8.14°時，風向不確定性對升阻比影響最大. 由圖12（c）可知，風速風向隨機耦合作用下，翼型升阻比3σ不確定度為39.022 3±13.708 5，即3σ置信區間相對不確定度為±35.13%，且其標準差σv，α分別是風速（σv）和風向（σα）單隨機因素時的4.76倍和1.08倍.

4.2? ?壓力分布

圖13為隨機風速、風向以及隨機風速和風向耦合作用下，攻角分別為1.99°、6.11°和8.14°時S809翼型上翼面壓力系數標準差分布. 隨機風速風向耦合作用時，壓力系數不確定程度大于單隨機因素影響程度;且攻角增大，壓力系數不確定度增加.在隨機風速和風向耦合作用下，攻角1.99°、6.11°和8.14°時壓力系數最大標準差分別為0.035 2、0.214 1和0.955 6，各壓力系數標準差最大值均出現在前緣.

5? ?結? ?論

基于NIPRCM方法和Transition SST轉捩模型，建立低雷諾數翼型氣動特性隨機數值分析模型，量化隨機風速大小和方向對風力機翼型升阻力和流場特性的影響，獲得自然來流條件下風力機翼型氣動特性不確定程度及流場非穩傳播規律.主要結論如下：

1）隨機風速風向對翼型升阻氣動因子不確定度影響顯著，且來流隨機性越大升阻比不確定程度越大，在計算攻角范圍內S809翼型升阻比3σ置信區間相對不確定度最大為±35.13%，該值遠大于常規設計標準靜態氣動因子評估不確定度，表明風力機設計應根據實際安裝位置氣流狀況，加強動態時變極限載荷和疲勞載荷風險評估.

2）隨機風速風向耦合作用時翼型升阻力系數不確定度比隨機風向或隨機風速時的不確定度明顯增大，且風向不確定影響更大.攻角8.14°隨機風速風向耦合作用且標準差為5%時，S809翼型升阻比不確定度分別是風速和風向單隨機因素下的4.76倍和1.08倍.

3）在隨機自然風速條件下，翼型對來流不確定性敏感區域為前緣，且中部0.6倍弦長位置轉捩處出現非穩流動加劇，可以考慮在翼型前緣部分進行氣動穩健性優化設計.

參考文獻

[1]? ? 黃宸武，楊科，廖猜猜，等. S809風力機專用翼型低雷諾數氣動特性實驗研究[J]. 工程熱物理學報，2014，35（11）：2197—2201.

HUANG C W，YANG K，LIAO C C，et al. Experimental study on low Reynolds number aerodynamic characteristics of special airfoil of S809 wind turbine[J]. Journal of Engineering Thermophysics，2014，35（11）：2197—2201. （In Chinese）

[2]? ? WANG W C，WANG J J，CHONG W T. The effects of unsteady wind on the performances of a newly developed cross-axis wind turbine：a wind tunnel study[J]. Renewable Energy，2019，131：644—659.

[3]? ? 李仁年，李銀然，王秀勇. 翼型后緣厚度對氣動性能影響的數值分析[J]. 空氣動力學學報，2011，29（2）：205—210.

LI R N，LI Y R，WANG X Y.Numerical analysis of aerodynamic performance of airfoils trailing edge thickness[J]. Acta Aerodynamica Sinica，2011，29（2）：205—210. （In Chinese）

[4]? ? 姜鑫，蘇猛，宋力，等. 柔性葉片風力機翼型氣動性能數值模擬[J]. 工程熱物理學報，2014，35（7）：1317—1320.

JIANG X，SU M，SONG L，et al. Numerical simulation of flexible blade wind turbine airfoil aerodynamic performance[J]. Journal of Engineering Thermophysics，2014，35（7）：1317—1320. （In Chinese）

[5]? ? 朱志斌，尚慶，白鵬，等. 翼型低雷諾數層流分離現象隨雷諾數的演化特征[J]. 航空學報，2019，40（5）：62—72.

ZHU Z B，SHANG Q，BAI P，et al. Evolution of laminar separation phenomenon on low Reynolds number airfoil at different Reynolds numbers[J].Acta Aeronautica et Astronautica Sinica，2019，40（5）：62—72.（In Chinese）

[6]? ? ABDALLAH I，NATARAJAN A，SORENSEN J D. Impact of uncertainty in airfoil characteristics on wind turbine extreme loads[J]. Renewable Energy，2015，75：283—300.

[7]? ? 董世充，王曉東，康順. 表面粗糙度對翼型氣動特性影響的數值模擬[J]. 工程熱物理學報，2015，36（6）：1238—1241.

DONG S C，WANG X D，KANG S. Numerical simulations on the effect of surface roughness on airfoil aerodynamic performance[J]. Journal of Engineering Thermophysics，2015，36（6）：1238—1241.（In Chinese）

[8]? ? 琚亞平，張楚華. 隨機粗糙度對風力機翼型氣動性能影響[J]. 工程熱物理學報，2012，33（5）：784—787.

JU Y P，ZHANG C H.Impact of surface roughness of wind turbine airfoil leading edge on aerodynamic performance[J].Journal of Engineering Thermophysics，2012，33（5）：784—787. （In Chinese）

[9]? ? 劉智益，王曉東，康順. 風輪失諧的氣動力不確定性研究[J]. 工程熱物理學報，2015，36（3）：517—521.

LIU Z Y，WANG X D，KANG S. Non-deterministic research on aerodynamic force of mistuned wind turbine[J].Journal of Engineering Thermophysics，2015，36（3）：517—521. （In Chinese）

[10]? 趙軻，高正紅，黃江濤，等. 基于PCE方法的翼型不確定性分析及穩健設計[J].力學學報，2014，46（1）：10—19.

ZHAO K，GAO Z H，HUANG J T，et al. Uncertainty quantification and robust design of airfoil based on polynomial chaos technique[J].Chinese Journal of Theoretical and Applied Mechanics，2014，46（1）：10—19.（In Chinese）

[11]? 鄔曉敬，張偉偉，宋述芳，等. 翼型跨聲速氣動特性的不確定性及全局靈敏度分析[J]. 力學學報，2015，47（4）：587—595.

WU X J，ZHANG W W，SONG S F，et al. Uncertainty quantification and global sensitivity analysis of transonic aerodynamics about airfoil[J].Chinese Journal of Theoretical and Applied Mechanics，2015，47（4）：587—595. （In Chinese）

[12]? ZHONG X P，DING J F，LI W J，et al. Robust airfoil optimization with multi-objective estimation of distribution algorithm[J]. Chinese Journal of Aeronautics，2008，21（4）：289—295.

[13]? 丁繼鋒，李為吉，張勇，等. 基于響應面的翼型穩健設計研究[J].空氣動力學學報，2007，25（1）：19—22.

DING J F，LI W J，ZHANG Y，et al. Robust airfoil optimization based on response surface method[J].Acta Aerodynamica Sinica，2007，25（1）：19—22. （In Chinese）

[14]? ANTONINI E G A，ROMERO D A，AMON C H.Improving CFD wind farm simulations incorporating wind direction uncertainty[J]. Renewable Energy，2019，133：1011—1023.

[15]? 唐新姿，黃軒晴，孫松峰，等. 考慮層流分離的低速風力機翼型氣動性能研究[J]. 動力工程學報，2017，37（1）：52—59.

TANG X Z，HUANG X Q，SUN S F，et al. Study on aerodynamic characteristics of low-speed wind turbine airfoil in consideration of laminar separation[J]. Journal of Chinese Society of Power Engineering，2017，37（1）：52—59. （In Chinese）

[16]? BUTTERFIELD C P，MUSIAL W P，SIMMS D A. Combined experiment phase I final report[R]. Golden，Colorado：National Renewable Energy Laboratory，1992：1—74.

[17]? LOEVEN G J A，WITTEVEEN J A S，BIJL H. Probabilistic collocation：an efficient non-intrusive approach for arbitrarily distributed parametric uncertainties[C]//45th AIAA Aerospace Sciences Meeting and Exhibit. Reston，Virginia：AIAA，2007.

[18]? 唐新姿，李鵬程，彭銳濤，等. 湍流工況小型風力機翼型氣動特性及穩健優化[J]. 機械工程學報，2020，56（2）：192—200.

TANG X Z，LI P C，PENG R T，et al. Aerodynamic characteristics and robust optimization of small wind turbine airfoil under turbulence condition[J]. Journal of Mechanical Engineering，2020，56（2）：192—200. （In Chinese）