{'fraction_non_alphanumeric': 0.028358774680403844, 'fraction_numerical': 0.015161434823059164, 'mean_word_length': 4.65510522213562, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
On turbulence models and lidar measurements for wind turbine control

To provide comprehensive information that will assist in making decisions regarding the adoption of lidar-assisted control (LAC) in wind turbine design, this paper investigates the impact of different turbulence models on the coherence between the rotor-effective wind speed and lidar measurement. First, the differences between the Kaimal and Mann models are discussed, including the power spectrum and spatial coherence. Next, two types of lidar systems are examined to analyze the lidar measurement coherence based on commercially available lidar scan patterns. Finally, numerical simulations have been performed to compare the lidar measurement coherence for different rotor sizes. This work confirms the association between the measurement coherence and the turbulence model. The results indicate that the lidar measurement coherence with the Mann turbulence model is lower than that with the Kaimal turbulence model. In other words, the potential value creation of LAC based on simulations during the wind turbine design phase, evaluated using the Kaimal turbulence model, will be diminished if the Mann turbulence model is used instead. In particular, the difference in coherence is more significant for larger rotors. As a result, this paper suggests that the impacts of different turbulence models should be considered uncertainties while evaluating the benefits of LAC.


Introduction
Turbine-mounted lidar sensors provide preview information about the inflow wind to be used for improving wind turbine control, which is referred to as wind-turbine-integrated lidar-assisted control (LAC). LAC is a promising technology for reducing wind turbine loads and the levelized cost of energy (LCOE) (Scholbrock et al., 2016;Simley et al., 2020;Schlipf et al., 2018). The potential benefits have been demonstrated in several works by simulation (Schlipf et al., 2010;Bossanyi, 2013;Schlipf et al., 2013b;Bossanyi et al., 2014) as well as in field experiments (Kumar et al., 2015;Fleming et al., 2014;Schlipf et al., 2014).
The topic of the optimal lidar scan pattern for wind energy applications is critical for the widespread deployment of LAC. Both practical considerations for overcoming the obstacles of LAC application and for optimizing li-dar scan patterns were discussed at the International Energy Agency (IEA) Wind Task 32 workshop (Simley et al., 2018). Three commonly used simulated measurement quality metrics for LAC application are defined in Simley et al. (2018): magnitude-squared coherence between the true rotoreffective wind speed (REWS) and the lidar-based estimate, mean square error (MSE) between the true REWS and the lidar-based estimate, and MSE between the generator speed and the rated generate speed. The REWS is commonly used to indicate the rotor-averaged wind condition. The correlation between the REWS measured by the lidar and experienced by the rotor has been discussed in Haizmann et al. (2015), Simley et al. (2012) and Schlipf et al. (2013a), in which the magnitude-squared coherence is suggested as a key metric to quantify the measurement quality. A fundamental component of simulation-based lidar measurement coherence is the theoretical spatial coherence of the turbulence: (a) the lateral-vertical spatial coherence is defined in the International Electrotechnical Commission (IEC) design standard (IEC, 2019); (b) wind evolution models (Bossanyi, 2013;Simley and Pao, 2015) are defined in terms of longitudinal spatial coherence. Note that the actual coherence in the field could be different from the theoretical coherence; thus experimental validation by field testing is important as well.
According to current wind turbine design requirements in the IEC standard (IEC, 2019), for the standard wind turbine classes the turbulence model must contain the following elements: the turbulence standard deviation, the longitudinal turbulence scale parameters and a recognized model for the coherence. The standards recommend the use of either the Kaimal turbulence model, together with a standard exponential lateral-vertical coherence model, or the Mann turbulence model to represent the random wind velocity field. Although extensive research has been carried out on evaluating lidar measurement coherence, there is a clear knowledge gap regarding the impact of different turbulence models on the lidar measurement coherence. The wind field model used in most of the above-mentioned studies consists of the Kaimal turbulence spectrum and the lateral-vertical spatial coherence model defined in the IEC standard (IEC, 2019). The impact of different turbulence models on the dynamic response of an offshore wind turbine has been evaluated by Nybø et al. (2020); the results showed that as the rotor size becomes larger, the variation of the wind in time and space also becomes increasingly important. There is a need to evaluate the load reduction potential of LAC using different turbulence models, which is critical for determining the value creation of LAC during the wind turbine design phase. Held and Mann (2019) extended the previous works by Haizmann et al. (2015), Simley et al. (2012) and Schlipf et al. (2013a) to analyze lidar measurement coherence with both the Mann turbulence model and Kaimal turbulence model. The theoreticalcoherence results were compared to field data from a nacelle lidar mounted on a Vestas V52 wind turbine. The results showed that the experimental data fit better to the coherence predicted by the Mann turbulence model, and the prediction based on the Kaimal turbulence model underestimates the coherence. However, the coherence analysis focused solely on a turbine with a small rotor diameter of 52 m; the impact of different rotor sizes and lidar scan patterns on coherence have not been investigated by Held and Mann (2019).
With the advent of larger rotor sizes and more flexible wind turbines, evaluating the value creation of LAC is becoming increasingly important. The analysis in this work is based on the framework proposed by Simley et al. (2018) and Held and Mann (2019). The specific objective of this study is to investigate the impact of different turbulence models recommended by the IEC standards on the lidar measurement coherence, especially for large rotor sizes (i.e., the Technical University of Denmark -DTU -10 MW reference turbine with a rotor diameter of 178 m; Bak et al., 2013), whereby the analysis can shed light on how to reasonably evaluate LAC benefits during the wind turbine design phase. First the differences between the Kaimal and Mann models are discussed. Then two types of commercial continuouswave (CW) lidar systems are examined to analyze the lidar measurement coherence, including a 4-beam lidar and 50beam circular scan lidar. The lidar measurement model has been created based on work by Simley et al. (2011), and numerical simulations have been performed to compare the lidar measurement coherence.
The remainder of this paper is organized as follows: Sect. 2 briefly describes the different turbulence models and compares the power spectra. The lidar measurement model is established in Sect. 3. In Sect. 4, numerical simulations for different lidar scan patterns and rotor sizes are performed. Conclusions and suggestions for future work are summarized in Sect. 5.

Preliminaries and evaluation of different turbulence models
Two different turbulence models are commonly used to evaluate the design loads in the IEC standard (IEC, 2019): the Kaimal spectrum with the exponential lateral-vertical coherence model (Kaimal model) and the Mann turbulence model (Mann model). The turbulence models use similar power spectra, and the major difference is the spatial distribution of the wind velocities.

Kaimal model
The advantage of the Kaimal model is that the onedimensional spectra are expressed as simple analytic expressions. The wind disturbance is described as turbulent velocity fluctuations and is assumed to be a stationary and random vector field with zero-mean Gaussian statistics. The power spectral densities (PSDs) of each wind component are given in a non-dimensional form: where f is the frequency in hertz, while the subscript k denotes the index of the velocity component in the longitudinal u, lateral v and upward w direction, respectively. The single-sided velocity component spectrum is denoted as S k , while σ k and L k represent the standard deviation and integral length scale parameters of the velocity component, respectively. The wind speed at hub height is denoted as V hub . For the longitudinal velocity component u, σ u is the representative value of the turbulence standard deviation, and L u is defined as L u = 8.1 u . For a modern wind turbine, the hub height is typically above z ≥ 60 m, and the longitudinal length scale parameter is u = 42 m.
The cross-power spectral density (CPSD) S u i ,u j (f ) between the wind at the two spatially separated points of u i , u j can be determined from the definition of spatial cocoherence γ i,j : where S u i ,u i and S u j ,u j are the PSDs of the wind speed at two different locations, i and j . The symbol denotes the real part of a complex number. Please note that the coherence can be split into a real part and an imaginary part, which are referred to as co-coherence and quad-coherence (Nybø et al., 2020). The coherence expressed in Eq.
(2) is in the real part form.
According to the IEC standard (IEC, 2019), the following exponential coherence model can be used in conjunction with the Kaimal PSD: where r is the magnitude of the distance between the two points projected onto a plane normal to the averaged wind direction and L c = L u is the coherence scale parameter. The definition in Eq.
(3) ignores the quad-coherence; thus the wind velocity fluctuations are assumed to be in phase. This assumption may be reasonable for small rotor sizes but can be questioned for larger rotor sizes (Eliassen and Obhrai, 2016).

Mann model
The Mann turbulence model (Mann, 1994) is a spectral tensor model based on von Kármán's model, which combines rapid distortion theory (RDT) with considerations about eddy lifetimes. The RDT in the Mann model gives an equation for the evolution or the "stretching" of the spectral tensor, and the tensor will be more and more "anisotropic" with time. The RDT will finally influence the lateral-vertical coherence in the rotor plane. The three-dimensional fluctuations around the mean wind speed u(x) can be represented by the vector field whereũ is the turbulent velocity field and U (x) is the mean wind field. Because of homogeneity, the covariance tensor is a function of the separation vector r between two points and is defined as follows: where denotes ensemble averaging.
All second-order statistics of turbulence, such as variances and cross spectra, can be derived from the covariance tensor.
The spectral tensor is given by is the non-dimensional spatial wavenumber for the three component directions k = 2πf/U and U is the mean wind speed. The resulting spectral tensor components can be found in Annex C of the IEC standard (IEC, 2019). For three-dimensional turbulent velocity vector u(x), the velocity components are determined from a decomposition of the spectral tensor and an approximation by discrete Fourier transform, following the procedure detailed in Mann (1998). Compared to the Kaimal spectrum and exponential coherence model, the advantage of using the Mann model to analyze lidar measurements is that it provides a threedimensional spectral tensor. The Mann model includes correlation between the components (u, v, w), whereas the Kaimal model has no correlation between different wind components.
The Mann model is based on three adjustable parameters: α 2/3 (the Kolmogorov constant multiplied with the rate of the viscous dissipation of specific turbulent kinetic energy raised to the power of two-thirds), the length scale l and the non-dimensional parameter related to the lifetime of the eddies.
The co-coherence γ ij for spatial separations (grid points i and j ) normal to the longitudinal direction is defined as where y is the lateral separation distance and z is the vertical separation distance. For the denominator in Eq. (7), when the two indices i = j , y = z = 0 and the wavenumber auto-spectrum ii (k 1 ) and jj (k 1 ) are expressed as where the subscript xx ∈ [ii, jj ].

Evaluation using different turbulence generators
The theoretical turbulence models are quite complicated, especially for the Mann model, although the application of the Mann model only requires three parameters (α 2/3 , l and ). Therefore, numerical simulations have been performed to compare the different turbulence models in this work.

Coordinate system
The coordinate system of the wind box as well as the lidar scan patterns is shown in Fig. 1. The size of the wind box should cover the entire rotor disk. The directions of the wind components (u, v, w) are aligned with the directions of the coordinate system axes (x, y, z). The lidar scan pattern will be elaborated in Sect. 3.2.

Turbulence generator
To generate the wind boxes for further analysis, two different turbulence simulators are used. The Kaimal model can be generated using the turbulence simulator TurbSim (Jonkman and Buhl, 2006), while the Mann model is generated by HAWC2 (Horizontal Axis Wind turbine simulation Code 2nd generation; Hansen et al., 2018).
All numerical simulations are performed for a wind field with mean wind speed U = 12 m s −1 and turbulence intensity given by the IEC Class A normal turbulence model (NTM). The parameters of the three-dimensional wind box are listed in Table 1. The grid size in the vertical and lateral directions is defined by the size of the wind box L grid and number of grid points N grids . Assuming Taylor's hypothesis of frozen turbulence, the grid size along the mean wind direction is defined as U T /N x , where U is the mean wind speed, T is the total time and N x is the number of longitudinal grid points.
Since the Mann turbulence fields are normally rescaled to the specified turbulence intensity inside HAWC2, the parameter α 2/3 is chosen to be 1 and the shear parameter should be approximately 3.9 for neutral conditions. The length scale l is recommended to be l = 0.7 u for normal conditions.
The method used in TurbSim is the Veers approach (Veers, 1988), wherein the PSDs in Eq. (1) and coherence function in Eq. (3) are used to correlate the Fourier components of different points in the y-z plane. Then the inverse fast Fourier transform (IFFT) is applied to obtain the correlated time series at each grid point. Although in the IEC standard the coherence function is only applied to the u component, the Veers approach is extended to apply the coherence to the components (v, w) in this work as well. It is assumed that the spatial coherence formula presented in Eq. (3) applies to all wind components (u, v, w), and the length scales for the different components are the same as defined for the PSDs. Otherwise, without the correlation of the v and w components the coherence between the REWS and its estimated value based on lidar measurements could be unrealistically high because the contribution of the v and w components could be close to zero after spatial averaging along the lidar beams (see Sect. 3.2). In contrast, the Mann model creates a turbulence field that is fully correlated in the directions of x, y and z.

Turbulence spectrum comparison
The differences between the Mann and Kaimal models are discussed in this section. Figure 2 shows the theoretical co-coherence γ i,j at different separation distances, in which the lateral separation distance y and vertical separation distance z are selected to be 10, 30 and 50 m. Some interesting findings are the following:   Fig. 2a, wherein the lateral co-coherence reduces as the lateral separation distance increases. With the small separation distance of 10 m, the coherence with the Mann model is higher than with the Kaimal model. Conversely, with increasing separation distance, the co-coherence with the Mann model falls sharply compared with the co-coherence with the Kaimal model; the co-coherence with the Mann model is far below the co-coherence with the Kaimal model for y = 50 m. Fig. 2b, the co-coherence with the Mann model is always higher than that with the Kaimal model for low wavenumbers. Unlike the lateral co-coherence, the vertical co-coherence does not drastically decrease with increasing separation distance.

For vertical separations in
3. The co-coherence with the Mann model is negative in some frequency ranges, which is not the case for the exponential coherence model with the Kaimal model expressed in Fig. 3. This implies an opposite phase of the wind components for some frequencies. Chougule et al. (2012) investigated the vertical cross-spectral phases in neutral atmospheric flow; the work demonstrated that the phase angle of the wind component u increases with stream-wise wavenumber and vertical separation distance.
With the advent of larger rotor sizes, lidar measurements must scan a larger area upstream of the rotor. So the findings above indicate that the choice of the turbulence model strongly influences the correlation between the lidar measurement and true REWS. This impact should be considered while evaluating the benefits of LAC.

Lidar coordinate system
Two different scan patterns based on commercial nacellemounted lidars are investigated here to illustrate the impact of different turbulence models on lidar measurement coherence: a 4-beam scan pattern (Fig. 1a) and a 50-beam circular scan pattern (Fig. 1b). The lidar is mounted on the nacelle, and the scan pattern may contain many different measurement points as shown in Fig. 1. Each scan pattern is further defined by the upstream preview distance d in the x direction and radial distance r between the scan point and the hub center in the y-z plane. The lidar is assumed to be installed at the hub center for simplicity. As suggested by Simley et al. (2018), the optimal lidar scan radius and preview distance used to achieve the best representation of the actual wind variables of interest that interact with the turbine can be expressed in non-dimensional units relative to the rotor radius. Coherence bandwidth is commonly used as a lidar measurement performance metric for LAC and will be described in detail in Sect. 4.2. The optimal scan parameters for maximizing the coherence bandwidth are summarized in Table 2, and the lidar scan parameters are defined accordingly in this work. For both lidar scan patterns, the scanning frequency for completing a full scan is 1 Hz, and the line-of-sight (LOS) measurement frequency is 4 and 50 Hz based on commercial examples.

Lidar simulator
The LOS velocity at one measurement point from a lidar system can be expressed as where l = [l x , l y , l z ] denotes the unit vector in the direction that the beam is oriented and [u, v, w] denotes the wind speed vector at the measurement point. Note that the sign of the upwind direction is negative. The velocity measured by a real scanning lidar is a spatial average of the LOS velocities along the lidar beam, which is described by the range-weighting function. The rangeweighting function for continuous-wave lidars is expressed as follows (Simley et al., 2014): where F denotes the lidar focal distance, denotes the distance from the focus position along the beam direction and K N is a normalizing factor so that the integral of W L from −∞ to ∞ gives unity. R R is the Rayleigh range and is given by where λ is the laser wavelength and a 2 is the beam radius at the output lens, which is calculated for the point at which the intensity has dropped to e −2 of its value at the beam center. The lidar beam radius a 2 takes the value 28 mm, which is broadly equivalent to the beam radius for current commercial lidar products (Pena et al., 2015). The wavelength λ is assumed to be the telecommunications wavelength of 1.55 × 10 −6 m. More details regarding lidar modeling can be found in Simley et al. (2014).

Reconstruction of the rotor-effective wind speed
The REWS is modeled as a sum of the u component wind speeds across the entire rotor disk area, assuming N p points on the rotor disk: The method of reconstructing the REWS from lidar measurements has been discussed by Schlipf et al. (2011). The lidar can only measure the wind speed component along the LOS; therefore, at least three beams are needed to estimate the three-dimensional wind vector at a single point. This limitation is referred to as the cyclops dilemma (Schlipf et al., 2011). Due to the cyclops dilemma and for the purpose of collective blade pitch control, the most common assumptions for reconstructing wind speeds from lidar measurements are the following: 1. no lateral v or vertical w wind components, 2. no shears or inflow angles.
The solution for estimating the REWS from LOS measurements is given by where N denotes the number of unique beams and l x,i denotes the x component of the orientation of beam i. The wind speed estimate u lid represents the average wind speed for a lidar measuring N points upstream of the turbine.

Numerical simulation settings
In order to investigate the impact of different turbulence models on lidar measurement coherence, numerical simulations have been performed. Apart from the Vestas V52 with a 52 m rotor diameter, two other reference wind turbines are used, including the National Renewable Energy Laboratory (NREL) 5 MW reference turbine with a 126 m rotor diameter (Jonkman et al., 2009) and the DTU 10 MW reference turbine with a 178 m rotor diameter. These two rotor sizes represent typical values for onshore and offshore turbines, respectively. The numerical simulations include 18 random turbulence boxes with different seeds for each turbulence model. The simulation time is 600 s. Therefore, the combination of two types of lidars, three different rotor sizes and two turbulence models results in 12 separate scenarios and 18 random realizations for each scenario.

Criteria for evaluating measurement quality and benefits
For indicating the measurement quality, the wavenumber k at which the magnitude-squared coherence γ 2 between u lid in Eq. (13) and u eff in Eq. (12) drops below 0.5 is commonly used as a performance metric (Schlipf et al., 2013b. This metric is referred to as the coherence bandwidth k 0.5 in this work. The wavenumber is inversely proportional to the eddy size in the longitudinal direction. So, the smallest detectable eddy size measured by a lidar is defined by the wavenumber k 0.5 . In other words, the smallest detectable eddy can be interpreted as the eddy size that can be captured with a correlation of 50 %. For reducing fatigue loads using LAC, detecting eddies with a length as small as 1 D (rotor diameter) in the longitudinal direction is important, because the thrust load induced by eddies with diameters of 1 D or larger across the rotor in the lateral and vertical directions can be mitigated using collective pitch control ; in turn, eddies covering the entire rotor disk in the lateral and vertical directions are expected to extend at least 1 D in the longitudinal direction. Thus, the magnitude-squared coherence γ 2 at k = 2π/D is the most critical metric.
By optimizing the lidar scan pattern, the measurement coherence bandwidth can be maximized, but the cost of the lidar will increase as well. Meanwhile, the benefits of fatigue load reduction may reach a plateau. Generally speaking, the lower the value of k 0.5 is, the lower the LAC benefits are. Integrating LAC into the turbine design phase involves a tradeoff optimization problem to consider the turbine cost and lidar cost simultaneously.

Coherence analysis
Based on the simulation results, the magnitude-squared coherence γ 2 between the lidar measurements and REWS are presented in Fig. 3 for the different scenarios investigated. For brevity, the dash-dot line labeled as 1 D represents the wavenumber corresponding to the rotor diameter D, whereas 2 D indicates the wavenumber corresponding to two rotor diameters. It can be clearly seen that the 50-beam circular scan lidar can achieve higher measurement coherence compared to the 4-beam lidar. For the NREL 5 MW turbine and the Kaimal model (see Fig. 3c and d), the maximum coherence bandwidth k 0.5 is approximately 0.03 rad m −1 for the 4-beam scan and 0.05 rad m −1 for the 50-beam scan. These results corroborate the findings of previous work by Simley et al. (2018).
The key findings of this study are included in the following discussion. For brevity, the magnitude-squared coherence with the Mann model is represented by γ 2 Mann , and the magnitude-squared coherence with the Kaimal model is represented by γ 2 Kaimal . Corresponding theoretical-coherence curves are also included in this figure following methods described in works by Held and Mann (2019) and Schlipf et al. (2013a).
1. For the Vestas V52 turbine in Fig. 3a and b, γ 2 Mann is higher than γ 2 Kaimal in the low-wavenumber region k ≤ 0.06 rad m −1 , which aligns with the findings of the work by Held and Mann (2019), in which the authors suggested that the Kaimal model gave a slight underestimation of the measurement coherence for a 52 m rotor diameter, and the coherence predicted from the Kaimal model is lower than the coherence predicted from the Mann model. Fig. 3c and d, γ 2 Mann is slightly higher than γ 2 Kaimal for low wavenumbers. Then, the coherence starts to separate around 2 D. Specifically, γ 2 Mann decreases more sharply than γ 2 Kaimal when k exceeds 2 D.

For the NREL 5 MW turbine in
3. For the DTU 10 MW turbine in Fig. 3e and f, the trend follows the trend with the NREL 5 MW turbine, but γ 2 Mann is considerably lower than γ 2 Kaimal . The coherence of γ 2 Mann drastically decreases before 2 D. For increasing wavenumbers, larger discrepancies are noticeable between γ 2 Mann and γ 2 Kaimal . 4. The additional measurement points with the circular scan provide an obvious improvement in measurement coherence in the frequency band k ∈ [ 2π 2 D , 2π 1 D ]. The maximum coherence bandwidth k 0.5 can reach 1 D with the 50-beam circular scan. The Kaimal model indicates that the 50-beam circular scan is a better scan pattern and can lead to realizing the full potential benefits of lidar-assisted collective pitch control. Surprisingly, the maximum coherence bandwidth k 0.5 with the Mann Figure 3. Magnitude-squared coherence γ 2 between lidar measurements and the REWS. Results (a, c, e) for the 4-beam scan pattern and (b, d, f) the 50-beam circular scan pattern. From top to bottom, the plots show results for the Vestas V52 turbine, NREL 5 MW turbine and DTU 10 MW turbine. The dash-dot lines labeled 1 and 2 D represent the wavenumbers k = 2π D and k = 2π 2 D , respectively. In the legend, "theory" denotes the theoretical coherence. model is far below 1 D, which will lead to lower benefits during the wind turbine design phase.

Conclusions
A novel finding in this work is that the coherence with the Mann model is lower that with the Kaimal model for large rotors, and this difference becomes larger with increasing rotor size. Conversely, for small rotor sizes, the coherence with the Mann model is higher than that with the Kaimal model. The differences between γ 2 Mann and γ 2 Kaimal are significant. These results are in accord with the theoretical coherence shown in Fig. 2, indicating lower coherence with the Mann model for larger separation distances.
In summary, these results provide important insights into the impact of different turbulence models on lidar measurement coherence. If the wind conditions at a site agree more closely with the Mann model, the lower coherence with the Mann model will diminish the advantages of LAC because inappropriate blade pitch actions in response to the lidar measurements will deteriorate the turbine structural loading. It can therefore be suggested that the turbulence model needs to be carefully considered while integrating the LAC solution with larger-rotor turbine designs.
This work confirms the association between lidar measurement coherence, the turbulence model and rotor size. Our results suggest that this impact should be considered an uncertainty when evaluating the benefits of LAC during the wind turbine design phase. Note that the impacts on the load reduction need to be further investigated using reference turbines and aero-elastic tools following the IEC standards. More broadly, further research should be undertaken to provide guidelines on how to determine the optimal scan pattern for different site-specific atmospheric conditions and rotor sizes. Field validation is strongly recommended to mitigate the risk induced by site-specific wind conditions if LAC is adopted, especially for large rotors.
Data availability. The turbulence box data could be made available on request.