An analytical linear two-dimensional actuator disc model and comparisons with computational fluid dynamics (CFD) simulations

• Published in: Wind Energy Science Vol. 8; no. 12; pp. 1853 - 1872
• Main Author: Madsen, Helge Aagaard
• Format: Journal Article
• Language: English
• Published: Göttingen Copernicus GmbH 13.12.2023; Copernicus Publications
• Subjects: Approximation; Fluid dynamics; Hydrodynamics; Simulation; Turbines; Velocity; Velocity distribution; Vortices; Wind power; Wind shear
• DOI: 10.5194/wes-8-1853-2023
• ISSN: 2366-7443

The continuous up-scaling of wind turbines enabled by more lightweight and flexible blades in combination with coning has challenged the assumptions of a plane disc in the commonly used blade element momentum (BEM)-type aerodynamic codes for the design and analysis of wind turbines. The objective with the present work is thus to take a step back relative to the integral 1-dimensional (1-D) momentum theory solution in the BEM model in order to study the actuator disc (AD) flow in more detail.We present an analytical, linear solution for a two-dimensional (2-D) AD flow with one equation for the axial velocity and one for the lateral velocity, respectively. Although it is a 2-D model, we show in the paper that there is a good correlation with axis-symmetric and three-dimensional (3-D) computational fluid dynamics (CFD) simulations on a circular disc. The 2-D model has thus the potential to form the basis for a simple and consistent rotor induction model.For a constant loading, the axial velocity distribution at the disc is uniform as in the case of the classical momentum theory for an AD. However, an important observation of the simulated flow field is that immediately downstream of the disc the axial velocity profiles change rapidly to a shape with increased induction towards the edges of the disc and less induction on the central part. This is typically what is seen at the disc in full non-linear CFD AD simulations, which is what we compare with in the paper.By a simple coordinate rotation the analytical solution is extended to a yawed disc with constant loading. Again, a comparison with CFD, now with a 3-D simulation on a circular disc in yaw, confirms a good performance of the analytical 2-D model for this more complicated flow.Finally, a further extension of the model to simulate a coned disc is obtained using a simple superposition of the solution of two yawed discs with opposite yaw angles and positioned so the two discs just touch each other. Now the validation of the model is performed with results from axis-symmetric CFD simulations of an AD with a coning of both 20 and -20∘. In particular, for the disc coned in the downwind direction there is a very good correlation between the simulated normal velocity to the disc, whereas some deviations are seen for the upwind coning.The promising correlation of the results for the 2-D model in comparison with 3-D simulations of a circular disc with CFD for complicated inflow like what occurs at yaw and coning indicates that the 2-D model could form the basis for a new, consistent rotor induction model. The model should be applied along diagonal lines on a rotor and coupled to an angular momentum model. This application is sketched in the outlook and is a subject for future research.