Research article 07 Jun 2021
Research article  07 Jun 2021
Wind turbine load validation in wakes using wind field reconstruction techniques and nacelle lidar wind retrievals
 ^{1}Department of Wind Energy, Technical University of Denmark, Frederiksborgvej 399, 4000 Roskilde, Denmark
 ^{2}Stuttgart Wind Energy (SWE), University of Stuttgart, Allmandring 5b, 70569 Stuttgart, Germany
 ^{1}Department of Wind Energy, Technical University of Denmark, Frederiksborgvej 399, 4000 Roskilde, Denmark
 ^{2}Stuttgart Wind Energy (SWE), University of Stuttgart, Allmandring 5b, 70569 Stuttgart, Germany
Correspondence: Davide Conti (davcon@dtu.dk)
Hide author detailsCorrespondence: Davide Conti (davcon@dtu.dk)
This study proposes two methodologies for improving the accuracy of wind turbine load assessment under wake conditions by combining nacellemounted lidar measurements with wake wind field reconstruction techniques. The first approach consists of incorporating wind measurements of the wake flow field, obtained from nacelle lidars, into random, homogeneous Gaussian turbulence fields generated using the Mann spectral tensor model. The second approach imposes wake deficit time series, which are derived by fitting a bivariate Gaussian shape function to lidar observations of the wake field, on the Mann turbulence fields. The two approaches are numerically evaluated using a virtual lidar simulator, which scans the wake flow fields generated with the dynamic wake meandering (DWM) model, i.e., the target fields. The lidarreconstructed wake fields are then input into aeroelastic simulations of the DTU 10 MW wind turbine for carrying out the load validation analysis. The power and load time series, predicted with lidarreconstructed fields, exhibit a high correlation with the corresponding target simulations, thus reducing the statistical uncertainty (realizationtorealization) inherent to engineering wake models such as the DWM model. We quantify a reduction in power and loads' statistical uncertainty by a factor of between 1.2 and 5, depending on the wind turbine component, when using lidarreconstructed fields compared to the DWM model results. Finally, we show that the number of lidarscanned points in the inflow and the size of the lidar probe volume are critical aspects for the accuracy of the reconstructed wake fields, power, and load predictions.
Wind turbines operating under wake conditions experience higher loading and lower power productions than those operating under wakefree conditions (Barthelmie et al., 2009; Larsen et al., 2013). The wakeinduced velocity deficit and its meandering are critical aspects in both loads and power analyses (Madsen et al., 2010; Doubrawa et al., 2017). The former reduces the inflow wind speed and causes unbalanced aerodynamic load distribution at the rotor, which in turn induces high load cycle amplitudes in the whole wind turbine structure (Lee et al., 2012). The latter is the main source of wakeadded turbulence (Madsen et al., 2010), affecting wind turbine responses and inducing high fatigue damage (Larsen et al., 2013). Moreover, small turbulence eddies that result from the breakdown of the tip vortices can cause small fatigue load cycles (Madsen et al., 2005). Thus, aeroelastic analysis of wind turbines operating under wake conditions requires detailed modeling of the wake flow fields.
To date, detailed predictions of wakegenerated turbulence can be achieved with large eddy simulation (LES); however, the computational cost is prohibitive when large numbers of simulations are required. This makes engineering wake models a practical alternative for certain applications. For design load evaluation, the IEC 614001 standard (IEC, 2019) recommends the dynamic wake meandering (DWM) model, among other loworder engineering wake models.
The DWM model considers wakes to act as passive tracers displaced in the lateral and vertical directions by the large eddies of the atmospheric flow (Madsen et al., 2010). The wake field is modeled as a “cascade” of quasisteady velocity deficits emitted by the source turbine that meander through a precalculated stochastic meandering path and that are advected in the streamwise direction adopting Taylor's hypothesis of frozen turbulence. These wake deficit time series are superposed on random threedimensional turbulence fields serving as input for aeroelastic simulations (Larsen et al., 2008; Madsen et al., 2010). The wake flow features simulated by the DWM model are conditional on both the ambient conditions, which can be measured from a local meteorological mast, and the operational conditions of the upstream wind turbines. In order to carry out load simulations, the 10 min statistical properties (mean and variance) of the simulated ambient inflow are set to match the measured ambient wind statistics (Dimitrov and Natarajan, 2017).
There are three primary sources of uncertainty intrinsic of engineering wake models that affect the accuracy in power and load predictions, which we here denote as the measurement, modeling, and statistical uncertainty. The measurement uncertainty includes deviations between the measured quantity of interest (e.g., the ambient wind field's characteristics or the power and load data) and its actual true values. The modeling uncertainty originates from the simplistic flowmodeling assumptions adopted to describe wake flow fields. This type of uncertainty can partly be reduced by improving the wake model (e.g., by adding further physical effects) (Keck et al., 2015) or by calibrating model parameters using measurements (Larsen et al., 2013; Reinwardt et al., 2020). Calibrating the DWM model with sitespecific observations improves the accuracy in power and load estimates; however, such calibrations do not hold at other sites (Madsen et al., 2010; Keck et al., 2012; Larsen et al., 2013; Reinwardt et al., 2020). As a result, DWMmodelbased power and load assessments might be highly uncertain at a given site unless highspatialresolution and hightemporalresolution measurements of the wake are available for model calibration.
The statistical uncertainty derives from the traditional method of performing aeroelastic simulations, for which the numerical wind fields are set to match the statistical properties (mean and variance) of the observed wind field on a 10 min basis. Since the numerical turbulence field and the wake meandering are stochastic processes, the instantaneous velocities of the simulated wake wind field and the resulting load prediction time series are uncorrelated with the observations. This can lead to simulation errors (Zwick and Muskulus, 2015) and introduces high statistical uncertainty into power and load predictions (Dimitrov and Natarajan, 2017; Pedersen et al., 2019). Further, to accurately predict power and loads in a validation analysis, it is essential to accurately reconstruct wake meandering time series.
Alternative load verification procedures are being explored to potentially reduce the statistical and modeling uncertainty in engineering wake models and replace measurements from masts with those from Doppler lidars (Dimitrov et al., 2019; Reinwardt et al., 2020; Conti et al., 2020). Lidars can provide highspatialresolution and hightemporalresolution inflow observations and extend (and eventually replace) traditional pointlike measurements such as those from cup and sonic anemometers. Further, as modern wind turbines have considerably increased in size, reaching rotor diameters of the order of 150–200 m, accurate measurements of the inflow wind field for aeroelastic calculations require multipoint and multiheight wind measurements within the entire rotor plane.
In particular, nacellemounted lidars have the advantage of being aligned with the rotor, which increases the number of validation data in contrast to a fixed mast where only a small wind direction sector is valid. The feasibility of nacellemounted lidar observations has been demonstrated for wake characterization (Trujillo et al., 2011; Fuertes et al., 2018; Herges and Keyantuo, 2019; Reinwardt et al., 2020), lidarassisted control (Schlipf et al., 2013; Simley et al., 2013, 2018), and power and load analysis in freestream conditions (Wagner et al., 2014; Dimitrov et al., 2019).
The recent work of Conti et al. (2020) proposed a lidarbased load validation procedure under wake conditions that describes wake flow fields by means of timeaveraged wind field characteristics estimated using nacelle lidar measurements. Although the quantified uncertainty in lidarbased power and load predictions was found to be comparable to estimates from IECrecommended practices that use the DWM model (Conti et al., 2020), the authors stated that lidarbased load validation procedures in wakes should account for a model of the wake deficit and its meandering dynamics to predict power and loads accurately.
Overall, developing lidarbased wake wind field reconstruction techniques that reduce the modeling and statistical uncertainties in the inflow inherent of loworder engineering wake models can improve loads and lifetime estimation accuracy (Rommel et al., 2020), enhance power curve testing in wind farms (Lydia et al., 2014; Wagner et al., 2015), and promote lidarassisted wind turbine and wind farm control strategies (Bossanyi et al., 2014; Raach et al., 2017; Simley et al., 2018; Schlipf et al., 2020).
The present work proposes two alternative approaches for wind turbine load validation under wake conditions using nacellemounted lidar retrievals combined with wake wind field reconstruction techniques. The first approach builds on the work of Dimitrov and Natarajan (2017), which incorporates multiple lidar retrievals in a turbulence field generated using the Mann spectral model (Mann, 1994) through a constrained Gaussian field algorithm. Incorporating nacellelidar measurements as constraints into turbulence fields can circumvent the DWM model's assumption to consider wakes passive tracers (Madsen et al., 2010), while reconstructing the actual observed inflow at a high spatial and temporal resolution.
The second approach reconstructs wake deficit characteristics including wake meandering by fitting a bivariate Gaussian shape function to lidar retrievals and superimposes these deficits on a random realization of the Mann turbulence field. This approach intends to minimize errors in wake deficit representations and introduce the observed wake meandering path directly into the simulations. Both lidarbased wake field reconstruction techniques can potentially decrease the modeling and statistical uncertainty inherent to the DWM model, thus predicting accurate power productions and loads.
We evaluate these lidarbased wake field reconstruction techniques in a tailored designed numerical framework that simulates a nacellemounted lidar scanning the synthetic wake fields generated with the DWM model. The main objective of this study is to verify that nacellemounted lidar measurements incorporated into wake field reconstruction methods improve the accuracy of power and load predictions when compared to wake field reconstruction using engineering wake models alone.
The work is structured as follows. In Sect. 2, we briefly formulate the load validation procedure. Section 3 introduces the methodology including the Mann spectral tensor model (Sect. 3.1) and the DWM model (Sect. 3.2). Section 3.3 describes the virtual lidar simulator and the analyzed scanning configurations. The wake field reconstruction techniques are formulated in Sect. 3.4. The results are provided in Sect. 4, including the uncertainty analysis of the lidarreconstructed fields in Sect. 4.1; a detailed analysis of the load validation results in Sect. 4.2; and the effects of the lidar specifications, e.g., probe volume size and sampling frequency, and those related to the atmospheric inflow conditions on the load prediction accuracy in Sect. 4.3. The last two sections are dedicated to the discussion of the findings and the conclusions from the study.
The design load cases (DLCs) and load verification procedure for wind turbines operating in wakes are described in the IEC standards (IEC, 2015, 2019). The present work covers the analysis of fatigue loads of wind turbines operating in wakes (see IEC 614001, DLC1.2). We apply the onetoone load validation procedure of IEC 6140013 (IEC, 2015), which consists of comparing simulated and targeted (e.g., measured) load statistics to assess the accuracy of aeroelastic simulations. As we carry out the load validation analysis numerically, we define a tailored designed load validation procedure, inspired by the approach of Dimitrov and Natarajan (2017) and illustrated in Fig. 1. The DTU 10 MW wind turbine (Bak et al., 2013) is used as reference in this study.
We use two sets of random turbulence field realizations, which we denote as set A and set B. These turbulence fields are generated using the model by Mann (1994); thus, they are defined as zeromean, homogeneous, uniformvariance Gaussian random fields. We simulate DWM modelbased wake fields using turbulence realizations from set A, which we denote as the target fields (see the rectangular black boxes in Fig. 1). In contrast, the DWM modelbased wake fields using turbulence field realizations from set B are denoted as the baseline (see the rectangular blue boxes in Fig. 1). Since the turbulence fields from set A and set B have the same turbulence characteristics, as they are generated using the same Mann parameters, but are statistically independent (i.e., the resulting wind fields time series are uncorrelated), we expect that the outcomes of load simulations with set A and set B will have the same statistical properties but will not be correlated (Dimitrov and Natarajan, 2017).
Hence, the result of a onetoone comparison of load statistics between the baseline and the target simulations is a direct measure of the statistical uncertainty (i.e., load scatter) that originates from both the random Mannbased turbulence realizations and the stochastic meandering process inherent to the DWM model. In a traditional load validation analysis, the target loads will be the measured loads, whereas the baseline loads will be the loads resulting from aeroelastic simulations using turbulence fields with the same properties as the measured inflow conditions (IEC, 2015).
To evaluate the lidarbased approaches, we use a virtual lidar simulator that scans the target wake fields and, through our proposed wake field reconstruction technique, incorporates these samples in a random turbulence field realization from set B (see Fig. 1). This numerical approach intends to imitate what we would eventually do when nacelle lidar measurements within wakes are available for load predictions.
Further, by incorporating lidar retrievals in the wind field reconstruction technique, we expect to reduce the amount of statistical uncertainty as the load time series resulting from this approach will have greater similarity with the load time series based on the target turbulence fields. Therefore, this procedure allows us to quantify the uncertainty in load predictions that results from lidarreconstructed wake fields (see the red elements in Fig. 1) against the target and, at the same time, to compare the associated statistical uncertainty with that of the baseline. To summarize, the following load simulation cases are defined:

Target. These are DWM modelbased wake fields imposed on random turbulence field realizations from set A.

Baseline. These are DWM modelbased wake fields imposed on random turbulence field realizations from set B.

Constrained simulations (CSs). These are lidarreconstructed wake fields, where lidar virtual measurements of the target fields are incorporated as constraints to random turbulence field realizations from set B.

Wake deficit simulations (WDSs). These are lidarreconstructed wake fields, where lidar virtual measurements of the target fields are fitted to a wake deficit shape function to compute wake deficits, which are then superimposed onto random turbulence field realizations from set B.
The load validation comprises a large number of simulations (we use 18 random turbulence field realizations for each individual 10 min statistic of the inflow wind) to quantify the statistical uncertainty in power and load predictions under inflow conditions measured at a site. More details on the load validation analysis are provided in Sect. 4.2. Eventually, we quantify the load uncertainties of the baseline and CS and WDS methods by comparison to the loads of the target simulations, and we define two main criteria to evaluate the proposed approaches:
 I.
The bias (here defined as the mean ratio between simulated and target loads) obtained with the lidarreconstructed CSs and WDSs is equal to that obtained with the baseline.
 II.
The statistical uncertainty (here defined as the standard deviation of the ratio between simulated and target loads) derived with the lidarreconstructed CSs and WDSs is lower than that obtained with the baseline.
Provided that these criteria are satisfied, the proposed lidarbased wake field reconstruction techniques will produce (I) power and load predictions in wakes that are statistically unbiased compared to the DWM model results and (II) a reduced statistical uncertainty in power and load predictions compared to the DWM model results, which is achieved by reconstructing wake fields with stronger similarities to the actual inflow.
3.1 Mann turbulence spectral model
The timedomain aeroelastic simulations require input of a threedimensional turbulence field that mimics atmospheric turbulence (Dimitrov et al., 2017). For this purpose, the IEC 614001 recommends, inter alia., the Mann uniform shear spectral tensor model (Mann, 1994) or the Kaimal model (Kaimal et al., 1972). The turbulence spectral properties of a threedimensional homogeneous wind field are described by the spectral velocity tensor Φ_{ij}(k) (Kristensen et al., 1989):
which is the Fourier transform of the covariance tensor R_{ij}(r); $\mathit{r}=(x,y,z)$ is the spatial separation vector defined in a righthanded coordinate system such that the longitudinal component of the wind field (u) is in the x direction, y and z are the directions of the transverse components (i.e., the v and wvelocity components), and $\mathit{k}=({k}_{\mathrm{1}},{k}_{\mathrm{2}},{k}_{\mathrm{3}})$ is the vector with the wavenumbers in the $(x,y,z)$ directions.
The model by Mann (1994) (hereafter referred to as the Mann model), assumes neutral atmospheric conditions and defines the spectral tensor as a function of three input parameters: ${\mathit{\alpha}}_{k}{\mathit{\u03f5}}^{\mathrm{2}/\mathrm{3}}$ is a product of the spectral Kolmogorov constant α_{k} and the turbulent energy dissipation rate ϵ; Γ is a parameter describing the anisotropy of the turbulence, and L is a length scale proportional to the size of turbulence eddies. From the spectral tensor, the cross spectra between two points located in a y–z plane and separated by a distance (Δ_{y},Δ_{z}) are calculated numerically by
Further, by inverse Fourier transforming the cross spectrum χ_{ij}, we can derive the auto and crosscorrelation structure of the turbulence field (Dimitrov and Natarajan, 2017) as
3.2 Dynamic wake meandering model
The DWM model is an engineering wake model that simulates wind field time series and includes three components: a quasisteady velocity deficit, the wakeadded turbulence, and the wake meandering (Madsen et al., 2010). Figure 2 illustrates these wake feature components qualitatively. The DWM model assumes wakes as passive tracers displaced in the lateral and vertical directions by the large eddies in the atmospheric flow. Further, the quasisteady wake deficits are advected in the streamwise direction adopting Taylor's assumption of frozen turbulence (Madsen et al., 2010). This set of assumptions allows decoupling the wake deficit and wakeadded turbulence components from the wake meandering model (Larsen et al., 2007). Hence, the three components of the DWM model are computed separately and subsequently superposed on random homogeneous turbulence field realizations (e.g., generated using the Mann model) to produce threedimensional wake field time series that are input into aeroelastic simulations (Larsen et al., 2013; Keck et al., 2014).
The velocity deficit definition is based on the work of Ainslie (1986, 1988), who applied a thin shearlayer approximation of the Navier–Stokes equations and a simple eddy viscosity formulation. The wake deficit expansion and recovery downstream of the generating turbine are driven by the turbulent mixing occurring due to the ambient turbulence and the turbulence generated by the wake shear field itself (Madsen et al., 2010; Keck et al., 2014, 2015). For a given wind turbine aerodynamic rotor design, a 10 min average inflow wind speed (${\overline{U}}_{\text{amb}}$), and ambient turbulence intensity (TI_{amb}), the DWM model calculates a twodimensional quasisteady velocity deficit defined in the meandering frame of reference (MFoR), which is a coordinate system with its origin in the center of symmetry of the deficit, as shown in Fig. 2a. Here, we use the numerical scheme of the standalone DWM model (Liew et al., 2020; Larsen et al., 2020) to compute the quasisteady velocity deficit.
The wakeadded turbulence originating from the breakdown of tip vortices and from the shear of the velocity deficit is accounted for by a semiempirical turbulence scaling factor (Madsen et al., 2010) as
where U_{def} is the axisymmetric velocity deficit in the MFoR (see also Fig. 2a) and k_{m1} and k_{m2} are calibration constants (Madsen et al., 2010). The twodimensional spatial distribution of k_{mt} is shown in Fig. 2b. As wake turbulence is both highly isotropic and characterized by a reduced turbulence length scale compared to that of the ambient turbulence (Madsen et al., 2005), k_{mt} of Eq. (4) scales the residual field of a Manngenerated turbulence field with a standard deviation of the longitudinal wind component equal to 1 m s^{−1} (IEC, 2019), assuming isotropic turbulence, i.e., Γ=0, and a small turbulence length scale (L≈10 %–25 % of the ambient turbulence length scale) (Madsen et al., 2010).
The wake meandering is assumed to be governed by the atmospheric turbulent structures of the order of two rotor diameters (D) or larger (Madsen et al., 2010). This assumption was verified using lidar observations of wakes (Bingöl et al., 2010; Trujillo et al., 2011). Thus, the simulated wake meandering time series is obtained by lowpass filtering atmospheric turbulence fluctuations (i.e., v and wvelocity components measured from a local mast or lidar or alternatively simulated by the Mann model) by a cutoff frequency ${f}_{\text{cut,off}}={\overline{U}}_{\text{amb}}/\left(\mathrm{2}D\right)$, which excludes contributions from smaller eddies to the meandering dynamics (Larsen et al., 2008).
As a result, the wake field simulated by the DWM model can be seen as a cascade of quasisteady velocity deficits that meander in the lateral and vertical directions and are advected downstream by the mean wind speed of the inflow using Taylor's assumption. These wake features are superposed on stochastic homogeneous turbulence field realizations to generate wake fields time series that are then input into aeroelastic simulations (see Fig. 2c).
Mathematically, a threedimensional synthetic wake flow field compliant with the DWM model formulation can be defined by a linear superposition of the ambient wind field and two inhomogeneous turbulence terms as
where ${\overline{U}}_{\text{amb}}\left(z\right)$ is the ambient wind speed including the atmospheric vertical wind shear profile, ${u}_{i{K}_{\text{def}}}^{\prime}(x,y,z)$ is a residual turbulence field with imposed wake deficits that follow the meandering path and ${u}_{j,{K}_{\text{turb}}}^{\prime}$ is a second turbulence field modeling wakeadded turbulence effects. Adopting Taylor's assumption, the wake field can be described by the spatial vector solely; thus, the time variable is disregarded in Eq. (5). The subscripts i and j indicate two random and uncorrelated turbulence field realizations. The ${u}_{i,{K}_{\text{def}}}^{\prime}$ field is computed as
where ${K}_{\text{def}}(x,y,z)$ denotes the DWM modelbased wake deficit time series including a precomputed stochastic meandering path, and ${u}_{i}^{\prime}$ is a random homogeneous turbulence field realization from the Mann model with the same Mann parameters as those of the ambient wind field. Note that ${K}_{\text{def}}(x,y,z)$ assumes values equal to unity when wake losses are not present. The Mann parameters, ${\mathit{\alpha}}_{k}{\mathit{\u03f5}}^{\mathrm{2}/\mathrm{3}}$, L, and Γ, are derived, e.g., from fitting the observed freestream turbulence velocity spectra to the Mann model with the use of precomputed lookup tables to speed up the fitting procedure (Peña et al., 2017). The wind field formulation of Eqs. (5) and (6) is consistent with the domain of wind fields typically input into aeroelastic simulations (Larsen and Hansen, 2007). Finally, ${u}_{j,{K}_{\text{turb}}}^{\prime}$ is obtained as
where ${K}_{mt}(x,y,z)$ denotes a time series of turbulence's scaling factors computed from Eq. (4) including a precomputed stochastic meandering path and ${u}_{j}^{\prime}$ is a random homogeneous turbulence field with σ_{u}=1 m s^{−1}, Γ=0, and L=10 % of the ambient turbulence length scale (Madsen et al., 2010; IEC, 2019).
3.3 Lidar simulator
We use the lidar simulator developed within the ViConDAR opensource numerical framework to virtually replicate lidar measurements (https://github.com/SWEUniStuttgart/ViConDAR, last access: 22 April 2021) (Pettas et al., 2020). The lidar simulator derives the lineofsight (LOS) velocities at each scanning location by transforming the u, v, and wvelocity components of the synthetic turbulence field into a LOS coordinate system. To simulate the probe volume of the lidars, a Gaussian weighting function W(F,r) is imposed along the LOS coordinate r and centered at the focal distance F:
The u velocity is computed from the projection of V_{LOS,eq} onto the longitudinal axis; i.e., the v and wvelocity components are neglected in the field reconstruction (Schlipf et al., 2013; Simley et al., 2013). This assumption leads to
where ϕ is the elevation and θ the azimuth angle of the scanning pattern, which refer to the rotations about the y and z axes, respectively. Neglecting the v and wvelocity components introduces uncertainty into the wind field reconstruction. However, the opening angles (ϕ,θ) relative to the scanning configurations of our work reach a maximum of 35^{∘} (see Sect. 3.3.1); thus, the errors introduced by Eq. (9) are marginal (Simley et al., 2013). Other sources of uncertainty in the radial velocity estimation inherent to lidars, e.g., from the optics and internal signal processing, are accounted for by adding Gaussian white noise. Here we add noise at a level that results in a signaltonoise ratio of −20 dB as in Pettas et al. (2020). We do not investigate the sensitivity of the noise level in the present work.
The lidar simulator can mimic any arbitrary scanning pattern and includes a time lag between each lidarsampled measurement to resemble the scanning frequency (see Fig. 3). In the present study, the virtual lidar data are computed from the synthetic wake flow fields generated using the DWM model. These wind fields are time series of the u, v, and wvelocity components defined over a turbulence box with a grid size of $\mathrm{8192}\times \mathrm{32}\times \mathrm{32}$ ($x\times y\times z$). A spatial resolution of 6.5 m is used for the grid in the rotor plane, which leads to a turbulence box with dimensions of 208 m×208 m in both lateral and vertical directions (y×z). The spatial resolution dx on the x axis depends on the simulated ambient wind speed at hub height: $\mathrm{d}x=\left({\overline{U}}_{\text{amb}}{T}_{\text{sim}}\right)/\mathrm{8192}$, where T_{sim} is the simulation time in seconds. These dimensions ensure an adequate turbulence field for a 10 min wind field simulation over a large rotor and a space–time resolution such that the probe volume effects can be captured by the virtual lidar (Dimitrov et al., 2017; Pettas et al., 2020).
3.3.1 Lidar scanning strategies
To evaluate currently available nacelle lidars' ability to perform wake characterization, we select a few standard scanning configurations and use them to perform load validation within wakes. These are a fourbeam lidar (4P) (Held and Mann, 2019a, b); an extended configuration with seven beams, six arranged at the corner of a hexagon and a central beam (7P) (Pettas et al., 2020); the conical scanning lidar (Cone) (Medley et al., 2014; Borraccino et al., 2017; Peña et al., 2017); the SpinnerLidar (SL) (Peña et al., 2019; Doubrawa et al., 2019); and a general grid pattern (Grid) covering the full turbulence box (see Fig. 4).
The lidar simulator is assumed to scan the selected patterns at the same single range upwind of the rotor. Pulsed and continuouswave (CW) lidar technologies apply different approaches at scanning multiple ranges (Peña et al., 2017). Pulsed lidars can scan multiple ranges along the LOS simultaneously within a single sample, while CW lidars typically sample much faster at a given range but need to refocus in order to change the sampling range. In the present paper, we only consider a single focusing range that is achievable with both lidar technologies. Further, a time lag between each sampling beam is simulated to mimic lidars' sampling frequency.
Although we do not optimize the scanning patterns, we use scan radii (defined as the radius between hub height and the location of the scanned points) of about 70 %–80 % of the rotor radius to estimate wind field characteristics based on previous recommendations (Dimitrov and Natarajan, 2017; Simley et al., 2018). Thus, we define the 4P, 7P, and Cone patterns accordingly as shown in Fig. 4. The SL trajectory is scaled to cover the full rotor area, and the positions over the plane of measurement are separated by 29 m in both vertical and transverse directions for the Grid pattern.
A preview distance of 0.7D (≈125 m) is assumed. Note that increasing the preview distance reduces the errors caused by the crosscontamination effects of the v and w components and reduces the induction effects but raises errors due to the wind evolution (Simley et al., 2012). These effects are not investigated in detail in this work, as we use DWM modelbased wake fields as target, which do not include induction effects or turbulence evolution as Taylor's assumption is applied.
We assume a 2 s scan period for all the simulated configurations, which refers to the time required for a beam to complete the full pattern. Given the finite resolution of the synthetic turbulence boxes (i.e., 6.5 m in both lateral and vertical directions), the Cone and SL scanned locations are binned within the box grid, as reported in Table 1.
A probe volume with an extension of 30 m in the LOS direction is assumed for all the analyzed patterns, which is comparable with the current CW lidar technology measuring at distances beyond 120 m (Peña et al., 2015). Further, a 30 m probe length is commonly used to model pulsed lidars (Schlipf, 2016). Here, we define the probe volume's length as the standard deviation of the Gaussian weighting function for convenience. Typically, Gaussian weighting functions are used to model pulsed lidars, whereas Lorentzian functions are used for CW lidars (Mann et al., 2010). Additionally, we define a case (Grid^{*}) that neglects probe volume averaging effects (see Table 1).
3.4 Wake field reconstruction techniques
By defining the DWM modelbased wake flow fields as the target fields, the underlying assumptions on which we define the lidarbased wake field reconstruction techniques are as follows:

The ambient wind conditions are known, including ${\overline{U}}_{\text{amb}}\left(z\right)$, the atmospheric turbulence intensity (TI_{amb}), and the atmospheric stability conditions (here implicitly prescribed through the Mann parameters – ${\mathit{\alpha}}_{k}{\mathit{\u03f5}}^{\mathrm{2}/\mathrm{3}}$, L, and Γ).

The lidarbased wake fields are reconstructed by incorporating lidar observations (e.g., in the form of constraints or lidarfitted velocity deficits) into a zeromean, homogeneous, and random Gaussian turbulence field generated by the Mann spectral tensor model.

The induction effects on lidar measurements are neglected, and Taylor's frozen turbulence hypothesis is assumed.

Only the uvelocity fluctuations are reconstructed from the lidar measurements of the target wake fields.
The corresponding random turbulence field realizations from set A (used for the target fields) and set B have similar spectral properties; however, these fields only describe the turbulence structures of the ambient wind field. The lidar measurements of the wake field, combined with the wake field reconstruction approach, should recover all the information regarding the wake characteristics, including velocity deficits, wakeadded turbulence, and meandering in lateral and vertical directions. Further, the first assumption is no longer needed if a second instrument is deployed at the site measuring the ambient conditions, for example, using a mast or a nacellemounted lidar (Borraccino et al., 2017; Peña et al., 2017).
The second and third assumptions are inherent in the modeling approach and limitations of the DWM model and other analytical wake models; however, in this study, the wake characteristics are extracted directly from the lidar observations rather than from a physically based deficit formulation. Eventually, wind turbine responses are mainly affected by the mean wind speed in the longitudinal direction (u velocity) and its variance (Dimitrov et al., 2018), while the effects of the v and w turbulence are generally marginal (Dimitrov and Natarajan, 2017).
3.4.1 Constrained Gaussian field simulations
The algorithm for applying constraints on a zeromean, homogeneous, and isotropic Gaussian random field was developed in Hoffman and Ribak (1991) and extended to Manngenerated turbulence fields for aeroelastic simulations in Nielsen et al. (2003) and Dimitrov and Natarajan (2017). The algorithm uses a set of constraints that are here derived from a virtual lidar simulator and an unconstrained random turbulence realization generated with the Mann spectral tensor model.
Following the notation in Dimitrov and Natarajan (2017), we denote $\stackrel{\mathrm{\u0303}}{g}\left(\mathit{r}\right)$, where $\mathit{r}=(x,y,z)$ is the spatial separation vector, an unconstrained random turbulence realization. The spectral property of $\stackrel{\mathrm{\u0303}}{g}\left(\mathit{r}\right)$ at each discrete lateral and vertical separation of the turbulence box can be computed from the Mann model in Eq. (2), given a set of parameters (${\mathit{\alpha}}_{k}{\mathit{\u03f5}}^{\mathrm{2}/\mathrm{3}},L,\mathrm{\Gamma}$). We denote a set of constraints as $\mathbf{H}=\mathit{\left\{}{h}_{i}\right(\mathit{r}){}_{{r}_{i}}={c}_{i},i,\mathrm{\dots},M\mathit{\}}$, where each constraint is a measured value of the wind speed for a particular spatial location r and M is the total number of constraints (i.e., the number of scanned points within a 10 min period). Note that the constraints are defined as a residual wind field; thus, we remove the mean ambient wind speeds from the lidar measurements of Eq. (9), i.e., ${c}_{i}\left({\mathit{r}}_{i}\right)={u}_{\text{lidar}}\left({\mathit{r}}_{i}\right){\overline{U}}_{\text{amb}}\left({\mathit{r}}_{i}\right)$, which are the values that are input into the algorithm.
The objective of the algorithm is to define a turbulence field g(r), subjected to the constraints in H that maintain the covariance and coherence properties of the unconstrained field $\stackrel{\mathrm{\u0303}}{g}\left(\mathit{r}\right)$. As demonstrated in Dimitrov and Natarajan (2017), the unknown points of the field can be defined by maximizing their conditional probability distribution on the constraint set H. We define the residual field $\mathit{\xi}\left(\mathit{r}\right)=g\left(\mathit{r}\right)\stackrel{\mathrm{\u0303}}{g}\left(\mathit{r}\right)$, which is the difference between the constrained and unconstrained fields. This residual field is also a random Gaussian field, where its values at the constraint locations are known: $\mathit{\xi}\left({\mathit{r}}_{\mathit{i}}\right)={c}_{i}\left({\mathit{r}}_{i}\right){\stackrel{\mathrm{\u0303}}{g}}_{c}\left({\mathit{r}}_{i}\right)$. The values of the residual field at unknown locations can be derived as (Dimitrov and Natarajan, 2017)
where 〈.〉 denotes ensemble averaging, ζ(r) is a vector of cross correlations between the constraints and the field, and Z is the symmetric correlation matrix of the constraints set. Both ζ(r) and Z can be computed from Eq. (3). Eventually, any constrained realization can be written as a sum of the unconstrained field and the mean of the residual field as
By denoting ${u}_{\text{CS},B,i}^{\prime}=g\left(\mathit{r}\right)$ as the constrained turbulence field that incorporates lidar measurements into a random turbulence realization i from set B (see Fig. 1), we can derive the reconstructed wake flow field to be input into aeroelastic simulations as
Note that the fidelity of the reconstructed wind field will depend on the accuracy of the nacelle lidar measurements used to characterize the wake field.
3.4.2 Wake deficit superposition simulations
The wake deficit superposition (WDS) approach assumes that velocity deficits can be described by a bivariate Gaussian shape function, which is fitted based on lidar measurements of the target wake flow field. Several studies have demonstrated the viability and robustness of the Gaussian curve fitting to track wake deficit displacements in the farwake region (Trujillo et al., 2011; Reinwardt et al., 2020).
In our study, the wake shape function not only tracks the wake meandering but also is used to quantify the depth and width of the wake at each quasiinstantaneous scan performed by the lidar. Traditionally, the normalized velocity deficit is defined as the difference between the ambient wind speed and that inside the wake as
where ${\overline{U}}_{\text{amb}}\left(z\right)$ is assumed to be known and the lidar measurements in the wake (u_{lidar}) are sampled by the lidar simulator using Eq. (9). Following the procedure of Trujillo et al. (2011), a bivariate Gaussian shape is used to describe the velocity deficit flow field as
where (μ_{y},μ_{z}) define the wake center location; (σ_{wy},σ_{wz}) are width parameters of the wake profile in the y and z directions, respectively; (y_{i},z_{i}) denote the spatial location of the LOS; and A is a scaling parameter dictating the depth of the wake. A leastsquares method is applied to fit the measured wind speed deficits from Eq. (13) to the bivariate Gaussian function in Eq. (14).
The optimal wake deficit parameters (${\mathit{\mu}}_{y},{\mathit{\mu}}_{z},{\mathit{\sigma}}_{wy},{\mathit{\sigma}}_{wz},A$) are obtained for each completed scanning period (i.e., ∼2 s as described in Table 1), resulting in approximately 300 lidarreconstructed deficits within a 10 min period. Finally, these lidarfitted wake deficits are superimposed on a random homogeneous turbulence field realization from set B, as shown in Fig. 1. A preliminary analysis showed that wide turbulence boxes (208 m×208 m) can present large turbulence structures within, i.e., broad regions across the box characterized by low wind speeds, whose sizes can alter the depth and width properties of the lidarfitted wake deficits in Eq. (14). As a result, the wake properties of the reconstructed field can considerably deviate from the actual imposed wake characteristics.
To compensate for these deviations and considering that the DWM modelbased wake fields can be defined as a linear summation of the ambient wind field ${\overline{U}}_{\text{amb}}$ scaled by the wake deficit function K_{def} and a random homogeneous turbulence realization term ${u}_{i}^{\prime}$, as reported in Eq. (6), we reformulate the leastsquares minimization problem as
where subscripts (m,n) indicate data points within the scanning configuration; the second term on the righthand side defines the velocity deficit as in Eq. (13), in which the reconstructed wake field is defined as ${\overline{U}}_{\text{amb}}(\mathrm{1}{K}_{\text{def,Gau}})+{u}_{B,i}^{\prime}$; and ${u}_{B,i}^{\prime}$ is the random homogeneous turbulence realization from set B. Note that when wake losses are present, (1−K_{def,Gau}) will reduce the ambient wind speed, as expected. As the sampling frequency of the lidar is lower than the sampling frequency of the synthetic wind field, we interpolate the fitted wake characteristics at each scan to the whole turbulence field by applying a nearestneighbor interpolation scheme. Finally, the reconstructed wake field input into aeroelastic simulations is defined by
where ${K}_{\text{def,Gau}}(x,y,z)$ is fitted using Eq. (15) for each completed scan by the nacelle lidar.
The results are divided into three parts. First, we assess the accuracy of lidarreconstructed wake fields against target fields in Sect. 4.1. Second, we carry out the load validation analysis in Sect. 4.2 and separately present the load prediction uncertainty of the CS approach in Sect. 4.2.2 and that of the WDS approach in Sect. 4.2.3. A detailed analysis of the predicted load time series and load spectral properties is conducted in Sect. 4.2.4 and 4.2.5. Finally, we evaluate the effects of both atmospheric turbulence conditions and the selected lidar technical specifications on the load prediction accuracy in Sect. 4.3.
4.1 Uncertainty in reconstructed wake fields
In this section, we evaluate the accuracy of the lidarreconstructed fields against the target fields. At first, we assess the accuracy of the reconstructed uvelocity time series across the turbulence box, by computing the root mean square error, $\text{RMSE}=\sqrt{\mathrm{1}/n{\sum}_{i}^{n}({\stackrel{\mathrm{\u0303}}{y}}_{i}{\widehat{y}}_{i}{)}^{\mathrm{2}}}/{\overline{y}}_{i}$, between the lidarreconstructed $\left(\stackrel{\mathrm{\u0303}}{y}\right)$ and target velocity $\left(\widehat{y}\right)$, where n=8192 is the grid size of the box in the longitudinal direction, normalized over the mean target velocity (${\overline{y}}_{i}$) at each grid point of the turbulence box. The normalized RMSE provides a measure of the quality of the lidarreconstructed fields with respect to the target fields; values closer to zero indicate a high precision and accuracy (see Fig. 5, top row).
Further, we compute the explained variance ratio across the turbulence box ${\mathit{\rho}}_{E}^{\mathrm{2}}=\left(\text{cov}\right(\stackrel{\mathrm{\u0303}}{y},\widehat{y})/{\mathit{\sigma}}_{\stackrel{\mathrm{\u0303}}{y}}{\mathit{\sigma}}_{\widehat{y}}{)}^{\mathrm{2}}$ (i.e., the square of Pearson's correlation coefficient; Achen, 1982), which defines the proportion of the variance in the inflow field that is transferred to the unconstrained turbulence field by imposing the constraints (Dimitrov and Natarajan, 2017). As the target and lidarreconstructed fields are based on two sets of random uncorrelated turbulence field realizations (see sets A and B in Fig. 1), ${\mathit{\rho}}_{E}^{\mathrm{2}}\sim \mathrm{0}$ is expected across the box if no lidar information is included. Contrarily, ${\mathit{\rho}}_{E}^{\mathrm{2}}=\mathrm{1}$ indicates that the reconstructed time series is fully correlated with the target; thus the two fields match completely.
Figure 5 also shows the spatial distribution of ${\mathit{\rho}}_{E}^{\mathrm{2}}$ derived from the CS and WDSreconstructed fields, with the 7P, Cone, and Grid configurations (see Table 1 for specifications). For this particular analysis, the turbine of interest is located 5D downstream of the upstream turbine, where D=179 m is the diameter of the DTU 10 MW turbine, and ambient conditions are characterized by U_{amb}=6 m s^{−1} and TI_{amb}=8 %. The inflow wind profile is defined by a powerlaw model with a shear exponent of 0.2.
As shown in Fig. 5, the locations of the imposed constraints are characterized by the lowest RMSE and highest ${\mathit{\rho}}_{E}^{\mathrm{2}}$. This effect is more pronounced for the CS results, as the algorithm imposes the actual observations directly on the synthetic field. The RMSE would tend to zero if the length of probe volume were neglected, the lidar's sampling frequency corresponded to the sampling frequency of the wind field, and crosscontamination effects were compensated for. The RMSE increases (and ${\mathit{\rho}}_{E}^{\mathrm{2}}$ decreases) for spatial regions that are farther from the lidar's beams. This occurs due to the covariance structure of the unconstrained turbulence field, for which the unknown points are nearly uncorrelated with the imposed constraints.
The errors introduced by the WDS fields are partly a consequence of an inaccurate estimation of the wake deficit characteristics (i.e., due to the limited spatial scanning configuration) and the smallscale turbulence structures contained in the turbulence box. Finally, the results in Fig. 5 confirm that the number of scanned positions by the lidar has a significant impact on the reconstructed fields' accuracy, affecting both the mean and variance of the reconstructed uvelocity component. Therefore, patterns that cover a larger region of the rotor lead to more accurate field representations (Dimitrov and Natarajan, 2017; Pettas et al., 2020).
In Fig. 6, we compare the lidarreconstructed uvelocity time series extracted at hub height, using the Grid pattern, with the target observations derived at the same location. The target wake field is simulated with U_{amb}=6 m s^{−1} and TI_{amb}=8 %. The time series of the virtual lidar measurements is also shown. We find that both field reconstruction approaches can predict the reduced wind speed within the wake region and recover the details of the wind speed fluctuations in the target field. However, uncertainty is introduced due to the limited lidar sampling frequency, the probe volume length (here assumed to be 30 m), and the adopted field reconstructing techniques. The results in Fig. 6 demonstrate that incorporating lidar data directly in the reconstructed field (i.e., the CS approach) leads to reproducing more accurate fields compared to the WDS approach.
In addition, we compute the power spectral density (PSD) of the aboveanalyzed time series of uvelocity fluctuations for a 10 min simulation and compare them in Fig. 7. We observe that the PSD of the reconstructed fields is comparable to that of the target for frequencies of up to ≈1 Hz, while the energy spectral content at higher frequencies is attenuated. According to Larsen et al. (2008), the dominant frequency of the wake meandering is defined as ${f}_{\text{cut,off}}={U}_{\text{amb}}/\left(\mathrm{2}D\right)=\mathrm{0.016}$ Hz (∼62 s period) for U_{amb}=6 m s^{−1}. As the lidar completes a full scan in 2 s, the largescale wake meandering dynamics are well captured. Further, as the wake meandering is the main source of wakeadded turbulence (i.e., ucomponent variance), the energy spectral content in the lowfrequency range is recovered, as shown in Fig. 7.
The enhanced turbulent energy content of the target field within the highfrequency range (>1 Hz) originates from the smallscale wakeadded turbulence (Madsen et al., 2010; Chamorro et al., 2012). These effects are not fully recovered in the reconstructed fields, mainly due to the lidar probe volume and limited sampling frequency and because the method fits the lidar measurements to a standard Mann turbulence field, without the smallscale wakeadded turbulence being explicitly included.
4.2 Load validation
The DTU 10 MW reference wind turbine is used for the load validation analysis (Bak et al., 2013). The load simulations are carried out using the aeroelastic code HAWC2 (Larsen and Hansen, 2007) and inflow wind conditions measured from an offshore site, as described in Sect. 4.2.1. Note that we run the analysis based on offshore wind conditions, which are characterized by low turbulence; thus, wake effects are more prominent. This work evaluates the load prediction accuracy at the main wind turbine structures, such as blades, shaft, and tower. Therefore, we neglect the modeling of the offshore substructures and foundations, and we use the onshore model of the DTU 10 MW.
Following the load validation procedure illustrated in Fig. 1, we quantify the uncertainty in power and load predictions resulting from the baseline, CSs, and WDSs against results obtained with the target fields. The CSs and WDSs are evaluated for the selected lidar configurations of Fig. 4, i.e., the 4P, 7P, Cone, SL, Grid, and Grid^{*} patterns with the parameters provided in Table 1. Two uncertainty indicators are defined to verify the load validation criteria I and II of Sect. 2:

Bias. ${\mathrm{\Delta}}_{R}=E\left(\stackrel{\mathrm{\u0303}}{y}\right)/E\left(\widehat{y}\right)$.

Uncertainty. ${X}_{R}=\sqrt{\langle (\stackrel{\mathrm{\u0303}}{y}/\widehat{y}E(\stackrel{\mathrm{\u0303}}{y})/E(\widehat{y}){)}^{\mathrm{2}}\rangle}$.
Here, the symbol E(.) denotes the mean value and 〈.〉 the ensemble average, $\widehat{y}$ is the quantity of interest (i.e., power or load statistics) derived from the target simulations, and $\stackrel{\mathrm{\u0303}}{y}$ corresponds to that produced by the reconstructed fields. We evaluate Δ_{R} and X_{R} based on the resulting 10 min power and load statistics and provide results in Sect. 4.2.2 for the CS fields and in Sect. 4.2.3 for the WDS fields.
The analyzed wind turbine responses include mean power production levels (Power_{mean}) and fatigue loads. We use the rainflowcounting algorithm to compute the 1 Hz damage equivalent fatigue loads with a Wöhler exponent of m=12 for blades and m=4 for steel structures such as tower and shaft. Thus, we compute fatigue loads at the blade root flapwise and edgewise moments, MxBR_{DEL} and MyBR_{DEL}, and towerbottom foreaft and side–side, MxTB_{DEL} and MyTB_{DEL}, and torsional loads at the tower top (also referred to as yaw moment), MzTT_{DEL}, and the drivetrain, MzSh_{DEL}.
Furthermore, we quantify the accuracy of the reconstructed wake fields based on estimates of the rotoreffective wind speed (U_{eff}), defined as the weighted sum of the u velocity measured across the rotor area, the explained variance ratio ${\mathit{\rho}}_{E}^{\mathrm{2}}$, and the uvelocity variance ${\mathit{\sigma}}_{u}^{\mathrm{2}}$ computed from the reconstructed turbulence fields. Finally, a load time series and spectral analysis is conducted in Sect. 4.2.4 and 4.2.5.
4.2.1 Site conditions
Load simulations are carried out using sitespecific observations collected from the FINO1 meteorological mast installed at the German offshore wind farm Alpha Ventus. The wind farm is situated in the North Sea and about 45 km north of the island of Borkum (Kretschmer et al., 2019). Data were collected over a period of 3 years from 2011 to 2014, and their details can be found in Kretschmer et al. (2019).
In the present work, we only use wind speeds and turbulence intensities measured under nearneutral conditions from a 90 m sonic anemometer installed at the mast. We extract 10 min average turbulence values binned for wind speeds ranging between 6 and 22 m s^{−1}; using wind speed bins with a 2 m s^{−1} bin width, we obtain nine bins with turbulence intensities of 8, 7, 7, 6, 6, 6, 6, 5, and 5 %, respectively. These are the statistics of the ambient wind field that we use as inputs for the load validation analysis.
For each 10 min sample of the inflow wind, we use 18 turbulence field realizations (IEC 614001 recommends at least 6 realizations), leading to 162 aeroelastic simulations for each scanning configuration analyzed. Simulations with ambient wind speeds below 6 m s^{−1} are disregarded, as the wind speed approaching the rotor drops below the turbine's cutin threshold due to wake deficit effects and the turbine shuts down.
Note that the recorded turbulence estimates at Alpha Ventus are considerably lower (by approximately a factor of 3) than values recommended by the lowturbulence IEC class C. Here, we perform the load validation analysis on more realistic turbulence estimates characterizing offshore sites, since IEC classC conditions would significantly attenuate the wakeinduced effects, as higher ambient turbulence leads to a faster recovery of the wake deficit.
We use standard IECrecommended turbulence parameters for the Mann model (i.e., L=29.4 m and Γ=3.9; IEC, 2019), whereas ${\mathit{\alpha}}_{k}{\mathit{\u03f5}}^{\mathrm{2}/\mathrm{3}}$ is tuned to obtain the target ambient turbulence levels of each simulation. The inflow is described by a power law with a fixed shear exponent of α=0.2, as recommended in the IEC standard. The spacing between the analyzed and upstream turbines is fixed at 5D.
The target wake field characteristics as a function of the ambient wind speed, which result from the 162 simulations, are shown in Fig. 8. The wake considerably reduces the inflow wind speed approaching the rotor (i.e., U_{eff}) by approximately 35 %, compared to the ambient wind speed (see Fig. 8a). This effect decreases for higher winds (>14 m s^{−1}) due to the low thrust coefficients of the turbine. However, the wake deficit does not fully recover at high wind speeds (${U}_{\text{eff}}/{U}_{\text{amb}}\approx \mathrm{0.93}$), as we simulate relatively low ambient turbulence levels, the spacing between the turbines is short (i.e., 5D), and the thrust coefficient of the turbine is nearly constant at high wind speeds. Further, the contribution of the wind shear to the ratio ${U}_{\text{eff}}/{U}_{\text{amb}}$ accounts for up to 3 % in freestream conditions; i.e., ${U}_{\text{eff}}/{U}_{\text{amb}}\approx \mathrm{0.97}.$
Turbulence levels within the wake region are nearly doubled at low wind speeds compared to the ambient conditions, as shown in Fig. 8b. Further, the wake meandering amplitudes, here computed as the standard deviation of the wake center displacements in the transverse directions normalized with the rotor diameter (${\mathit{\sigma}}_{{\mathit{\mu}}_{y}}/D$ and ${\mathit{\sigma}}_{{\mathit{\mu}}_{z}}/D$), are also shown in Fig. 8c and d. As expected, larger wake displacements occur in the lateral direction compared to in the vertical direction (Keck et al., 2014; Machefaux et al., 2016).
4.2.2 Load uncertainty of constrained simulations
The uncertainties (Δ_{R} and X_{R}) of load predictions obtained with the CS fields as a function of the ambient wind speed are shown in Fig. 9. We find that the biases largely vary depending on the simulated scanning pattern and analyzed load sensor. First, we observe that the patterns with fewer “points” (i.e., 4P, 7P, and Cone) overestimate U_{eff} by 2 %–10 % (see Fig. 9a). This is because (1) these patterns scan an insufficient number of positions within the inflow area to characterize the wake flow fully and (2) the autocorrelation structure of the unconstrained turbulence box is such that the spatial regions that are not scanned by the lidar are nearly uncorrelated with the locations of the imposed constraints, as also shown in Fig. 5. Thus, in the regions that are not scanned by the lidar, the reconstructed wind speed approaches the ambient wind speed values.
As a result, lower deficits are simulated or equivalently higher rotoreffective wind speeds are predicted. Consequently, the power predictions are overestimated (Δ_{R}>10 %), as seen for ambient wind speeds below 14 m s^{−1} in Fig. 9b. Patterns with high spatial resolution, such as those of SL, Grid, and Grid^{*}, provide rotoreffective wind speed and power production estimates in good agreement with the baseline.
The statistics of ${\mathit{\rho}}_{E}^{\mathrm{2}}$ in Fig. 9c indicate that increasing the number of points scanned by the lidar (see SL, Grid, and Grid^{*}) leads to a more accurate reconstruction of the wake turbulence. The biases of both ${\mathit{\rho}}_{E}^{\mathrm{2}}$ and U_{eff} decrease for high wind speeds due to the attenuated wakeinduced effects (see Fig. 8). The improved performance of SL and Grid (and Grid^{*}) is also confirmed by estimates of ${\mathit{\sigma}}_{u}^{\mathrm{2}}$ in Fig. 9d, which show that the SL and Grid configurations can match the target variance with an accuracy of up to 98 %, compared to 40 %–60 % accuracy inherent of the 4P, 7P, and the Cone configurations. Nevertheless, the observed biases of U_{eff}, ${\mathit{\rho}}_{E}^{\mathrm{2}}$, and ${\mathit{\sigma}}_{u}^{\mathrm{2}}$ reveal that the 4P, 7P, and Cone patterns lead to inaccurate wake field representations and do not satisfy the criteria of the load validation (see Criterion I in Sect. 2).
The results from simulations with the SL, Grid, and Grid^{*} patterns provide fatigue load statistics of MxBR_{DEL}, MxTB_{DEL}, MzTT_{DEL}, and MzSh_{DEL} in good agreement with the results of the baseline (see Fig. 9e–h). However, the calculated biases indicate a consistent underprediction at all wind speeds. This gap is largely compensated for when probe volume effects are neglected, as seen for Grid^{*} (green lines). Overall, the observed deviations in the load predictions are due to the uncertainty in lidar measurements (i.e., size of the probe volume, crosscontamination effects, limited sampling frequency) and the limited scanning coverage of the patterns.
Figure 10 shows the statistics of Δ_{R} and X_{R} including all wind speeds. As expected, the baseline leads to Δ_{R}∼1 for all the analyzed load sensors, which indicates that the adopted 18 turbulence seeds are sufficient for the load statistics to converge. The large biases from simulations with the 4P, 7P, and Cone patterns (Δ_{R}∼0.87–1.37) follow from the inaccurate wind field reconstruction discussed above. The load predictions with the SL and the Grid configurations provide biases closer to the baseline, although turbulencedriven load sensors are underpredicted by 2 %–7 %. These deviations decrease as probe volume effects are neglected (i.e., Δ_{R}∼1 % for Grid^{*} in Fig. 10, left).
The statistics of X_{R} are shown in Fig. 10, right. The baseline's X_{R} is a direct measure of the statistical uncertainty intrinsic to the DWM model, which is due to the stochastic properties of the synthetic turbulence field and wake meandering. Thus, the turbine responses largely affected by wakeinduced effects are identified by high X_{R} values (see baseline in Fig. 10, right). The power predictions and the majority of fatigue loads show a relatively high statistical uncertainty (X_{R}∼0.05–0.09), resulting in a large load scatter.
The X_{R} values of MyTB and MzSh are significantly higher than other load sensors. The cause of the former is structural resonance occurring at low wind speeds for which the 3P frequency (≈0.3 Hz) excites the tower's natural frequency (≈0.25 Hz) (Bak et al., 2013). This effect originates from a design aspect of the DTU 10 MW turbine and is amplified under wake conditions due to the unbalanced aerodynamic load distribution induced at the rotor. Nevertheless, structural resonance is independent of the wake field reconstruction approach. The high X_{R} values of MzSh originate from the intense controller activity to regulate the generator torque under highly variable inflow conditions. Significantly lower X_{R} values characterize the CSbased load predictions compared to the statistics obtained with the baseline. X_{R} values are reduced by a factor of between 1.4–5 for the main winddriven turbine responses such as Power_{mean} and fatigue loads (i.e., MxBR_{DEL}, MxTB_{DEL}, MzTT_{DEL}, and MzSh_{DEL}). The CS fields, reconstructed using scanning patterns with a sufficient number of scanned positions and limited lidar probe volume, satisfy both the load validation criteria I and II in Sect. 2.
4.2.3 Load uncertainty of wake deficit simulations
We present the results relative to the WDSs in the same fashion as for the CSs in Sect. 4.2.2. Thus, we plot the load prediction uncertainty as a function of the ambient wind speeds in Fig. 11. The 4P, 7P, and Cone patterns lead to improved biases of U_{eff} and consequently Power_{mean} (see Fig. 11a and b) compared to the results obtained with the CS fields (shown in Fig. 10). The Power_{mean} predictions computed with the WDS approach and the 7P pattern are comparable with the baseline, while the results using the CS fields overpredicted it by ≈10 %. Also, improved estimates of both ${\mathit{\rho}}_{E}^{\mathrm{2}}$ and ${\mathit{\sigma}}_{u}^{\mathrm{2}}$ are found in Fig. 11c and d, including for low wind speeds, which indicates a more accurate reconstruction of the wake turbulence by the WDS than by the CS approach. In contrast, we find lower values of ${\mathit{\rho}}_{E}^{\mathrm{2}}$ and ${\mathit{\sigma}}_{u}^{\mathrm{2}}$ under higher wind speeds compared to lower wind speeds because of the less pronounced wake deficits of the target fields, as also shown in Fig. 8a.
These findings suggest that more details on the wake characteristics are better recovered by fitting a wake deficit function rather than by incorporating lidar measurements directly into the turbulence boxes, when looking at patterns where the inflow is scanned at few positions. Overall, simulations with the 7P, SL, Grid, and Grid^{*} patterns can produce power predictions comparable with the baseline (see Fig. 11b), whereas the 4P and Cone patterns lead to inaccurate predictions. Figure 11e–h show that the fatigue loads obtained with the 7P, SL, Grid, and Grid^{*} configurations are generally lower than those from the baseline.
We quantify the statistics of Δ_{R} and X_{R}, including all the wind speeds using WDS fields, and present the results in Fig. 12. As discussed above, the 4P and Cone patterns overpredict the rotoreffective wind speed and underpredict the wake turbulence; these effects counteract each other, leading to fictitious lack of biases in fatigue loads. Similar conclusions can be made for the 7P configuration, although it provides reliable power estimates. As seen for the CS results, the SL, Grid, and Grid^{*} configurations provide biases in good agreement with the baseline, although fatigue loads are underpredicted by up to 6 %. By neglecting volumeaveraging effects (i.e., Grid^{*}), only a marginal improvement of the biases is achieved. Simulations with the WDS fields can reduce the statistical uncertainty in Power_{mean} by a factor of 5 and that of the main load components (i.e., MxBR_{DEL}, MxTB_{DEL}, MzTT_{DEL}, and MzSh_{DEL}) by a factor of 1.2–2 compared to the baseline (see X_{R} in Fig. 12, right).
4.2.4 Timeseries analysis of load predictions
In this section, we investigate the accuracy of lidarreconstructed load time series against target observations. An illustrative example is provided in Fig. 13, where the lidarbased power and load timeseries predictions are compared with the target simulations. As shown, both CS and WDS approaches recover to a large extent wakeinduced effects on and the instantaneous events of the wind turbine responses, leading to load time series that are highly correlated with the target observations. This finding explains the reductions in X_{R} observed in Figs. 10 and 12.
In order to quantify the accuracy of the predicted load time series, we evaluate the cross correlations $\mathit{\rho}(\stackrel{\mathrm{\u0303}}{y},\widehat{y})=\text{cov}(\stackrel{\mathrm{\u0303}}{y},\widehat{y})/{\mathit{\sigma}}_{\stackrel{\mathrm{\u0303}}{y}}{\mathit{\sigma}}_{\widehat{y}}$ between the lidarbased results $\left(\stackrel{\mathrm{\u0303}}{y}\right)$ and the target simulations $\left(\widehat{y}\right)$ (ρ=1 means perfect correlation). We focus the analysis on the SL, Grid, and Grid^{*} configurations, which provide the most promising results, as demonstrated in the previous sections. We compute ρ for all the 162 simulations and for each load component, and provide average estimates in Fig. 14. We find that both the CS and WDSpredicted Power time series reach a nearly perfect correlation with the actual target observations (ρ=0.96–0.99). Note that Power is a lowfrequency signal (see Fig. 13), which is marginally affected by the local turbulence fluctuations. A high correlation value is also obtained for MxBR (ρ=0.89–0.98) and for the tower top and shaft load components (ρ=0.60–0.90).
The correlation relative to MxTB drops to ≈0.33 with the WDSs, while higher values are achieved by the CS results. It should be noted that the structural resonance occurring at low wind speeds, which excites the tower, can potentially affect the correlation results. Figure 13 shows that the MxTB time series presents a nearly periodic signal, where the wind turbulence imprint is marginal. Overall, the accuracy of lidarreconstructed load time series shows a significantly higher degree of correlation with the target observations compared to that achieved by the baseline. Furthermore, the CS approach can predict more accurately the observed load fluctuations compared to the WDS approach.
4.2.5 Spectral coherence analysis of load predictions
We conduct a spectral analysis on the time series of MxBR, MxTB, and MzTT, which are highly correlated with the wake meandering (Muller et al., 2015; Moens et al., 2019; Ning and Wan, 2019) and are primarily affected by wake turbulence. The PSD analysis is provided in Appendix A and shows that neither of the wake field reconstruction methods shifts the energy content among frequencies or introduce instabilities (i.e., artificial artifacts).
The spectral coherence analysis provides more insight into the accuracy of reconstructed blade and tower loads. Here, we compute the coherence as ${\mathit{\gamma}}^{\mathrm{2}}=\leftS\right(\stackrel{\mathrm{\u0303}}{y},\widehat{y}\left)\right(f){}^{\mathrm{2}}/(S\left(\stackrel{\mathrm{\u0303}}{y}\right)\left(f\right)S\left(\widehat{y}\right)\left(f\right))$, where $S\left(\stackrel{\mathrm{\u0303}}{y}\right)$ and $S\left(\widehat{y}\right)$ are the autospectra of the CS (or WDS) and target load estimates and $S(\stackrel{\mathrm{\u0303}}{y},\widehat{y})$ is their cross spectrum. We compare the coherence resulting from the load time series produced by either CSs or WDSs with the target observations for U_{amb}=6 m s^{−1} and TI_{amb}=8 % (see Fig. 15).
It is observed that both field reconstruction techniques lead to high coherence in the proximity of the principal load frequencies, such as the rotational frequency (1P for the blade and 3P for the tower; see Fig. A1 for more details), the natural frequency of the tower (≈0.25 Hz, which is close to that of 3P at 6 m s^{−1}), and the dominant wake meandering frequency (≈0.016 Hz). In general, the coherence from the CSs is nonzero at frequencies of up to 0.7 Hz (6P) and is higher than that from WDSs. This confirms that higherfrequency fluctuations can be reconstructed more accurately using the CS approach.
By increasing the scanning pattern's temporal resolution and the number of scanned points and by neglecting volumeaveraging effects, the CS approach could potentially reconstruct the whole spectrum of the loads. With the WDS approach, we can only reconstruct turbulence structures corresponding to the size of the wake deficit. Finally, given the limitation of the reconstruction techniques to recover smallscale turbulence structures, as discussed in Fig. 7, the accuracy of tower loads, which are driven by highfrequency fluctuations (see Fig. A1), is lower compared to that of the blades. This can partly explain the larger deviations of Δ_{R}, X_{R}, and ρ, inherent in MxTB_{DEL} and relative to MxBR_{DEL}, observed in Figs. 10, 12, and 14, as well as explain why Δ_{R} for MxTB_{DEL} improves the most when the probe volume size is neglected, as seen in Fig. 10.
4.3 Sensitivity analysis
The load validation of Sect. 4.2 is carried out using statistics collected under nearneutral conditions at Alpha Ventus, i.e., low atmospheric turbulence. Nevertheless, atmospheric turbulence conditions have a strong impact on the wake development (Kumer et al., 2017; Zhan et al., 2020) and wind turbine loads (Sathe et al., 2013; Kretschmer et al., 2018). Further, the lidar measuring characteristics can impact the accuracy of reconstructed fields (Lundquist et al., 2015) and thus that of load predictions. In the next subsections, we investigate the effect of atmospheric turbulence conditions as well as selected lidar specifications on the accuracy of lidarbased load predictions using the Grid pattern as an example.
4.3.1 Effect of atmospheric turbulence conditions on load prediction accuracy
Figure 16a shows the effect of the lidarbased load predictions bias on TI_{amb} within the range of 4 %–20 %. The high TI_{amb} leads to faster recovery of the velocity deficit (Doubrawa et al., 2019), amplifies the wake meandering (Machefaux et al., 2016), and affects the accuracy of lidarreconstructed fields (Pettas et al., 2020). This has a negligible effect on the accuracy of load predictions obtained with the CS fields, while larger deviations are observed for the WDS results. This is partly due to the limited area scanned by the lidar combined with the large wake displacements. The fitting procedure intrinsic to the WDS approach can lead to an inaccurate estimation of the wake shape parameters when the wake moves out of the scanned area (Trujillo et al., 2011).
We investigate the influence of the atmospheric turbulence length scale on the load prediction accuracy in Fig. 16b, by varying L between 5 and 70 m. Earlier studies have shown the strong dependency of load statistics on the turbulence length scales (Sathe et al., 2013; Dimitrov et al., 2017; Conti et al., 2020). Further, L provides a measure of the scanning configuration resolution useful for performing constraints (Dimitrov and Natarajan, 2017). The turbulence length scale affects the predicted statistics of the explained variance ratio of the CS fields, which decreases from ${\mathit{\rho}}_{E}^{\mathrm{2}}\sim \mathrm{0.8}$ for L=29 m to ${\mathit{\rho}}_{E}^{\mathrm{2}}\sim \mathrm{0.6}$ for L=5 m (not shown). This indicates that when L is low, the turbulence structure sizes fall below the sampling fidelity of the CS approach (note that the scanned points of the Grid configuration are separated by 29 m as described in Sect. 3.3). Also, the turbulence structure sizes become small relative to the lidar probe volume, causing the lidar measurements to average out more of the turbulence. The CSbased load predictions' biases show a dependency on the turbulence length scales, while the WDS fields are not significantly affected (see Fig. 16b).
4.3.2 Effect of lidar probe volume and scanning period on load prediction accuracy
One of the main limitations of continuouswave lidars is that the probe volume size increases proportionally with the square of the focal distance (Sathe and Mann, 2013). As the diameter of modern wind turbines has reached 150–200 m, measuring at farther distances upstream of the rotor becomes an issue due to the larger probe volumes. Hence, we investigate the effect of the lidar probe volume on the load prediction accuracy in Fig. 17a, by varying the probe volume length between 0 to 210 m. As shown, the magnitude of Δ_{R} decreases almost linearly with increasing probe volume lengths. Further, the probe volume effects are more pronounced for the CS approach, which directly incorporates the lowpass filtered wind speed fluctuations into the reconstructed field.
Another limitation inherent in the pulsed lidar technology is the reduced sampling frequency compared to continuouswave lidars (Peña et al., 2015). The lidar sampling frequency effect on the load prediction accuracy is assessed by varying the scanning period, which is defined as the time to complete a full scan (1–30 s). For this particular analysis, the target simulations are run for U_{amb}=6 m s^{−1} and TI_{amb}=16 %. Although the scanning period does not make an important contribution to the load prediction accuracy, as shown in Fig. 17b, this outcome is conditional on the dominant frequency of the wake meandering, which in turn decreases with larger rotors (${f}_{\text{cut,out}}={U}_{\text{amb}}/\left(\mathrm{2}D\right)$), and the number of scanned points in the inflow. The CS results show that a bias lower than 2 % in power predictions is found for scanning periods of up to ≈20 s, which corresponds to onethird of the wake meandering dominant period.
One of the main elements used in the study is to regard as target the wake flow fields generated by the DWM model. The DWM model is a simplified engineering wake model subjected to modeling uncertainties. Although the mean wind velocity and turbulence fields in the farwake region can deviate from highfidelity simulations (e.g., computational fluid dynamics, CFD) or field data, the calibration of the DWM model coefficients can considerably improve the accuracy and provide wake fields in good agreement with lidar observations (Reinwardt et al., 2020) and CFD simulations (Keck et al., 2012, 2014, 2015). The modeling uncertainty originating from an inaccurate calibration of the DWM model is not expected to significantly alter this study's findings, as we demonstrate the robustness of the lidarbased approaches under a large variety of inflow wind and operational conditions.
The wake turbulence spectral properties are described, to the extent needed for the load analysis, by an isotropic Manngenerated turbulence field with a low length scale (Madsen et al., 2005). A more realistic modeling choice to accurately simulate the turbulence structures within the wake fields, which can also affect aeroelastic load simulations, is found in LES (Churchfield et al., 2015; Nebenführ and Davidson, 2017). Further, lidarbased wind field reconstruction techniques applied to LES fields have been recently developed (Bauweraerts and Meyers, 2020, 2021). Nevertheless, the computational burden of highfidelity simulations, such as LES, would make the statistical load analysis of this work unfeasible.
Another limitation stems from the lidar simulator used in the study, which replaces fullfield lidar measurements. Real lidar data taken upstream of the rotor should be corrected for induction (Borraccino et al., 2017; Mann et al., 2018), blade blockage effects, and wind evolution (Bossanyi, 2013; de Mare and Mann, 2016). These effects are not simulated due to the modeling assumptions of the DWM model and should be further investigated, e.g., using LES fields. Despite the limitations mentioned above, the numerical framework developed within this work is useful to assess the influence of several uncertainty sources on power and load predictions and to evaluate different lidar scanning strategies in an idealized yet fully controllable environment.
Characterizing the smallscale wakeadded turbulence poses a challenge given the limitations of the lidar's sampling frequency and probe volume size (Peña et al., 2017). The smallscale wakeadded turbulence enhances the energy spectral content in the highfrequency range, 0.4–20 Hz (Madsen et al., 2010; Chamorro et al., 2012; Singh et al., 2014), and its contribution to the fatigue damage varies according to the load component and turbine operational strategy (Tibaldi et al., 2015). The CS algorithm could potentially reconstruct highfrequency wakeadded turbulence. However, this would require sampling the wind field at a high temporal frequency and without probe volume effects. In contrast, the WDS method cannot explicitly reconstruct the highfrequency wakeadded turbulence.
Bergami and Gaunaa (2014) demonstrated that the most serious fatigue damage on the blades occurs at frequencies of around 1P (0.1–0.16 Hz for the DTU 10 MW), whereas structures such as the tower top (nacelle) and tower bottom are mainly affected by the tower eigenfrequency (≈0.25 Hz) and the 3P frequency (0.3–0.48 Hz). As the PSD of tower loads exhibits large energy spectral content at high frequencies (see Fig. A1), the accuracy of tower load predictions decreases compared to that achieved by blade loads, as found in Sect. 4.2.4 and 4.2.5.
We demonstrate that a high number of lidarscanned positions of the inflow are required to ensure an acceptable level of accuracy in the reconstructed wake fields. The results reveal that the current commercially available nacellemounted lidars (e.g., the 4P, 7P, and Cone patterns) will not provide sufficient information to reconstruct the wake fields accurately for the load assessments. In contrast, the scanning requirements are fulfilled by the SpinnerLidar and any arbitrary lidar that can potentially scan a greater region of the rotor, e.g., a Gridlike configuration. Although we do not optimize the scanning strategies, it is inferred that the required number of positions scanned by the lidar depends on the size of the turbulence structures in the wake field.
Incorporating a sufficient number of lidar measurements directly in the turbulence fields leads to more accurate load predictions than assuming a wake deficit's generic shape function. The CS algorithm can also be extended to reconstruct the v and wturbulence fluctuations (Dimitrov and Natarajan, 2017). Additionally, the CS method finds direct application for reconstructing more complex flow fields occurring in wind farms, e.g., multiple wakes.
On the other hand, the accuracy of the WDSpredicted loads is conditional on the selected shape function's capability to represent velocity deficits. The wake deficit can deviate from a Gaussian shape as the atmosphere becomes more unstable (Ning and Wan, 2019); it exhibits a doublepeak shape in the nearwake region (Keck et al., 2014) and a more complex geometry in a multiplewake scenario. Overall, reproducing the actual observed wake meandering path in the wake field simulations can potentially reduce the statistical uncertainty in power and load predictions.
The fitting procedure of the WDS approach is relatively fast and can provide realtime spatial and temporal characteristics of the wake flow field, which are useful for power and load predictions, wind farm monitoring, and control strategies. The computational cost of the CS algorithm considerably increases with the number of constraints simulated and the dimension of the turbulence boxes. For reference, a single wind field with 27 900 constraints (i.e., using the SL configuration) and a turbulence box with a grid size of $\mathrm{8192}\times \mathrm{32}\times \mathrm{32}$ points currently requires 1.5 h of simulation time on a single CPU.
This study proposed two alternative wind turbine load validation procedures under wake conditions that reconstruct synthetic wake fields from time series of lidar retrievals. The first approach consisted of incorporating nacelle lidar measurements of the wake as constraints into random Mann turbulence field realizations. The second approach relied on the superposition of lidarfitted bivariate Gaussian wake deficit time series on the Mann turbulence fields. The two approaches were numerically evaluated, adopting a tailored designed framework that uses a virtual lidar simulator to scan threedimensional wake fields simulated by the DWM model (i.e., the target fields).
We demonstrated that lidarreconstructed wake fields recovered the main wake flow features affecting wind turbine power and load predictions, such as the spatial distribution of the velocity deficit and its meandering dynamics. However, the accuracy of power and load estimates was highly conditional on the number of scanned points by the lidar, the probe volume size, and the ambient turbulence intensity that in turn affected the wake evolution.
The load validation analysis showed that the current commercially available nacellemounted lidars would not provide sufficient spatial resolution to characterize wakes for power and load assessments, whereas research lidars, e.g., the SpinnerLidar and the Gridlike configuration, fulfilled these requirements.
Provided that a sufficient number of wind measurements were taken upwind of the rotor (e.g., using the SpinnerLidar or the Grid), incorporating them as constraints into turbulence fields was the most robust and accurate procedure for reconstructing wake fields and predicting power and loads. The lidarreconstructed wake fields produced power and load time series that were highly correlated with the target turbine responses; thus, reducing the statistical uncertainty (realizationtorealization) by a factor of 1.2–5 when compared to the traditional load validation procedure (i.e., using the DWM model). Although unbiased power productions were predicted, the SpinnerLidar and Gridbased reconstructed wake fields underpredicted fatigue load estimates by 1 %–8 % depending on the load component and the size of the probe volume. The biases in fatigue load predictions were further reduced to less than 2 % when neglecting probe volume effects.
Further investigations should evaluate the effects of rotor induction and turbulence evolution on the accuracy of lidarreconstructed wake fields. Besides, the proposed wake field reconstruction techniques should be validated using fullfield data collected in operating wind farms.
Figure A1 shows a comparison of the PSD of MxBR, MxTB, and MzTT between the lidarreconstructed and target simulations for U_{amb}=6 m s^{−1} and TI_{amb}=8 %. Figure A1a displays the PSD of MxBR, where the first three peaks correspond to the subsequent rotor harmonics (1P, 2P, and 3P). The highest observed peak is at 1P (∼0.1 Hz), which indicates that the greatest load cycle amplitude is due to asymmetric blade loading conditions. This effect is amplified by the inhomogeneous wake field approaching the rotor.
Compared to the rotating blades, the PSD of the tower loads MxTB and MzTT exhibits the largest energy content at higher frequencies of up to 3P (∼0.3 Hz). Further, the natural frequency of the tower (0.25 Hz) nearly corresponds to the 3P frequency at 6 m s^{−1}. This explains the very high peak seen for MxTB. Overall, the PSD produced by the simulations with lidarreconstructed fields (CSs and WDSs) shows good agreement with that of the target simulations, meaning that the energy content is not being shifted between frequencies. However, it is observed that the energy content at high frequencies (>1 Hz), induced by the wakeadded turbulence, is not fully recovered due to the lidar probe volume and limited sampling frequency.
The ViConDAR framework can be accessed at https://github.com/SWEUniStuttgart/ViConDAR (https://doi.org/10.5281/zenodo.4889772, Vaspe, 2021), and the constrained simulation framework can be accessed at https://gitlab.windenergy.dtu.dk/nkdi/constrainedsimulation (https://doi.org/10.5281/zenodo.4896514, Conti and Dimitrov, 2021).
DC, VP, ND, and AP participated in the conception and design of the work. VP developed the ViConDAR framework. ND developed the constrained simulation framework. DC developed the numerical framework that gets ViConDAR to scan wakegenerated fields and use virtual lidar measurements as inputs to reconstruct wake field via CS and WDS approaches. DC also carried out aeroelastic simulations and wrote the draft manuscript. VP, ND, and AP provided key elements of the programming code, supported the overall analysis, and critically revised the manuscript.
The authors declare that they have no conflict of interest.
Vasilis Pettas is partially funded by the German Federal Ministry for Economic Affairs and Energy (BMWi) in the framework of the national joint research project RAVE – OWP Control (ref. 0324131B).
This paper was edited by Raúl Bayoán Cal and reviewed by two anonymous referees.
IEC: International Standard IEC6140013: Wind turbines – Part 13: Measurement of mechanical loads, Standard, IEC, 2015. a, b, c
IEC: International Standard IEC614001: wind turbines – Part 1: design guidelines, Fourth; 2019, Standard, IEC, 2019. a, b, c, d, e
Achen, C. H.: Interpreting and Using Regression, Sage Publications, Beverly Hills, https://doi.org/10.4135/9781412984560, 1982. a
Ainslie, J.: Calculating the flow field in the wake of wind turbines, J. Wind Eng. Ind. Aerod., 27, 213–224, https://doi.org/10.1016/01676105(88)900372, 1988. a
Ainslie, J. F.: Wake modelling and the prediction of turbulence properties, in: Proceedings of the Bwea Wind Energy Conference, british Wind Energy Association, 20–24 October 1986, Cambridge, 115–120, 1986. a
Bak, C., Zahle, F., Bitsche, R., et al.: Description of the DTU 10 MW reference wind turbine, technical report, DTU Wind Energy, Roskilde, Denmark, 2013. a, b, c
Barthelmie, R. J., Hansen, K. S., Frandsen, S. T., Rathmann, O., Schepers, J., Schlez, W., Phillips, J., Rados, K., Zervos, A., Politis, E., and Chaviaropoulos, P.: Modelling and Measuring Flow and Wind Turbine Wakes in Large Wind Farms Offshore, Wind Energy, 12, 431–444, https://doi.org/10.1002/we.348, 2009. a
Bauweraerts, P. and Meyers, J.: Bayesian based estimation of turbulent flow fields from lidar observations in a conventionally neutral atmospheric boundary layer, J. Phys. Conf. Ser., 1618, 032047, https://doi.org/10.1088/17426596/1618/3/032047, 2020. a
Bauweraerts, P. and Meyers, J.: Reconstruction of turbulent flow fields from lidar measurements using largeeddy simulation, J. Fluid Mech., 906, A17, https://doi.org/10.1017/jfm.2020.805, 2021. a
Bergami, L. and Gaunaa, M.: Analysis of aeroelastic loads and their contributions to fatigue damage, J. Phys. Conf. Ser., 555, 012007, https://doi.org/10.1088/17426596/555/1/012007, 2014. a
Bingöl, F., Mann, J., and Larsen, G. C.: Light detection and ranging measurements of wake dynamics Part I: Onedimensional Scanning, Wind Energy, 13, 51–61, https://doi.org/10.1002/we.352, 2010. a
Borraccino, A., Schlipf, D., Haizmann, F., and Wagner, R.: Wind field reconstruction from nacellemounted lidar shortrange measurements, Wind Energ. Sci., 2, 269–283, https://doi.org/10.5194/wes22692017, 2017. a, b, c
Bossanyi, E.: Unfreezing the turbulence: application to LiDARassisted wind turbine control, IET Renew. Power Gen., 7, 321–329, https://doi.org/10.1049/ietrpg.2012.0260, 2013. a
Bossanyi, E. A., Kumar, A., and HuguesSalas, O.: Wind turbine control applications of turbinemounted LIDAR, J. Phys. Conf. Ser., 555, 012011, https://doi.org/10.1088/17426596/555/1/012011, 2014. a
Chamorro, L. P., Guala, M., Arndt, R. E., and Sotiropoulos, F.: On the evolution of turbulent scales in the wake of a wind turbine model, J. Turbul., 13, 1–13, https://doi.org/10.1080/14685248.2012.697169, 2012. a, b
Churchfield, M. J., Moriarty, P. J., Hao, Y., Lackner, M. A., Barthelmie, R., Lundquist, J. K., and Oxley, G. S.: A comparison of the dynamic wake meandering model, largeeddy simulation, and field data at the egmond aan Zee offshore wind plant, in: 33rd Wind Energy Symposium, 5–9 January 2015, Kissimmee, Florida, 20 pp., 2015. a
Conti, D. and Dimitrov, N.: Constrained Gaussian turbulence field simulations (Version 0.1.0), Zenodo, https://doi.org/10.5281/zenodo.4896514, 2021. a
Conti, D., Dimitrov, N., and Peña, A.: Aeroelastic load validation in wake conditions using nacellemounted lidar measurements, Wind Energ. Sci., 5, 1129–1154, https://doi.org/10.5194/wes511292020, 2020. a, b, c, d
de Mare, M. T. and Mann, J.: On the SpaceTime Structure of Sheared Turbulence, Bound.Lay. Meteorol., 160, 453–474, https://doi.org/10.1007/s105460160143z, 2016. a
Dimitrov, N., Borraccino, A., Peña, A., Natarajan, A., and Mann, J.: Wind turbine load validation using lidarbased wind retrievals, Wind Energy, 22, 1512–1533, https://doi.org/10.1002/we.2385, 2019. a, b
Dimitrov, N. K. and Natarajan, A.: Application of simulated lidar scanning patterns to constrained Gaussian turbulence fields for load validation, Wind Energy, 20, 79–95, https://doi.org/10.1002/we.1992, 2017. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p
Dimitrov, N. K., Natarajan, A., and Mann, J.: Effects of normal and extreme turbulence spectral parameters on wind turbine loads, Renew. Energ., 101, 1180–1193, https://doi.org/10.1016/j.renene.2016.10.001, 2017. a, b, c
Dimitrov, N., Kelly, M. C., Vignaroli, A., and Berg, J.: From wind to loads: wind turbine sitespecific load estimation with surrogate models trained on highfidelity load databases, Wind Energ. Sci., 3, 767–790, https://doi.org/10.5194/wes37672018, 2018. a
Doubrawa, P., Barthelmie, R. J., Wang, H., and Churchfield, M. J.: A stochastic wind turbine wake model based on new metrics for wake characterization, Wind Energy, 20, 449–463, https://doi.org/10.1002/we.2015, 2017. a
Doubrawa, P., Debnath, M., Moriarty, P. J., Branlard, E., Herges, T. G., Maniaci, D. C., and Naughton, B.: Benchmarks for Model Validation based on LiDAR Wake Measurements, J. Phys. Conf. Ser., 1256, 012024, https://doi.org/10.1088/17426596/1256/1/012024, 2019. a, b
Fuertes, F. C., Markfort, C. D., and PorteacuteAgel, F.: Wind Turbine Wake Characterization with NacelleMounted Wind Lidars for Analytical Wake Model Validation, Remote Sens.Basel, 10, 668, https://doi.org/10.3390/rs10050668, 2018. a
Held, D. P. and Mann, J.: Detection of wakes in the inflow of turbines using nacelle lidars, Wind Energ. Sci., 4, 407–420, https://doi.org/10.5194/wes44072019, 2019a. a
Held, D. P. and Mann, J.: Lidar estimation of rotoreffective wind speed – an experimental comparison, Wind Energ. Sci., 4, 421–438, https://doi.org/10.5194/wes44212019, 2019b. a
Herges, T. G. and Keyantuo, P.: Robust Lidar Data Processing and Quality Control Methods Developed for the SWiFT Wake Steering Experiment, J. Phys. Conf. Ser., 1256, 012005, https://doi.org/10.1088/17426596/1256/1/012005, 2019. a
Hoffman, Y. and Ribak, E.: Constrained realizations of Gaussian fields – A Simple algorithm, Astrophys. J., 380, L5–L8, https://doi.org/10.1086/186160, 1991. a
Kaimal, J., Izumi, Y., Wyngaard, J., and Cote, R.: Spectral characteristics of surfacelayer turbulence, Q. J. Roy. Meteor. Soc., 98, 563, https://doi.org/10.1002/qj.49709841707, 1972. a
Keck, R.E., Veldkamp, D., Aagaard Madsen, H., and Larsen, G. C.: Implementation of a Mixing Length Turbulence Formulation Into the Dynamic Wake Meandering Model, J. Sol. Energy Eng., 134, 021012, https://doi.org/10.1115/1.4006038, 2012. a, b
Keck, R.E., de Mare, M. T., Churchfield, M. J., Lee, S., Larsen, G. C., and Aagaard Madsen, H.: On atmospheric stability in the dynamic wake meandering model, Wind Energy, 17, 1689–1710, https://doi.org/10.1002/we.1662, 2014. a, b, c, d, e
Keck, R. E., De Maré, M., Churchfield, M. J., Lee, S., Larsen, G., and Madsen, H. A.: Two improvements to the dynamic wake meandering model: Including the effects of atmospheric shear on wake turbulence and incorporating turbulence buildup in a row of wind turbines, Wind Energy, 18, 111–132, https://doi.org/10.1002/we.1686, 2015. a, b, c
Kretschmer, M., Schwede, F., Faerron Guzmán, R., Lott, S., and Cheng, P. W.: Influence of atmospheric stability on the load spectra of wind turbines at alpha ventus, J. Phys. Conf. Ser., 1037, 052009, https://doi.org/10.1088/17426596/1037/5/052009, 2018. a
Kretschmer, M., Pettas, V., and Cheng, P. W.: Effects of wind farm downregulation in the offshore wind farm Alpha ventus, in: ASME 2019 2nd International Offshore Wind Technical Conference, Iowtc 2019, 3–6 November 2019, St. Julian's, Malta, https://doi.org/10.1115/IOWTC20197554, 2019. a, b
Kristensen, L., Lenschow, D., Kirkegaard, P., and Courtney, M.: The Spectral Velocity Tensor for Homogeneous Boundary Layer Turbulence, Bound.Lay. Meteorol., 47, 149–193, https://doi.org/10.1007/BF00122327, 1989. a
Kumer, V. M., Reuder, J., and Eikill, R. O.: Characterization of turbulence in wind turbine wakes under different stability conditions from static Doppler LiDAR measurements, Remote Sens.Basel, 9, 242, https://doi.org/10.3390/rs9030223, 2017. a
Larsen, G., Ott, S., Liew, J., van der Laan, M., Simon, E., R.Thorsen, G., and Jacobs, P.: Yaw induced wake deflection – a fullscale validation study, J. Phys. Conf. Ser., 1618, 062047, https://doi.org/10.1088/17426596/1618/6/062047, 2020. a
Larsen, G. C., Madsen Aagaard, H., Bingöl, F., Mann, J., Ott, S., Sørensen, J., Okulov, V., Troldborg, N., Nielsen, N. M., Thomsen, K., Larsen, T. J., and Mikkelsen, R.: Dynamic wake meandering modeling, Risø National Laboratories, Roskilde, Denmark, 2007. a
Larsen, G. C., Madsen Aagaard, H., Thomsen, K., and Larsen, T. J.: Wake meandering: A pragmatic approach, Wind Energy, 11, 377–395, https://doi.org/10.1002/we.267, 2008. a, b, c
Larsen, T. J. and Hansen, A. M.: How 2 HAWC2, the user's manual, Risø National Laboratory, Roskilde, Denmark, 2007. a, b
Larsen, T. J., Aagaard Madsen, H., Larsen, G. C., and Hansen, K. S.: Validation of the dynamic wake meander model for loads and power production in the Egmond aan Zee wind farm, Wind Energy, 16, 605–624, https://doi.org/10.1002/we.1563, 2013. a, b, c, d, e
Lee, S., Churchfield, M., Moriarty, P., Jonkman, J., and Michalakes, J.: Atmospheric and wake turbulence impacts on wind turbine fatigue loadings, in: 50th Aiaa Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, AIAA 2012–0540, 9–12 January 2012, Nashville, Tennessee, https://doi.org/10.2514/6.2012540, 2012. a
Liew, J., Raimund Pirrung, G., and Meseguer Urbán, A.: Effect of varying fidelity turbine models on wake loss prediction, J. Phys. Conf. Ser., 1618, 062002, https://doi.org/10.1088/17426596/1618/6/062002, 2020. a
Lundquist, J. K., Churchfield, M. J., Lee, S., and Clifton, A.: Quantifying error of lidar and sodar Doppler beam swinging measurements of wind turbine wakes using computational fluid dynamics, Atmos. Meas. Tech., 8, 907–920, https://doi.org/10.5194/amt89072015, 2015. a
Lydia, M., Kumar, S. S., Selvakumar, A. I., and Prem Kumar, G. E.: A comprehensive review on wind turbine power curve modeling techniques, Renew. Sust. Energ. Rev., 30, 452–460, https://doi.org/10.1016/j.rser.2013.10.030, 2014. a
Machefaux, E., Larsen, G. C., Koblitz, T., Troldborg, N., Kelly, M. C., Chougule, A. S., Hansen, K. S., and Rodrigo, J. S.: An experimental and numerical study of the atmospheric stability impact on wind turbine wakes, Wind Energy, 19, 1785–1805, https://doi.org/10.1002/we.1950, 2016. a, b
Madsen, H. A., Larsen, G. C., and Thomsen, K.: Wake flow characteristics in low ambient turbulence conditions, in: Proceedings (cdrom), Copenhagen Offshore Wind Conference, Copenhagen, Denmark, 2005. a, b, c
Madsen, H. A., Larsen, G. C., Larsen, T. J., Troldborg, N., and Mikkelsen, R. F.: Calibration and Validation of the Dynamic Wake Meandering Model for Implementation in an Aeroelastic Code, J. Sol. Energy Eng., 132, 041014, https://doi.org/10.1115/1.4002555, 2010. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p
Mann, J.: The spatial structure of neutral atmospheric surfacelayer turbulence, J. Fluid Mech., 273, 141–168, 1994. a, b, c, d
Mann, J., Pena Diaz, A., Bingöl, F., Wagner, R., and Courtney, M.: Lidar Scanning of Momentum Flux in and above the Atmospheric Surface Layer, J. Atmos. Ocean. Tech., 27, 959–976, https://doi.org/10.1175/2010jtecha1389.1, 2010. a
Mann, J., Peña, A., Troldborg, N., and Andersen, S. J.: How does turbulence change approaching a rotor?, Wind Energ. Sci., 3, 293–300, https://doi.org/10.5194/wes32932018, 2018. a
Medley, J., Barker, W., Harris, M., Pitter, M., Slinger, C., Mikkelsen, T., and Sjöholm, M.: Evaluation of wind flow with a nacellemounted, continuous wave wind lidar, in: Proceedings of EWEA 2014, EWEA, Barcelona, Spain, 2014. a
Moens, M., Coudou, N., and Philippe, C.: A numerical study of correlations between wake meandering and loads within a wind farm, J. Phys. Conf. Ser., 1256, 012012, https://doi.org/10.1088/17426596/1256/1/012012, 2019. a
Muller, Y. A., Aubrun, S., and Masson, C.: Determination of realtime predictors of the wind turbine wake meandering, Exp. Fluids, 56, 1–11, https://doi.org/10.1007/s0034801519239, 2015. a
Nebenführ, B. and Davidson, L.: Prediction of windturbine fatigue loads in forest regions based on turbulent LES inflow fields, Wind Energy, 20, 1003–1015, https://doi.org/10.1002/we.2076, 2017. a
Nielsen, M., Larsen, G. C., Mann, J., Ott, S., Hansen, K. S., and Pedersen, B.: Wind Simulation for Extreme and Fatigue Loads, Risø National Laboratory, Roskilde, Denmark, 2003. a
Ning, X. and Wan, D.: LES study of wake meandering in different atmospheric stabilities and its effects on wind turbine aerodynamics, SustainabilityBasel, 11, 6939, https://doi.org/10.3390/su11246939, 2019. a, b
Pedersen, M. M., Larsen, T. J., Madsen, H. A., and Larsen, G. C.: More accurate aeroelastic windturbine load simulations using detailed inflow information, Wind Energ. Sci., 4, 303–323, https://doi.org/10.5194/wes43032019, 2019. a
Pettas, V., García, F. C., Kretschmer, M., Rinker, J. M., Clifton, A., and Cheng, P. W.: A numerical framework for constraining synthetic wind fields with lidar measurements for improved load simulations, AIAA Scitech 2020 Forum, Orlando, Florida, https://doi.org/10.2514/6.20200993, 2020. a, b, c, d, e, f
Peña, A., Hasager, C. B., Badger, M., Barthelmie, R. J., Bingöl, F., Cariou, J.P., Emeis, S., Frandsen, S. T., Harris, M., Karagali, I., Larsen, S. E., Mann, J., Mikkelsen, T., Pitter, M., Pryor, S., Sathe, A., Schlipf, D., Slinger, C., and Wagner, R.: Remote Sensing for Wind Energy, DTU Wind Energy, Roskilde, Denmark, 2015. a, b
Peña, A., Mann, J., and Dimitrov, N.: Turbulence characterization from a forwardlooking nacelle lidar, Wind Energ. Sci., 2, 133–152, https://doi.org/10.5194/wes21332017, 2017. a, b, c, d, e
Peña, A., Mann, J., and Rolighed Thorsen, G.: SpinnerLidar measurements for the CCAV52, DTU Wind Energy, Roskilde, Denmark, 2019. a
Raach, S., Schlipf, D., and Cheng, P. W.: Lidarbased wake tracking for closedloop wind farm control, Wind Energ. Sci., 2, 257–267, https://doi.org/10.5194/wes22572017, 2017. a
Reinwardt, I., Schilling, L., Dalhoff, P., Steudel, D., and Breuer, M.: Dynamic wake meandering model calibration using nacellemounted lidar systems, Wind Energ. Sci., 5, 775–792, https://doi.org/10.5194/wes57752020, 2020. a, b, c, d, e, f
Rommel, D. P., Di Maio, D., and Tinga, T.: Calculating wind turbine component loads for improved life prediction, Renew. Energ., 146, 223–241, https://doi.org/10.1016/j.renene.2019.06.131, 2020. a
Sathe, A. and Mann, J.: A review of turbulence measurements using groundbased wind lidars, Atmos. Meas. Tech., 6, 3147–3167, https://doi.org/10.5194/amt631472013, 2013. a
Sathe, A., Mann, J., Barlas, T. K., Bierbooms, W., and van Bussel, G.: Influence of atmospheric stability on wind turbine loads, Wind Energy, 16, 1013–1032, https://doi.org/10.1002/we.1528, 2013. a, b
Schlipf, D.: Lidarassisted control concepts for wind turbines, PhD thesis, Universitat Stuttgart, Stuttgart, Germany, 2016. a
Schlipf, D., Schlipf, D. J., and Kuehn, M.: Nonlinear model predictive control of wind turbines using LIDAR, Wind Energy, 16, 1107–1129, https://doi.org/10.1002/we.1533, 2013. a, b
Schlipf, D., Guo, F., and Raach, S.: Lidarbased Estimation of Turbulence Intensity for Controller Scheduling, J. Phys. Conf. Ser., 1618, 032053, https://doi.org/10.1088/17426596/1618/3/032053, 2020. a
Simley, E., Pao, L. Y., Kelley, N., Jonkman, B., and Frehlich, R.: LIDAR wind speed measurements of evolving wind fields, in: 50th Aiaa Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, AIAA, 9 January–12 January 2012, Nashville, Tennessee, 2012–0656, https://doi.org/10.2514/6.2012656, 2012. a
Simley, E., Y. Pao, L., Frehlich, R., Jonkman, B., and Kelley, N.: Analysis of Wind Speed Measurements using Continuous Wave LIDAR for Wind Turbine Control, 49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Orlando, Florida, 2013. a, b, c
Simley, E., Fürst, H., Haizmann, F., and Schlipf, D.: Optimizing lidars for wind turbine control applicationsResults from the IEA Wind Task 32 workshop, Remote Sens.Basel, 10, 863, https://doi.org/10.3390/rs10060863, 2018. a, b, c
Singh, A., Howard, K. B., and Guala, M.: On the homogenization of turbulent flow structures in the wake of a model wind turbine, Phys. Fluids, 26, 025103, https://doi.org/10.1063/1.4863983, 2014. a
Tibaldi, C., Henriksen, L. C., Hansen, M. H., and Bak, C.: Wind turbine fatigue damage evaluation based on a linear model and a spectral method, Wind Energy, 19, 1289–1306, https://doi.org/10.1002/we.1898, 2015. a
Trujillo, J.J., Bingöl, F., Larsen, G. C., Mann, J., and Kühn, M.: Light detection and ranging measurements of wake dynamics. Part II: twodimensional scanning, Wind Energy, 14, 61–75, https://doi.org/10.1002/we.402, 2011. a, b, c, d, e
Vaspe, A.: SWEUniStuttgart/ViConDAR: ViConDAR V1.0 (Version V1.0), Zenodo, https://doi.org/10.5281/zenodo.4889772, 2021. a
Wagner, R., Friis Pedersen, T., Courtney, M., Antoniou, I., Davoust, S., and Rivera, R.: Power curve measurement with a nacelle mounted lidar, Wind Energy, 17, 1441–1453, https://doi.org/10.1002/we.1643, 2014. a
Wagner, R., Courtney, M. S., Friis Pedersen, T., and Davoust, S.: Uncertainty of power curve measurement with a twobeam nacellemounted lidar, Wind Energy, 19, 1269–1287, https://doi.org/10.1002/we.1897, 2015. a
Zhan, L., Letizia, S., and Valerio Iungo, G.: LiDAR measurements for an onshore wind farm: Wake variability for different incoming wind speeds and atmospheric stability regimes, Wind Energy, 23, 501–527, https://doi.org/10.1002/we.2430, 2020. a
Zwick, D. and Muskulus, M.: The simulation error caused by input loading variability in offshore wind turbine structural analysis, Wind Energy, 18, 1421–1432, https://doi.org/10.1002/we.1767, 2015. a