基于Theodorsen理论的新型动态失速预测模型及其应用

风力机叶片动态失速时的非定常气动特性及严重的迟滞现象使得风力机功率实测值严重偏离其静态预测值。鉴于此,基于Theodorsen理论、基尔霍夫势流理论,在忽略低阶附加质量引起的下洗气流加速度项及状态变量转换后,提出一种包括翼型附着流和后缘动态分离流的新型动态失速模型。利用该模型分析NREL 5 MW海上风力机叶片6种翼型的非定常动态失速特性得出：通过翼型的气流在完全附着流与完全分离流之间不断转换,受附着流脱落尾诱导的动态下洗气流影响及边界层动态分离产生的压力滞后的双重作用,动态升力系数变化曲线和静态升力现象曲线偏差较大,6种翼型动态升力系数变化曲线均呈非常明显的迟滞环现象。DU40、DU35、DU30、DU25、DU21和NACA64这6种翼型动态升力系数增幅明显,分别达17.6%、60.9%、60.7%、55.1%、63.7%和40.8%。动态失速攻角极大地超过静态失速攻角,分别增大到36.53°、21.40°、20.20°、17.68°、16.97°和21.42°。6种翼型动态失速预测结果与公开实验数据结论一致,证实所提出的动态失速气动模型计算结果准确可信,具有较强通用性。     

Computational fluid dynamics simulation of the aerodynamics of a high solidity,small‐scale vertical axis wind turbine

A computational fluid dynamics simulation was performed for a small‐scale, high solidity (σ = 0.48) H‐type Darrieus vertical axis wind turbine. Two‐dimensional unsteady Reynolds‐averaged Navier–Stokes equations were solved for the turbine numerical model, which has a large stationary domain and smaller rotating subdomain connected by a sliding mesh interface. The simulation results were first validated against steady‐state airfoil data. The model was then used to solve for three rotating blades with constant ambient flow velocity (Re = 360,000) over numerous blade speed ratios. The high solidity and the associated low blade speed ratio and rotational speed of the turbine result in complex flow–blade interaction mechanisms. These include dynamic stall resulting in vortex shedding, vortex impingement on the source blade and significant flow momentum extraction causing reduced power production from the downstream blade pass. Copyright © 2011 John Wiley & Sons, Ltd.     

定常吸气对垂直轴风力机性能影响研究

定常吸气装置可有效提高垂直轴风力机气动性能,改善风轮流场结构及翼型动态失速特性。基于CFD方法对垂直轴风力机进行数值模拟,研究不同叶尖速比（TSR）下定常吸气对风力机气动及流场特性的影响,对比分析原始风力机及定常吸气作用下的风能利用率、整机转矩系数及涡量分布。结果表明：不同尖速比下定常吸气均可显著提高风力机气动性能,减小风轮载荷波动,降低最佳叶尖速比,提高风力机运行稳定性;叶尖速比为2.51时,风能利用系数增加34.69%;定常吸气削弱了风轮叶片间尾涡脱落的影响,抑制叶片前缘涡的形成,减缓了叶片的动态失速现象,对风轮流场有良好的改善效果。     

Simulating variable pitch crossflow water turbines: A coupled unsteady ONERA-EDLIN model and streamtube model

This article describes a new method for simulating unsteady hydrodynamics forces and moments on the blades of a crossflow ‘Darrieus’ turbine with active pitch variation. This method is based on the ONERA-EDLIN dynamic stall model, coupled with a momentum streamtube model to take into account the turbine interference on the flow. Both models are presented, and compared separately with experimental results for a pitching airfoil for the ONERA-EDLIN model; and for Darrieus turbine for the momentum theory. The model coupling is then detailed and compared with experimental data taken from the open literature [1] The turbine has 2 straight blades with a NACA 0012 section operating in water at a mean chord Reynolds number of 4 × 10⁴ for tip speed ratio λ = 2.5, 5 and 7.5. Good agreement was found for average λ = 5, and qualitative agreement could be obtained at low and high λ, where dynamic stall effects and interference effects respectively are predominant. This is positive because λ = 5 is the closest value from the optimal power production point. Variable pitch is finally introduced in the model and several functions are tested in order to increase efficiency. A maximum increase of 53% on the power coefficient was found to occur with a sinusoidal law.     

Wind tunnel quantification of dynamic stall on an S817 airfoil and its control using synthetic jet actuators

Wind tunnel experiments were performed to quantify the aerodynamic characteristics of the S817 airfoil in dynamic stall conditions, and the subsequent application of active flow control to modify the manner by which dynamic stall incepts. Both quasi‐2D and cantilevered finite span configurations were tested. Surface pressure, six‐component force‐torque sensor, and stereoscopic particle image velocimetry (SPIV) were used to quantify the baseline flow and the benefits of actuating synthetic jets (installed at x/c = 0.35, angled 45° into the flow, and at a momentum coefficient C_μ = 0.012). The airfoil was pitched at reduced frequencies of k_f = 0.025 and 0.05 and at shallow and deep stall. Vortex induced lift from dynamic stall was observed and was eliminated by the use of synthetic jets for nearly all conditions; pitching moment deviation was also observed to be significant, and was eliminated at shallow stall and significantly reduced during deep dynamic stall when the synthetic jets were actuated. Moreover, the activation of synthetic jets resulted in significant reduction in the hysteresis (area within the pitching up and pitching down load history) of the lift and pitching moment through all experimental conditions, as much as 41% and 85%, respectively. SPIV flow fields in shallow dynamic stall demonstrated that actuation of synthetic jets confined the separated region to the trailing edge, in both the instantaneous and time averaged sense. To further reduce the lift and pitching moment hysteresis at high angles of attack, a pulse modulation technique was used and showed a marked increase in synthetic jet performance compared with the continuously actuated case and achieved this result with approximately 65% less power consumption.     

Effects of nonuniform incident velocity on a dynamic wind turbine airfoil

For further insight into the performance of a horizontal axis wind turbine blade section under yaw loads, a 2D numerical simulation of a pitching S809 airfoil under dynamic stall with an unsteady incident velocity is presented. The streamwise incident velocity and pitch angle of incidence oscillated with the same frequency but with a range of phase differences, ? π ≤ Φ ≤ π. Changing Φ caused variation of the results, which can be highlighted as significantly augmented and dramatically damped dynamic stall loads, both increasing and decreasing trends for vortex growth time during Φ increase, a shifted location of the maximum loads and a change in the order of the vortex pair circulation in each cycle. The results showed strong dependency on the velocity and acceleration of the freestream during dynamic stall, which categorized the results in four individual subdomains with different behaviors. Copyright © 2014 John Wiley & Sons, Ltd.     

吸气控制策略对垂直轴风力机气动性能影响研究

针对垂直轴风力机复杂气动特性,将吸气孔置于风力机翼型上下表面,提出不同吸气控制策略以改善其气动性能。基于CFD方法,研究不同叶尖速比下吸气策略对风力机风能利用率、叶片切向力系数及流场特性的影响,综合考虑能量消耗与风力机输出功率。结果表明：提出的3种控制策略在低叶尖速比下均能大幅提升整机气动效率。效果最佳的迎、背风区交替吸气策略可显著推迟分离点,延缓翼型动态失速发生,并减少分离涡周期性脱落造成的损失。此外,该策略对动态尾迹效应有良好的控制效果,同时降低整机转矩波动幅值,消除中低叶尖速比下风轮负转矩,从而提高获能效率且延长风力机使用寿命。     

A vortex panel model for the simulation of the wake flow past a vertical axis wind turbine in dynamic stall

A 2D vortex panel model with a viscous boundary layer formulation has been developed for the numerical simulation of a vertical axis wind turbine (VAWT), including the operation in dynamic stall. The model uses the ‘double wake’ concept to reproduce the main features of the unsteady separated flow, including the formation and shedding of strong vortical structures and the wake–blade interaction. The potential flow equations are solved together with the integral boundary layer equations by using a semi‐inverse iterative algorithm. A new criterion for the reattachment of the boundary layer during the downstroke of a dynamically stalled aerofoil is implemented. The model has been validated against experimental data of steady aerofoils and pitching aerofoils in dynamic stall at high and low Reynolds numbers (Re = 1.5 × 10⁶ and Re = 5 × 10⁴). For the low Reynolds number case, time‐resolved 2D particle image velocimetry (PIV) measurements have been performed on a pitching NACA 0012 aerofoil in dynamic stall. The PIV vorticity fields past the oscillating aerofoil are used to test the model capability of capturing the formation, growth and release of the strong leading edge vortex that characterizes the dynamic stall. Furthermore, the forces extracted from the PIV velocity fields are compared with the predicted ones for a quantitative validation of the model. Finally, the model is applied to the computation of the wake flow past a VAWT in dynamic stall; the predicted vorticity fields and forces are in good agreement with phase‐locked PIV data and CFD‐DES available in the literature. Copyright © 2012 John Wiley & Sons, Ltd.     

Dynamic stall control via adaptive blowing

An aerodynamic load control concept termed “adaptive blowing” was successfully tested on a NACA 0018 airfoil model at Reynolds numbers ranging from 1.5·10⁵ to 5·10⁵. The global objective was to eliminate lift oscillations typically encountered on wind turbine blade sections. Depending on the jet momentum flux, steady blowing from a control slot in the leading-edge region can be utilized to either enhance or reduce lift by suppressing or inducing boundary layer separation respectively. Furthermore, high momentum blowing effectively eliminated the dynamic stall vortex during deep dynamic stall conditions. Based on these previous findings, the present work explores the feasibility of controlling unsteady aerodynamic loads by dynamically varying the jet momentum flux to compensate for transient changes of the inflow. Various scenarios including high amplitude pitching, rapid freestream oscillations and combinations of both were investigated in a custom-built unsteady wind tunnel facility. An iterative control algorithm was implemented which successfully identified the momentum coefficient time profiles required to minimize the lift excursions. The combination of fully suppressing dynamic stall and dynamically adjusting the lift coefficient provided an unprecedented control authority, producing virtually constant phase averaged lift in all cases.     

凹槽-襟翼对翼型动态失速特性影响研究

风力机复杂运行环境使叶片常处于失速环境,导致翼型升力骤降,严重影响风力机气动性能.为改善翼型流动分离,延缓失速,对凹槽-襟翼对翼型动态失速特性作用效果开展研究,并利用计算流体力学方法分析不同折合频率与翼型厚度时凹槽-襟翼对翼型气动性能的影响.结果表明:俯仰振荡过程中,凹槽-襟翼可有效提升翼型吸力面流速,降低失速攻角下逆...     

Effect of acceleration on dynamic stall of airfoil in unsteady operating conditions

In this paper the effect of accelerated flow over a moving airfoil is considered and based on the flow field around the airfoil the dynamic stall is evaluated. In contrast to ordinary pitching motion, the dynamic stall evaluation in this study is performed with a different motion pattern, in which the airfoil has a heaving motion in one direction. This motion pattern is also similar to rotation of an element of blade in horizontal axis wind turbines (HAWTs). In present investigation, the Reynolds number is changed during simulation time and variations of this parameter from initial to final values are shown by acceleration parameter. The operating Reynolds number is more than 10⁶, and a S809 airfoil is selected to move with accelerations of 1, 4 and 6 m/s² in normal direction to free stream. To resolve accelerated flow filed in the two‐dimensional computational domain and to achieve results within a reasonable computation time, the unsteady Reynolds‐Averaged Navier–Stokes (URANS) equations are employed. The governing equations are discretized based on the finite volume approach and semi‐implicit method for pressure linked equations algorithm is used for pressure–velocity coupling. Furthermore, turbulence effect on flow field is accounted using shear stress transport (SST) k‐ω model. It is shown that the accelerated flow can significantly influence on the aerodynamic loads and dynamic stall trend. This study may introduce a new concept regarding dynamic stall and aerodynamic loads when the rotational acceleration is involved in HAWTs. Copyright © 2014 John Wiley & Sons, Ltd.     

Validation of the Beddoes–Leishman dynamic stall model for horizontal axis wind turbines using MEXICO data

The aim of this study is to assess the load predicting capability of a classical Beddoes–Leishman dynamic stall model in a horizontal axis wind turbine environment, in the presence of yaw misalignment. The dynamic stall model was tailored to the horizontal axis wind turbine environment and validated against unsteady thick airfoil data. Subsequently, the dynamic stall model was implemented in a blade element‐momentum code for yawed flow, and the results were compared with aerodynamic measurements obtained in the MEXICO (Model Rotor Experiments under Controlled Conditions) project on a wind turbine rotor placed in a large scale wind tunnel. In general, reasonable to good agreement was found between the blade element‐momentum model and MEXICO data. When large yaw misalignments were imposed, poor agreement was found in the downstroke of the movement between the model and the experiment. Still, over a revolution, the maximum normal force coefficient predicted was always within 8% of experimental data at the inboard stations, which is encouraging especially when blade fatigue calculations are being considered. Copyright © 2012 John Wiley & Sons, Ltd.     

Analysis of vortex‐induced and stall‐induced vibrations at standstill conditions using a free wake aerodynamic code

Vortex‐induced and stall‐induced vibrations of a 2D elastically mounted airfoil at high angles of attack in the vicinity of 90° are investigated using a vortex type model. Such conditions are encountered in parked or idling operation at extreme yaw angles provoked by control system failures. At very high angles of attack, massive flow separation takes place over the entire blade span, and vortex shedding evolves downstream of the blade giving rise to periodically varying loads at frequencies corresponding to the Strouhal number of the vortices shed in the wake. As a result, vortex‐induced vibrations may occur when the shedding frequency matches the natural frequency of the blade. A vortex type model formulated on the basis of the ‘double wake’ concept is employed for the modelling of the stalled flow past a 2D airfoil. By tuning the core size of the vortex particles in the wake, the model predictions are successfully validated against averaged 2D measurements on a DU‐96‐W‐180 airfoil at high angles of attack. In order to assess the energy fed to the airfoil by the aerodynamic loads, the behaviour under imposed sinusoidal edgewise motions is analysed for various oscillation frequencies and amplitudes. Moreover, stall‐induced and vortex‐induced vibrations of an elastically mounted airfoil section are assessed. The vortex model predicts higher aeroelastic damping as compared with that obtained using steady‐state aerodynamics. Excessive combined vortex‐induced and stall‐induced edgewise vibrations are obtained beyond the wind speed of 30 m s^?1. Copyright © 2014 John Wiley & Sons, Ltd.     

潮流能水轮机叶片翼型水动力性能的数值模拟研究

鉴于潮流能水轮机叶轮叶片性能的好坏关系到叶轮的捕能效果,而叶片翼型的水动力性能是叶片设计的基础,基于CFD软件ANSYS CFX对部分翼型的水动力性能进行了分析。首先在不同湍流模型下对NACA2412翼型的部分水动力特性进行数值模拟并与试验值做对比分析;然后在SST湍流模型下,对NACA64421翼型在攻角为0°时,探讨了边界层网格对数值模拟结果的影响;并对其在攻角为-2°~20°情形下的水动力性能进行了数值模拟。结果表明,湍流模型的选用、边界层网格划分及近壁尺寸的选取对数值模拟结果的影响均较大,并通过数值模拟得到了NACA64421翼型的升力系数、阻力系数、最佳攻角、失速角等水动力性能参数,可为水轮机叶片设计时翼型的选取提供参考。     

尾缘非定常射流襟翼对垂直轴风力机气动特性影响研究

为提升垂直轴风力机气动性能并改善其动态失速特性,将射流襟翼布置于翼型尾缘压力面,并提出5种射流控制策略,采用计算流体力学方法研究不同策略对垂直轴风力机气动性能影响,从而确定最佳控制策略。结果表明：在180°~360°相位角范围内施加射流控制可使风力机风能利用系数在最佳尖速比下提升31.31%,并有效抑制吸力面尾缘涡形成与发展,增大翼面两侧压差;射流越靠近尾缘,垂直轴风力机气动性能提升效果越好。     

Airfoil optimization for wind turbine application

An airfoil optimization method for wind turbine applications that controls the loss in performance due to leading edge contamination is developed and tested. The method uses the class‐shape‐transformation technique to parametrize the airfoil geometry and uses an adjusted version of the panel code XFOIL to calculate the aerodynamic performance. To find optimal airfoil shapes, the derivative‐free Covariance Matrix Adaptation Evolution Strategy is used in combination with an adaptive penalty function. The method is tested for the design of airfoils for the outer part of a megawatt‐class wind turbine rotor blade, and the results are compared with airfoils from Delft University. It is found that the method is able to automatically create airfoils with equal or improved performance compared with the Delft designs. For the tested application, the adjustments performed to the XFOIL code improve the maximum lift, post stall, and the overall drag predictions.     

Vortex‐induced vibrations on a modern wind turbine blade

This article investigates the aero‐elastic response of the DTU 10‐MW RWT blade in deep stall conditions with angles of attack in the vicinity of 90 degrees. The simulations were conducted with the high‐fidelity fluid–structure interaction simulation tool HAWC2CFD employing the multi‐body‐based structural model of HAWC2 and the incompressible computational fluid dynamics solver EllipSys3D. The study utilizes detached eddy simulation computations and considers the three‐dimensional blade geometry including blade twist and taper. A preliminary frequency analysis of the load variations on a stiff blade showed that an inclined inflow with a velocity component along the blade axis can trigger a spanwise correlated vortex shedding over large parts of the blade. Moderate wind speeds were sufficient to generate vortex shedding with frequencies close to the first edgewise eigenfrequency of the blade. Aero‐elastic computations of the elastic blade confirmed the findings of the frequency analysis. Inflow conditions with inclination angles between Ψ = 20° and Ψ = 55° and relatively low to moderate wind speeds between V = 16 and V = 26ms^?1 were sufficient to trigger severe edgewise blade vibrations with blade tip amplitudes of several metres. The investigated inflow conditions are considered realistic and might occur when the wind turbine is idling or standing still and the yaw system is unable to align the wind turbine with the incoming wind. Copyright © 2016 John Wiley & Sons, Ltd.     

考虑转捩的风力机翼型动态失速数值模拟

以风力机专用翼型的动态失速为对象,采用一种基于流场当地变量的Gamma-Theta转捩模型配合SSTk-ω湍流模型进行数值模拟,研究转捩对动态失速性能的影响和动态失速下的转捩规律。结果表明,使用考虑转捩效应,能够使动态失速过程中上仰段大迎角状态下失速和下俯段气流再附的模拟得到改善。在动态失速上仰段,上表面转捩由后缘分离泡向前缘分离泡的转变过程较快,导致转捩点迅速前移;而在下俯段,前缘分离泡向后缘分离泡的转变过程中经过了自然转捩和再层流化的过渡,因此转捩点的移动较上仰段平滑。     

Spectral analysis of New MEXICO standstill measurements to investigate vortex shedding in deep stall

Spectral analysis was performed on the time series data computed from pressure measurements on the New MEXICO (Model Rotor Experiments under Controlled Conditions) rotor in standstill conditions. As a priori, 3D airfoil polars were recreated from standstill measurements and compared against 2D airfoil polars and flat plate theory results to verify the measurements. The spectral analysis revealed the presence of dominant shedding frequencies for certain ranges of the geometric angle of attack. Two dominant shedding modes were identified: One was associated with bluff body vortex shedding, and the other was associated with low Strouhal number shedding. No dominant shedding frequencies were observed for angles of attack beyond 50°. The research improves on our current understanding of the unsteady nature of the stall regime, along with providing insight into the existence of vortex‐induced vibrations on a wind turbine in standstill conditions.     

On the S809 airfoil's unsteady aerodynamic characteristics

The S809 airfoil dynamic characteristics, which are based on the airfoil dynamic tests at the University of Glasgow, are presented in this paper. The airfoil tests include static, ramp‐type (ramp‐up and ramp‐down) and oscillatory motions at Reynolds numbers of 1.0 × 10⁶ and 1.5 × 10⁶ with and without the sand‐tripped leading edge. This study aimed to explore the unsteady aerodynamic features of the S809 airfoil, such as the progression of separation from leading edge to trailing edge, the large trailing edge separation before stalling, the stall onset inception and the re‐attachment convection, and to provide some useful data for tuning/refining the semi‐empirical dynamic stall (DS) models, such as the Leishman–Beddoes DS model or its variations. Experimental results show that the S809 airfoil has a complicated DS process that renders this airfoil a challenge for any modeller of the unsteady airloads. The leading‐edge roughness has small effect on the static features, but significantly invokes earlier stall onset inception under dynamic conditions, while has small influence on the convective phase of the re‐establishment of fully attached flow. Copyright © 2009 John Wiley & Sons, Ltd.