Wake modeling and simulation of an experimental wind turbine using large eddy simulation coupled with immersed boundary method alongside a dynamic adaptive mesh refinement

Abstract

Accurate numerical models of the flow around wind turbines is highly relevant for optimizing the energy efficiency of wind farms. The main goal of this work is to apply large eddy simulation (LES) along with an immersed boundary (IB) method to generate spatial and temporal information of the velocity field and forces acting on a fully-resolved NREL S809 airfoil and NREL Phase VI wind turbine. The numerical framework used in the simulations performs LES under a Cartesian block-structured mesh that is dynamically refined via adaptive mesh refinement (AMR) to increase accuracy and reduce computational costs. The mesh refinement criteria is based on Lagrangian points and vorticity. Solid surfaces of the wind turbine are captured using the Lagrangian formulation of the IB method. The simulations are based on the NREL Phase VI wind turbine, where blade analysis of lift and drag coefficients were performed, followed by time-averaged streamwise velocity and vorticity evaluations. For different angles of attack, the wake centerline velocity analysis demonstrates high energy extraction in the near wake of the wind turbine blade. In terms of validation, the drag coefficient of the NREL S809 airfoil achieves better agreement with the experimental data than the lift coefficient. Validation from the NREL Phase VI shows that the power performance difference between the simulation and experiment is lower than 6%. Wake analysis of the wind turbine scenario are based on downstream profiles of the time-averaged streamwise velocity. The numerical simulation shows a higher loss of kinetic energy in the near wake region, mainly in the centerline, compared to benchmark data, but achieves very similar results to literature in the far-wake region.