A Novel Analytical Wake Model with a Cosine-Shaped Velocity Deficit

Abstract: A novel analytical model is proposed and validated in this paper to predict the velocity deficit in the wake downwind of a wind turbine. The model is derived by employing mass and momentum conservation and assuming a cosine-shaped distribution for the velocity deficit. In this model, a modified wake growth rate rather than a constant one is chosen to take into account the effects of the ambient turbulence and the mechanical turbulence generated. The model was tested against field observations, wind-tunnel measurements in different thrust operations and high-resolution large-eddy simulations (LES) for two aerodynamic roughness lengths. It was found that the normalized velocity deficit predicted by the proposed model shows good agreement with experimental and numerical data in terms of shape and magnitude in the far wake region (x/d0 > 3). Based on the proposed model, predictions from multiple views and at different locations are demonstrated to show the spatial distribution of streamwise velocity downwind of a wind turbine. The result shows that the model is suitable for predicting streamwise velocity fields and thus could provide some references for the selection of wind turbine spacing.


Introduction
Wind energy is one of the most profitable renewable resources and is expected to develop substantially all over the world [1]. Wind turbines in wind farms have to operate in the downwind wake flow which is characterized by lower mean velocity and higher turbulence intensity than those under unaffected conditions. As a result, wind-turbine wakes lead to overall power losses in large wind farms and increased fatigue loading on turbines [2][3][4]. It is reported that the average power losses due to the wake effect could be of the order 10-20% of overall power output in wind farms [3]. Due to the need for large wind farm installation, wind-turbine wake has become a significant topic in wind energy [5] and many researchers have conducted extensive research on turbine wakes via approaches such as field observations, wind-tunnel tests and numerical simulations. Some non-exhaustive works related to wind-turbine wake are outlined here. For instance, Virtual Wakes Laboratory conducted a series of measurements on offshore wind farms, such as Horns Rev, Middelgrunden, Nysted and Vindeby [6]. Zhang et al. [7], Howland et al. [8] and Bastankhah and Porté-Agel [9] carried out systematic wind-tunnel experiments to study the development of wake structures. Abkar et al. [10] used the large-eddy simulations (LES) technique to investigate the wake flow in a wind farm during a diurnal cycle and Qian and Ishihara [11] employed the actuator disk method to study the wake deflection caused by a yawed wind turbine. Although experimental and numerical approaches are powerful tools in the research of wind-turbine wakes, physically-based analytical models are still extensively used to estimate wake deficit and power production because of their simplicity, good accuracy and low computational cost [12,13]. Various analytical models were proposed (e.g., [4,[14][15][16][17][18][19][20]) and applied in practical engineering applications. For example, Sun and Yang [21] adopted three different wake models to estimate the wind energy output in Hong Kong. Shao et al. [22] combined the well-known Jensen model with a new interaction model to evaluate the power losses in the Lillgrund and Horn Rev offshore wind farms with different inflow directions. Gao et al. [19] employed a two-dimensional Jensen-Gaussian wake model to optimize wind turbine layout within a wind farm based on the multiple population genetic algorithm. Many studies have informed us of the importance of wake model selection and refinement in terms of accuracy. In summary, the analytical wake model is an efficient and useful tool with which to evaluate wake-turbine flows and wake effects on power generation. Thus, it has played a significant role in the wind energy community.

Review of Some Previous Wake Models
The Jensen model is one of the pioneering analytical models which assumes a top-hat distribution for the velocity deficit and employs mass conservation [14]: where U ∞ is the incoming velocity at hub height, U W is the velocity in the streamwise direction, C T is the thrust coefficient of the turbine, k t is the constant rate of wake expansion, d 0 is the diameter of the rotor and x is the downwind distance from the wind turbine. In this model, the wake width r W is characterized by a linear law (see Figure 1): r W is the wake width which indicates the distance between the rotor axis and the first point where the velocity resumes the freestream value. Jensen [14] suggested a constant value k t = 0.01 and Barthelmie et al. [23] considered that the values are 0.075 and 0.05 in onshore and offshore cases respectively. Alternatively, k t can be estimated by the following relation [24]: where z h and z 0 are the turbine hub height and the roughness length respectively. The famous Jensen model has gained its popularity in many pieces of commercial software (e.g., WAsP, WindFarmer, WindPRO and OpenWind) [6,23,25,26]. Later, Frandsen et al. [27] proposed a wake model by employing mass and momentum conservation and the same top-hat assumption: where A 0 and A W are the cross-sectional areas of the rotor and the wake respectively. In this model, the wake radius r W grows nonlinearly: where α is of order 10k t and β is given as: Without considering radial variation of the wake, the 1D top-hat models may give rise to considerable errors in the power estimation of a wind farm with varying wind directions [17]. In addition, underestimation of the maximum velocity deficit is a frequently mentioned deficiency [28]. Hence, developing a model with a more authentic velocity deficit profile is still appealing. Many previous studies (e.g., [29][30][31][32][33]) reported the velocity deficit in the wake exhibits a Gaussian or trigonometric behavior in the far wake region. Based on those findings, Bastankhah and Porté-Agel [17] proposed a 2D analytical model with a Gaussian shape distribution for the velocity deficit by using the mass and momentum conservation principle: where σ is standard deviation of the velocity deficit, which is given by: where r is the radial distance from the center of the wake, k * is the wake growth rate and is the initial wake width. It should be noted that k * is different from k t . Niayifar and Porte-Agel [34] proposed an empirical linear equation to estimate k * based on numerical simulation data [30]: where I 0 is the streamwise turbulence intensity in the incoming flow. Compared to the 1D model, this Gaussian model greatly improves the prediction of wake flows and has shed light on new theoretical research. For example, Xie and Archer [35] made some modifications to the Bastankhah and Porté-Agel model by introducing an elliptic Gaussian distribution to take into consideration the anisotropic wake expansion. The elliptic Gaussian function was also applied by Abkar and Porté-Agel [36] to characterize the profile of wake deficit for different atmospheric stability. Ishihara and Qian [4] proposed a modified Gaussian model based on the Bastankhah and Porté-Agel model, which improves the prediction of velocity deficit in the near wake region. Other Gaussian-like models can be found in these papers [16,19,20]. Instead of using a Gaussian distribution, Tian et al. [18] chose a cosine function to redistribute the mass flux predicted by the Jensen model: In this model, the wake decay coefficient k w is associated with turbulence intensities: This model could give a good prediction in terms of both shape and velocity magnitude [18]. It also successfully tackled multi-wake problems and provided a reasonable estimation of power production in a large offshore wind farm [33]. From the aforementioned introduction, it can be seen that Gaussian or trigonometric functions are better approximations for wake deficit modeling compared with single top-hat shape. However, there are some challenges when using the Gaussian-like models, which prevent them from being widely used in real-world wind farm projects. For example, in the Bastankhah and Porté-Agel model (see Equation (7)), the key empirical parameter k * = ∂σ/∂x has to be determined by data fitting from experimental or numerical data on a case-by-case basis [37], which could be quite difficult for engineers to put into practice because they know little prior information about the standard deviation of wind-turbine wake. Standard deviation σ is often used as an abstract characterization of wake width for the Gaussian-like models (e.g., σ [17], 2σ [38], 2.58σ [19] and 2 √ ln 2σ [4]) and many complicated definitions make it hard to determine the appropriate range of expansion rates for guidelines. In fact, the Gaussian function only approaches to zero when infinitely far away from the center, while the authentic profile of the velocity deficit tends toward zero much faster than the Gaussian shape [38]. Rather than using the above abstract concepts, the wake width in the trigonometric model proposed by Tian et al. [18] acts as an interface where velocity recovers to freestream value, which is physically more intuitive. Although the accuracy of interface prediction depends on the chosen combination of parameters (e.g., I 0 and C T ), the model conforms to the physical law of wake development and good results can still be obtained [18,33]. It should be mentioned that the model is derived by an equivalence to the mass flux predicted by the Jensen model. Thus, the balance of momentum is not taken into account. In general, analytical wake models are categorized into two types: those that rely on mass conservation, and others that are based on both mass and momentum theory. We consider that employing the momentum principle to derive wake models is a more advanced approach due to the fact that the turbine-induced thrust is directly linked to the momentum deficit in the flow governing equation of momentum, which better captures the physical essence behind wind-turbine aerodynamics. In order to take advantage of the superiority of momentum principle and the intuitive characterization for wake expansion, there is a need to propose a new trigonometric model employing both mass and momentum conservation, which is a research gap that deserves to be further investigated and studied.

Current Study
In this paper, a new analytical wake model is proposed and validated. The comparison between the proposed model and the previous wake models is explicitly listed in Table 1. In Section 2, a cosine-shaped distribution is considered for the velocity deficit in the wake, and mass and momentum conservation are employed to derive the model. An nonlinear rate is considered for the development of the wake expansion in this model. In Section 3, the proposed model is then validated against the field measurements [39], the wind-tunnel tests in two different thrust operations [40] and the numerical simulations under two different atmospheric turbulence conditions [30]. Based on the proposed model, a series of predictions at different locations are demonstrated in Section 4 to show the spatial distribution of the normalized velocity downwind of a wind turbine. Finally, conclusions and future research are presented in Section 5.

Derivation of the Wake Model
We apply the mass equation to the control volume (see Figure 2): Figure 2. Schematic of the control volume around the wind turbine.
Similarly, the momentum theory is employed to the same volume, which is extended upwind and downwind far enough so that the pressure at the surface becomes equal to the freestream pressure. If it is appropriate to neglect pressure and viscous terms, the following relation can be obtained [41]: where T is the thrust force exerted by the wind turbine on the flow, which indicates the momentum transfer between the turbine and flow. It can be characterized by the thrust coefficient [42]: Multiplying both sides of Equation (12) by U ∞ and then substituting it into Equation (13) gives: In order to derive the wake model, the following assumption is made: Assumption: A trigonometric distribution is considered for the velocity deficit and the normalized velocity deficit can be expressed as: where A, K and B need to be determined. To appropriately characterize the profile of the velocity deficit, the period of cosine function should be twice the wake radius (see Figure 3): Figure 3. Assumption of the trigonometric distribution for velocity deficit.
When r tends to the wake radius r W , the wake velocity recovers to the freestream velocity: By solving the Equations (17) and (18), we can obtain: Inserting Equation (19) into Equation (16) yields: Substituting Equations (14) and (20) into Equation (15) and integrating over the cross-section area of wake region A W (x) gives: There are two solutions to Equation (21), but only one of them is physically meaningful, which estimates a smaller value of velocity deficit at larger downwind distance. Hence, the appropriate solution is expressed as: Inserting Equation (22) into Equation (20) and rearranging yields: From Equation (23) it can be seen that the model preserves the self-similarity of velocity deficit and the maximum normalized velocity deficit occurs at the center of the wake: The only parameter that has not yet been determined is the wake radius r W and it plays an important role in predicting the velocity deficit in the wake. Hence, we need to estimate it as accurately as possible. A common practice is to introduce a modified wake growth rate for the wake expansion in consideration of the effects of both ambient turbulence and rotor generated turbulence [18][19][20]33,43]: where k W is the modified wake growth rate and I W is the streamwise turbulence intensity. It should be noted that Equation (25) is an empirical formula. There are several turbulence intensity models (e.g., [4,[44][45][46]) proposed to evaluate the added turbulence intensity I + which is defined as a function of I 0 and I W : Based on both experimental and numerical methods, Crespo and Hernández [46] suggested an empirical relation for the parameter ranges 5 < x/d < 15, 0.07 < I 0 < 0.14 and 0.1 < a < 0.4: where a is the induced velocity factor, equal to (1 − √ 1 − C T )/2. Recently, a more sophisticated turbulence intensity model with a dual-Gaussian shape was proposed by Ishihara and Qian [4]. Turbulence intensity at the top tip height is estimated by these two models in Figure 4. It can be seen that the Ishihara and Qian model predicts slightly larger I W for small thrust cases of C T = 0.3, 0.5, while the discrepancies between the models decrease with downwind distance. The Crespo and Hernández model only estimates larger I W in the very far wake region for a high thrust case of C T = 0.7. The discrepancy between the turbulence intensity estimated by the models exists, but it is not very evident. In our research, the Crespo and Hernández model is employed to couple with the proposed wake model (Equation (23)).  Figure 5 and the procedure of using this model is divided into the following three steps: 1. Estimate the streamwise turbulence intensity I W at a specific downwind distance according to Equations (26) and (27)

Validation of the Proposed Model
In order to validate the proposed wake model, published measurement data and LES results are cited in our study.

Field Measurement
In the 1980s, field measurements of Nibe B wind turbine were conducted by Taylor [39]. The layout of the Nibe site is shown in Figure 6. Four meteorological masts (M1, M2, M3, M4) provided wake measurements data at different downstream locations behind the Nibe B wind turbine (x = 2.5d 0 , 4d 0 , 6d 0 , 7.5d 0 ). Nibe A wind turbine located 5d 0 downstream of Nibe B turbine was not in operation during the measurements. Due to the obstacle effect of this standstill turbine Nibe A, the field measurement data from M3 (x = 6.0d 0 ) were excluded in this discussion. All the measured data downwind of Nibe B turbine were recorded using an average period of 1 min. The rotor diameter and the hub height of the turbine were 40 m and 45 m respectively. The mean wind velocity at hub height was 8.55 m/s with streamwise turbulence intensity I 0 = 10% and the roughness length 0.07 m. Based on this situation, the thrust coefficient was estimated to be C T = 0.82. For comparison, various models and the LES results published in Troldborg et al. [47] were tested against the field measurement data at two downwind distances (x = 4d 0 , 7.5d 0 ); see Figure 7. The LES data and the models except for the top-hat ones overpredict the maximum velocity deficit at downwind distance 4d 0 . It may be the reason that in reality Nibe A turbine located 1d 0 downstream of M2 would have a obstacle blocking effect on upstream wind fields although it was not in operation during the measurement. In contrast, the Bastankhah and Porté-Agel model, Ishihara and Qian model and the proposed model as well as LES data show a good agreement in predicting the maximum velocity deficit at 7.5d 0 . It can be seen that both Gaussian and cosine shapes are much more suitable for describing the profiles of the velocity deficit compared with top-hat shape. From Equations (10) and (23) we find that Tian model and the proposed model are directly related to the variable wake decay rate k W which takes into account the effects of the ambient turbulence and the added turbulence, while the proposed one offers some improvement over the Tian model by considering momentum theory.

Wind-Tunnel Measurement
In 1989, Garrad Hassan conducted a wind-tunnel measurement toward turbine wakes at the Marchwood Engineering Laboratory's atmospheric boundary layer wind tunnel. The corresponding wind-tunnel measurement data were published by Schlez et al. [40]. The wind-tunnel test chose a 1/160 scale horizontal axis wind turbine, which corresponds to a full-scale turbine with a 43.2 m rotor and a 50 m hub. The model site was flat with a uniform roughness length of 0.075 m and the free wind speed at hub height was 5.3 m/s. Two different operating conditions C T = 0.62 and C T = 0.85 were chosen in the study.
In Figure 8, the predictions obtained from different wake models are compared with measurement data for a wind turbine operating at two different thrust coefficients. In the case of the lower thrust coefficient, the proposed model exhibits a better performance than other models in terms of amplitude and shape of the velocity deficit at selected downwind locations. In the case of higher coefficient, the mean velocity at the center of the wake predicted by the proposed model is slightly underestimated by 6% at downwind distance 5d 0 , while the discrepancy becomes less at 7d 0 . It was also found that a narrower wake is estimated by this model in both cases. However, the top-hat series models underestimate the velocity deficit at the center of the wake and overestimate it in proximity to the boundary of the wake in both cases. The result obtained from the Bastankhah and Porté-Agel model does not compare well with experimental data either. This can be explained by the fact that the wake expansion parameter k * (Equation (9)) is obtained by linear fitting using just a few cases and it is not calibrated correctly. More research should be carried out so that generalization of k * can incorporate a wider range of operating conditions and incoming flow characteristics. In general, the wind turbine operating at a higher thrust coefficient is responsible for a larger velocity deficit. By comparing the results of two different operations, it can be found that the proposed model is in accordance with this property. The recovery rates of velocity at the center of the wake predicted by the proposed model are about 20% and 25% from distance 5d 0 to 7d 0 corresponding to cases of C T = 0.62 and C T = 0.85 respectively. This is attributed to the fact that mechanical generated turbulence strengthens the turbulent mixing process in the far wake region, resulting in a faster wake recovery when the wind turbine operates in larger thrust condition.

Numerical Simulation
Results are also compared with Wu and Porté-Agel [30], who performed a numerical simulation to investigate the atmospheric turbulence effects on wind-turbine wakes. A Lagrangian scale-dependent dynamic model in combination with actuator-disk model with rotation was applied in their work. Two cases of different roughness lengths, including grasslands (z 0 ≈ 0.005 m) and snow-covered flats (z 0 ≈ 0.00005 m) corresponding to a streamwise turbulence intensity 6.9% and 4.8% respectively, were chosen to verify the accuracy of the proposed wake model. The Vestas V80-2MW wind turbine with a 80 m rotor and a 70 m hub was adopted in the simulation. The incoming wind speed at hub height was 9 m/s. Under this situation, the wind turbine had a thrust coefficient C T ≈ 0.8.
Horizontal and vertical profiles of the velocity deficit in two roughness lengths for different analytical models and the LES data are presented in Figure 9. The figure shows that the results obtained from the proposed model are in acceptable agreement with the LES data not only in the hub-height plane but also at different heights, although a minor overestimation of the velocity deficit was found near the ground. The Ishihara and Qian model moderately underestimates the velocity deficit in both cases. The Bastankhah and Porté-Agel model performs well in the case of z 0 = 0.005 m but it markedly overestimates deficit in case of z 0 = 0.00005 m. The discrepancies are attributable to the incorrect calibration of k * . The maximum normalized velocity deficit predicted by various models and the LES data are plotted as a function of the normalized downwind distance (x/d 0 > 3) in Figure 10. To make a quantitative analysis, a relative error δ is defined as: where ∆U i and ∆U i represent the maximum velocity deficit predicted by the models and LES respectively. The relative errors of different models are presented in Table 2, and it was found that the proposed model exhibits a more acceptable agreement with the LES data compared with other models in both cases. In Figure 10, the proposed model slightly underestimates the minimum velocity, with the largest discrepancies being up to about 8% and 13% at location 4d 0 in the cases of z 0 = 0.005 m and z 0 = 0.00005 m respectively. It was also found that the Jensen model predicts the maximum normalized velocity deficit with accuracy in the very far wake region (x/d 0 > 8), but its underestimation is quite obvious with a smaller downwind distance (3 < x/d 0 < 8), mainly due to the unrealistic assumption of the top-hat shape . Table 2. Relative errors of the models in two cases with different roughness lengths. The figure also shows that the wake recovers faster in the case of stronger incoming turbulence. This can be explained by the fact that higher ambient turbulence enhances the mixing process in the flows and finally leads to a faster wake recovery. It can be found that this feature is shown in the proposed model.

Predictions of the Proposed Model
In order to illustrate the predictions of the proposed model, contour plots from three different views will be discussed below. As shown in Figure 11, when a wind turbine operates in the turbulent boundary layer, the wakes behind the rotor do not show axisymmetric characteristics, whereas the velocity deficit is nearly axisymmetric after some downwind distance [29]. Therefore, the predicted velocity at a point (x, y, z) is determined by means of superposition of upstream velocity and velocity deficit: where U upstream is the logarithmic incoming wind speed, which is expressed as: There are three steps to estimate the velocity at a point (x, y, z). The first step is to determine the parameters such as the operating conditions and the inflow characteristics. The second is to calculate the wake width according to Equation (25) and the distance between the point and the rotor axis. Finally, the wind speed can be obtained by Equation (23). It should be noted that the point is within the wake-influenced zone if r < r W . High-fidelity numerical simulations with four different incoming turbulence intensities (cases 1-4) were performed by Wu and Porté-Agel [30] and we made a comparison against case 3 (C T = 0.8, I 0 = 6.9%, z 0 = 0.005 m, U ∞ =9 m/s, r 0 = 40 m and z h = 70 m). Figure 11. Schematic of the vertical profile of (a) streamwise velocity and (b) velocity deficit.

Wake Profiles from the X -Y View
The X -Y views of the wind velocity at different height positions are shown in Figure 12: (a) Z = z h + 0.5d 0 plane; (b) Z = z h plane; (c) Z = z h − 0.5d 0 plane.
Because the incoming wind speed varies with the vertical coordinate, it can be seen that velocity distributions at different heights are not symmetrical to the turbine-hub plane. The wake-influenced region gradually expands as the distance to the turbine increases, while it is only severely affected near the wake center in the upstream zone (x/d 0 < 6).

Wake Profiles from the X -Z View
In Figure 13, the proposed model is compared against high-fidelity LES performed by Wu and Porté-Agel [30] in terms of the streamwise velocity in the far wake region. This region is far less influenced by the detailed features of the wind turbine (including tower, nacelle and blade profile). The result shows that the model is able to capture the characteristics of streamwise velocity field, including both velocity magnitude and wake development trend. It is reasonably deduced that wake velocity could also be estimated with acceptable accuracy in other X -Z view planes. However, it should be noted that the analytical solution (Equation (23)) will diverge in the wake region (x/d 0 < 2-3) which prevents the model's prediction near the rotor. The X -Z views of the wind velocity at different positions along the Y distance are shown in Figure 14: (a) Y = 0 plane; (b) Y = 0.25d 0 plane; (c) Y = 0.5d 0 plane; (d) Y = 0.75d 0 plane; (e) Y = 1d 0 plane. Because the wind velocity predicted by the model is symmetric about the Y = 0 plane, sections of one side are able to show the velocity distribution. In the Y = 0 plane, the wake effect is serious within 6d 0 downwind of the wind turbine. As the distance from the Y = 0 plane increases, the wake-affected zone gradually shrinks and becomes less obvious beyond Y = 0.5d 0 . Hence, when designing a wind farm, the distance between the two in-line wind turbines should not be less than 6d 0 so that the downwind turbine can be exposed to larger wind speed, thereby extracting more energy from the wind-farm flows.
It should be mentioned that the presence of the ground will hinder the wake growing downward in reality. Due to the intrinsic difficulty of analytical models assuming axial symmetry for the velocity deficit, it is nearly impossible to incorporate the ground effect into the derivation of analytical models. The image technique is applied to tackle this problem [17,48]. Based on this method, a symmetrical turbine is added below the ground that it assumed to be inviscid and the velocity deficits of both the real and the image turbines are added together, so that the drag conservation is satisfied [12]. If the image technique is employed to consider the ground effect, the normalized velocity deficit differs by not more than 0.01 at hub height for x/d 0 < 20 compared with a stand-alone turbine wake deficit without considering the ground effect. It seems that the ground has a weak impact on the wake region of our interest. Therefore, it is rational to neglect the ground effect for simplicity. result shows that a slower recovery rate is noticed in the model compared with LES; this discrepancy becomes minor in further downwind distance (X = 15d 0 ). Velocity distribution in Y -Z view obtained by the model was still in appropriate agreement with LES, and its non-axisymmetric characteristics were found because of the presence of the boundary layer. The figure shows that the wake circle is still clearly visible around 5d 0 and it gradually deforms due to the collision with the ground.
When the wind flows through a wind turbine, the torque exerted by the wind makes the turbine blades rotate. Due to the conservation of angular momentum, the wake rotates in the opposite direction to that of the blades. Given the theoretical difficulty in analytical models, there is no appropriate way to deal with the rotation effect. However, experimental [7] and numerical [49] studies found out flow rotation is not noticeable in the far wake (x/d 0 > 5). Hence, the model is still competent to predict the streamwise velocity in the far wake region without consideration of the rotation effect.

Conclusions and Further Work
In this paper, a new analytical model is proposed to predict the velocity deficit in the wake downwind of a wind turbine. A cosine-shaped distribution is considered for the velocity deficit and mass and momentum conservation are employed to derive the model, which are expressed as: where k t and I + are obtained by Equations (3) and (27) respectively. The effects of the ambient turbulence and the rotor generated turbulence are considered in the variable rate for wake expansion. Compared to the Tian model, which only employs mass conservation, this model makes use of the superiority of momentum conservation, filling a research gap in the development of cosine-series models. The proposed model is tested against field observations, wind-tunnel measurements and high-resolution LES data, and the result shows that it can predict the velocity deficit downwind of a wind turbine with acceptable accuracy in terms of both shape and magnitude. By contrast, the top-hat models generally underestimate the velocity deficit at the wake center and overestimate it in proximity to the wake edge. Based on the new model, a series of predictions from different views are presented to demonstrate the wake-influenced region. The result indicates that wake effect is inconspicuous in the very far wake region and provides some references for the selection of wind turbines' spacing. It should be noted that a major drawback of the model is that the analytical solution (Equation (23)) will diverge in the near wake region (x/d 0 < 2-3), which prevents the model's estimation of wake deficit near the rotor. However, due to its low computational cost, strong physics-based principles and good accuracy, this new model is still quite appealing for use in wind energy programs.
Future research will extend the capability of the proposed model to estimate the velocity deficit for different atmospheric stabilities. In addition, the model will be used to evaluate the power outputs of large wind farms for different layout configurations and inflow conditions. Funding: This study was supported by Chinese National Natural Science Foundation (51678452), which is gratefully acknowledged.