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

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

复杂工况下的大型风力机气动性能和尾迹研究Investigation of aerodynamic performance and wake of the large scale wind turbine in complicated conditions

在复杂工况下,大型风力机非定常特性会更严重,导致风力机气动性能变化和尾迹预测更加复杂。本文主要针对稳态偏航、动态偏航、风剪切和随机风速场等复杂工况,基于自由涡尾迹方法,嵌入复杂工况的模块,加入了动态失速模型和三维旋转效应模型修正,实现了复杂工况数值模拟计算,比较了不同复杂工况的气动载荷和尾迹形状。最后,得出了风力机在复杂工况下的气动性能、载荷和尾迹叶尖涡线特性,并计算出风力机在复杂工况下的气动载荷超调量和迟滞时间。对推进自由涡尾迹方法应用于风力机工程的大批工况载荷计算,提高大型风力机的载荷计算精度和设计水平等具有重要意义。     

Reduced order model of offshore wind turbine wake by proper orthogonal decomposition

A wide range of previously designed methods for faster parametrization of partial differential equations requires them to be solved using existing finite volume, finite element, and finite difference solvers. Due to the requirement of high degrees of freedom to accurately model the physical system, computational costs often becomes a bottle-neck. It poses challenges to conducting efficient repeated parametric sampling of the input parameter that disrupts the whole design process. Model reduction techniques adopted to high fidelity systems provide a basis to accurately represent a physical system with a lower degree of freedom. The present work focuses on one such method for high-fidelity simulations that combines finite volume strategy with proper orthogonal decomposition and Galerkin projection to test reduced-order models for high Reynolds number flow applications. The model is first benchmarked against flow around a cylinder for which extensive numerical and experimental data is available in the literature. The models are then tested to full-scale NREL 5MW offshore wind turbines to evaluate wake evolution in the downstream direction. The simulations results show relative errors of wind turbines for the first seventy modes approach 4.7% in L²-norm for velocities.     

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.     

大型海上风力机尾迹区域风场分析

随着世界范围内海上风电场的不断兴起及海上风电场设计规模的日趋庞大,大型海上风力机尾迹区域风场特征的研究对于海上风电场的优化布局有着重要的指导意义。本文基于三维Navier-Stocks控制方程和适用于旋转流场分析的RNGk-ε湍流模型,采用滑移网格技术对美国可再生能源实验室(NREL)的5 MW海上风力机的性能及其尾迹区域的风场特征进行了较为系统的数值模拟。通过将不同风速下风力机输出功率的数值结果与NREL的设计参考数据进行对比,三维数值模型的有效性得到了很好地验证。此外,在此基础上进一步研究了大型海上风力机额定风速下及不同风轮转速下尾迹区域平均风速的分布特征,并得到了一系列具有参考价值的重要结论。     

Research status and trend of wind turbine aerodynamic noise

The main components of the wind turbine aerodynamic noise are introduced. A detailed review is given on the theoretical prediction, experimental measurement, and numerical simulation methods of wind turbine noise, with specific attention to applications. Furthermore, suppression techniques of wind turbine aerodynamic noise are discussed. The perspective of future research on the wind turbine aerodynamic noise is presented.     

Determination of real-time predictors of the wind turbine wake meandering

The present work proposes an experimental methodology to characterize the unsteady properties of a wind turbine wake, called meandering, and particularly its ability to follow the large-scale motions induced by large turbulent eddies contained in the approach flow. The measurements were made in an atmospheric boundary layer wind tunnel. The wind turbine model is based on the actuator disc concept. One part of the work has been dedicated to the development of a methodology for horizontal wake tracking by mean of a transverse hot wire rake, whose dynamic response is adequate for spectral analysis. Spectral coherence analysis shows that the horizontal position of the wake correlates well with the upstream transverse velocity, especially for wavelength larger than three times the diameter of the disc but less so for smaller scales. Therefore, it is concluded that the wake is actually a rather passive tracer of the large surrounding turbulent structures. The influence of the rotor size and downstream distance on the wake meandering is studied. The fluctuations of the lateral force and the yawing torque affecting the wind turbine model are also measured and correlated with the wake meandering. Two approach flow configurations are then tested: an undisturbed incoming flow (modelled atmospheric boundary layer) and a disturbed incoming flow, with a wind turbine model located upstream. Results showed that the meandering process is amplified by the presence of the upstream wake. It is shown that the coherence between the lateral force fluctuations and the horizontal wake position is significant up to length scales larger than twice the wind turbine model diameter. This leads to the conclusion that the lateral force is a better candidate than the upstream transverse velocity to predict in real time the meandering process, for either undisturbed (wake free) or disturbed incoming atmospheric flows.     

The stability and development of tip and root vortices behind a model wind turbine

An inclined turbulent jet discharging a passive scalar into a turbulent crossflow is investigated under conditions of favorable, zero and adverse streamwise pressure gradient. Experiments are conducted in water by means of magnetic resonance velocimetry and magnetic resonance concentration measurements. The velocity ratio and density ratio are equal to one for all cases. The flow configuration is relevant to film cooling technology, the molecular properties of the fluid being immaterial in the fully turbulent regime. Under favorable pressure gradient (FPG), the streamwise acceleration tilts the jet trajectory toward the wall, which would be beneficial for the film cooling performance. However, the counter-rotating vortex pair is strengthened in the accelerating flow by streamwise stretching. Also, the crossflow boundary layer is significantly thickened by increasingly adverse pressure gradient, which affects the mass transfer from the jet. Overall, the more intense counter-rotating vortices and the thinner boundary layer associated with increasingly FPG enhance the scalar dispersion into the main flow, hampering the film cooling performance in terms of surface effectiveness.     

The wake flow control behind a circular cylinder using ion wind

Active and passive flow control methods have been studied for decades, but there have been only a few studies of flow control methods using ion wind, which is the bulk motion of neutral molecules driven by locally ionized air from a corona discharge. This paper describes an experimental study of ion wind wake control behind a circular cylinder. The experimental conditions consisted of a range of electrohydrodynamic numbers—the ratio of an electrical body force to a fluid inertial force—from 0 to 2 and a range of Reynolds numbers from 4×10³ to 8×10³. Pressure distributions over the cylinder surface were measured and flow visualizations were carried out using a smoke-wire method. The flow visualizations confirmed that ion wind significantly affects the wake structure behind a circular cylinder, and that the pressure drag can be dramatically reduced by superimposing ion wind.List of symbols BR blockage ratio - C _d coefficient of the pressure drag - C _p coefficient of the surface pressure, 2(p–p ₀)/(U ₀ ²) - C _pb coefficient of the base surface pressure, 2(p _b–p ₀)/(U ₀ ²) - D diameter of the cylinder - D _P pressure drag - d _p diameter of particle - E the electric field - F _e Coulombian force (qE) - F _v viscous force - H wire-to-cylinder spacing - I total electric current (A) - L the axial length of cylinder (m) - N _EHD electrohydrodynamic number - p _b base pressure of cylinder at =180° - p ₀ reference static pressure at 10D upstream - q the charge on the particle - R radius of the cylinder - V applied voltage (kV) - U ₀ mean flow velocity (m/s) - ion mobility in air (m²/(s V)) - ₀ permittivity of free space - viscosity of fluid (kg/ms) - density of fluid (kg/m³) - installation angle of a wire electrode (°)     

Dynamic stall model for wind turbine airfoils

A model is presented for aerodynamic lift of wind turbine profiles under dynamic stall. The model combines memory delay effects under attached flow with reduced lift due to flow separation under dynamic stall conditions. The model is based on a backbone curve in the form of the static lift as a function of the angle of attack. The static lift is described by two parameters, the lift at fully attached flow and the degree of attachment. A relationship between these parameters and the static lift is available from a thin plate approximation. Assuming the parameters to be known during static conditions, nonstationary effects are included by three mechanisms: a delay of the lift coefficient of fully attached flow via a second-order filter, a delay of the development of separation represented via a first-order filter, and a lift contribution due to leading edge separation also represented via a first-order filter. The latter is likely to occur during active pitch control of vibrations. It is shown that all included effects can be important when considering wind turbine blades. The proposed model is validated against test data from two load cases, one at fully attached flow conditions and one during dynamic stall conditions. The proposed model is compared with five other dynamic stall models including, among others, the Beddoes–Leishman model and the ONERA model. It is demonstrated that the proposed model performs equally well or even better than more complicated models and that the included nonstationary effects are essential for obtaining satisfactory results. Finally, the influence of camber and thickness distribution on the backbone curve are analysed. It is shown that both of these effects are adequately accounted for via the static input data.     

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

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

Investigation of horizontal-axis wind turbine (HAWT) blade threedimensional rotational effect based on field experiments

Field experiments are performed on a two-bladed 33 kW horizontal-axis wind turbine (HAWT). The pressures are measured with 191 pressure sensors positioned around the surfaces of seven spanwise section airfoils on one of the two blades. Three-dimensional (3D) and two-dimensional (2D) numerical simulations are performed, respectively, on the rotor and the seven airfoils of the blade. The results are compared with the experimental results of the pressure distribution on the seven airfoils and the lift coefficients. The 3D rotational effect on the blade aerodynamic characteristics is then studied with a numerical approach. Finally, some conclusions are drawn as follows. From the tip to the root of the blade, the experimental differential pressure of the blade section airfoil increases at first and then decreases gradually. The calculated 3D result of the pressure distribution on the blade surface is closer to that of the experiment than the 2D result. The 3D rotational effect has a significant impact on the blade surface flow and the aerodynamic load, leading to an increase of the differential pressure on the airfoils and their lift coefficient than that with the 2D one because of the stall delay. The influence of the 3D rotational effect on the wind turbine blade especially takes place on the sections with flow separation.     

Motion of triple configuration of shock waves with formation of wake behind branching point

If three developing shock-wave fronts come together at one point the laws of conservation connecting the parameters of the gas in the vicinity of this point give an overdetermined system of equations. To remove the possible contradiction it is necessary to increase the number of initial parameters. As a rule, the assumption of the presence of a contact discontinuity emerging from the branching point is sufficient. It is also possible for two contact discontinuities to develop which form two shock waves with respect to the branching point opposite the boundary of the isobaric region filled with gas in a state of rest. Such a region is called the wake of the triple point by analogy with the aerodynamic wake for flow around bodies with flow separation. A closed system of five simultaneous differential equations which describes approximately the dynamics of a stream containing a branched system of shock waves with a developing wake behind the triple point is derived and discussed in the report.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 67–75, November–December, 1973.     

基于涡尾迹方法的风力机气动性能分析

基于(势流)涡尾迹方法开发了水平轴风力机叶片气动性能分析程序,采用固定尾迹涡模型和自由尾迹模型分别对气动设计性能进行计算分析,得到风力机设计工况下的涡位置、诱导因子、功率系数及扭矩系数等气动性能参数,并与设计结果对比。结果表明,涡尾迹方法能够快速准确地计算风力机叶片气动性能参数,对风力机叶片气动分析,固定尾迹涡模型较自由尾迹模型计算时间短,具有较好的实用性。     

A general model for trailing edge serrations simulation on wind turbine airfoils

Trailing edge serrations(TESs) are capable of noticeably suppressing the turbulent trailing edge noise induced by rotating wind turbine blades and become an integral part of a blade. However, the challenges involved in the dimensional design of serration height 2 h, wavelength λ and flap angle are Φ yet to be dealt with in a satisfactory manner. To address the problem, a general model for simulating the effects of serrations on the hydrodynamic and aeroacoustic performance is proposed due to its ease of use and relatively low requirements for user input. The solid serrations are replicated by momentum sources calculated by its aerodynamic forces. Then, a case relevant to wind turbine airfoil is examined, a hybrid improved delay detached eddy simulation(IDDES) method coupled with FW-H integration is deployed to obtain the flow features and far-field sound pressure level. It is found that the modeling method could reproduce the flow field and noise as serrated airfoil.     

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

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

A minimum-domain impulse theory for unsteady aerodynamic force with discrete wake

We extend the impulse theory for unsteady aerodynamics, from its classic global form to finite-domain formulation, then to a minimum-domain version for discrete wake. Each extension has been confirmed numerically. The minimum-domain theory indicates that the numerical finding of Li and Lu(2012) is of general significance: The entire force is completely determined by only the time rate of impulse of those vortical structures still connecting to the body, along with the Lamb-vector integral thereof that captures the contribution of all the rest disconnected vortical structures.     

“Blind Test 3” calculations of the performance and wake development behind two in-line and offset model wind turbines

This is a report on the comparison between measurement data and predictions presented at the “Blind Test 3” Workshop organized jointly by Nowitech and Norcowe in Bergen, 10 and 11 December, 2013. A number of researchers were invited to predict the performances and the wake development behind two model wind turbines that have been extensively tested at the Department of Energy and Process Engineering, NTNU. The turbines were arranged in-line, but slightly offset so that the wake of the upstream turbine only interacted with roughly half the area swept by the second rotor. This is a common event in most wind parks and produces flow fields that are both complicated and harmful for the downstream turbine. As expected it turned out to be a difficult flow to predict.Contributions were received from five different groups using a range of methods, from fully resolved Reynolds averaged Computational Fluid Dynamics (CFD) models to Large Eddy Simulations (LES). The range of results was large but the overall trend is that the current methods predict the power generation as well as the thrust force reasonably well. But there is a large uncertainty in the prediction of the turbulence field in the wake. Hence, the LES method consistently performed better than the others.