Actuator Line Wall-Modeled Immersed Boundary Method for Predicting the Aerodynamic Performance of Wind Turbines

This study addresses the trade-off between accuracy and efficiency in predicting the aerodynamics and wakes of large wind turbines. We developed a unified immersed boundary–actuator line framework with large-eddy simulation. The actuator line efficiently represents blade loading, while the immersed boundary method (IBM) with a wall model resolves near-blade turbulence. The solver uses a staggered Cartesian discretization and is accelerated by a hybrid CPU/GPU implementation. An implicit signed-distance geometry treatment and a ghost cell wall function based on Spalding’s law reduce near-wall grid requirements and eliminate body-fitted meshing. Flow past a three-dimensional cylinder at Re = 3900 validates the accuracy and good grid convergence of the IBM. For the wind turbine, three meshes show converged thrust and torque, with differences below 1% between the two finer grids. At the rated condition (U∞ = 11.4 m/s), thrust and torque agree with STAR-CCM+ and FAST, with deviations of 6.3% and 1.2%, respectively. Parametric cases at 4–10 m/s show thrust and torque increasing nonlinearly with inflow, approximately quadratically, in close agreement with reference models. As wind speed rises, the helical pitch tightens, the wake broadens, and breakdown occurs earlier, consistent with stronger shed vorticity. The framework delivers high fidelity and scalability without body-fitted meshes, offering a practical tool for turbine design studies and extensible wind plant simulations.