Thek-ε-fPmodel applied to double wind turbine wakes using different actuator disk force methods

The newly developed k-ε-fP eddy viscosity model is applied to double wind turbine wake configurations in a neutral atmospheric boundary layer, using a Reynolds-Averaged Navier–Stokes solver. The wind turbines are represented by actuator disks. A proposed variable actuator disk force method is employed to estimate the power production of the interacting wind turbines, and the results are compared with two existing methods: a method based on tabulated airfoil data and a method based on the axial induction from 1D momentum theory. The proposed method calculates the correct power, while the other two methods overpredict it. The results of the k-ε-fP eddy viscosity model are also compared with the original k-ε eddy viscosity model and large-eddy simulations. Compared to the large-eddy simulations-predicted velocity and power deficits, the k-ε-fP is superior to the original k-ε model.