黃國慶 劉廷林 閆渤文 莫蕊瑜 鄧曉蔚 信志強(qiáng)

DOI：?10.11835/j.issn.2096-6717.2022.019

基金項(xiàng)目：國家自然科學(xué)基金（51878104）；外國專家重點(diǎn)支撐計(jì)劃（zcsf2021008）

作者簡介：黃國慶（1976-?），男，博士，教授，主要從事高層建筑、風(fēng)機(jī)和輸電線塔等結(jié)構(gòu)抗風(fēng)以及隨機(jī)振動(dòng)研究，E-mail：ghuang1001@gmail.com。

通信作者：閆渤文（通信作者），男，博士，副教授，E-mail：bowenyancq@cqu.edu.cn。

Received： 2021?11?30

Foundation items： National Natural Science Foundation of China （No. 51878104）; Key Support Program of China for Foreign Experts （No. zcsf2021008）

Author brief： HUANG Guoqing （1976-?）， PhD， professor， main research interests： wind resistance of tall building， wind turbine and transmission tower; random vibration， E-mail： ghuang1001@gmail.com.

corresponding author：YAN Bowen （corresponding author），?PhD，?associate professor，?E-mail：?bowenyancq@cqu.edu.cn.

（1. 重慶大學(xué)?土木工程學(xué)院，重慶?400045；?2. 中國船舶集團(tuán)有限公司中船海裝風(fēng)電股份有限公司，重慶?401122；?3. 香港大學(xué)?土木工程系，香港?999077；?4. 河海大學(xué)?力學(xué)與材料學(xué)院，南京?211100）

摘要：隨著大型風(fēng)電基地建設(shè)，上游風(fēng)機(jī)在運(yùn)行時(shí)會(huì)使下游風(fēng)場風(fēng)速下降，湍流度增大，造成下游風(fēng)機(jī)發(fā)電功率降低，加劇風(fēng)機(jī)的疲勞破壞并縮短其服役周期。因此，亟需開展風(fēng)機(jī)尾流研究，明確其特性及演化規(guī)律。為了揭示不同入流及偏航角下的單風(fēng)機(jī)尾流特性，基于單風(fēng)機(jī)尾流風(fēng)洞試驗(yàn)，驗(yàn)證基于大渦模擬（Large Eddy Simulation， LES）結(jié)合致動(dòng)線模型（Actuator Line Model， ALM）數(shù)值模擬方法的準(zhǔn)確性；基于LES-ALM模擬方法研究入流風(fēng)場（包括風(fēng)速及湍流度）及偏航角對風(fēng)機(jī)尾流特性的影響，闡明正負(fù)偏航角下單風(fēng)機(jī)尾流的對稱性。結(jié)果表明：隨著背景湍流度的增大，風(fēng)機(jī)尾流恢復(fù)速度加快；當(dāng)入流條件相同時(shí)，風(fēng)機(jī)設(shè)置正負(fù)對稱偏航角，其尾流風(fēng)速也表現(xiàn)出一定的對稱性；風(fēng)機(jī)偏航角越大，風(fēng)機(jī)尾流膨脹寬度會(huì)逐漸減小，并降低尾流風(fēng)速的虧損程度。

關(guān)鍵詞：單風(fēng)機(jī)尾流；偏航角；大渦模擬-致動(dòng)線模型；來流條件

中圖分類號：TK89 ????文獻(xiàn)標(biāo)志碼：A ????文章編號：2096-6717（2024）03-0115-09

Wake characteristics of single wind turbine under different inflow conditions and yaw angles

HUANG Guoqing¹，?LIU Tinglin¹，?YAN Bowen¹，?MO Ruiyu²，?DENG Xiaowei³，?XIN Zhiqiang⁴

（1. School of Civil Engineering， Chongqing University， Chongqing 400045， P. R. China；?2. CSSC Haizhuang Wind Power Co.， Ltd.， Chongqing 401122， P. R. China；?3. Department of Civil Engineering， The University of Hong Kong， Hong Kong 999077， P. R. China；?4. College of Mechanics and Materials， Hohai University， Nanjing 211100， P. R. China）

Abstract： With the rapid construction of large-scale wind farms， the upstream wind turbine will deteriorate the wind velocity in the wind turbine wake region and increase the wake turbulence， thereby reducing the power efficiency of the downstream wind turbine and aggravating the fatigue damage of the wind turbine to shorten its service life. Therefore， it is urgent to carry out the study of wind turbine wake characteristics. The main purpose of this study is to reveal the wake characteristics of a standalone wind turbine under different upstream flow conditions and with various yaw angles. Firstly， the accuracy of Large Eddy Simulation （LES） combined with Actuator Line Model （ALM） was verified by the corresponding wind tunnel test. And then， based on the LES-ALM simulation method， the effects of inlet （including velocity and turbulence intensity） and yaw angle on the turbine wake characteristics were studied， and the symmetry of the horizontal profiles of the wind turbine wake at the hub height was clarified with positive and negative yaw angles. The results show that with the increase of the background turbulence intensity， the recovery of the wind turbine wake will be accelerated. Moreover， under the same upstream flow conditions and symmetrical yaw angle， the certain symmetry of wind turbine wake velocity was also identified. When the yaw angle of the wind turbine increases， both the expansion rate of the turbine wake width and the deficit degree of the wake wind velocity would decrease.

Keywords： wind turbine wake；?yaw angle；?LES-ALM；?upstream flow conditions

傳統(tǒng)化石能源環(huán)境污染較大，且逐漸枯竭，亟須發(fā)展環(huán)境友好的清潔可再生能源。風(fēng)能由于其商業(yè)化程度高，是全世界可再生能源開發(fā)的重點(diǎn)。在大型風(fēng)電場中，上游風(fēng)機(jī)在運(yùn)行中會(huì)造成下游流場出現(xiàn)風(fēng)速降低且湍流增加的尾流區(qū)域^[1]。據(jù)統(tǒng)計(jì)，海上風(fēng)電場中風(fēng)機(jī)尾流引起的發(fā)電量損失達(dá)到10%～20%^[2]。因此，對風(fēng)機(jī)尾流特性的研究成為近年來風(fēng)能領(lǐng)域的熱點(diǎn)之一。

近年來，由于計(jì)算性能的大幅提升和計(jì)算流體力學(xué)（Computational Fluid Dynamics， CFD）方法的快速發(fā)展，基于數(shù)值仿真方法對風(fēng)機(jī)尾流進(jìn)行的研究越來越多^[3-5]。Xie等^[6]使用大渦模擬的方法研究了單臺(tái)風(fēng)機(jī)的尾流特性，他們發(fā)現(xiàn)尾流風(fēng)速虧損在水平方向和垂直方向的上部保持自相似的高斯特性，在接近地面的尾流區(qū)域內(nèi)無該特性。Vermeer等^[7]研究了均勻來流下的風(fēng)機(jī)近尾流特性，并研究了風(fēng)機(jī)尾流在遠(yuǎn)尾流區(qū)域的衰減規(guī)律。Churchfield等^[8]揭示了大氣穩(wěn)定性和地面粗糙度對風(fēng)機(jī)響應(yīng)、發(fā)電量以及尾流發(fā)展有顯著影響。Na等^[9]的研究表明，風(fēng)機(jī)尾流干擾對下游風(fēng)機(jī)的功率輸出及風(fēng)機(jī)載荷有較大影響。但目前系統(tǒng)開展不同來流風(fēng)速及湍流度對尾流發(fā)展演化規(guī)律影響的研究相對較少。

風(fēng)電場中通常使用尾流主動(dòng)控制的方法減少尾流之間的相互影響。尾流主動(dòng)控制的核心是通過上游風(fēng)機(jī)的偏航以減少下游風(fēng)機(jī)所在區(qū)域的風(fēng)速虧損^[10]。風(fēng)機(jī)的偏航角越大，其尾流向側(cè)邊偏移的程度也就越大^[11]。仿真研究表明，通過優(yōu)化風(fēng)電場的偏航角，兩串列式風(fēng)機(jī)的發(fā)電量至少可以提升5%^[12]。Adaramola等^[13]的研究也表明偏航控制在提高風(fēng)電場的發(fā)電量上有巨大潛力。因此，研究不同偏航角對風(fēng)機(jī)尾流發(fā)展的影響具有重要意義。

Dou等^[14]通過風(fēng)洞試驗(yàn)將偏航后的風(fēng)機(jī)輪轂高度平面的遠(yuǎn)尾流（風(fēng)機(jī)下游5倍風(fēng)輪直徑及以上位置處尾流）與基于高斯分布的尾流模型進(jìn)行對比，發(fā)現(xiàn)兩者符合程度較好。該結(jié)論可驗(yàn)證偏航后風(fēng)機(jī)輪轂高度平面的遠(yuǎn)尾流同樣符合高斯分布，證明偏航后的遠(yuǎn)尾流分布在統(tǒng)計(jì)上具有一定的對稱性。而在實(shí)際工程中，下游風(fēng)機(jī)一般處于上游風(fēng)機(jī)的遠(yuǎn)尾流影響之中，且風(fēng)機(jī)偏航包括正向偏航（風(fēng)輪面逆時(shí)針旋轉(zhuǎn)）和負(fù)向偏航（風(fēng)輪面順時(shí)針旋轉(zhuǎn)），風(fēng)機(jī)尾流也具有旋轉(zhuǎn)性。在偏航與尾流旋轉(zhuǎn)兩個(gè)因素耦合作用下，正負(fù)對稱偏航角下的遠(yuǎn)尾流是否具有對稱性值得進(jìn)一步研究。

筆者主要研究風(fēng)機(jī)尾流在不同入流、不同湍流度和不同偏航角的情況下的發(fā)展規(guī)律，并探究風(fēng)機(jī)在正負(fù)對稱偏航時(shí)遠(yuǎn)尾流是否具有對稱性。

1 數(shù)值方法

1.1 大渦模擬

大渦模擬根據(jù)不同渦旋的尺度對湍流中的渦分別進(jìn)行處理，利用瞬時(shí)的N-S方程直接求解大尺度渦旋，而小尺度渦旋對大尺度渦旋的影響則是通過亞格子尺度模型來近似。結(jié)合實(shí)際情況，在不可壓縮N-S方程中加入科里奧利力項(xiàng)以及浮力項(xiàng)來模擬地球自轉(zhuǎn)效應(yīng)以及由地表溫度差異引起的垂向浮力。聯(lián)立連續(xù)性方程得到的控制方程為 （2）

式中：Ⅰ為由地球自轉(zhuǎn)引起的科里奧利力項(xiàng)；Ⅱ?yàn)閴毫π拚?xiàng)；Ⅲ為背景壓力梯度項(xiàng)，作為推動(dòng)風(fēng)形成的驅(qū)動(dòng)力；Ⅳ為應(yīng)力張量項(xiàng)，其包括黏性應(yīng)力和?；瘉喐褡討?yīng)力，本研究采用Smagorinsky亞格子模型^[15]，其表達(dá)式為 （3） （4）（5）

式中：為可解尺度湍流變形張量；為亞格子渦粘模型；為濾波尺寸；Smagorinsky常數(shù)取值0.14。式（2）中Ⅴ為浮力項(xiàng)；Ⅵ為風(fēng)機(jī)葉片對流場的作用力源項(xiàng)，該項(xiàng)來自致動(dòng)線模型。

1.2 致動(dòng)線模型

數(shù)值仿真中，使用致動(dòng)線模型代替風(fēng)機(jī)。致動(dòng)線模型最早由So?rensen等^[16]提出。該模型思想在于將真實(shí)風(fēng)機(jī)葉片沿徑向離散，每個(gè)離散點(diǎn)用虛擬體積力代替。該方法不需要建立風(fēng)機(jī)葉片實(shí)體模型，無需求解葉片表面邊界層流動(dòng)，大大降低了建模難度，提高了計(jì)算效率。致動(dòng)線模型在模擬葉片旋轉(zhuǎn)及捕捉葉尖與葉根的渦脫落具有較高精度。相關(guān)研究表明，致動(dòng)線模型是除實(shí)體建模以外的一種精度較高的風(fēng)機(jī)流場計(jì)算方法^[17]。對于致動(dòng)線模型中的每個(gè)致動(dòng)點(diǎn)上的體積力，可以用葉素動(dòng)量理論計(jì)算，其表達(dá)式為 （6） （7）

式中：為相對于葉素的空氣流速；與分別為攻角對應(yīng)的升力與阻力系數(shù)；為攻角；為半徑位置處葉素弦長，為葉素寬度；為空氣密度，取值1.225 kg/m³。 （8）

式中：為葉片的旋轉(zhuǎn)角速度；和為軸向速度和切向速度。

2 數(shù)值方法驗(yàn)證

使用LES-ALM數(shù)值計(jì)算方法，為驗(yàn)證其計(jì)算準(zhǔn)確性，將該方法計(jì)算結(jié)果與挪威科技大學(xué)（Norges Teknisk-Naturvitenskapelige Universitet， NTNU）研究團(tuán)隊(duì)開展的單風(fēng)機(jī)尾流風(fēng)洞試驗(yàn)的風(fēng)速虧損結(jié)果^[18]進(jìn)行對比。

2.1 風(fēng)洞試驗(yàn)與數(shù)值仿真參數(shù)設(shè)置

風(fēng)洞試驗(yàn)^[18]中的風(fēng)洞長×寬×高為11.15 m×2.72 m×1.8 m。風(fēng)機(jī)模型的翼型型號為S826，直徑D為0.894 m，輪轂高度H_hub為0.817 m。風(fēng)機(jī)模型的弦長與扭轉(zhuǎn)角沿徑向分布圖以及風(fēng)機(jī)的升阻力系數(shù)隨攻角的變化如圖1所示。試驗(yàn)中來流為均勻定常流，輪轂中心的參考風(fēng)速=10 m/s，風(fēng)機(jī)葉尖速比為6。對應(yīng)在數(shù)值仿真中，設(shè)置長×寬×高為10 m×4.6 m×2 m的計(jì)算域，網(wǎng)格使用blockmesh生成，通過網(wǎng)格無關(guān)性驗(yàn)證，確保網(wǎng)格精度達(dá)到要求，計(jì)算域共960萬網(wǎng)格。入流參數(shù)和風(fēng)機(jī)參數(shù)與風(fēng)洞試驗(yàn)保持一致，數(shù)值仿真中使用的邊界條件如表1所示。

2.2 風(fēng)洞試驗(yàn)與數(shù)值仿真對比驗(yàn)證

圖2為數(shù)值仿真中計(jì)算的風(fēng)速虧損與風(fēng)洞試驗(yàn)測量的風(fēng)速虧損值對比，從圖中可以明顯看出，在x/D=1、3位置處，數(shù)值仿真和風(fēng)洞試驗(yàn)結(jié)果整體趨勢符合較好，但是在輪轂中心位置處，數(shù)值仿真與風(fēng)洞試驗(yàn)的結(jié)果有一定偏差，這是由致動(dòng)線模型忽略輪轂對流場的影響造成的。而x/D=5位置距離風(fēng)機(jī)較遠(yuǎn)，輪轂對該位置已無影響，從圖中也可看出，無論是數(shù)值大小還是尾流趨勢，數(shù)值仿真與風(fēng)洞試驗(yàn)都符合較好。因此可認(rèn)為LES-ALM方法能較好地計(jì)算風(fēng)機(jī)尾流流場特征。

3 數(shù)值結(jié)果分析

3.1 入流生成與數(shù)值仿真基本參數(shù)

控制好入流特性是研究不同入流參數(shù)下風(fēng)機(jī)尾流發(fā)展規(guī)律的前提。使用能較好控制入流特性的CDRFG（Consistent Discrete Random Flow Generation）方法^[19]來生成入流，該方法是一種大渦模擬入口湍流生成方法^[18]。研究中考慮由3種風(fēng)速（U_h=6、10.4、15 m/s）和3種湍流度（I=2%、6%、10%）組合成的9種入流。

選用風(fēng)輪直徑D為152 m的5 MW風(fēng)機(jī)，其基本參數(shù)如表2所示。上述入流生成方法雖然能較好地模擬出目標(biāo)的脈動(dòng)風(fēng)速譜，并合理保持速度之間的相干性，但也具有一定的局限性。該方法生成的湍流風(fēng)場能滿足連續(xù)性方程，卻較難滿足動(dòng)量方程。為使入流在計(jì)算域內(nèi)能較好保持，避免出現(xiàn)湍流衰減過快現(xiàn)象，本研究將原型風(fēng)機(jī)進(jìn)行1：250縮尺，風(fēng)機(jī)縮尺后參數(shù)如表2所示。足尺風(fēng)機(jī)與同型號縮尺風(fēng)機(jī)的尾流特性主要受雷諾數(shù)影響，研究中風(fēng)機(jī)縮尺后的雷諾數(shù)為2×10⁴，Chamorro等^[20]的研究表明，該量級及以上的雷諾數(shù)對風(fēng)機(jī)的遠(yuǎn)尾流特性的影響不大。因此，雷諾數(shù)對數(shù)值仿真結(jié)果影響較小，可利用縮尺后的風(fēng)機(jī)進(jìn)行尾流特性分析。

數(shù)值仿真中計(jì)算域的長×寬×高為8 m×2.4 m×1.8 m，其網(wǎng)格劃分如圖3所示，網(wǎng)格最大尺寸為0.05 m，最小尺寸為0.012 m。且該網(wǎng)格劃分已通過無關(guān)性驗(yàn)證，最終網(wǎng)格數(shù)目達(dá)到587萬。計(jì)算時(shí)間步長取入流風(fēng)速跑過整個(gè)計(jì)算域所需時(shí)間的1/1 500。計(jì)算的所有工況如表3所示。

3.2 入流自保持性驗(yàn)證結(jié)果分析

將生成的入流在空風(fēng)場（計(jì)算域中不加入致動(dòng)線模型）中計(jì)算的風(fēng)速和湍流度結(jié)果與目標(biāo)設(shè)定值進(jìn)行對比，驗(yàn)證入流的自保持性。圖4、圖5分別為不同工況輪轂中心位置處的平均風(fēng)速和湍流度沿順風(fēng)向的變化圖。從圖4可知，風(fēng)速的保持性較好，輪轂中心處的風(fēng)速值基本穩(wěn)定在目標(biāo)值附近。從圖5可知，2%和10%湍流度對應(yīng)工況的湍流自保持性整體也較好，6%湍流度工況與目標(biāo)湍流度的誤差在±1%之內(nèi)。因此，認(rèn)為入流在計(jì)算域中的自保持性較好，可進(jìn)行后續(xù)尾流流場計(jì)算。

3.3 不同入流下風(fēng)機(jī)尾流發(fā)展規(guī)律

在實(shí)際工程中，風(fēng)機(jī)在無偏航狀態(tài)下，尾流的發(fā)展主要受入流的影響。圖6是入流風(fēng)速和入流湍流度為6 m/s和2%、6 m/s和6%、6 m/s和10%、10.4 m/s和6%、15 m/s和6%下的單風(fēng)機(jī)尾流輪轂高度平面風(fēng)速歸一化后的云圖，云圖所描述的物理量為平均風(fēng)速歸一化（）后的結(jié)果。圖7為對應(yīng)工況的不同位置的風(fēng)速歸一化后的廓線。其中圖6（a）、（b）、（c）與圖7（a）、（b）、（c）為對比不同湍流度對風(fēng)機(jī)尾流的影響，圖6（b）、（d）、（e）為對比不同風(fēng)速對風(fēng)機(jī)尾流的影響。

從圖6可以看出，在風(fēng)機(jī)后都存在一個(gè)風(fēng)速減小的區(qū)域，并且隨著離風(fēng)機(jī)的位置越來越遠(yuǎn)，尾流風(fēng)速與來流風(fēng)速差距越小。而在圖7中可以更加清晰地看出，距離風(fēng)機(jī)越遠(yuǎn)，尾流的影響半徑逐漸增大。這就是尾流風(fēng)速的兩個(gè)重要特征：隨著與風(fēng)機(jī)距離的增加，尾流風(fēng)速會(huì)逐漸恢復(fù)為來流風(fēng)速，并且尾流風(fēng)速的膨脹半徑會(huì)逐漸增加。從圖6（a）、（b）、（c）和圖7（a）、（b）、（c）可看出，隨著入流湍流度的增加，尾流風(fēng)速的虧損逐漸減小，尾流風(fēng)速的恢復(fù)速度逐漸加快。從圖6（b）、（d）、（e）和圖7（b）、（d）、（e）可看出，在相同湍流度下，當(dāng)入流風(fēng)速小于額定風(fēng)速時(shí)，尾流對風(fēng)速變化不敏感。當(dāng)風(fēng)速大于額定風(fēng)速時(shí)，即在風(fēng)速為15 m/s時(shí)，風(fēng)機(jī)后x/D=4位置處尾流的最低速度為來流風(fēng)速的87.1%，并且不同位置的風(fēng)速相差較小。因此，可認(rèn)為在風(fēng)速大于額定風(fēng)速情況下，尾流風(fēng)速恢復(fù)快。

3.4 正負(fù)對稱偏航下尾流對稱性分析

尾流的發(fā)展在風(fēng)機(jī)偏航與尾流旋轉(zhuǎn)兩個(gè)因素的耦合作用下是否依然具有對稱性，該問題值得研究。若對稱性較好，則后續(xù)研究中只需研究正或負(fù)偏航工況或者即可。因LES-ALM數(shù)值方法沒有考慮機(jī)艙與塔筒對風(fēng)機(jī)尾流的影響，所以主要研究遠(yuǎn)尾流處輪轂高度平面尾流風(fēng)速對稱性。

以入流風(fēng)速6 m/s和入流湍流度6%為例，研究單風(fēng)機(jī)偏航角為±10°、±20°、±30°時(shí)對輪轂高度平面遠(yuǎn)尾流的對稱性影響。圖8為輪轂高度平面處經(jīng)歸一化處理之后的平均風(fēng)速云圖。從圖8可以定性分析，風(fēng)機(jī)在正負(fù)對稱偏航情況下，其尾流也表現(xiàn)出一定的對稱性。為更好地比較，圖9給出不同對稱偏航角下，風(fēng)機(jī)后x/D=3～10位置處輪轂高度平面風(fēng)速歸一化后的廓線圖，其橫坐標(biāo)的相鄰刻度間距表示范圍為0～1。圖中對負(fù)偏航工況的風(fēng)速廓線進(jìn)行對稱處理。負(fù)偏航風(fēng)速廓線對稱處理后的結(jié)果與正偏航的風(fēng)速廓線相比，其尾流風(fēng)速的數(shù)值大小及發(fā)展趨勢都符合得很好。

上述現(xiàn)象在入流風(fēng)速和入流湍流度分別為6 m/s和10%、10.4 m/s和2%以及10.4 m/s和6%條件下也較為明顯。因此可得到結(jié)論，當(dāng)來流情況相同時(shí)，單風(fēng)機(jī)在正負(fù)對稱偏航角下，其輪轂高度平面的遠(yuǎn)尾流處風(fēng)速也表現(xiàn)出較好的對稱性。

3.5 不同偏航角下風(fēng)機(jī)尾流發(fā)展規(guī)律

偏航控制是工程中常用的減少上游風(fēng)機(jī)尾流對下游風(fēng)機(jī)影響的方法。單風(fēng)機(jī)正負(fù)對稱偏航下，其尾流發(fā)展也呈現(xiàn)較好的對稱性，因此主要對比偏航角為10°、20°和30°與無偏航工況的尾流結(jié)果。圖10（a）、（b）、（c）分別是入流風(fēng)速和入流湍流度為6 m/s和6%、10.4 m/s和6%、15 m/s和6%條件下風(fēng)機(jī)后6D位置處有無偏航的尾流風(fēng)速廓線對比圖。圖11為對應(yīng)工況x/D=6位置處的不同風(fēng)速和不同偏航角下尾流風(fēng)速云圖。從圖10、圖11中可看出，單風(fēng)機(jī)采用偏航控制后，尾流風(fēng)速會(huì)向一側(cè)偏移，尾流風(fēng)速虧損的最大值減小，并且隨著偏航角的增大，上述效果也越來越明顯。當(dāng)風(fēng)速大于額定風(fēng)速時(shí)，偏航對其尾流的影響會(huì)越來越小，尾流的偏轉(zhuǎn)效應(yīng)較為不明顯。因此可得到結(jié)論：偏航會(huì)使尾流發(fā)生偏轉(zhuǎn)，偏航角越大，尾流偏轉(zhuǎn)程度越大。并且偏航后會(huì)使得尾流虧損幅值減小。隨著風(fēng)速增大，尾流的偏航效應(yīng)會(huì)減弱。

4 結(jié)論

基于縮尺風(fēng)機(jī)模型，使用LES-ALM方法研究不同入流風(fēng)速、入流湍流度、不同偏航角條件下的單風(fēng)機(jī)尾流特性規(guī)律?？傻玫揭韵陆Y(jié)論：

1）來流為均勻定常流時(shí)，將LES-ALM方法計(jì)算尾流流場的結(jié)果與NTNU開展的風(fēng)洞試驗(yàn)對比，符合情況較好，說明該方法能較好模擬均勻流下的流場特征。

2）利用CDRFG方法生成的入流風(fēng)速及湍流度在基于LES-ALM方法的數(shù)值計(jì)算中具有較好的自保持性，從而為縮尺后風(fēng)機(jī)的尾流數(shù)值仿真提供了支撐。

3）隨著尾流的發(fā)展，距離風(fēng)機(jī)越遠(yuǎn)處尾流的虧損越小，尾流風(fēng)速的膨脹半徑越大。當(dāng)來流風(fēng)速小于額定風(fēng)速時(shí)，尾流發(fā)展對于風(fēng)速較為不敏感，當(dāng)來流風(fēng)速大于額定風(fēng)速時(shí)，風(fēng)機(jī)的尾流效應(yīng)相對較弱，尾流風(fēng)速虧損較小。并且隨著湍流度的增大，尾流的恢復(fù)速度越快，對應(yīng)相同位置的速度虧損越小。

4）在同一入流下，單風(fēng)機(jī)設(shè)置正負(fù)對稱偏航后，其輪轂高度平面的遠(yuǎn)尾流處風(fēng)速發(fā)展也呈現(xiàn)出較好的對稱性，為簡化計(jì)算工作量，后續(xù)在進(jìn)行單風(fēng)機(jī)尾流研究時(shí)，可只研究正偏航工況或者只研究負(fù)偏航工況即可。

5）當(dāng)單風(fēng)機(jī)設(shè)置主動(dòng)偏航控制時(shí)，尾流會(huì)發(fā)生一定的偏轉(zhuǎn)，尾流膨脹寬度會(huì)相對減小，且尾流虧損幅值低于同一入流下的無偏航工況。當(dāng)風(fēng)速小于額定風(fēng)速時(shí)，偏航角越大，上述效應(yīng)越明顯。當(dāng)風(fēng)速大于額定風(fēng)速時(shí)，由于尾流虧損本身較小，偏航后的尾流偏轉(zhuǎn)效應(yīng)較為不明顯。

參考文獻(xiàn)

[1] ?ARCHER C L， VASEL-BE-HAGH A， YAN C， et al. Review and evaluation of wake loss models for wind energy applications [J]. Applied Energy， 2018， 226： 1187-1207.

[2] ?BARTHELMIE R J， HANSEN K， FRANDSEN S T， et al. Modelling and measuring flow and wind turbine wakes in large wind farms offshore [J]. Wind Energy， 2009， 12（5）： 431-444.

[3] ?IVANELL S， S?RENSEN J N， MIKKELSEN R， et al. Analysis of numerically generated wake structures [J]. Wind Energy， 2009， 12（1）： 63-80.

[4] ?TROLDBORG N， S?RENSEN J N， MIKKELSEN R. Actuator line simulation of wake of wind turbine operating in turbulent inflow [J]. Journal of Physics： Conference Series， 2007， 75： 012063.

[5] ?WU Y T， PORTE-AGEL F. Large-eddy simulation of wind-turbin： Evaluation of turbine parametrisations [J]. Boundary-Layer Meteorology， 2011， 138： 345-366.

[6] ?XIE S B， ARCHER C. Self-similarity and turbulence characteristics of wind turbine wakes via large-eddy simulation [J]. Wind Energy， 2015， 18（10）： 1815-1838.

[7] ?VERMEER L J， S?RENSEN J N， CRESPO A. Wind turbine wake aerodynamics [J]. Progress in Aerospace Sciences， 2003， 39（6/7）： 467-510.

[8] ?CHURCHFIELD M J， LEE S， MICHALAKES J， et al. A numerical study of the effects of atmospheric and wake turbulence on wind turbine dynamics [J]. Journal of Turbulence， 2012， 13： N14.

[9] ?NA J S， KOO E， KO S C， et al. Stochastic characteristics for the vortical structure of a 5-MW wind turbine wake [J]. Renewable Energy， 2019， 133： 1220-1230.

[10] ?ARCHER C L， VASEL-BE-HAGH A. Wake steering via yaw control in multi-turbine wind farms： Recommendations based on large-eddy simulation [J]. Sustainable Energy Technologies and Assessments， 2019， 33： 34-43.

[11] ?BROMM M， ROTT A， BECK H， et al. Field investigation on the influence of yaw misalignment on the propagation of wind turbine wakes [J]. Wind Energy， 2018， 21（11）： 1011-1028.

[12] ?AHMAD T， BASIT A， AHSAN M， et al. Implementation and analyses of yaw based coordinated control of wind farms [J]. Energies， 2019， 12（7）： 1266.

[13] ?ADARAMOLA M S， KROGSTAD P ?. Experimental investigation of wake effects on wind turbine performance [J]. Renewable Energy， 2011， 36（8）： 2078-2086.

[14] ?DOU B Z， QU T M， LEI L P， et al. Optimization of wind turbine yaw angles in a wind farm using a three-dimensional yawed wake model [J]. Energy， 2020， 209： 118415.

[15] ?GERMANO M， PIOMELLI U， MOIN P， et al. A dynamic subgrid-scale eddy viscosity model [J]. Physics of Fluids A： Fluid Dynamics， 1991， 3（7）： 1760-1765.

[16] ?SO?RENSEN J N， SHEN W Z. Numerical modeling of wind turbine wakes [J]. Journal of Fluids Engineering， 2002， 124（2）： 393-399.

[17] ?VRANA T K， FLYNN D， GOMEZ-LAZARO E， et al. Wind power within European grid codes： Evolution， status and outlook [J]. WIREs Energy and Environment， 2018， 7（3）： e285.

[18] ?KROGSTAD P ?， ERIKSEN P E. “Blind test”?calculations of the performance and wake development for a model wind turbine [J]. Renewable Energy， 2013， 50： 325-333.

[19] ?ABOSHOSHA H， ELSHAER A， BITSUAMLAK G T， et al. Consistent inflow turbulence generator for LES evaluation of wind-induced responses for tall buildings [J]. Journal of Wind Engineering and Industrial Aerodynamics， 2015， 142： 198-216.

[20] ?CHAMORRO L P， ARNDT R E A， SOTIROPOULOS F. Reynolds number dependence of turbulence statistics in the wake of wind turbines [J]. Wind Energy， 2012， 15（5）： 733-742.

（編輯??王秀玲）