On Modelling Wind-Farm Wake Turbulence Autospectra and Coherence from a Database

This study addresses the feasibility of modeling wind-farm wake-turbulence autospectra and coherences from a database: flow velocity points from experimental and computational fluid dynamics (CFD) investigations. Specifically, it first applies an earlier-exercised framework to construct the autospectral models from a database and then it adopts a recently proposed framework to construct the coherence models from a database. While this proposed framework has not been tested against a database, the methodology has been completely formulated with a theoretical basis. These models of autospectrum and coherence are interpretive, and in closed form. Both frameworks basically involve the perturbation series expansion of the autospectra and coherences. The framework for modeling autospectra is tested against a demanding database of wake turbulence inside a wind farm over a complex terrain from a full-scale test. The suitability of these autospectral models for simulation through white-noise driven filters is also demonstrated. Finally, coherence models are generated for assumed values of the perturbation series constants, and these coherence models are used to demonstrate how the coherence models of homogeneous isotropic turbulence deviate from the coherence models of non-homogeneous non-isotropic turbulence such as wind-farm wake turbulence. This feasibility of extracting both the one-point statistics of autospectral models and the two-point statistics of coherence models from a database represents a research avenue that is new and promising in the treatment of wind-farm wake turbulence. This paper also demonstrates the feasibility of fruitfully exploiting the wake treatment methods developed in other fields.


Introduction
During the past thirty years, wake turbulence and its effects on wind turbines and wind farms have been extensively investigated, primarily analytically and to some extent experimentally. The extensive literature up to 2010 has been well covered in the widely used text of Manwell et al. [1]. As for the extensive analytical investigations since, suffice it to mention, as representative samples, Keck et al. [2] for wake-turbulence modeling from the low-fidelity CFD treatment of the Navier-Stokes (NS) equations, and Carrion et al. [3] for wake-turbulence modeling from the high-fidelity treatment of NS equations. The work of Carrion et al. [3] for example, includes a concise account of the state of the art of modeling wind-farm wake turbulence.
An overview of these investigations [1][2][3] is included here; although extremely brief, this should help appreciate how the present work serves as a desirable adjunct of experimental and high-fidelity CFD based investigations, and why it represents a new and promising avenue of wake-turbulence modeling. Wake-turbulence modeling falls into three categories:

1.
Semi-empirical models: These models are based on conservation of momentum and such simplified assumptions; and they typically contain an empirical constant. For example, the widely used model due to Katic et al. [1] belongs to this category; therein, the empirical constant is referred as wake decay constant. A recent work of Ge M. et al. [4] merits mention: it gives a thorough account of the key features of this class of models and their continual evolution as well as their utility base in the treatment of wind energy applications such as design of wind-farm layout. 2.
Low-fidelity CFD solutions: Several numerical schemes have been continually proposed (e.g., dynamic wake meandering model [2]); these are based on a wide range of physics-based approximations to NS equations. As more and more experimental databases become available, they are being continually updated and they merit further validation. 3.
High-fidelity CFD solutions: Despite the severe CPU-hour constraint, these high-fidelity approaches are indispensable in generating a database that serves as a reference point and in supplementing databases from experiments.
While major strides have been made in generating databases and in providing a much-improved understanding of wake turbulence, the current capability for modeling the one-point statistics of autospectrum, much more so, for modeling the two-point statistics of coherence, merits significant improvement. In fact, empirical exponential coherence functions are still being used [1]. Accordingly, the present study explores the feasibility of extracting the one-point statistics of autospectrum and the two-point statistics of coherence from a database.
The autospectral model extraction from a database of the present study is based on the framework due to Schau, Gaonkar and Polsky [5]. This framework guarantees that the extracted model and the autospectral data points have the same mean square value (a measure of turbulence energy), time scale and the Kolmogorov −5/3 spectral decay. For completeness, an earlier study by Gaonkar [6] should be mentioned as well; therein, the autospectral model extraction from a database is based on the framework of reference [7], in which the −5/3 spectral law is bypassed. Now it is expedient to address the development of a framework for extracting the two-point statistics models from a database. This can be approached either through cross-spectrum, which is a complex quantity involving the magnitude and phase or through coherence, which, as a spectral correlation coefficient, is a real quantity. The first approach generally leads to modeling the magnitude and not the phase, as was the case in Ref. [7]. The second approach through coherence is relatively more convenient and provides a means of capturing the two-point statistics from a database completely. The recent study due to Krishnan and Gaonkar [8] follows this second approach; although not tested against a database, the framework is formulated with a mathematical basis and the present study adopts this framework [8].
By design, these autospectral and coherence models are in closed form and they have a simple analytical structure to facilitate interrogation and interpretation of voluminous data points on autospectra and coherences. And they lend themselves well to routine use as a predictive tool. Compared to a description through such voluminous, numerically generated, autospectral and coherence data points, they describe wake turbulence analytically with better transparency and bring better understanding. Thus, these interpretive models broaden the scope and utility base of the database that invariably involves enormous resources. While the extracted models are database-specific (thus they are not predictive by themselves), the framework can be applied to any database and the model extraction is a routine exercise.
In the treatment of coherence for homogeneous isotropic turbulence for which the frozen turbulence hypothesis is applicable (HIT), the present study is motivated by and built on the earlier studies of Burton et al. [9], Houbolt and Sen [10], Frost et al. [11] and Irwin [12]. This treatment of coherence for HIT is found to show differences among these studies [9][10][11][12], and the present study, after an in-depth examination, follows Frost et al. [10] for cross-spectra and Irwin [12] for coherences. Given this background, the present study seems to provide a unified account of coherence for HIT in the treatment of wake turbulence.
To sum up: These interpretive models complement the experimental and CFD-based investigations as surrogate analytical models for both the one-point statistics of autospectrum and the two-point statistics of coherence. Moreover, this paper also demonstrates the feasibility of fruitfully exploiting the methodologies from other fields to the treatment of wind-turbine wake turbulence. And these methodologies offer promise towards providing a foothold on a formidably complex flow field inside a windfarm for engineering analysis.

Basic of Modeling
A comparison of the measured autospectra of ambient atmospheric boundary layer turbulence (ABL) and wake turbulence shows that the ABL autospectrum has gone through changes in energy distribution with respect to frequency. Figure 1 [13] should help bring a better understanding of this comparison; specifically, it shows measured dimensionless longitudinal autospectrum fŠ uu /σ 2 u versus dimensionless frequency f z/U, where z is the mast height and U is the mean wind speed. These autospectra were experimentally generated at the same location in a wind farm over a complex terrain. Figure 1a refers to ABL with a turbulence intensity of 0.103, when the turbines were under stand-still conditions. Furthermore, Figure 1b refers to wake turbulence with a turbulence intensity of 0.204, when the turbines were fully operational. This change in the shape of the wake turbulence autospectrum cannot be realized through a linear superposition of a series of independently occurring changes at different frequencies on the ABL autospectrum. Thus, the autospectral morphing must be due to a nonlinear transformation of ABL. Stated otherwise, wind-farm wake turbulence could be idealized as nonlinearly transformed ABL and in turn, an earlier-developed mathematical framework for autospectral modeling of airwake-downwash turbulence could be adapted to modeling wind-farm wake turbulence as well [5,14]. (Airwake-downwash turbulence refers to the coupled flow-field of ship's airwake shed from the superstructure and the helicopter downwash. Therein [5,14], the mathematical framework "posits" that airwake-downwash turbulence is nonlinearly transformed ABL. Regarding the coherence, the framework of Ref. 8 is adopted with the same justification that is used for the autospectrum.)

Basic of Modeling
A comparison of the measured autospectra of ambient atmospheric boundary layer turbulence (ABL) and wake turbulence shows that the ABL autospectrum has gone through changes in energy distribution with respect to frequency. Figure 1 [13] should help bring a better understanding of this comparison; specifically, it shows measured dimensionless longitudinal autospectrum / versus dimensionless frequency / , where is the mast height and is the mean wind speed. These autospectra were experimentally generated at the same location in a wind farm over a complex terrain. Figure 1a refers to ABL with a turbulence intensity of 0.103, when the turbines were under stand-still conditions. Furthermore, Figure 1b refers to wake turbulence with a turbulence intensity of 0.204, when the turbines were fully operational. This change in the shape of the wake turbulence autospectrum cannot be realized through a linear superposition of a series of independently occurring changes at different frequencies on the ABL autospectrum. Thus, the autospectral morphing must be due to a nonlinear transformation of ABL. Stated otherwise, wind-farm wake turbulence could be idealized as nonlinearly transformed ABL and in turn, an earlier-developed mathematical framework for autospectral modeling of airwake-downwash turbulence could be adapted to modeling wind-farm wake turbulence as well [5,14]. (Airwake-downwash turbulence refers to the coupled flow-field of ship's airwake shed from the superstructure and the helicopter downwash. Therein [5,14], the mathematical framework "posits" that airwake-downwash turbulence is nonlinearly transformed ABL. Regarding the coherence, the framework of Ref. 8 is adopted with the same justification that is used for the autospectrum.)

Methodology of AutoSpectra
The Lateral component v(t) is selected for providing details of the autospectral model extraction methodology [5]. The methodology remains the same for the vertical component with minor changes in the parameters used in the constraint equations. However, for the Longitudinal component, additional changes in the constraint equations are also required [14]. One-sided autospectrum is used throughout.
Statistical independence of velocity components is assumed [9]. The mean square value, time scale and autospectral asymptotic limit law are different for each component of non-homogeneous turbulence according to Kolmogorov's −5/3. The framework combines four elements: (1) A mathematical formulation based on a perturbation series expansion of the autocorrelation/autospectrum functions; (2) Extraction of time scale and autospectral asymptotic limit from the database; (3) Development of constraint equations in closed form to ensure that the developed model satisfies the requirements related to normalization, time scale, and autospectral asymptotic limit; and (4) Evaluation of the constants in the series expansion subject to satisfying the constraint equations and fitting a curve on a set of selected autospectral data points in a least squares sense.

Lateral Wake Turbulence
The perturbation series for the autocorrelation of lateral wake turbulence velocity v(t) can be expressed as in Equation (1).
The calculated autocorrelation as well as series expansion autocorrelation follow the properties of normalized autocorrelations, that is, R vv (0) = R vv (0) = 1. The Fourier transform of Equation (1) gives the series expansion for the autospectrum S vv ( f ): where S vvn ( f ) is the Fourier transform of R n vv (τ).
The autospectrum is typically normalized with respect to dimensional time scale T v , which is traditionally defined as v , the mean square value, Equations (4) and (5) typify the normalization: According to Kolmogorov's spectral law, the autospectrum model should decay as given in Equation (6), where A v is a scaling parameter determined from the data, and f h f represents high frequencies.
The basis function R vv (τ) on the right-hand side of Equation (1) is the von Karman lateral correlation function as given in Equation (7).
The scaling parameter α v in Equation (7) ensures that the relation in Equation (4) is satisfied.

Database
The database is generated from a full-scale experimental study of wake turbulence by Mofiadakis et al. [13]. Specifically, data is collected at several points along a complex windy terrain at an altitude of 320-330 m with seven Vestas (27-225 kW) installed in a row. The thirty cases of autospectra that were generated from the database represent free stream to fully wake affected conditions. The autospectral decay was found to be in the range of −1.36 to −1.75 for free stream conditions (all wind turbines at stand-still condition). Moreover, for notational simplicity, β iu , β iv and β iw are simply referred to as β coefficients in this section.
The database typically comprises the temporal flow velocity points. In this case, however, the database [13] has already been transformed from the temporal to the frequency domain and normalized to dimensionless form f S ii ( f )/σ 2 i . Moreover, the temporal autocorrelation is not available, and the dimensionless autospectra are presented on a log-log scale and in turn S ii ( f = 0) is not available. This limitation is overcome by assuming that S ii (0) = S ii ( f ) as f approaches zero. It is emphasized that the framework is designed to develop autospectral models from a database; thus the lack of a temporal database ceases to be a major issue. Having extracted σ 2 i from Equation (5), key information to be extracted from the database is T i , the time scale. As for T i , it is calculated using S ii ( f ) at the lowest frequency in a log-log plot; see Equation (4). Finally, the autospectral decay constant A i is graphically evaluated from Equation (6).
Having thus generated time scale T i and autospectral asymptotic limit A i , the numerical scheme now focuses on computing the series expansion β coefficients. It is emphasized that these β coefficients determine the scaling parameter α i ; see Equation (9). As an iterative procedure, the scheme involves selecting the β coefficients, beginning with the von Karman model (e.g., β 1v = 1 for the lateral component) and strictly enforcing the constraint as typified by Equation (8). For completeness, the gist of the iterative procedure is included; for details see [14].
The numerical scheme minimizes the sum of two errors in a least squares sense: model's deviation from the autospectral data and the A i constraint error. That is, a selected set of β coefficients gives a model with a value of A i ; stated otherwise, these β coefficients carry a least squares error with respect to the measured autospectral data and an error with respect to the graphically measured A i value. The resulting A i error is expected to be within acceptable limits for wind turbine applications (<<10%). This error is perhaps acceptable, after all, A i is not rigorously defined with respect to the data sets, nor is there a standard method of determining when a computed autospectrum has reached its point of asymptotic decay. This lack of precision also means A i will vary somewhat from user to user. The A i constraint typified by Equation (10) also merits one final comment. For some isolated cases of data sets, the high frequency limit does not exhibit accurately the −5/3 spectral decay [13]. For these cases, it is sensible to exclude this constraint. The next section elaborates this scenario.
Computationally, the above numerical scheme is found to be inexpensive. For example, a gradient based search algorithm (MATLAB SQP) starting with the von Karman model does not require large iteration counts for a converged model. The reason is that the computational cost, now, depends only on the length of the discrete data array, as an error between the model and this data array, and this error has to be calculated at each iteration of the search algorithm. Given the thoughtful selection of search parameters such as step size in β-space and convergence criteria, this numerical scheme should prove inexpensive computationally.

Result of Autospectra
For illustration, just one example of the vertical component w(t) is selected. Modeling is presented based on both first-order (two-term series) and second-order (three-term series) correction. And in each case, modeling covers two approaches. In the first approach, the Kolmogorov −5/3 law is not enforced; that is, by satisfying only the first two constraints, typified by Equations (8) and (9). In the second approach, all three constraints are satisfied; that is, in addition to satisfying these two constraints, the model also satisfies the Kolmogorov −5/3 law (see Equation (10)) by a minimized error. This enforcement is identified in the respective figures by "A error = xx%", where xx% indicates the percentage error in satisfying the −5/3 law. Typically, an "A error" of less than 10% is considered satisfactory. Throughout, the dimensionless autospectrum f S ii ( f )/σ 2 i is presented against frequency f (Hz). Furthermore, in each figure, the corresponding von Karman model (e.g., β 1v = 1 for the lateral) is also included; this helps assess how far the developed model is an improvement over the von Karman, a widely used model for the free-stream case [13]. For additional results, see Schau [14].
For the vertical component in Figure 2, the corresponding first-order-correction (a two-term series) models represent appreciable improvement over the von Karman, particularly for f > 10 −1 Hz. Overall, modeling still merits further improvements for f > 10 −1 Hz. For the vertical component (Figure 2b), the enforcement of the Kolmogorov −5/3 law in a least squares sense involves "A error = 23.56%", well above the stipulated "A error" of 10%. These two features, the feasibility of improving the correlation and reducing the "A error", is explored in the next figure based on the second-order-correction (a three-term series). The results of Figure Figure 3b) shows that the three-term series model improves the correlation throughout, particularly for f > 10 −1 Hz. Second, the "A error", which is 23.56% for the two-term series model comes down to 7.64%. Thus, this comparison shows that the three-term series model is a noteworthy improvement over the two-term series model, without or with the enforcement of the −5/3 law. To sum up: modeling based on first-order correction (a two-term series) is generally adequate, and further improvement in correlation and further reduction in "A error = x" can be achieved through modeling based on second-order correction (a three-term series). Figure 4 shows how the autospectrum from the white-noise-driven filter for the developed vertical model compares with the one from the database and the developed model itself (specifically, Figure 4 refers to Figure 3a). As seen from this Figure 4, the developed model and simulation are almost indistinguishable. (The filter represents a single-input, single-output system driven by white noise; the design is routine and thus the details are omitted [14]).

Methodology of Coherence
As done for autospectra, for coherences also, a mathematical framework is developed for extracting interpretive coherence models from a database of flow velocity points from experimental and CFD investigations. Here as well, each velocity component is considered statistically independent of the other two. For each velocity component, the framework begins with a perturbation series expansion of the coherence; therein, the basis function or the first term of the series is represented by the corresponding coherence for HIT. The perturbation coefficients are evaluated by satisfying the theoretical constraints and fitting a curve on a set of numerically generated coherence points from a database.
In the literature, the development of the cross-spectra and coherences for the longitudinal, vertical and lateral components is scattered and piecemeal; what is more, the expressions for these cross-spectra and coherences show difference among these studies. Accordingly, this section first presents the cross-spectrum, after all, coherence is cross-spectrum that is normalized by the corresponding autospectrum (details to follow). Then it presents the coherences and finally a perturbation theory scheme for the wind-farm wake-turbulence coherence.

Construction of the Vertical Cross Spectrum
For illustration, vertical turbulence w(t) is considered under headwind conditions. Given , the mean wind velocity, , the elapsed time ( − ) and the correlation distance = , the von Karman correlation function ( ) for vertical turbulence ( ) is given by Equation (11) [15]: where = /1.339 , is the scale length and is the modified Bessel function of the second kind. Now consider the cross-correlation between vertical turbulence ( ) at Point 1 and ( ) at Point 2, where these two points are separated by the across-wind distance , as typified by Figure 5. For this scenario, Figure 5 shows that the correlation distance changes to the expression given in Equation

Methodology of Coherence
As done for autospectra, for coherences also, a mathematical framework is developed for extracting interpretive coherence models from a database of flow velocity points from experimental and CFD investigations. Here as well, each velocity component is considered statistically independent of the other two. For each velocity component, the framework begins with a perturbation series expansion of the coherence; therein, the basis function or the first term of the series is represented by the corresponding coherence for HIT. The perturbation coefficients are evaluated by satisfying the theoretical constraints and fitting a curve on a set of numerically generated coherence points from a database.
In the literature, the development of the cross-spectra and coherences for the longitudinal, vertical and lateral components is scattered and piecemeal; what is more, the expressions for these cross-spectra and coherences show difference among these studies. Accordingly, this section first presents the cross-spectrum, after all, coherence is cross-spectrum that is normalized by the corresponding autospectrum (details to follow). Then it presents the coherences and finally a perturbation theory scheme for the wind-farm wake-turbulence coherence.

Construction of the Vertical Cross Spectrum
For illustration, vertical turbulence w(t) is considered under headwind conditions. Given V, the mean wind velocity, τ, the elapsed time (t 2 − t 1 ) and the correlation distance x = Vτ, the von Karman correlation function R ww (x) for vertical turbulence w(t) is given by Equation (11) [15]: where u = x/1.339L, L is the scale length and K n is the modified Bessel function of the second kind. Now consider the cross-correlation between vertical turbulence w 1 (t) at Point 1 and w 2 (t) at Point 2, where these two points are separated by the across-wind distance S, as typified by Figure 5. For this scenario, Figure 5 shows that the correlation distance changes to the expression given in Equation (12a): where σ = S/L. Now, the cross-correlation R w 1 w 2 (x) can be expressed as in Equation (12b) [10]: Given this background, it is emphasized that in the present study, the expressions of crossspectra as typified by Equations (12c), (17a) and (18a) for the vertical, longitudinal and lateral components agree with those of Frost et al. [11]. As for coherence, the corresponding expressions given by Equations (16), (17b) and (18b) agree with those of Irwin [12]. Figures 6-8 Figure 9, which is a recasting of Figure 7 for a much expanded vertical scaling (1 to 10 in Figure 9 in comparison to 1 to 10 in Figure 7). The crosses (*) in Figure 9 indicate the termination of the curve to avoid generating negative coherence values. Given the state of the art and one's initiation into cross-spectra and coherence, it is difficult to come up with a basis for these negative values of cross-spectrum and coherence for HIT; a resolution of this difficulty would require further research [11]. The Fourier transform of R w 1 w 2 (τ) is the cross-spectrum S w 1 w 2 (ν) [10]: where, In Equation (12d), ν represents the dimensionless frequency ν = ωL/V and σ = S/L, the dimensionless distance.

Coherence for HIT
Coherence is also referred to as spectral correlation coefficient in that it quantifies the normalized cross-correlation between the turbulence velocities at two points as a function of frequency. For illustration, consider the vertical turbulence velocities at two points which are separated by a distance S, as typified by Figure 5. By definition, coherence is given by: where S w 1 w 2 (ν) is the cross-spectrum between vertical turbulence w 1 (t) at Point 1 and w 2 (t) at Point 2, and similarly S w 1 w 1 (ν) and S w 2 w 2 (ν) are the corresponding autospectra of w 1 (t) and w 2 (t). For HIT, cross-spectrum is real and S w 1 w 1 (ν) ≈ S w 2 w 2 (ν). Therefore, coherence from Equation (13) simplifies to Equation (14).
As seen from Equation (14), C w 1 w 2 (σ, ν) is a ratio of the cross-spectrum from Equation (12c) and the autospectrum from Equation (15). It is expedient to reiterate that this autospectrum is due to von Karman [15] and that the cross-spectrum is due to Houbolt and Sen [10], as an extended version of the von Karman spectral equations that accounts for the cross-correlation between vertical turbulence velocities at two points; also see Figure 5. After some algebra, Equation (14) simplifies to Equation (16) [8].
Given this background, it is emphasized that in the present study, the expressions of cross-spectra as typified by Equations (12c), (17a) and (18a) for the vertical, longitudinal and lateral components agree with those of Frost et al. [11]. As for coherence, the corresponding expressions given by Equations (16), (17b) and (18b) agree with those of Irwin [12]. Figures 6-8, respectively, show vertical, longitudinal and lateral coherence between Points 1 and 2 as a function of dimensionless frequency ν = ωL/V for σ = S/L = 0, 0.1, 0.2, . . . 1. For σ = 0, Point 2 merges into Point 1 and in turn the cross-spectra become the respective autospectra and thus the coherence represents the perfect coherence. For example, as seen from Figure 6 for the vertical coherence, C w 1 w 2 (σ, ν) → 1 . Similarly, as seen from Figures 7 and 8, C u 1 u 2 (σ, ν) → 1 and C v 1 v 2 (σ, ν) → 1 . Exactly the opposite happens with increasing σ = S/L. That is, with increasing σ, the correlation between these two points decreases and so does the corresponding coherence. For example, as seen from Figures 6-8 Moreover, as seen from these figures, the coherence decreases rapidly for ν > 1 or so.
The longitudinal cross-spectrum S u 1 u 2 (ν) and coherence C u 1 u 2 (σ, ν) are typified by Equations (17a) and (17b), respectively, and Figure 7 shows coherence C u 1 u 2 (σ, ν) as a function of dimensionless frequency ν = ωL/V; all of this merits revisiting. The reason is that S u 1 u 2 (ν) and in turn the corresponding coherence can become negative at high frequencies. As seen from Equations (17a) and (17b), respectively, S u 1 u 2 (ν) and C u 1 u 2 (σ, ν) can take on negative values for K 5/6 (z) ≤ z/2K 1/6 (z). See Figure 9, which is a recasting of Figure 7 for a much expanded vertical scaling (1 to 10 −6 in Figure 9 in comparison to 1 to 10 −2 in Figure 7). The crosses (*) in Figure 9 indicate the termination of the curve to avoid generating negative coherence values. Given the state of the art and one's initiation into cross-spectra and coherence, it is difficult to come up with a basis for these negative values of cross-spectrum and coherence for HIT; a resolution of this difficulty would require further research [11].

Coherence Modeling for Wind-Farm Wake Turbulence
Wind farm wake turbulence deviates from HIT. Accordingly, the framework for coherence modeling from a database accounts for this deviation based on perturbation theory. Here as well, the framework assumes the same topology that was assumed in the development of the basis functions; for illustration vertical coherence ( , ) is selected.
The basis function or the first term of the series is given by Equation (16). Since ( , ) = 1 for = 0, Equation (17) is subject to the constraint: The second condition that ( , ) = 0 for = ∞ is automatically satisfied since ( , ) = 0 for = ∞ . The coefficients in the series are evaluated by satisfying the theoretical constraint of Equation (20) and fitting a curve on a set of selected numerically generated coherence points in a least squares sense.
For illustrations, longitudinal coherence of wake turbulence ( , ) is selected with a twoterm perturbation series (also see Equation (19)): Specifically, consider = 0.7 and = 0.3 (also see constraint Equation (20)). Descriptively stated, this case represents wake turbulence, which deviates weakly from HIT. It is plausible that this case belongs to wake turbulence at locations that are downwind of the first two rows or so. Therein, wake turbulence is expected to deviate only weakly from HIT as depicted in Figure 10.

Coherence Modeling for Wind-Farm Wake Turbulence
Wind farm wake turbulence deviates from HIT. Accordingly, the framework for coherence modeling from a database accounts for this deviation based on perturbation theory. Here as well, the framework assumes the same topology that was assumed in the development of the basis functions; for illustration vertical coherence C w 1 w 2 (σ, ν) is selected.
Let C w 1 w 2 (σ, ν) represent the vertical coherence of wake turbulence. The framework begins with a perturbation series expansion of C w 1 w 2 (σ, ν) (essentially the same procedure applies to the other two components): The basis function or the first term of the series is given by Equation (16). Since C w 1 w 2 (σ, ν) = 1 for σ = 0, Equation (17) is subject to the constraint: The second condition that C w 1 w 2 (σ, ν) = 0 for σ = ∞ is automatically satisfied since C w 1 w 2 (σ, ν) = 0 for σ = ∞. The coefficients in the series C iw are evaluated by satisfying the theoretical constraint of Equation (20) and fitting a curve on a set of selected numerically generated coherence points in a least squares sense.
For illustrations, longitudinal coherence of wake turbulence C u 1 u 2 (σ, ν) is selected with a twoterm perturbation series (also see Equation (19)): Specifically, consider C 1u = 0.7 and C 2u = 0.3 (also see constraint Equation (20)). Descriptively stated, this case represents wake turbulence, which deviates weakly from HIT. It is plausible that this case belongs to wake turbulence at locations that are downwind of the first two rows or so. Therein, wake turbulence is expected to deviate only weakly from HIT as depicted in Figure 10.

Conclusions
This study has shown that an earlier-exercised mathematical framework lends itself well to extracting interpretive autospectral models of wind farm wake turbulence from a database. While these models are database specific, the framework can be applied to any database and the model construction is straightforward. As to the two-point statistics of wake-turbulence, this study first presents a unified account of cross-spectrum and coherence for HIT; this account is of considerable utility in that, in the literature, the expressions of cross-spectrum and coherence show differences from one study to the other. Given these expressions of coherences, this study then builds a framework for extracting wake turbulence coherence models from a database. The frameworks for autospectra and coherences do not follow classical perturbation theory approach of a solution for a linearized problem along with successively added corrections. Both frameworks represent a practical combination of a series expansion, exploitation of a database, and theoretical constraints in closed form.
This study also leads to following specific findings: 1. Generally, no more than a three-term series (second-order correction) is necessary to develop an autospectral model; in most cases, a two-term series (first-order correction) is found to be adequate for wind engineering applications. 2. The addition of a third term to the series has significant power in reducing the " " between the model and the data. Recall that the " " refers to minimizing the sum of the errors in a least squares sense: the model deviation from the autospectral data points and from the measured high frequency autospectral decay level (when applicable). 3. These developed models lend themselves well to design of filters driven by white noise; that is, the filter design is as routine as the currently used procedure for the von Kármán models. 4. While this framework to constructing the coherence models from a database has not been tested against a database, it has been formulated from first principles and with a theoretical basis. 5. This study has shown the feasibility of constructing both the one-point statistics of autospectral models and the two-point statistics of coherence models from a database. These models could serve as surrogate analytical models in the experimental and CFD investigations; thus, this feasibility offers promise in providing an improved understanding of wake turbulence. 6. The two frameworks for the autospectrum and coherence increase the utility base of the database, involving enormous resources. Given the simple analytical structure of these models, they bring better understanding and transparency to a dataset.

Conclusions
This study has shown that an earlier-exercised mathematical framework lends itself well to extracting interpretive autospectral models of wind farm wake turbulence from a database. While these models are database specific, the framework can be applied to any database and the model construction is straightforward. As to the two-point statistics of wake-turbulence, this study first presents a unified account of cross-spectrum and coherence for HIT; this account is of considerable utility in that, in the literature, the expressions of cross-spectrum and coherence show differences from one study to the other. Given these expressions of coherences, this study then builds a framework for extracting wake turbulence coherence models from a database. The frameworks for autospectra and coherences do not follow classical perturbation theory approach of a solution for a linearized problem along with successively added corrections. Both frameworks represent a practical combination of a series expansion, exploitation of a database, and theoretical constraints in closed form.
This study also leads to following specific findings: 1. Generally, no more than a three-term series (second-order correction) is necessary to develop an autospectral model; in most cases, a two-term series (first-order correction) is found to be adequate for wind engineering applications.

2.
The addition of a third term to the series has significant power in reducing the "A error" between the model and the data. Recall that the "A error" refers to minimizing the sum of the errors in a least squares sense: the model deviation from the autospectral data points and from the measured high frequency autospectral decay level (when applicable).

3.
These developed models lend themselves well to design of filters driven by white noise; that is, the filter design is as routine as the currently used procedure for the von Kármán models.

4.
While this framework to constructing the coherence models from a database has not been tested against a database, it has been formulated from first principles and with a theoretical basis.

5.
This study has shown the feasibility of constructing both the one-point statistics of autospectral models and the two-point statistics of coherence models from a database. These models could serve as surrogate analytical models in the experimental and CFD investigations; thus, this feasibility offers promise in providing an improved understanding of wake turbulence. 6.
The two frameworks for the autospectrum and coherence increase the utility base of the database, involving enormous resources. Given the simple analytical structure of these models, they bring better understanding and transparency to a dataset.