Probabilistic estimation of the Dynamic Wake Meandering model parameters using SpinnerLidar-derived wake characteristics
- 1Department of Wind Energy, Technical University of Denmark, Frederiksborgvej 399, 4000 Roskilde, Denmark
- 2Wind Energy Technologies, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
Correspondence: Davide Conti (firstname.lastname@example.org)
We study the calibration of the Dynamic Wake Meandering (DWM) model using high-spatial- and high-temporal-resolution SpinnerLidar measurements of the wake field collected at the Scaled Wind Farm Technology (SWiFT) facility located in Lubbock, Texas, USA. We derive two-dimensional wake flow characteristics including wake deficit, wake turbulence, and wake meandering from the lidar observations under different atmospheric stability conditions, inflow wind speeds, and downstream distances up to five rotor diameters. We then apply Bayesian inference to obtain a probabilistic calibration of the DWM model, where the resulting joint distribution of parameters allows for both model implementation and uncertainty assessment. We validate the resulting fully resolved wake field predictions against the lidar measurements and discuss the most critical sources of uncertainty. The results indicate that the DWM model can accurately predict the mean wind velocity and turbulence fields in the far-wake region beyond four rotor diameters as long as properly calibrated parameters are used, and wake meandering time series are accurately replicated. We show that the current DWM model parameters in the IEC standard lead to conservative wake deficit predictions for ambient turbulence intensities above 12 % at the SWiFT site. Finally, we provide practical recommendations for reliable calibration procedures.
Wake effects are perceived as one of the largest sources of uncertainty in energy production and load estimates of onshore and offshore wind farms (Walker et al., 2016). Within an iterative design process and/or optimization study, wake effects on aeroelastic turbine responses are predicted using engineering wake models, e.g. the Dynamic Wake Meandering (DWM) (Madsen et al., 2010) and Frandsen (Frandsen, 2007) models, which can be used within simple and fast design tools (Braunbehrens and Segalini, 2019). Their main limitation is their reduced ability to fully resolve the turbulence structures of the wake field, which often leads to an inaccurate representation of the flow field and biased power and load predictions (Reinwardt et al., 2018). To minimize the modelling uncertainty, it is a common practice to calibrate engineering wake models using field measurements when available or using higher-fidelity simulations like computational fluid dynamics (CFD).
Wind lidars have become popular for studying wind turbine wakes due to their higher spatial resolution and ease of installation compared to traditional anemometers mounted on meteorological masts (Machefaux et al., 2016). The use of lidar measurements to calibrate low-order wake models has already been successfully adopted (Trabucchi et al., 2017; Reinwardt et al., 2020; Zhan et al., 2020a). Although high-quality lidar observations of the wake field are available (Käsler et al., 2010; Iungo et al., 2013; Aitken et al., 2014), the spatial and temporal resolution required to characterize wake deficit, wake turbulence, and meandering characteristics is rarely achieved. Such resolution is a key characteristic for the development and evaluation of dynamic wake models.
The Scaled Wind Farm Technology (SWiFT) experiment, conducted at Sandia National Laboratories between 2016 and 2017 (Herges et al., 2017, 2018; Herges and Keyantuo, 2019), provides a fairly complete and suitable dataset for the calibration and evaluation of wake models (Doubrawa et al., 2019, 2020; Conti et al., 2020a). The SWiFT dataset consists of concurrent measurements of inflow conditions from a heavily instrumented meteorological mast, high-spatial- and high-temporal-resolution measurements of a single wake flow field behind a turbine from a nacelle-mounted SpinnerLidar, and power and load measurements from a second turbine operating in the waked field. The detailed instrumentation of the site allows the investigation of the wake field variability under different atmospheric-stability conditions as well as the analysis of the wake-induced effects on the waked-turbine operation (i.e. power and load predictions).
Here, we analyse the SWiFT dataset aiming at calibrating and evaluating the DWM model. This model is recommended in the IEC 61400-1 standard (IEC, 2019) for the purpose of wind turbine and wind farm design certification, and it is widely used in load assessments under wake conditions (Larsen et al., 2013; Galinos et al., 2016; Reinwardt et al., 2018, 2020; Dimitrov, 2019). The DWM model simulates wind field time series and is divided into three parts: a wake deficit component, which simulates the velocity deficit; a wake-added turbulence component; and a wake meandering component, which is a stochastic meandering process. These three components are presumed to affect wind turbine loading conditions (Keck et al., 2012; Galinos et al., 2016; Larsen et al., 2013; Dimitrov, 2019). Although several studies have demonstrated the superior performance of the DWM model compared to other engineering wake models that only predict steady wake features (Thomsen et al., 2007; Larsen et al., 2013; Reinwardt et al., 2018), the accuracy of both the DWM-simulated wake flow fields and the resultant turbine power and load predictions is still to be assessed.
1.1 A review of the DWM model
The underlying hypothesis of the DWM model is to consider the wake as a passive tracer of the large incoming turbulence structures. The so-called split-in-scales assumption (Larsen et al., 2008) states that the large-scale turbulent eddies contained in the atmospheric boundary layer are the main drivers of the wake meandering, whereas the smaller turbulent eddies govern the wake deficit evolution downstream of the rotor. Further, wake deficits from upstream turbines are transported in the streamwise direction, assuming Taylor's hypothesis of frozen turbulence (Larsen et al., 2015). This set of assumptions allows for the decoupling of the wake deficit and wake-added turbulence formulations from the wake meandering process (Larsen et al., 2007). Therefore, the three components of the DWM model can be computed separately and successively superimposed on turbulence fields to generate wake time series, which can be used as inputs to aeroelastic simulations (Larsen et al., 2013; Keck et al., 2014a).
The wake deficit formulation of the DWM model is mainly based on the work of Ainslie (1987) and solves the axisymmetric Navier–Stokes (N–S) equations with an eddy viscosity term and a set of calibration parameters. Initially, the DWM model was calibrated with CFD simulations performed by Madsen et al. (2010). Keck et al. (2012) derived a two-dimensional model of the eddy viscosity term and updated the calibration parameters based on CFD simulations. Larsen et al. (2013) found that the calibration parameters of the former two studies were not suitable for predicting power and loads at the Egmond aan Zee offshore wind farm. To match the measured power, they introduced an artificial filtering function in the eddy viscosity term and re-calibrated the deficit model; however, this calibration was not based on the spatial description of the wake flow field but on power production data. The eddy viscosity model to predict velocity deficits in the current IEC standard (IEC, 2019) is inspired by the work of Larsen et al. (2013).
Keck et al. (2014a, 2015) proposed a correction factor to the eddy viscosity term, which includes the effects of atmospheric stability and shear on the turbulence mixing occurring in the wake, and re-calibrated the model parameters. Although these improvements were verified against large-eddy simulations (LESs), the influence of atmospheric stability on the wake deficit evolution was hardly observed during a lidar campaign (Machefaux et al., 2016; Larsen et al., 2015), in which it was argued that atmospheric stability affects to a large extent the meandering process. A load validation study using the DWM model with calibrated parameters from both Madsen et al. (2010) and Keck et al. (2012) as well as the IEC standard (IEC, 2019) was conducted by Reinwardt et al. (2018), who collected load measurements at the ECN wind turbine test site in Germany and at the Technical University of Denmark (DTU) test site in Høvsøre in Denmark. They found fatigue load biases within the range of 11 %–15 % for the tower bottom and 8 %–21 % for the blade-root flapwise bending moments. Reinwardt et al. (2020) derived a new set of calibration parameters based on full-field lidar observations of the wake field from a wind farm in the south-east of Hamburg, Germany. They demonstrated that improved wake deficit predictions can be obtained by calibrating the DWM model with nacelle-mounted lidars.
The fidelity of the simulated wake meandering dynamics also affects the accuracy of load predictions (Larsen et al., 2013; Conti et al., 2021). Modelling of the meandering process relies on a suitable stochastic turbulence field and definition of the large-scale turbulence structures. Larsen et al. (2008) and Trujillo et al. (2011) demonstrated that the large-scale eddies can be extracted from the incoming atmospheric turbulence field from local mast measurements. Alternatively, the wake meandering process can be simulated through synthetic wind fields generated using stochastic turbulence models (i.e. the turbulence model by Mann, 1994) and a definition of the large-scale eddies or by means of LESs. Machefaux et al. (2015) showed that inconsistencies between a Mann-based and LES-based meandering process can arise due to differences in the input turbulence fields. Larsen et al. (2008) and Trujillo et al. (2011) defined the large-scale eddies on the order of two rotor diameters (D) or larger as responsible for wake meandering, whereas other studies defined scales larger than 3–4 D as dominant (Espana et al., 2011; Muller et al., 2015; Yang and Sotiropoulos, 2019). Despite the severe impact of the wake meandering dynamics on load predictions, its uncertainty has not been assessed in load validation studies due to lack of data (Larsen et al., 2013; Churchfield et al., 2015; Reinwardt et al., 2018). However, aeroelastic simulations with constrained wake meandering dynamics can potentially decrease the uncertainty in load predictions under wake conditions (Conti et al., 2021).
Further, the added turbulence formulation in the DWM model accounts for additional mechanically generated turbulence caused by the wake shear and the breakdown of tip and root vortices. These contributions are modelled by a semi-empirical formulation that uses parameters, which were calibrated against CFD simulations (Madsen et al., 2010). To our knowledge, no further development has been made on this subject.
1.2 Problem statement
As described above, there is no consensus for the values of the DWM model parameters when studying load predictions at any given site. Also, and perhaps most importantly, we do not know the sources of uncertainty observed in previous studies that used the model (Larsen et al., 2013; Churchfield et al., 2015; Reinwardt et al., 2018), which need to be addressed to provide reliable load predictions. The common practice has been to derive optimized sets of model parameters based on limited synthetic or experimental data. This has led to an unknown confidence in the overall model prediction ability; incorrect calibration of the model parameters may impact significantly the model performance and lead to suboptimal wind turbine designs.
To address this issue, we estimate uncertainties in the calibration parameters of the DWM model by applying Bayesian inference (Box and Tiao, 1973), which consists of updating any related prior information on model parameters by incorporating new knowledge obtained from wake flow characteristics derived through lidar measurements. Further, the Bayesian calibration provides a systematic approach to include various types of uncertainty such as physical variability as well as measurement and modelling errors. This paper focuses on improving and validating the calibration of DWM model parameters using lidar-derived wake features and has a fourfold primary purpose:
Derive wake flow features such as the two-dimensional velocity deficit and wake-added turbulence profiles as well as time series of the wake meandering in both lateral and vertical directions from the SpinnerLidar measurements under different inflow wind speeds and atmospheric-stability conditions.
Calibrate the DWM-model-based wake deficit and wake-added turbulence predictions using the SpinnerLidar-derived wake flow features and the Bayesian inference framework.
Propagate modelling uncertainties in fully resolved wake flow fields for robust predictions that take into account the calibrated uncertainties.
Conduct a sensitivity analysis to determine the most significant sources of uncertainty in simulated wake fields that are typically inputs to aeroelastic simulations.
This study contributes to the ongoing discussion regarding the accuracy of power and load predictions of wind turbines operating under wake situations (Conti et al., 2020b, 2021) by quantifying uncertainties in wake simulations performed with the DWM model under a variety of inflow wind conditions. The outcomes of this study are useful for improving currently adopted wake simulation procedures for load analysis in the IEC standards as well as to provide practical recommendations for wake model calibration studies based on measurements from nacelle-mounted lidars.
The work is organized as follows. Section 2 describes the DWM model. The SWiFT layout and relative wind site conditions are described in Sect. 3. In Sect. 4, we present the wind field retrieval assumptions used to derive wake features from SpinnerLidar measurements. The Bayesian calibration of the DWM model is performed in Sect. 5. We carry out the validation of the wind turbine wake simulations and conduct a sensitivity analysis to investigate the most influential parameters in Sect. 6. Finally, the last two sections are dedicated to the discussions and conclusions.
The DWM model resolves three main wake features: the quasi-steady velocity deficit, the wake-added turbulence, and the wake meandering. Each model component is described separately in the following subsections.
2.1 Quasi-steady velocity deficit
The quasi-steady velocity deficit component describes the wake expansion and recovery caused partly by the recovery of the rotor pressure field and partly by turbulence diffusion moving farther downstream of the rotor (Larsen et al., 2013). The wake deficit is formulated in the meandering frame of reference (MFoR), which is a coordinate system with origin in the centre of symmetry of the deficit.
In the far-wake region, i.e. distances larger than two rotor diameters (Sanderse, 2015), the deficit evolution is assumed to be governed by turbulent mixing and is described by the thin shear layer approximation of the rotational symmetric N–S equations with the pressure term disregarded (Madsen et al., 2010). To account for the neglected pressure gradient effects, an initial wake deficit is analytically formulated based on the turbine's axial induction derived from blade element momentum (BEM) theory (Madsen et al., 2010). The turbulence closure of the N–S equations is obtained by means of an eddy viscosity term, and the momentum equation is solved numerically using a finite difference scheme with the artificial initial deficit as a boundary condition (Madsen et al., 2010). Here, we use the numerical scheme of the stand-alone DWM model (Liew et al., 2020; Larsen et al., 2020). We refer to the generalized definition of the non-dimensional eddy viscosity term by Keck et al. (2012), who considered two major drivers to the turbulence mixing: the ambient turbulence (TIamb) and turbulence induced by the wake shear layer:
where νT is the eddy viscosity, Uamb is the ambient wind speed at hub height, and R is the rotor radius. The first term on the right-hand side of Eq. (1) describes the contribution of the ambient turbulence and the second the self-generated turbulence by the wake shear layer. Madsen et al. (2010) proposed the wake radius , where is the downstream distance normalized by R, and the maximum velocity difference (Uamb−Umin), where Umin is the minimum wind speed in the wake, as the turbulent length and velocity scales, respectively, that govern turbulent mixing due to the wake shear layer.
Based on classical mixing length theory, Keck et al. (2012) defined the turbulence stresses to be proportional to the local velocity gradient , which provides a two-dimensional eddy viscosity formulation that is a function of the axial and radial coordinates, and r, respectively. The max operator is included to avoid underestimating the turbulent stresses at locations where the velocity gradient of the deficit approaches zero. Both terms in Eq. (1) include a filter function ( and ) and a model constant (k1 and k2). The filter functions are required to model the turbulence development behind the rotor and have values in the range of 0–1, depending on the downstream distance only (Keck et al., 2012). F1 accounts for the delay of the ambient turbulence entrainment into the wake and is assumed to “activate” ambient turbulence effects at downstream distances where the pressure has recovered (i.e. 2 D downstream, where the far wake begins; Sanderse, 2015).
F2 compensates for the initial non-equilibrium between the mean velocity field and the turbulent energy content created due to the rapid change in mean flow gradients close to the rotor. We refer to Eqs. (17) and (18) in Keck et al. (2015) for the mathematical formulation of F1 and F2; k1 and k2 are calibration parameters that govern the turbulence mixing and presumably do not change with wind turbine design and ambient conditions (Keck et al., 2012).
2.2 Wake turbulence
where TIm denotes the turbulence induced by the meandering of the wake deficit, and TIadd is the wake-added turbulence. TIm is commonly denoted as the apparent turbulence (Madsen et al., 2005) as the stochastic meandering of the wake deficit induces additional velocity fluctuations into time series taken at fixed locations in the wake. This term is considered the main source of added turbulence in the far wake (Madsen et al., 2010), while its spatial distribution can be computed by the convolution of the wake deficit in the MFoR and the probability distribution function (PDF) of the wake meandering in the lateral and vertical directions (Keck et al., 2014a). TIadd accounts for the shear- and mechanically generated turbulence due to blade tip and root trailing vortices. The inhomogeneity of the wake-added turbulence is modelled by scaling the local turbulence using the factor kmt (Madsen et al., 2010) as
where Udef, MFoR is the velocity deficit in the MFoR, and km1 and km2 are constants calibrated based on CFD results (Madsen et al., 2010). The wake-added turbulence derived from Eq. (3) is presumed to meander together with the wake deficit, thus being displaced by the large-scale eddies in the atmosphere.
2.3 Meandering model
Here, the meandering model is confined to a single wake scenario, whereas multiple wake dynamics are described in Machefaux (2015). The wake field is modelled by considering a cascade of consecutive wake deficits that are displaced by the large-scale lateral- and vertical-velocity fluctuations, i.e. the wake transport velocities (vc and wc), corresponding to the lateral (y) and the vertical axis (z), respectively. Adopting Taylor's hypothesis, the downstream advection of these deficits is assumed to be controlled by the mean wind speed of the ambient wind field. Larsen et al. (2008) estimated vc and wc by low-pass filtering of atmospheric turbulence fluctuations. They defined a filtering cut-off frequency , thus excluding contributions from smaller eddies to the meandering dynamics. This assumption was verified using full-scale lidar-based measurements collected behind an operating turbine (Bingöl et al., 2010). The wake displacements are computed as
where defines the time for an air particle to move from the rotor to the downstream distance in the wake region. An appropriate choice of the transport velocity of the wake advection lies between the ambient wind speed and the centre velocity of the wake deficit (Keck et al., 2014b; Machefaux et al., 2015). The contribution from the yaw misalignment, which can redirect wakes in the lateral direction, is accounted for by , where is the yaw offset at the specific time (Machefaux, 2015; Vollmer et al., 2016). The contribution of the rotor tilt is considered by (Machefaux, 2015).
The SWiFT facility is a research site located in Lubbock, Texas, operated by Sandia National Laboratories (Herges et al., 2017). The site includes three Vestas V27 wind turbines, two meteorological towers, and a SpinnerLidar (Peña et al., 2018) mounted on the nacelle of one of the turbines and looking backwards. The entire site is on a fibre optic data acquisition and control network that synchronizes recordings from masts, turbines, and the SpinnerLidar (Herges et al., 2017, 2018). The measurement campaign took place between 2016 and 2017 with the main objective of characterizing wake fields and investigating wake steering control strategies (Herges et al., 2017).
Figure 1 provides an overview of the test layout together with the notation used throughout the paper. In this study, we analyse data collected at the meteorological mast (METa1), the turbine (WTGa1), and the SpinnerLidar mounted on the nacelle of the WTGa1. The METa1 (hereafter referred to as the mast) is 60 m tall and instrumented with sonic anemometers at 10, 18, 32, 45, and 58 m, sampling at 100 Hz. Other instruments installed on the mast are reported in Herges et al. (2017). The mast is placed 2.5 D south of WTGa1 in compliance with the IEC standard guidelines (IEC, 2015, 2017). As southerly winds are prevalent at the site (see Fig. 1, right), this layout allows the retrieval of concurrent incoming wind conditions from METa1, wake measurements behind the WTGa1 performed by the SpinnerLidar, and power and load measurements on the waked WTGa2 installed 5 D downstream. The WTGa1 and WTGa2 are variable-speed and pitch-regulated turbines with a hub height of 32.1 m, D=27 m, a cut-in wind speed of 3 m/s, and a maximum power output of 192 kW reached at the rated wind speed of 12 m/s (Berg et al., 2014). The supervisory control and data acquisition (SCADA) is available for both turbines, providing records of the rotor speed, pitch and yaw angles, and power production, among other information, at 50 Hz.
The SpinnerLidar is a research Doppler wind lidar developed at DTU based on a continuous-wave (CW) laser system (Peña et al., 2018). Hereafter, SpinnerLidar and lidar denote the same system. The SpinnerLidar has been mounted either in the spinner or on top of the nacelle of a wind turbine (Angelou and Sjöholm, 2015; Peña et al., 2018). In this study, the SpinnerLidar was installed on the nacelle of the WTGa1 and scanned the rotor wake at a high temporal and spatial resolution so that wake features could be derived. For the SWiFT campaign, the SpinnerLidar scanned continuously in a rose-pattern every 2 s (see Fig. 2), and the system internally subdivided the rose into 984 sections. The accumulated Doppler-shifted spectra at each of the sections was also recorded (Herges et al., 2017).
Once a scan was completed, the SpinnerLidar refocused at a different range, and this process took about 2 s (Herges et al., 2018). For the SWiFT campaign, several scanning strategies were adopted as described below.
Strategy I. The SpinnerLidar scanned seven downstream distances: 1, 1.5, 2, 2.5, 3, 4, and 5 D. A full cycle (i.e. from 1 to 5 D) took 30–42 s. This dataset is suitable for investigating the wake deficit evolution and recovery behind the rotor; however, the frequency is too low to properly derive turbulence estimates or meandering dynamics.
Strategy II. The SpinnerLidar scanned at the fixed distance of 2.5 D, ensuring both high spatial and temporal resolution; ≈298 rosette scans were generated within a 10 min period. This dataset is suited for turbulence and meandering investigations.
Strategy III. The SpinnerLidar scanned at the fixed distance of 5 D behind the rotor, generating about 298 scans each 10 min. During this period, power and load measurements were recorded on WTGa2. This dataset is suitable for load validation analysis. Since it provides a description of the wake flow field, including velocity deficits, turbulence, and meandering at a distance that corresponds to typical spacings in wind farms, it is a valuable dataset for validating fully resolved wake flow predictions as long as induction effects are accounted for.
3.2 Site conditions
For extended periods of the campaign, WTGa1 operated under large yaw misalignment as wake steering strategies were being investigated (Herges et al., 2017). To consider periods where WTGa1 is nearly aligned with the mean inflow, we filtered out 10 min periods characterized by an average yaw offset larger than compared to the free-stream wind direction (Conti et al., 2020a). The yaw offset is here defined as the difference between the nacelle orientation and the wind direction measured at the mast. Further, we focus the analysis on periods for which the free-stream wind direction is within 90–270∘ (thus southern winds; see Fig. 1, right). This leads to about 850 available 10 min periods.
Figure 3 shows 10 min statistics of the hub-height turbulence intensity (TIamb), the power-law shear exponent (α), and the power production of WTGa1 as a function of the hub-height mean wind speed (Uamb) based on the mast inflow measurements; α is computed from the sonic measurements at 18 and 45 m. As shown, the site is characterized by a wide range of turbulence and shear conditions, which are a consequence of the varying atmospheric stability (Doubrawa et al., 2019; Conti et al., 2020a). Further, relatively low wind speeds are recorded (3–10 m/s); thus WTGa1 operates below rated power, as seen in Fig. 3c. Because of this range of operating conditions, high rotor thrust coefficients that induce strong wake deficits characterize this dataset.
3.2.1 Atmospheric stability
Here, we investigate the variability in the wake flow characteristics under varying stability and inflow wind speed conditions. We classify each 10 min sonically derived statistic into atmospheric-stability classes defined by ranges of the dimensionless stability parameter (), where L is the Obukhov length (Monin and Obukhov, 1954) computed from the sonic measurements as
where is the friction velocity, is the local kinematic momentum flux, k=0.4 is the von Kármán constant, g is the acceleration due to gravity, T is the mean surface-layer temperature, the vertical velocity component is denoted by w, and Θv is the virtual potential temperature (which we approximate by the sonic temperature). The prime denotes fluctuations around the mean value, and the overbar is a time average. We define three main atmospheric-stability classes based on ranges by Peña (2019): unstable (), near-neutral (), and stable () atmospheric conditions. We use the measurements at the 18 m sonic anemometer to derive the stability within each 10 min period. As shown in Conti et al. (2020a), the sonic measurements at 18 m provide the best fit to the polynomial form of Högström (1988), which describes the relation between the dimensionless wind shear ϕm and the dimensionless stability parameter in the surface layer (see middle panel of Fig. 3 in Conti et al., 2020a).
3.3 Data statistics
The statistics of the inflow wind parameters are presented separately in Tables 1, 2, and 3 according to the relative SpinnerLidar scanning strategy. Table 1 presents data collected during Strategy I. There is a fair number of 10 min periods to characterize the variability in the wake deficit with respect to atmospheric stability, inflow wind speeds, and downstream distances. The table shows increasing turbulence levels under unstable compared to stable cases, whereas relatively high vertical wind shears are found under stable conditions, as expected. The dataset is thus suitable for analysing the effects of atmospheric stability on the wake recovery. The dataset collected during Strategy II is reported in Table 2 and is used to characterize wake turbulence and meandering under different stability conditions. For Strategy III, represented in Table 3, the dataset is characterized by stable conditions mainly as the records correspond to night hours within three consecutive nights in July 2017.
As lidars only measure the line-of-sight (LOS) velocity (vlos), assumptions are needed to reconstruct the three-dimensional wind field , where u is the longitudinal, v the lateral, and w the vertical velocity component. If we neglect any probe volume averaging along the beam, vlos depends on the unit directional vector , which describes the scanning geometry through the elevation (ϕ) and azimuth (θ) angles and the wind field u,
Considering the small elevation angles and the typical low values of w, we assume w=0 (Doubrawa et al., 2019, 2020). This assumption may introduce an error of up to 3 % in the reconstructed horizontal wind speed at short distances (1–2 D) (Debnath et al., 2019). Following the approach of Doubrawa et al. (2020), we can combine the u and v velocity components into a total horizontal wind vector, U, and Eq. (6) becomes
where is the yaw offset, and the overbar indicates a smoothed signal as we apply a moving average operator with a 15 s window to the yaw misalignment to account for any temporal delay from the spatial distances among the mast, turbine's nacelle, and SpinnerLidar measurements (Conti et al., 2020a). With Eq. (7), we can reconstruct horizontal wind velocity measures at each individual scanned point within the rosette pattern. Further, we linearly interpolate the reconstructed wind speeds across the rosette pattern into a two-dimensional regular grid with a 2 m resolution, which is sufficient to characterize the spatial characteristics of the wind field in wakes (Fuertes et al., 2018; Conti et al., 2020a).
4.1 Lidar-estimated wake deficit
To perform comparisons with predicted velocity deficits from the DWM model, we aim at isolating the contribution of the wake deficit from that of the vertical wind shear in lidar measurements. As defined in Trujillo et al. (2011), the quasi-instantaneous wake deficit profile can be obtained by subtracting the mean vertical shear profile (Uamb(z)) from the quasi-instantaneous wake recording as
where is estimated from lidar measurements using Eq. (7), and Uamb(z) is the relative 10 min average inflow vertical wind speed profile measured at the mast. The deficit is then normalized with respect to the ambient wind speed profile. The vlos measurements and also the reconstructed U wind velocities are defined on a coordinate system that is attached either to the nacelle (nacelle frame of reference, NFoR), which rotates with the yawing of the turbine, or to the ground (fixed frame of reference, FFoR). To perform direct comparisons with the DWM model predictions, the lidar-estimated deficits obtained from Eq. (8) need to be computed in the MFoR. Here, this is performed by tracking the wake centre position through the method of Trujillo et al. (2011), where a bivariate Gaussian shape is fitted to the velocity deficit flow field, and the wake centre is the geometric centroid of the Gaussian function:
where μy and μz define the wake centre location; σwy and σwz are width parameters of the wake profile in the y and z directions, respectively; yi and zi denote the spatial locations of the lidar measurements; and A is a scaling parameter. Each scanned point of the quasi-instantaneous wake recording can be translated into the MFoR using the estimated μy and μz from Eq. (9) (Reinwardt et al., 2020). Therefore, we can compute the multiple wake recordings within a 10 min period in the MFoR and subsequently compute flow statistics such as the ensemble-average deficit profile as well as the spatial distribution of the wake turbulence in the MFoR. To ensure a high-quality fit, we reject scans where the estimated wake centre location is within ≈10 % of the lateral bounds of the scanning area and at more than 0.75 D from the hub height in the vertical direction (Doubrawa et al., 2020; Conti et al., 2020a).
Figure 4 illustrates ensemble-average measured deficit profiles in the MFoR at 2, 3, 4, and 5 D behind the rotor obtained from all 10 min periods characterized by an incoming wind speed of 7 m/s and under varying stability regimes during Strategy I (see Table 1 for reference). We can clearly observe the impact of the atmospheric stability and in particular of the associated turbulence levels on the wake recovery behind the rotor. A strong and well-defined symmetric wake deficit shape is seen under stable conditions (top row), whereas the deficits recover faster moving downstream as the atmosphere becomes more unstable (bottom row).
4.2 Lidar-estimated wake turbulence
Turbulence measures derived from lidar radial velocity measurements are “filtered” because of their relatively large probe volume (Peña et al., 2017), and so they are generally lower than those obtained from sonic observations. Nevertheless, if the Doppler spectrum of the vlos is available, we can potentially circumvent the averaging effects and estimate the unfiltered variance of vlos (Peña et al., 2017; Mann et al., 2010). Mann et al. (2010) assume that the ensemble-averaged Doppler spectrum over a time period 〈S(vlos)〉 is related to the probability distribution of the vlos at the focus distance and can be computed as
where φ(s) is the spatial averaging function of the lidar that depends on the position along the beam s, and p(vlos|s) denotes the PDF of vlos at the location s. If we assume that the PDF of vlos is independent of s, (i.e. there is no velocity gradient along the beam), then Eq. (10) reduces to 〈S(vlos)〉=p(vlos). As a result, the vlos statistics (i.e. mean and variance) can be computed from the first and second central moments of p(vlos) as
where and denote the mean and unfiltered variance of vlos, respectively. Nevertheless, velocity gradients along the lidar beam may appear when measuring at the wake edges, which can introduce errors in the estimated turbulence (Meyer Forsting et al., 2017). Following the procedure of Peña et al. (2019), we compute the ensemble-averaged normalized Doppler spectrum within 10 min periods by thresholding the noise-flattened spectra with a value of 1.2 and correcting them by subtracting the background spectrum. We accumulate the LOS Doppler spectra onto the regular grid of the scanned area and estimate and for each grid cell using Eq. (11). As discussed in Herges and Keyantuo (2019), invalid measurements occur due to the boresight and ground return as well as the return from the rotating rotor of WTGa2, if in operation. These invalid observations appear as a very high return signal in the Doppler spectrum in proximity to low wind speeds (i.e. at approximately 1 m/s) and are removed. The filtering effects due to the probe volume can be quantified by computing the ratio between filtered and unfiltered LOS variances across the rosette pattern; we find ratios in the range 0.8–0.9 at 2.5 D, which vary according to stability conditions (not shown).
Examples of 10 min ensemble-averaged Doppler spectra obtained at three fixed locations across the scanned area – a wake centre, a wake edge, and a wake-free position – are shown in Fig. 5 for an incoming wind speed of 7 m/s and ambient turbulence of 6 %. A narrow spectrum with a single-peak distribution centred at about 7 m/s for the wake-free location (green) is seen, whereas spectrum-broadening effects induced by small-scale generated turbulence are noticeable for the positions within the wake. The wake centre (red) shows a wider spectrum with a peak at a significantly lower wind speed than the incoming flow, whereas the wake edge (cyan) shows a double-peak distribution that may be partially due to the inhomogeneity of the wind field along the beam (Herges and Keyantuo, 2019) and also due to the meandering occurring within the analysed 10 min period.
To characterize the spatial distribution of the wake turbulence within the scanned area, we derive estimates directly by applying the variance operator to Eq. (7):
where is the variance of the horizontal wind speed, and as shown, covariance terms are neglected. As the LOS is almost never aligned with the u velocity component across the rosette, except at the centre of the pattern, can be “contaminated” by the variances and covariances of the other velocity components (Peña et al., 2017). Therefore, the relation in Eq. (12) can lead to inaccurate estimations of the longitudinal-velocity variances. Peña et al. (2019) estimated the contamination of different components on the LOS variances for the SpinnerLidar and showed that the ratio of the unfiltered LOS velocity variance to the variance of the longitudinal velocity component is generally lower than 1 across the scanned area, except at the centre, where the ratio is 1, and within an area above the centre, where it can be higher than unity. Although the adopted retrieval assumption in Eq. (12) introduces uncertainties in the turbulence measures, we can account for the expected errors in the Bayesian inference framework.
Figure 6 illustrates the spatial distribution of the unfiltered computed in the MFoR, normalized with the u-velocity variance of the ambient wind field measured at the 32 m sonic anemometer (). Under stable conditions and for a downstream distance of 2.5 D, we can observe an enhancement in turbulence levels in proximity to the rotor tips, especially in the upper part of the rotor (see Fig. 6a). The observed added turbulence is caused by the breakdown of the rotor tip vortices. These features are no longer noticed as the atmosphere becomes more unstable where a more uniform and less prominent distribution of the turbulence is found (see Fig. 6b and c).
The calibration of the wake deficit and wake-added turbulence components are conducted in the MFoR using a Bayesian inference framework. We describe the Bayesian model in Sect. 5.1 and provide calibration results for the wake deficit in Sect. 5.2 and for the wake-added turbulence in Sect. 5.3. We investigate wake meandering dynamics separately in Sect. 5.4.
5.1 Bayesian inference formulation
The basis of the Bayesian inference is to estimate the probability distribution of the model parameters based on available observations. Let be a set of model parameters to be estimated using lidar-derived wake features (i.e. wake deficit and wake-added turbulence profiles in the MFoR) denoted by , where n is the number of available observations. We consider that the experimental data and the model predictions satisfy the prediction error equation:
where denotes the DWM model predictions obtained from a particular set of model parameters (θm) and a set of observable variables (Xm). Here includes the inflow wind conditions measured at the mast (); the rotor thrust coefficient of the turbine (CT), which is derived from the BEM model (Madsen et al., 2010); and the spatial locations of the scanning pattern (). denotes a random prediction error composed of two terms: the measurement error ϵy and the model prediction error ϵm. The former is described by a zero mean normal distribution with standard deviation , which is determined from field observations. The latter is assumed to have zero mean, which implies unbiased model predictions, and a standard deviation to be determined by the Bayesian estimation along with the model parameters. To facilitate statistical inference, we assume Xm as deterministic inputs (i.e. free of uncertainty) and that the model error ϵm is independent of the set of input variables Xm and described by a normal distribution. This implies that the model predictions are normally distributed for a given Xm, which is a reasonable choice for wake deficit profiles. The Bayesian approach for model calibration deals with updating the combined parameter set , given a set of observations (yd,Xm) by applying the Bayes theorem:
where is the updated posterior distribution of the model parameters, is the prior distribution that is typically assigned based on subjective or previous information, denotes the likelihood of observing the data yd from a model with corresponding θm parameters, and f(yd) is the prior predictive distribution that is defined as the marginal distribution . By using the prediction error in Eq. (13) and assuming that the error terms are jointly normal with a zero mean vector and covariance matrix and , the measured quantities follow the normal distribution , where the covariance matrix takes the form . As a result, the likelihood function of observing the data follows the multi-variable normal distribution defined as
where the denotes the determinant. The analytical and differentiable solution of the posterior distribution of the parameters in the N–S equations with the eddy viscosity term of Eq. (1) is not readily available. Therefore, we employ a numerical sampling method to approximately evaluate the posterior distribution and its first and second moments. Here, the adaptive no-U-turn Markov chain Monte Carlo (MCMC) sampler is employed to generate samples from the posterior distribution (Hoffman and Gelman, 2014; Salvatier et al., 2016).
The outcome of the calibration is a joint probability distribution of the inferred model parameters. From this joint PDF, we can estimate the posterior PDF of any wake feature simulated by the DWM model, i.e. the wake deficit and wake-added turbulence profiles in the MFoR or the fully resolved wakes in the FFoR, among others, which we denote by q:
The posterior distribution of the wake feature q in Eq. (16) can be solved numerically using sampling methods (e.g. Monte Carlo simulations), so its first and second moment can be estimated.
5.2 Wake deficit parameter estimation
We use lidar-derived wake deficit profiles in the MFoR collected during Strategy I and Strategy II and employ the Bayesian model to infer uncertainty in the k1 and k2 parameters of the eddy viscosity term in Eq. (1). These parameters were found to be the most sensitive to the resulting wake deficit predictions (Keck et al., 2012). The prediction model of Eq. (13) is constructed as follows. The experimental data (yd) comprise two-dimensional ensemble-average lidar-estimated deficit profiles binned according to downstream distances (3, 4, and 5 D), atmospheric stability (i.e. stable, near-neutral, and unstable), and wind speed bins of 1 m/s in the range of 3–9 m/s. Note that we discard measurements in the near wake (at 2 D) for improving the quality of the fitting in the far-wake region (see Fig. 8).
The set of observable variables comprises , where the inflow parameters () are provided in Tables 1 and 2; the rotor thrust coefficient CT is derived from the BEM model implemented in the aeroelastic code HAWC2 (Larsen and Hansen, 2007) and based on the aerodynamics and airfoil inputs of the SWiFT turbine (Doubrawa et al., 2020) (CT=0.84, that is nearly constant for wind speeds below 9 m/s); and x,y, and z refer to the spatial coordinates of the deficits resolved in the MFoR. The uncertainties in measured deficit profiles are computed as , where σ(r) is the standard deviation of all 10 min deficits within the analysed case at the radial position r, and n is the number of 10 min periods (also referred to as samples in Table 1).
We select uniform prior distributions on model parameters and on intervals that consider physical constraints, ensuring convergence of results and covering previous calibrations reported in the literature. Thus, we employ an MCMC algorithm to sample from the posterior PDFs of the calibration parameters using the Bayesian framework. The inferred joint and marginal posterior distributions of the model parameters (k1 and k2) are shown in Fig. 7, together with point values from earlier studies (Madsen et al., 2010; Keck et al., 2012; Larsen et al., 2013; IEC, 2019; Reinwardt et al., 2020). As shown, the lidar-based wake deficits are informative, and we obtain well-defined posteriors that follow a normal distribution with and . The negative correlation between k1 and k2 seen in Fig. 7 indicates the interdependence of the physically induced effects as both parameters contribute to turbulence diffusion. It is found that the posterior means of the informed parameters k1 and k2 differ from those recommended in the IEC standard (IEC, 2019). Generally, low k1 and k2 values attenuate the degree of turbulence mixing in the wake, which lead to strong velocity deficits persisting at distances farther downstream. We provide the statistical properties (mean, variance, and coefficient of variation) of the inferred parameters in Table 4 and correlation measures in Table 5.
5.2.1 Wake deficit predictions
We propagate the uncertainties in k1 and k2 to predict wake deficits in the MFoR and compare them with the ensemble-average lidar-derived profiles in Fig. 8. First, we observe that the lidar-estimated deficits exhibit a faster wake recovery as the atmosphere becomes more unstable compared to stable regimes. This effect is mainly caused by the enhanced turbulence mixing occurring under unstable conditions as they are characterized by ambient turbulence levels 2–3 times higher than those of the stable cases (see Table 1). The lidar-observed maximum deficit varies between 30 % and 60 % within the first five rotor diameters, depending on the inflow turbulence conditions. Similar behaviours were reported in recent lidar measurement campaigns (Iungo and Porté-Agel, 2014; Machefaux et al., 2016; Fuertes et al., 2018; Zhan et al., 2020b).
The lidar-estimated deficits are approximately Gaussian under stable to near-neutral conditions, whereas the Gaussian shape is lost under more unstable conditions. This may result from errors in the wake tracking procedure due to the larger meandering amplitudes and also due to the presence of large-scale turbulence structures in the inflow (Conti et al., 2020a). The rotor thrust is another factor governing the variability in the wake recovery (Zhan et al., 2020b); however, its influence is secondary for the dataset analysed here due to the relatively low incoming wind speeds and relatively constant thrust coefficients (Conti et al., 2020a).
The DWM-model-predicted deficit profiles with parameters specified by their posterior distributions are in good agreement with the lidar observations for distances beyond 4 D (see Fig. 8). For these distances, the turbulence mixing effects dictated by the ambient and self-generated wake turbulence on the deficit recovery are fairly well captured by the inferred parameters. The nominal model predictions generally fit the observations, whereas an overlap between the measurements and the region of modelling uncertainty is found. The largest deviations between predicted and measured deficits are found at shorter distances (2–3 D) and are mainly due to the model inadequacy to simultaneously fit all the experimental measurements and experimental uncertainties. The assumptions adopted to describe the near-wake region also introduce uncertainty in the deficit predictions at short distances (Keck et al., 2015; Machefaux et al., 2016).
It can be observed that the uncertainties in k1 and k2 parameters primarily influence the depth of the wake (i.e. the maximum deficit), while the sensitivity to these parameters decreases significantly with the outer radial distance. It is also noticed that the uncertainty in the deficit predictions increases for high ambient turbulence and far downstream distances. This is because k1 is proportional to TIamb in Eq. (1), and the sensitivity of the model parameters increases as the wake recovers. To provide a measure of the uncertainty in wake deficit predictions, we compute the coefficient of variation COV = , where σ is the standard deviation, and μ is the mean value of the maximum deficit, obtained by propagating the PDFs of k1 and k2, as in Eq. (16). We find COV=3 % under stable conditions (Tamb=0.07), which increases to 6 % under unstable conditions (Tamb=0.14), for an incoming wind speed of 7 m/s. This result confirms that the uncertainty in wake deficit predictions increases for higher turbulence, but it also shows that uncertainties in k1 and k2 parameters do not lead to uncertainty of the same magnitude in deficit predictions resolved in the MFoR (for reference, % and %, as reported in Table 4).
We provide comparisons between measured and predicted wake deficit profiles using the calibration from this work as well as those reported in early studies in Fig. 9. For this particular analysis, we analyse predictions at 5 D behind the rotor for an inflow wind speed of 7 m/s under stable, near-neutral, and unstable regimes. The main discrepancy among the models is the relative sensitivity of the wake recovery to the ambient turbulence. This is primarily governed by k1; however, in the eddy viscosity model of Larsen et al. (2013), IEC (2019), and Reinwardt et al. (2020), it also depends on a nonlinear coupling function Famb(TIamb) that attenuates the wake recovery for turbulence above ≈12 % (see Fig. 6 in Larsen et al., 2013). This function was introduced to fit the power productions at the Egmond aan Zee offshore wind farm, and it is not based on observations of the wake field (Larsen et al., 2013). The wake recovery predicted with the models of Larsen et al. (2013) and IEC (2019) is practically insensitive to ambient turbulence rising from 7 % to 16 % at downstream distances up to 5 D. Similar outcomes are reported in Fig. 13 in Reinwardt et al. (2020), who showed that the model of Larsen et al. (2013) provided conservative deficits for ambient turbulence up to 16 %.
By considering k1 and k2 as universal constants (Keck et al., 2012), the quality of the dataset utilized for the model calibration is essential to ensure reliable parameter estimation. As previous calibrations were carried out on larger rotors than those of the SWiFT turbines and utilized either power production data (Larsen et al., 2013; IEC, 2019) or limited CFD simulations (Madsen et al., 2010; Keck et al., 2015) and one-dimensional scans of the wake by a nacelle lidar (Reinwardt et al., 2020), these aspects may explain the observed deviations in Fig. 9.
5.3 Improved wake-added turbulence formulation
The wake-added turbulence model (Eq. 3) assumes that turbulent structures (i.e. tip and root vortices) are unaffected by atmospheric turbulence. The rotor-induced vortices are rapidly disrupted under high-turbulence conditions, causing the breakdown within the first 2 D (Madsen et al., 2005). However, vortices can persist and extend at farther distances under low to moderate turbulence combined with stable stratification conditions (Ivanell et al., 2009; Subramanian et al., 2018; Conti et al., 2020a). Thus, the wake-added turbulence profiles can exhibit a more pronounced double-peak feature in the proximity of the rotor tips or a more uniform distribution depending on the atmospheric turbulence conditions (this effect is also seen in Fig. 6). Equation (3) also assumes radially symmetric wake-added turbulence profiles. However, the inflow vertical wind shear re-distributes the wake turbulence. Enhanced turbulence levels are actually observed in the proximity of the upper tip of the rotor blade (Vermeer et al., 2003; Chamorro and Porte-Agel, 2009; Conti et al., 2020a). Figure 6a shows this effect.
Due to these two assumptions, we propose an improved semi-empirical formulation of the wake-added turbulence scaling factor , which produces wake profiles in better agreement with the lidar observations. This is achieved by relating both the depth and the velocity gradient terms to the ambient turbulence and by including the effect of the inflow vertical wind shear on the vertical-velocity deficit gradient as
where (U(z)Udef(y,z)), and , and are parameters to be determined using Bayesian inference. Figure 10 illustrates the two-dimensional profiles of the depth and velocity deficit gradient terms of Eqs. (3) and (17). As illustrated, by including the term , we obtain a turbulence field that mimics qualitatively well the observed enhanced turbulence within the upper wake region.
5.3.1 Estimation of wake-added turbulence parameters
The calibration parameters in Eq. (17) are inferred based on the lidar-estimated wake-added turbulence profiles in the MFoR collected during Strategy II. For this particular dataset, the SpinnerLidar scans at a fixed distance of 2.5 D, ensuring about 298 scans for every 10 min period. The inflow characteristics are reported in Table 2. As the Doppler LOS velocity spectrum is available, we derive unfiltered LOS variances as described in Sect. 4.2 and subsequently the turbulence intensity as the ratio of the standard deviation to the mean of the horizontal wind speed. From Eq. (2), we can isolate TIadd (wake-added turbulence term) by firstly resolving the wake recordings in the MFoR, which eliminates the contribution of TIm, and then by subtracting TIamb, which is derived from the 32 m sonic observations in front of the rotor.
As a first step, we fit the “original” analytical formulation of kmt from Eq. (3) to each individual lidar-estimated wake-added turbulence profile to estimate the values of km1 and km2 using a simple least-squares optimization algorithm. Figure 11a and d show the relation between the estimated optimal parameters (in markers) and the ambient turbulence. It is shown that km1 increases, and km2 decreases almost linearly for increasing turbulence intensity. This indicates the strong effect of the deficit gradient term (proportional to km2 in Eq. 3) under low-turbulence conditions, which both amplifies the double-peak feature of the wake turbulence profile at the rotor tips and mimics the wake vortices. As the ambient turbulence increases, km1 and lower km2 become larger, which indicates a more uniform distribution of the wake turbulence. These effects are observed in Fig. 11b and c (see markers), which show the spanwise distribution of the lidar-derived wake-added turbulence for two 10 min periods with relatively low and high values of ambient turbulence intensity, 7 % and 12 %, respectively. We also compare the measured and predicted vertical distribution of TIadd in Fig. 11e and f, which shows that slightly improved predictions can be obtained using Eq. (17), i.e. considering the effects of the atmospheric shear on wake turbulence.
To infer the posterior PDFs of , and , we assign prior distributions based on the linear dependencies observed in Fig.11a and d, namely , , , and as well as and previously estimated in Sect. 5.2. The uncertainties in lidar-derived wake-added turbulence profiles account for errors introduced by the flow modelling assumptions of Eq. (12), which neglects cross-contamination effects on the LOS variance. We define these errors as zero-mean normally distributed with standard deviation , which leads to a coefficient of variation of 10 % (Peña et al., 2019).
The resulting posterior PDFs of the parameters follow a normal distribution (not shown) with statistical properties tabulated in Table 4. We provide the mean values of the inferred , and values in the form of a linear regression model that relates the depth and deficit gradient terms of Eq. (17) to the ambient turbulence in Fig. 11a and d. The shaded area represents the 95 % confidence interval, which is obtained by propagating the PDFs of the turbulence-related and wake deficit parameters. As shown, the predictions are characterized by a relatively high degree of uncertainty. This is because, for example, there are few available observations to characterize turbulence; the measurement uncertainties are relatively high; and the modelling simplifications, such as the semi-empirical , introduce further uncertainties. We provide the correlations among all the inferred parameters obtained from the joint PDF in Table 5. As expected, the wake deficit parameters k1 and k2 are negatively correlated to and as and are proportional to the wake deficit term in Eq. (17), which is the main driver to the intensity of the wake-added turbulence.
5.4 Wake meandering
Here, we investigate the relationship between the inflow turbulence fluctuations and the lidar-tracked wake positions to characterize the large-scale eddies responsible for the meandering. The analysis is carried out by comparing the spectra of the lidar-tracked meandering time series, which are derived by means of the tracking algorithm in Eq. (9), with that simulated from the meandering model in Eq. (4), where vc and wc are obtained from the sonic observations at 32 m. Yaw misalignment from SCADA is also included. We conduct the analysis using the dataset collected at 2.5 D behind the rotor during Strategy II and classify all the available 10 min periods according to stability to derive ensemble-average spectra from multiple observations with similar inflow conditions. Note that the measured inflow wind speeds are below rated, and turbulence levels from all 10 min periods within each stability class are similar to those reported in Table 2.
Results of the spectral analysis for both lateral and vertical meandering are shown in Fig. 12. Here, the ensemble-average spectra from the SpinnerLidar observations are compared to those from the meandering model without low-pass filtering the incoming turbulence fluctuations (denoted as DWMwf). The spectra are normalized with their relative variances and plotted as a function of the commonly used Strouhal number , where f denotes frequency, and U∞ is the aggregated wind speed from the ensemble-average statistics. As shown, the slope of the lidar-based spectra (red lines) matches that of DWMwf (blue lines) up to St = 0.3–0.5, which corresponds to three and two rotor diameters, respectively. For St> 0.5, the energy content of the lidar-estimated spectra remarkably decreases compared to that of DWMwf. These observations indicate that large-scale turbulent structures (>2 D) are dominant in the wake meandering (Trujillo et al., 2011; Heisel et al., 2018). When compared to the stable case, the spectra under unstable conditions show a higher energy content at large turbulence scales or equivalently low Strouhal number.
Figure 12 shows that a stochastic description of the large-scale eddy size might be appropriate. Thus, we describe the large-scale eddies responsible for the wake meandering by introducing the stochastic variable Dm, which is normally distributed with a mean equal to 2.5 D (it corresponds to the wake diameter at 5 D behind the rotor) and a standard deviation of 0.3 D based on the observations in Fig. 12. The resulting 95 % confidence intervals are shown as shaded areas in Fig. 12. The uncertainty in Dm is found negligible when computing wake meandering time series; this is shown in Appendix A.
The validation of the DWM model is performed by resolving wake fields in the FFoR; thus the simulated wakes include the combined effects of the velocity deficit, added turbulence, and wake meandering dynamics in both lateral and vertical directions. This analysis is carried out using data from Strategy III, i.e. at a fixed distance of 5 D behind the rotor, ensuring a sufficient number of scans to derive unfiltered turbulence estimates as well as wake meandering time series within a 10 min period. The analysed dataset is primarily characterized by stable conditions, as seen in Table 3, i.e. low turbulence intensities (6 %–8 %) and strong vertical shears (α = 0.18–0.38).
6.1 Correction for rotor induction effects
The SpinnerLidar measurements collected during Strategy III were taken in the induction zone of the WTGa2 (Herges et al., 2018). For induction correction, we employ the two-dimensional induction model of Troldborg and Meyer Forsting (2017), which accounts for both longitudinal and radial variation in the induced wind velocity:
where a0 is the induction factor at the rotor centre area, ; γa=1.1, is the distance in front of the rotor normalized by the rotor radius; denotes the radial distance from the rotor centre axis; and , , λa=0.587, and ηa=1.32 (Troldborg and Meyer Forsting, 2017; Dimitrov, 2019). Note that only lidar measurements taken across the rotor area of WTGa2 are corrected by the induction factor in Eq. (18). The estimated induction factors indicate that the wind speed can be reduced at hub height by up to 12 % upstream of the WTGa2 and below rated power (see also Fig. 14).
6.2 Uncertainty propagation of simulated wake fields
We derive the two-dimensional spatial distribution of the mean wind speed in the wake region (UFFoR) by the convolution between the wake deficit in the MFoR (Udef, MFoR) and the PDF of the meandering path (fm) (Keck et al., 2015):
where ym and zm denote the spatial coordinates of the wake meandering time series, and ym, ϵ and zm, ϵ are measures of their relative uncertainties. We introduce these errors to account for incorrect wake tracking positions that can arise due to the adopted wake tracking algorithm; ym, ϵ and zm, ϵ are assumed to be uncorrelated and to follow a normal distribution with zero mean and standard deviation such that the 95 % percentile corresponds to approximately 4 m, which is twice the resolution adopted to interpolate SpinnerLidar measurements onto the regular grid (see Sect. 4). Note that the atmospheric shear profile Uamb(z) is superposed after the wake deficit calculation (Madsen et al., 2010). Similarly, the two-dimensional spatial distribution of the u-velocity variance () can be computed as
where , UFFoR is derived from Eq. (19), and is the variance of the ambient u velocity component. Alternatively, UFFoR and can be equivalently derived by superposing the wake deficit and the wake-added turbulence factor on stochastic turbulence fields with constrained meandering path and subsequently by computing the first- and second-order statistics of the synthetic wind fields. However, the analytical forms in Eqs. (19) and (20) can be easily used to propagate the posterior PDFs of the calibration parameters to predict profiles of UFFoR and with relative uncertainties. The inflow parameters (Uamb, α, and ) in Eqs. (19) and (20) are required inputs for the DWM model and are estimated from mast measurements.
We compute UFFoR and from Eqs. (19) and (20) by constraining the meandering path (fm) using either the lidar-tracked meandering time series, which we denote as DWM*, or by using the meandering model in Eq. (4) with low-pass-filtered vc- and wc-velocity fluctuations obtained from the 32 m sonic anemometer and the yaw misalignment of WTGa1 obtained from SCADA, which we denote as DWM. In the latter model, the low-pass-filtered frequency is defined as a function of the stochastic variable associated with the size of the large-scale eddies (Dm), as discussed in Sect. 5.4. Figure 13 shows the comparison between observed and predicted (with both DWM* and DWM models) profiles of the mean wind speed and its standard deviation obtained from two different 10 min periods.
A good agreement between measurements and predictions is found. The vertical profile of the mean wind speed exhibits a single-peak shape resulting from the combined effects of the inflow vertical shear (modelled by a power law) and the wake-induced Gaussian-like deficit shape. The wake turbulence in the lateral direction exhibits a double-peak shape with larger values near the locations associated with strong velocity gradients that are further enhanced by the wake meandering. Enhanced turbulence levels in proximity of the upper wake region are observed from lidar measurements. The deviations between DWM* and model predictions are more pronounced for the simulated turbulence than for the wind speed fields and are exclusively due to differences in the meandering representations.
The wake simulation uncertainties shown in Fig. 13 are determined by propagating uncertainties in model parameters (Table 4) and by accounting for relative correlations (Table 5) using Monte Carlo simulations. The uncertainties in lidar-measured wind speeds account for volume averaging effects and for errors introduced by the retrieval assumptions (e.g. wsin (ϕ)=0 m/s) (Debnath et al., 2019). The uncertainties in lidar-derived turbulence (here defined as the standard deviation of the wind speed) account for errors introduced by neglecting cross-contamination effects in Eq. (12) (Peña et al., 2019).
To evaluate the performance of the DWM model, we calculate two flow metrics that are relevant in aeroelastic simulations and compare them to relative measured quantities: the rotor-effective wind speed (Ueff), defined as the weighted sum of wind speeds across the rotor area, and the effective wake turbulence (σu, eff), which is derived as the weighted sum of turbulence estimates across the rotor. Figure 14 shows the one-to-one comparison between measured and DWM*- and DWM-model-predicted flow metrics for all 10 min periods. We find a slope that deviates <1 % from unity, with R2=0.95 for Ueff and a bias of ≈ 4 % and R2≈0.93 for σu, eff. The observed scatter is explained by the large measurement uncertainties, by the uncertainties in inflow wind parameters, and by those of the model predictions. The latter is estimated by propagating uncertainties in model parameters provided in Table 4, which result in a COV of 1 % for Ueff, and of 3 % for σu, eff. These findings indicate that the inferred uncertainties in model parameters do not propagate into an error of the same magnitude on the fully resolved wakes that are eventually input to aeroelastic simulations.
Further, a sensitivity analysis indicates that the variations in Ueff and σu, eff are mostly explained by the uncertainties in k1 and k2 and the wake tracking position ; this is shown in Appendix B. The contribution of the wake-added turbulence from Eq. (17) to the total wake turbulence σu, eff in the FFoR is marginal and accounts for approximately 2 %–7 %. The turbulence induced by the meandering of the wake deficit is thus the major source of added turbulence (Madsen et al., 2010).
The SpinnerLidar measurements show that the wake recovers faster under unstable compared to stable conditions primarily due to the high turbulence levels of the former (Iungo et al., 2013; Zhan et al., 2020b). The wake recovery can be predicted accurately by the DWM model when using appropriate calibration parameters. Further, accurate reproduction of the wake meandering dynamics in both lateral and vertical directions is key for accurate wake simulations.
Although the SWiFT campaign provides a comprehensive dataset, including a wide range of inflow wind conditions and detailed two-dimensional lidar measurements of the wake field, the calibration parameters obtained from the 192 kW turbine with 32 m diameter should be further evaluated for multi-megawatt turbines with larger rotors. While this study cannot explicitly demonstrate the transferability of the obtained results for modern-size turbines, it highlights the need for datasets that include observations of the wake deficit profiles under varying stability conditions, inflow wind speeds, and downstream distances for the reliable and robust calibration of engineering wake models. Here, we find that the velocity deficit's recovery rate for increasing turbulence is the main difference among calibrations reported in the literature (see Fig. 9), which were based on either power production data (Larsen et al., 2013; IEC, 2019) or limited CFD simulations (Madsen et al., 2010; Keck et al., 2015) and one-dimensional observations of the wake field from a nacelle lidar (Reinwardt et al., 2020).
We also demonstrate that uncertainties in calibration parameters (e.g. COV % and COV %) do not lead to uncertainty of the same magnitude in the wind speed and turbulence fields (e.g. COV % and COV %) that are eventually used as input for the aeroelastic simulations. Note that the uncertainty in the simulated wake fields increases for increasing ambient turbulence due to the proportionality of k1 to TIamb in Eq. (1). Although a different dataset might produce other calibration parameters, it does not necessarily imply an improved accuracy of the simulated wake fields in the FFoR. Particularly, we expect our calibration and that of Madsen et al. (2010), Larsen et al. (2013), and IEC (2019) to provide similar wake fields in the FFoR for low ambient turbulence (i.e. ≤7 %), as illustrated by the blue lines in Fig. 9.
We recommend conducting DWM model calibrations using two-dimensional high-temporal- and high-spatial-resolution measurements of the wake field. Such resolution can be achieved by research-based nacelle lidars today, which allow the resolution of wake flow features including the wake deficit and wake-added turbulence profiles in the MFoR for model calibration. When using power production data for such calibrations (the IEC standard values are based on this approach), we are unable to distinguish between uncertainties from inaccurate wake deficit predictions and those from erroneous wake meandering representations. The latter play an important role in the accuracy of the fully resolved wake fields, which are inputs to aeroelastic simulations.
As wind turbines are typically spaced 5 D and beyond, it is recommended to focus future measurement campaigns and model validation studies on including those regions. Additional high-resolution datasets will allow a more comprehensive validation analysis to verify that the model proposed in this work is applicable for a wide range of conditions (e.g. turbine types and site conditions). To this extent, Bayesian inference is a valuable approach for updating the PDFs of model parameters estimated within this work by directly including data from future observations while retaining information from the earlier observations. Further, power and load validation analyses using the proposed calibration parameters at multiple sites can increase the confidence in our calibration methodology.
Characterizing wake turbulence using lidars is challenging due to the limited sampling frequency and probe volume effects (Peña et al., 2017). We demonstrate the usefulness of Doppler radial velocity spectra to compute unfiltered LOS variances in wake conditions. In addition to the enhanced turbulence intensity, the wake turbulence is characterized both by being highly isotropic and reduced turbulence length scale compared to the ambient turbulence (Madsen et al., 2005). These turbulence characteristics were not investigated in this work.
We analysed high-spatial- and high-temporal-resolution SpinnerLidar measurements of the wake field collected at the SWiFT facility and derived wake features such as the wake deficit, wake-added turbulence, and wake meandering under varying atmospheric-stability conditions, inflow wind speeds, and downstream distances. The SpinnerLidar-estimated wake characteristics computed in the MFoR were used to determine uncertainties in the DWM model parameters using Bayesian inference. The uncertainties in model parameters were propagated to predict fully resolved wake flow fields in the FFoR. This approach allowed us to quantify uncertainties in the DWM-model-simulated wake fields and to investigate the sensitivity of the model parameters to flow features that primarily affect power and load predictions.
The SpinnerLidar-derived wake deficit profiles revealed the strong impact of atmospheric stability on wake evolution. In particular, we observed the faster recovery of the deficit under unstable compared to stable regimes as higher turbulence intensities characterized the former. These effects were accurately reproduced by the eddy viscosity term of the DWM model with the inferred parameters for distances beyond 4 D. Our results indicate that the currently adopted parameters in the IEC standard lead to conservative velocity deficit predictions (up to 18 % for moderate to high ambient turbulence TIamb≥12 %) at distances up to 5 D behind the rotor.
We proposed and verified an improved semi-empirical formulation of the wake-added turbulence model that captured the effects of the atmospheric shear and the ambient turbulence on the spatial re-distribution of the wake turbulence observed at 2.5 D. We also demonstrated that the wake meandering is the major source of added turbulence in the wake region.
The underlying hypothesis of the DWM model (i.e. wakes are advected passively by the large eddies in the incoming wind field) was verified by means of the SpinnerLidar-tracked meandering time series. The spectral analysis indicated that large eddies associated with sizes larger than 2 D are responsible for the wake meandering; however, the large-eddy “definition” had only marginal effects on the predicted wake fields. Accurate tracking of the wake centre position was the most influential factor in simulating wake flow fields accurately. We expect it to also play a central role in the accuracy of power and load predictions.
In a future study, we will quantify uncertainties in power and load predictions based on the proposed calibration at two different sites, the SWiFT facility and the Nørrekær Enge wind farm in Denmark (Conti et al., 2020b).
Figure A1 compares measured and predicted time series of the wake meandering in the lateral direction under varying stability conditions. The DWM model predictions are computed by applying the filtering cut-off frequency with the stochastic definition of the large-scale eddies (Dm). A reasonable agreement between the two signals is found, where the major wake displacements are captured by the meandering model. Note that improved correlation can be achieved by utilizing a reduced advection wind speed in the time-lag parameter in Eq. (4) compared to the ambient wind velocity. Nevertheless, the largest observed movements are induced by the yawing of the WTGa1, which is fairly frequent within the analysed dataset. Further, the 95 % confidence interval of the meandering predictions obtained by propagating uncertainties in the Dm definition is nearly negligible. This indicates that uncertainties in the filtering cut-off frequency used in the DWM model formulation have a marginal effect on the accuracy of wake simulations.
Sensitivity analysis is conducted to identify the most important parameters affecting the accuracy of wake simulations. Here, we only investigate the uncertainties in the model parameters, which are listed in Table 4, and assume that the inflow conditions are perfectly prescribed. We employ a variance-based sensitivity method and compute total Sobol indices (Sobol, 2001; Saltelli et al., 2010). The Sobol sensitivity decomposes the variance of the response (e.g. Ueff and σU, eff) into contributions from input parameters and associated interactions.
The Sobol indices computed from wake simulations at five rotor diameters behind the rotor are illustrated in Fig. B1. Note that uncertainties in wake centre locations (ym, ϵ and zm, ϵ) have a similar influence as for the calibration parameters (i.e. k1 and k2) on the predicted flow features such as Ueff and σu, eff. Overall, tracking the wake meandering in both lateral and vertical directions is of primary importance in wake field representations and therefore in power and load validations. The sensitivity of the wake-added turbulence model parameters as well as the large-scale-eddy definition (Dm) is marginal. It is inferred that these parameters can be considered to be deterministic without compromising the accuracy of wake simulations.
The SpinnerLidar and meteorological mast data can be requested from the authors.
DC, ND, AP, and TH participated in the conception and design of the work. DC conducted the data analysis and wrote the draft manuscript. ND and AP supported the overall analysis and critically revised the manuscript. TH contributed with the acquisition of the dataset and critically revised the manuscript.
The authors declare that they have no conflict of interest.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Sandia National Laboratories is a multi-mission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC., a wholly owned subsidiary of Honeywell International, Inc., for the US Department of Energy’s National Nuclear Security Administration under contract DE-NA-0003525.
This paper was edited by Sandrine Aubrun and reviewed by two anonymous referees.
International Standard IEC61400-13: Wind turbines – Part 13: Measurement of mechanical loads, Standard, International Electrotechnical Commission (IEC), 2015. a
International Standard IEC61400-12-1: Wind energy generation systems – Part 12-1: Power performance measurements of electricity producing wind turbines, Standard, International Electrotechnical Commission (IEC), 2017. a
Ainslie, J. F.: Calculating the flowfield in the wke of wind turbines, J. Wind. Eng. Ind. Aerod., 27, 213–224, 1987. a
Aitken, M. L., Banta, R. M., Pichugina, Y. L., and Lundquist, J. K.: Quantifying wind turbine wake characteristics from scanning remote sensor data, J. Atmos. Ocean. Technol., 31, 765–787, https://doi.org/10.1175/JTECH-D-13-00104.1, 2014. a
Angelou, N. and Sjöholm, M.: UniTTe WP3/MC1: Measuring the inflow towards a Nordtank 500 kW turbine using three short-range WindScanners and one SpinnerLidar, DTU Wind Energy, Roskilde, Denmark, 2015. a
Berg, J., Bryant, J., LeBlanc, B., Maniaci, D. C., Naughton, B., Paquette, J. A., Resor, B. R., White, J., and Kroeker, D.: Scaled Wind Farm Technology Facility Overview, 32nd Asme Wind Energy Symposium, 13–17 January 2014, National Harbor, Maryland, https://doi.org/10.2514/6.2014-1088, 2014. a
Box, G. and Tiao, G.: Bayesian inference in statistical analysis, Addison Wesley, Voorburg, the Netherlands, 1973. a
Chamorro, L. P. and Porte-Agel, F.: A Wind-Tunnel Investigation of Wind-Turbine Wakes: Boundary-Layer Turbulence Effects, Bound.-Lay. Meteorol., 132, 129–149, https://doi.org/10.1007/s10546-009-9380-8, 2009. a
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, large-eddy simulation, and field data at the egmond aan Zee offshore wind plant, 33rd Wind Energy Symposium, 5–9 January 2015, Kissimmee, Florida, 20 pp., 2015. a, b
Conti, D., Dimitrov, N., Peña, A., and Herges, T.: Wind turbine wake characterization using the SpinnerLidar measurements, J. Phys. Conf. Ser., 1618, 062040, https://doi.org/10.1088/1742-6596/1618/6/062040, 2020a. a, b, c, d, e, f, g, h, i, j, k, l
Conti, D., Dimitrov, N., and Peña, A.: Aeroelastic load validation in wake conditions using nacelle-mounted lidar measurements, Wind Energ. Sci., 5, 1129–1154, https://doi.org/10.5194/wes-5-1129-2020, 2020b. a, b
Conti, D., Pettas, V., Dimitrov, N., and Peña, A.: Wind turbine load validation in wakes using wind field reconstruction techniques and nacelle lidar wind retrievals, Wind Energ. Sci., 6, 841–866, https://doi.org/10.5194/wes-6-841-2021, 2021. a, b, c
Debnath, M., Doubrawa, P., Herges, T., Martínez-Tossas, L. A., Maniaci, D. C., and Moriarty, P.: Evaluation of Wind Speed Retrieval from Continuous-Wave Lidar Measurements of a Wind Turbine Wake Using Virtual Lidar Techniques, J. Phys. Conf. Ser., 1256, 012008, https://doi.org/10.1088/1742-6596/1256/1/012008, 2019. a, b
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/1742-6596/1256/1/012024, 2019. a, b, c
Doubrawa, P., Quon, E. W., Martinez-Tossas, L. A., Shaler, K., Debnath, M., Hamilton, N., Herges, T. G., Maniaci, D., Kelley, C. L., Hsieh, A. S., Blaylock, M. L., van der Laan, P., Andersen, S. J., Krueger, S., Cathelain, M., Schlez, W., Jonkman, J., Branlard, E., Steinfeld, G., Schmidt, S., Blondel, F., Lukassen, L. J., and Moriarty, P.: Multimodel validation of single wakes in neutral and stratified atmospheric conditions, Wind Energy, 23, 2027–2055, https://doi.org/10.1002/we.2543, 2020. a, b, c, d, e
Espana, G., Aubrun, S., Loyer, S., and Devinant, P.: Spatial study of the wake meandering using modelled wind turbines in a wind tunnel, Wind Energy, 14, 923–937, https://doi.org/10.1002/we.515, 2011. a
Frandsen, S. T.: Turbulence and turbulence-generated structural loading in wind turbine clusters, Risø National Laboratory, Roskilde, Denmark, 2007. a
Fuertes, F. C., Markfort, C. D., and Porteacute-Agel, F.: Wind Turbine Wake Characterization with Nacelle-Mounted Wind Lidars for Analytical Wake Model Validation, Remote Sens., 10, 668, https://doi.org/10.3390/rs10050668, 2018. a, b
Galinos, C., Dimitrov, N. K., Larsen, T. J., Natarajan, A., and Hansen, K. S.: Mapping Wind Farm Loads and Power Production – A Case Study on Horns Rev 1, J. Phys. Conf. Ser., 753, 032010, https://doi.org/10.1088/1742-6596/753/3/032010, 2016. a, b
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/1742-6596/1256/1/012005, 2019. a, b, c
Herges, T. G., Maniaci, D. C., Naughton, B. T., Mikkelsen, T. K., and Sjöholm, M.: High resolution wind turbine wake measurements with a scanning lidar: Paper, J. Phys. Conf. Ser., 854, 012021, https://doi.org/10.1088/1742-6596/854/1/012021, 2017. a, b, c, d, e, f, g
Herges, T. G., Berg, J. C., Bryant, J. T., White, J. R., Paquette, J. A., and Naughton, B. T.: Detailed analysis of a waked turbine using a high-resolution scanning lidar, J. Phys. Conf. Ser., 1037, 072009, https://doi.org/10.1088/1742-6596/1037/7/072009, 2018. a, b, c, d
Hoffman, M. D. and Gelman, A.: The no-U-turn sampler: Adaptively setting path lengths in Hamiltonian Monte Carlo, J. Mach. Learn. Res., 15, 1593–1623, 2014. a
Iungo, G. V. and Porté-Agel, F.: Volumetric lidar scanning of wind turbine wakes under convective and neutral atmospheric stability regimes, J. Atmos. Ocean. Technol., 31, 2035–2048, https://doi.org/10.1175/JTECH-D-13-00252.1, 2014. a
Ivanell, S. S. A., Henningson, D. P., and Joachim, P. P.: Numerical computations of wind turbine wakes, KTH, Stockholm, 2009. 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. Energ., 134, 021012, https://doi.org/10.1115/1.4006038, 2012. a, b, c, d, e, f, g, h, i, j
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, 2014a. a, b, c
Keck, R.-E., Mikkelsen, R. F., Troldborg, N., de Maré, M., and Hansen, K. S.: Synthetic atmospheric turbulence and wind shear in large eddy simulations of wind turbine wakes, Wind Energy, 17, 1247–1267, https://doi.org/10.1002/we.1631, 2014b. a
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 build-up in a row of wind turbines, Wind Energy, 18, 111–132, https://doi.org/10.1002/we.1686, 2015. a, b, c, d, e, f
Käsler, Y., Rahm, S., Simmet, R., and Kühn, M.: Wake measurements of a multi-MW wind turbine with coherent long-range pulsed doppler wind lidar, J. Atmos. Ocean. Technol., 27, 1529–1532, https://doi.org/10.1175/2010JTECHA1483.1, 2010. a
Larsen, G., Ott, S., Liew, J., van der Laan, M., Simon, E., R.Thorsen, G., and Jacobs, P.: Yaw induced wake deflection – a full-scale validation study, J. Phys. Conf. Ser., 1618, 062047, https://doi.org/10.1088/1742-6596/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, DTU Wind Energy, Roskilde, Denmark, 2007. a
Larsen, G. C., Machefaux, E., and Chougule, A.: Wake meandering under non-neutral atmospheric stability conditions – Theory and facts, J. Phys. Conf. Ser., 625, 012036, https://doi.org/10.1088/1742-6596/625/1/012036, 2015. a, b
Larsen, T. J. and Hansen, A. M.: How 2 HAWC2, the user's manual, Risø National Laboratory, Roskilde, Denmark, 2007. a
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, f, g, h, i, j, k, l, m, n, o, p, q, r, s
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/1742-6596/1618/6/062002, 2020. a
Machefaux, E., Larsen, G. C., Troldborg, N., Gaunaa, M., and Rettenmeier, A.: Empirical modeling of single-wake advection and expansion using full-scale pulsed lidar-based measurements, Wind Energy, 18, 2085–2103, https://doi.org/10.1002/we.1805, 2015. a, b
Machefaux, E., Larsen, G. C., Koblitz, T., Troldborg, N., Kelly, M. C., Chougule, A., 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, c, d
Madsen, A. H., Larsen, G. C., and Thomsen, K.: Wake flow characteristics in low ambient turbulence conditions, Proceedings (cd-rom), Copenhagen Offshore Wind 2005, 26–28 October 2005, 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. Energ., 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, q, r
Mann, J.: The spatial structure of neutral atmospheric surface-layer turbulence, J. Fluid Mech., 273, 141–168, 1994. a
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. Technol., 27, 959–976, https://doi.org/10.1175/2010jtecha1389.1, 2010. a, b
Meyer Forsting, A. R., Troldborg, N., and Borraccino, A.: Modelling lidar volume-averaging and its significance to wind turbine wake measurements: Paper, J. Phys. Conf. Ser., 854, 012014, https://doi.org/10.1088/1742-6596/854/1/012014, 2017. a
Monin, A. and Obukhov, A.: Osnovnye zakonomernosti turbulentnogo peremeshivanija v prizemnom sloe atmosfery (Basic laws of turbulent mixing in the atmosphere near the ground), Trudy Geofiz. Inst, 24, 163–187, 1954. a
Peña, A., Sjöholm, M., Mikkelsen, T. K., and Hasager, C. B.: Inflow characterization using measurements from the SpinnerLidar: the ScanFlow experiment, J. Phys. Conf. Ser., 1037, 052027, https://doi.org/10.1088/1742-6596/1037/5/052027, 2018. a, b, c
Reinwardt, I., Gerke, N., Dalhoff, P., Steudel, D., and Moser, W.: Validation of wind turbine wake models with focus on the dynamic wake meandering model, J. Phys. Conf. Ser., 1037, 072028, https://doi.org/10.1088/1742-6596/1037/7/072028, 2018. a, b, c, d, e, f
Reinwardt, I., Schilling, L., Dalhoff, P., Steudel, D., and Breuer, M.: Dynamic wake meandering model calibration using nacelle-mounted lidar systems, Wind Energ. Sci., 5, 775–792, https://doi.org/10.5194/wes-5-775-2020, 2020. a, b, c, d, e, f, g, h, i
Saltelli, A., Annoni, P., Azzini, I., Campolongo, F., Ratto, M., and Tarantola, S.: Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index, Comput. Phys. Commun., 181, 259–270, https://doi.org/10.1016/j.cpc.2009.09.018, 2010. a
Subramanian, B., Chokani, N., and Abhari, R.: Impact of atmospheric stability on wind turbine wake evolution, J. Wind Eng. Ind. Aerod., 176, 174–182, https://doi.org/10.1016/j.jweia.2018.03.014, 2018. a
Thomsen, K., Madsen Aagaard, H., Larsen, G. C., and Larsen, T. J.: Comparison of methods for load simulation for wind turbines operating in wake, J. Phys. Conf. Ser., 75, 012072, https://doi.org/10.1088/1742-6596/75/1/012072, 2007. a
Trabucchi, D., Trujillo, J.-J., and Kuehn, M.: Nacelle-based Lidar Measurements for the Calibration of a Wake Model at Different Offshore Operating Conditions, 14th Deep Sea Offshore Wind R&d Conference, Eera Deepwind '2017, 137, 77–88, https://doi.org/10.1016/j.egypro.2017.10.335, 2017. 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: two-dimensional scanning, Wind Energy, 14, 61–75, https://doi.org/10.1002/we.402, 2011. a, b, c, d, e
Vollmer, L., Steinfeld, G., Heinemann, D., and Kühn, M.: Estimating the wake deflection downstream of a wind turbine in different atmospheric stabilities: an LES study, Wind Energ. Sci., 1, 129–141, https://doi.org/10.5194/wes-1-129-2016, 2016. a
Walker, K., Adams, N., Gribben, B., Gellatly, B., Nygaard, N. G., Henderson, A., Marchante Jimémez, M., Schmidt, S. R., Rodriguez Ruiz, J., Paredes, D., Harrington, G., Connell, N., Peronne, O., Cordoba, M., Housley, P., Cussons, R., Håkansson, M., Knauer, A., and Maguire, E.: An evaluation of the predictive accuracy of wake effects models for offshore wind farms, Wind Energy, 19, 979–996, https://doi.org/10.1002/we.1871, 2016. a
Yang, X. and Sotiropoulos, F.: Wake characteristics of a utility-scale wind turbine under coherent inflow structures and different operating conditions, Physical Review Fluids, 4, 024604, https://doi.org/10.1103/PhysRevFluids.4.024604, 2019. 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, 2020b. a, b, c