the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Changing the rotational direction of a wind turbine under veering inflow: a parameter study
Antonia Englberger
Julie K. Lundquist
Andreas Dörnbrack
All currentday windturbine blades rotate in clockwise direction as seen from an upstream perspective. The choice of the rotational direction impacts the wake if the wind profile changes direction with height. Here, we investigate the respective wakes for veering and backing winds in both hemispheres by means of largeeddy simulations. We quantify the sensitivity of the wake to the strength of the wind veer, the wind speed, and the rotational frequency of the rotor in the Northern Hemisphere. A veering wind in combination with counterclockwiserotating blades results in a larger streamwise velocity output, a larger spanwise wake width, and a larger wake deflection angle at the same downwind distance in comparison to a clockwiserotating turbine in the Northern Hemisphere. In the Southern Hemisphere, the same wake characteristics occur if the turbine rotates counterclockwise. These downwind differences in the wake result from the amplification or weakening or reversion of the spanwise wind component due to the effect of the superimposed vortex of the rotor rotation on the inflow's shear. An increase in the directional shear or the rotational frequency of the rotor under veering wind conditions increases the difference in the spanwise wake width and the wake deflection angle between clockwise and counterclockwiserotating actuators, whereas the wind speed lacks a significant impact.
 Article
(8922 KB)  Fulltext XML
 BibTeX
 EndNote
The copyright of the authors Antonia Englberger and Andreas Dörnbrack for this publication are transferred to Deutsches Zentrum für Luft und Raumfahrt e. V., the German Aerospace Center. This work was authored (in part) by the National Renewable Energy Laboratory, operated by Alliance for Sustainable Energy, LLC, for the US Department of Energy (DOE) under contract no. DEAC3608GO28308. Funding provided by the US Department of Energy Office of Energy Efficiency and Renewable Energy Wind Energy Technologies Office. The views expressed in the article do not necessarily represent the views of the DOE or the US Government. The US Government retains and the publisher, by accepting the article for publication, acknowledges that the US Government retains a nonexclusive, paidup, irrevocable, worldwide license to publish or reproduce the published form of this work, or allows others to do so, for US Government purposes.
Most modern industrialscale wind turbines rotate clockwise, as seen from a viewer looking downwind. Traditional Danish windmills turned counterclockwise (Maegaard et al., 2013), as they were built by righthanded millers who preferred the thin end of the laths to be pointing towards the left on the blades. This rotational direction was adapted by the windturbine pioneer Christian Riisager and subsequently by the company Tvind. In 1978, Erik GroveNielsen designed the first 5 m fiberglass blades. He and his wife Tove chose a clockwise rotational direction of the blades to distinguish their product from Tvind. Descendants of the Riisager wind turbine (Windmatic and Tellus) rotate counterclockwise, while those of GroveNielsen (Vestas, Bonus – now Siemens, Nordtank, and Enercon) rotate clockwise. Three of the four clockwiserotating blade manufacturers became market leaders in the international wind power industry, and the clockwiserotating blades, eventually, became the global standard (Maegaard et al., 2013). The clockwise blade rotation is, therefore, a historical coincidence without physical motivation.
Rotating blades encounter a variety of wind conditions. In a convective regime during the daytime above the surface layer, there is no significant change in the incoming wind direction or wind speed with height and the inflow conditions are uniform over the whole rotor area. A nocturnal stably stratified regime, however, often generates wind profiles with changing magnitude (vertical wind shear) and direction (wind veer; Lindvall and Svensson, 2019). Vertical variations in both quantities reflect the balance between Coriolis force and friction. Friction affects the lowest part of the wind profile and contributes as internal friction in the flow while the rotational direction of the wind vector in the Ekman spiral aloft depends on the hemisphere. In the Northern Hemisphere (NH; Southern Hemisphere, SH), winds tend to rotate clockwise (counterclockwise) with height (Stull, 1988). Veer occurs on many nights both onshore (Walter et al., 2009; Rhodes and Lundquist, 2013; Sanchez Gomez and Lundquist, 2020) and offshore (Bodini et al., 2020, 2019). According to 2 years of meteorological tower measurements in Lubbock (Texas; Walter et al., 2009) and 3 months of lidar observations in northcentral Iowa (Sanchez Gomez and Lundquist, 2020), veer occurs in well over 70 % of those stable boundary layer (SBL) occurrences (≈ 76 % in Walter et al., 2009, and ≈ 78 % in Sanchez Gomez and Lundquist, 2020). In the remaining 22 % (Sanchez Gomez and Lundquist, 2020) to 24 % (Walter et al., 2009), a backing wind occurs. A backing wind is characterized by a counterclockwise wind direction change with height in the NH.
The frequency of occurrence of a veering wind depends on many criteria. A wind direction change with height occurs mainly at night. Secondly, seasonal differences occur. Thirdly, the frequency of occurrence of veering or backing is location specific. In Lubbock (Texas; Walter et al., 2009) and in northcentral Iowa (Sanchez Gomez and Lundquist, 2020), a veering wind occurs on three out of four nights. In their global climatology of veer based on radiosonde data, Lindvall and Svensson (2019) find stronger veer (or backing for the SH) in midlatitudes (see their Fig. 3). Of course, topography can change the frequency of occurrence significantly. Each location has its own percentage values of the occurrence of a veering wind. In this study, a directional shear of 0.08${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$ is applied in the reference case, as it corresponds to the mean of the frequency of occurrence of a veering inflow in Walter et al. (2009). Further, seasonal dependence occurs (Bodini et al., 2019, 2020). The longerlasting nights during winter are characterized by smaller mean values of the directional shear (minimum winter values in December of 0.03${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$, according to 13 months of offshore lidar measurements in Massachusetts; Bodini et al., 2020, Fig. 4). The shorter nights during summer, however, are characterized by larger mean values of the directional shear (maximum summer values in June of 0.095${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$; Bodini et al., 2020, Fig. 4). The occurrence of a specific directional shear (e.g., ds = 0.08${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$) for a specific wind speed (between 9 and 11 m s^{−1}, similar to the reference case used here) is much larger during winter (85 % of the veering cases are characterized by ds ≤ 0.08${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$) in comparison to the summer (45 %; Bodini et al., 2020, Fig. 5).
The wind turbine's wake characteristics in a veering wind regime differ for counterclockwise and clockwiserotating blades (Englberger et al., 2020). The induced vortex component of the near wake's flow is determined by the rotation of the blades. The wake rotates opposite to the blade rotation due to aerodynamics and design of the windturbine blades (Zhang et al., 2012). In contrast, the rotational direction of the far wake is determined by the Ekman spiral. If a northern hemispheric Ekman spiral interacts with clockwiserotating blades, the spanwise flow component in the wake weakens or even reverses due to a superposition with the vortex of the near wake and attenuates to the inflow in the far wake. In this case, the near wake's counterclockwise rotation diminishes and becomes clockwise. After this reversion or likewise in the case of a reduction in the spanwise wake component, the wake's rotation strength intensifies downwind. Conversely, if the same inflow interacts with counterclockwiserotating blades, the spanwise flow component is amplified in the near wake, because the rotational direction persists in the whole wake and the wake's rotation strength weakens downwind.
The modification of the spanwise flow component also impacts the streamwise velocity in the wake. It affects the velocity deficit in the near wake, the streamwise wake elongation of the wake, the spanwise wake width, and the deflection angle of the wake (Englberger et al., 2020). There also exists a rotationaldirection impact on the streamwise velocity component in a flow regime without wind veer (Vermeer et al., 2003; Shen et al., 2007; Sanderse, 2009; Kumar et al., 2013; Hu et al., 2013; Yuan et al., 2014; Mühle et al., 2017). VaselBeHagh and Archer (2017) simulated a wind farm similar to Lillgrund in Sweden with alternative rotational direction of the rotors, starting with clockwise in the first row. Including wind shear but no wind veer in the inflow conditions, the power output was 1.4 % larger in comparison to only clockwiserotating rotors in the wind farm. However, compared to the wake differences for clockwise and counterclockwiserotating rotors in a flow regime with wind veer, the differences are small in the case of no wind veer (Englberger et al., 2020). Therefore, the 1.4 % in VaselBeHagh and Archer (2017) can be considered a lower limit; the consideration of veer amplifies this difference.
In this study, we investigate the relationship between the upstream wind profile and the direction of the turbine rotation using largeeddy simulations (LESs). Both clockwise and counterclockwiserotating actuators are embedded in a veering as well as a backing inflow for both hemispheres. In the case of a veering inflow in the NH, we carry out a parameter study investigating the impact of the magnitude of the geostrophic wind, the directional shear, and the rotational frequency of the rotor. The results of the rotationaldirection impact on the wake are interpreted for all simulations with a theoretical analysis considering a Rankine vortex representation of the wake. To our knowledge, this is the first parameter study which investigates the interactions of wake rotational direction in combination with an Ekman spiral.
Our previous study (Englberger et al., 2020) lays the groundwork for this investigation, describing in detail the rotationaldirection impact in a stably stratified regime under veered inflow conditions and in an evening boundary layer regime under nonveered conditions (in the Northern Hemisphere). Further, that work explains the physical mechanism responsible for the rotationaldirection impact of the blades on the wake by simple analysis of a linear superposition of the veering inflow wind field with a Rankine vortex.
This paper is organized as follows. The numerical model EULAG, the windturbine simulation setup, and the metrics applied in this work are described in Sect. 2. The analysis predictions are introduced in Sect. 3. The corresponding idealized simulations investigating the rotationaldirection impact on the wake follow in Sect. 4. A comparison of the simulation results to the analysis predictions is given in Sect. 5, and a conclusion follows in Sect. 6.
2.1 The numerical model EULAG
The windturbine simulations, with prescribed wind and turbulence conditions, are conducted with the flow solver EULAG (Prusa et al., 2008). For a comprehensive description and discussion of EULAG we refer to Smolarkiewicz and Margolin (1998) and Prusa et al. (2008).
The Boussinesq equations for a flow with constant density ρ_{0} = 1.1 kg m^{−3} are solved for the Cartesian velocity components u, v, and w and for the potentialtemperature perturbations Θ^{′} = Θ − Θ_{e} (Smolarkiewicz et al., 2007):
with Θ_{e} representing the environmental or background state and Θ_{0} representing the constant reference value of 300 K. In Eqs. (1)–(3), d∕dt, ∇ and ∇ ⋅ represent the total derivative, the gradient, and the divergence, respectively. The quantity p^{′} represents the pressure perturbation with respect to the environmental state. Further, g represents the vector of acceleration due to gravity. The subgridscale terms 𝓥 and ℋ symbolize viscous dissipation of momentum and diffusion of heat. F_{WT} corresponds to the turbineinduced force and β_{v} to the rotational direction. All following simulations are performed without an explicit subgridscale closure as implicit LESs (Grinstein et al., 2007), to remove any question of the influence of the subgridscale closure on the results. Further, we apply a freeslip vertical boundary condition.
The turbineinduced forces (F_{WT}) in Eq. (1) are parametrized with the blade element momentum (BEM) method as the actuator disc, including a nacelle and excluding the tower. The BEM method enables the calculation of the steady loads, thrust, and power for different wind speeds and rotational speeds of the blades. The airfoil data of the 10 MW reference wind turbine from DTU (Bak et al., 2013) are applied, whereas the radius of the rotor as well as the chord length of the blades is scaled to a rotor with a diameter of 100 m. For a more detailed description of the windturbine parametrization and all values used in the windturbine parametrization, we refer to parametrization B of Englberger and Dörnbrack (2017).
The actuator disc rotates in a clockwise or counterclockwise direction, depending on the choice of β_{v}∈{−1, 1}. The rotor rotation is not directly simulated; instead, the rotor forces are exerted directly on the velocity fields in Eq. (1). A clockwiserotating rotor initiates a counterclockwise wake rotation and vice versa, following conservation of angular momentum (Zhang et al., 2012). In this work, a common clockwise rotor rotation “cr” is defined as β_{v} = 1 and β_{w} = −1 and a counterclockwise rotor rotation “ccr” as β_{v} = −1 and β_{w} = 1, with β_{u} = 1 in both cases.
2.2 Setup of the windturbine simulations
Windturbine simulations on 512 × 64 × 64 grid points with a horizontal and vertical resolution of 5 m and open boundaries are performed for veering and nonveering inflow lasting 40 min. The rotor of the wind turbine has a diameter D as well as a hub height z_{h} of 100 m and is located at 300 m downwind from the inflow boundary and centered in the spanwise y direction.
A veering wind profile can be described by the Ekman spiral:
following Stull (1988), with a geostrophic wind u_{g} and
representing a Coriolis parameter f = 1.0 × 10^{−4} s^{−1} and an eddy viscosity coefficient κ.
Wind direction change between two heights is defined as directional shear. As we assume the directional shear to be an impact factor for the interaction process of a rotating system with veering inflow, a modified version of the Ekman spiral is applied as v_{f} in the simulations in this work. Further, the negative vertical gradient of the streamwise velocity in the supergeostrophic component of the Ekman spiral is not considered in u_{f}.
The simulations are initialized with the streamwise velocity profile
with an eddy viscosity coefficient κ = 0.06 m^{2} s^{−1}. The corresponding spanwise velocity profile is
in the case of no veering inflow with $\frac{\partial {v}_{f}}{\partial z}$ = 0 and
in the case of veering inflow with $\frac{\partial {v}_{f}}{\partial z}$ ≠ 0, with a given directional shear
with Δϕ = ϕ_{150 m}−ϕ_{50 m} and
in the lowest 200 m and constant above. The influence of the Coriolis force on the flow field is only included in the simulations via Eqs. (6) and (8). Note that no Coriolis force is applied in the numerical model (Eq. 1).
For u_{f} and v_{f}, we consider the NH (f>0) and the SH (f<0) and a veering (Δϕ>0 in NH; Δϕ<0 in SH) and a backing (Δϕ<0 in NH; Δϕ>0 in SH) wind. In the reference simulation (with a veering wind in the NH), the directional shear is 0.08${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$ with v_{f}(z_{h}) = 0. The initial vertical velocity is
in all simulations. The flow components u_{f}, v_{f}, and w_{f} are used for specifying the initial conditions. The pressure solver in EULAG further applies u_{f} and v_{f} as boundary conditions. The potential temperature is
in the lowest 200 m and 303 K above.
For a veering wind in the NH, we modify the geostrophic wind, the directional shear, and the rotational frequency of the rotor. We sample winds with a geostrophic wind component of u_{g} = 6 m s^{−1}, u_{g} = 10 m s^{−1} (reference simulation), and u_{g} = 14 m s^{−1}, the first and last of which are referred to with the acronyms u6 and u14 in the simulation nomenclature. Further, we apply a directional shear of 0.04, 0.08, 0.12, 0.16, and 0.20${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$ corresponding to weak (ds4), moderate (reference simulation), moderate to strong (ds12), strong (ds16), and very strong (ds20) shear. As an additional parameter, the rotational frequency ranges from Ω = 0.058 s^{−1}, Ω = 0.12 s^{−1}, and Ω = 0.175 s^{−1} to Ω = 0.23 s^{−1}, corresponding to low (Ωl), moderate (reference simulation), high (Ωh), and very high (Ωvh) in the simulation nomenclature.
This work is a parameter study investigating the impacts of the inflow (directional shear, wind speed) and the rotating system (rotational frequency) on the wake. The wake's impact by the rotational direction of the rotor depends on the mean wind profile, which is determined by the geostrophic wind (Eqs. 6 and 8) and the directional shear (Eq. 9). Turbulence modifies the strength of the wakes but not the occurrence (Appendix and Sect. 5). Therefore, we perform the simulations as implicit LESs with no explicit subgridscale closure model. Moreover, we apply the turbulence parametrization by Englberger and Dörnbrack (2018b) to perturb the flow field during the numerical integration. This turbulence parametrization provides a computationally fast method for windturbine simulations with open horizontal boundary conditions in a small domain. It includes stabilitydependent atmospheric characteristics in the inflow. This makes the method very suitable for parameter studies. We superimpose upon the inflow wind field turbulent fluctuations in a neutral boundary layer precursor simulation (Englberger and Dörnbrack, 2017), as represented by the term I in Eq. (13), where ${\mathit{u}}_{p}{}_{{i}^{\ast},j,k}$ is the velocity vector of a neutral boundary layer equilibrium state at each grid point i, j, and k.
The indices of the grid points are denoted by i = 1 … n, j = 1 … m, and k = 1 … l in the x, y, and z directions, respectively. The star refers to a streamwise shift by one grid point at every time step δ with i^{∗} = i+δ, whereas i^{∗} = [1,n] and δ^{∗} represents the passed number of time steps. The prefactor α_{0} represents the amplitude of the turbulence perturbations and ${\mathit{\alpha}}_{{i}^{\ast},j,k}$ represents adjustable stratificationdependent parameters for convective and stable regimes as well as the transitions between them. The stratificationdependent parameters were retrieved from a 30 h diurnalcycle simulation from Englberger and Dörnbrack (2018a).
In the following simulations we apply nighttime representations using values of α = 0.3, α_{u} = 0.15, α_{v} = 0.24, and α_{w} = 0.13 (Englberger and Dörnbrack, 2018b, Table 1). A rather similar setup including the turbulence parametrization has been applied in Englberger and Lundquist (2020).
2.3 Metrics
For the investigation of the rotationaldirection impact on the wake, the following characteristics are calculated from the simulation results: the spatial distribution of the timeaveraged discrete streamwise velocity $\stackrel{\mathrm{\u203e}}{{u}_{i,j,k}}$, the timeaveraged discrete spanwise velocity $\stackrel{\mathrm{\u203e}}{{v}_{i,j,k}}$, and the streamwise velocity deficit
The characteristics are averaged over the last 30 min of the 40 min windturbine simulation. The 30 min temporal average is calculated online in the numerical model according to the method of Fröhlich (2006, Eq. 9.1).
In the following, the quantities $\stackrel{\mathrm{\u203e}}{{u}_{i,j,k}}$ and $\stackrel{\mathrm{\u203e}}{{v}_{i,j,k}}$ are evaluated and discussed for top and bottom sectors. They result from a division of the rotor area into four sections of 90^{∘}, as shown in Fig. 1, including all grid points with a distance r from the rotor center 0 m < r ≤ R. The left and right sectors are defined from a view looking downwind towards the wind turbine on the disc.
In this section, we construct a simple analytical model of the interaction of the rotating wake of a wind turbine with a sheared ambient flow. The rotating wake is prescribed by a Rankine vortex, whereas the ambient flow is described by three different inflow conditions (no veer, veering wind, backing wind). In the case of a veering inflow, the relations are also evaluated for three different parameters (wind speed, directional shear, rotational velocity). The approach follows Englberger et al. (2020) and is modified to allow different directionalshear values.
A rotating system can be described by a Rankine vortex with the radial dependence r, the rotational velocity ω, and the angle ϑ in a y–z plane:
The veering inflow is described by Eqs. (4) and (5), whereas no wind veer is described by Eq. (7). In this analysis we apply the simplified Eqs. (6) and (8) for the veering inflow, as they allow a variety of directionalshear values. Both inflow cases result in a superposition of the spanwise components v_{f} (Eqs. 7 and 8) and v_{v} (Eq. 15) in Eq. (17):
In Eq. (17), a linear decrease in v_{v}(z,x_{down}) is assumed for a given downwind distance x_{down} from the rotating system x_{WT} up to x_{ξ} with v(z,x_{ξ}) = 0. In this work, we only consider the spanwise flow component, as w_{f} = 0 (Eq. 11).
Figure 2 represents the spanwise velocity component v resulting from Eq. (17) at z = 125 m and at z = 75 m at the rotor center in a lateral direction. The rotating system has a rotor center z_{h} = 100 m and a rotor radius R = 50 m. The vertical positions are centered in the top and the bottom sectors of Fig. 1.
The wake resulting from a clockwise (cr) or counterclockwise (ccr) rotating rotor interacting with no wind veer are represented in Fig. 2b. Following Eq. (17), the rotational direction of the rotor determines the sign of v(z,x_{down}). Therefore, the spanwise velocity component has the opposite sign in the top and the bottom rotor part of both cr and ccr in Fig. 2b. Approaching x_{ξ}, v(z,x_{ξ}) = v_{f}(z) = 0.
In the case of veering inflow, however, the spanwise flow component impacts the wake (Eq. 17). The spanwise flow component results from the Ekman spiral, which is hemispheric dependent. In the NH, f > 0, and, therefore, the spanwise flow component ${v}_{f}({z}_{\mathrm{h}}R/\mathrm{2})$ > 0 in the lower rotor half and ${v}_{f}({z}_{\mathrm{h}}+R/\mathrm{2})$ < 0 in the upper rotor half with v_{f}(z_{h}) = 0 (Eq. 8). This situation corresponds to a flow from right to left in the lower rotor half and from left to right in the upper rotor half, looking from upwind towards downwind (Fig. 1). If “+” is applied in Eq. (17), ${v}_{v}({z}_{\mathrm{h}}R/\mathrm{2},{x}_{\text{down}})=+\mathit{\omega}r\mathrm{sin}\left(\mathrm{270}{}^{\circ}\right)=\mathit{\omega}r$ < 0 in the lower rotor half and ${v}_{v}({z}_{\mathrm{h}}+R/\mathrm{2},{x}_{\text{down}})=+\mathit{\omega}r\mathrm{sin}\left(\mathrm{90}{}^{\circ}\right)=\mathit{\omega}r$ > 0 in the upper rotor half. However, if “−” is applied in Eq. (17), ${v}_{v}({z}_{\mathrm{h}}R/\mathrm{2},{x}_{\text{down}})=\mathit{\omega}r\mathrm{sin}\left(\mathrm{270}{}^{\circ}\right)=\mathit{\omega}r$ > 0 in the lower rotor half and ${v}_{v}({z}_{\mathrm{h}}+R/\mathrm{2},{x}_{\text{down}})=\mathit{\omega}r\mathrm{sin}\left(\mathrm{90}{}^{\circ}\right)=\mathit{\omega}r$ < 0 in the upper rotor half. The sign “+” corresponds to a counterclockwise wake rotation which arises from a clockwise rotor rotation (Zhang et al., 2012), whereas the sign “−” corresponds to a clockwise wake rotation arising from a counterclockwise rotor rotation ccr.
In the case of a clockwiserotating rotor, the rotor competes against the veer effect with ${v}_{v}({z}_{\mathrm{h}}R/\mathrm{2},{x}_{\text{down}})$ < 0 superpositioning ${v}_{f}({z}_{\mathrm{h}}R/\mathrm{2})$ > 0 and ${v}_{v}({z}_{\mathrm{h}}+R/\mathrm{2},{x}_{\text{down}})$ > 0 superpositioning ${v}_{f}({z}_{\mathrm{h}}+R/\mathrm{2})$ < 0. In both the top and bottom half of the rotor, the spanwise component of the inflow v_{f} is weakened by the vortex component v_{v} or even reversed if $\mid {v}_{v}\mid $ > $\mid {v}_{f}\mid $. Approaching downwind, the impact of v_{v} decreases and v(z,x_{down}) approaches v_{f}(z) at x_{down} = x_{ξ} with v_{v}(z,x_{ξ}) = 0.
In the case of a counterclockwiserotating rotor, the wake vortex intensifies the inflow with ${v}_{v}({z}_{\mathrm{h}}R/\mathrm{2},{x}_{\text{down}})$ > 0 superpositioning ${v}_{f}({z}_{\mathrm{h}}R/\mathrm{2})$ > 0 and ${v}_{v}({z}_{\mathrm{h}}+R/\mathrm{2},{x}_{\text{down}})$ < 0 superpositioning ${v}_{f}({z}_{\mathrm{h}}+R/\mathrm{2})$ < 0. The vortex intensifies the inflow v_{f}(z_{h}) at all rotor heights. Approaching downwind, the impact of v_{v} decreases and v(z,x_{down}) approaches v_{f}(z) at x_{down} = x_{ξ} with v_{v}(z,x_{ξ}) = 0. At x_{down} = x_{ξ}, the situation is independent of the vortex and the wake has completely recovered.
The different behavior of the spanwise wake component is presented in Fig. 2e. In the case of cr, the vortex component weakens the spanwise inflow component, resulting in a reversion of the sign of v(z,x_{down}) behind the rotor at x_{down} < x_{ξ}. In the case of ccr, however, the vortex component intensifies the spanwise inflow component. At x_{ξ}, both rotational directions show the same result, approaching the inflow conditions.
Figure 2h represents the situation for a backing wind. Only ϕ(z) (Eq. 10) and, therefore, the flow component v_{f}(z) (Eq. 8) change sign in both the top and bottom half of the rotor. The vortex component ${v}_{v}(z\pm R/\mathrm{2},{x}_{\text{down}})$ is not inflow dependent. Therefore, the wake behavior of ccr (cr) in the case of backing wind is comparable to cr (ccr) under veering inflow, resulting in a decrease (intensification) of v_{f}(z) in the wake, following Eq. (17) with a “−” in Eq. (10).
This analysis shows a rotationaldirectiondependent downwind behavior of the spanwise flow component in the case of $\frac{\partial {v}_{f}}{\partial z}$ ≠ 0. The superposition of the Rankine vortex with a veering inflow (and likewise the backing wind) has three impacts. The veering inflow is determined by the geostrophic wind u_{g} and the directional shear ds over the rotor height. The vortex component is determined by the rotational velocity ω of the rotor. The impact of u_{g}, ds, and ω on the expected mean behavior of the spanwise wake component is presented in Fig. 2 for low values of the parameters in the left column (panels a, d, g) and high values in the right column (panels c, f, i), whereas Fig. 2e represents the veering case for moderate parameter values.
A decrease in u_{g} (Fig. 2a) or ds (Fig. 2d) and likewise an increase in u_{g} (Fig. 2c) or ds (Fig. 2f) impacts the mean value of the spanwise wake field. When decreasing the atmospheric parameter values (Fig. 2a and d), the values of ${v}_{f}({z}_{\mathrm{h}}\pm R/\mathrm{2})$ also decrease, leading to a downwind shift in the signchanging point of v(z,x_{down}) = 0 (compare Fig. 2a, d to Fig. 2e). A further decrease in ${v}_{f}(z\pm R/\mathrm{2})$ approaching ${v}_{f}(z\pm R/\mathrm{2})$ = 0 of the nonveering inflow case results in Fig. 2b. An increase in u_{g} (Fig. 2c) or ds (Fig. 2f) results in an increase in ${v}_{f}(z\pm R/\mathrm{2})$ and an upward shift in the signchanging point. If the atmospheric parameter values increase (decrease), the difference in the slope of the spanwise component between cr and ccr also increases (Fig. 2c and f; decreases in Fig. 2a and d). Likewise, the slope in the case of cr increases for high values (Fig. 2c and f in comparison to low values Fig. 2a and d), whereas in the case of ccr, the slope decreases for larger values of the inflow parameters (Fig. 2c and f vs. Fig. 2e vs. Fig. 2a and d). This behavior can be interpreted as an increase in the difference in the wake between cr and ccr if the atmospheric parameters increase.
The rotational velocity ω controls the magnitude of the spanwise vortex component. A decrease in ω (Fig. 2g) and likewise an increase (Fig. 2i) also influence the mean value of the spanwise wake field, especially in the near wake. A decrease in ω decreases v(z,x_{down}) directly behind the rotor (Fig. 2g), whereas an increase results in an increase in v(z,x_{down}) in the near wake (Fig. 2i). Larger values of ω lead to a less rapid wake recovery in the near wake.
As a veering wind in the NH is comparable to a backing wind in the SH (following the definition via Eq. 10), all panels in Fig. 2 are also valid for the SH with red lines representing ccr_SH, blue lines representing cr_SH, dashed lines referring to the top rotor part, and solid lines referring to the bottom rotor part.
4.1 Veering vs. no veering inflow
The analysis of the preceding section predicts a rotationaldirection impact on the spanwise velocity component v (and likewise the vertical component w) in the wake under veering (or backing) inflow, whereas the wake characteristics in the case of no veer are independent of the rotational direction of the rotor. This rotationaldirection impact is investigated by LESs with veering and no veering inflow, with the simulation CR, CCR, CR_NV, and CCR_NV conducted with the parameters as listed in Table 1. The interactions between the wake rotation and the inflow are embodied in Fig. 3 in the cross stream and vertical velocities at x = 3 D. The first two columns represent CR and CCR in the case of veering inflow, whereas the last two columns correspond to no wind veer in the incoming flow field. The top row (Fig. 3a–d) represents the vectors (v, w). The evolution of v and w is represented in the second row (Fig. 3e–h) for v and in the third row (Fig. 3i–l) for w. In the case of no wind veer, the sign of v is opposite in the upper and the lower rotor half for CR_NV and CCR_NV (Fig. 3g and h), as predicted by the analysis (Eq. 17). The same is valid for the sign of w (Fig. 3k and l). The numerical model shows a clockwiserotating wake in the case of CCR_NV (while looking from upwind towards downwind; Fig. 1) and a counterclockwiserotating wake in the case of CR_NV.
Under veering inflow, the simulated wake rotates clockwise in the case of CCR (Fig. 3b) and counterclockwise in the case of CR (Fig. 3a), similar to the noveer case (Fig. 3d and c). However, in comparison to the noveer case, the strength of rotation differs and is much more pronounced in the case of CCR (Fig. 3b) in comparison to CR (Fig. 3a). This rotation arises from the spanwise velocity component, as the vertical velocity (Fig. 3i and j) is comparable to the nonveering cases (Fig. 3k and l). The positive and negative perturbations in v have the same positive and negative patterns in CR (Fig. 3e vs. Fig. 3g) and CCR (Fig. 3f vs. Fig. 3h) as in CR_NV and CCR_NV in the corresponding rotor sector at x = 3 D, although with smaller $\mid \phantom{\rule{0.125em}{0ex}}v\phantom{\rule{0.125em}{0ex}}\mid $ values in the upper and lower rotor sector in the case of CR and larger $\mid \phantom{\rule{0.125em}{0ex}}v\phantom{\rule{0.125em}{0ex}}\mid $ values in the case of CCR. This simulated amplification of the spanwise flow component in the case of CCR (Fig 3f) and weakening up to a reversion of the sign in the wake region at x = 3 D in the case of CR (Fig. 3e) is in agreement with the predictions of the analysis (Eq. 17) and Fig. 2.
A downwind distance of x = 3 D is visualized in Fig. 3 because a significant spanwise vortex impact on the spanwise flow component in the wake can be expected. In the following, special emphasis is placed at x = 7 D, which is often considered a typical downwind distance for a hypothetical waked wind turbine in numerical simulation studies (e.g., Gaumond et al., 2014; Abkar et al., 2016). At x = 7 D, the vortex impact is much smaller compared to x = 3 D (Fig. 3), resulting in an increase in the impact of the atmospheric flow.
As the rotational direction has a significant impact on the spanwise flow component at x = 3 D (Fig. 3), an impact on the streamwise flow component is also expected. The numerical results for the streamwise velocity component are presented for veering (CR, CCR) and nonveering (CR_NV, CCR_NV) inflow by x–y cross sections of the streamwise velocity in the top half of the rotor disc at z = 125 m (Fig. 4a–d), at hub height at z = 100 m (Fig. 4e–h), and in the bottom half of the rotor disc at z = 75 m (Fig. 4i–l).
The effect of wind veer on the streamwise velocity component of clockwiserotating wind turbines is investigated by comparing CR to CR_NV, in Fig. 4a vs. Fig. 4c at z = 125 m, Fig. 4e vs. Fig. 4g at z = 100 m, and Fig. 4i vs. Fig. 4k at z = 75 m. Inflow veer causes a more rapid wake recovery at all heights, based on comparison of the velocity deficit contours. Because enhanced $\frac{\partial {v}_{f}}{\partial z}$ ≠ 0 in the case of veering wind, it provides a source of resolved turbulence resulting in higher entrainment in comparison to the noveer case. Further, inflow wind veer causes wake deflection in relation to a vertical plan through the nacelle at y = 0 in both the top half (Fig. 4a vs. Fig. 4c) and the bottom half (Fig. 4i vs. Fig. 4k) of the rotor disc. The wake in the veered simulation CR is deflected towards the right (y > 0 D; towards the left is y < 0 D) in the upper (lower) rotor part (Fig. 4a and i). In the nonveered simulation CR_NV, the wake is only slightly deflected towards the left in the toptip sector (Fig. 4c) and towards the right in the bottomtip sector (Fig. 4k). This effect is caused by the rotation of the rotor, which transports highermomentum air counterclockwise, resulting in a wake deflection to the left at z = 125 m (Fig. 4c). Consequently, the opposite situation prevails at z = 75 m (Fig. 4k). As the inflow veer contribution to wake deflection is much larger compared to the effect of a clockwiserotating rotor, the wake deflection changes from the left in CR_NV (Fig. 4c) to the right in CR (Fig. 4a) in the upper rotor half and vice versa in the lower rotor half.
As a next step, the rotationaldirection impact in the nonveered simulations CR_NV and CCR_NV is investigated (Fig. 4c vs. Fig. 4d, Fig. 4g vs. Fig. 4h, and Fig. 4k vs. Fig. 4l). The impact of the rotational direction on the wake is limited to the wake deflection differences at the upper (Fig. 4c and d) and the lower (Fig. 4k and l) rotor height, which are nearly axissymmetric to y = 0 D and result from the rotational direction of the rotor. These differences in the nonveered simulations agree with results of Vermeer et al. (2003), Shen et al. (2007), Sanderse (2009), Kumar et al. (2013), Hu et al. (2013), Yuan et al. (2014), Mühle et al. (2017), and Englberger et al. (2020).
The rotationaldirection impact on the wake structure under veering inflow is investigated by a comparison of CCR to CR (Fig. 4b vs. Fig. 4a, f vs. Fig. 4e, and Fig. 4j vs. Fig. 4i). In CCR, the wake recovers more rapidly (Fig. 4f vs. Fig. 4e) and the wake deflection angle is larger (Fig. 4b vs. Fig. 4a and Fig. 4j vs. Fig. 4i) in comparison to CR. Further, the wake width is larger in the spanwise direction in CCR in comparison to CR (Fig. 4b vs. Fig. 4a, Fig. 4f vs. Fig. 4e, and Fig. 4 vs. Fig. 4i).
The differences in the spanwise wake width and the wake deflection angle are investigated in more detail with the y–z cross sections at x = 7 D in Fig. 5 for veering inflow (CR in Fig. 5a, CCR in Fig. 5b) and no wind veer (CR_NV in Fig. 5c, CCR_NV in Fig. 5d) with both rotational directions of the actuator. In the case of no veering inflow, the simulated wake at x = 7 D retains the shape of the rotor (Fig. 5c). In the case of a veering inflow, however, the wake in the lower rotor half is shifted to the left and in the upper rotor half to the right (Fig. 5a). The striking difference between veering and nonveering inflow simulations in combination with a clockwiserotating actuator corresponds to the inflow profile (Eqs. 6 and 8), where a veering inflow is characterized by a wind component from right to left for z < 100 m and from left to right for z > 100 m, whereas the spanwise inflow velocity is zero in the case of no veer at all rotor heights. The skewed wake structure under veering inflow resembles those of the simulations of Abkar and PortéAgel (2016), Vollmer et al. (2017), Bromm et al. (2017), Churchfield and Sirnivas (2018), and Englberger and Dörnbrack (2018a).
Further, we compare the differences between a clockwise and a counterclockwiserotating actuator for nonveering and veering inflow. In the case of no wind veer, the simulated wake structures of CCR_NV (Fig. 5d) and CR_NV (Fig. 5c) show no striking difference. In the case of veering inflow, however, the skewed wake structure differs in CR and CCR (Fig. 5a and b). Whereas the wake is elliptical in CR, this shape is stretched in the rotor region in CCR. This difference in shape explains the difference in the spanwise wake width at hub height (Fig. 4f) and also in the lower (Fig. 4j) and the upper (Fig. 4b) rotor part. The wake structure outside the rotor region also differs between Fig. 5a and b. Due to the elongation of the elliptical structure in CCR in the rotor region (Fig. 4b) and approximately the same vertical wake extension in CCR and CR, the wake deflection angle increases in the case of CCR (Fig. 5b vs. Fig. 5a), as shown in the lower rotor half in Fig. 4j vs. Fig. 4i and also in the upper rotor half in Fig. 4b vs. Fig. 4a.
A quantitative description of the streamwise velocity differences is presented in Fig. 6 at the downwind position x = 7 D. Figure 6 represents the vertical profiles at y = 0 D in a and spanwise profiles of u at z = 75 m in b, at z = 100 m in c, and at z = 125 m in d for both rotational directions CR and CCR. The heights correspond to Fig. 4. Figure 6e–h represent the nonveering inflow simulations CR_NV and CCR_NV. Whereas the vertical and spanwise profiles of CR and CCR in the case of no inflow veer (_NV) are almost overlapping (Fig. 6e–h), there is a difference in the case of veering inflow (Fig. 6a–d). Firstly, the streamwise wake elongation difference in Fig. 4f vs. Fig. 4e is represented by larger u values in the lower and the upper rotor half in the case of CCR in Fig. 6b and d. The larger wake deflection angle in CCR in comparison to CR (Fig. 4b vs. Fig. 4a and Fig. 4j vs. Fig. 4i) is represented by a larger spanwise distance of the minimum of u from y = 0 D in the case of CCR in the lower (Fig. 6b) and the upper (Fig. 6d) rotor half. This spanwise difference in u_{min} is accompanied by larger u values in the case of CCR for y < $\mathrm{1}/\mathrm{2}$ D in the lower rotor part (Fig. 6b) and for y < 1/2 D in the upper rotor part (Fig. 6d). Secondly, the difference in the spanwise wake width is represented at all three heights by a larger ΔL_{y} with smaller u values in CCR in the outermost region of the left and the right sectors in comparison to CR (Fig. 6b–d).
As a final step, the difference in the wake is summarized by the 30 min timeaveraged and rotorareaaveraged streamwise velocity $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$. Figure 7a represents the difference between clockwise and counterclockwiserotating rotors for a veering inflow and in the case of no wind veer from x = 4 D to 10 D. At x = 7 D, $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ is 0.24 m s^{−1} larger in the counterclockwiserotating rotor simulation CCR in comparison to CR, whereas there is no difference between CCR_NV and CR_NV. According to Fig. 6a–d, these larger $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ values in the case of CCR result from larger u values in the upper and lower sector related to the larger wake deflection angle in the case of CCR, which compensates for the larger u values in the outer region of the left and right sectors resulting from a larger spanwise wake width in the case of CCR.
The previous investigations show a striking dependence of the rotational direction of the rotor on the wake under veering inflow, which is qualitatively well explained by the analysis. A schematic illustration of the deceleration or even reversion of the spanwise flow if a clockwiserotating rotor CR interacts with a veering wind is presented in Fig. 8a. The amplification of the spanwise flow in the case of a counterclockwiserotating rotor CCR interacting with veering inflow is presented in Fig. 8b.
4.2 Veering wind vs. backing wind
According to the analytical results (Fig. 2e vs. Fig. 2h), the spanwise component v in the wake is expected to be comparable for a clockwiserotating rotor in veering inflow and a counterclockwiserotating rotor in backing inflow, as well as for a clockwiserotating rotor in backing inflow and a counterclockwiserotating rotor in veering inflow. The wake characteristics resulting from a backing wind with both rotational directions are investigated in the simulations CR_b and CCR_b and compared to the veering wind cases CR and CCR in Fig. 9. The parameters applied in the corresponding simulations are listed in Table 1.
The behavior in the upper and the lower rotor part in Fig. 9 can directly be compared after mirroring at y = 0 D, an effect resulting from the opposite sign of the directional shear and Δϕ in Eqs. (9) and (10). A strong similarity is prevalent in the streamwise velocity component at hub height (Fig. 9f, g and e, h), in the lower rotor half (Fig. 9i, l and j, k) as well as in the upper rotor half (Fig. 9a, d and b, c). The more rapid wake recovery and the larger spanwise wake width for CCR and CR_b in comparison to CR and CCR_b are present at all rotor heights. The larger wake deflection angle in the upper and the lower rotor half in CCR (Fig. 9b and j) is also comparable to CR_b (Fig. 9c and k), whereas the smaller wake deflection angle in CR (Fig. 9a and i) is comparable to CCR_b (Fig. 9d and l).
The qualitative comparison in Fig. 6i–l shows the differences from Fig. 9 between CR_b and CCR_b in the case of streamwise wake elongation, spanwise wake width, and the wake deflection angle. Further, comparing the backingwind situation (Fig. 6i–l) to the veeringwind situation (Fig. 6a–d), CR_b corresponds to CCR and CCR_b to CR in the vertical profiles (Fig. 6i vs. Fig. 6a) and at hub height (Fig. 6k vs. Fig. 6c). After mirroring at y = 0 D, CR_b corresponds to CCR and likewise CCR_b to CR in the lower and the upper rotor part (compare Fig. 6j to Fig. 6b and Fig. 6l to Fig. 6d).
Expressing the differences between a backing and a veering wind from both rotational directions of the rotor by the quantity $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ in Fig. 7a, the $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ values are 0.24 m s^{−1} larger if a backing wind (CR_b) interacts with a clockwiserotating rotor in comparison to a veering wind (CR). Similarly, the $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ values are larger if a backing wind interacts with a counterclockwiserotating rotor (CCR_b). Therefore, $\mathrm{\Delta}\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ is the same for CR and CCR_b and likewise for CCR and CR_b.
The northern hemispheric results of CR and CCR are comparable to southern hemispheric CR_b and CCR_b situations, whereas the northern hemispheric results of CR_b and CCR_b correspond to CR and CCR in the SH. The schematic illustration of a backing wind interacting with both rotational directions is presented in Fig. 8c and d with an amplification of the spanwise wind component in the case of a backing wind and a clockwiserotating rotor CR_b (Fig. 8d) and a weakening or reversion in the case of a counterclockwiserotating rotor CCR_b (Fig. 8c).
4.3 Wind speed
Wind speed may also affect the veering inflow (Eq. 8 via Eq. 6), modifying the spanwise velocity component. There is no significant impact of u_{g} = (6, 10, 14) m s^{−1} on the wake elongation, the spanwise wake width, and the wake deflection angle between clockwise and counterclockwiserotating actuators. Therefore, the contour plots are not shown. Only a qualitative comparison is presented in Fig. 6m–p for u_{g} = 6 m s^{−1} and in Fig. 6q–t for u_{g} = 14 m s^{−1}. The occurrence of a wake width as well as the wake deflection angle difference between clockwise and counterclockwiserotating actuators from the reference case u_{g} = 10 m s^{−1} (Fig. 6b–d) is independent of u_{g}. Only for smaller velocity values (u_{g} = 6 m s^{−1}; Fig. 6o) are the hub height differences less pronounced.
The similarity of the vertical and spanwise profiles for all geostrophicwind values in Fig. 6 results in no remarkable difference in $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ in Fig. 7b between clockwise and counterclockwiserotating simulations. Independent of u_{g}, the values of $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ are slightly larger in the case of counterclockwiserotating simulations in comparison to clockwise ones, except the difference in $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ between clockwise and counterclockwiserotating simulations increases for decreasing u_{g}.
4.4 Directional shear
The directional shear is the second contributing parameter of the veering inflow (Eq. 8), modifying the spanwise velocity component resulting from analysis. The impact of all five directionalshear values from Table 1 on the wake is investigated at hub height (Fig. 10) and in the upper (Fig. 11) and lower (Fig. 12) rotor half. In the clockwiserotating as well as the counterclockwiserotating actuator simulations (Figs. 10–12) the wake recovers more rapidly if directional shear increases. A larger directional shear represents a larger turbulence source due to an increase in $\frac{\partial {v}_{f}}{\partial z}$. Therefore, the simulations with larger directionalshear values result in higher entrainment rates and a more rapid wake recovery. Our simulated dependence of the wake recovery on the amount of wind veer for clockwiserotating simulations is comparable to the numerical results in Fig. 11 of Bhaganagar and Debnath (2014).
The magnitude of directional shear affects the wake elongation in dependence on the rotor direction but not to the same extent in clockwiserotating simulations in comparison to counterclockwiserotating ones. The wake elongation in CR_ds4 is much longer in comparison to CCR_ds4 (panel a vs. panel b in Figs. 10–12). It is still larger in CR in comparison to CCR (panel c. vs. panel d in Figs. 10–12). A further increase in the directional shear finally results in a similar wake recovery of CR_ds12 and CCR_ds12 (panel e vs. panel f in Figs. 10–12) and a slightly more rapid wake recovery of CR_ds16 in comparison to CCR_ds16 (panel g vs. panel h in Figs. 10–12). Comparing the very strong directionalshear cases CR_ds20 and CCR_ds20 (panel i vs. panel j in Figs. 10–12), the wake recovery is significantly faster in CR_ds20. Further, the difference between CR_ds4 and CR is larger in comparison to CCR_ds4 and CCR (panels a, c and panels b, d in Figs. 10–12). This trend continues for increasing directional shear.
Another difference between clockwise and counterclockwiserotating actuators is the spanwise wake width (Fig. 10c vs. d). The impact of the directional shear on the spanwise wake width at hub height results in an increase in the difference in the spanwise wake width between a clockwise and a counterclockwiserotating simulation (Fig. 10a, c, e, g, i vs. Fig. 10b, d, f, h, j). In addition, the wake deflection angle increases for increasing values of the directional shear. The difference in larger wake deflection angles in the case of a counterclockwiserotating actuator in the upper and the lower rotor part is also prevalent for all directionalshear values (right column of Figs. 11 and 12, panels b, d, f, h, j). Further, large values of the directional shear in combination with a counterclockwiserotating actuator leads to a breakup of the wake (panels h and j in Figs. 11 and 12). This erosion could be related to high spanwisevelocity values in the case of CCR_ds16 and CCR_ds20 due to amplification of the inflow and the vortex spanwise component, which is not the case in CR_ds16 and CR_ds20 (panels g and i in Figs. 12 and 11).
For a quantitative investigation of the directionalshear impact on the differences in the wake between clockwise and counterclockwiserotating actuators, vertical and horizontal profiles at x = 7 D for all five cases of different directionalshear values are presented in Fig. 13. Considering the vertical profile through y = 0 D (left column of Fig. 13), the vertical wake extension decreases if the directional shear increases, as the wake deflection is influenced by the incoming wind direction at each height (Churchfield and Sirnivas, 2018; Tomaszewski et al., 2018; Bodini et al., 2017; Englberger and Lundquist, 2020). This dependency of wake veer on wind veer is also represented at z = 75 m (Fig. 13b, f, j, n, r) and at z = 125 m (Fig. 13d, h, l, p, t), where the wake deflection angle is additionally influenced by the rotational direction of the actuator. The wake deflection angle is larger if the actuator rotates counterclockwise, independent of the values of directional shear.
The directionalshear impact on the spanwise wake width is investigated via the profiles of Fig. 13. Especially at hub height (Fig. 13c, g, k, o, s), the wake width decreases if the directional shear increases. This effect can be related to the increase in skewness in the wake for an increasing directional shear. Comparing clockwise and counterclockwiserotating actuators, the spanwise wake width is larger in the case of a counterclockwiserotating actuator, independent of the directionalshear value.
Considering the rotoraveraged values $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ in Fig. 7c, the rotoraveraged wind speeds are larger for a weak wind veer in the counterclockwiserotating actuator simulations (CCR_ds4) in comparison to the clockwiserotating ones (CR_ds4). As the wind veer increases, the difference in $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ between clockwise and counterclockwiserotating disc simulations decreases. In the case of a moderate to strong wind veer, $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ is only slightly larger for the clockwiserotating rotor CR_ds12. Approaching an even higher directional shear in the strong and very strong wind shear cases, this difference between CR_ds16 and CCR_ds16 and likewise between CR_ds20 and CCR_ds20 increases, whereas now the rotoraveraged wind speeds are larger for clockwiserotating actuators in comparison to counterclockwise ones.
Independent of the directional shear, the streamwisevelocity values at x = 7 D are larger in the case of a counterclockwiserotating actuator in the top and the bottom sector (Fig. 13 second and fourth column). The difference between counterclockwise and clockwise rotation increases in the radial direction away from the nacelle (not shown). This is related to the larger wake deflection angle in the counterclockwiserotating case in comparison to the clockwise case. Also independent of the directional shear, the streamwisevelocity values are larger in the case of a clockwiserotating actuator in the left and right sectors (Fig. 13 third row). This is an effect of the narrower wake width in the case of a clockwiserotating actuator. If the directional shear is small, the larger u values of counterclockwiserotating actuators in the top and bottom sectors are compensating for the larger u values in the case of a clockwiserotating actuator in the right and left sectors. If the rotational direction is very high, the opposite is the case.
4.5 Rotational frequency
The rotational frequency contributes to the windturbine forces in Eq. (1) and modifies the spanwise velocity component (Eq. 17). The wake impact of the four rotationalfrequency values from Table 1 is presented at hub height (Fig. 14), at z = 125 m (Fig. 15), and at z = 75 m (Fig. 16). In comparison to the impact of a change in the atmospheric parameters, which was mainly limited to the far wake, the rotational frequency also significantly impacts the near wake. An increase in the rotational frequency results in a larger minimum value of the velocity deficit (panels g, h vs. panels a, b in Figs. 15 and 16) and a less rapid wake recovery (panels g, h vs. panels a, b in Figs. 15 and 16). As the rotational frequency increases, the wake structure differs more between clockwise and counterclockwiserotating actuators. The difference in the spanwise wake width increases for an increasing rotational frequency at all heights. Further, an increase in the rotational frequency results in a slightly larger downwind wake extension in the case of a clockwiserotating actuator at all heights. In the case of a counterclockwiserotating actuator, however, an increase in the rotational frequency results in a similar downwind wake extension of the velocity deficit (Eq. 14).
A wakesplitting pattern exists in the upper (Fig. 15f and h), as well as in the lower (Fig. 16f and h) rotor part for large rotationalfrequency values. The pattern is similar to the breakup of the wake for large directionalshear values interacting with a counterclockwiserotating actuator in CCR_ds16 and CCR_ds20 (panels h and j in Figs. 11 and 12). The occurrence of the pattern in combination with high rotationalfrequency values could also be related to a very large spanwise flow component, now resulting from a large contribution of the vortex. An additional simulation (not shown) with u_{g} = 10 m s^{−1}, ds = 0.20${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$, and Ω = 0.23 s^{−1} reinforces the slitting pattern of Figs. 12 and 11h and j, supporting our assumption as to why the splitting occurs only for counterclockwiserotating actuators. Further, a similar but less distinctive wakesplitting pattern for a counterclockwiserotating actuator was observed in the veeraffected lower rotor half with ds = 0.28${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$ and Ω = 0.12 s^{−1}; see Fig. 10f by Englberger et al. (2020).
For a quantitative investigation of the rotationalfrequency impact on the wake differences between clockwise and counterclockwiserotating rotors, the vertical and horizontal profiles at x = 7 D are presented for all four cases in Fig. 17. Considering the vertical and spanwise profiles at z = 100 m (Fig. 17 first and third column), the rotationaldirection decrease in u results from larger windturbine forces due to an increase in Ω. The difference in the wake defection angle (Fig. 17 second and fourth column) and in the spanwise wake width (Fig. 17 second, third, and fourth column) between clockwise and counterclockwiserotating discs increases for increasing Ω. An increase in Ω further results in two u minima in the lower (Fig. 17j and n) and the upper (Fig. 17l and p) rotor half and a larger decrease in u approaching r = R at hub height (Fig. 17k and o) in the counterclockwiserotating simulations.
The increase in u in the lower and upper sector compensates for the decrease in the left and right sector for increasing Ω in the case of counterclockwiserotating discs (Fig. 17), resulting in larger values of $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ in Fig. 7d for all counterclockwiserotating simulations. The $\stackrel{\mathrm{\u203e}}{{u}_{\mathrm{A}}}$ difference between clockwise and counterclockwiserotating actuators increases for an increasing rotational frequency, which is related to the splitting of the wake.
The idealized numerical simulations investigated the impact of the rotational direction of the actuator in combination with veering inflow, no wind veer, and backing inflow on the wake of a single wind turbine. The parameter study investigated the streamwise dependency of the wake on wind speed, directional shear, and rotational frequency. For a comparison of the simulated results with the expected results from analysis in Fig. 2, Fig. 18 is plotted for 90^{∘} bottom and top sectors in the range 0 m < r ≤ 50 m (Fig. 1). Figure 2 represents the results at y = 0 D and z = 75 m as the bottom part of the rotor disc and at y = 0 D and z = 125 m as the top of the rotor disc. The general structure (slope, signchanging point), however, is independent of the vertical location in the analysis. Only the magnitude of the spanwise inflow v_{f} at x_{down} > x_{ξ} and of the spanwise vortex component v_{v} at x_{down} < x_{ξ} are affected in Fig. 2, as v_{f} is height dependent and asymmetric to the rotor center and v_{v} has a radial dependency. Therefore, the panels of Fig. 18 are directly comparable to those of Fig. 2 regarding the difference in Δu_{g}, Δds, and ΔΩ between low, moderate, and highvalue cases.
Comparing the nonveering simulations CR_NV and CCR_NV (Fig. 18b) to the analysis prediction (Fig. 2b), v > 0 in the top sector and v < 0 in the bottom sector both with a clockwiserotating rotor. In the case of a counterclockwiserotating rotor, v in the top sector corresponds to v in the bottom sector of a clockwiserotating simulation and vice versa. Only the downwind slope for x_{down} < x_{ξ} is much smaller. This results from a different radial distribution and a smaller absolute value of the windturbine forces applied in the numerical simulations (Eq. 1) in comparison to the Rankine vortex applied in the theoretical analysis (Eqs. 15 and 16). Further, the smaller slope in the numerical simulations can be related to the resolved turbulence and the resulting wake recovery.
Comparing the simulation with moderateveering inflow CR and CCR (Fig. 18e) to the analysis predictions (Fig. 2e), the acceleration in the case of a counterclockwiserotating rotor and the weakening in the case of a clockwiserotating rotor up to x_{down} ≈ 10 D are prevalent. The smaller slope values are for the same reason as in the nonveering case. The different slope values between the simulations and the analysis predictions simply result in an upwind shift in the signchanging location in the wake. An increase in Ω (Fig. 18i) approaches the structure predicted by analysis (steeper slope, flow reversion behind the rotor, downwind shift in the signchanging point) of Fig. 2b. The values of Ω applied in the BEM method to calculate the windturbine forces (Eq. 1) are the same as the values of ω applied in the Rankine vortex (Eq. 15). Due to the differences in the calculation of the spanwise flow field between the BEM method and the Rankine vortex, the nearwake absolute values are not comparable in Figs. 2 and 18. Considering the backing inflow simulation (Fig. 18h) and exchanging clockwise and counterclockwise and likewise top and bottom, it corresponds to the veering inflow in Fig. 18e, which is predicted by Fig. 2h and e.
The analytic model predicts the same impact of the geostrophic wind and the directional shear on the spanwise wake structure (Fig. 2a vs. Fig. 2d and Fig. 2c vs. Fig. 2f). The general structure at x_{down} < x_{ξ} is comparable in the numerical simulations with u_{g} = 6 m s^{−1} and ds = 0.08${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$ (Fig. 18a) and u_{g} = 10 m s^{−1} and ds = 0.04${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$ (Fig. 18d) and likewise with u_{g} = 14 m s^{−1} and ds = 0.08${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$ (Fig. 18c) and u_{g} = 10 m s^{−1} and ds = 0.12${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$ (Fig. 18f). Minor differences exist; e.g., compare the difference between clockwise and counterclockwiserotating simulations at x = 7 D in Fig. 18a and d. The larger difference between CR_ds4 and CCR_ds4 in Fig. 18d can be related to a decrease in $\frac{\partial {v}_{f}}{\partial z}$ and a smaller amount of resolved turbulence generated by the inflow with ds = 0.04${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$ (Fig. 18d) in comparison to ds = 0.08${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$ (Fig. 18a), whereas the change in u_{g} has no influence on $\frac{\partial {v}_{f}}{\partial z}$. Changing the wind speed and the directional shear has a significant impact on $\mid v\mid $ further downwind at x_{down} > x_{ξ}; e.g., v_{20D} in Fig. 18f ≈ 2 ⋅ v_{20D} in Fig. 18d. This corresponds to the differences between Fig. 2f and d at x_{down} > x_{ξ}.
Changing the rotational frequency of the vortex has its largest impact on the spanwise wake velocity directly behind the rotor, whereas an increase in the spanwise vortex component v_{v} results in a larger amplification of $\mid v\mid $ in the case of a counterclockwiserotating rotor and a larger weakening of $\mid v\mid $ in the case of a clockwiserotating rotor (Fig. 18g, e, and i). This behavior corresponds to the nearwake differences in Fig. 2g, e, and i. For large enough values of the rotational frequency, the spanwise wake component reverses sign in the simulation CR_Ωh directly behind the rotor (Fig. 18i).
The comparison of the simulation results with the analysis predictions can be summarized as follows:

The simulated amplification or weakening or reversion of the spanwise inflow wind component in the wake follows the theoretical analysis in the case of no wind veer, a veering inflow, and a backing inflow for clockwise and counterclockwiserotating discs. It can be understood and described by the superposition of the rotational flow induced by the disc and the vertical shear of the incoming wind.

The agreement of the simulation results with the analysis predictions proves that the impact of the rotational direction on the spanwise wake field is determined by the mean values of the inflow wind field and influenced by the resulting turbulence.

The inflow parameters (wind speed and directional shear) and the rotation rate of the rotor are two counteracting processes. The individual magnitudes determine the differences in the spanwise wake component between clockwise and counterclockwiserotating actuators. The difference increases (decreases) for increasing (decreasing) values of u_{g}, ds, and Ω.
The rotationaldirection impact on the spanwise velocity component in the wake also modifies the streamwise flow component. The streamwise velocity components in the wake are shown in Fig. 19 for the same sectors and simulations as in Fig. 18. The larger values in the top sector in comparison to the bottom sector result from the heightdependent streamwise velocity (Eq. 6) with larger values in the upper rotor half. In the case of no wind direction change with height (_NV, Fig. 19b), the $\stackrel{\mathrm{\u203e}}{u}$ values are independent of the rotational direction of the actuator. In the case of a veering (backing) inflow (Fig. 19e and b), the streamwise wake velocity is larger in the case of CCR (CR_b) in comparison to CR (CCR_b) in both sectors. The parameters under veering inflow impact the difference $\mathrm{\Delta}\stackrel{\mathrm{\u203e}}{u}$ between counterclockwise and clockwiserotating actuators in the corresponding sector (Fig. 19a, d, g and c, f, i). An increase in u_{g}, ds, or Ω results in larger $\mathrm{\Delta}\stackrel{\mathrm{\u203e}}{u}$ values (Fig. 19c, f, i), whereas smaller parameter values decrease $\mathrm{\Delta}\stackrel{\mathrm{\u203e}}{u}$ (Fig. 19a, d, g). The difference $\mathrm{\Delta}\stackrel{\mathrm{\u203e}}{u}$ in the left and right sectors is much less distinct in comparison to the top and bottom sectors (not shown).
The rotationaldirection impact on the wake can be summarized as follows:

A rotationaldirection impact on the streamwise velocity components in the wake exists only in the case of a veering (or backing) inflow.

In the case of veering inflow in the NH the spanwise wake width (ΔL_{y}) as well as the wake deflection angle (Δϵ) is larger in the case of a counterclockwiserotating actuator CCR (Fig. 20). This behavior is independent of the magnitude of the parameters.

Increasing the magnitudes of the directional shear ds or the rotation rate Ω increases the spanwise wake width difference (ΔL_{y}) and the wake deflection angle difference (Δϵ) between a counterclockwise and a clockwiserotating actuator. The impact of the geostrophic wind u_{g} is much less pronounced.

An increase in u_{g} or likewise a decrease in Ω results in a more rapid wake recovery. Increasing ds, there is no wake difference between a clockwise and a counterclockwiserotating actuator for a specific directional shear ds_{c}. If ds < ds_{c}, the streamwise velocity is larger in the case of a counterclockwiserotating rotor, whereas for ds > ds_{c}, the streamwise velocity is larger in the case of a clockwiserotating rotor. Approaching smaller or larger values of the directional shear, $\mathrm{\Delta}\stackrel{\mathrm{\u203e}}{u}$ increases between a clockwise and a counterclockwiserotating actuator.
We investigate the impact of the rotational direction on the wake of a wind turbine for veering and backing inflow conditions, as well as in the case of no wind veer in both hemispheres, using idealized LESs in comparison to a simple analytic model. In addition, the impact of the geostrophic wind and the directional shear as well as the impact of the rotational frequency on the wake differences between clockwise and counterclockwiserotating wind turbines was investigated in the case of veering inflow.
The rotational direction of a wind turbine has only a minor impact on the wake in the case of no wind veer. This result from the numerical experiments is consistent with previous investigations by Vermeer et al. (2003), Shen et al. (2007), Sanderse (2009), Kumar et al. (2013), Hu et al. (2013), Yuan et al. (2014), Mühle et al. (2017), and Englberger et al. (2020). An inflow without wind veer is a typical daytime situation and also occurs in the evening transition of the diurnal boundary layer evolution, when the flow is still influenced by daytime turbulence.
In the case of veering or backing inflow, however, the wake characteristics (streamwise wake elongation, spanwise wake width, wake deflection angle) depend significantly on the rotational direction. Veering and backing inflow are characteristic nighttime situations of the boundary layer flow if no other processes such as topographically induced circulations or largescale weather systems prevent the establishment of an SBL regime. Veer within the windturbine rotor layer has been observed with several field campaigns with towers and lidars (Walter et al., 2009; Sanchez Gomez and Lundquist, 2020; Bodini et al., 2019, 2020), and veer throughout the boundary layer has been observed globally using radiosonde datasets (Lindvall and Svensson, 2019).
Under veering inflow in the NH (backing inflow in the SH), the spanwise wake width and the wake deflection angle are larger for a counterclockwiserotating (clockwiserotating) actuator in comparison to a clockwiserotating one. An increase (decrease) in the directional shear in the atmospheric flow or in the rotational frequency of the rotor increases (decreases) the differences in the spanwise wake width and the wake deflection angle. The wind speed does not impact these wake characteristics significantly. In locations with veering inflow in the NH (backing inflow in the SH) and directionalshear values ds < ds_{c} with 0.12${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$ < ds_{c} < 0.16${}^{\circ}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$, the streamwise velocity is larger in the case of a counterclockwiserotating (clockwiserotating) rotor. These differences apply to the wake ranging from x = 4 D to at least x = 10 D downwind. For less common higher values of the directional shear ds > ds_{c}, the streamwise velocity is larger in the case of a clockwiserotating (counterclockwiserotating) rotor in the NH (SH).
Different operating conditions (e.g., yaw control) of upwind turbines are already applied to mitigate downwind impacts in wind parks (Fleming et al., 2019). This work suggests that counterclockwiserotating blades in the case of veering inflow and clockwiserotating blades in the case of backing inflow in the NH (and vice versa in the SH) could have benefits as well. The wake deflection angle becomes larger if the spanwise flow component is amplified by the vortex induced by the rotating wind turbine. This process occurs independent of the magnitude of the parameter values applied in the numerical simulations.
As the numerical results of this study arise from an idealized parameter study employing specific assumptions, they have limitations. For example, the turbulent perturbations applied in the numerical simulations are incorporated via a simple turbulence parametrization. The imposed turbulence parameters were retrieved from a precursor LES and no realtime wind and potentialtemperature profiles are applied. Particularly, the impact of the rotational direction on a windturbine wake under veering (or backing) inflow resulted from basic analytical predictions and was compared with the numerical model. However, the impact of rotational directions has never been measured, as no counterclockwiserotating wind turbines currently exist. Despite the limitations of this numerical study, the simple analysis as well as the idealized parameter study shows a consistent and clear impact of the rotational direction of a wind turbine on the wake flow during conditions for which the wind direction turns with height.
To explore a more comprehensive assessment of the wake impact, further investigations would be interesting: the investigation of the nonlinearity of the interaction process, numerical simulations applying the turbulence of an SBL precursor simulation for different strengths of stratification and directional shears, or even considering a lowlevel jet at the rotor height. Topography could influence the wake dynamic explored here. We have assessed the wake of an individual turbine, but these results could be extended to a large farm in which the presence of upwind turbines could affect turbulence intensity, which probably affects the magnitude. However, an important point will be to prove the theoretically predicted effect resulting from superposition of inflow veer with the vortex component on the wake with measurements.
Finally, the overall assessment of the impact of these results depends on the frequency of occurrence of veering inflow. Only limited sets of longterm observations provide an assessment of the frequency of veering (Walter et al., 2009; Sanchez Gomez and Lundquist, 2020; Bodini et al., 2019, 2020) in the windturbine rotor layer. The global climatology of veer throughout the atmospheric boundary layer based on radiosonde data (Lindvall and Svensson, 2019) suggests that veer occurs broadly in midlatitudes and polar regions, but further investigation is required to assess if that boundary layer veer broadly affects wind energy generation.
The turbulence parametrization of Englberger and Dörnbrack (2018a) is applied in the simulations of this work. The main part is conducted with α = 0.3, α_{u} = 0.15, α_{v} = 0.24, and α_{w} = 0.13. These values are nighttime representations following Table 1 of Englberger and Dörnbrack (2018b). Figure A1 presents the reference CR and CCR windturbine simulation applied in this work at z = 125 m in a and b, at z = 100 m in e and f, and at z = 125 m at i and j. Further, simulation results α = 0.3 and ${\mathit{\alpha}}_{{i}^{\ast},j,k}$ = 0 are presented at z = 125 m in c and d, at z = 100 m in g and h, and at z = 125 m in k and l. Panels a and b are the reference simulations CR and CCR with veering inflow. Panels c and d correspond to α = 0.3, panels e and f to α = 0.5, and panels g and h to α = 0.7, with ${\mathit{\alpha}}_{{i}^{\ast},j,k}$ = 0 in all three cases. The streamwise velocity at hub height, as well as in the lower and the upper rotor half, shows similar characteristics of the nearwake velocity deficit maximum, the streamwise wake elongation, the spanwise wake width, and the wake deflection angle. Only the strength of occurrence of these wake characteristics depends on the turbulent intensity, which is larger in the case of ${\mathit{\alpha}}_{{i}^{\ast},j,k}$ = 0. This reinforces the assumption that wake characteristic differences depend on the mean wind profile, which is the same in all simulations of Fig. A1 and is not an effect of the applied turbulence parametrization.
Currently, the data are not publicly available.
All authors conceived the idea. AE performed the simulations and prepared the manuscript with contributions from both coauthors.
The authors declare that they have no conflict of interest.
The authors gratefully acknowledge the Gauss Center for Supercomputing e. V. (http://www.gausscentre.eu, last access: 19 November 2020) for funding this project by providing computing time on the GCS supercomputer SuperMUC at the Leibniz Supercomputing Center (LRZ; http://www.lrz.de, last access: 19 November 2020).
The article processing charges for this openaccess publication were covered by a Research Center of the Helmholtz Association.
This paper was edited by Sandrine Aubrun and reviewed by Paul van der Laan and two anonymous referees.
Abkar, M. and PortéAgel, F.: Influence of the Coriolis force on the structure and evolution of wind turbine wakes, Physical Review Fluids, 1, 063701, https://doi.org/10.1103/PhysRevFluids.1.063701, 2016. a
Abkar, M., Sharifi, A., and PortéAgel, F.: Wake flow in a wind farm during a diurnal cycle, J. Turbul., 17, 420–441, https://doi.org/10.1080/14685248.2015.1127379, 2016. a
Bak, C., Zahle, F., Bitsche, R., Kim, T., Yde, A., Henriksen, L. C., Hansen, M. H., Blasques, J. P. A. A., Gaunaa, M., and Natarajan, A.: The DTU 10MW reference wind turbine, in: Danish Wind Power Research 2013, DTU Trinity, Fredericia, Denmark from 27 May 2013–28 May 2013. a
Bhaganagar, K. and Debnath, M.: Implications of Stably Stratified Atmospheric Boundary Layer Turbulence on the NearWake Structure of Wind Turbines, Energies, 7, 5740–5763, https://doi.org/10.3390/en7095740, 2014. a
Bodini, N., Zardi, D., and Lundquist, J. K.: Threedimensional structure of wind turbine wakes as measured by scanning lidar, Atmos. Meas. Tech., 10, 2881–2896, https://doi.org/10.5194/amt1028812017, 2017. a
Bodini, N., Lundquist, J. K., and Kirincich, A.: US East Coast Lidar Measurements Show Offshore Wind Turbines Will Encounter Very Low Atmospheric Turbulence, Geophys. Res. Lett., 46, 5582–5591, https://doi.org/10.1029/2019GL082636, 2019. a, b, c, d
Bodini, N., Lundquist, J., and Kirincich, A.: Offshore Wind Turbines Will Encounter Very Low Atmospheric Turbulence, J. Phys. Conf. Ser., 1452, National Renewable Energy Lab. (NREL), Golden, CO, USA, 012023, https://doi.org/10.1088/17426596/1452/1/012023, 2020. a, b, c, d, e, f, g
Bromm, M., Vollmer, L., and Kühn, M.: Numerical investigation of wind turbine wake development in directionally sheared inflow, Wind Energy, 20, 381–395, 2017. a
Churchfield, M. J. and Sirnivas, S.: On the effects of wind turbine wake skew caused by wind veer, in: 2018 Wind Energy Symposium, p. 0755, AIAA American Institute of Aeronautics and Astronautics, Reston, 2018. a, b
Englberger, A. and Dörnbrack, A.: Impact of Neutral BoundaryLayer Turbulence on WindTurbine Wakes: A Numerical Modelling Study, Bound.Lay. Meteorol., 162, 427–449, https://doi.org/10.1007/s105460160208z, 2017. a, b
Englberger, A. and Dörnbrack, A.: Impact of the diurnal cycle of the atmospheric boundary layer on windturbine wakes: a numerical modelling study, Bound.Lay. Meteorol., 166, 423–448, https://doi.org/10.1007/s1054601703093, 2018a. a, b, c
Englberger, A. and Dörnbrack, A.: A Numerically Efficient Parametrization of Turbulent WindTurbine Flows for Different Thermal Stratifications, Bound.Lay. Meteorol., 169, 505–536, https://doi.org/10.1007/s105460180377z, 2018b. a, b, c
Englberger, A. and Lundquist, J. K.: How does inflow veer affect the veer of a windturbine wake?, J. Phys. Conf. Ser., 1452, 012068, https://doi.org/10.1088/17426596/1452/1/012068, 2020. a, b
Englberger, A., Dörnbrack, A., and Lundquist, J. K.: Does the rotational direction of a wind turbine impact the wake in a stably stratified atmospheric boundary layer?, Wind Energ. Sci., 5, 1359–1374, https://doi.org/10.5194/wes513592020, 2020. a, b, c, d, e, f, g, h
Fleming, P., King, J., Dykes, K., Simley, E., Roadman, J., Scholbrock, A., Murphy, P., Lundquist, J. K., Moriarty, P., Fleming, K., van Dam, J., Bay, C., Mudafort, R., Lopez, H., Skopek, J., Scott, M., Ryan, B., Guernsey, C., and Brake, D.: Initial results from a field campaign of wake steering applied at a commercial wind farm – Part 1, Wind Energ. Sci., 4, 273–285, https://doi.org/10.5194/wes42732019, 2019. a
Fröhlich, J.: Large Eddy Simulation turbulenter Strömungen, Teubner Verlag/GWV Fachverlage GmbH, Wiesbaden, 414 pp., 2006. a
Gaumond, M., Réthoré, P.E., Ott, S., Pena, A., Bechmann, A., and Hansen, K. S.: Evaluation of the wind direction uncertainty and its impact on wake modeling at the Horns Rev offshore wind farm, Wind Energy, 17, 1169–1178, 2014. a
Grinstein, F. F., Margolin, L. G., and Rider, W. J.: Implicit Large Eddy Simulation, Cambridge University Press, New York, USA, 546 pp., 2007. a
Hu, H., Yuan, W., Ozbay, A., and Tian, W.: An experimental investigation on the effects of turbine rotation directions on the wake interference of wind turbines, in: 51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Grapevine, Texas, 7–10 January 2017, p. 607, 2013. a, b, c
Kumar, P. S., Abraham, A., Bensingh, R. J., and Ilangovan, S.: Computational and experimental analysis of a counterrotating wind turbine system, NISCAIRCSIR, India, 300–306, available at: http://hdl.handle.net/123456789/17441 (last access: 19 November 2020), 2013. a, b, c
Lindvall, J. and Svensson, G.: Wind turning in the atmospheric boundary layer over land, Q. J. Roy. Meteor. Soc., 145, 3074–3088, 2019. a, b, c, d
Maegaard, P., Krenz, A., and Palz, W.: Wind power for the world: the rise of modern wind energy, Jenny Stanford, Boca Raton, Florida, 2013. a, b
Mühle, F., Adaramola, M. S., and Sætran, L.: The effect of rotational direction on the wake of a wind turbine rotor – a comparison study of aligned coand counter rotating turbine arrays, Energ Proced., 137, 238–245, https://doi.org/10.1016/j.egypro.2017.10.346, 2017. a, b, c
Prusa, J. M., Smolarkiewicz, P. K., and Wyszogrodzki, A. A.: EULAG, a computational model for multiscale flows, Comput. Fluids, 37, 1193–1207, https://doi.org/10.1016/j.compfluid.2007.12.001, 2008. a, b
Rhodes, M. E. and Lundquist, J. K.: The effect of windturbine wakes on summertime US Midwest atmospheric wind profiles as observed with groundbased doppler lidar, Bound.Lay. Meteorol., 149, 85–103, https://doi.org/10.1007/s105460139834x, 2013. a
Sanchez Gomez, M. and Lundquist, J. K.: The effect of wind direction shear on turbine performance in a wind farm in central Iowa, Wind Energ. Sci., 5, 125–139, https://doi.org/10.5194/wes51252020, 2020. a, b, c, d, e, f, g
Sanderse, B.: Aerodynamics of wind turbine wakes, Energy Research Center of the Netherlands (ECN), ECNE–09016, Petten, the Netherlands, Tech. Rep., 5, 153, 2009. a, b, c
Shen, W. Z., Zakkam, V. A. K., Sørensen, J. N., and Appa, K.: Analysis of counterrotating wind turbines, J. Phys. Conf. Ser., 75, 012003, https://doi.org/10.1088/17426596/75/1/012003, 2007. a, b, c
Smolarkiewicz, P. K. and Margolin, L. G.: MPDATA: A FiniteDifference Solver for Geophysical Flows, J. Comput. Phys., 140, 459–480, https://doi.org/10.1006/jcph.1998.5901, 1998. a
Smolarkiewicz, P. K., Sharman, R., Weil, J., Perry, S. G., Heist, D., and Bowker, G.: Building resolving largeeddy simulations and comparison with wind tunnel experiments, J. Comput. Phys., 227, 633–653, https://doi.org/10.1016/j.jcp.2007.08.005, 2007. a
Stull, R. B.: An Introduction of Boundary Layer Meteorology, Kluwer Academic, Dordrecht, 1988. a, b
Tomaszewski, J. M., Lundquist, J. K., Churchfield, M. J., and Moriarty, P. J.: Do wind turbines pose roll hazards to light aircraft?, Wind Energ. Sci., 3, 833–843, https://doi.org/10.5194/wes38332018, 2018. a
VaselBeHagh, A. and Archer, C. L.: Wind farms with counterrotating wind turbines, Sustainable Energy Technologies and Assessments, 24, 19–30, 2017. a, b
Vermeer, L., Sørensen, J. N., and Crespo, A.: Wind turbine wake aerodynamics, Prog. Aerosp. Sci., 39, 467–510, https://doi.org/10.1016/S03760421(03)000782, 2003. a, b, c
Vollmer, L., Lee, J. C., Steinfeld, G., and Lundquist, J.: A wind turbine wake in changing atmospheric conditions: LES and lidar measurements, J. Phys. Conf. Ser., 854, 012050, https://doi.org/10.1088/17426596/75/1/012003, 2017. a
Walter, K., Weiss, C. C., Swift, A. H., Chapman, J., and Kelley, N. D.: Speed and direction shear in the stable nocturnal boundary layer, J. Sol. Energ.T. ASME, 131, 011013, https://doi.org/10.1115/1.3035818, 2009. a, b, c, d, e, f, g, h
Yuan, W., Tian, W., Ozbay, A., and Hu, H.: An experimental study on the effects of relative rotation direction on the wake interferences among tandem wind turbines, Sci. China Phys. Mech., 57, 935–949, 2014. a, b, c
Zhang, W., Markfort, C. D., and PortéAgel, F.: Nearwake flow structure downwind of a wind turbine in a turbulent boundary layer, Exp. Fluids, 52, 1219–1235, https://doi.org/10.1007/s0034801112508, 2012. a, b, c
 Abstract
 Copyright statement
 Introduction
 Numerical model framework
 Theoretical analysis
 Idealized simulations – rotationaldirection impact on the wake
 Comparison to analytic model
 Conclusions
 Appendix A: Turbulence parametrization
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References
 Abstract
 Copyright statement
 Introduction
 Numerical model framework
 Theoretical analysis
 Idealized simulations – rotationaldirection impact on the wake
 Comparison to analytic model
 Conclusions
 Appendix A: Turbulence parametrization
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References