Research article 16 Aug 2018
Research article  16 Aug 2018
From standard wind measurements to spectral characterization: turbulence length scale and distribution
 Wind Energy Department, Risø Lab./Campus, Danish Technical University, Roskilde 4000, Denmark
 Wind Energy Department, Risø Lab./Campus, Danish Technical University, Roskilde 4000, Denmark
Correspondence: Mark Kelly (mkel@dtu.dk)
Hide author detailsCorrespondence: Mark Kelly (mkel@dtu.dk)
In wind energy, the effect of turbulence upon turbines is typically simulated using wind “input” time series based on turbulence spectra. The velocity components' spectra are characterized by the amplitude of turbulent fluctuations, as well as the length scale corresponding to the dominant eddies. Following the IEC standard, turbine load calculations commonly involve use of the Mann spectraltensor model to generate time series of the turbulent threedimensional velocity field. In practice, this spectraltensor model is employed by adjusting its three parameters: the dominant turbulence length scale L_{MM} (peak length scale of an undistorted isotropic velocity spectrum), the rate of dissipation of turbulent kinetic energy ε, and the turbulent eddylifetime (anisotropy) parameter Γ. Deviation from “ideal” neutral sheared turbulence – i.e., for nonzero heat flux and/or heights above the surface layer – is, in effect, captured by setting these parameters according to observations.
Previously, sitespecific $\mathit{\{}{L}_{\text{MM}},\mathit{\epsilon},\mathrm{\Gamma}\mathit{\}}$ values were obtainable through fits to measured threedimensional velocity component spectra recorded with sample rates resolving the inertial range of turbulence (≳1 Hz); however, this is not feasible in most industrial wind energy projects, which lack multidimensional sonic anemometers and employ loggers that record measurements averaged over intervals of minutes. Here a form is derived for the shear dependence implied by the eddylifetime prescription within the Mann spectraltensor model, which leads to derivation of useful forms of the turbulence length scale. Subsequently it is shown how L_{MM} can be calculated from commonly measured sitespecific atmospheric parameters, namely mean wind shear (dU∕dz) and standard deviation of streamwise fluctuations (σ_{u}). The derived L_{MM} can be obtained from standard (10 min average) cup anemometer measurements, in contrast with an earlier form based on friction velocity.
The new form is tested across several different conditions and sites, and it is found to be more robust and accurate than estimates relying on friction velocity observations. Assumptions behind the derivations are also tested, giving new insight into rapiddistortion theory and eddylifetime modeling – and application – within the atmospheric boundary layer. The work herein further shows that distributions of turbulence length scale, obtained using the new form with typical measurements, compare well with distributions P(L_{MM}) obtained by fitting to spectra from researchgrade sonic anemometer measurements for the various flow regimes and sites analyzed. The new form is thus motivated by and amenable to sitespecific probabilistic loads characterization.
Of the atmospheric parameters which are generally input into (or required by) wind turbine load calculation codes, several stand out due to their prominence in load contributions: the “mean” wind speed U, the standard deviation of streamwise turbulent velocity σ_{u}, the shear dU∕dz or shear exponent α, and the characteristic turbulence length scale L corresponding to the most energetic turbulent motions (Wyngaard, 2010). Dimitrov et al. (2015) explored the importance of shear (α); Dimitrov et al. (2017) found that both fatigue and extreme turbine loads can be sensitive to L in addition to the dominant influences of mean wind speed U and streamwise turbulence “strength” σ_{u}^{1}. These are also consistent with the earlier finding of Sathe et al. (2013) that stability could affect fatigue loads through α and σ_{u}.
Within the context of obtaining sitedependent statistics of the most crucial loaddriving parameters ($U,{\mathit{\sigma}}_{u},\mathit{\alpha},L$) from conventional industrial wind measurements, this work focuses on the one parameter which has thus far been most difficult to measure: the turbulence length scale L. The turbulence length scale corresponds to the “energycontaining subrange” of turbulent velocity fluctuations associated with the peak of the streamwise velocity spectrum, which contribute most to turbulent kinetic energy (and σ_{u}) – and which can dominate the turbulence contribution to wind turbine loads. Measurements used in wind energy are usually stored as 10 min statistics (average and standard deviation of wind speed and direction), so one cannot obtain turbulence spectra from them, nor can one calculate integral time or length scales from such observations.
Because of its widespread use in the wind industry and its inclusion in the IEC 61400–1, Edition 3 (2005) standard on design requirements for wind turbines, here we consider the spectral turbulence model of Mann (1994) and L as prescribed for this model. Within the “Mann model”, which uses rapiddistortion theory (RDT) to account for shearinduced distortion of isotropic turbulence Pope (2000); Savill (1987), there is also a prescription for the scaledependent time over which turbulent eddies of a given size are distorted. This timescale is key to proper representation of atmospheric turbulence and reproduction of component spectra via RDT. However, the eddy lifetime was not directly derived, but rather cleverly prescribed, by Mann (1994). Concurrent to and independent of the work herein, de Mare and Mann (2016) also derived some relations to create a model for timevarying eddy lifetime. The present article provides direct derivation of the eddy lifetime, which results in a relation between the three (spectral) parameters of the Mann model and measurable quantities. More importantly, the derivations here include connection of the turbulence length scale to routinely available quantities from typical 10 min industrial wind records. The turbulence length scale is in fact that corresponding to the von Kármán (1948) spectral form, and thus the relation here is applicable to other turbulence models used in wind engineering, such as those relying on the Kaimal et al. (1972) spectrum.
After deriving the eddy lifetime and giving subsequent expressions for the turbulence length scale, this article proceeds to validation of the underlying assumptions. Constraints implied by fitting the Mann model to measured spectra in nonneutral conditions, given eddy lifetime and mixinglength relations, are also tested. This includes dependence of predicted velocity variance on model anisotropy parameter (Γ), as well as implications in the surface layer and connection to previous findings in boundarylayer meteorology. Finally, the length scale obtained from conventional 10 min wind measurements via the new expression is compared to the length scale found from fits of Mannmodel output to measured component spectra; this is done using data from multiple sites, representing several types of site conditions.
Relation of the turbulence length (spectral “peak”) scale to measurable statistics is possible through the eddylifetime form of Mann (1994), where the latter is defined in terms of the isotropic von Kármán spectrum that is distorted using RDT.
2.1 Eddy lifetime
A number of forms exist to estimate eddy lifetime τ_{e}, though these can be generally expressed as the ratio of a length scale (taken as the reciprocal of wavenumber, k^{−1}) to a velocity scale which follows from some integrated form of the (scalar) kinetic energy spectrum E(k):
where the characteristic velocity scale can be generically described by
In contrast to the “coherencedestroying diffusion time” of ComteBellot and Corrsin (1971) and the reciprocal of eddydamping rates from Lesieur (1990), for use with rapiddistortion theory Mann (1994) chose an eddy lifetime that depends on eddy size (wavenumber) according to
i.e., equivalent to p=0 in terms of Eq. (1). The choice of Eq. (2) for eddy lifetime was found to behave more reasonably than both the ComteBellot and Corrsin (1971) “diffusion time” (where p=1)^{2}, as well as the timescale $\left[{k}^{\mathrm{3}}E\right(k){]}^{\mathrm{1}/\mathrm{2}}$ (which in the inertial range is equivalent to $p=\mathrm{1}$)^{3} implicit in eddydamped quasinormal Markovian models (Andre and Lesieur, 1977; Lesieur, 1990); both of the latter lifetime models do not (reliably) integrate to give finite ${\mathit{\sigma}}_{u}^{\mathrm{2}}$.
Mann (1994) rewrites τ_{M} as
where _{2}F_{1} is Gauss' hypergeometric function (Abramowitz and Stegun, 1972)^{4}, and L_{MM} is the turbulence length scale associated with the peak of the turbulent kinetic energy spectrum E(k) as in Eq. (2). The eddy lifetime definition (3) is used in practical implementation of the spectral tensor model (Mann, 2000), and it notably defines a parameter of this model: the eddy lifetime factor Γ, also known as the anisotropy factor. The Mann (1994) spectraltensor model employs RDT, whereby the shear dU∕dz distorts turbulence from an isotropic state, based on an initial turbulent kinetic energy spectrum of the von Kármán form
where α=1.7 (von Kármán, 1948). This in effect defines the length scale L_{MM} through the peak of the initial spectrum^{5}. Using Eq. (4) in the proportionality expression (2) produces
where we have introduced the proportionality constant c_{τ} to write the result of integrating the proportionality relation (2) as an equation. Now τ_{M} can be seen to depend upon k, L_{MM}, and ε. The eddy lifetime can be reduced and clarified via ${}_{\mathrm{2}}\phantom{\rule{0.125em}{0ex}}{F}_{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}\left\{\frac{\mathrm{1}}{\mathrm{3}},\frac{\mathrm{17}}{\mathrm{6}};\frac{\mathrm{4}}{\mathrm{3}},(k{L}_{\text{MM}}{)}^{\mathrm{2}}\right\}\simeq [\mathrm{1}+\mathrm{3.07}(k{L}_{\text{MM}}{)}^{\mathrm{2}}{]}^{\mathrm{1}/\mathrm{3}}$ to give the more transparent von Kármánlike form^{6}
Since Eqs. (3) and (5) are equal, we have an expression relating the Mannmodel parameters to the shear dU∕dz:
The expression (7) can be made yet more useful to relate the turbulent length scale to measurable parameters, as shown in Sect. 2.2.
Eddy lifetime and equilibrium
The parameters $\mathit{\{}\mathit{\epsilon},\mathrm{\Gamma},{L}_{\text{MM}}\mathit{\}}$ are sitedependent and in practice have been obtained from measurements through fits of the model output to observed spectra (Mann, 2000), relying on (at least three of) F_{11}, F_{13}, F_{33}, and F_{22} (Dimitrov et al., 2017; Sathe et al., 2013). The model starts with an (undistorted) isotropic incompressible turbulence spectral tensor:
where E(k) is taken to be E_{vK}(k) (shown in Eq. 4), then the Φ_{ij} are distorted – i.e., the rapiddistortion equations are solved – per (threedimensional) wavenumber over a time τ_{M}(k) via RDT.
The rapiddistortion equations do not explicitly solve for production of normal stresses (which sum to twice the turbulent kinetic energy) or shear stress, though they do include terms that perturb the stresses^{7} to account for the (anisotropic) effect of a constant shear dU∕dz. Further, the RDT discussed here does not include dissipation (Mann, 1994; Pope, 2000); instead, in the spectraltensor model the dissipation rate of turbulent kinetic energy ε is a parameter giving the amplitude of the undistorted (initial) spectrum via Eq. (4). In practice ε is obtained via fits of precalculated Mannmodel output to measured spectra. So ε in effect gives the inertialrange amplitudes of the distorted velocity component spectra, which have been distorted for a time τ_{M}(k). From Eq. (3) one sees that the parameter Γ serves as a factor that determines the amount of distortion and associated anisotropy: increasing Γ corresponds to longer distortion time τ_{M} and thus more anisotropy, with Γ=0 corresponding to isotropy (zero distortion of the initial isotropic Φ_{ij}). The separation between the peaks of the different component spectra increases with Γ; the spectral peak of F_{33} is at higher wavenumbers (smaller scales) than the F_{22} peak, which is at higher wavenumbers than the peak of F_{11} (Mann, 1994).
A stationary equilibrium result is achieved via the eddylifetime prescription together with rapid distortion of the isotropic spectral tensor – with τ_{M} and (initial) inertialrange amplitudes depending on ε via Eqs. (3)–(4) and (7), whereby shearproduction of TKE is in effect balanced by dissipation. That is, the resultant shear stress 〈uw〉 (expressible now in terms of ε) can be multiplied by $\mathrm{2}\partial U/\partial z$ to give the implied production rate of 〈uu〉, which with vv and ww (through Γ) gives the implied TKE production rate, amounting to P=ε; such an equilibrium, enforced by τ_{M}, can also be inferred from de Mare and Mann (2016).
2.2 Characteristic length scale
Noting that the spectrum of a variable integrates to the variance of said variable, then invoking Eq. (8) with the isotropic von Kármán form Eq. (4) for E(k) and exploiting ${F}_{\mathrm{11}}\left({k}_{\mathrm{1}}\right)=\int \phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}\int {\mathrm{\Phi}}_{\mathrm{11}}\mathrm{d}{k}_{\mathrm{2}}\mathrm{d}{k}_{\mathrm{3}}$, one obtains the isotropic streamwise turbulence variance
which is the undistorted streamwise variance. The factor 0.69 is the numerical value of $\frac{\mathrm{9}}{\mathrm{55}}\sqrt{\mathit{\pi}}{f}_{\mathrm{\Gamma}}\left(\frac{\mathrm{1}}{\mathrm{3}}\right)/{f}_{\mathrm{\Gamma}}\left(\frac{\mathrm{5}}{\mathrm{6}}\right)$, and f_{Γ}(x) is the Euler gamma function (Abramowitz and Stegun, 1972). Then using Eq. (9) in Eq. (7) we get a relation for the isotropic (undistorted) turbulence length scale implied by the lifetime model (3),
where the leading term in parentheses is expected to be of the order of 1.
2.2.1 Relation to observations
Peña et al. (2010) suggested that the Mannmodel length scale is proportional to the classic mixing length ${\mathrm{\ell}}_{*}\equiv {u}_{*}/(\mathrm{d}U/\mathrm{d}z)$ multiplied by an empirical constant,
where they assign c_{m}=1.7. However, we find from observations that on average c_{m}≈2.3 over flat land, i.e., $\langle {L}_{\text{MM}}/{\mathrm{\ell}}_{*}\rangle =\mathrm{2.3}$ (see next section). Combining Eqs. (10)–(11) one sees that c_{τ} decreases with the relative magnitude of measured shear stress (as ${\mathit{\sigma}}_{\text{iso}}/{u}_{*,\text{obs}}$); this is also expressed usefully through the measured ratio of streamwise fluctuation amplitude to friction velocity:
From the above and Eq. (10) one subsequently then finds
For constant (${\mathit{\sigma}}_{u,\text{obs}}/{u}_{*,\text{obs}}$), Eq. (13) implies that the turbulence scale L_{MM} can be expressed independently of Γ, given σ_{u,obs} and dU∕dz.
Caughey et al. (1979) reported the mean profile of ${\mathit{\sigma}}_{u}^{\mathrm{2}}\left(z\right)$ from the seminal “Kansas experiment”, showing that $({\mathit{\sigma}}_{u}/{u}_{*}{)}_{\mathrm{0}}^{\mathrm{2}}\approx \mathrm{5}$–6 in the homogeneous atmospheric surface layer (their Fig. 5). The corresponding value of $({\mathit{\sigma}}_{u}/{u}_{*}{)}_{\mathrm{0}}$ is approximately 2.3; thus, if c_{m}≈2.3 as well, then Eq. (12) reduces to
Given the definition of c_{τ} through Eq. (7), c_{τ} is a constant; since Eq. (9) shows σ_{iso} is independent of Γ, then ${\mathit{\sigma}}_{u,\text{obs}}\propto \mathrm{\Gamma}$. Consistent with this argument, Eq. (13) reduces to
which is also evident inserting Eq. (14) into Eq. (10). Using Eq. (15), L_{MM} can simply be diagnosed from typical measurements, e.g., 10 min average cupanemometer output at two (or more) heights. The length L_{MM} can also be cast in terms of variables commonly used in wind engineering, notably the turbulence intensity I_{u} and shear exponent α. Invoking $\mathrm{d}U/\mathrm{d}z=\mathit{\alpha}U/z$ (Kelly et al., 2014a) and defining ${I}_{\text{obs}}\equiv {\mathit{\sigma}}_{u,\text{obs}}/U$, then Eq. (15) becomes
2.2.2 Modeled spectra: covariances, anisotropy, and Γ
The spectral Mann model (“MM”) distorts the isotropic von Kármán spectral tensor (Φ_{ij}(k), Eq. 4), per wavenumber via rapiddistortion theory over the wavenumberdependent eddy lifetime τ_{M}, such that the component spectra become anisotropic at wavenumbers outside (lower than) the inertial range; the degree of distortion – and thus anisotropy – are consequently represented by the eddylifetime parameter Γ. Above we showed via mixinglength arguments that L_{MM} is independent of Γ, resulting in Eq. (15). Possible Γ dependences can also be examined by considering the shear stress
obtained from the modeled spectral tensor component ${F}_{\mathrm{13}}\left({k}_{\mathrm{1}}\right)=\int \phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}\int {\mathrm{\Phi}}_{\mathrm{13}}\mathrm{d}{k}_{\mathrm{2}}\mathrm{d}{k}_{\mathrm{3}}$, which is expected to be a function of Γ. Indeed Mann (1994) shows this to be the case, with modeled stress $\langle uw{\rangle}_{\text{MM}}/{\mathit{\sigma}}_{\text{iso}}^{\mathrm{2}}$ varying almost linearly between 0 and −1 for $\mathrm{0}<\mathrm{\Gamma}<\mathrm{5}$; then ${u}_{*,\text{MM}}^{\mathrm{2}}/{\mathit{\sigma}}_{\text{iso}}^{\mathrm{2}}\approx \mathrm{\Gamma}/\mathrm{5}$. Subsequently from Eq. (12) one has
for c_{m}≈2.3, in analogy with Eq. (14); thus we expect ${u}_{*,\text{obs}}\propto \mathrm{\Gamma}$, similar to the expected behavior of ${\mathit{\sigma}}_{u,\text{obs}}\propto \mathrm{\Gamma}$ following Eq. (14).
In addition to the approximate expression (18), which is based on the simplified relation ${u}_{*,\text{MM}}^{\mathrm{2}}/{\mathit{\sigma}}_{\text{iso}}^{\mathrm{2}}\approx \mathrm{\Gamma}/\mathrm{5}$, it is possible to derive an exact relation based on the Mannmodel shear stress (Eq. 17) – but this is cumbersome and analytically intractable. Though de Mare and Mann (2016) derived implicit expressions toward relating {Γ, dU∕dz, L_{MM}} to the eddy lifetime and integral of the modeled stress spectrum (Eq. 17), these must be evaluated numerically or graphically. An explicit expression corresponding to ${c}_{m}^{\mathrm{1}}={\mathrm{\ell}}_{*}/{L}_{\text{MM}}$ (like Eq. 11 here) was derived by de Mare and Mann (2016), but it depends on numerically integrating the stress spectrum.
As spectra fitted to Mannmodel outputs correspond to distorted anisotropic turbulence, and noting the Γ dependence of ${u}_{*,\text{MM}}$ discussed above, we expect σ_{u,MM} to also depend on Γ. From Fig. 4 of Mann (1994) we find ${\mathit{\sigma}}_{u,\text{MM}}^{\mathrm{2}}/{\mathit{\sigma}}_{\text{iso}}^{\mathrm{2}}\simeq (\mathrm{1}+\mathrm{0.14}{\mathrm{\Gamma}}^{\mathrm{2}})$, which for Γ≳2, the range corresponding to atmospheric boundary layer (ABL) observations (Sathe et al., 2013), becomes roughly ${\mathit{\sigma}}_{u,\text{MM}}\approx {\mathit{\sigma}}_{\text{iso}}(\mathrm{0.61}+\mathrm{0.3}\mathrm{\Gamma})$.
2.3 Ideal neutral surfacelayer implications
Within the atmospheric surface layer (ASL), in the homogeneous stationary limit under neutral conditions, $\mathrm{d}U/\mathrm{d}z\phantom{\rule{0.125em}{0ex}}\to {u}_{*}/\left(\mathit{\kappa}z\right)$ so that Eq. (11) reduces to ${L}_{\text{MM}}\to {c}_{m}\mathit{\kappa}z\approx \mathrm{0.92}z$. Similarly, in this “loglaw regime” ${\mathit{\epsilon}}_{\text{ASL,N}}={u}_{*}^{\mathrm{3}}/\left(\mathit{\kappa}z\right)$ so that Eq. (7) becomes ${\mathrm{\Gamma}}_{\text{ASL,N}}={c}_{\mathit{\tau}}(\mathrm{3}\mathit{\alpha}/\mathrm{2}{)}^{\mathrm{1}/\mathrm{2}}({L}_{\text{MM}}/\mathit{\kappa}z{)}^{\mathrm{2}/\mathrm{3}}$, or equivalently ${L}_{\text{MM}}{}_{\text{ASL,N}}=(\mathrm{3}\mathit{\alpha}/\mathrm{2}{)}^{\mathrm{3}/\mathrm{4}}\mathit{\kappa}z[\mathrm{\Gamma}/{c}_{\mathit{\tau}}{]}^{\mathrm{3}/\mathrm{2}}$, which via Eq. (12) can be written
Thus for ${c}_{m}={\mathit{\sigma}}_{\text{u,obs}}/{u}_{\text{*,obs}}$, we see that the Mann (1994) eddylifetime formulation (3) implies ${L}_{\text{MM}}\to \mathrm{1.1}\mathit{\kappa}z({\mathit{\sigma}}_{u,\text{obs}}/{\mathit{\sigma}}_{\text{iso}}{)}^{\mathrm{3}/\mathrm{2}}$ in the neutral ASL. Meanwhile, as noted just above, the mixinglength form (11) implies L_{MM}→c_{m}κz; this is consistent with Eq. (19) under the condition that $({\mathit{\sigma}}_{u,\text{obs}}/{\mathit{\sigma}}_{\text{iso}})\simeq ({c}_{m}/\mathrm{1.1}{)}^{\mathrm{2}/\mathrm{3}}$ or roughly ${\mathit{\sigma}}_{u,\text{obs}}\approx \mathrm{1.6}{\mathit{\sigma}}_{\text{iso}}$ for c_{m}≃2.3.
Since the choice of eddy lifetime form (3) leads to a sheardependent relation (7) between the spectraltensor model parameters, one obtains Eq. (10) for the undistorted (isotropic) length scale, with ${L}_{\text{MM}}\propto \phantom{\rule{0.125em}{0ex}}{\mathit{\sigma}}_{\text{iso}}(\mathrm{d}U/\mathrm{d}z{)}^{\mathrm{1}}$; further invoking a mixinglength argument then leads to a relation (15) for L_{MM} in terms of quantities that are directly measurable via standard windindustry (onedimensional cup) anemometers. Here we test Eq. (15) as well the assumptions leading to it, through measured wind speed, shear, and turbulent velocity component spectra. We also find a form for the distribution of L_{MM} over all conditions – as would be needed in practice to represent the turbulence length scales of flows experienced by wind turbines at a given site.
For the assumption testing in this section, the spectra used are measured via threedimensional sonic anemometers on the primary meteorological mast located at the Danish National Test Centre for Large Wind Turbines (Høvsøre), 1.75 km from the western coast of Denmark (Mann et al., 2005; Peña et al., 2016). The anemometers give 20 Hz samples of all three velocity components and temperature^{8} at heights of 10, 20, 40, 60, 80, and 100 m. This allows calculation of mean speeds, directions, and vertical shear of mean speed over individual 10 min records; in particular we focus on heights of z=20 m and z=80 m, as we are able to calculate shear at (across) these heights using the measurements at 10, 40, 60, and 100 m while also using the measured wind speed components and subsequent spectra at $z=\mathit{\{}\mathrm{20},\mathrm{80}\mathit{\}}$ m. The parameters $\mathit{\{}{L}_{\text{MM}},\mathrm{\Gamma},\mathit{\u03f5}\mathit{\}}$ are obtained via fits of precalculated Mannmodel spectra to the measured velocitycomponent and stress spectra F_{11}(f), F_{22}(f), F_{33}(f), and F_{13}(f); this is done via Taylor's hypothesis (${k}_{\mathrm{1}}=\mathrm{2}\mathit{\pi}f/U$) along with combined leastsquares fits (Chougule et al., 2017; Mann, 1994).
3.1 Testing of assumptions and predicted constraints
The implications of Eqs. (12)–(15) included the independence of L_{MM} and c_{τ} from Γ, as well as (for example) the expected dependence ${\mathit{\sigma}}_{u,\text{obs}}\propto \mathrm{\Gamma}{\mathit{\sigma}}_{\text{iso}}$. Indeed we find that L_{MM} is independent of Γ, with no significant statistical correlation: $\langle {L}_{\text{MM}}\mathrm{\Gamma}\rangle /\sqrt{\langle {L}_{\text{MM}}^{\mathrm{2}}\rangle \langle {\mathrm{\Gamma}}^{\mathrm{2}}\rangle}<\mathrm{0.15}$ for land or sea sectors at any given height. We also confirm that ${\mathit{\sigma}}_{u,\text{obs}}\propto \mathrm{\Gamma}{\mathit{\sigma}}_{\text{iso}}$, which is demonstrated by Figs. 1–2. The first figure displays the joint probability density $P({\mathit{\sigma}}_{u,\text{obs}},{\mathit{\sigma}}_{\text{iso}})$, where σ_{u,obs} is the streamwise turbulent variance measured in 10 min intervals, and σ_{iso} is calculated using Eq. (9) with L_{MM} and ϵ from spectral fits corresponding to the same intervals. One can see from Fig. 1 that σ_{u,obs} generally follows σ_{iso}, and we find ${\mathit{\sigma}}_{u,\text{obs}}/{\mathit{\sigma}}_{\text{iso}}\approx \mathrm{5}/\mathrm{3}$. Such evidence corresponds closely to the predicted constraint following Eq. (19) that ${\mathit{\sigma}}_{u,\text{obs}}/{\mathit{\sigma}}_{\text{iso}}$ should have a value of roughly 1.6 in the neutral surface layer; this is reasonable in the mean, since conditions on average are essentially neutral due to the shape of the stability distribution at Høvsøre (Kelly and Gryning, 2010). Figure 2 further shows that ${\mathit{\sigma}}_{u,\text{obs}}/{\mathit{\sigma}}_{\text{iso}}\propto \mathrm{\Gamma}$, consistent with c_{τ} being a constant independent of Γ following Eq. (14). The slope of the line in Fig. 2 also corresponds to the approximate Mannmodel behavior ${\mathit{\sigma}}_{u,\text{MM}}\approx {\mathit{\sigma}}_{\text{iso}}(\mathrm{0.61}+\mathrm{0.3}\mathrm{\Gamma})$ for Γ≳2, outlined at the end of Sect. 2.2.2 above.
Considering wind speeds in the typical turbine operating range of 4–25 m s^{−1}, the Høvsøre data also confirm that $\langle {\mathit{\sigma}}_{u}/{u}_{*}{\rangle}_{\text{obs}}\approx \mathrm{2.3}$, consistent with the findings of Caughey et al. (1979). Further, the data also show that 〈c_{m}〉≈2.3 so that Eq. (13) reduces approximately to Eq. (15). It is also found that the same approximate trends are seen when considering only U>7 m s^{−1} (not shown), but with slightly less scatter (narrower joint distributions) away from the predicted ${\mathit{\sigma}}_{u,\text{obs}}/{\mathit{\sigma}}_{\text{iso}}$ behaviors shown in Figs. 1–2 (dashed/dotted lines) and discussed above.
The data also show that ${\mathit{\sigma}}_{u}/{u}_{*}$ is not correlated with L_{MM}, whether we include all speeds or limit the wind speed range to 7–25 or 4–25 m s^{−1}. Thus this ratio can be treated as a constant in Eq. (13) for a given height (or throughout the surface layer), using Eq. (13) over a range of wind speeds.
3.2 Turbulence lengthscale distributions P(L_{MM})
The efficacy of using Eq. (15) to estimate the spectral length scale L_{MM} can be seen by considering Fig. 3. The figure displays the joint distribution of turbulence length scale at a height of z=80 m, i.e., $P\left({L}_{\text{MM,obs}},{\mathit{\sigma}}_{u,\text{obs}}\mathrm{d}U/\mathrm{d}z{}^{\mathrm{1}}\right)$; this is obtained through Eq. (15) from 10 min measurements and via fitting observed spectra. Figure 3 is usefully interpreted as the probabilityweighted performance of Eq. (15) for predicting L_{MM} (from σ_{u,obs} measured at z=80 m and the shear dU∕dz observed over z=60–100 m) versus the L_{MM} obtained from fits of the spectraltensor model to corresponding 10 min spectra. One sees a 1 : 1 relationship, particularly for the most commonly found values of the length scale; these L_{MM} values range ∼15–50 m^{9}. Compared to the scales calculated from observed spectra, there is some misprediction of L_{MM} calculated by Eq. (15), but it is relatively rare; this is shown by the low probabilities in Fig. 3 away from the wellpredicted, most commonly occurring L_{MM}.
To demonstrate the statistical character of Eq. (15), as well as its potential for probabilistic use (e.g., as input to probabilistic load calculations), Fig. 4 shows the probability density P(L_{MM}). As in Fig. 3, L_{MM} is again calculated from fits to 10 min spectra and also estimated by ${\mathit{\sigma}}_{u,\text{obs}}/(\mathrm{d}U/\mathrm{d}z)$, i.e., Eq. (15). Additionally Fig. 4 displays P(L_{MM}) for L_{MM} calculated through Eq. (11), i.e., ${c}_{m}{u}_{*}/(\mathrm{d}U/\mathrm{d}z)$; this is done both using the value of c_{m}=1.7 reported by Peña et al. (2010) as well as using the approximate mean of 2.3 found to be consistent with measurements and theory in Sects. 3.1 and 2.2 above. From Fig. 4 one sees that, for values of turbulent peak scale greater than the mode (∼20 m) up to roughly 150 m, there is a match between the distribution of the diagnosed L_{MM} and distributions of length scale estimated from the forms (15) based on σ_{u,obs} and (11) based on u_{*} with c_{m}=2.3; these are roughly equivalent for this case over relatively simple homogeneous terrain. It is found that the Peña et al. (2010) value of c_{m}=1.7 leads to overprediction of L_{MM} by a factor of 2 or more at scales smaller than 10 m and underprediction by 50 % or more at scales larger than 50 m. The u_{*}based form (11) using c_{m}=2.3 matches the spectrally fit diagnosed distribution P(L_{MM}) slightly better than the σ_{u}based form (15), with predicted peak (mode) values of L_{MM} being about 3–4 m smaller than the diagnosed peak L_{MM}.
For the homogeneous land case in Fig. 4 the probability density function of $\mathrm{2.3}{u}_{*,\text{obs}}/(\mathrm{d}U/\mathrm{d}z)$ matches P(L_{MM}) observed from the spectral fits to within 10 %, over the range 10 m ≲ L_{MM}≲75 m, and the probability density function of ${\mathit{\sigma}}_{u,\text{obs}}/(\mathrm{d}U/\mathrm{d}z)$ also matches within 10 % over the range 15 m ≲ L_{MM}≲50 m. This is consistent with the darkly colored 1 : 1 patch evident in Fig. 3 and also shows that Eq. (15) (and also Eq. 11 with c_{m}=2.3) is sufficient for probabilistic wind load simulations, for two reasons. First, the wellmatched range of scales corresponds to the most commonly found L_{MM}. Secondly, although scales smaller than ∼15 m are not rare (with an occurrence of roughly 1 in 6), they will have a diminishing effect on turbine loads. More specifically, L_{MM} is more than 70 % likely to fall in the 15–75 m range, i.e., P(15 m $<{L}_{\text{MM}}<\mathrm{75}$ m)>0.7, and L_{MM} has more than 86 % likelihood of occurrence between 0 and 75 m, for this homogeneous land case at z=80 m. The relatively common shorter scales correspond to weaker turbulent fluctuations (thus loads), because on average ${\mathit{\sigma}}_{u,\text{obs}}\propto {L}_{\text{MM}}^{\mathrm{2}/\mathrm{3}}$, as implied by Fig. 1 and Eqs. (9)–(15). Further, turbine loads are less influenced by fluctuations characterized by spatial scales significantly smaller than the blade lengths; thus the error in predicted probability for these shorter scales, and the slight underprediction of the most common L_{MM}, should not significantly influence probabilistic load calculations relying on sitespecific L_{MM} obtained via measurements and Eq. (15).
While Eq. (15) is useful to estimate L_{MM} and P(L_{MM}) as shown above, one expects Eq. (13) to perform better, as it does not rely on the approximation ${c}_{m}={\mathit{\sigma}}_{u}/{u}_{*}$. Indeed $\langle {c}_{m}{u}_{*}/{\mathit{\sigma}}_{u}\rangle $ is actually 1.11 (or 1.13 if considering winds only down to 7 m s^{−1}) due to ${\mathit{\sigma}}_{u}/{u}_{*}$ being slightly smaller and c_{m} slightly larger than 2.3; using these values in Eq. (13) gives estimates of L_{MM} closer to the spectrally diagnosed L_{MM}, and within 10 % of P(L_{MM}) over a range of L_{MM} from below 10 m to beyond 100 m. It should also be noted that ignoring speeds below 7 m s^{−1} can lead to slightly smaller L_{MM}, since these low wind speeds are more influenced by unstable conditions. Indeed for L_{MM}≳50 m, including the lower wind speeds causes both diagnosed and predicted L_{MM} to increase roughly 10 %; this is consistent with larger turbulent eddies being created under unstable conditions.
3.2.1 Estimating P(L_{MM}) in coastal/offshore conditions
To demonstrate the (probabilistic) use of Eq. (13) or Eq. (15) for L_{MM} in somewhat different conditions, we now consider flow from offshore, using data from the same mast and height as above (Høvsøre, z=80 m) but for wind directions between 240 and 300^{∘}. The mast is roughly 1.75 km east of the coastline and subsequently 1.65 km east of a 16–17 m high sand dune that lies 100 m inland, where both are locally oriented in the N–S direction (i.e., for the range of wind directions considered). The dune causes enhanced/accelerated transition of the flow from an offshore (water roughness) to an overland flow regime (Berg et al., 2015); this results in winds which reflect onshore and coastal conditions at low heights (below ∼40–80 m depending on stability) and offshore conditions at higher z.
Figure 5 displays the distribution P(L_{MM}) of spectralpeak (Mann model) length scales for coastal/offshore winds (from west ±30^{∘}), again using Eq. (15) to estimate L_{MM} along with L_{MM} diagnosed through spectral fits. For comparison the corresponding P(L_{MM}) for easterly winds from Fig. 4 is also included. Just as for the homogeneous land case shown in Fig. 4, one sees in Fig. 5 that, for inhomogeneous coastal conditions, again Eq. (15) gives P(L_{MM}) basically matching the spectrally fit observations for scales beyond ∼15 m; in this coastal regime the range of wellpredicted L_{MM} extends further, to ∼150 m. While one sees that the distribution of L_{MM} is a bit different for the (western) inhomogeneous coastal case than for the (eastern) homogeneous land case, the simple expression (15) functions similarly for both flow regimes, with the arguments presented in Sect. 3.2 again applying here.
The u_{*}based Eq. (11) also behaves similarly (not shown) as in the homogeneous land case of Fig. 4, i.e., with gross overpredictions at small scales and underpredictions at large scales. One difference between the coastal and land cases is that, for small L_{MM}, Eq. (15) overestimates the distribution P(L_{MM}) a bit more for the coastal regime than for the homogeneous land regime (L_{MM}<20 m); as explained above for the land case, an overprediction at the smallest scales is not expected to significantly impact load calculations, due to the relatively small length scales involved.
3.2.2 Estimation of P(L_{MM}) in more complex conditions
To further show the behavior of L_{MM} and the utility of Eq. (15) at a site with more complex conditions, we examine data from the inhomogeneous forested Danish National Test Centre for Large Wind Turbines site near Østerild in Denmark (Hansen et al., 2014). Here sonic anemometer data are available at heights of 10 and 44 m, with concurrent data from three lidars at $z=\mathit{\{}\mathrm{45},\phantom{\rule{0.125em}{0ex}}\mathrm{80},\phantom{\rule{0.125em}{0ex}}\mathrm{140},\phantom{\rule{0.125em}{0ex}}\mathrm{200},\phantom{\rule{0.125em}{0ex}}\mathrm{300}\mathit{\}}$ m. In this study we consider data from the site's “western lidar”^{10}, to measure winds that flow over the forest more than 70 % of the time, where the canopy height is 10–20 m (Hansen et al., 2014; Sogachev et al., 2017). The analysis here uses one year (May 2010–May 2011) of wind speeds U≥5 m s^{−1} from the lidar at 45 and 80 m heights along with the “fast” (20 Hz) data from the sonic anemometer at 44 m. The shear dU∕dz is measured across 45–80 m; the spectra and subsequent turbulence/Mannmodel parameters {${L}_{\text{MM}},\mathrm{\Gamma},\mathit{\epsilon}$}, as well as and measured quantities $\mathit{\{}{\mathit{\sigma}}_{u,\text{obs}},{u}_{*,\text{obs}}\mathit{\}}$, are obtained from the sonic anemometer. The measurements are significantly higher than twice the forest canopy height, and thus above the roughness sublayer (Garratt, 1980; Raupach et al., 1980) and amenable to similarity and mixinglength theory (Sogachev and Kelly, 2016) as well as Mannmodel use (Chougule et al., 2015).
Just as Fig. 4 showed for flow over homogeneous land at Høvsøre in Sect. 3, here Fig. 6 displays the probability density of turbulence (Mannmodel) length scale L_{MM} observed via spectral fits at z=44 m for Østerild, along with predictions based on both Eq. (11) via ${u}_{*,\text{obs}}$ and Eq. (15) via σ_{u,obs}.
As in the cases above (homogeneous land and inhomogeneous coastal), the new form (15) predicts the distribution rather well, particularly for scales ∼10–100 m – despite the shape of P(L_{MM}) being different due to the trees. For the forest case of Fig. 6 the σ_{u}based form captures both the peak (most likely L_{MM}) and magnitude of P(L_{MM}), while the u_{*}based form grossly underpredicts L_{MM}, more so than for the previous cases. The latter is likely due to ${u}_{*,\text{obs}}$ being predominantly affected by the canopy (via larger effective roughness) more so than σ_{u,obs}, which tends to be more characteristic of the entire ABL (Wyngaard, 2010). There is, however, a curious minor peak (with a probability ∼1 % as large as the main peak) around scales of $\sim \mathrm{300}\pm \mathrm{50}$ m in the lengthscale distribution obtained from spectral fits shown in Fig. 6; this is captured by neither the u_{*}based form (11) nor σ_{u}based formulations (13) and (15). Although this peak falls spectrally at small wavenumbers that are more difficult to capture when spectrally fitting the Mann model, it actually corresponds to the distance to the next upwind edge of the forest (orchard segment) in the predominant wind directions.
Towards concluding, we first revisit the motivation for (and thus context of) this work: (1) to “close” the Mann (1994) eddylifetime (τ_{M}) formulation as implemented in rapiddistortion theory – allowing relation between Mannmodel parameters (${L}_{\text{MM}},\mathit{\epsilon},\mathrm{\Gamma}$) and the shear (dU∕dz) taken to distort the modeled turbulence; (2) to connect the parameters of the Mann (1994) spectral turbulence and eddylifetime models with atmospheric statistics, both in theory and in practice; (3) to provide a formulation for the turbulence length scale L_{MM} in terms of quantities commonly measured in wind energy; and (4) to demonstrate that the “measurable” form developed for L_{MM} is robust and amenable to use in (probabilistic) wind turbine load calculations. These four motivating goals have basically been realized, as shown in the previous sections, and this work has a number of implications.
Implications and application
A previously suggested form (11) for L_{MM}, based on friction velocity u_{*} and (10 min) mean wind shear dU∕dz (Peña et al., 2010), was confirmed here to be sensitive to its proportionality constant c_{m}. But this constant can vary from site to site (and possibly with height), and the published value of c_{m}=1.7 (Peña et al., 2010) leads to significant error in prediction of L_{MM} for the different conditions (land and sea directions) at Høvsøre and at the forested site of Østerild. Finding c_{m} from sonic anemometer observations via L_{MM} from fits to spectra and friction velocity measurements, Eq. (11) may perform slightly better over uniform flat terrain compared to the σ_{u}based form (15) – but this can be considered a sitedependent fit in itself, as was the case when using a diagnosed value of c_{m}=2.3 for the homogeneous flat land sectors at Høvsøre. However, obtaining c_{m} is generally not possible in industrial practice; where it can be obtained, it relies on L_{MM} – which is the quantity desired – thus negating the purpose of Eq. (11). While u_{*} can also in principle be estimated from wind speeds taken at multiple heights by cup anemometers, this too is difficult in practice: one must account for stability, not to mention the need for measurements at multiple heights in the surface layer (or worse, the limited validity of similarity theory above the ASL). Furthermore, it is expected that c_{m} is a function of the (local) surface roughness, as demonstrated by the different results found over the forested Østerild site. Thus the form (15) is preferable, since it requires only the commonly measured quantities σ_{u} and dU∕dz. This simple form also gave good estimates of P(L_{MM}) in the forested case – without the need for tuning, whereas the u_{*}based form (11) requires a recalculation of its coefficient c_{m} for such cases.
Since Eq. (13) gave yet better performance than both its simplified form (15) and the u_{*}based relation (11), one might suggest its use. But Eq. (13) requires ${c}_{m}/({\mathit{\sigma}}_{u}/{u}_{*})$, where c_{m} is difficult to obtain, as discussed in the previous paragraph. However, although c_{m} might vary from site to site (or perhaps with height), it was found that the ratio ${c}_{m}/({\mathit{\sigma}}_{u}/{u}_{*})$ did not vary appreciably – consistent with the good performance of the simplified form (15), which assumes ${c}_{m}\approx {\mathit{\sigma}}_{u}/{u}_{*}$, across sites and regimes.
One interesting implication of the testing of assumptions then follows from the finding that $\langle {\mathit{\sigma}}_{u}/{u}_{*}{\rangle}_{\text{obs}}\approx \mathrm{2.3}$, consistent in the surface layer with Caughey et al. (1979). Examining the joint behavior of ${\mathit{\sigma}}_{u}/{u}_{*}$ and the stability parameter (inverse Obukhov length) L^{−1}, the sonic anemometer data available at multiple heights in this study show no correlation between these two quantities. The dimensionless profiles ${\mathit{\sigma}}_{u}^{\mathrm{2}}\left(z\right)/{u}_{*\mathrm{0}}^{\mathrm{2}}$ and ${u}_{*}^{\mathrm{2}}\left(z\right)/{u}_{*\mathrm{0}}^{\mathrm{2}}$ shown by Caughey et al. (1979) also imply
with the ratio converging to a constant above the surface layer (z≳0.1h, where the atmospheric boundarylayer depth h typically ranges from ∼200 m in stable conditions to 1 km or more in convective conditions). The flatterrain Høvsøre data in fact show the mean value $\langle {\mathit{\sigma}}_{u}/{u}_{*}{\rangle}_{\text{obs}}$ to be independent of z. If one knew how c_{m} varied with height (and stability), then one could also use Eqs. (20) and (13) from measurements at one height range, to estimate L_{MM} at higher z (for a given stability range). Over flat terrain, on average the peak spectral scale for streamwise fluctuations (λ_{u}) grows with z (Caughey et al., 1979; Peltier et al., 1996)^{11}. Therefore, if we take L_{MM}∝λ_{u}, then with Eq. (20) one expects the ratio ${c}_{m}/({\mathit{\sigma}}_{u}/{u}_{*})$ to increase with z as well. Thus from (13) the Mannmodel length scale L_{MM} will increase with height relative to the mixing length ${\mathrm{\ell}}_{*}\equiv {u}_{*}/(\mathrm{d}U/\mathrm{d}z)$, so at higher z one would expect the general form (13) to be yet more accurate than its approximate form (15); however, this is not likely for wind turbine rotor heights, except in very stable conditions (Kelly et al., 2014b; Liu and Liang, 2010). Unfortunately the sonicanemometer measurements available for this study did not include heights well beyond the surface layer, so such variation was difficult to detect.
It is also notable that Fig. 3 appears to imply the relative error (e.g., in %) in estimating L_{MM} with Eq. (15) grows for less common values of L_{MM}, particularly very large scales (and also at very small scales if including U<7 m s^{−1}). Thus Eq. (15) is recommended first for estimation of P(L_{MM}). However, the error at large scales is in part dependent on the limited (10 min) sample lengths and the fitting routine, as there are very few points to fit at the lowest frequencies. Use of 30 min samples can reduce such scatter, and modification of the fitting algorithms may also improve estimations of the larger scales.
Ongoing work includes windspeeddependent prediction of L_{MM}, particularly the conditional statistics P(L_{MM}U). Further concurrent work also entails systematic accounting for the rotor size (shear distance) relative to height (i.e., Δz∕z) within the distribution of length scales; following Kelly and Gryning (2010) and Kelly et al. (2014a) a semiempirical derivation of P(L_{MM}) including Δz∕z has been obtained but demands more data for validation and publication. Understanding of the latter facilitates “vertical extrapolation” of L_{MM} and measured turbulence and shear statistics, as well as accounting for the effect of rotor size or shear measurement span.

The eddy lifetime of Mann (1994), which is part of commonly used turbulence modeling for wind turbine design load cases (IEC 61400–1, Edition 3, 2005), leads to a relation for turbulence (spectralpeak) length scale L_{MM} of
$$}{L}_{\text{MM}}\simeq \phantom{\rule{0.33em}{0ex}}{\displaystyle \frac{{c}_{m}}{({\mathit{\sigma}}_{\text{u,obs}}/{u}_{*,\text{obs}})}}\phantom{\rule{0.33em}{0ex}}{\displaystyle \frac{{\mathit{\sigma}}_{u}}{\mathrm{d}U/\mathrm{d}z}},$$where c_{m} and ${\mathit{\sigma}}_{u,\text{obs}}/{u}_{*,\text{obs}}$ are essentially constants for a given height z, and ${c}_{m}/({\mathit{\sigma}}_{u,\text{obs}}/{u}_{*,\text{obs}})$ is found to fall between 1 and 1.11 for the three flow regimes analyzed.

Theory and measurements support the assumption that ${c}_{m}/({\mathit{\sigma}}_{\text{u,obs}}/{u}_{\text{*,obs}})\approx \mathrm{1}$, roughly constant for different atmospheric flow regimes; the turbulence length scale can consequently be approximated by
$$}{L}_{\text{MM}}\simeq \phantom{\rule{0.33em}{0ex}}{\displaystyle \frac{{\mathit{\sigma}}_{u}}{\mathrm{d}U/\mathrm{d}z}}.$$Thus typical 10 min mean cup anemometer measurements can be used to estimate L_{MM}.

L_{MM} is affected by atmospheric stability; this effect is contained within σ_{u} and dU∕dz.

In terms of the classic mixinglength form ${u}_{*}/\mathrm{d}U/\mathrm{d}z$, the turbulence length scale L_{MM} in the spectraltensor model is observed to be larger (by ca. 30–40 %) than previously reported by Peña et al. (2010).
The data are within an SQL database at DTU and are not publicly available.
The author declares that he has no conflict of interest.
The author thanks the reviewers for their time and effort towards
constructive criticism of the present article, and thanks are owed to Nikolay
Dimitrov for discussions around probabilistic loads. This work was partly
supported by the DTU Wind Energy internally funded
crosssectional project “Wind to Loads”.
Edited by: Horia Hangan
Reviewed by: two anonymous referees
Abramowitz, M. and Stegun, I. A.: Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; 9th printing, Dover, New York, 1972. a, b
Andre, J. C. and Lesieur, M.: Influence of helicity on the evolution of isotropic turbulence at high Reynolds number, J. Fluid Mech., 81, 187–207, https://doi.org/10.1017/S0022112077001979, 1977. a
Berg, J., Vasiljevic, N., Kelly, M. C., Lea, G., and Courtney, M.: Addressing Spatial Variability of SurfaceLayer Wind with LongRange WindScanners, J. Atmos. Ocean. Tech., 32, 518–527, https://doi.org/10.1175/JTECHD1400123.1, 2015. a
Caughey, S., Wyngaard, J. C., and Kaimal, J.: Turbulence in the evolving stable boundarylayer, J. Atmos. Sci., 36, 1041–1052, 1979. a, b, c, d, e
Chougule, A., Mann, J., Segalini, A., and Dellwik, E.: Spectral tensor parameters for wind turbine load modeling from forested and agricultural landscapes, Wind Energy, 18, 469–481, 2015. a
Chougule, A. S., Mann, J., Kelly, M., and Larsen, G. C.: Modeling Atmospheric Turbulence via Rapid Distortion Theory: Spectral Tensor of Velocity and Buoyancy, J. Atmos. Sci., 74, 949–974, https://doi.org/10.1175/JASD160215.1, 2017. a
ComteBellot, G. and Corrsin, S.: Simple Eulerian time correlation of full and narrowband velocity signals in gridgenerated, “isotropic” turbulence, J. Fluid Mech., 48, 273–337, 1971. 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, b, c, d
Dimitrov, N. K., Natarajan, A., and Kelly, M.: Model of Wind Shear Conditional on Turbulence and its impact on Wind Turbine Loads, Wind Energy, 18, 1917–1931, https://doi.org/10.1002/we.1797, 2015. a
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, 2017. a, b
Garratt, J. R.: Surface Influence Upon Vertical Profiles in the Atmospheric NearSurface Layer, Q. J. Roy. Meteor. Soc., 106, 803–819, https://doi.org/10.1002/qj.49710645011, 1980. a
Hansen, B. O., Courtney, M., and Mortensen, N. G.: Wind Resource Assessment – Østerild National Test Centre for Large Wind Turbines, Tech. Rep. DTU Wind Energy E0052(EN), Risø Lab/Campus, Danish Tech. Univ. (DTU), Roskilde, Denmark, 2014. a, b, c
IEC 61400–1, Edition 3: Wind turbine generator systems – Part 1: Safety requirements, International Electrotechnical Comission, Geneva, Switzerland, 2005. a, b, c
Kaimal, J. C., Wyngaard, J. C., Izumi, Y., and Coté, O. R.: Spectral characteristics of surfacelayer turbulence, Q. J. Roy. Meteor. Soc., 98, 563–589, 1972. a
Kelly, M. and Gryning, S.E.: LongTerm Mean Wind Profiles Based on Similarity Theory, Bound.Lay. Meteorol., 136, 377–390, 2010. a, b, c
Kelly, M., Larsen, G., Dimitrov, N. K., and Natarajan, A.: Probabilistic Meteorological Characterization for Turbine Loads, J. Phys. Conf. Ser., 524, 012076, https://doi.org/10.1088/17426596/524/1/012076, 2014a. a, b
Kelly, M., Troen, I., and Jørgensen, H. E.: Weibullk revisited: “tall” profiles and height variation of wind statistics, Bound.Lay. Meteorol., 152, 107–124, 2014b. a
Lesieur, M.: Turbulence in fluids: stochastic and numerical modelling, Springer, Dordrecht, NL, https://doi.org/10.1007/9789400905337, 1990. a, b
Liu, S. and Liang, X.Z.: Observed Diurnal Cycle Climatology of Planetary Boundary Layer Height, J. Climate, 23, 5790–5809, 2010. a
Mann, J.: The spatial structure of neutral atmospheric surfacelayer turbulence, J. Fluid Mech., 273, 141–168, 1994. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p
Mann, J.: The spectral velocity tensor in moderately complex terrain, J. Wind Eng., 88, 153–169, 2000. a, b
Mann, J., Ott, S., Jørgensen, B. H., and Frank, H. P.: WAsP Engineering 2000, Tech. Rep. RisøR1356(EN), Risø National Laboratory, Roskilde, Denmark, 2002. a
Mann, J., Astrup, P., Jensen, N., Landberg, L., and Jørgensen, H.: The meteorology of “the very large wind turbines”, in: Proc. of the 2004 European Wind Energy Conference, European Wind Energy Association (EWEA), London, 2005. a
Peltier, L. J., Wyngaard, J. C., Khanna, S., and Brasseur, J. G.: Spectra in the Unstable Surface Layer, J. Atmos. Sci., 53, 49–61, 1996. a, b, c
Peña, A., Gryning, S.E., and Mann, J.: On the lengthscale of the wind profile, Q. J. Roy. Meteor. Soc., 136, 2119–2131, 2010. a, b, c, d, e, f, g
Peña, A., Floors, R. R., Sathe, A., Gryning, S.E., Wagner, R., Courtney, M., Larsén, X. G., Hahmann, A. N., and Hasager, C. B.: Ten Years of BoundaryLayer and WindPower Meteorology at Høvsøre, Denmark, Bound.Lay. Meteorol., 158, 1–26, https://doi.org/10.1007/s1054601500798, 2016. a
Pope, S. B.: Turbulent Flows, Cambridge University Press, Cambridge, UK, 2000. a, b, c
Raupach, M., Thom, A., and Edwards, I.: A WindTunnel Study of TurbulentFlow Close to Regularly Arrayed Rough Surfaces, Bound.Lay. Meteorol., 18, 373–397, https://doi.org/10.1007/BF00119495, 1980. 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, 2013. a, b, c
Savill, A. M.: Recent developments in rapiddistortion theory, Annu. Rev. Fluid Mech., 19, 531–575, 1987. a
Sogachev, A. and Kelly, M.: On Displacement Height, from Classical to Practical Formulation: Stress, Turbulent Transport and Vorticity Considerations, Bound.Lay. Meteorol., 158, 361–381, https://doi.org/10.1007/s105460150093x, 2016. a
Sogachev, A., Cavar, D., Kelly, M., and Bechmann, A.: Effective roughness and displacement height over forested areas, via reduceddimension CFD, Tech. Rep. DTU Wind Energy E0161(EN), Risø Lab/Campus, Danish Technical University (DTU), Roskilde, Denmark, 2017. a
von Kármán, T.: Progress in the Statistical Theory of Turbulence, P. Natl. Acad. Sci. USA, 34, 530–539, https://doi.org/10.1073/pnas.34.11.530, 1948. a, b
Wyngaard, J. C.: Turbulence in the Atmosphere, Cambridge University Press, Cambridge, UK, 2010. a, b
To a lesser extent, some sensitivity to the Mannmodel anisotropy parameter Γ has also been found.
The Mann (1994) expression is also equivalent (or at least proportional) to the “convection time” of ComteBellot and Corrsin (1971).
The reciprocal of eddydamping rate, $\left[{k}^{\mathrm{3}}E\right(k){]}^{\mathrm{1}/\mathrm{2}}$, is equal in the inertial range to Eq. (1) with $p=\mathrm{1}$ since $E\left(\mathit{\kappa}\right)\to {\mathit{\kappa}}^{\mathrm{5}/\mathrm{3}}$ there. This expression is also similar to the “rotation time” or “strain time” given by ComteBellot and Corrsin (1971), but it should be noted that such expressions integrate from 0 to k, i.e., over eddies larger than 1∕k.
The hypergeometric function ${}_{\mathrm{2}}{F}_{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}\left[\frac{\mathrm{1}}{\mathrm{3}},\frac{\mathrm{17}}{\mathrm{6}};\frac{\mathrm{4}}{\mathrm{3}};(k{L}_{\text{MM}}{)}^{\mathrm{2}}\right]$ approaches 1 for kL_{MM}≫1 (the inertial range) and simplifies to a_{HG}(kL_{MM})^{2∕3} for kL_{MM}≪1, where ${a}_{\mathrm{HG}}\equiv (\mathrm{3}\sqrt{\mathit{\pi}}/\mathrm{4}){f}_{\mathrm{\Gamma}}(\mathrm{4}/\mathrm{3})/{f}_{\mathrm{\Gamma}}(\mathrm{17}/\mathrm{6})\simeq \mathrm{0.69}$ and f_{Γ}(x) is the Euler gamma function.
The peak of the von Kármán isotropic TKE spectrum E_{vK}(k) occurs at $k{L}_{\text{MM}}=\sqrt{\mathrm{12}/\mathrm{5}}$, i.e., ${L}_{\text{MM}}\simeq \mathrm{1.55}/{k}_{\text{peak}}$.
Note $\sqrt{\mathrm{2}/\mathrm{3}}\simeq \mathrm{0.82}$ and $\mathrm{3.07}={a}_{\mathrm{HG}}^{\mathrm{3}}$; cf. footnote 4. In Eq. (6), αε^{2∕3} is kept together for comparison with Eq. (4) and because αε^{2∕3} is commonly used as an input to the spectraltensor model instead of ε (IEC 61400–1, Edition 3, 2005; Mann et al., 2002).
Assuming a constant mean shear dU∕dz, the spectraltensor model solves Fouriertransformed versions of rapiddistortion equations for streamwise normal stress 〈u_{1}u_{1}〉 and shear stress 〈u_{1}u_{3}〉; multiplying these by dU∕dz one obtains the corresponding production rates: ${P}_{\mathrm{11}}=\mathrm{2}\langle {u}_{\mathrm{1}}{u}_{\mathrm{3}}\rangle \mathrm{d}U/\mathrm{d}z$ and ${P}_{\mathrm{13}}=\langle {u}_{\mathrm{3}}{u}_{\mathrm{3}}\rangle \mathrm{d}U/\mathrm{d}z$ (Pope, 2000).
The sonic anemometers actually give a temperature very close to the virtual temperature, i.e., the temperature including buoyant effects of water vapor.
The spectral fits were done using spectraltensor model output over the parameter ranges of $\mathrm{5}<{L}_{\text{MM}}<\mathrm{500}$ m and $\mathrm{0}\le \mathrm{\Gamma}\le \mathrm{5}$. Some spectra were poorly fitted; since these occurred when Γ=5, cases with Γ>4.95 were excluded from the analysis here. As justification, I note that only a small fraction of the cases (<10 %) had such Γ, and we only consider wellfit spectra for reliable comparison of parameters.
The “western lidar” at Østerild is located ∼1 km west of the northernmost turbines but less than 100 m east of a forest patch and 5–20 km from the North Sea coastline in the prevailing (W–NW) wind directions (Hansen et al., 2014).
The peak length scale also grows with boundarylayer depth h in convective conditions and thus with increasingly negative inverse Obukhov length L^{−1} (Peltier et al., 1996). But over all stability conditions, which are dominated by neutral conditions (Kelly and Gryning, 2010), and over an expected distribution of h at a given site, the basic growth of λ_{u} with z is consistent with Peltier et al. (1996) reporting λ_{u}∝z for neutral conditions.