*Article* **A Spectrum Correction Method Based on Optimizing Turbulence Intensity**

**Wenwu Yi 1,2, Ziqi Lu 1,2, Junbo Hao 1,2, Xinge Zhang 1,2, Yan Chen 1,2,\* and Zhihong Huang 1,2**


**\*** Correspondence: ychen@stu.edu.cn; Tel.: +86-138-2968-2901

**Abstract:** Based on the classical spectral representation method of simulating turbulent wind speed fluctuation, a harmonic superposition algorithm was introduced in detail to calculate the homogeneous turbulence wind field simulation in space. From the view of the validity of the numerical simulation results in MATLAB and the simulation efficiency, this paper discussed the reason for the bias existing between three types of turbulence intensity involved in the whole simulation process: simulated turbulence intensity, setting reference turbulence intensity, and theoretical turbulence intensity. Therefore, a novel spectral correction method of a standard deviation compensation coefficient was proposed. The simulation verification of the correction method was carried out based on the Kaimal spectrum recommended by IEC61400-1 by simulating the uniform turbulent wind field in one-dimensional space at the height of the hub of a 15 MW wind turbine and in two-dimensional space in the rotor swept area. The results showed that the spectral correction method proposed in this paper can effectively optimize the turbulence intensity of the simulated wind field, generate more effective simulation points, and significantly improve the simulation efficiency.

**Keywords:** uniform wind field simulation; turbulence intensity; deviation of standard deviation; spectral representation

### **1. Introduction**

The random turbulent wind is one of the critical factors causing the fatigue load of wind turbine blades. Aerodynamic analysis of wind turbines under the influence of turbulence is a vital basis for designing the ultimate load and fatigue load [1]. Meanwhile, with the trend of large-scale and flexible blades, the coupling of nonlinear aerodynamic and structural problems with the environment becomes more complex. Therefore, it is of great significance to establish a pulsating wind field model with high adaptability to multiple design conditions and complicated incoming flow environments in line with engineering applications.

In engineering practice, considering the speed of the solution and data processing is important. While using the computational fluid dynamics (CFD) method to solve Navier-Stokes equations results in a high accuracy, this method is complex and heavily dependent on computer performance. Accordingly, the current research on turbulent wind field simulation tends to be based on classical stochastic process theory. The PSD (power spectral density) function is employed to simulate the time history of pulsating wind speed. Regularly, Harris spectrum, Von Karman spectrum, Simiu spectrum, and Kaimal spectrum are applied to power spectral density models [2–6]. Shinozuka proposed the harmonic synthesis method to settle the matter of the stationary Gaussian random process and non-stationary random process of wind speed time history simulation, introducing the double index frequency combined with FFT technology to achieve the ergodic properties of each state of the simulation curve [7,8]. Lagrange interpolation, Hermite interpolation,

**Citation:** Yi, W.; Lu, Z.; Hao, J.; Zhang, X.; Chen, Y.; Huang, Z. A Spectrum Correction Method Based on Optimizing Turbulence Intensity. *Appl. Sci.* **2022**, *12*, 66. https:// doi.org/10.3390/app12010066

Academic Editors: Luis Hernández-Callejo, Sara Gallardo Saavedra and Sergio Nesmachnow

Received: 20 October 2021 Accepted: 6 December 2021 Published: 22 December 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

PoD-based interpolation, and other interpolation methods [9–11] have significantly reduced the number of analog signal decompositions and have further optimized the operation efficiency. In essence, the power spectrum simulation method is a Monte Carlo numerical statistical method. To date, most studies have mainly concentrated on how to optimize the algorithm to improve the computational efficiency in order to solve more dense simulation points, or to verify the fitting degree between the power spectrum simulating fluctuating wind speed and the original spectrum [12–14]. However, few scholars have paid attention to the turbulence intensity of the simulated fluctuating wind field and the validity of the time history of the simulating point fluctuating wind speed.

Based on applying the spectral representative method to simulate wind speed, the significant standard deviation's bias, generated due to valuing a truncated simulated frequency interval and the simulation processing itself, respectively, is proposed to shorten the time required to complete the method by creatively introducing a compensation coefficient. The corrective method modifies the original power spectral density function to produce more effective simulated points in the uniform simulated wind field that can generate the simulated turbulence intensity, meeting the set requirements. Moreover, this paper also discussed the correction method's further development by expanding its usage to the two-dimensional level. The modified method would be verified by simulating the uniform fluctuating wind field of a 15 MW referenced wind turbine, which may lay the foundation for simulating a turbulent field that is applied to a larger scale wind turbine.

### **2. Materials and Methods**

### *2.1. Turbulent Wind Field Model*

In a turbulent wind field, the wind speed can be decomposed into average wind speed and fluctuating wind speed. Thus, the simulated three-dimensional wind speed at any simulation point is the linear superposition of the average wind speed and fluctuating wind speed in longitudinal, transverse, and vertical directions [2], which can be calculated by:

$$\mathcal{V}(z,t) = \Box(z) + \upsilon(z,t) \tag{1}$$

where, *U*(*z*) denotes the average wind speed, *v*(*z*, *t*) denotes the fluctuating wind speed of a simulated point, z denotes simulated height, and t denotes time.

When setting the average wind speed, the influence of wind shear must be considered for large wind turbines whose installation height is usually higher than 100 m. Commonly used wind shear models include exponential model and logarithmic model, and the stable modified logarithmic model was used in this paper [8], which can be calculated by:

$$u(z) = \frac{u^\*}{\kappa} \left( \ln \frac{z}{z\_0} - \psi \right) \tag{2}$$

where, *u*<sup>∗</sup> denotes friction velocity, it can be represented as *u*<sup>∗</sup> = 0.045*Vref* − 0.012, *Vref* , which means the reference mean wind speed; *κ* denotes Von Karman constant, under neutral atmospheric conditions, *κ* = 0.4; *z*<sup>0</sup> denotes terrain roughness parameter, which can be calculated as *z*<sup>0</sup> = *Ac* (*u*∗) 2 *<sup>g</sup>* , where *Ac* = 0.034 when the terrain type is an offshore area, g denotes the acceleration of gravity; and *ψ* denotes the stability function when in the neutral condition, its value is 0.

Pulsating wind speed is regarded as a stationary Gaussian random process. According to the stochastic process theory, the power spectral density function combines with the coherence function, which is set to describe the wind speed correlation between two different points that may generate different wind speeds during a period in a stochastic wind field while a smaller separation distance of any two points will bring a greater correlation, is used to simulate the pulsating wind speed time history [15].

### *2.2. Basic Algorithm for Simulating Fluctuating Wind Speed Time History*

Considering the simulation process at a certain point *vj*(*t*),(*j* = 1, 2, 3, ... , *m*) in the turbulent wind field in space, the basic algorithm for simulating and solving the one-dimensional *M*-variable turbulent wind field ' *vj*(*t*) ( is described, and its algorithm flow is shown as follows:

	- Sampling frequency *fs*, in order to ensure that the analog signal can be reconstructed accurately without aliasing, *fs* is required to be greater than two times the value of *Fmax*, *Fmax* indicates the upper limit of the analog frequency range, which denotes also the cut-off frequency in this paper;
	- The initial frequency *Fmin* and cutoff frequency *Fmax*, which are selected by considering the dimensionless power spectral density function image and the influence of truncation error in the simulation frequency range on the simulation variance [17], the frequency step size is denoted by *d f* = *Fmax*−*Fmin <sup>N</sup>* , where *N* represents the frequency sampling number.
	- Sampling interval *dt*, *dt* = 1/ *fs*;
	- Analog time points *M*, generally *M* = 2*N*.

$$dr = \left(\frac{z}{m-1}\right) \cdot \sqrt{\left(x\_k - x\_j\right)^2 + \left(y\_k - y\_j\right)^2} \tag{3}$$

where, *xk*, *xj*, *yk*, *yj* represent the abscissa and ordinate values of simulation point *k* and simulation point *j*, respectively.

(4) Calculate the cross-spectral density matrix of each frequency sampling point {*S*(*fn*)}, (*n* = 1, 2, 3, ··· , *N*):

$$S(f\_n) = \begin{bmatrix} S\_{11} & S\_{12} & \cdots & S\_{1m} \\ S\_{21} & S\_{22} & \cdots & S\_{2m} \\ \vdots & \vdots & \ddots & \vdots \\ S\_{m1} & S\_{m2} & \cdots & S\_{mm} \end{bmatrix} \tag{4}$$

where, *S*(*f*) denotes the power spectral density function, *Sjj* denotes the autocorrelation spectrum of corresponding points, and *Skj*(*k* = *j*) denotes the cross-correlation spectrum between two simulated points, (*k* = 1, 2, 3, ··· , *j*), which can be calculated as:

$$S\_{kj}(f\_n) = \sqrt{S\_k(f\_n) S\_j(f\_n)} \cdot \operatorname{col}(k, j; f) \tag{5}$$

Since the cross-correlation function of the fluctuating wind speed in one-dimensional space does not involve the change of height, the cross-correlation spectrum between the simulated point j and other simulated points is consistent, and the cross-correlation spectrum of the fluctuating wind speed in one-dimensional space can be obtained by the following formula:

$$S\_{kj}(f\_n) = S\_j(f\_n) \cdot \text{coh}(k, j; f) \tag{6}$$

(5) Combined with the double-indexed frequency method, the Cholesky decomposition method is performed on the cross-spectral density matrix of sampling points at each frequency, and the decomposed lower triangular matrix is obtained, that is:

$$S(f\_{kn}) = H(f\_{kn})H^T(f\_{kn})\tag{7}$$

$$f\_{kn} = (n-1)df + (k/m)df\tag{8}$$

where, *fkn* denotes the double-indexed frequency; *H*(*fkn*) denotes the decomposed lower triangular matrix.

(6) Random phase obedience *φkn* is introduced and subjected to independent random distribution among intervals (0, 2*π*), carrying out the Fast Fourier transformation, namely:

$$G\_{ku}(t) = FFT(B\_{ku})\tag{9}$$

$$B\_{kn}(ndf) = \begin{cases} \ \ \end{cases} \begin{aligned} H(ndf + (k/m)df)e^{i\phi\_{kn}}, (0 \le n < N) \\ 0, (N \le n < M) \end{aligned} \tag{10}$$

(7) Employing harmonic superposition method to generate the wind speed time history of the simulated point *j* samples, and *vj*(*t*) can be calculated by Equation (11):

$$w\_j(p \cdot dt) = \sqrt{2df} \operatorname{Re} \sum\_{k=1}^{j} G\_{jk} \exp(i \frac{k}{m} \cdot 2\pi df \cdot p \cdot dt) \tag{11}$$

where, *p* = 1, 2, 3, ··· ,(*M* · *m* − 1).

### *2.3. Power Spectrum Correction Method*

In the practical application of simulating a uniform pulsating wind field, since the dimensionless power density spectrum of one-dimensional space accords with Gaussian distribution (Figure 1 shows the longitudinal dimensionless wind speed power density spectrum curve of several classical spectra), the relationship between the simulating pulsating wind speed variances *σ*<sup>2</sup> and the power spectrum density function *S*(*f*) is defined as below [18,19]:

$$
\sigma^2 = \int\_0^\infty S(f) df \tag{12}
$$

Therefore, when the analog frequency range [*F*max, *Fmin*] is a segment interval and is cut out from the spectral density function originally defined interval (0, ∞), theoretically, the variance of the simulated fluctuating wind speed *σtheory*<sup>2</sup> is defined as:

$$\left| \sigma\_{\text{theory}} \right|^2 = \int\_{F\_{\text{min}}}^{F\_{\text{max}}} S(f) df \tag{13}$$

The defined standard deviation of the fluctuating wind speed *σ* can be determined by an empirical formula in the simulation algorithm (e.g., Equation (20) in Section 3.1), according to the calculation formula of turbulence intensity *IT* = *σ*/*v*(*t*), choosing the segment interval [*F*max, *Fmin*] will cause the simulated results to produce truncated standard deviation, which will inevitably affect the numerical error of the simulated turbulence intensity due to the existence of the deviation of *σ* and *σtheory*. The deviation between the two can be calculated as:

$$
\varepsilon = \left| \frac{\sigma\_{\text{theory}} - \sigma}{\sigma} \right| \tag{14}
$$

The standard deviation's deviation *ε* affects the simulated turbulence intensity to a certain extent, and the smaller the deviation is, the smaller the influence is.

**Figure 1.** The longitudinal dimensionless velocity of the wind power density spectrum curve.

It can be seen from Equations (12)–(14) that when increasing the proportion of the interval length in the simulated frequency interval, the theoretically calculated standard deviation's deviation will decrease. However, it can be seen from the simulation algorithm introduced in Section 2.2 that the frequency step has a significant impact on the simulation results. On the other hand, according to the formula *d f* = *Fmax*−*Fmin <sup>N</sup>* , after determining the basic range of the simulation frequency interval, a slight change caused by the minor increase in cutoff frequency relative to the magnitude of *N* is almost negligible. When the ratio of the turbulence integral length to the average wind speed *Xk* = *Lk*/*U* is constant, the upper limit of the interval expands to the right indefinitely, and the deviation will not significantly decrease [10].

This is also consistent with the simulation conclusion in this paper: under the same simulation conditions, when the initial frequency changes by one step unit, the impact on the standard deviation's deviation *ε* is significantly greater than that when the cutoff frequency changes by one step unit (see Section 3.2 for details). In addition, the simulation results also show that the theoretical truncation bias *ε* caused by the truncation of the simulation frequency range will affect the numerical value of the simulated object in the simulation process, which leads to the simulation bias *ε*ˆ. Simulation deviation is not only related to truncation deviation, but also to the setting of other simulation parameters and the application of interpolation methods. The simulated deviation is a manifestation of the error generated by simulating the turbulence intensity, which can be expressed by the following calculation formula:

$$
\mathfrak{E} = \left| \frac{\partial \mathfrak{T} - \sigma}{\sigma} \right| \tag{15}
$$

where, *σ*ˆ denotes the standard deviation of actual simulated wind speed.

Considering df = *Fmax*−*Fmin <sup>N</sup>* , if the magnitude of the frequency sample *N* is increased by an order (usually *N* is an exponential form with base 2), the simulation bias can indeed be reduced, but at the same time the computational memory will be doubled and the simulation speed will be seriously slowed down. Take the analog frequency range [*F*max, *Fmin*]=[10<sup>−</sup>5, 12], and the other simulation conditions are consistent with the case of one-dimensional simulation in Section 3. The simulation experiment results are shown in Table 1.

Table 1 and Figure 2 show that increasing the scale *N* causes the whole standard deviation bias of the turbulent wind field simulation to almost linearly reduce, or even when increased to *N* = 214, the overall standard deviation's bias of the turbulent wind field has reached the truncation error, but with a significant disadvantage of longer running time, and a great computing cost if more simulated points are set in the turbulent wind field space. On the other hand, this result is not very ideal. When increasing the magnitude of *N*, the simulated turbulence intensity of more simulation points falls outside the range of reference turbulence intensity and the defined turbulence intensity and needs to be removed, resulting in an unnecessary waste of computing resources.

(**a**) (**b**)

**Figure 2.** (**a**) The deviation of standard deviation at a different frequency sampling magnitude; (**b**) Simulated turbulence intensity at a different frequency sampling magnitude.



To sum up, the correction method should realize the minimum truncation deviation *εmin* in theory and the minimum simulation deviation *ε*ˆ*min* generated in the actual simulation process, without adding more frequency sampling points to reduce the simulation efficiency. In addition, the method should generate as many effective simulation points as possible (effective simulation points are defined as the simulation points whose simulated turbulence intensity ˆ*IT* falls within the range of reference turbulence intensity *Iref* and defined turbulence intensity *IT*).

Consider the error values of truncation bias *ε* and simulation bias *ε*ˆ:

$$
\Delta \varepsilon = \pounds - \varepsilon \tag{16}
$$

Introduce compensation coefficient *β*:

$$
\beta = \left(1 + \Delta \varepsilon\right)^2 \tag{17}
$$

If the variance *σ*<sup>2</sup> in the PSD function is substituted into the compensation coefficient *β*, then the original factor will be corrected *<sup>σ</sup>*<sup>2</sup> *β* .

### **3. Results and Discussion**

Given the IEA 15 MW (NREL) wind turbine [20], its basic parameters are shown in Table 2, and frequency domain simulation parameter settings as shown in Table 3.


**Table 2.** 15 WM wind turbine parameters.

**Table 3.** Frequency domain simulation parameters.


### *3.1. Power Density Spectrum and Coherence Function Selection*

Kaimal spectrum is widely used to describe the spectral density during wind speed [21], and its simulation results are more conservative when applied to the fatigue load design of wind turbines. The function expression of Kaimal spectrum in IEC61400-1 4th Edition [18] is as below, *Xk* = *Lk*/*U* represents the ratio of the turbulence integral length to average wind speed (*k* = 1, 2, 3 represents the longitudinal, transverse and vertical component spectra, respectively):

$$S(f) = \sigma\_k^2 \frac{4X\_k}{\left(1 + 6fX\_k\right)^{5/3}}\tag{18}$$

IEC61400-1 uses an exponential form of coherence function, which is expressed as follows:

$$\text{Coh}(r,f) = \exp[-12((\frac{fr}{u})^2 + (0.12\frac{r}{L\_k})^2)^{0.5}]\tag{19}$$

In Equation (18), according to IEC61400-1, the parameters of the longitudinal component are calculated as follows:

(1) Wind speed definition standard deviation *σ*1:

$$
\sigma\_1 = I\_{ref} (0.75 \Pi + b); \ b = 5.6 m/s \tag{20}
$$

where, *Iref* = 0.14, determined by 15 MW wind turbine design grade Class IB; *U* = 10.59 *m*/*s*.

(2) Integral scale parameter *L*1:

$$L\_1 = 8.1\Lambda\_1\tag{21}$$

The longitudinal integral scale parameters at the height of the hub *Z* are:

$$\Lambda\_1 = \begin{cases} \ 0.7Z\_\prime (Z \le 60 \,\, m) \\ \ 42 \,\, m\_\prime (Z \ge 60 \,\, m) \end{cases} \tag{22}$$

The horizontal and vertical components of turbulence spectrum parameters simulating a three-dimensional fluctuating wind speed are shown in Table 4.

**Table 4.** Turbulence spectrum parameters of horizontal and vertical components.


### *3.2. Determine the Analog Frequency Range* [*Fmax*,*Fmin*]

By setting the initial frequency *Fmin* <sup>∈</sup> <sup>1</sup> 10−5, 10−<sup>1</sup> 2 , setting the cut-off frequency *Fmax* <sup>∈</sup> [1, 10] (step by 101 units and 1 unit respectively to form an 5 <sup>×</sup> 10 analog frequency interval matrix), and performing integral operations one by one according to Equation (13), the ratio of the standard deviation of the truncated interval integral to the standard deviation of the fully defined interval integral is *η* = *σtheory*/*σ*, the result of which is shown in Figure 3a. By combining this with the 3*σ* principle, it can be found that when the initial frequency *Fmin* < 10−<sup>4</sup> and cut-off frequency is given any value, the standard deviation ratio *η* will always fall within a confidence interval of +2*σ*, as the cutoff frequency increases, *η* infinity goes to +3*σ*, therefore, it achieves a value of *Fmin* = 10−4.

Under the same simulation conditions, the simulation frequency interval element from the matrix described above was substituted, one by one, into the uniform turbulent wind field simulation algorithm, for simulating fluctuating wind speed of time history in one-dimensional space (at hub height), accompanied with the standard deviation's deviation and truncation deviation. The comparison between the two is shown in Figure 4b. The figure shows that when *Fmin* = 10−4, the simulation bias *ε*ˆ does not decrease significantly with the increase of the simulation frequency range, the truncated bias step rises when *Fmax* > 6, instead. To weaken the distortion of the analog signal, the cutoff frequency *Fmax* = 5 was selected to simulate the uniform turbulent wind field in one-dimensional space (at the height of the hub) and two-dimensional space (at the swept surface of the wind wheel), respectively, (the distribution of simulated points of these two types of spatial wind field simulation are shown in Figure 4a,b, as above) considering the truncation error and simulation error comprehensively, and the compensation coefficient was determined.

**Figure 3.** (**a**) The ratio of cut-off standard deviation and defined standard deviation along with the cut-off frequency change curve; (**b**) Standard deviation scatter distribution under different cut-off frequency intervals.

**Figure 4.** (**a**) One-dimensional spatial simulation point distribution diagram; (**b**) Two-dimensional spatial simulation point distribution diagram.

n

### *3.3. Compensation Coefficient β Correction Method*

According to Section 3.2, the time history simulation of the fluctuating wind speed in a turbulent wind field is performed under the condition that the simulated frequency range is determined as [*F*max, *Fmin*]=[10−4, 5]. The main parameters are shown in Table 5, and the distribution of the standard deviation error and compensation coefficient at each simulated height in space is obtained, as shown in the Figure 5.

The longitudinal Kaimal spectrum is modified by considering the compensation coefficient *β*:

$$S^0(f) = \frac{\sigma\_1^2}{\beta} \frac{4X\_1}{(1 + 6fX\_1)^{5/3}}\tag{23}$$

where, *X*<sup>1</sup> = *L*1/*U*.

**Table 5.** The compensation coefficient of relevant parameters.


**Figure 5.** The distribution of correction factor of simulated points for two-dimensional space.

After compensation and correction, the main parameters of one-dimensional space are shown in Table 6.



The modified simulated point fluctuating wind speed diagram is shown as Figure 6.

**Figure 6.** 10 min fluctuating wind speed time history diagram of points: (**a**) One-dimensional space for points 20, 21, 40; (**b**) Two-dimensional space for points 1, 25, 49.

A comparison of the power spectral density at simulated points and power spectral density related to simulated wind speed before and after modification is shown in Figure 7. Among which, the line of the 'Corrected Spectrum Simulated' represents the power density spectrum of the fluctuating wind speeds that are generated by the corrected Kaimal spectrum; the line of the 'Corrected Spectrum' represents the corrected Kaimal power density spectrum itself; the line of the 'Original Simulated Spectrum' represents the power density spectrum of the fluctuating wind speed that are generated by the original Kaimal spectrum; the line of the 'Original Spectrum' represents the original Kaimal power density spectrum itself.

Comparison of the standard deviation, standard deviation's bias, and turbulence intensity of the simulated wind speed before and after modification is shown in Figure 8. Among which, in Figure 8a,c,e, the line of the 'Corrected Simulation' represents the standard deviation/standard deviation's deviation/turbulence intensity of the fluctuating wind speed that is simulated by corrected Kaimal spectrum and generated by 40 points in a one-dimensional simulated wind field; the line of the 'Original Simulation' represents the standard deviation/standard deviation's deviation/turbulence intensity of the fluctuating wind speed that is simulated by the Kaimal spectrum and generated by 40 points in a one-dimensional simulated wind field; the line of the 'Corrected mean value' represents the mean value of the simulation results from the corrected Kaimal spectrum among 40 points in a one-dimensional simulated wind field; the line of the 'Original mean value' represents the mean value of the simulation results from the Kaimal spectrum among 40 points in a one-dimensional simulated wind field; in Figure 8e, the line named 'Defined' represents the defined turbulence intensity *IT*; the line named 'Referenced' represents the reference turbulence intensity *Iref* ; in Figure 8b,d,f, the surface named 'Corrected'

represents the standard deviation/standard deviation's deviation/turbulence intensity of the fluctuating wind speed of the points varying along with the simulated heights generated by the corrected Kaimal spectrum in a two-dimensional wind field; the surface named 'Uncorrected' represents the standard deviation/standard deviation's deviation/turbulence intensity of the fluctuating wind speed of the points varying along with the simulated heights generated by the Kaimal spectrum in a two-dimensional wind field; and the surface named 'Defined' represents the defined standard deviation/turbulence intensity of the fluctuating wind speed of the points varying along the simulated heights.

**Figure 7.** The power density spectrum fitting: (**a**) One-dimensional simulated point 20 being corrected and uncorrected; (**b**) Two-dimensional simulated point 1 being corrected.; (**c**) One-dimensional simulated point 21 being corrected and uncorrected; (**d**) Two-dimensional simulated point 25 being corrected; (**e**) One-dimensional simulated point 40 being corrected and uncorrected; (**f**) Two-dimensional simulated point 49 being corrected.

Concerning the longitudinal one-dimensional space spectrum correction method and the lateral spectrum and the vertical spectrum parameters in Table 4, 20 evenly distributed simulation points are set to verify whether the compensate correction method suggested is applicable for the lateral spectrum and the vertical spectrum. A one-dimensional simulation of three-dimensional wind speed is acquired, and the three-dimensional fluctuating wind velocity contour map and 3 d surface figure of simulation point 10 are shown in Figure 9.

**Figure 8.** The contrast of standard deviation of simulated points between being corrected and uncorrected: (**a**) One-dimensional space; (**b**) Two-dimensional space; The contrast of standard deviation deviation of simulated points between being corrected and uncorrected: (**c**) One-dimensional space; (**d**) Two-dimensional space; The contrast of turbulence intensity of simulated points between being corrected and uncorrected: (**e**) One-dimensional space; (**f**) Two-dimensional space.

**Figure 9.** (**a**) The three-dimensional fluctuating wind velocity contour map of point 10; (**b**) Standard The three-dimensional fluctuating wind velocity 3D surface of point 10.

The above results show that:


### **4. Conclusions**


**Author Contributions:** Conceptualization, W.Y., Z.L.; Methodology, W.Y.; Data curation, W.Y., J.H.; Formal analysis, W.Y.; Investigation, W.Y., X.Z.; Project administration, Y.C.; Software, W.Y.; Supervision, Z.H.; Validation, W.Y.; Visualization, W.Y., J.H.; Writing—original draft, W.Y., Z.L.; Writing—Review and editing, W.Y., X.Z., Z.L., J.H.; All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the National Natural Science Foundation of China, grant number 51976113.

**Institutional Review Board Statement:** Not Applicable.

**Informed Consent Statement:** Not Applicable.

**Data Availability Statement:** Not Applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.

### **References**

