Research article 05 Mar 2019
Research article  05 Mar 2019
The aerodynamics of the curled wake: a simplified model in view of flow control
 National Renewable Energy Laboratory, Golden, CO, USA
 National Renewable Energy Laboratory, Golden, CO, USA
Correspondence: Luis A. MartínezTossas (luis.martinez@nrel.gov)
Hide author detailsCorrespondence: Luis A. MartínezTossas (luis.martinez@nrel.gov)
When a wind turbine is yawed, the shape of the wake changes and a curled wake profile is generated. The curled wake has drawn a lot of interest because of its aerodynamic complexity and applicability to wind farm controls. The main mechanism for the creation of the curled wake has been identified in the literature as a collection of vortices that are shed from the rotor plane when the turbine is yawed. This work extends that idea by using aerodynamic concepts to develop a controloriented model for the curled wake based on approximations to the Navier–Stokes equations. The model is tested and compared to timeaveraged results from largeeddy simulations using actuator disk and line models. The model is able to capture the curling mechanism for a turbine under uniform inflow and in the case of a neutral atmospheric boundary layer. The model is then incorporated to the FLOw Redirection and Induction in Steady State (FLORIS) framework and provides good agreement with power predictions for cases with two and three turbines in a row.
This work was authored by the National Renewable Energy Laboratory, operated by Alliance for Sustainable Energy, LLC, for the U.S. Department of Energy (DOE) under Contract No. DEAC3608GO28308. Funding provided by the U.S. 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 U.S. Government. The U.S. Government retains and the publisher, by accepting the article for publication, acknowledges that the U.S. Government retains a nonexclusive, paidup, irrevocable, worldwide license to publish or reproduce the published form of this work, or allow others to do so, for U.S. Government purposes.
A curled wake is a phenomenon observed in the wake of a wind turbine when the turbine is yawed relative to the freestream velocity. When a wind turbine is yawed, the wake is not only deflected in a direction opposite to the yaw angle, but its shape changes. The mechanism behind this effect has drawn attention not only from fluid dynamicists because of the interesting physics phenomena happening in the wake but also from the wind turbine controls community who intends to use it to control wind farm flows (Fleming et al., 2017).
It has been shown that wake steering (redirection of the wake through yaw misalignment) can lead to an increase in power production of wind turbine arrays (Adaramola and Krogstad, 2011; Park et al., 2013; Gebraad et al., 2016). Previous studies have used largeeddy simulations (LESs) and analytical tools to show the effects of yawing in the redirection of the wake (Jiménez et al., 2010; Bastankhah and PortéAgel, 2016; Shapiro et al., 2018). Wind tunnel experiments have also been used to study the wake of a wind turbine in yaw (Medici and Alfredsson, 2006; Bartl et al., 2018). The curled wake mechanism, in the context of wind turbine wakes, was first identified by Howland et al. (2016) during a porous disk experiment, and in LESs using an actuator disk model (ADM) and actuator line model (ALM). This mechanism was described by a pair of counterrotating vortices that are shed from the top and bottom of the rotor when the rotor is yawed. Further, these vortices move the wake to the side and create a curled wake shape. This mechanism was later confirmed and elaborated by the work of Bastankhah and PortéAgel (2016) by performing experiments using particle image velocimetry of a scaled wind turbine and by doing a theoretical analysis using potential flow theory. Shapiro et al. (2018) show that the vorticity distribution shed from the rotor because of yaw has an elliptic shape as opposed to only a pair of counterrotating vortices. The curled wake has also been observed in LESs with different atmospheric stabilities (Vollmer et al., 2016). Berdowski et al. (2018) were able to reproduce curled wake profiles by using a vortex method. Fleming et al. (2017) show that the generated vortices affect the performance of wake steering and motivate the development of engineering models (like the one in this paper), which include wake curling physics. Controllers based on such models would pursue different, and likely more effective, wind farm control strategies.
FLOw Redirection and Induction in Steady State (FLORIS) is a software framework used for wind plant performance optimization (Gebraad et al., 2016; Annoni et al., 2018). A wake model is used in FLORIS to compute the effect of wind turbine wakes on downstream turbines. Different models can be used inside of FLORIS to compute the wind turbine wakes, including the Jensen and Gaussian wake models (Bastankhah and PortéAgel, 2016; Jensen, 1983). A review of controloriented models can be found in Annoni et al. (2018). The model proposed in this work is a new wake model and it has been incorporated to the FLORIS framework. The model is used to compute the wake from each turbine in a wind farm. In the curled wake model, the V and W velocity components generated by each turbine are superposed linearly. After computing the individual wakes, they are added using the sum of squares method (Katic et al., 1987).
In this work, we describe the aerodynamics of the curled wake, and propose a new, simple, and computationally efficient model for wake deficit, based on a linearized version of the Navier–Stokes equations with approximations. The model is tested and compared to LESs using actuator disk and line models.
Here, we develop a simplified model of the wake deficit considering the aerodynamics of the curled wake. We start by writing the Reynoldsaveraged Navier–Stokes streamwise momentum equation for an incompressible flow:
where u, v, and w are the timeaveraged velocity components in the streamwise, spanwise, and wallnormal directions. p is the timeaveraged pressure, ρ is the fluid density, and ν_{eff} is the effective (turbulent and molecular) viscosity. The flow field is now decomposed into a base solution and a perturbation about this solution. The base solution includes what we consider to be the main effects that convect the wake, including, but not limited to,

the streamwise velocity profile,

the rotational velocity from the shed vortices caused by yawing, and

the rotational velocity due to the blade rotation.
The velocity components can be expanded as
where U, V, and W are the streamwise, spanwise, and wallnormal velocity components from the base solution, and u^{′}, v^{′}, and w^{′} are the perturbation velocities. U, V, and W represent a base flow that is responsible for convecting the wakes. And u^{′}, v^{′}, and w^{′} are the perturbation velocities about the base solution which represent the wake deficits. We use a reference frame where x is the streamwise direction, y is the spanwise direction, and z is the wallnormal direction. Assuming that the inflow is fully developed, the base solution is only a function of the spanwise and wallnormal directions.
Linearizing the Euler momentum equation for the streamwise component leads to
For this initial test of the model, the pressure gradient is assumed to be zero, which is a valid approximation in the far wake but not necessarily near the rotor. Also, we assume that the streamwise velocity, U, is not a function of the streamwise (x) and spanwise (y) coordinates. We also assume that the viscous force from the boundary layer is balanced by its pressure gradient. With these simplifications, we are left with the equation
Equation (4) describes the downstream evolution of the wake deficit, u^{′}. We note that the only velocity component being solved is the streamwise component, u^{′}. The base solution for the flow (U, V, W) includes the spanwise effects from different features, such as rotation and curl. Because these effects are the main drivers for wake deformation, the spanwise perturbations, v^{′} and w^{′}, are assumed to be zero. This assumption reduces the model to a single partial differential equation. Equation (4) will be solved numerically in the next section. The model does not use the continuity equation so mass conservation is not strictly enforced for the perturbation velocities.
2.1 Curled wake
The curled wake effect is added to the model by adding a distribution of counterrotating vortices to the base flow solution. Figure 1 shows a schematic of the rotor and a collection of vortices being shed at the rotor plane. The superposition of all the vortices leads to a spanwise velocity distribution that creates the curled wake shape.
Each vortex is described as a Lamb–Oseen vortex with a tangential velocity distribution given by
where u_{t} is the tangential component of the velocity, r is the distance from the vortex core, Γ is the circulation strength, and σ determines the width of the vortex core. The circulation related to the strength of each vortex is still unknown. We assume that the problem is symmetric and that the entire circulation leaves the disk through the shed vortices. In the current implementation, the vortices are assumed not to decay as they are convected downstream; however, it would be possible to incorporate this effect into the model. It is important to note that the vorticity shed has an elliptic distribution according to the shape of the rotor (Shapiro et al., 2018). The total circulation is then Γ=ρDU_{∞}F_{L}, where ρ is the fluid density, D is the rotor diameter, U_{∞} is the inflow velocity at hub height, and F_{L} is the force perpendicular to the inflow wind. The force component perpendicular to the inflow wind generates circulation. If we use the definition of the thrust coefficient, it is now possible to redefine the strength of the circulation as a function of thrust coefficient and yaw angle (Shapiro et al., 2018) as
where C_{T} is the thrust coefficient for the given inflow velocity, and γ is the yaw angle. The discrete elliptic distribution of shed vortices added to match the shape of the disk is
where i is the index denoting each of the vortices distributed on a line between the top and bottom of the rotor diameter, N is the total number of vortices, and the coordinates y_{i} and z_{i} are centered at the location of the shed vortex. The size of the vortex core is set to $\mathit{\sigma}=D/\mathrm{5}$. The choice of σ intends to represent a realistic vortex core size for an elliptic distribution. This value represents $\mathit{\sigma}/D\sim \mathrm{0.2}$, which is on the order of the optimal size for flow over an airfoil (MartínezTossas et al., 2017). It is possible to use other values; however, in the simulations presented, this value is a good compromise between a physical width and numerical stability. The strength of each vortex is associated with the elliptical distribution by
where r_{i} is the radial location of the ith vortex in the rotor coordinate system and ${\mathrm{\Gamma}}_{\mathrm{0}}=\mathrm{4}/\mathit{\pi}\mathrm{\Gamma}$ is used to ensure that the total amount of circulation Γ is conserved. This means that each vortex has a unique amount of circulation. When the distribution of circulation in Eq. (9) is integrated, the total circulation is obtained. The results use N=200 discrete vortices, which has shown that this has little effect on the solutions, and values around N=20 provide the same results. Note that, when using just the two tip vortices at the top and bottom of the rotor, the shape of the curled wake does not match the simulation results. For this reason, an elliptic distribution of shed vorticity, which matches the shape of the rotor, should be used.
2.2 Wake rotation
It is important to include wake rotation in the model, because the rotation will move the wake in a preferred direction. Wake rotation is taken into account by adding a tangential velocity distribution that is caused by the rotation inside the rotor area. The tangential induction factor is defined as
where a is the induction factor based on the thrust coefficient from standard actuator disk theory, R is the rotor radius, r is the radial distance from the center of the rotor, and λ is the tip speed ratio (Burton et al., 2002). From this equation, the tangential component of the velocity u_{t} is a singular vortex with 1∕r behavior:
A Lamb–Oseen vortex is used to desingularize the behavior near the center of the rotor. The circulation strength for the wake rotation vortex based on Eq. (10) is now
We assume that the wake rotation vortex does not decay or deform as it moves downstream. This is not necessarily true, as turbulence mixing will decrease the wake rotation. However, the present model does not diffuse the spanwise velocity components and some of the errors in the model can be attributed to this. The current implementation uses a Lamb–Oseen vortex with a core size $\mathit{\sigma}=D/\mathrm{5}$ to eliminate numerical instabilities caused by high velocities near the vortex core, but other values could be used.
2.3 Atmospheric boundary layer
The atmospheric boundary layer can be specified as part of the background flow. A profile including streamwise and spanwise velocity components can be specified. The streamwise profile is described by using a power law:
where U_{h} is the velocity at hub height, z_{h} is the hub height, and α is the shear exponent. Also, it is possible to add a profile for the spanwise velocity component to take veer into account. In order to avoid numerical instabilities, the minimum wind speed is set to 20 % of U_{h}. This happens close to the bottom wall where the results do not affect the solution in the wake.
2.4 Turbulence modeling
The turbulent viscosity in Eq. (4) is determined by using a mixing length model. The turbulent viscosity in the atmospheric boundary layer is dependent on a mixing length, ℓ_{m}, and a velocity gradient (Pope, 2001). The mixing length for flows in the atmospheric boundary layer is defined by
where κ is the von Kàrmàn constant, z is the distance from the wall, and λ=15 m is the value reached by ℓ_{m} in the free atmosphere (Blackadar, 1962; Sun, 2011). Now the turbulent viscosity is given by
where $\frac{\text{d}u}{\text{d}z}$ is the streamwise velocity gradient in the wallnormal direction.
2.5 Ground effect
The presence of the ground will have an effect on the shed vortices. The ground effect is incorporated by applying a symmetry boundary condition at the ground (Bastankhah and PortéAgel, 2016). This is done by using Eqs. (7) and (8), with the y and z coordinates placed below the ground and inverting the sign of the circulation. This condition has a more dominant effect on the vortices close to the ground as they interact with the boundary.
2.6 Superposition of solutions
Superposing all the effects mentioned earlier leads to a base flow that includes all the features presented. The linearized equation allows us to add features by superposing them onto the velocity components. Notice that, in this implementation of the model, the base solutions are a function of only the spanwise directions, y and z, and there is no dependency on the streamwise coordinate, x. However, dependency on the streamwise direction could also be included in the base solution. The vortices do induce motion on each other and the spanwise component of the momentum equation should be used. However, these motions are smaller than the streamwise motions and solving one equation as opposed to three reduces computational cost significantly.
2.7 Initial and boundary conditions
The initial condition for the perturbation velocity, u^{′}, is specified as the yawed disk projected onto a plane normal to the streamwise direction. This shape represents the shape of the wake downstream right after the rotor. The exact shape of the wake near the rotor is much more complicated and is not taken into account in the present model. The initial profile is set as a uniform distribution of wake deficit (${u}^{\prime}=\mathrm{2}a{U}_{\mathrm{\infty}}$) inside the rotor projected area, where a is the induction. This step function is smoothed using a Gaussian filter to avoid numerical oscillations in the spanwise directions. The lateral boundaries are set to zero perturbation (${u}^{\prime}=\mathrm{0}$) because there is no wake in that region.
It is now possible to discretize Eq. (4) and solve it numerically. Because the timederivative and pressuregradient terms were dropped, the equation is parabolic, and it can be solved as a marching problem. The equation is solved by starting from an initial condition at the rotor plane and marching downstream. This is done by using a firstorder upwind discretization for the streamwise derivative and secondorder finite differencing for the spanwise derivatives. The discrete equation is now
where i, j, and k are the indices denoting the grid points in the x, y, and z coordinates, ∇^{2} is the wallnormal and spanwise components in the Laplacian operator, and ν_{eff} is the effective viscosity. Because of the explicit form of the equation, it is possible to include the nonlinear term ${u}^{\prime}\frac{\partial {u}^{\prime}}{\partial x}$, which is not present in Eq. (4). We note that this term has a small effect on the solution, and not including it provides similar results. The resolution used is on the order of 30–40 grid points per diameter. The typical domain size needed to capture the wake is on the order of 3D in the spanwise and wallnormal directions, and as long as needed for the downstream location. The computational expense of the algorithm without any optimization is small (∼1–3 s) and can be used to generate curled wake profiles quickly.
3.1 Numerical stability
The proposed numerical method uses a forwardtime, centeredspace method (Hoffman and Frankel, 2001). This algorithm is explicit, meaning that it can become unstable for certain conditions. Using Eq. (16), we can establish the numerical criteria for stability from the forwardtime, centeredspace method (Hoffman and Frankel, 2001). Equation (17) shows a guideline for the stability requirement of the algorithm:
This stability criterion is based on a twodimensional equation and the equation we are solving is threedimensional. However, after testing various conditions, this criterion served as a good guideline for the threedimensional version of the equation. After some algebraic manipulation, it is possible to show that the maximum grid spacing in the streamwise direction, Δx, is independent of the spanwise grid resolution. Equation (17) can be written as
After testing the model with several grid spacings, it was found that a resolution on the order of $D/\mathrm{\Delta}\sim \mathrm{30}$–40 provided converged results for the model without numerical oscillations. In the case of laminar inflow, an arbitrary viscosity needs to be added to stabilize the numerical solution. After experimenting with different values, an effective viscosity based on a Reynolds number
proved to be sufficient to stabilize the solution.
In this section, we compare the proposed model to LESs with an actuator disk/line model. Different simulations are used to test the proposed model: (1) a simulation of an isolated turbine using the ADM under uniform inflow, (2) a simulation of an isolated turbine using the ALM under uniform inflow, and (3) a simulation of a turbine using the ALM inside the atmospheric boundary layer under neutral stability conditions.
4.1 Actuator disk/line model under uniform inflow
Here, we compare the results from the model to LESs of a wind turbine under uniform inflow of a turbine using an actuator disk/line model under uniform inflow from Howland et al. (2016). The yaw angle for these simulations is γ=30^{∘}. First, the model is compared to a simulation using an actuator disk without rotation. Figure 2 shows downstream planes with contours of streamwise velocity normalized by the inflow velocity, for the case of an ADM without rotation (thrust only). The streamlines are based on the crossstream components of velocity. The overall shape of the curled wake is well captured by the model. The overall streamline shapes between the LES and proposed wake model are similar. The pair of counterrotating vortices is clearly visible in both cases; however, the LES computes streamlines with a more complex shape than the simpler model is capable of capturing. The resulting effect is that, in both cases, the wake deficit cross sections are deformed and curled in a similar fashion. The model does not contain the tower, which is present in the simulation; however, this has a small effect on the wake.
Figure 3 shows downstream velocity contours for the case of an LES using an actuator line model under uniform inflow (Howland et al., 2016) and the proposed model including curl and rotation. The difference between the model implementation in this case and the case of the actuator disk model is that rotation of the wake has been added. Again, the streamlines are similar in both cases, but the LES produces more complex patterns. Further, the resultant wake deficit deformation is also similar in both cases, with more deficit remaining at the top of the wake. In this case, asymmetry is observed with respect to the centerline across the y axis. This asymmetry is caused by the wake rotation induced by the torque applied to the fluid by the rotor. Interestingly, the combination of curl and rotation pushes most of the deficit in a preferred direction. In this case, it is pushed upward in the positive z and negative y directions. This insight can be used to steer the wake accordingly. We also observe that, in the LES, the wake diffuses faster than in the model. This is because the model does not yet take into account the turbulence generated by the wake.
Figure 4 shows the axial velocity along a horizontal line at hub height for different downstream locations from the simulations in Figs. 2 and 3. There is good agreement between the simulations and the model (in general), although some differences can be observed. Near the edges of the wake there is an acceleration in the LES, which is caused by the blockage effect, and the model is not able to capture this. In the case of the ALM, the main difference can be attributed to the different initial condition. The model assumes a step function, which is different from the wake resolved by the LES. However, the general shape of the wake and its deviation to the sides are wellcaptured by the model.
Figure 5 shows the axial velocity along a vertical line that passes through the center of the rotor for different downstream locations from the simulations in Figs. 2 and 3. Good agreement can be observed from the profiles between the model and LES. The general behavior of the profiles is wellcaptured by the model. In the case of the ADM, differences in the near wake are present due to the inclusion of the tower model. The profiles from the model have sharper gradients near the edges, showing that a turbulence model should be present to incorporate the diffusion due to turbulent mixing.
4.2 Largeeddy simulations using the actuator line model in the atmospheric boundary layer
The framework presented can easily be extended by adding more features. As an example, we present a comparison of the model with an LES of a wind turbine inside a neutral atmospheric boundary layer with a yaw angle, γ=20^{∘}. The largeeddy simulation was performed using the Simulator fOr Wind Farm Applications (SOWFA) from the National Renewable Energy Laboratory (Churchfield and Lee, 2012). We also include the Gaussian wake model from Bastankhah and PortéAgel (2016) in the comparison.
To add the effects of the atmosphere to the curled wake model, a vertical profile of velocity in the streamwise direction is added to the base solution. Also, a linear spanwise velocity component is added to the base solution to take veer into account, although this had little effect on the results presented. The veer profile was chosen as a linear profile that matched the inflow from the LES results. Figure 6 shows the inflow wind vertical profile three diameters upstream of the rotor for the LES and the one specified from the model, and the resulting turbulent viscosity as a function of height. A power law with a shear exponent of α=0.15 was chosen in the model to match the LES inflow condition. The turbulence intensity at hub height from the LES is 10 %.
Figure 7 shows the mean streamwise velocity contours for the LES of a wind turbine inside a neutral atmospheric boundary layer and results from the model. In general, there is good agreement between the model and the simulation. We can see that the main difference comes from the wake in the LES diffusing more than in the model. This is expected because the turbulence model does not take into account the turbulence generated by the turbine wake, only the turbulence caused by the velocity gradients in the atmospheric boundary layer. There are also differences resulting from the lateral motion of the vortices. The present model does not take into account the convection of the vortices. This is shown in Fig. 7, where the top and bottom vortices stay at the same place when using the model. In reality, these vortices are convected to the sides, as shown in the LES.
Figure 8 presents velocity along horizontal and vertical lines passing through the center of the rotor comparing the model presented, the Gaussian wake model from Bastankhah and PortéAgel (2016) and largeeddy simulations. The Gaussian model follows the implementation by Bastankhah and PortéAgel (2016) with the wake growth rate of ${k}_{x}={k}_{y}=\mathrm{0.022}$. There are differences present in the proposed model, which we attribute to the simplifications of the proposed model and the turbulence model. However, there is good agreement in terms of nearwake predictions. From the vertical profiles, we can see that further downstream the agreement between the curled wake model and the LES deteriorates. This is because the turbulent diffusion from the model does not provide enough dissipation, and because the vortices do not decay nor are they convected to the side. This means that the vortices will convect the wake deficit, thereby providing unrealistic stronger wake deficits further downstream.
In the LES and proposed model, the curled wake shape is produced in the near wake. As the wake evolves, the turbulence diffuses the curled wake shape, and eventually it becomes more similar to a Gaussian wake, as observed by Bastankhah and PortéAgel (2016). In contrast to the LES and curled wake model, the Gaussian model predicts a very symmetric wake at every downstream location and is not able to predict the complex shape of the curled wake. Further downstream, a selfsimilar profile seems to be a good description of the wake, since the turbulence has diffused most of the structures caused by the curling mechanism.
It is difficult to track a wake centerline in the curled wake model and the LES. The curled wake is characterized by a complex threedimensional structure and a wake center is not really descriptive of this mechanism, especially in the near wake. Figure 9 shows lateral displacement of the center of the wake for LESs, the curled wake model, the Gaussian model, and the model from Shapiro et al. (2018). The curled wake and LES cases compute this displacement by averaging a collection of tracers around the center of the rotor inside a radius of 0.2D. We see that all the models agree relatively well in predicting the lateral displacement. We argue that the displacement is not a proper measure of the wake because it cannot track the nonsymmetric and more complex shape of the curled wake. The tracer is not only displaced laterally but also in the vertical direction. To properly represent the curled wake and its displacement, a more robust and threedimensional model, such as LES and/or the curled wake model, should be used.
We now test the proposed model inside the FLORIS framework (Gebraad et al., 2016; Annoni et al., 2018). We compare the new curled wake model to the twodimensional Gaussian steering model from Bastankhah and PortéAgel (2016). In the Gaussian steering model (Bastankhah and PortéAgel, 2016), the wake deflection due to yaw misalignment of turbines is defined by doing budget analysis on the Reynoldsaveraged Navier–Stokes equations. In the curled wake model, the wake steering is computed by solving a linearized version of the Navier–Stokes equations.
First, we run a case of two turbines aligned with the flow and the upstream turbine is yawed by 25^{∘}. Figure 10 shows the streamwise velocity profiles for a FLORIS simulation with the curled wake model and the Gaussian model from the Bastankhah and PortéAgel (2016) model. We observe that, in the curled wake model, the wake of the second turbine is affected by the curl of the first turbine. It was observed by Fleming et al. (2017) that the curled wake mechanism does affect the wake of the second turbine, but current yaw steering models are not able to take this effect into account. This is one of the main advantages of the curled wake model, because secondary wake steering has recently been observed and can be used to optimize wind farm performance (Fleming et al., 2017).
Now, we present results for three turbines aligned with the upstream turbine yawed 25^{∘}. The curled wake model predicts deflections up to the third turbine's wake. This approach addresses previous concerns about models not being able to capture wake deflection from downstream turbines (Fleming et al., 2017).
Table 1 shows a comparison between the power predictions and performance of the Gaussian and curled wake models and simulations performed in SOWFA. In the case of two turbines, both models agree very well in terms of power predictions. However, we notice that the power predictions from the curled wake overpredict results from LESs in the case of three aligned turbines. This outcome is expected because the vortices resulting from curl do not decay as they travel downstream. Without a decay model, the spanwise velocities from the yawed turbine would never decay. In reality, these vortices decay due to the turbulence in the atmosphere and in the wake.
The curled wake model provides improvements in predicting power gains for more than two turbines in a row. This outcome is because the vortices from the first turbine are propagated downstream. However, because the vortices do not decay in time, the power may be overpredicted.
The key differences between the model and simulations can be summarized as follows:

The vortices caused by the curl effect in the model do not change their position and do not decay. In reality, these vortices induce motion on each other and are advected by the freestream flow, which may have a lateral component.

The turbulence model does not take into account the wind turbine wake. It can only take into account the turbulence from the atmospheric boundary layer background flow. This is why the wake decays faster in the largeeddy simulations compared to the model.

The vortices in the model do not decay with downstream distance. In reality, vortices decay because of the radial diffusion of tangential momentum.

The model does not take into account all the nonlinear interactions present in the simulation. For this reason, the model is only able to capture the behavior of the larger scales, and hence, not all the details of the flow (such as the deformation of the vortices) can be captured.
The model can be further improved by taking some of these factors into account. However, the present model is able to capture the main dynamics of the curled wake with a reduced computational cost. Further improvements are part of future work.
A new model has been proposed to study the aerodynamics of the curled wake. The model solves a linearized version of the Navier–Stokes momentum equation with the curl effect added as a collection of vortices with an elliptic distribution shed from the rotor plane. The main difference between the model presented and the Gaussian models from Bastankhah and PortéAgel (2016) and Shapiro et al. (2018) is that there is no assumption on the shape of the wake, apart from its initial condition. The wake profile generated by the curled wake model is the solution to the linearized momentum equation with some assumptions. This allows the wake to take any shape, which in many cases was observed to differ significantly from a Gaussian wake.
The model has the ability to include several features of the wake including effects due to yaw (“curl”), wake rotation, a boundary layer profile, and turbulence modeling. The model has been implemented and tested to reproduce curled wake profiles. Good agreement is observed when comparing the model to largeeddy simulations of flow past a yawed turbine using an actuator disk/line model. The model was implemented and tested using the FLORIS framework. Good agreement was observed in predicting power extraction by yawing the first turbine in a row of two and three turbines. We observe that the effects of the vortices shed by a yawed turbine propagate for downstream distances longer than the separation between two turbines. This means that a yawed turbine can be used to redirect not only its own wake but the wake of other downstream turbines as well. Also, we note that the shed vortices allow for spanwise velocity components, which are vital when considering wake redirection and wind farm controls. The vortices generated are not limited to only yawing, as they can also be used for tilt and combinations of tilt and yaw. This work sets a foundation for a simplified wake steering model to be used in a more general wind farm controloriented framework. Future work consists of improving the curled wake model with emphasis on implementing a robust decay model for the vortices and comparing the model to experimental data.
This code is currently under development and not publicly available yet. Information can be obtained in the meantime from the corresponding author.
The authors declare that they have no conflict of interest.
The authors would like to acknowledge Charles Meneveau and Patrick Hawbecker
for suggestions on turbulent modeling in the atmospheric boundary layer and
Patrick Moriarty for providing wind turbine aerodynamics insights.
Simulations for the NREL code SOWFA were performed using the National
Renewable Energy Laboratory's Peregrine highperformance computing system.
Numerical implementation and plots were done using Python with NumPy,
Matplotlib, and SciPy libraries.
Edited by: Jens Nørkær Sørensen
Reviewed by: three anonymous referees
Adaramola, M. and Krogstad, P.A.: Experimental investigation of wake effects on wind turbine performance, Renew. Energ., 36, 2078–2086, 2011. a
Annoni, J., Fleming, P., Scholbrock, A., Roadman, J., Dana, S., Adcock, C., PorteAgel, F., Raach, S., Haizmann, F., and Schlipf, D.: Analysis of controloriented wake modeling tools using lidar field results, Wind Energ. Sci., 3, 819831, https://doi.org/10.5194/wes38192018, 2018. a, b, c
Bartl, J., Mühle, F., Schottler, J., Sætran, L., Peinke, J., Adaramola, M., and Hölling, M.: Wind tunnel experiments on wind turbine wakes in yaw: effects of inflow turbulence and shear, Wind Energ. Sci., 3, 329–343, https://doi.org/10.5194/wes33292018, 2018. a
Bastankhah, M. and PortéAgel, F.: Experimental and theoretical study of wind turbine wakes in yawed conditions, J. Fluid Mech., 806, 506–541, https://doi.org/10.1017/jfm.2016.595, 2016. a, b, c, d, e, f, g, h, i, j, k, l
Berdowski, T., Ferreira, C., van Zuijlen, A., and van Bussel, G.: ThreeDimensional FreeWake Vortex Simulations of an Actuator Disc in Yaw and Tilt, in: AIAA 2018 Wind Energy Symposium, 8–12 January 2018, 0513, 2018. a
Blackadar, A. K.: The vertical distribution of wind and turbulent exchange in a neutral atmosphere, J. Geophys. Res., 67, 3095–3102, 1962. a
Burton, T., Sharpe, D., Jenkins, N., and Bossanyi, E.: Aerodynamics of HorizontalAxis Wind Turbines, John Wiley and Sons Ltd, 41–172, https://doi.org/10.1002/0470846062.ch3, 2002. a
Churchfield, M. J. and Lee, S.: SOWFA – NWTC Information Portal, available at: https://nwtc.nrel.gov/SOWFA (last access: 26 February 2019), 2012. a
Fleming, P., Annoni, J., Churchfield, M., MartinezTossas, L. A., Gruchalla, K., Lawson, M., and Moriarty, P.: A simulation study demonstrating the importance of largescale trailing vortices in wake steering, Wind Energ. Sci., 3, 243–255, https://doi.org/10.5194/wes32432018, 2018. a, b, c, d, e
Gebraad, P. M. O., Teeuwisse, F. W., van Wingerden, J. W., Fleming, P. A., Ruben, S. D., Marden, J. R., and Pao, L. Y.: Wind plant power optimization through yaw control using a parametric model for wake effects – a CFD simulation study, Wind Energ., 19, 95–114, https://doi.org/10.1002/we.1822, 2016. a, b, c
Hoffman, J. D. and Frankel, S.: Numerical methods for engineers and scientists, CRC press, 2001. a, b
Howland, M. F., Bossuyt, J., MartínezTossas, L. A., Meyers, J., and Meneveau, C.: Wake structure in actuator disk models of wind turbines in yaw under uniform inflow conditions, J. Renew. Sustain. Ener., 8, 043301, https://doi.org/10.1063/1.4955091, 2016. a, b, c, d, e, f, g
Jensen, N. O.: A note on wind generator interaction, Tech. rep., Risø National Laboratory Roskilde, 1983. a
Jiménez, Á., Crespo, A., and Migoya, E.: Application of a LES technique to characterize the wake deflection of a wind turbine in yaw, Wind Energ., 13, 559–572, 2010. a
Katic, I., Højstrup, J., and Jensen, N. O.: A simple model for cluster efficiency, European Wind Energy Association, Conference and Exhibition, Rome Italy, 7–9 October 1986, 1987. a
MartínezTossas, L. A., Churchfield, M., and Meneveau, C.: Optimal smoothing length scale for actuator line models of wind turbine blades based on Gaussian body force distribution, Wind Energ., 20, 1083–1096, https://doi.org/10.1002/we.2081, 2017. a
Medici, D. and Alfredsson, P.: Measurements on a wind turbine wake: 3D effects and bluff body vortex shedding, Wind Energ., 9, 219–236, 2006. a
Park, J., Kwon, S., and Law, K. H.: Wind farm power maximization based on a cooperative static game approach, in: Proc. SPIE, vol. 8688, https://doi.org/10.1117/12.2009618, 2013. a
Pope, S.: Turbulent Flows, Cambridge University Press, 2001. a
Shapiro, C. R., Gayme, D. F., and Meneveau, C.: Modelling yawed wind turbine wakes: a lifting line approach, J. Fluid Mech., 841, R1, https://doi.org/10.1017/jfm.2018.75, 2018. a, b, c, d, e, f
Sun, J.: Vertical variations of mixing lengths under neutral and stable conditions during CASES99, J. Appl. Meteorol. Clim., 50, 2030–2041, 2011. a
Vollmer, L., Steinfeld, G., Heinemann, D., and Kühn, M.: Estimating the wake deflection downstream of a wind turbine in different atmospheric stabilities: an LES study, Wind Energ. Sci., 1, 129–141, https://doi.org/10.5194/wes11292016, 2016. a
 Abstract
 Copyright statement
 Introduction
 Aerodynamics of the curled wake: a controloriented model
 Numerical solution
 Comparison between the model and largeeddy simulations
 Controloriented modeling
 Possible improvements for the model
 Conclusions
 Code availability
 Competing interests
 Acknowledgements
 References
 Abstract
 Copyright statement
 Introduction
 Aerodynamics of the curled wake: a controloriented model
 Numerical solution
 Comparison between the model and largeeddy simulations
 Controloriented modeling
 Possible improvements for the model
 Conclusions
 Code availability
 Competing interests
 Acknowledgements
 References