CN112347611A - 一种风力机远场尾流流向湍流度计算方法 - Google Patents
一种风力机远场尾流流向湍流度计算方法 Download PDFInfo
- Publication number
- CN112347611A CN112347611A CN202011102957.9A CN202011102957A CN112347611A CN 112347611 A CN112347611 A CN 112347611A CN 202011102957 A CN202011102957 A CN 202011102957A CN 112347611 A CN112347611 A CN 112347611A
- Authority
- CN
- China
- Prior art keywords
- wind turbine
- turbulence
- flow direction
- wake
- turbine generator
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000012937 correction Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 6
- 238000011156 evaluation Methods 0.000 claims description 2
- 239000004744 fabric Substances 0.000 claims description 2
- 238000005457 optimization Methods 0.000 abstract description 2
- 238000010248 power generation Methods 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 description 14
- 238000004364 calculation method Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 3
- 230000007935 neutral effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000009141 biological interaction Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 238000013524 data verification Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 210000001944 turbinate Anatomy 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Computer Hardware Design (AREA)
- General Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了属于新能源风力发电技术领域的一种风力机远场尾流流向湍流度计算方法。包括步骤:1,获取风电机组入流情况、风电机组参数和风电机组运行状态的基本数据;2,根据获取的基本数据计算轮毂高度水平面内尾流速度损失剖面标准差σ随下游距离x的线性变化函数;在尾流速度损失分布已知的情况下,或通过高斯速度损失剖面拟合直接获取,相当于已知的输入条件;3,将轮毂高度水平面内尾流速度损失剖面标准差σ作为输入值输入到附加湍流度模型中,结合风电机组入流数据得出风电机组远场尾流流向湍流度的预测结果。本发明可以实现风电机组远场尾流流向湍流度的准确预测,对机组排布优化具有重要的指导意义。
Description
技术领域
本发明属于新能源风力发电技术领域,特别涉及一种风力机远场尾流流向湍流度计算方法。
背景技术
流向湍流度是导致风力发电机组疲劳失效的主要原因,以此用来计算风力发电机组的疲劳载荷。因此,利用结构简单、计算效率高的解析模型准确预测风力发电机组尾流区的流向湍流度成为风电企业的迫切要求;有助于提高风电场微观选址的技术水平,降低风力机疲劳载荷、延长风力机使用寿命。
理论上风力发电机组尾流区的流向湍流度可以用以下公式表达:
其中,Iu为风电机组尾流区中的流向湍流度;I0为环境湍流度;ΔIu为风电机组运行时产生的附加流向湍流度;它们都是无量纲量。
在上式中,通常环境湍流度I0分布是已知的,重点是通过对附加流向湍流度的预测,实现尾流区内流向湍流度的计算。而尾流区的附加流向湍流度分布,受到风力发电机组推力系数CT、风力发电机组轮毂高度位置环境湍流度Ia、与机组的轴向距离x、径向距离r等多个因素影响;建立精确的附加流向湍流度计算模型,是提高尾流区流向湍流度计算精度的核心。
在附加湍流度模型中,主要分为三类(Porté-Agel F*,Bastankhah M,ShamsoddinS.Wind-Turbine and Wind-Farm Flows:A Review[J].Boundary Layer Meteorology,2019,174(1):1-59):第一类是横截面附加流向湍流度呈顶帽分布的Frandsen湍流度模型(Frandsen S,Thgersen ML.Integrated fatigue loading for wind turbines in windfarms by combining ambient turbulence and wakes[J].Multi-Science PublishingCo.Ltd.1999,23(6):327–339);第二类是只考虑上叶尖位置附加流向湍流度变化的Crespo及其类似的湍流度模型(Crespo A,J.Herna′ndez J.Turbulence characteristics inwind-turbine wakes[J].Journal of Wind Engineering and IndustrialAerodynamics,1996,61(1):71-85);第三类是横截面附加流向湍流度呈双高斯分布,附加流向湍流度最大值在叶尖处的Takeshi湍流度模型(Ishihara T,Qian GW*.A newGaussian-based analytical wake model for wind turbines considering ambientturbulence intensities and thrust coefficient effects[J].Journal of WindEngineering and Industrial Aerodynamics,2018,177:275-292)。
其中,Frandsen湍流度模型未考虑与环境湍流度的映射关系,附加流向湍流度的预测局限于轮毂中心位置处,尾流截面处的附加流向湍流度呈顶帽分布与工程实际明显不符。
Crespo及其类似的湍流度模型大多局限于机组上叶尖位置处的附加流向湍流度,并未对尾流截面处的湍流分布情况加以考虑。
Takeshi湍流度模型的思路是:先假定附加流向湍流度沿轮毂中心线呈双高斯轴对称分布,最大值位置位于叶尖处,再考虑来流剪切层的影响对垂直方向上的不对称性进行修正。其对附加流向湍流度的分析较为全面,但存在以下问题:
(1)未考虑尾流发展过程中膨胀特性对附加流向湍流度最大值位置的影响,只在小湍流度、小推力系数,尾流膨胀较小,可近似为不膨胀的情况下适用;在其它情况下,则存在明显误差,由此进行的高斯分布也将因为顶点位置的不准确出现较大偏差。
(2)在对垂直方向上附加流向湍流度的不对称性进行修正时,未考虑轮毂高度以上来流剪切对不对称性的影响,且地面修正范围未在展向y方向加以限制,认为只与垂直高度z有关,与地面对风电机组附加流向湍流度的影响仅在机组尾流区中的特性明显不符。
(3)数据验证时所采用的入流湍流度为3.5%、机组推力系数为0.36等部分实验过于理想;在工程应用中具有很大的局限性。
综上所述,上述附加湍流度模型能全面引起风电机组流向湍流度的预测结果不准确问题。
发明内容
本发明的目的是提出一种风力机远场尾流流向湍流度计算方法,其特征在于,所述方法包括:
步骤1:获取风电机组入流情况、风电机组参数和风电机组运行状态的基本数据;
步骤2:根据步骤1获取的基本数据计算轮毂高度水平面内尾流速度损失剖面标准差σ随下游距离x的线性变化函数;在尾流速度损失分布已知的情况下,或通过高斯速度损失剖面拟合直接获取,相当于已知的输入条件;
步骤3:将步骤1的基本数据和步骤2中轮毂高度水平面内尾流速度损失剖面标准差σ作为输入值输入到附加湍流度模型中,结合风电机组入流数据得出风电机组远场尾流流向湍流度的预测结果。
所述步骤1获取风电机组基本数据具体包括:风电机组来流环境湍流度I0的分布、风电机组轮毂高度位置环境湍流度Ia、轮毂高度zh、风轮直径D、机组推力系数CT。
所述步骤2获取的基本数据计算轮毂高度水平面内尾流速度损失剖面标准差σ随下游距离x的线性变化函数具体包括:
步骤201:将风电机组轮毂高度位置环境湍流度Ia、风轮直径D、机组推力系数CT输入到尾流速度损失剖面标准差σ的求解公式中,计算σ=kwx+σ0,
所述步骤3将轮毂高度水平面内尾流速度损失剖面标准差σ作为输入值输入到附加湍流度模型中,结合风电机组入流数据得出风电机组远场尾流流向湍流度的预测结果具体包括:
步骤301:根据附加流向湍流度模型中双高斯分布假设,以轮毂高度水平面内尾流速度损失剖面标准差σ作为输入,求取附加流向湍流度最大值位置与轮毂中心线的径向距离r1/2、附加流向湍流度呈高斯分布时的标准偏差σT,其中
步骤302:将附加流向湍流度最大值位置与轮毂中心线的径向距离r1/2、附加流向湍流度呈高斯分布时的标准偏差σT、风电机组轮毂高度位置环境湍流度Ia、轮毂高度zh、风轮直径D、机组推力系数CT代入到附加湍流度模型,得出风电机组尾流附加流向湍流度的预测结果;
所述附加湍流度模型表达式为:
所述修正函数表达式为:
α表示尾流截面内任意位置坐标与轮毂中心轴的连线相对于轮毂高度水平面的夹角;P为权重分布,表示尾流截面内不同位置处修正函数的权重变化;z表示距地面的垂直距离,与轮毂高度zh单位一致。
所述权重分布P表达式为:
步骤304:根据风电机组远场尾流流向湍流度的预测结果,评估风电机组之间尾流效应对疲劳载荷的影响,基于相关的疲劳载荷评估规范,判断风电机组是否满足疲劳载荷要求。
步骤305:当风电机组的疲劳载荷要求不满足时,对风电机组的布机位置进行调整;再重复步骤304,直到要求满足为止;从而延长风力机使用寿命。
本发明的有益效果在于在不考虑热浮力效应影响中性大气条件为发生在早晨、傍晚或多云的天气下、平坦地形、风电机组不偏航情况下,风电机组轮毂高度位置环境湍流度0.06<Ia<0.14、机组推力系数0.6<CT<0.84,可以实现风电机组远场尾流流向湍流度的准确预测,对机组排布优化具有重要的指导意义。
附图说明
图1是附加流向湍流度双高斯分布假设示意图。
图2是一种风力机远场尾流流向湍流度计算方法的应用流程图。
图3是各湍流度解析模型与刘惠文进行风洞实验所得数据的对比图:其中3D、4D、5D、6D、7D、8D、9D、10D分别表示相应的机组下游位置,即x=3D、4D、5D、6D、7D、8D、9D、10D。
图4是各湍流度解析模型与Takeshi所得实验数据的对比图:(a)case1、(d)case2、(g)case3轮毂高度水平面附加流向湍流度最大值随下游距离x的变化;(b)case1、(e)case2、(h)case3经过轮毂中心竖直面x-z内不同下游位置处(x/D=2,4,6,8)的流向湍流度;(c)case1、(f)case2、(i)case3轮毂高度处x-y水平面内不同下游位置处(x/D=2,4,6,8)的流向湍流度。
图5是各湍流度解析模型与Shengbai Xie进行风洞实验所得实验数据的对比图:其中5D、7D、10D、14D、20D分别表示相应的机组下游位置,即x=5D、7D、10D、14D、20D。
图6是各湍流度解析模型与Wu YT所得大涡模拟实验数据的对比图:(a)case1,(c)case2,(e)case3,(g)case4在经过轮毂中心x-z竖直面内的对比;(b)case1,(d)case2,(f)case3,(h)case4在轮毂高度处x-y水平面内的比较。其中5D、7D、10D、15D分别表示相应的机组下游位置,即x=5D、7D、10D、15D。
具体实施方式
本发明提出一种风力机远场尾流流向湍流度计算方法,包括如下步骤:
步骤1:获取风电机组入流情况、风电机组参数和风电机组运行状态的基本数据;
步骤2:根据步骤1获取的基本数据计算轮毂高度水平面内尾流速度损失剖面标准差σ随下游距离x的线性变化函数;在尾流速度损失分布已知的情况下,或通过高斯速度损失剖面拟合直接获取,相当于已知的输入条件;
步骤3:将步骤1的基本数据和步骤2中轮毂高度水平面内尾流速度损失剖面标准差σ作为输入值输入到附加湍流度模型中,结合风电机组入流数据得出风电机组远场尾流流向湍流度的预测结果。下面结合附图和优选实施例对本发明作详细说明。
首先根据图1所示的附加流向湍流度模型中双高斯分布假设,以轮毂高度水平面内尾流速度损失剖面标准差σ作为输入,求取附加流向湍流度最大值位置与轮毂中心线的径向距离r1/2、附加流向湍流度呈高斯分布时的标准偏差σT,其中
图2为一种风力机远场尾流流向湍流度计算方法的应用流程图,以下结合实施例对所述方法的具体过程进行详细说明。
实施例1
在风洞实验风力机模型尾流中进行一种风力机远场尾流流向湍流度计算方法的应用,包括以下步骤:
步骤1:确立参考坐标系,以机组安装位置为坐标原点,平行于风轮旋转轴为x轴(平行于来流方向),垂直于来流方向为y轴,竖直方向为z轴。
步骤2:获取风电机组入流情况、风电机组参数、风电机组运行状态等基本数据。具体包括:
风力机来流环境湍流度I0的分布、风力机轮毂高度位置环境湍流度Ia=0.105、轮毂高度zh=125mm、风轮直径D=120mm、机组推力系数CT=0.5。
步骤3:根据步骤2获取的基本数据计算轮毂高度水平面内尾流速度损失剖面标准差σ随下游距离x的线性变化函数。具体包括:
步骤32:将系数kw=0.044、σ0=0.220D、以及风轮直径D输入到尾流速度损失剖面标准差σ的求解公式得到:σ=0.044x+0.220D;该公式可对风电机组下游任意位置x处的尾流速度损失剖面标准差σ进行求解。
步骤4:将步骤2的基本数据和步骤3中轮毂高度水平面内尾流速度损失剖面标准差σ作为输入值输入到附加湍流度模型中,结合风力机入流数据得出风力机远场尾流流向湍流度的预测结果。具体包括:
步骤42:将附加流向湍流度最大值位置与轮毂中心线的径向距离附加流向湍流度呈高斯分布时的标准偏差风力机轮毂高度位置环境湍流度Ia=0.105、轮毂高度zh=125mm、风轮直径D=120mm、机组推力系数CT=0.5代入到附加湍流度模型,即可计算得到整个尾流区内任意位置(x,y,z)处的附加流向湍流度ΔIu。
所述的附加湍流度模型表达式为:
其中:
r表示尾流截面内任意位置坐标到轮毂中心线的直线距离,它与下游距离x、轮毂高度zh、风轮直径D单位一致;δ(r)表示修正函数,对垂直方向上附加流向湍流度不对称性进行修正;α表示尾流截面内任意位置坐标与轮毂中心轴的连线相对于轮毂高度水平面的夹角,P为权重分布,表示尾流截面内不同位置处修正函数的权重变化。
步骤43:以风力机来流环境湍流度I0(z)、风力机附加流向湍流度的预测结果ΔIu(x,y,z)作为输入;判断ΔIu(x,y,z)≥0是否成立,如果成立,则令否则,则令得出整个尾流区内任意位置(x,y,z)处风力机流向湍流度的预测结果。
步骤5:将计算得到的结果与刘惠文进行风洞实验的结果、Frandsen湍流度模型、Crespo湍流度模型、Takeshi湍流度模型进行对比,为便于对比,Frandsen、Crespo湍流度模型均采用顶帽分布且以rw=2r1/2=2.35σ作为边界;对比结果如图3所示。
实施例2
本实施例采用的分析步骤与实施例1相同;尾流速度损失剖面标准差σ随下游距离x的线性变化函数在实施例1中已给出,σ当作已知条件,其余的分析步骤与实施例1相同;以下在实施例2、实施例3和实施例4中分别进行结果分析。
本实施例是在Takeshi流向湍流度风洞实验与数值计算实验结果中(case1-3)进行应用;其中Exp表示风洞实验,WT-M、WT-P为参考Exp实验条件采用旋转制动盘机组模型计算得到的大涡模拟数值结果;WT-M里的机组尺寸与Exp一致,轮毂高度zh=0.7m、风轮直径D=0.57m;WT-P采用的是容量为2.4MW的海上风电机组,zh=80m、风轮直径D=92m。同一案例中各实验入流条件基本一致:case1中,Ia=0.035、CT=0.81;case2中,Ia=0.137、CT=0.37;case3中,Ia=0.137、CT=0.81。
机组尺寸的影响按风轮直径D归一化之后可忽略,因此所有模型均以WT-M中的机组尺寸zh=0.7m、D=0.57m以及各案例中相应的入流条件作为输入,对比结果如图4所示。
实施例3
在Shengbai Xie风洞实验风力机模型尾流中进行应用;风洞实验中,风力机轮毂高度位置环境湍流度Ia=0.070、轮毂高度zh=125mm、风轮直径D=150mm、机组推力系数CT=0.42,尾流速度损失剖面标准差σ=0.023x+0.219D。
所有模型均以Ia=0.070、zh=125mm、D=150mm、CT=0.42、σ=0.023x+0.219D作为输入,对比结果如图5所示。
实施例4
在Wu YT采用旋转制动盘机组模型进行的大涡模拟附加流向湍流度数值计算实验中(case1-4)进行应用。各实验风电机组参数、运行状态完全一致,zh=70m、D=80m、CT=0.8。区别在于各实验的入流条件和尾流速度损失剖面标准差:case1中,Ia=0.134、σ=0.055x+0.238D;case2中,Ia=0.094、σ=0.040x+0.253D;case3中,Ia=0.069、σ=0.030x+0.272D;case4中,Ia=0.048、σ=0.031x+0.257D。在附加流向湍流度已知的情况下,直接对附加流向湍流度进行对比,对比结果更加直观,更容易观测出模型预测的准确性。
以相应的已知条件作为模型输入,各模型与实验数据的对比结果如图6所示。
Claims (4)
1.一种风力机远场尾流流向湍流度计算方法,其特征在于,所述方法包括:
步骤1:获取风电机组入流情况、风电机组参数和风电机组运行状态的基本数据;
步骤2:根据步骤1获取的基本数据计算轮毂高度水平面内尾流速度损失剖面标准差σ随下游距离x的线性变化函数;在尾流速度损失分布已知的情况下,或通过高斯速度损失剖面拟合直接获取,相当于已知的输入条件;
步骤3:将步骤1的基本数据和步骤2中轮毂高度水平面内尾流速度损失剖面标准差σ作为输入值输入到附加湍流度模型中,结合风电机组入流数据得出风电机组远场尾流流向湍流度的预测结果。
2.根据权利要求1所述一种风力机远场尾流流向湍流度计算方法,其特征在于,所述步骤1获取风电机组基本数据具体包括:风电机组来流环境湍流度I0的分布、风电机组轮毂高度位置环境湍流度Ia、轮毂高度zh、风轮直径D、机组推力系数CT。
4.根据权利要求1所述一种风力机远场尾流流向湍流度计算方法,其特征在于,所述步骤3将轮毂高度水平面内尾流速度损失剖面标准差σ作为输入值输入到附加湍流度模型中,结合风电机组入流数据得出风电机组远场尾流流向湍流度的预测结果具体包括:
步骤301:根据附加流向湍流度模型中双高斯分布假设,以轮毂高度水平面内尾流速度损失剖面标准差σ作为输入,求取附加流向湍流度最大值位置与轮毂中心线的径向距离r1/2、附加流向湍流度呈高斯分布时的标准偏差σT,其中
步骤302:将附加流向湍流度最大值位置与轮毂中心线的径向距离r1/2、附加流向湍流度呈高斯分布时的标准偏差σT、风电机组轮毂高度位置环境湍流度Ia、轮毂高度zh、风轮直径D、机组推力系数CT代入到附加湍流度模型,得出风电机组尾流附加流向湍流度的预测结果;
所述附加湍流度模型表达式为:
所述修正函数表达式为:
α表示尾流截面内任意位置坐标与轮毂中心轴的连线相对于轮毂高度水平面的夹角;P为权重分布,表示尾流截面内不同位置处修正函数的权重变化;z表示距地面的垂直距离,与轮毂高度zh单位一致;
所述权重分布P表达式为:
步骤304:根据风电机组远场尾流流向湍流度的预测结果,评估风电机组之间尾流效应对疲劳载荷的影响,基于相关的疲劳载荷评估规范,判断风电机组是否满足疲劳载荷要求;
步骤305:当风电机组的疲劳载荷要求不满足时,对风电机组的布机位置进行调整;再重复步骤304,直到要求满足为止;从而延长风力机使用寿命。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011102957.9A CN112347611B (zh) | 2020-10-15 | 2020-10-15 | 一种风力机远场尾流流向湍流度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011102957.9A CN112347611B (zh) | 2020-10-15 | 2020-10-15 | 一种风力机远场尾流流向湍流度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112347611A true CN112347611A (zh) | 2021-02-09 |
CN112347611B CN112347611B (zh) | 2024-05-31 |
Family
ID=74360669
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011102957.9A Active CN112347611B (zh) | 2020-10-15 | 2020-10-15 | 一种风力机远场尾流流向湍流度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112347611B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113111608A (zh) * | 2021-05-10 | 2021-07-13 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种新型局部湍流脉动生成方法 |
CN113239648A (zh) * | 2021-06-22 | 2021-08-10 | 华北电力大学 | 风电场远场尾流流向湍流度确定方法及装置 |
CN113536487A (zh) * | 2021-08-02 | 2021-10-22 | 华能新能源股份有限公司 | 一种基于双变量高斯函数的尾流计算方法、装置及存储介质 |
CN113627097A (zh) * | 2021-07-22 | 2021-11-09 | 中国长江三峡集团有限公司 | 一种利用风力机scada数据校正尾流评估模型的方法 |
CN113705118A (zh) * | 2021-08-20 | 2021-11-26 | 南京航空航天大学 | 一种风力机尾流湍流强度的计算方法 |
CN116050287A (zh) * | 2022-12-12 | 2023-05-02 | 中广核风电有限公司 | 海上漂浮式风机尾流解析的建模方法和装置 |
CN117272869A (zh) * | 2023-11-15 | 2023-12-22 | 南京航空航天大学 | 一种兼顾风力机近尾流和远尾流特性的全尾流解析方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2840258A1 (en) * | 2013-08-20 | 2015-02-25 | General Electric Company | System and method for preventing excessive loading on a wind turbine |
CN106897486A (zh) * | 2017-01-12 | 2017-06-27 | 华北电力大学 | 考虑湍流强度影响的抛物线形风电机组尾流模型计算方法 |
CN108108562A (zh) * | 2017-12-27 | 2018-06-01 | 华北电力大学 | 一种基于高斯分布的风电机组尾流的解析建模方法 |
CN109885845A (zh) * | 2017-12-06 | 2019-06-14 | 北京普华亿能风电技术有限公司 | 一种风力发电机的尾流分析方法及*** |
CN110009736A (zh) * | 2019-05-09 | 2019-07-12 | 华北电力大学(保定) | 三维尾流模型的建立方法、装置、设备及存储介质 |
CN110046442A (zh) * | 2019-04-22 | 2019-07-23 | 国电联合动力技术有限公司 | 风机与尾流投影面积叠加计算方法及装置 |
EP3527817A1 (de) * | 2018-02-19 | 2019-08-21 | Siemens Aktiengesellschaft | Verfahren und vorrichtungen zur automatischen ermittlung und/oder kompensation des einflusses einer wirbelschleppe auf eine windkraftanlage |
CN110159488A (zh) * | 2019-06-26 | 2019-08-23 | 天津超算科技有限公司 | 提高发电效率的风机*** |
DE102018108858A1 (de) * | 2018-04-13 | 2019-10-17 | Wobben Properties Gmbh | Windenergieanlage, Windpark sowie Verfahren zum Regeln einer Windenergieanlage und eines Windparks |
-
2020
- 2020-10-15 CN CN202011102957.9A patent/CN112347611B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2840258A1 (en) * | 2013-08-20 | 2015-02-25 | General Electric Company | System and method for preventing excessive loading on a wind turbine |
CN106897486A (zh) * | 2017-01-12 | 2017-06-27 | 华北电力大学 | 考虑湍流强度影响的抛物线形风电机组尾流模型计算方法 |
CN109885845A (zh) * | 2017-12-06 | 2019-06-14 | 北京普华亿能风电技术有限公司 | 一种风力发电机的尾流分析方法及*** |
CN108108562A (zh) * | 2017-12-27 | 2018-06-01 | 华北电力大学 | 一种基于高斯分布的风电机组尾流的解析建模方法 |
EP3527817A1 (de) * | 2018-02-19 | 2019-08-21 | Siemens Aktiengesellschaft | Verfahren und vorrichtungen zur automatischen ermittlung und/oder kompensation des einflusses einer wirbelschleppe auf eine windkraftanlage |
DE102018108858A1 (de) * | 2018-04-13 | 2019-10-17 | Wobben Properties Gmbh | Windenergieanlage, Windpark sowie Verfahren zum Regeln einer Windenergieanlage und eines Windparks |
CN110046442A (zh) * | 2019-04-22 | 2019-07-23 | 国电联合动力技术有限公司 | 风机与尾流投影面积叠加计算方法及装置 |
CN110009736A (zh) * | 2019-05-09 | 2019-07-12 | 华北电力大学(保定) | 三维尾流模型的建立方法、装置、设备及存储介质 |
CN110159488A (zh) * | 2019-06-26 | 2019-08-23 | 天津超算科技有限公司 | 提高发电效率的风机*** |
Non-Patent Citations (4)
Title |
---|
J.M. PROSPATHOPOULOS ETC: "Evaluation of the effects of turbulence model enhancements on wind turbine wake predictions", WIND ENERGY, vol. 14, no. 02, pages 285 - 300 * |
刘永前 等: "大型风电机组尾流效应非定常数值模拟研究", 工程热物理学报, vol. 38, no. 03, pages 541 - 547 * |
吴正人 等: "风电场运行对局地边界层气象因素影响", 太阳能学报, no. 02, pages 118 - 125 * |
孟航: "单台风电机组远场尾流CFD计算模型研究", 《中国优秀硕士学位论文全文数据库 工程科技II 辑》, no. 01, pages 042 - 534 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113111608A (zh) * | 2021-05-10 | 2021-07-13 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种新型局部湍流脉动生成方法 |
CN113111608B (zh) * | 2021-05-10 | 2022-06-28 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种新型局部湍流脉动生成方法 |
CN113239648A (zh) * | 2021-06-22 | 2021-08-10 | 华北电力大学 | 风电场远场尾流流向湍流度确定方法及装置 |
CN113239648B (zh) * | 2021-06-22 | 2023-11-21 | 华北电力大学 | 风电场远场尾流流向湍流度确定方法及装置 |
CN113627097A (zh) * | 2021-07-22 | 2021-11-09 | 中国长江三峡集团有限公司 | 一种利用风力机scada数据校正尾流评估模型的方法 |
CN113627097B (zh) * | 2021-07-22 | 2022-10-04 | 中国长江三峡集团有限公司 | 一种利用风力机scada数据校正尾流评估模型的方法 |
CN113536487A (zh) * | 2021-08-02 | 2021-10-22 | 华能新能源股份有限公司 | 一种基于双变量高斯函数的尾流计算方法、装置及存储介质 |
CN113705118A (zh) * | 2021-08-20 | 2021-11-26 | 南京航空航天大学 | 一种风力机尾流湍流强度的计算方法 |
CN116050287A (zh) * | 2022-12-12 | 2023-05-02 | 中广核风电有限公司 | 海上漂浮式风机尾流解析的建模方法和装置 |
CN116050287B (zh) * | 2022-12-12 | 2023-12-08 | 中广核风电有限公司 | 海上漂浮式风机尾流解析的建模方法和装置 |
CN117272869A (zh) * | 2023-11-15 | 2023-12-22 | 南京航空航天大学 | 一种兼顾风力机近尾流和远尾流特性的全尾流解析方法 |
CN117272869B (zh) * | 2023-11-15 | 2024-02-09 | 南京航空航天大学 | 一种兼顾风力机近尾流和远尾流特性的全尾流解析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112347611B (zh) | 2024-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112347611A (zh) | 一种风力机远场尾流流向湍流度计算方法 | |
CN107194097B (zh) | 基于风电场气动模拟和风速风向数据的分析方法 | |
Schepers | ENDOW: Validation and improvement of ECN's wake model | |
CN112784509B (zh) | 基于改进高斯分布尾流模型的风力发电机功率计算方法 | |
CN110009736B (zh) | 三维尾流模型的建立方法、装置、设备及存储介质 | |
CN114169614B (zh) | 一种基于风力机尾流模型优化的风电场优化调度方法及*** | |
CN104794287A (zh) | 一种风工程尾流计算方法 | |
CN112966454B (zh) | 一种风电场风机尾流动态耦合模拟方法 | |
Park et al. | Wind farm layout optimization using genetic algorithm and its application to Daegwallyeong wind farm | |
CN113627101A (zh) | 一种基于改进型ad/rsm模型的风力机尾流模拟方法 | |
Chen et al. | High-order redesign method for wind turbine blade optimization in model test considering aerodynamic similarity | |
Sugathapala et al. | Aerodynamic modeling of simplified wind turbine rotors targeting small-scale applications in Sri Lanka | |
Rodrigo et al. | Validation of meso-wake models for array efficiency prediction using operational data from five offshore wind farms | |
CN109766611B (zh) | 一种考虑地形坡度的风电场功率简化预测方法 | |
Tang et al. | An on-site measurement coupled CFD based approach for wind resource assessment over complex terrains | |
Rafiei et al. | Optimized configuration with economic evaluation for shrouded vertical axis wind turbines applicable for urban structures | |
CN113627097B (zh) | 一种利用风力机scada数据校正尾流评估模型的方法 | |
Maniaci et al. | Wind energy high-fidelity model verification and validation roadmap | |
Hwang et al. | Wind farm layout optimization using multidisciplinary model | |
Bekken et al. | Performance enhancement of an induced draught axial flow fan through pressure recovery | |
Ivanell et al. | Validation of methods using EllipSys3D | |
CN117313399B (zh) | 一种适用于复杂地形的水平轴风力机三维各向异性超高斯全尾流模型的建立及应用方法 | |
Höning et al. | Influence of rotor blade flexibility on the near-wake behavior of the NREL 5 MW wind turbine | |
Porter et al. | Comparative study of numerical modelling techniques to estimate tidal turbine blade loads | |
CN117556170A (zh) | 一种偏航风机的非高斯尾流风速损失分布预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |