Next Article in Journal
A Control Scheme without Sensors at the PV Source for Cost and Size Reduction in Two-Stage Grid Connected Inverters
Next Article in Special Issue
Simulated Annealing Algorithm for Wind Farm Layout Optimization: A Benchmark Study
Previous Article in Journal
The Effect of Renewable Energy Consumption on Sustainable Economic Development: Evidence from Emerging and Developing Economies
Previous Article in Special Issue
Aerodynamically Interacting Vertical-Axis Wind Turbines: Performance Enhancement and Three-Dimensional Flow
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Wake Modeling Paradigm for Wind Farm Design and Control

Department of Mechanical Engineering, Johns Hopkins University, Baltimore, MD 21218, USA
*
Author to whom correspondence should be addressed.
Energies 2019, 12(15), 2956; https://doi.org/10.3390/en12152956
Submission received: 28 May 2019 / Revised: 5 July 2019 / Accepted: 23 July 2019 / Published: 1 August 2019
(This article belongs to the Special Issue State of the Art of Wind Farm Optimization)

Abstract

:
Wake models play an integral role in wind farm layout optimization and operations where associated design and control decisions are only as good as the underlying wake model upon which they are based. However, the desired model fidelity must be counterbalanced by the need for simplicity and computational efficiency. As a result, efficient engineering models that accurately capture the relevant physics—such as wake expansion and wake interactions for design problems and wake advection and turbulent fluctuations for control problems—are needed to advance the field of wind farm optimization. In this paper, we discuss a computationally-efficient continuous-time one-dimensional dynamic wake model that includes several features derived from fundamental physics, making it less ad-hoc than prevailing approaches. We first apply the steady-state solution of the model to predict the wake expansion coefficients commonly used in design problems. We demonstrate that more realistic results can be attained by linking the wake expansion rate to a top-down model of the atmospheric boundary layer, using a super-Gaussian wake profile that smoothly transitions between a top-hat and Gaussian distribution as well as linearly-superposing wake interactions. We then apply the dynamic model to predict trajectories of wind farm power output during start-up and highlight the improved accuracy of non-linear advection over linear advection. Finally, we apply the dynamic model to the control-oriented application of predicting power output of an irregularly-arranged farm during continuous operation. In this application, model fidelity is improved through state and parameter estimation accounting for spanwise inflow inhomogeneities and turbulent fluctuations. The proposed approach thus provides a modeling paradigm with the flexibility to enable designers to trade-off between accuracy and computational speed for a wide range of wind farm design and control applications.

1. Introduction

Improving wind farm design and control is at the core of wind energy research aimed at facilitating greater wind energy penetration. Design problems focus on maximizing wind farm power production and reducing the levelized cost of energy by optimizing wind farm layouts [1,2,3,4] and wind turbine set points for yaw, tilt, and thrust [5,6]. Active control attempts to actuate turbines dynamically to reduce farm level power fluctuations [7], track a power output reference signal to provide power grid services [8,9,10,11], or maximize power production [12,13].
Design and control studies rely on engineering models that can be easily computed and optimized because high-fidelity simulations are too computationally expensive to be practical. Traditional engineering design models focus on wake growth rates, the shape of the velocity deficit in the transverse directions downstream of the turbine (the wake profile), and wake interactions. However, the details and capabilities of these models vary considerably. For example, the wake profile shape can be modeled as a top-hat [14,15], Gaussian [16,17], or a more complicated step-wise function [18]. Wake areas are often modeled as growing linearly with downstream distance [14,15], but have also been represented with other functional forms [19]. Velocity deficits arising from wake interactions may be computed as a linear sum [20,21] or as a sum of squares [14,15]. Related modeling approaches include the effects of the atmospheric boundary layer (ABL) on the wind turbine array through one-way [22] or two-way coupling [23,24] with ABL models [19,25,26]. Design-related models may also include certain dynamic effects such as yaw and tilt actuation [27,28,29] or the turbulent meandering of wakes [30]. These design models are typically either static or account for limited temporal variations and therefore cannot capture changing wind farm conditions or power output levels.
Active and closed-loop control applications, on the other hand, require time-dependent models that can capture the effects of control actions. Dynamic modeling approaches also vary considerably in fidelity, modeling philosophy, and underlying assumptions. The simplest models are discrete-time extensions of steady-state wake models [31] and continuous-time one-dimensional models [8,32]. Higher fidelity, yet more computationally-expensive, models include two-dimensional Reynolds-averaged Navier–Stokes (RANS) simulations [33] and very coarse large-eddy simulations (LES) [34].
Design and control-oriented models also differ in their implementation. In design problems, the goal is to estimate the effect of wind turbines on the fluid flow prior to implementation or construction of the farm. Therefore, steady-state model coefficients cannot be directly measured and are instead estimated empirically [20,21,35] based on site-specific data, or more infrequently, via physics-based methods [24]. In closed-loop control applications, parameter estimation is enabled through real-time data regarding power output or from flow field sensors that provide a means of fitting model coefficients to local site-specific attributes based on the current flow state [4,36,37].
In other ways, dynamic models have a range of complications that are not relevant to steady-state models. For example, the temporal effect of large-scale turbulent structures and the non-linearity of advection affect inter-turbine travel times. These interactions impact the instantaneous power output of individual turbines within a farm, but may not change the time-averaged power output. Therefore, these effects have not been extensively discussed in studies that focus on steady-state applications.
In this paper, we build on our previously-proposed continuous-time one-dimensional dynamic wake model [32] to provide a unified modeling framework for steady-state design and active control applications. We demonstrate how the selection of modeling parameters, including the wake mixing and expansion model, the wake profile shape, the advection scheme, and the inflow model, can improve model fidelity. For the dynamic model, we introduce a filtering approach to incorporate spanwise inhomogeneity in the inlet conditions to improve the inflow model, as well as state and parameter estimation to improve power output prediction. We then demonstrate that these estimation techniques can correct model errors, such as those introduced when using a linear advection model.
The remainder of this paper is organized as follows. We introduce the model proposed in [8] in Section 2, along with a new definition of the wake profile that improves the model accuracy. In Section 3, we discuss how to parameterize the model for design applications that employ the steady-state version of the model, where improved performance is achieved by linking the wake expansion rates to top-down atmospheric models [26]. In Section 4, we apply the dynamic model to a range of problems starting with the prediction of the power trajectory of a wind farm at start-up. We then introduce measurement-based state and parameter estimation approaches to account for spanwise inhomogeneities in the inlet and modeling errors. Finally, we summarize and discuss the conclusions in Section 5.

2. Wake Model

This paper builds on the dynamic wake model introduced in [8] and extended to include the effect of yaw in [32]. The model presented here differs from previous versions in the following ways. Velocity deficits are linearly superposed; the wake growth rate is linked to the friction velocity; and the wake profile model incorporates a smooth transition from a top-hat to a Gaussian wake profile.
We begin by considering a wind farm with N turbines each with diameter D = 2 R . The coordinate system is defined such that x is the streamwise direction, y is the spanwise direction, and z is the vertical direction. The center of the n th turbine’s rotor is located at s n = ( s x , n , s y , n , z h ) . The inflow velocity to the farm may vary in time and space U ( x , t ) ; however, in some cases, the time-averaged U ( x ) or time- and space-averaged velocity U may be the more appropriate, and the arguments of the averaged quantities will be omitted.
The streamwise velocity deficit behind the n th turbine δ u n is governed by the one-dimensional partial differential equation (PDE) (no Einstein summation convention over n is implied):
δ u n t + U a , n ( x , t ) δ u n x = w n ( x , t ) δ u n ( x , t ) + S n ( t ) G ( x s x , n ) .
The advection velocity U a , n can be specified is several ways. Here, we focus on two cases: the linear advection model of the form U a , n = U , which was used in [37] for ease of implementation, and the nonlinear form U a , n ( x , t ) = U δ u n ( x , t ) . The function G ( x ) is a normalized Gaussian function with characteristic width R:
G ( x ) = 1 R 2 π exp x 2 2 R 2 .
The function w n ( x ) defines the wake expansion rate as:
w n ( x , t ) = U a , n ( x , t ) 1 A n ( x ) d A n d x ,
which depends on the wake area A n ( x ) that is defined in Section 2.1.
The source term S n ( t ) = U a , n ( s x , n , t ) δ u 0 , n is specified by the initial velocity deficit δ u 0 , n , which is obtained from an inviscid model [14,38]. In the case where the thrust is pointing in the direction of the upstream velocity, the appropriate inviscid model that relies on actuator disk theory [38] gives an initial velocity deficit of:
δ u 0 , n = u , n ( 1 1 C T , n ) = 1 2 C T , n u d , n ,
where C T , n is the thrust coefficient and u , n is the upstream velocity at turbine n. The initial velocity deficits can also be described in terms of the local thrust coefficient C T , n and the disk-averaged velocity across the turbine rotor u d , n .
For the dynamic model, we follow the convention of LES [25,39] and use C T and u d in the model. Assuming that the resulting velocity field is consistent under Galilean transformations, we adopt a linear superposition of velocity deficits. The resulting velocity field is found by distributing the velocity deficits using the wake profile W n ( x ) defined in Section 2.2 and subtracting the result from the inlet flow velocity U ( x , t ) :
u ( x , t ) = U ( x , t ) n δ u n ( x , t ) W n ( x ) .
The linear superposition approach simplifies the implementation within a model-based control framework and agrees well with the data, as shown in Section 3. This disk-averaged velocity for turbine n, u d , n , is found by sampling the velocity field according to:
u d , n = 0 u ( x , t ) H ( R r n ) δ ( x s n , x ) 2 π r n d r n d x ,
where r n 2 = ( y s y , n ) 2 + ( z z h ) 2 is the squared radial coordinate of the n th turbine.
For the steady-state case, it is easier to use C T and u . The corresponding steady-state solution of the dynamic model [32] is:
δ u n ( x ) = δ 0 , n A n ( x ) π R 2 2 1 + erf x R 2 .
Linear superposition of the other wakes can be used to find the upstream velocity of each turbine:
u , n = 0 U ( x ) m n δ u m ( x ) W m ( x ) H ( R r n ) δ ( x s n , x ) 2 π r n d r n d x .

2.1. Turbulent Wake Mixing and Expansion

The wake model in Section 2 captures the dynamics of wake advection, inter-turbine wake interactions (through velocity deficit superposition), and the effect of wake growth on the magnitude of the velocity deficit in the wake. A functional form for the wake area A n ( x ) = π R 2 d w , n 2 ( x ) , defined using the normalized wake diameter d w , n ( x ) , that properly captures the effect of turbulent mixing of the wake with the surrounding turbulent atmospheric boundary layer is needed to specify the model fully. We now derive the required form for the wake area through the self-similarity assumption applied to the Reynolds-averaged mean momentum equation associated with a turbine in the ABL. For convenience, we dispense with the turbine indexing in this section, and the origin of the coordinate system is placed at the center of the turbine rotor.
After linearizing the advective term, assuming the pressure gradient vanishes in the region sufficiently downstream from the turbine, neglecting the streamwise Reynolds stress, and using the eddy viscosity model, the Reynolds-averaged mean momentum equation for the streamwise velocity u in the region away from the turbine forcing becomes:
U u x = 1 r r r ν T u r ,
where ν T is the eddy viscosity, and we have assumed a constant inflow velocity U for the turbine of interest. The far wake of a turbine has been shown empirically to have a self-similar profile [16]; therefore, the velocity may be written as:
u ( x , r ) = U δ u ( x ) W r ( x ) ,
where W ( ξ ) is a self-similar function of ξ = r ( x ) and ( x ) D d w ( x ) is proportional to the diameter of the wake.
For a free wake, the eddy viscosity scales with the wake width and the velocity deficit ν T δ u [40]. In the atmospheric boundary layer, however, the appropriate velocity scale is the friction velocity u * , giving the following eddy viscosity:
ν T ( x ) = u * ( x ) .
Substituting (10) and (11) into (9) and expanding the derivatives, we obtain:
U d δ u d x W + U δ u ( x ) r ( x ) 2 d d x W = u * δ u ( x ) r W u * δ u ( x ) ( x ) W ,
where primes denote derivatives with respect to ξ . Rearranging (12) in dimensionless form leads to:
U u * ( x ) δ u ( x ) d δ u d x W = U u * d d x ξ + 1 ξ W + W .
The expression (13) shows that in order for the profile to be self-similar, the following two conditions must hold in the wake:
U u * d d x = constant
U u * ( x ) δ u ( x ) d δ u d x = constant .
Integrating (14), we find that ( x ) x , and the wake diameter expands linearly with downstream distance, as in the Jensen model [14,15]. This yields a normalized wake diameter of:
d w ( x ) = 1 + 2 k w x D ,
where the wake expansion coefficient:
k w = α u * U
is proportional to the ratio of the friction velocity and hub-height velocity, validating the empirical suggestion of Frandsen [35]. The parameter α will remain as the only free parameter available to the model.
While the solution to (15) only requires δ u x n , we can find the scaling using momentum flux conservation in the far field limit using:
0 U U u ( x , r ) 2 π r d r = 0 U δ u ( x ) W ξ 2 π 2 ( x ) d ξ = U δ u ( x ) 2 ( x ) C ,
where C is a constant that depends on the initial velocity deficit in the wake δ u 0 and the integral 0 W ( ξ ) d ξ . Therefore, the velocity deficit scales inversely with the square of the downstream distance δ u x 2 . Solving the steady-state version of the PDE (1) also provides this scaling.
While the linear growth rate is valid in the wake region of the turbine, it is ill-posed upstream of the turbine, where it can vanish or become negative. Therefore, we instead use the following modified function [8] that smoothly approximates the linear expansion in the far wake while avoiding becoming less than unity close to the turbine:
d w ( x ) = 1 + k w ln 1 + e x / R .
Defined this way, the normalized wake diameter specifies the wake expansion rate (3) needed to define the velocity deficit PDE (1) fully.

2.2. Wake Profile

We now propose a new functional form for the wake profile W ( x ) in (5). Experimental observations have revealed that the shape of the wake profile as a function of spanwise distance slowly transitions from a top-hat profile (cf. the left panel of Figure 1) immediately following the rotor disk to a self-similar Gaussian profile farther downstream in the far wake [16,20,28]. For simplicity, the Jensen model traditionally used top-hat profiles for the entire downstream evolution of the wake [14,15], but such profiles poorly describe the far wake profile that dominates the interaction of upstream wakes on downstream turbines. Other models use the Gaussian profile for the entire wake [16], which overestimates the centerline deficit near the turbine, or model the transition from the top-hat as a decaying potential core [28].
Here, we propose a super-Gaussian function that more accurately approximates the known transition from a top-hat profile to a Gaussian profile. As in Section 2.1, we remove the turbine indexing for convenience and set the origin of the coordinate at the center of the turbine rotor. The proposed wake profile function is of the form:
W ( x ) = C ( x ) exp D 2 8 σ 0 2 2 r D d w ( x ) p ( x ) ,
where σ 0 = D / 4 is a parameter defining the width of the wake relative to the wake diameter D d w ( x ) and p ( x ) is a parameter that characterizes the smoothness of the wake. In order to transition from a top-hat profile at the rotor disk to a Gaussian in the far wake, we require that:
lim x 0 p =
lim x p = 2 .
Furthermore, p ( x ) should be specified such that the profile becomes nearly Gaussian at x / D 4 [28]. To satisfy these conditions, we define:
p ( x ) = 2 1 + D max ( x , 0 ) ,
where p = 2.5 at x / D = 4 .
The wake function must also be defined to conserve mass, which means that for a single turbine subject to the inflow velocity U :
0 2 π r ( U u ) d r = A ( x ) δ u ( x , t ) .
Alternatively, the integral can be expressed as:
0 2 π r W ( x ) d r = A ( x ) ,
the solution of which [41] (p. 342 #3.478.1) is:
0 2 π r W ( x , r ) d r = C ( x ) A ( x ) 2 p ( x ) 8 σ 0 2 D 2 2 / p ( x ) Γ 2 p ( x ) .
Therefore, the scaling value C ( x ) is:
C ( x ) = p ( x ) 2 Γ 2 / p ( x ) D 2 8 σ 0 2 2 / p ( x )
with lim p 2 C = D 2 / 8 σ 0 2 and lim p C = 1 .
The resulting profiles for the velocity in the wake, assuming k w = 0.1 , and C T = 2 , are shown in Figure 1. These profiles begin as a top-hat with the expected velocity of u d / U = a = 1 / 3 , at the rotor disk, where a is the induction factor. Around x / D = 1 , the velocity is a smoothed top-hat with u / U = 2 a = 2 / 3 in the middle of the wake. Further downstream, the deficit decays, transitioning to a Gaussian by mimicking the potential core region noted in [28].
This functional form captures the mean velocity in both the near wake, where the velocity deficit of an actuator disk has a top-hat profile, and the far wake, where the velocity deficit has a Gaussian profile. Small-scale features, such as trailing root and tip vortices, and turbine-specific thrust profiles are not captured by this functional form. In cases where these specifics are deemed important, the functional form can only be assumed to be valid for x / D > 2 . With the wake profile specified, the model presented in this section is fully defined with a single free parameter α that specifies the wake recovery rate inside the farm.

3. The Steady-State Model for Design Applications

For design problems, we require a method to specify the unknown coefficient α for the wake expansion coefficient k w (17). In typical models, it is obtained using empirical tables based on surface and atmospheric conditions. The recently-proposed coupled wake boundary layer (CWBL) model [23,24,42] eliminates the need for empirical tables by coupling the wake expansion coefficient of the Jensen model to the top-down model [19,25], which provides information about the turbulence characteristics of the wind turbine array boundary layer.
The CWBL determines the single free parameter of the Jensen wake model, the fully-developed wake expansion coefficient k w , , by matching the power prediction of the top-down and Jensen wake model. For a given wake expansion coefficient, the top-down wake model planform thrust coefficient is determined by finding the wake area of the farm; i.e. c f t = ( N π D 2 C T ) / ( 4 w f A p f ) , where A p f is the planform area of the farm and w f is the effective wake area coverage. The wake expansion coefficient in the developing region is then assumed to transition from k w , 0 = u * / U to k w , using a heuristic equation.
Here, we use a related approach, but make several important modifications. First, we replace the Jensen wake model with the steady-state solution of the wake model described above. The remaining differences in the approaches are highlighted in the remainder of this section. Most notably, we dispense with the rather strong assumption of the standard top-down model that the disk-averaged velocity is equal to the planar-averaged hub-height velocity u d u ¯ h to calculate the planform-averaged thrust coefficient and power production. We represent entrance effects using the developing top-down model [26,43] and replace the heuristic wake growth parameter development with the classic equation for the developing internal boundary layer.
We detail our approach by first introducing the fully-developed top-down model [25]. This model assumes constant stress regions above and below the wind turbine layer, with friction velocities of u * , h i and u * , l o and roughness heights of z 0 , h i and z 0 , l o , respectively. The planar-averaged momentum balance results in the following relationship:
u * , h i 2 = u * , l o 2 + 1 2 c f t u ¯ h 2 ,
where u ¯ h is the planar-averaged velocity at hub height z = z h . Within the wind turbine wake region z h R z z h + R , where R is the rotor radius, the wind turbine wakes are assumed to generate an additional eddy viscosity ν w . The resulting eddy viscosity in the wake region is ( 1 + ν w * ) κ u * , l o below the hub height z h R z z h and ( 1 + ν w * ) κ u * , h i above the hub height z h z z h + R , where ν w * 28 1 2 c f t .
Applying momentum conservation and velocity continuity across the layers, the roughness height of the wind farm is:
z 0 , h i = z h 1 + R z h β exp c f t 2 κ 2 + ln z h z 0 , l o 1 R z h β 2 1 / 2 ,
where β = ν w * / ( 1 + ν w * ) . The friction velocities are:
u * , h i = u * ln ( δ / z 0 , l o ) ln ( δ / z 0 , h i )
u * , l o = u * , h i ln z h z 0 , h i 1 + R z h β ln z h z 0 , l o 1 R z h β
where u * is the free-stream friction velocity and δ is the height of the boundary layer. The planar-averaged hub-height velocity is:
u ¯ h = u * , h i κ ln z h z 0 , h i 1 + R z h β .
In the CWBL model, the disk-averaged velocity is then assumed to be approximately equal to the planar-averaged hub-height velocity u d u ¯ h to calculate the power production of the farm.
In contrast to the CWBL model, we assume that the top-down model only provides useful predictions of the planar-average hub-height velocity, but not of the individual turbine power, which depends on the more local velocity at each turbine rather than the planar-averaged velocity. As a result, the proper matching condition to find the free parameter α of the present model, which scales the wake expansion coefficients as a function of friction velocity, is to require consistency between the planar averaged velocities from the top-down and the wake model. In practice, we find α that minimizes the difference between the modeled planar-averaged velocities:
min α n D u ¯ h , n w m u ¯ h , n t d 2 ,
where the superscripts t d and w m correspond to the top-down and wake models, respectively.
In order to match the planar-averaged velocity from the wake model u ¯ h , n w m corresponding to each turbine, the wind farm’s planform area must be divided into regions that “belong” to each turbine. Here, we propose to use Voronoi tessellation to identify these areas. Difficulties arise at the edge of the farm, where the edge turbines are closest to infinite sections spanning the entire planform area. We address this issue by performing an initial Voronoi tessellation and reflecting each edge turbine across the adjacent line segments of the tessellation to create ghost points. Using the ghost points, Voronoi tessellation is used again to find the planform area of each turbine. The area sections obtained through this algorithm are shown in Figure 2 and Figure 3 for the Horns Rev wind farm and a randomly-generated wind farm, respectively. The left panels show the original Voronoi tessellation with the ghost point projections. The right panels show the final tessellation of the wind farm.
We model entrance effects using the developing top-down model [26,43], where the boundary transitions from a single log-law profile with surface roughness height z 0 , l o to the fully-developed top-down model described above. In the transition region, an internal boundary layer grows from the height of the wind farm at z = z h until it reaches the final boundary layer height at z = δ . The classical internal boundary layer growth model dictates that:
δ i b l ( x ) = z h + z 0 , h i x z 0 , h i 4 / 5
until δ i b l = δ . The friction velocities (30) and (31) and planar-averaged velocity (32) evolve with streamwise distance x, where δ is replaced with the internal boundary layer height δ i b l ( x ) .
In our application of the model, the planform thrust coefficient needed for the top-down model is calculated using the force determined by the wake model:
c f t , n = C T , n u d , n 2 π R 2 A p f , n u ¯ h , n w m 2 ,
where A p f , n is the planform area of the Voronoi cell around the n th turbine. This planform thrust coefficient also provides the surface roughness z 0 , h i . Although this value should strictly represent the roughness height of the fully-developed regime, in smaller wind farms, it is difficult to find the roughness height from matching. Therefore, we assume that within each cell, the planform thrust coefficient represents the average coefficient from upstream turbines. Furthermore, the starting location of the internal boundary layer is assumed to coincide with the most upstream turbine, the wake of which, defined by the cone with diameter 2 D d w , n ( x ) , intersects the rotor of the turbine.
Using the friction velocities at each turbine, we can determine the wake expansion coefficients that the effective friction velocity at each turbine is the average of the low and high friction velocities ( u * , l o + u * , h i ) / 2 , and the advection velocity is the planar-averaged velocity u ¯ h . As a result, the wake expansion coefficient (17), which is the input to the wake model, can be written as:
k w = α 2 u * , l o + u * , h i u ¯ h .
For a single turbine, this reduces to the expected k w = α u * / U . This improvement replaces the heuristic development function with a physically-based function and eliminates the difficulty of identifying a fully-developed regime to determine the effective wake area coverage that is needed in the CWBL approach.
We validate our proposed steady state model by applying it to the Horns Rev wind farm and compared the results to LES, as well as the Jensen and CWBL models. The Horns Rev wind farm is composed of 10 rows of eight turbines with rotor diameters of D = 80 m and hub heights of z h = 70 m. The turbines are arranged such that the spacings are 7 D and 6.95 D in the streamwise and spanwise directions when the wind direction is at 270 . A friction velocity of u * = 0.45 m/s, boundary layer height of 500 m, and surface roughness of 0.002 m are selected to match the LES data. The thrust coefficient is C T = 0.78 ( C T = 1.45 ). Details about the LES setup and CWBL and Jensen model implementations can be found in [24,44].
The wind farm power production predicted by the model is compared to the LES, Jensen, and CWBL models in Figure 4. The current model shows good agreement with the LES data for all angles, with particularly good agreement for the 270 wind angle and in the early rows of the farm for all cases. The results are comparable to the CWBL model, but our approach outperforms the Jensen model. Figure 5 illustrates that the flow field obtained is more realistic than that generated by models using simple top-hat or Gaussian profiles (the latter are shown in [21]), due to the use of the wake profile function in Equation (20). In particular, the figures show that the wake profiles transitions from top-hats to Gaussians, and the expansion rate of downstream turbines is larger than the first turbines in each row, as seen in the LES data. As a result, at closer streamwise spacings, we expect the model to perform well compared to actuator disk model simulations; however, actuator line simulations with radially-varying induction factors and small-scale tip vortices will not be fully described by the model at close streamwise spacings.

4. A Control-Oriented Dynamic Model

Models for active and closed-loop control applications need to account for the effect of wind turbine control actions on other turbines. In Section 4.1, we apply the dynamic model described in Section 2 to an extreme example of such a problem, a wind farm at start-up. We simulate the model with both linear and non-linear advection U a and compare the results to ensemble-averaged LES. In Section 4.2, we provide measurement-based state and parameter estimation algorithms. The first provides a filter to obtain estimates of the time-varying, spanwise inhomogenous inflow velocities. We then apply an ensemble Kalman filter (EnKF) to correct for modeling errors, such as the linear advection assumption and additional turbulence effects, which are not fully accounted for in the model. Finally, the dynamic model with state and parameter estimation is compared to LES data.

4.1. Dynamic Thrust Modulation

In this section, we test the model’s capability to capture the dynamics of wake advection in response to dynamic changes in a wind turbine’s thrust control actions. Specifically, we examine the power output trajectory in each row of the farm using the LES results from a wind farm at start-up [8]. Both the LES and the model are applied to a wind farm composed of 14 rows of 10 turbines with streamwise and spanwise spacings of 7 D and 5 D , respectively. Each turbine has a rotor diameter of D = 100 m, hub height of z h = 100 m, and operating thrust coefficient of C T = 4 / 3 . The incoming atmospheric boundary layer flow has a friction velocity of u * = 0.517 m/s and roughness height z 0 = 0.1 m. The inflow velocity is therefore set as U ( x , t ) = ( u * / κ ) log ( z h / z 0 ) , where κ = 0.4 is the von Karman constant. The wind farm start-up is modeled as a step change in the local thrust coefficient from C T = 0 to C T = 4 / 3 . The model is simulated using the linear and nonlinear advection velocities discussed in Section 2. The wake expansion coefficients are determined using the top-down coupling approach discussed in Section 3, with C T = 4 / 3 .
The dynamic change in power for the first four rows, where the total power for the entire row is plotted, is shown in Figure 6. An ensemble of 20 LES realizations is averaged to reduce the effect of turbulent fluctuations and focus on the dynamic effects during start-up. The results demonstrate that the model is able to accurately predict the power levels of the wind farm using only the layout of the farm, the inflow friction velocity and the roughness height. The linear model results are similar to those seen in [8]. In this case, however, the wake expansion coefficients are not set by fitting the steady-state simulation data because the extended framework described in Section 3 provides a means of calculating the wake expansion coefficients. As expected, the nonlinear advection model provides a better prediction of the dynamics, particularly where the first wake hit the downstream turbine at t 1.5 min.

4.2. State and Parameter Estimation

In this section, we use the power measurements that are readily available in physical wind farms as feedback for a state and parameter estimation. This type of estimation scheme would be used in a control application in order to correct for modeling errors and unmodeled effects such as turbulent fluctuations that are clearly visible in the spectrum of wind turbine power production [45]. To demonstrate the effectiveness of this approach, we introduce a known error into the model by using the linear advection model.
We first present an empirical model of the local inlet velocity that accounts for the known large-scale streaks and rolls [46] that generate spanwise variability in the inflow velocity. The model is based on estimates of the freestream velocity from the power measurements P n ( t ) collected from unwaked turbines in the farm; i.e., those turbines whose power output is not affected by wakes of upstream turbines. The velocity magnitudes collected from these unwaked turbines are then time filtered using a linear temporal relaxation towards the instantaneous values, according to:
d u , n d t = 1 τ 4 + C T , n ( t ) 4 8 P n ( t ) ρ π D 2 C T , n ( t ) 1 / 3 u , n ( t ) n W .
The set W denotes the set of unwaked turbines in the farm, and τ = 5 min is the relaxation time constant by which the measurements are effectively filtered.
These filtered inflow velocities are used in two ways. First, the linear advection velocity model is used with the advection velocity in (1) set to the mean of the filtered inflow velocities U a = 1 N f n W u , n , where N f is the number of unwaked turbines. While this approach is known to be less accurate than using the nonlinear advection model, it introduces a known error that enables us to demonstrate the effectiveness of this approach. In particular, errors due to this assumption are mitigated through the updated measurements of the power at each turbine and at each time-step. A linear model may also be preferred to simplify computations and analysis. Second, a spanwise profile of the incoming velocity U ( y , t ) is generated by linearly interpolating between the freestream velocity points found above. This spanwise profile is used to find the velocity field U ( y , t ) for (5).
An EnKF is used to reduce the error in the velocity deficits δ u n (the states of the model) and the wake expansion coefficients k w , n (the parameters of the model) based on power measurements at each turbine. In particular, we extend the approach of Shapiro et al. [37] to incorporate the modified dynamic model presented in this paper. This extension allows the estimation method to be applied to arbitrary turbine configurations and includes spanwise variations in inflow velocities. We recall that the previous approach of Shapiro et al. [37] was limited to regular wind farm configurations in which all rows were averaged and the farm could be treated as a single column of representative turbines.
The EnKF assumes that the resulting modeled wind farm system is governed by the discrete update equations:
ψ k + 1 = f ( ψ k , C T k ) + B χ k ,
ξ k = h ( ψ k , C T k ) + ϵ k
where ψ k + 1 = f ( ψ k , C T k ) and ξ k = h ( ψ k , C T k ) are temporal and spatial discretizations of the wake model, ψ k is the vector of the model states, C T k is a vector of the local thrust coefficients at each turbine, and the function h ( ψ k , C T k ) maps the states to the measured quantities (the power at each turbine). Measurement and modeling errors are represented by the zero-mean white noise processes ϵ and χ , respectively. The assumed root-mean-square (RMS) of the measurement, state, and parameter noises must be prescribed. In the examples shown, they are specified as σ δ u = 0.1 m/s, σ k = 0.00025 , and σ P = 4.83 kW.
An ensemble of 256 wake model implementations, governed by (38) and (39), is used to estimate the error covariance of the wake model compared to measurements. The state and parameter estimates are computed in two steps. First, the wake model is used to perform an intermediate forecast of the wake model states. Then, power measurements are used in conjunction with the ensemble estimates in the measurement analysis step to produce states and parameter estimates. Further details of this method can be found in [37].
The estimation algorithm is tested using an LES of a 38-turbine wind farm composed of four staggered rows of six turbines followed by 14 irregularly-arranged turbines at the back of the farm, as shown in Figure 7. This arrangement allows us to demonstrate the utility of this approach on wind farms with both regularly- and irregularly-spaced turbines. The simulations are performed using the pseudo-spectral code LESGO [47], which is open source and publicly available at https://lesgo.me.jhu.edu/, using the concurrent precursor inflow method [48], the Lagrangian-averaged scale-dependent subgrid scale model [49], and actuator disk wind turbine models [25,39]. The domain extends 9 × 6 × 1 km in the streamwise, spanwise, and vertical directions with 384 × 256 × 128 grid points in these directions, respectively. The friction velocity is set to be u * = 0.45 m/s, and the surface roughness height is specified as z 0 = 0.1 m. The turbines have a rotor diameter of D = 100 m, a hub height of z h = 100 m, and a local thrust coefficient of C T = 4 / 3 .
The measured and estimated power production for six turbines, which are representative of the farm, are shown in Figure 8. The estimated power production show very good agreement with the measurements. The RMS error has a value of 10.4 kW, which corresponds to a percent error of 0.8 % . The modeled velocity field at t = 60 min is compared to a time-averaged LES velocity field averaged between t = 57.5 min and t = 62.5 min in Figure 9. As demonstrated by this plot, the state-estimated dynamic wake model generates realistic velocity fields at a computational cost that is significantly lower than higher fidelity approaches like LES or RANS. Specifically, for this simulation comprised of an ensemble of 256 perturbed wake models, the ensemble calculations are spread over 32 processors, which means that each processor handles eight of the ensemble members. The average time to advance these models is 0.9 s for a simulation time step of 2.1 s, which means that these calculations can be computed in real time based on typical wind farm control loops.

5. Conclusions

We introduced a continuous-time dynamic wake model and its steady-state counterpart, which provided a unifying framework for studying wind farm design and control problems. Our previously-proposed dynamic model was extended by deriving the linear wake expansion scaling from self-similarity, showing rigorously that k w u * / U and developing a realistic wake profile function that transitioned from a top-hat profile to a Gaussian. For steady-state design problems, the developing top-down model was coupled to the wake model to specify the wake expansion coefficient k w properly and closely match LES data. The dynamic model was shown to reproduce the behavior of a wind farm at start-up. State and parameter estimation was used to incorporate unmodelled turbulence effects, as well as variations in spanwise inflow velocity and streamwise advection velocities. These extensions lay the groundwork for improving wind farm design and control optimization by providing physics-based estimates of the wake expansion rate, more accurately predicting the wind farm advection time, and generating more realistic wind farm velocity fields through realistic wake profiles and state estimation.

Author Contributions

C.R.S. formulated the problem, analyzed results and wrote the first draft of the paper. G.M.S. formulated the problem, performed the simulations and analyzed results in Section 4.2. C.M. and D.F.G. helped formulate the problem, interpret the results and contributed to writing the paper.

Funding

This material is based on work supported by the National Science Foundation under Grant No. CMMI 1635430 and the National Science Foundation Graduate Research Fellowship under Grant No. DGE 1746891. Computations were made possible by Maryland Advanced Research Computing Center (MARCC).

Acknowledgments

The authors thank R.J.A.M. Stevens and J. Meyers for useful conversations at the initial stages of this research and P. Bauweraerts for the start-up LES data.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Lackner, M.A.; Elkinton, C.N. An Analytical Framework for Offshore Wind Farm Layout Optimization. Wind Eng. 2007, 31, 17–31. [Google Scholar] [CrossRef]
  2. Chowdhury, S.; Zhang, J.; Messac, A.; Castillo, L. Unrestricted wind farm layout optimization (UWFLO): Investigating key factors influencing the maximum power generation. Renew. Energy 2012, 38, 16–30. [Google Scholar] [CrossRef]
  3. Bokharaie, V.S.; Bauweraerts, P.; Meyers, J. Wind-farm layout optimisation using a hybrid Jensen–LES approach. Wind Energy Sci. 2016, 1, 311–325. [Google Scholar] [CrossRef]
  4. Adcock, C.; King, R.N. Data-driven wind farm optimization incorporating effects of turbulence intensity. In Proceedings of the 2018 Annual American Control Conference (ACC), Milwaukee, WI, USA, 27–29 June 2018; pp. 701–706. [Google Scholar] [CrossRef]
  5. Fleming, P.A.; Gebraad, P.M.; Lee, S.; van Wingerden, J.W.; Johnson, K.; Churchfield, M.; Michalakes, J.; Spalart, P.; Moriarty, P. Evaluating techniques for redirecting turbine wakes using SOWFA. Renew. Energy 2014, 70, 211–218. [Google Scholar] [CrossRef]
  6. Annoni, J.; Scholbrock, A.; Churchfield, M.; Fleming, P. Evaluating tilt for wind plants. In Proceedings of the 2017 Annual American Control Conference (ACC), Seattle, WA, USA, 24–26 May 2017; pp. 717–722. [Google Scholar] [CrossRef]
  7. De Rijcke, S.; Driesen, J.; Meyers, J. Power smoothing in large wind farms using optimal control of rotating kinetic energy reserves. Wind Energy 2015, 18, 1777–1791. [Google Scholar] [CrossRef]
  8. Shapiro, C.R.; Bauweraerts, P.; Meyers, J.; Meneveau, C.; Gayme, D.F. Model-based receding horizon control of wind farms for secondary frequency regulation. Wind Energy 2017, 20, 1261–1275. [Google Scholar] [CrossRef]
  9. Bay, C.J.; Annoni, J.; Taylor, T.; Pao, L.; Johnson, K. Active power control for wind farms using distributed model predictive control and nearest neighbor communication. In Proceedings of the 2018 Annual American Control Conference (ACC), Milwaukee, WI, USA, 27–29 June 2018; pp. 701–706. [Google Scholar] [CrossRef]
  10. Boersma, S.; Doekemeijer, B.; Siniscalchi-Minna, S.; van Wingerden, J. A constrained wind farm controller providing secondary frequency regulation: An LES study. Renew. Energy 2019, 134, 639–652. [Google Scholar] [CrossRef]
  11. Vali, M.; Petrović, V.; Boersma, S.; van Wingerden, J.W.; Pao, L.Y.; Kühn, M. Adjoint-based model predictive control for optimal energy extraction in waked wind farms. Control Eng. Pract. 2019, 84, 48–62. [Google Scholar] [CrossRef]
  12. Goit, J.P.; Meyers, J. Optimal control of energy extraction in wind-farm boundary layers. J. Fluid Mech. 2015, 768, 5–50. [Google Scholar] [CrossRef] [Green Version]
  13. Munters, W.; Meyers, J. Dynamic strategies for yaw and induction control of wind farms based on large-eddy simulation and optimization. Energies 2018, 11, 177. [Google Scholar] [CrossRef]
  14. Jensen, N.O. A Note on Wind Generator Interaction; Technical Report Risø-M-2411; Risø National Laboratory: Roskilde, Denmark, 1983. [Google Scholar]
  15. Katic, I.; Højstrup, J.; Jensen, N.O. A simple model for cluster efficiency. In Proceedings of the European Wind Energy Association Conference and Exhibition, Rome, Italy, 7–9 October 1986; pp. 407–410. [Google Scholar]
  16. Bastankhah, M.; Porté-Agel, F. A new analytical model for wind-turbine wakes. Renew. Energy 2014, 70, 116–123. [Google Scholar] [CrossRef]
  17. Ge, M.; Wu, Y.; Liu, Y.; Yang, X.I. A two-dimensional Jensen model with a Gaussian-shaped velocity deficit. Renew. Energy 2019, 141, 46–56. [Google Scholar] [CrossRef]
  18. Gebraad, P.M.O.; van Wingerden, J.W. A Control-Oriented Dynamic Model for Wakes in Wind Plants. J. Phys. Conf. Ser. 2014, 524, 012186. [Google Scholar] [CrossRef]
  19. Frandsen, S.; Barthelmie, R.; Pryor, S.; Rathmann, O.; Larsen, S.; Højstrup, J.; Thøgersen, M. Analytical modelling of wind speed deficit in large offshore wind farms. Wind Energy 2006, 9, 39–53. [Google Scholar] [CrossRef]
  20. Niayifar, A.; Porté-Agel, F. A new analytical model for wind farm power prediction. J. Phys. Conf. Ser. 2015, 625, 012039. [Google Scholar] [CrossRef]
  21. Niayifar, A.; Porté-Agel, F. Analytical Modeling of Wind Farms: A New Approach for Power Prediction. Energies 2016, 9, 741. [Google Scholar] [CrossRef]
  22. Peña, A.; Rathmann, O. Atmospheric stability-dependent infinite wind-farm models and the wake-decay coefficient. Wind Energy 2013, 17, 1269–1285. [Google Scholar] [CrossRef]
  23. Stevens, R.J.A.M.; Gayme, D.F.; Meneveau, C. Coupled wake boundary layer model of wind-farms. J. Renew. Sustain. Energy 2015, 7, 023115. [Google Scholar] [CrossRef] [Green Version]
  24. Stevens, R.J.A.M.; Gayme, D.F.; Meneveau, C. Generalized coupled wake boundary layer model: Applications and comparisons with field and LES data for two wind farms. Wind Energy 2016, 19, 2023–2040. [Google Scholar] [CrossRef]
  25. Calaf, M.; Meneveau, C.; Meyers, J. Large eddy simulation study of fully developed wind-turbine array boundary layers. Phys. Fluids 2010, 22, 015110. [Google Scholar] [CrossRef] [Green Version]
  26. Meneveau, C. The top-down model of wind farm boundary layers and its applications. J. Turbul. 2012, 13, N7. [Google Scholar] [CrossRef]
  27. Martínez-Tossas, L.A.; Annoni, J.; Fleming, P.A.; Churchfield, M.J. The aerodynamics of the curled wake: A simplified model in view of flow control. Wind Energy Sci. 2019, 4, 127–138. [Google Scholar] [CrossRef]
  28. Bastankhah, M.; Porté-Agel, F. Experimental and theoretical study of wind turbine wakes in yawed conditions. J. Fluid Mech. 2016, 806, 506–541. [Google Scholar] [CrossRef]
  29. Jiménez, Á.; Crespo, A.; Migoya, E. Application of a LES technique to characterize the wake deflection of a wind turbine in yaw. Wind Energy 2010, 13, 559–572. [Google Scholar] [CrossRef]
  30. Larsen, G.C.; Madsen, H.A.; Thomsen, K.; Larsen, T.J. Wake meandering: a pragmatic approach. Wind Energy 2008, 11, 377–395. [Google Scholar] [CrossRef]
  31. Gebraad, P.M.O.; Teeuwisse, F.W.; van Wingerden, J.W.; Fleming, P.A.; Ruben, S.D.; Marden, J.R.; Pao, L.Y. Wind plant power optimization through yaw control using a parametric model for wake effects-a CFD simulation study. Wind Energy 2016, 19, 95–114. [Google Scholar] [CrossRef]
  32. Shapiro, C.R.; Gayme, D.F.; Meneveau, C. Modelling yawed wind turbine wakes: a lifting line approach. J. Fluid Mech. 2018, 841. [Google Scholar] [CrossRef] [Green Version]
  33. Boersma, S.; Doekemeijer, B.; Vali, M.; Meyers, J.; van Wingerden, J.W. A control-oriented dynamic wind farm model: WFSim. Wind Energy Sci. 2018, 3, 75–95. [Google Scholar] [CrossRef] [Green Version]
  34. Bauweraerts, P.; Meyers, J. On the Feasibility of Using Large-Eddy Simulations for Real-Time Turbulent-Flow Forecasting in the Atmospheric Boundary Layer. Bound.-Layer Meteorol. 2019, 171, 213–235. [Google Scholar] [CrossRef]
  35. Frandsen, S. On the wind speed reduction in the center of large clusters of wind turbines. J. Wind Eng. Ind. Aerodyn. 1992, 39, 251–265. [Google Scholar] [CrossRef]
  36. Doekemeijer, B.M.; Boersma, S.; Pao, L.Y.; van Wingerden, J.W. Ensemble Kalman filtering for wind field estimation in wind farms. In Proceedings of the 2017 Annual American Control Conference (ACC), Seattle, WA, USA, 24–26 May 2017; pp. 19–24. [Google Scholar] [CrossRef]
  37. Shapiro, C.R.; Meyers, J.; Meneveau, C.; Gayme, D.F. Dynamic wake modeling and state estimation for improved model-based receding horizon control of wind farms. In Proceedings of the 2017 Annual American Control Conference (ACC), Seattle, WA, USA, 24–26 May 2017; pp. 709–716. [Google Scholar] [CrossRef]
  38. Burton, T.; Jenkins, N.; Sharpe, D.; Bossanyi, E. Wind Energy Handbook; John Wiley & Sons: Chichester, West Sussex, UK, 2011. [Google Scholar]
  39. Meyers, J.; Meneveau, C. Large Eddy Simulations of Large Wind-Turbine Arrays in the Atmospheric Boundary Layer. In Proceedings of the 48th AIAA Aerospace Sciences Meeting, Orlando, FL, USA, 4–7 January 2010; p. 827. [Google Scholar] [CrossRef]
  40. Tennekes, H.; Lumley, J.L. A First Course in Turbulence; MIT: Cambridge, MA, USA, 1972. [Google Scholar]
  41. Gradshteyn, I.; Ryzhik, I. Table of Integrals, Series, and Products: Corrected and Enlarged Edition; Academic Press: New York, NY, USA, 1980. [Google Scholar]
  42. Stevens, R.J.A.M.; Gayme, D.; Meneveau, C. Using the coupled wake boundary layer model to evaluate the effect of turbulence intensity on wind farm performance. J. Phys. Conf. Ser. 2015, 625, 012004. [Google Scholar] [CrossRef] [Green Version]
  43. Crespo, A.; Frandsen, S.; Gómez-Elvira, R.; Larsen, S.E.; Petersen, E.L.; Hjuler Jensen, P.; Rave, K.; Helm, P.; Ehmann, H. Modelization of a large wind farm, considering the modification of the atmospheric boundary layer. In Proceedings of the European Wind Energy Conference, Nice, France, 1–5 March 1999. [Google Scholar]
  44. Porté-Agel, F.; Wu, Y.T.; Chen, C.H. A Numerical Study of the Effects of Wind Direction on Turbine Wakes and Power Losses in a Large Wind Farm. Energies 2013, 6, 5297–5313. [Google Scholar] [CrossRef]
  45. Stevens, R.J.A.M.; Meneveau, C. Temporal structure of aggregate power fluctuations in large-eddy simulations of extended wind-farms. J. Renew. Sustain. Energy 2014, 6, 043102. [Google Scholar] [CrossRef] [Green Version]
  46. Kline, S.J.; Reynolds, W.C.; Schraub, F.A.; Runstadler, P.W. The structure of turbulent boundary layers. J. Fluid Mech. 1967, 30, 741. [Google Scholar] [CrossRef]
  47. Stevens, R.J.A.M.; Martínez, L.A.; Meneveau, C. Comparison of wind farm large eddy simulations using actuator disk and actuator line models with wind tunnel experiments. Renew. Energy 2018, 116, 470–478. [Google Scholar] [CrossRef]
  48. Stevens, R.J.; Graham, J.; Meneveau, C. A concurrent precursor inflow method for Large Eddy Simulations and applications to finite length wind farms. Renew. Energy 2014, 68, 46–50. [Google Scholar] [CrossRef] [Green Version]
  49. Bou-Zeid, E.; Meneveau, C.; Parlange, M. A scale-dependent Lagrangian dynamic model for large eddy simulation of complex turbulent flows. Phys. Fluids 2005, 17, 025105. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Profiles of streamwise velocity u ( x , r ) at various downstream distances based on the wake profile function W ( x ) in (20).
Figure 1. Profiles of streamwise velocity u ( x , r ) at various downstream distances based on the wake profile function W ( x ) in (20).
Energies 12 02956 g001
Figure 2. Voronoi tessellation of the Horns Rev wind farm showing wind turbine locations (black dots), ghost point locations (blue dots), and Voronoi cells (black lines). The left panel shows the original tessellation with ghost point projections. The right panel shows the final tessellation with the ghost points removed.
Figure 2. Voronoi tessellation of the Horns Rev wind farm showing wind turbine locations (black dots), ghost point locations (blue dots), and Voronoi cells (black lines). The left panel shows the original tessellation with ghost point projections. The right panel shows the final tessellation with the ghost points removed.
Energies 12 02956 g002
Figure 3. Voronoi tessellation of a randomly-generated wind farm showing wind turbine locations (black dots), ghost point locations (blue dots), and Voronoi cells (black lines). The left panel shows the original tessellation with ghost point projections. The right panel shows the final tessellation with the ghost points removed.
Figure 3. Voronoi tessellation of a randomly-generated wind farm showing wind turbine locations (black dots), ghost point locations (blue dots), and Voronoi cells (black lines). The left panel shows the original tessellation with ghost point projections. The right panel shows the final tessellation with the ghost points removed.
Energies 12 02956 g003
Figure 4. Wind farm power production of Columns 2–4 (counting from the south) at various wind angles from large-eddy simulations (LES) (black), the Jensen model (red), the coupled wake boundary layer (CWBL) (blue), and the present model (green). LES data were digitally extracted from [24].
Figure 4. Wind farm power production of Columns 2–4 (counting from the south) at various wind angles from large-eddy simulations (LES) (black), the Jensen model (red), the coupled wake boundary layer (CWBL) (blue), and the present model (green). LES data were digitally extracted from [24].
Energies 12 02956 g004
Figure 5. Velocity field of Horns Rev wind farm at various wind angles. Results compare the current model (left) to time-averaged LES velocity fields (right). LES fields reproduced from [44], which is licensed under the Creative Commons Attribution CC-BY 3.0 license.
Figure 5. Velocity field of Horns Rev wind farm at various wind angles. Results compare the current model (left) to time-averaged LES velocity fields (right). LES fields reproduced from [44], which is licensed under the Creative Commons Attribution CC-BY 3.0 license.
Energies 12 02956 g005
Figure 6. Power production of a wind farm at start-up comparing LES (symbols), the model with linear advection (red), and the model with nonlinear advection (blue). Only the first four rows of the farm are shown.
Figure 6. Power production of a wind farm at start-up comparing LES (symbols), the model with linear advection (red), and the model with nonlinear advection (blue). Only the first four rows of the farm are shown.
Energies 12 02956 g006
Figure 7. Locations of the 38 turbines in the wind farm simulated in LES.
Figure 7. Locations of the 38 turbines in the wind farm simulated in LES.
Energies 12 02956 g007
Figure 8. Measured (——)and estimated power production (– – –) for the 38-turbine wind farm. The turbine locations are shown in Figure 7.
Figure 8. Measured (——)and estimated power production (– – –) for the 38-turbine wind farm. The turbine locations are shown in Figure 7.
Energies 12 02956 g008
Figure 9. Comparison of state-estimated dynamic wake model velocity fields at t = 60 min and LES averaged between t 57.5 min and t = 62.5 min.
Figure 9. Comparison of state-estimated dynamic wake model velocity fields at t = 60 min and LES averaged between t 57.5 min and t = 62.5 min.
Energies 12 02956 g009

Share and Cite

MDPI and ACS Style

Shapiro, C.R.; Starke, G.M.; Meneveau, C.; Gayme, D.F. A Wake Modeling Paradigm for Wind Farm Design and Control. Energies 2019, 12, 2956. https://doi.org/10.3390/en12152956

AMA Style

Shapiro CR, Starke GM, Meneveau C, Gayme DF. A Wake Modeling Paradigm for Wind Farm Design and Control. Energies. 2019; 12(15):2956. https://doi.org/10.3390/en12152956

Chicago/Turabian Style

Shapiro, Carl R., Genevieve M. Starke, Charles Meneveau, and Dennice F. Gayme. 2019. "A Wake Modeling Paradigm for Wind Farm Design and Control" Energies 12, no. 15: 2956. https://doi.org/10.3390/en12152956

APA Style

Shapiro, C. R., Starke, G. M., Meneveau, C., & Gayme, D. F. (2019). A Wake Modeling Paradigm for Wind Farm Design and Control. Energies, 12(15), 2956. https://doi.org/10.3390/en12152956

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop