基于CFD技术的垂直轴风力机动态尾流特性研究

利用CFD软件对麦克马斯特大学垂直轴风力机进行不同叶尖速比下的数值模拟,计算结果与风洞试验数据吻合良好。近场尾流中,与单叶片的风力机模拟结果比较,上游叶片产生并向下游延伸的旋涡影响下游运行轨道上叶片的升阻力特性,不仅使叶片扭矩输出峰值降低,而且峰值产生的时间延迟。对垂直轴风力机叶片叶梢进行修改,模拟结果显示,叶片扭矩输出峰值不变,但是谷值有所降低,修改后风力机沿风向推力幅值降低明显;远场尾流中,采用风速轮廓线原理,以瑞典的法尔肯贝里市200kW垂直轴风力机为原型,按照真实的空间排布进行数值模拟。模拟结果显示,上游风力机上下两端处产生较为集中的远场尾流,影响下游风力机叶片下半段的气动性能,下游风力机功率输出降低明显。     

浮式海上风力机叶片气动性能的流固耦合分析

随着深水浮式海上风电场在世界范围内的兴起,浮式平台运动性能对风力机稳定运行及叶片气动载荷影响的研究具有重要意义。基于三维粘性不可压缩Navier-Stokes方程和适用于旋转流场分析的重整化群k-ε(RNG)湍流模型,数值模拟美国可再生能源实验室(NREL)5MW海上风力机的气动性能,并将数值模拟结果与NREL的设计参考数据进行对比分析,较好地验证了该数值模拟方法的有效性。进一步利用滑移网格技术模拟风力机叶片随浮式平台的典型周期性运动,实现浮式风力机叶片与周围流场的复杂非线性流固耦合分析,分别研究浮式平台不同运动幅值和运动周期对风力机叶片气动性能的影响规律,并从物理机理角度进行阐明分析。本文的主要研究成果,将对未来大型深水浮式海上风力机的气动性能分析及浮式平台系统的运动性能设计起到积极的指导作用。     

浮式海上风力机叶片气动性能的流固耦合分析

随着深水浮式海上风电场在世界范围内的兴起,浮式平台运动性能对风力机稳定运行及叶片气动载荷影响的研究具有重要意义。基于三维粘性不可压缩Navier-Stokes方程和适用于旋转流场分析的重整化群k-ε（RNG）湍流模型,数值模拟美国可再生能源实验室（NREL）5MW海上风力机的气动性能,并将数值模拟结果与NREL的设计参考数据进行对比分析,较好地验证了该数值模拟方法的有效性。进一步利用滑移网格技术模拟风力机叶片随浮式平台的典型周期性运动,实现浮式风力机叶片与周围流场的复杂非线性流固耦合分析,分别研究浮式平台不同运动幅值和运动周期对风力机叶片气动性能的影响规律,并从物理机理角度进行阐明分析。本文的主要研究成果,将对未来大型深水浮式海上风力机的气动性能分析及浮式平台系统的运动性能设计起到积极的指导作用。     

基于Navier-Stokes／Kirchhoff方程的悬停旋翼跨声速气动声学计算

研究了Kirchhoff积分面是否有盖有底,以及是否计及旋翼网格上的流场值,这两个因素对噪声预测结果的影响.发展了一种基于重叠网格的计算悬停旋翼远场噪声的数值方法.数值计算过程分为流场模拟和声场模拟两部分.悬停旋翼流场的数值模拟是在两个相互重叠的网格上进行的:在高质量的旋翼网格上求解Navier-Stokes方程,用于模拟旋翼附近的粘性流动和近场尾涡的捕捉;在远离粘性区域处布置符合悬停流场物理特征的圆柱形背景网格,控制方程为Euler方程,用于远场尾涡的捕捉.计算得到的流场信息插值到用于声场计算的Kirchhoff积分面上.观测点处的噪声可以认为是由这个完全包含桨叶的Kirchhoff积分面上的面元(声源)发声得到.远场声波的传播由Kirchhoff积分公式描述.计算结果表明:采用有盖有底的Kirchhoff积分面并且同时计及旋翼网格流场值时,计算得到的HSI噪声与实验值吻合最好.     

一种风力机气动计算的全自由涡尾迹模型

采用全自由方式建立风力机尾流场的涡尾迹模型,引入“虚拟周期”的概念,并发展一种自适应松弛因子方法,从而改善了自由尾迹迭代的稳定性,提高了迭代收敛速度。利用建立的自由涡尾迹模型,计算了风力机叶片的尾流场结构、气动性能及叶片载荷,并与实验结果进行了对比分析。结果表明,尖速比越大,自适应松弛因子方法对缩小模型计算时间越有效;全自由涡尾迹模型能准确给出风力机尾流场的结构,包括尾迹的扩张以及叶尖涡和叶根涡的产生、发展和耗散的过程,风轮扭矩与实验数据吻合;叶片载荷分布的计算结果在低风速下与实验值基本一致,但是在大风速下差别较大,说明需要一个准确的失速模型。     

非结构动网格在三维可动边界问题中的应用

研究用于非结构动网格的弹簧近似方法，采用顶点弹簧描述，导出并讨论了弹簧倔强系数的取值。通过引入边界修正和扭转效应修正，对标准弹簧近似方法进行了改进，转动翼型算例的结果表明，改进后的方法大大提高了网格变形能力和网格质量，应用该动网格方法耦合求解基于（Arbitrary Lagrangian-Eulerian,ALE）描述的三维Euler方程，模拟了作俯抑振动的矩形机翼绕流，计算结果与实验数据及文献计算结果十分一致，作为多个自由刚与流体耦合运动问题的简单例证，耦合刚体动力学方程，模拟了激波与双立方体的相互作用，得到了非定常流场结构，研究表明，基于弹簧近似的非结构动网格与有限体积流式流场解算器相结合，是模拟包含运动边界的非定常流动问题的有效方法。     

串列叶片式前向离心风机气动与噪声特性的优化研究

对采用串列叶片的某前向离心风机内部三维非定常流动进行了数值计算,重点研究了串列叶片不同叶片相对长度和不同叶片相对周向位置两个参数对风机气动性能及气动噪声的影响.通过响应面方法对数值结果进行二次回归拟合,得到两个参数与风机效率和A声级间的函数关系,并进行了优化分析.数值结果表明:两个参数对串列叶片式前向离心风机效率和A声级均有较大影响,合理的串列叶片设计能够在保持气动性能基本不变的情况下降低风机的气动噪声.将可靠的CFD数值技术与响应面方法结合起来用于指导离心风机的改进及试验设计是可行的,本文的研究结果可为串列式离心风机在节能与降噪的总体设计方面提供参考.     

方柱尾流流态对床面冲淤形态的影响

应用平面二维悬沙数学模型对方柱尾流区的泥沙输运及床面调整进行了数值模拟，并对重要物理参数做了分析说明。数值计算采用时间分裂一全隐式有限差分格式，流场的计算基于沿水深平均的RANS方程。通过系统的数值模拟，揭示了不同来流情况下，方柱尾流区的流态变化，及与之相应的床面变化规律，并对不同流态下的泥沙运动作了分析。计算结果表明在方柱尾流区，不同的流场流态对尾流区的床面调整有明显影响。流场较弱的情况下，尾流区中湍动强度相应较弱，此时悬浮泥沙由于流速减小而普遍落淤，床面应力的减小也致使床面冲刷量减小。随着流场强度的逐渐增大，尾流区中湍动强度相应增强，床面应力增大，同时湍流的强烈交换作用增强了对泥沙的输运作用，床面变化不再是简单的淤积状态，部分区域出现了冲刷。     

非结构动网格在多介质流体数值模拟中的应用

采用非结构动网格方法对含多介质的流场进行数值模拟.采用改进的弹簧方法来处理由于边界运动而产生的网格变形.采用基于格心的有限体积方法求解守恒型的ALE(Arbitrary Lagrangiall-Eulerian)方程,控制面通量的计算采用HLLC(Hartem,Lax,van Leer,Contact)方法,采用几何构造的方法使空间达到二阶精度,时间离散采用四阶Runge-Kutta方法.物质界面的处理采用虚拟流体方法.本文对含动边界的激波管、水下爆炸等流场进行数值模拟,取得较好的结果,不同时刻界面的位置和整个扩张过程被准确模拟.     

计算含动边界非定常流动的无网格算法

在无网格算法中考虑了含动边界的流动问题,研究了可以计算处理包含一定位移及扭转动边界非定常流动的算法.创建了无网格算法的动点法则,并引入抗扭方法对弹簧方法进行改进来处理离散点运动,提高了方法的可用度及精度.发展了求解基于无网格的ALE方程组的算法,在点云离散的基础上采用曲面逼近计算空间导数及HLLC格式计算数值通量,运用四步龙格-库塔法进行时间推进.在跨、超音速条件下,计算模拟了典型翼型简谐振动流场,计算结果与实验结果及文献对比吻合,验证了该算法的正确性.     

Effects of a rooftop wind turbine on the dispersion of air pollutant behind a cube-shaped building

The concentration distribution of urban air pollutants is closely related to people's health. As an important utilization form of urban wind power, rooftop wind turbines have been widely used in cities. The wake effect of the rooftop wind turbines will change the flow behind buildings and then affect the pollutant dispersion. To this end, the pollutant dispersion behind the building is studied via the computational fluid dynamics method. The actuator disk model and idealized cube are adopted to model the wind turbine and the building, respectively. The study shows that the rooftop wind turbine can reduce the pollutant mass fraction near the ground and the pedestrian level. Due to the wake effect of the rooftop wind turbine, the turbulent fluctuation behind the building is weakened, and the spanwise pollutant dispersion is suppressed. Besides, the rooftop wind turbine weakens the downwash movement of the building, which enhances the vertical pollutant dispersion.     

Progresses in application of computational fluid dynamic methods to large scale wind turbine aerodynamics

The computational fluid dynamics (CFD) methods are applied to aerodynamic problems for large scale wind turbines. The progresses including the aerodynamic analyses of wind turbine profiles, numerical flow simulation of wind turbine blades, evaluation of aerodynamic performance, and multi-objective blade optimization are discussed. Based on the CFD methods, significant improvements are obtained to predict two/three dimensional aerodynamic characteristics of wind turbine airfoils and blades, and the vortical structure in their wake flows is accurately captured. Combining with a multi-objective genetic algorithm, a 1.5 MW NH-1500 optimized blade is designed with high efficiency in wind energy conversion.     

单一水平轴风电机组尾迹的模拟方法与流动机理研究综述

为实现碳达峰、碳中和“3060”目标, 风能将在我国能源体系发挥重要作用. 风力机尾迹是影响风电性能和度电成本的关键因素, 需在风力机布置和控制设计中充分考虑. 本文首先介绍风力机尾迹的数值模拟方法, 包括解析模型、低阶模型、大涡模拟和来流湍流生成方法. 解析模型和低阶模型可快速计算风力机尾迹, 但依赖于模型参数, 且不能或不能准确预测尾迹湍流特性. 结合风力机参数化模型的大涡模拟可准确预测尾迹蜿蜒等湍流特征, 是流动机理研究的有力工具, 可为发展快速预测模型提供数据和理论支撑. 接着, 本文介绍了叶尖涡、中心涡和尾迹蜿蜒并讨论其产生机理. 对于湍流来流, 叶尖涡主要存在于近尾迹. 蜿蜒是远尾迹的主要特征, 影响下游风力机的来流特征. 尾迹蜿蜒的产生有两种机制: 来流大尺度涡和剪切层失稳. 数值和观测结果显示两种机制共同存在. 机舱和中心涡对尾迹蜿蜒有重要影响. 采用叶片和机舱的致动面模型可准确预测尾迹蜿蜒. 研究显示不同风力机尾迹间的湍流特征存在相似性, 为发展尾迹湍流的快速预测模型提供了理论依据. 当前研究多关注平坦地形上的风力机尾迹, 复杂地形和海洋环境下的大气湍流和风力机尾迹的机理复杂, 现有工程模型无法准确预测, 有待深入研究.      

Dynamic wind loads and wake characteristics of a wind turbine model in an atmospheric boundary layer wind

An experimental study was conducted to characterize the dynamic wind loads and evolution of the unsteady vortex and turbulent flow structures in the near wake of a horizontal axis wind turbine model placed in an atmospheric boundary layer wind tunnel. In addition to measuring dynamic wind loads (i.e., aerodynamic forces and bending moments) acting on the wind turbine model by using a high-sensitive force-moment sensor unit, a high-resolution digital particle image velocimetry (PIV) system was used to achieve flow field measurements to quantify the characteristics of the turbulent vortex flow in the near wake of the wind turbine model. Besides conducting “free-run” PIV measurements to determine the ensemble-averaged statistics of the flow quantities such as mean velocity, Reynolds stress, and turbulence kinetic energy (TKE) distributions in the wake flow, “phase-locked” PIV measurements were also performed to elucidate further details about evolution of the unsteady vortex structures in the wake flow in relation to the position of the rotating turbine blades. The effects of the tip-speed-ratio of the wind turbine model on the dynamic wind loads and wake flow characteristics were quantified in the terms of the variations of the aerodynamic thrust and bending moment coefficients of the wind turbine model, the evolution of the helical tip vortices and the unsteady vortices shedding from the blade roots and turbine nacelle, the deceleration of the incoming airflows after passing the rotation disk of the turbine blades, the TKE and Reynolds stress distributions in the near wake of the wind turbine model. The detailed flow field measurements were correlated with the dynamic wind load measurements to elucidate underlying physics in order to gain further insight into the characteristics of the dynamic wind loads and turbulent vortex flows in the wakes of wind turbines for the optimal design of the wind turbines operating in atmospheric boundary layer winds.     

Wake impacts on aerodynamic and aeroelastic behaviors of a horizontal axis wind turbine blade for sheared and turbulent flow conditions

The magnitude and temporal variations of wind speed considerably influence aerodynamic and structural responses of MW-sized horizontal axis wind turbines. Thus, this paper investigates the variations in airloads and blade behavior of a wind turbine blade resulting from operations in sheared and turbulent flow conditions. First, in order to validate the present methods, comparisons of aerodynamic results were made among the blade element momentum method, free-wake method, and numerical results from the previous studies. Then, the validated methods were applied to a national renewable energy laboratory 5 MW reference wind turbine model for fluid–structure interaction analyses. From the numerical simulations, it can be clearly seen that unfavorable airloads and blade deformations occur due to the sheared and turbulent flow conditions. In addition, it is clear that wake impacts are not as substantial at those of high wind speeds; however, the effects obviously affect the aerodynamic and structural behaviors of the blade at lower wind speeds. Therefore, it is concluded that the numerical results markedly indicate the demand for accurate assessment of wake dynamics for accurate estimations of the aerodynamic and structural responses for sheared and turbulent flow environments.     

Reviews of fluid dynamics researches in wind energy engineering

中国风电行业发展十分迅速, 总装机容量和新增装机容量已达到世界第一, 主力机型将从兆瓦级发展到多兆瓦级(> 5 MW), 风资源已从陆地扩展到海上, 这些变化将带来一系列新的问题. 本文将就目前风力发电行业中的流体力学问题进行综述和探讨, 总结其研究现状及取得的进展. 主要包括大气边界层风特性、局部复杂地形影响及微观选址、风力机尾流干扰、风力机专用翼型族、风力机性能计算、气动弹性、海上风力发电的基础及水力作用等问题. 此外, 还针对目前发展迅速的计算流体力学问题进行了讨论. 最后指出中国风电技术的不足及措施.     

Recent improvements of actuator line-large-eddy simulation method for wind turbine wakes

In a large wind farm, the wakes of upstream and downstream wind turbines can interfere with each other, affecting the overall power output of the wind farm. To further improve the numerical accuracy of the turbine wake dynamics under atmosphere turbulence, this work proposes some improvements to the actuator line-large-eddy simulation(AL-LES) method. Based on the dynamic k-equation large-eddy simulation(LES), this method uses a precursor method to generate atmospheric inflow turbulence,models th...     

Modeling of lifting‐device aerodynamics using the actuator surface concept

The actuator surface (AS) concept and its implementation within a differential, Navier–Stokes control volume finite‐element method (CVFEM) are presented in this article. Inspired by vortex and actuator disk methods, the AS concept consists of using porous surfaces carrying velocity and pressure discontinuities to model the action of lifting surfaces on the flow. The underlying principles and mathematics associated with AS are first reviewed, as well as their implementation in a CVFEM. Results are presented for idealized 2D cases with analytical solutions, as well as for the 3D cases of a finite wing and an experimental wind turbine. In the case of the finite wing, wake induction is well handled by the model with accurate predictions of induced angles and drag when compared with the Prandtl lifting line model. Comparisons with volume force approaches, often used to model the action of propellers or wind turbine blades in a simplified analysis, show that the AS concept has some interesting advantages in terms of accuracy and respect of flow physics. This new approach is easy and rapid to embed in most computational fluid dynamics (CFD) methods. It is applicable to a wide range of problems involving thin lifting devices like finite wings, propellers, helicopter or wind turbine blades. Copyright © 2009 John Wiley & Sons, Ltd.