1. Introduction
Offshore wind energy is emerging as the next large-scale energy production device, with better economic conditions, technological development, and favorable policies opening new markets and areas for production. According to the most recent estimates of the International Energy Agency (IEA), total offshore wind capacity is set to more than triple by 2026, reaching almost 120 GW of installed capacity, which would account for one fifth of the total capacity of installed wind energy [
1].
When more wind farms with bigger turbines are being built offshore, the occurrence and importance of different atmospheric conditions are increasing. This makes understanding how this affects the flow dynamics of wind turbines more important in order to optimize energy production and the operation of wind farms.
As mentioned in a review paper by Porté-Agel et al. [
2], it has been shown in numerous papers that the convective boundary layer displays stronger wake meandering and faster wake recovery compared to the neutral and stable atmospheric boundary layers [
3,
4,
5,
6,
7,
8,
9,
10]. Wake meandering is associated with one of the main impacts that stability has on the flows, namely changes in length and velocity scales of the atmospheric turbulence. Faster wake recovery is attributed to the increase in turbulence intensity, another consequence of decreasing stability, as found by Abkar and Porté-Agel [
9]. They found that not only the magnitude but also the spatial distribution of the mean velocity deficit, turbulence intensity, and turbulent momentum fluxes were affected. In the present study, the Mann turbulence model is used to replicate different sized eddies in the different inflows and the corresponding spatial coherence.
There are two primary theories that attempt to explain the meandering of wind turbine wakes. One is that the wake is passively advected by large-scale turbulent motion originating from the inflow [
11]. The other is an observation that wake meandering behind wind turbines exhibits similar behavior to that of vortex shedding behind bluff bodies in high Reynolds number flows. The non-dimensional Strouhal number represents the frequency associated with this vortex shedding, and its expected value for flow around a solid circular disk is around
, the same as for a turbine with a very high tip speed ratio. As the tip speed ratio decreases,
increases. Medici et al. [
12] found experimentally that for a variety of turbine configurations,
was in the range of
. Trivellato et al. [
13] found
to be a recurring number both in their computational fluid dynamics (CFD) study of their wind turbine and in many other studies of seemingly unrelated flows presented in their literature review.
Wake steering shows potential in improving wind farm efficiency and reducing wake losses, both in wind tunnel studies [
14] and field campaigns [
15]. The research into wake steering under different atmospheric stabilities has not been as extensive as for the non-yawed conditions, with a few notable exceptions. Churchfield et al. [
16] performed experimental and numerical measurements at the The Scaled Wind Farm Technology (SWiFT) facility with the aim of measuring wakes and wake deflection resulting from yaw misalignment under a variety of atmospheric conditions. Actuator line large-eddy simulation (AL-LES) simulations are initially run in order to better plan and predict subsequent experiments using light detection and ranging (LiDAR) measurements. These experiments indicated that wake deflection will be strongest under stable to neutrally stratified conditions, and the enhanced mixing of the unstably stratified conditions decreases the amount of expected deflection. Under stable and neutral conditions, a maximum wake deflection of about one third of a rotor diameter is expected 5 rotor diameters (5D) downstream, while this is reduced by roughly half under convective conditions.
In a study by Vollmer et al. [
17], an actuator disc large-eddy simulation (AD-LES) representation of the 5 MW NREL turbine was simulated under different atmospheric conditions. This study suggested that it might not be reasonable to deflect the turbine wake through yaw in unstable conditions, as no correlation was found between the wake position and turbine yaw angle under convective conditions. In a later PhD thesis by Vollmer [
18], wind farm control under different atmospheric conditions using LES was addressed. It was found that parameters such as atmospheric stability, wind veer, shear, and turbulence intensity are important parameters to predict the wake deflection of wind turbines, and that wake deflection can increase the energy yield of a two turbine array in a neutral and stable atmosphere, while the same could not be said about a convective atmosphere. Similar results have also been reported in field campaigns by Fleming et al. [
19] and in LES simulations by Wei and Wan [
20] and Wei et al. [
21].
The Actuator Line model is particularly well-suited for simulating wind turbine wakes under varying atmospheric conditions due to its ability to accurately represent the aerodynamic forces on rotating blades without the computational overhead of fully resolving the blade geometry. This approach enables high-resolution simulations while using fewer grid cells and allowing for larger time steps, which is beneficial for studying large-scale wind farm scenarios with multiple turbines [
22]. By applying the Actuator Line model in this study, we can closely examine the influence of atmospheric stability on wake behavior, particularly under yawed conditions, providing valuable insights for advancing wind farm control strategies.
In the research by Ning and Wan [
23], large-eddy simulation (LES) studies were conducted using OpenFOAM v2206 to examine wind turbine wakes under neutral and unstable atmospheric boundary layer conditions. They presented a comprehensive analysis of wake meandering, including amplitude, spectrum, and probability distribution in both horizontal and vertical planes at various downstream locations. Additionally, they explored the impact of these wake dynamics on the structural loading stability of downstream turbines through correlation analysis and power spectrum evaluation. Rivera-Areba et al. [
24] conducted a comparative study evaluating the engineering model for the DWM against LES in the context of horizontal and vertical steering of wind turbine wakes. Their analysis involved performing LES on the IEA 15 MW wind turbine under various yaw and tilt angles to assess the accuracy and performance of the DWM model.
The allure of capitalizing on the potential energy production gains of wake steering has prompted the construction of several analytical models. Analytical models for calculating the wake trajectory of yawed wind turbines have been proposed by Jiménez et al. [
25], Bastankhah and Porté-Agel [
26], Qian and Ishihara [
27], and Shapiro et al. [
28], in chronological order. Each model was based on previous work, attempting to improve the accuracy of wake deflection estimates. A mixture of theoretical derivation with different approaches, numerical simulations, and wind tunnel experiments were used for development and testing. One of the main distinctions between models is the assumed shape of the distribution of velocity deficit and skew angle, those being a top-hat and Gaussian distribution. The Jimenez model uses the top-hat, which is often used to explain why its predictions differ greatly from the others that use the Gaussian distribution.
Previous work done by the authors tested the effects of the yaw misalignment on two different wind turbines. RANS simulations were verified and validated and compared to the analytical wake deflection models described previously [
29]. Selected figures and work from that study are included in
Appendix C.
With the announced increase in energy production from wind farms, understanding the flow dynamics involved allows for more accurate modeling to help with design and optimization. Wake deflection is likely to be viewed as an attractive option within wind farm control. For this technique to mature, however, one of the questions that needs to be addressed further is how the wake behind yawed wind turbines behaves in non-neutrally stratified atmospheres, which are dominant offshore, and how well analytical models are able to predict this wake deflection.
This paper intends to investigate the wake dynamics of wind turbines with and without yaw under different atmospheric stabilities. This will be done through large-eddy simulations of the Vestas V80 wind turbine with inflow representing stable, neutral, and unstable atmospheres generated using the Mann turbulence generator superimposed onto the wind shear profile. A fast Fourier transform will be performed on the wake meandering, and the energy-containing frequencies and wake center position will be studied 8D downstream of the wind turbine. Finally, the wake deflection trajectories will be compared with the predictions made by analytical models.
The novelty of our research lies in its comprehensive approach to examining wind turbine wake dynamics by integrating multiple factors that have previously been addressed individually. While existing LES studies have explored the impact of atmospheric stability and wake deflection separately, this work uniquely combines these elements. Additionally, the study employs synthetic turbulence generation as the inflow condition, diverging from the conventional LES precursor approach. This choice is particularly significant as it enables the creation of a dataset tailored for direct comparison with existing engineering models, specifically the Dynamic Wake Meandering (DWM) models implemented in tools such as Fast.Farm, HAWC2, and DIWA. By doing so, the research not only enhances the understanding of wake behavior under varying atmospheric conditions and yaw but also provides a valuable benchmark for validating and improving current engineering models used in wind turbine design and analysis.
Since these numerical simulations are strictly aerodynamic representations of the wind turbines, the second big question that needs to be addressed for adoption of yaw steering strategies cannot be answered here, namely the structural loads on the wind turbine subject to different wake conditions. Different turbulence generators are not investigated either in the development of this model, and the inflow parameters are based on findings from Riverra-Arreba et al. [
30].
5. Numerical Simulations
CFD simulations are performed using OpenFOAM in order to analyze the flow dynamics behind wind turbines. The flow field is solved using the OffWindSolver, which is a solver under development by Balram Panjwani at SINTEF Industry [
42]. Both RANS and LES simulations are used in this work. The RANS simulations use a RNG (Re-Normalization Group)
k-
turbulence model, based on the authors’ previous experience of this model’s suitability for solving wind turbine flows [
29]. It models the effect of turbulence as an added viscosity, the effect of which is introduced into the NSEs. Separate transport equations for turbulent kinetic energy
k and turbulent dissipation rate
are solved to find the turbulent viscosity. The RNG
k-
model is an extension of the standard
k-
model, incorporating additional modifications that work to provide better accuracy in capturing the flow physics and reduce the dependency on user-defined constants. The LES simulations use the Smagorinsky turbulence model to close the NSE equations [
43]. The Smagorinsky turbulence model is a method used in LES to approximate the effects of the unresolved small-scale turbulence on the simulated larger-scale flow. It is also an eddy viscosity model and therefore assumes the turbulent stresses are proportional to the local deformation of the flow. Their magnitude is determined using the Smagorinsky coefficient. Its value is usually based on empirical or theoretical considerations, and it affects the accuracy of the model.
5.1. RANS Simulations
RANS simulations were initially run in order to verify and validate the wind turbine models used for this study, in addition to being a building block towards the higher-fidelity LES simulations.
5.1.1. Mesh Description
Mesh parameters are given in
Table 5 along with pictures of the mesh in
Figure 4a,b for the Samsung turbine.
is the domain length and
the number of cells in the
directions.
refers to the cell-size-lengths, which are equal in the
x,
y, and
z directions, in the region of the mesh where the turbine is located. Near the ground, the cells are refined enough to capture and maintain the inflow shear profile, while at the top and back, cells are expanded. This reduces computational cost in the regions that are of less importance.
5.1.2. Verification
A total run-time of 700 and 1500 s for each simulation proved to give a converged solution with respect to time for the Vestas and Samsung turbines, respectively. A case with a reduced time step was also run to check the simulations sensitivity to a reduction in residuals. It showed a minimal effect on both wake, thrust, and power for both turbines. Residuals along with results for velocity deficit, turbine thrust, power, and time convergence can be found in
Appendix A.
A mesh refinement study was performed to ensure that grid convergence for the RANS simulations was achieved. Three meshes of increasing refinement named coarse, orig, and fine were initially run for both turbines with minimum cell length
in the refined region. The cell length
was changed by a factor of
between each mesh for the Vestas turbine and
for the Samsung turbine. Based on these results, a fourth mesh, named medium, with a refinement level in between coarse and orig was proposed and verified for the V80 turbine, while the orig mesh was deemed adequate for the Samsung turbine.
Figure 5a,b show the results of a GCI study performed using the medium and fine meshes, while output for power and thrust from the rotors is presented in
Table 6.
Shear in the inflow profile proved to have a pronounced effect on the velocity deficit in the wake, as can be seen in
Figure 6. The shear profile was therefore included in subsequent iterations of the models. A no-slip boundary condition and wall functions were applied at the bottom wall, and, in order to capture the steep gradient, cells were refined in this lower region. A sufficiently low value for the non-dimensional wall distance y+ needed for an accurate near-wall solution was never achieved. However, a plot of the wind velocity profile for locations downstream of the inlet and upstream of the turbine shows good agreement with the input velocity profile, especially in the region where the turbine operates; see
Figure 7. This analysis was done for the Vestas turbine because of the availability of inflow data from Keck et al. [
7]. Changes based on the analysis were made to both turbine models.
5.1.3. Validation
RANS simulations in uniform inflow were run at different wind speeds to validate turbine thrust and power. The results were compared to performance data given by the manufacturer and found in literature and are given in
Figure 8 and
Figure 9. It can be seen in these figures that both thrust and power are modeled accurately. In
Figure 9b, it can be seen that at a wind speed of 10 m/s, the thrust is overestimated. This is likely because the turbine would approach rated wind speed and go into a different operating region, which would mean altering the TSR. Tip speed ratio for the V80 turbine is
, and for the Samsung it is as described earlier in
Section 2.
LES data from Keck et al. [
7] of a Vestas V80 turbine was used to validate the wake of the RANS case. Keck et al. [
7] used a hub-height wind speed of
m/s in an unstable atmosphere with turbulence intensity
. A shear profile for the RANS simulation was generated with the power law using
as found by Rivera-Arreba et al. [
30] for an unstable atmosphere; see
Table 4. Another parameter with significant impact is the turbulent dissipation rate
prescribed at the inlet. It is based on a turbulent length scale
. Setting
gave the best agreement with validation data (see
Figure 10a) and will be used in further simulations. The Samsung turbine was validated with LiDAR data from the TotalControl project [
31], where measurements are done on the Levenmouth turbine operated by ORE Catapult. A rather large turbulent length scale of
m, or
(1Radius), gave a solution most similar to that of the validation data, as shown in
Figure 10b. It should be noted that the data used to validate this turbine is quite noisy, and there is an absence of information describing the inflow conditions when the measurements were made. With a basis in the information given from the TotalControl project, a wind speed of 8 m/s and a turbulence intensity of 10% were therefore assumed, while the shear profile was imposed in a similar manner to that of the Vestas turbine.
5.2. LES
A higher-fidelity simulation was run using LES for the 2 MW Vestas V80 turbine. It was verified and then validated using the same data as previously from Keck et al. [
7] before being used to generate results. The Samsung turbine was not used for the LES simulations primarily due to lack of validation data, as explained in
Section 5.3. For the results, six simulations were run, namely with and without a yaw angle of
for each of the three stability conditions presented earlier in
Section 4.4. The
yaw angle was chosen based on previous research by the authors and is expected to produce a clear wake displacement with downstream flow structures, which has been observed in literature. This angle has also been among the yaw angles studied in the development of the analytical models by Bastankhah and Porté-Agel [
26], Shapiro et al. [
28], and Jiménez et al. [
25].
5.2.1. Mesh Description
The new mesh has the same domain and cell size as the one used for RANS, only with added refinement regions around the turbine and wake. The reason being that, whilst RANS is designed for coarser meshes, LES needs a finer mesh in order to resolve eddies at smaller length scales. Running an LES simulation with the RANS solutions interpolated onto the finer mesh makes for a more accurate result. The mesh was also changed to be uniform in the
x,
y-, and
z-directions. The refinement regions were created using the openFOAM utility
snappyHexMesh, and the resulting cubic cells are described in
Table 7. Pictures of the mesh are presented in
Figure 11. For the simulations of yawed turbines, both wake refinement regions were moved 30 m in the negative
y-direction, and the turbine refinement region was expanded 13 m in front of and behind the turbine. This ensures that both the turbine and the wake are still well captured by the finer cells. Due to the angle of the turbine, the wind velocity experienced by the rotor is decomposed by the cosine of the yaw angle. With the new wind speed assumed to be the component normal to the plane of the turbine, a new rotational speed was calculated keeping the tip speed ratio constant.
5.2.2. Verification
A simulation was run on a finer mesh with the number of cells increased by a factor of
in the
x,
y, and
z directions. The two solutions were then used to calculate the grid convergence index for the velocity deficit downstream of the turbine. Along with values for turbine output, the GCI study shows that a finer mesh is not needed; see
Figure 12 and
Table 8.
A total run time of 3000 s proved to give a converged solution with respect to wake, thrust, and power. A case with a halved time step of
was run to test the model’s sensitivity to this parameter. This also had the effect of lowering the residuals. Results show that
and the corresponding residuals are sufficiently low, which is shown in
Appendix A.
Both the shear profile and turbulence intensity of the simulation were checked for positions upstream of the turbine. It is shown in
Figure 13 that shear profiles for all positions and stabilities gave good agreement with the desired shear-profile shape. Upon inspection, turbulence intensity was found to decrease substantially from the inlet to the turbine. To correct for this, intensity at the inlet was increased.
Table 9 gives an overview of the achieved and target values for turbulence intensities.
In addition to the parameters presented above, a number of other simulations were run in order to examine the numerical models sensitivity and reliance on these parameters. This included comparing the Smagorinsky and a dynamic one equation eddy-viscosity turbulence models, testing for different values for the Smagorinsky coefficient, examining different damping functions within the Smagorinsky model and running a simulation with a larger domain. Of these, it was found that the turbulence model, the Smagorinsky coefficient, and the extended domain influenced the solution to a very small degree. The choice of damping function provided in openFOAM did impact the solution more, and among the Prandtl, van Driest, and cube-root volume damping functions, the simpler cube-root volume damping functions were shown to give results most representative of the validation case.
5.2.3. Validation
Simulated turbine thrust and power for all stabilities at 0 degrees yaw are compared with values from the turbine manufacturer given in
Table 10. The wake deficit is then compared with Keck et al. [
7], following the same procedure as in
Section 5.1.3, and presented in
Figure 14. Ideally, field measurement data or a higher fidelity model would be used to validate this LES model, thereby avoiding one LES model being used to validate another LES model. Time restrictions and data availability made this unachievable, but the validation results are an indication of a realistic turbine wake nonetheless.
5.3. Samsung Turbine
As previously stated in
Section 2 and
Section 5.1.3, the numerical Samsung turbine matches the performance expected from literature well. Even though the LiDAR measurements are noisy, the wake of the Samsung turbine was able to replicate the measured data to a reasonable degree with a turbulent length scale of one rotor radius for the RANS simulations (
Figure 10b). For future implementations of the turbine, the tilt angle should be smaller. The value of this parameter followed from the 5 MW NREL turbine and was not altered. It can be seen in
Figure 4b that the wake and wake center of the turbine have a downward trajectory, which is likely due to this high tilt angle. The wake center not going parallel to the horizon could also be a reason as to why this turbine needs such a high turbulent length scale to match the validation data.
Ultimately, however, the choice was made to use the V80 Vestas turbine in favor of the Samsung S7.0-171 turbine. This came down to the fact that for the Vestas turbine, the validation data were given for well-described inflow conditions, while the same could not be said for the Samsung turbine. Seeing that in many cases the wind speed was not well defined and that the turbulence intensity was not specified for any of the measurements, the authors could not validate the turbine with sufficient confidence to pursue high-fidelity LES simulations aiming at investigating the influence of atmospheric conditions on the flow.
However, the rotor performance of the turbine compares well with literature, and the RANS simulations give reasonable agreement with the TotalControl data. Other publications may therefore find the information useful, especially with access to validation data with a better description of the inflow.
7. Conclusions
Understanding the impact of atmospheric conditions on wind turbine flow dynamics is crucial for optimizing energy production and operational efficiency. Offshore wind farms introduce challenges related to larger turbines and different operating conditions. The effect of atmospheric stability on wind turbine wake dynamics has been examined in this work, in particular wake deflection and meandering. High-fidelity numerical AL-LES simulations with Mann generated turbulent inflows were used to investigate three different stabilities for turbines with and without yaw.
An attempt was made to reproduce and use the Samsung S7.0-171 in the study. A successful creation of a numerical representation of the turbine has been achieved, which can be useful for later studies. However, due to a high degree of uncertainty in some parts of the validation data that are important for this work, sufficient accuracy of the AL-LES model could not be guaranteed, and the Vestas V80 turbine was used instead.
As anticipated, with increasing atmospheric stability and a consequent decrease in turbulence intensity, the wake deficit also increases. Despite these differences observed for wake deficit between the neutral and stable cases, they have very similar deflected wake trajectories when the turbines are in yaw. Both exhibit a greater initial wake deflection than what is expected from the analytical wake deflection models, but this initial sharp gradient flattens out, and at 8D, good agreement is found between the numerical and analytical results. At 5D, the wake center position was found to be 0.29D, which is in excellent agreement with Churchfield et al. [
16], who reported a deflection of 0.33D for the same downstream distance. In the furthest wake position examined, at 8D, deflection increases marginally to 0.375D, exposing a new potential turbine at this location to partial wake conditions. The effect of this is made more substantial when considering that the oscillating motion from meandering has a magnitude equal to 0.25D. It was also found that the wake trajectories were very similar if the mean velocity was used to find wake centers or if the wake center position was found using the mean value of the meandering wake. This would suggest that the meandering behind yawed turbines occurs around the deflection axis.
Meandering analysis was conducted at 8D for neutral and stable stabilities, while the unstable case used 4D due to the absence of a clearly defined wake. The Welch method with specific parameters provided results for cutoff frequency, peak frequency, and Strouhal number. In the neutral atmosphere, meandering behavior was similar between yawed and non-yawed turbines, while slightly larger amplitudes were observed for yaw compared to non-yaw in the stable case. Neutral simulations exhibited larger amplitudes compared to stable, non-yawed simulations. Meandering frequency appeared higher in both stable cases, supported by the peak frequency results. Results at 4D in the neutral and stable cases were similar to those at 8D, indicating some independence from downstream position. The results for stable and neutral inflows aligned well with the analytical cutoff frequency using the wake diameter. Strouhal number findings were consistent with previous studies for stable and neutral cases, but significantly lower for the unstable case. This could suggest different dynamics; however, considering the model limitations and challenges in determining the wake center, further inspection and verification are needed to draw conclusive insights.