Development of an interface to introduce stationary LES data to the URANS solver THETA for HAWT performance prediction

Abstract

The impact of different sheared velocity profiles on the performance prediction of a horizontal axis wind turbine in the atmospheric boundary layer is investigated. Firstly, the wall roughness in the analytical logarithmic description of the atmospheric boundary layer is varied to obtain different velocity profiles. Subsequently, it is proposed to replace the analytical logarithmic description of the atmospheric boundary layer by the time‐averaged velocity data of a precursor large eddy simulation (LES) and to reconstruct the turbulence of the velocity fluctuations. The LES data are introduced as inflow condition through a LES‐RANS interface in a one‐way coupling approach. Three different methods to reconstruct URANS turbulence values out of the velocity fluctuations are investigated. It is shown that the reconstruction method has an impact on the development of the velocity profile, turbulent kinetic energy, and the turbulent dissipation during the transport through the URANS domain. The different inflow data, which the horizontal axis wind turbine experiences, are responsible for changes in the overall rotor thrust (up to 2.7%) and rotor torque (up to 2.4%). Conversely, the induction factors and effective angles of attack hardly change and can well be compared with a blade element momentum method. Finally, the results of both approaches to prescribe the atmospheric boundary layer are compared. The thrust and power coefficients, and wake recovery are close to each other. Simulations are carried out on an industrial 900 kW wind turbine with the incompressible URANS solver THETA.