Doppler Lidar Investigations of Wind Turbine Near-Wakes and LES Modeling with New Porous Disc Approach

: Many bottom-mounted offshore wind farms are currently planned for the coastal areas of Japan, in which wind speeds of 6.0–10.0 m/s are extremely common. The impact of such wind speeds is very relevant for the realization of bottom-mounted offshore wind farms. In evaluating the feasibility of these wind farms, therefore, strict evaluation at wind speeds of 6.0–10.0 m/s is important. In the present study, the airﬂow characteristics of 2 MW-class downwind wind turbine wake ﬂows were ﬁrst investigated using a vertically proﬁling remote sensing wind measurement device (lidar). The wind turbines used in this study are installed at the point where the sea is just in front of the wind turbines. A ground-based continuous-wave (CW) conically scanning wind lidar system (“ZephIR ZX300”) was used. Focusing on the wind turbine near-wakes, the detailed behaviors were considered. We found that the inﬂuence of the wind turbine wake, that is, the wake loss (wind velocity deﬁcit), is extremely large in the wind speed range of 6.0–10.0 m/s, and that the wake loss was almost constant at such wind speeds (6.0–10.0 m/s). It was additionally shown that these results correspond to the distribution of the thrust coefﬁcient of the wind turbine. We proposed a computational ﬂuid dynamics (CFD) porous disk (PD) wake model as an intermediate method between engineering wake models and CFD wake models. Based on the above observations, the wind speed range for reproducing the behavior of the wind turbine wakes with the CFD PD wake model we developed was set to 6.0–10.0 m/s. Targeting the vertical wind speed distribution in the near-wake region acquired in the “ZephIR ZX300”, we tuned the parameters of the CFD PD wake model (C RC = 2.5). We found that in practice, when evaluating the mean wind velocity deﬁcit due to wind turbine wakes, applying the CFD PD wake model in the wind turbine swept area was very effective. That is, the CFD PD wake model can reproduce the mean average wind speed distribution in the wind turbine swept area.


Introduction
Generally, as the blade of a wind turbine rotates, a velocity loss region is formed downstream of the turbine (see Figure 1). In a large-scale offshore wind farm, two major problems have been highlighted regarding the wind turbines located on the downstream side due to the wind turbine wakes located on the upstream side: (i) a decrease in the amount of power generated by the wind turbine located on the downstream side and (ii) the occurrence of failures inside the wind turbine nacelle and wind turbine blades.
A review of the state-of-the-art in wind turbine wake modeling is organized in detail in the literature [1][2][3]. The literature also describes in detail various approaches to the behavior of wind turbine wakes [4][5][6][7][8][9]. The most common engineering wake model/analytical wake model for wind farms is the so-called Park wake model by Jensen [8] and Katic [9]. The Park wake model implements a simple formula for the size of the wake deficit and its expansion downstream with a single adjustable parameter (wake decay constant). Therefore, the Park wake model is still widely used for evaluating business feasibility in the wind industry in Japan and overseas. Engineering wake models are extremely difficult to accurately predict the mutual interference of wind turbine wakes.
We have already proposed the CFD PD wake mode [10]. In the present study, the validity of the CFD PD wake model was verified by the wake of a 2 MW-class downwind turbine. The airflow characteristics of 2 MW-class downwind wind turbine wake flows were first investigated using a vertically profiling remote sensing wind measurement device (lidar). The wind turbines used in this study are installed at the point where the sea is just in front of the wind turbines. A ground-based CW conically scanning wind lidar system ("ZephIR ZX300") [11,12] was used for lidar. Particular attention was concentrated on the behavior of the wind turbine near-wakes. Next, targeting the near-wake flow characteristics acquired in the "ZephIR ZX300", we tuned the parameters of the wake model we proposed and considered the numerical results, including those of the conventional engineering wake models.

Overview of Omonogawa Wind Power Station in Akita Prefecture
The Omonogawa Wind Power Station, owned by Hitachi Zosen Corporation, is located on the north side of the mouth of the Omonogawa River in Akita City, Akita Prefecture, Japan. As shown in Figure 2, two downwind-type wind turbines with a rated output of 2-MW, manufactured by Hitachi, Ltd., are in operation: the first, since 1 March 2015 and the second, since 1 February 2016. As shown in Figure 2a The predominant wind direction in this area is west-northwest, and it is a major feature that the two wind turbines are arranged almost in a straight line. Furthermore, wind turbine No.2 is located about 80 m from the coastline. When a winter monsoon (main wind direction: west-northwest to northwest) occurs, the wind environment in this area is similar

Overview of Omonogawa Wind Power Station in Akita Prefecture
The Omonogawa Wind Power Station, owned by Hitachi Zosen Corporation, is located on the north side of the mouth of the Omonogawa River in Akita City, Akita Prefecture, Japan. As shown in Figure 2, two downwind-type wind turbines with a rated output of 2-MW, manufactured by Hitachi, Ltd., are in operation: the first, since 1 March 2015 and the second, since 1 February 2016. As shown in Figure 2a The predominant wind direction in this area is west-northwest, and it is a major feature that the two wind turbines are arranged almost in a straight line. Furthermore, wind turbine No.2 is located about 80 m from the coastline. When a winter monsoon (main wind direction: west-northwest to northwest) occurs, the wind Energies 2021, 14, 2101 3 of 33 environment in this area is similar to the wind environment that would be found above the ocean. In this study, the wind field from the sea (sea breeze) was the subject of research. Figure 3 shows the various aerodynamic characteristics of downwind wind turbines with a rated output of 2-MW targeted in this study. The wind speed at which the power coefficient shows the maximum value (about 0.48) is around 8.0 m/s, and the thrust coefficient corresponding to this is about 0.8. In this study, we focused on the airflow characteristics of the wind turbine wake at 6.0-10.0 m/s. Energies 2021, 14, x FOR PEER REVIEW 3 of 33 to the wind environment that would be found above the ocean. In this study, the wind field from the sea (sea breeze) was the subject of research. Figure 3 shows the various aerodynamic characteristics of downwind wind turbines with a rated output of 2-MW targeted in this study. The wind speed at which the power coefficient shows the maximum value (about 0.48) is around 8.0 m/s, and the thrust coefficient corresponding to this is about 0.8. In this study, we focused on the airflow characteristics of the wind turbine wake at 6.0-10.0 m/s.     Figure 4 shows the positional relationship between the two wind turbines, the Doppler lidar, and the wind measurement by the 50 m observation mast. In this study, a British-made Doppler lidar "ZephIR ZX300" was used. The "ZephIR ZX300" can measure a total of 10 points of wind direction, wind speed and turbulence intensity at any altitude from about 40 m to about 200 m.  Figure 4 shows the positional relationship between the two wind turbines, the Doppler lidar, and the wind measurement by the 50 m observation mast. In this study, a Britishmade Doppler lidar "ZephIR ZX300" was used. The "ZephIR ZX300" can measure a total of 10 points of wind direction, wind speed and turbulence intensity at any altitude from about 40 m to about 200 m.

Verification of Accuracy of Doppler Lidar Measurement
First, we verified the accuracy of the Doppler lidar measurement. As shown in Figure  4, the Doppler lidar system was installed downstream of wind turbine No.2 at 3D and 4D, where D (=86 m) means the rotor diameter of the wind turbine No.2. All the measurement data shown below were 10-minutes mean statistics. Figure 5 shows the situation when the Doppler lidar is installed at 3D and 4D distance from wind turbine No.2. Only plants are present between the wind turbine and the Doppler lidar, as shown in Figure 5a.
In the case where the Doppler lidar system was installed downstream of wind turbine No.2 at 3D, we compared the nacelle-measured wind speed and turbulence intensity of turbine No.2 with the wind speed and turbulence intensity by the Doppler lidar for the westerly wind (270.0 ± 45.0 deg), which is not affected by the wind turbine wakes (see Figure 4). The hub height (78 m) of wind turbine No.2 was considered as the reference altitude. Considering the balance of measurement points, 46, 57, 78, 88, 98, 108, 118, 120, 128, 138 and 207 m were set as the measurement points.
On the other hand, in the case where the Doppler lidar system was installed downstream of wind turbine No.2 at 4D, we also compared the 50 m observation mast data and the Doppler lidar data with respect to the above wind direction (270.0 ± 45.0 deg). The 50 m observation mast was considered as the reference altitude. Considering the balance of

Verification of Accuracy of Doppler Lidar Measurement
First, we verified the accuracy of the Doppler lidar measurement. As shown in Figure 4, the Doppler lidar system was installed downstream of wind turbine No.2 at 3D and 4D, where D (=86 m) means the rotor diameter of the wind turbine No.2. All the measurement data shown below were 10-minutes mean statistics. Figure 5 shows the situation when the Doppler lidar is installed at 3D and 4D distance from wind turbine No.2. Only plants are present between the wind turbine and the Doppler lidar, as shown in Figure 5a.
In the case where the Doppler lidar system was installed downstream of wind turbine No.2 at 3D, we compared the nacelle-measured wind speed and turbulence intensity of turbine No.2 with the wind speed and turbulence intensity by the Doppler lidar for the westerly wind (270.0 ± 45.0 deg), which is not affected by the wind turbine wakes (see Figure 4). The hub height (78 m) of wind turbine No.2 was considered as the reference altitude. Considering the balance of measurement points, 46, 57, 78, 88, 98, 108, 118, 120, 128, 138 and 207 m were set as the measurement points.
On the other hand, in the case where the Doppler lidar system was installed downstream of wind turbine No.2 at 4D, we also compared the 50 m observation mast data and the Doppler lidar data with respect to the above wind direction (270.0 ± 45.0 deg). The 50 m observation mast was considered as the reference altitude. Considering the balance  Figure 6 shows the verification results of the accuracy of the Doppler lidar measurement. In this study, the comparison of wind speed and turbulence intensity were performed. As shown in Figure 6a (that is, correlation between the measurements made by the nacelle-mounted wind vane and anemometer and by the Doppler lidar for the wind speed and turbulence intensity), no clear correlation was found for the turbulence intensity due to the influence of the nacelle (correlation coefficient R = 0.59, coefficient of determination R 2 = 0.35). On the other hand, a good correlation was confirmed for the wind speed (R = 0.98, R 2 = 0.96). On the other hand, in comparison with the 50 m observation mast shown in Figure 6b (that is, correlation between the measurements made by the 50 m observation mast and by the Doppler lidar for the wind speed and turbulence intensity), it was confirmed that there is a very high correlation between the wind speed (R = 0.99, R 2 = 0.98) and the turbulence intensity (R = 0.94, R 2 = 0.88). Therefore, the Doppler lidar "ZephIR ZX300" was found suitable to measure the wind speed and the turbulence intensity.  Figure 6 shows the verification results of the accuracy of the Doppler lidar measurement. In this study, the comparison of wind speed and turbulence intensity were performed. As shown in Figure 6a (that is, correlation between the measurements made by the nacellemounted wind vane and anemometer and by the Doppler lidar for the wind speed and turbulence intensity), no clear correlation was found for the turbulence intensity due to the influence of the nacelle (correlation coefficient R = 0.59, coefficient of determination R 2 = 0.35). On the other hand, a good correlation was confirmed for the wind speed (R = 0.98, R 2 = 0.96). On the other hand, in comparison with the 50 m observation mast shown in Figure  6b (that is, correlation between the measurements made by the 50 m observation mast and by the Doppler lidar for the wind speed and turbulence intensity), it was confirmed that there is a very high correlation between the wind speed (R = 0.99, R 2 = 0.98) and the turbulence intensity (R = 0.94, R 2 = 0.88). Therefore, the Doppler lidar "ZephIR ZX300" was found suitable to measure the wind speed and the turbulence intensity.

Doppler Lidar Investigation of Wind Turbine Near-Wakes and Consideration
The main purpose of this study was to understand the airflow characteristics of the wind turbine near-wake in a pseudo-offshore environment. We focused on the period in which the monsoon is generated from the sea in winter (28 November 2016-2 December 2016 and 17 November 2017-21 December 2017), the Doppler lidar "ZephIR ZX300" was installed downstream of wind turbine No.2 at 3D (D = 86 m). As shown in Figure 4, we analyzed the measured data with a wind direction of 332.0 ± 10.0 deg (prevailing wind direction at this time), and all outliers and missing values were removed. Figure 7 shows the vertical distribution of the mean wind speed downstream of wind turbine No.2 at 3D (D = 86 m); in this figure, the results are organized in increments of 1 m/s (see Table 1) and the data are averaged at each altitude. For the wind speed shown in Table  1, the nacelle-measured wind speed was corrected with the Doppler lidar output at the hub height (78 m) of wind turbine No.2 as the true value. Furthermore, for example, for a wind speed of 5 m/s, alle wind speeds in the range of 5.00 m/s to 5.99 m/s were counted. Figure 7

Doppler Lidar Investigation of Wind Turbine Near-Wakes and Consideration
The main purpose of this study was to understand the airflow characteristics of the wind turbine near-wake in a pseudo-offshore environment. We focused on the period in which the monsoon is generated from the sea in winter (28 November 2016-2 December 2016 and 17 November 2017-21 December 2017), the Doppler lidar "ZephIR ZX300" was installed downstream of wind turbine No.2 at 3D (D = 86 m). As shown in Figure 4, we analyzed the measured data with a wind direction of 332.0 ± 10.0 deg (prevailing wind direction at this time), and all outliers and missing values were removed. Figure 7 shows the vertical distribution of the mean wind speed downstream of wind turbine No.2 at 3D (D = 86 m); in this figure, the results are organized in increments of 1 m/s (see Table 1) and the data are averaged at each altitude. For the wind speed shown in Table 1, the nacelle-measured wind speed was corrected with the Doppler lidar output at the hub height (78 m) of wind turbine No.2 as the true value. Furthermore, for example, for a wind speed of 5 m/s, alle wind speeds in the range of 5.00 m/s to 5.99 m/s were counted. Figure 7 reveals that there is a clear difference in the formation of wind turbine wakes at each wind speed. That is, there is a significant deficit in each wind speed due to wind turbine wakes. Considering this together with Figure 3, it can be seen that at wind speeds of about 8.0-9.0 m/s, where the power coefficient is maximized, the amount of wake loss near the height of the wind turbine hub (78 m above the ground) is extremely large compared to the other wind speed classes. The main reasons for this are as follows: at wind speeds of about 8.0-9.0 m/s, the wind turbine operates most efficiently in obtaining wind energy, so the amount of wake loss (wind velocity deficit) is considered to be large. When the wind speed is about 10.0 m/s or more, the amount of wake loss near the height of the wind turbine hub (78 m above the ground) gradually decreases as the wind speed increases. We consider that this is because the pitch angle is controlled so that the wind turbine blade does not become a source of resistance to the wind (so-called pitch angle control).
Wake loss was clearly confirmed in the wind speed range of 6.0-10.0 m/s ( Figure 7); Figure 8 shows the results for this wind speed range normalized by the wind speed near the wind turbine hub height (ground clearance of 78 m). It should be noted that the measurement results for wind speeds of 6.0-10.0 m/s in Figure 8 have almost the same vertical distribution in the range of the wind turbine swept area (z = 35-121 m); in this range, there is no significant difference in the amount of wake loss regardless of the wind speed. This tendency corresponds to the distribution of the thrust coefficient of the wind turbine shown in Figure 3b. Based on these results, in this study, the wind speed range targeted by the CFD PD wake model was set to 6.0-10.0 m/s (see Figure 3b). Figure 9a shows the vertical distribution of the turbulence intensity at 3D (D = 86 m) distance from wind turbine No.2 corresponding to Figure 7. It can be seen that there is a clear difference in the vertical distribution of the turbulent intensity at each wind speed. When the wind speed was 11.0 m/s or higher, the turbulence intensity gradually decreased as the wind speed increased. Figure 9b shows the results for wind speeds of 6.0-10.0 m/s, which were of particular interest in this study. In this range, the turbulence intensity near the height of the wind turbine hub (78 m above the ground) is very large compared to other wind speeds. Furthermore, as shown in Figure 8, the vertical distribution is almost the same in the range of the wind turbine swept area (z = 35-121 m), regardless of the wind speed. In this study, it was shown that the amount of wake loss was very large in the wind speed range of 6.0-10.0 m/s, compared to other wind speeds, and that the turbulence intensity was also large in the range of the wind turbine swept area. Figure 10 shows a comparison of the measured turbulence intensity with the International Electrotechnical Commission (IEC) model. Measurement point is 3D (D = 86 m) distance from wind turbine No.2 and wind turbine hub height 78 m above the ground. Measurement period is 18 November 2017-20 December 2017. It can be seen that the turbulence intensity is high at wind speeds of 4.0-10.0 m/s and exceeds the most severe turbulence category (Category A). As shown in Figure 3a, the wind turbine starts at a wind speed of 4.0 m/s and operates at a rated speed at wind speeds beyond 10.0 m/s. The turbulence intensity decreases at wind speeds higher than 11.0 m/s; as mentioned earlier, it is presumed that this is because the pitch angle is controlled so that the wind turbine blade does not become a source of resistance.

Comparison of Wind Measurement Results for Wind Turbine Near-Wake According to Doppler Lidar and Existing Engineering Wake Model
Here, we describe a comparison between the wind measurement results obtained by Doppler lidar and existing engineering wake models such as the Park model [8,9] and Ishihara model [13] derived from the Park model. Details of the Park model [8,9] formulation can be found in Appendix A. Figure 11 shows the visualization results for the Park model (wind speed in the streamwise direction) at the wind turbine hub height (78 m above the ground) when the wind speed is 8.0 m/s. It is confirmed that the wake effect extends to the vicinity of the Doppler lidar installed 3D downstream of wind turbine No.2. Figure 12 shows the wind speed prediction from the Park model (wind speed in the streamwise direction) at the wind turbine hub height (78 m above the ground) when the wind speed is 5.0-15.0 m/s. Now, please refer to the results of wind speeds of 6.0-10.0 m/s in Figure 12, which are our particular focus in this study. The wind speed ratios for the speeds 6.0-10.0 m/s are very similar, being the thrust coefficient of the wind turbine (see Figure 3b) almost constant in this range.

Comparison of Wind Measurement Results for Wind Turbine Near-Wake According to Doppler Lidar and Existing Engineering Wake Model
Here, we describe a comparison between the wind measurement results obtained by Doppler lidar and existing engineering wake models such as the Park model [8,9] and Ishihara model [13] derived from the Park model. Details of the Park model [8,9] formulation can be found in Appendix A. Figure 11 shows the visualization results for the Park model (wind speed in the streamwise direction) at the wind turbine hub height (78 m above the ground) when the wind speed is 8.0 m/s. It is confirmed that the wake effect extends to the vicinity of the Doppler lidar installed 3D downstream of wind turbine No.2. Figure 12 shows the wind speed prediction from the Park model (wind speed in the streamwise direction) at the wind turbine hub height (78 m above the ground) when the wind speed is 5.0-15.0 m/s. Now, please refer to the results of wind speeds of 6.0-10.0 m/s in Figure 12, which are our particular focus in this study. The wind speed ratios for the speeds 6.0-10.0 m/s are very similar, being the thrust coefficient of the wind turbine (see Figure 3b) almost constant in this range.    . A detailed examination of the obtained results revealed that there was a significant difference between the prediction results from the Park and Ishihara models and the measurement results from the Doppler lidar. In other words, it was shown that the prediction results from both the Park and Ishihara models gave a larger amount of wake loss (wind velocity deficit) compared to those from the Doppler lidar.   . A detailed examination of the obtained results revealed that there was a significant difference between the prediction results from the Park and Ishihara models and the measurement results from the Doppler lidar. In other words, it was shown that the prediction results from both the Park and Ishihara models gave a larger amount of wake loss (wind velocity deficit) compared to those from the Doppler lidar.  . A detailed examination of the obtained results revealed that there was a significant difference between the prediction results from the Park and Ishihara models and the measurement results from the Doppler lidar. In other words, it was shown that the prediction results from both the Park and Ishihara models gave a larger amount of wake loss (wind velocity deficit) compared to those from the Doppler lidar.   The wind speed range is 5.0-15.0 m/s. Note the measurement results obtained by Doppler lidar shown by the solid red line. As already mentioned, the wind velocity deficit is almost the same at wind speeds of 6.0-10.0 m/s, to which we paid particular attention in this study. Next, we focus on the prediction results from the Park and Ishihara models. In both cases, it was clear that there was a larger amount of wake loss (wind velocity deficit) near the wind turbine hub height (78 m above the ground), compared with the wind measurement results obtained by the Doppler lidar. Comparing the predicted value from the Park model with that from the Ishihara model, it is apparent that the wind velocity deficit from the Ishihara model becomes large in the wind speed range of 6.0 m/s or more.
Energies 2021, 14, x FOR PEER REVIEW 13 of 33 Figure 14 shows a comparison of the wind velocity deficit. Here, measurement point is 3D distance from wind turbine No.2 (D = 86 m) and hub height 78 m above the ground. The wind speed range is 5.0-15.0 m/s. Note the measurement results obtained by Doppler lidar shown by the solid red line. As already mentioned, the wind velocity deficit is almost the same at wind speeds of 6.0-10.0 m/s, to which we paid particular attention in this study. Next, we focus on the prediction results from the Park and Ishihara models. In both cases, it was clear that there was a larger amount of wake loss (wind velocity deficit) near the wind turbine hub height (78 m above the ground), compared with the wind measurement results obtained by the Doppler lidar. Comparing the predicted value from the Park model with that from the Ishihara model, it is apparent that the wind velocity deficit from the Ishihara model becomes large in the wind speed range of 6.0 m/s or more.

Overview of LES-Based CFD Approach
In the present study, we used the "RIAM-COMPACT" real urban version [14] based on a mixed-time-scale (MTS) subgrid-scale (SGS) model [15] implemented in LES (Equations (1)-(11)). Please refer to the previous papers [16,17] for details of numerical simulation methods.

Overview of LES-Based CFD Approach
In the present study, we used the "RIAM-COMPACT" real urban version [14] based on a mixed-time-scale (MTS) subgrid-scale (SGS) model [15] implemented in LES (Equations (1)-(11)). Please refer to the previous papers [16,17] for details of numerical simulation methods.

Inflow Turbulence Generation Methods and Consideration
Inflow turbulence was artificially generated and imposed as an inflow airflow condition in order to verify the effectiveness of the CFD PD wake model [10] in a real field environment. Here, we explain the methods for artificially generating the inflow turbulence. Figure 15 shows a conceptual diagram of the inflow turbulence generation methods. The recycling method by Lund [18] was simplified and used. The grid resolution in all directions was 2 m (∆x = ∆y = ∆z = 2 m). Therefore, the total number of grids was approximately 41 million points. Concerning the boundary conditions, a power law distribution following N = 10 was assigned to inflow conditions (Equation (12)). Other boundary conditions are the same as in previous paper [14]. The time step was set to ∆t = (1 × 10 -3 ) R/U in , where the Courant number was less than about 1.0. In addition, the variable R represents the rotor radius.
Turbulent fluctuations (fluctuating velocity: u (t, x, y, z), v (t, x, y, z) and w (t, x, y, z)) on a sampling plane were recycled back to the inlet as the inflow boundary condition. After that, time integration was performed until the flow field was fully developed, and turbulent flow data (instantaneous velocity: u(t, x, y, z), v(t, x, y, z) and w(t, x, y, z)) in the sampling plane at a position different from the above were temporarily saved to disk space. These datasets were used as inflow conditions for the main domain that simulated the wind turbine wake. A single two-dimensional (2D) roughness element was placed on the ground in order to obtain turbulent flows in a short period of time. In other words, it could trigger the laminar-turbulent boundary-layer transition. Figure 16 shows the spatial distribution of the instantaneous dimensionless U-velocity in the central cross section of the span, from immediately after the start of the CFD simulation to the time when the flow field was fully developed. From the dimensionless time t* (=t U in /R, R represents the rotor radius) = 0.5 shown in Figure 16a to the dimensionless time 100.0 shown in Figure 16e, the single 2D roughness element described above was installed on the ground. From the dimensionless time 200.0 shown in Figure 16f to the dimensionless time 300.0 shown in Figure 16g, the single 2D roughness element was removed. By observing these figures, it is clear that a separation flow is formed from the single 2D roughness element, moves to the downstream side, and flows again from the inflow boundary cross section. It is observed that the laminar-turbulent boundary-layer transition occurs with the passage of time and that inflow turbulence is generated over the entire computational domain.  Figure 17 shows a rear view of the spatial distribution of the instantaneous U-velocity at the sampling plane. Here, the dimensionless time is t* = 200.0. The position of the swept area of the wind turbine is also described in the figure. It can be seen that the flow field is sufficiently developed in both the y and z directions. Focusing especially on the z direction, it can be seen that the height of the turbulent boundary layer generated in this study is sufficiently higher than the blade top of the wind turbine; that is, the wind turbine is completely buried in the turbulent boundary layer. Figure 18 shows a top view of the spatial distribution of the instantaneous U-velocity at z = 2, 34, 78 and 122 m at the time corresponding to Figure 17. At z = 2 m near the  Figure 17 shows a rear view of the spatial distribution of the instantaneous U-velocity at the sampling plane. Here, the dimensionless time is t* = 200.0. The position of the swept area of the wind turbine is also described in the figure. It can be seen that the flow field is sufficiently developed in both the y and z directions. Focusing especially on the z direction, it can be seen that the height of the turbulent boundary layer generated in this study is sufficiently higher than the blade top of the wind turbine; that is, the wind turbine is completely buried in the turbulent boundary layer. Figure 18 shows a top view of the spatial distribution of the instantaneous U-velocity at z = 2, 34, 78 and 122 m at the time corresponding to Figure 17. At z = 2 m near the ground, it can be seen that very fine vortex structures are formed. A relatively large vortex structures occur, as the altitude from the ground gradually increases. ground, it can be seen that very fine vortex structures are formed. A relatively large vortex structures occur, as the altitude from the ground gradually increases.   Figure 19 shows the mean velocity profile and corresponding three standard deviations evaluated by the sampling plane. Focusing on the mean velocity profile shown in Figure 19a, it can be seen that the numerical results represented by the solid black line reproduce the target profile well (power law distribution with N = 10) shown by the purple circle symbol. On the other hand, in the vertical distribution of the three standard deviations shown in Figure 19b, the velocity fluctuations are in the order of σ u > σ v > σ w in the wind turbine swept area. The figure shows the ratio to the standard deviation in the streamwise direction. Here, the position in which they were evaluated is the height of the wind turbine hub (z = 78 m).
As described above, it was shown that a well-developed turbulent boundary layer can be formed temporally and spatially by using a simplification of the Lund recycling method. In the next section, we explain the wind turbine wake simulation when the turbulent flow data (instantaneous velocity: u (t, x, y, z), v (t, x, y, z) and w (t, x, y, z)) generated here are imposed as inflow conditions.  Figure 19 shows the mean velocity profile and corresponding three standard deviations evaluated by the sampling plane. Focusing on the mean velocity profile shown in Figure 19a, it can be seen that the numerical results represented by the solid black line reproduce the target profile well (power law distribution with N = 10) shown by the purple circle symbol. On the other hand, in the vertical distribution of the three standard deviations shown in Figure 19b, the velocity fluctuations are in the order of σu > σv > σw in the wind turbine swept area. The figure shows the ratio to the standard deviation in the streamwise direction. Here, the position in which they were evaluated is the height of the wind turbine hub (z = 78 m).
As described above, it was shown that a well-developed turbulent boundary layer can be formed temporally and spatially by using a simplification of the Lund recycling method. In the next section, we explain the wind turbine wake simulation when the turbulent flow data (instantaneous velocity: u (t, x, y, z), v (t, x, y, z) and w (t, x, y, z)) generated here are imposed as inflow conditions.

Large Eddy Simulation (LES) of the Wind Turbine Near-Wake Flow Field Using CFD PD Wake Model and Consideration
The numerical simulation method is almost the same as that in the previous section, but in order to reproduce the behavior of the 2 MW-class downwind turbine wake, we used the CFD PD wake model that we recently developed [10]. The effect of the CFD PD wake model, described in Equations (13) and (14), was implemented as an external force in Equation (2).
The only parameter of the CFD PD wake model, CRC, is the resistance coefficient. The purpose of the simulation in this study was to reproduce the Doppler lidar-based wind field measurement in the wind turbine near-wake region (x/D = 3, rotor diameter D = 86

Large Eddy Simulation (LES) of the Wind Turbine Near-Wake Flow Field Using CFD PD Wake Model and Consideration
The numerical simulation method is almost the same as that in the previous section, but in order to reproduce the behavior of the 2 MW-class downwind turbine wake, we used the CFD PD wake model that we recently developed [10]. The effect of the CFD PD wake model, described in Equations (13) and (14), was implemented as an external force in Equation (2).
The only parameter of the CFD PD wake model, C RC , is the resistance coefficient. The purpose of the simulation in this study was to reproduce the Doppler lidar-based wind field measurement in the wind turbine near-wake region (x/D = 3, rotor diameter D = 86 m). For this purpose, we conducted a thorough preliminary investigation of the resistance coefficient (C RC = 1.0 to 10 in 0.5 increments), based on which we adopted 2.5 for the coefficient (C RC ), with a cosine-shaped distribution function in the swept area, as we have already shown in our previous paper [10]. Figure 20 shows the computational domain, computational grid, coordinate system, etc. employed in this study. Of particular note is the fact that the spatial grid resolution in each direction was set at equal intervals of 2 m. As a result, the total number of grid points was about 41 million. In this study, parallel simulations were performed using the supercomputer "SX-Aurora TSUBASA", owned by the RIAM, Kyushu University.
In this study, four patterns of simulations were performed as shown in Table 2. As already explained, the CFD PD wake model was applied to the swept area in all cases. Case 1 considered the effects of nacelle, tower and inflow turbulence. Case 2 considered the effects of nacelle and inflow turbulence. In Case 3, only the effect of inflow turbulence was considered. Case 4 did not consider the effects of nacelle, tower and inflow turbulence. By comparing the results of these simulations, the effect on the airflow characteristics in the wind turbine near-wake region (x/D = 3) was examined in detail.
Energies 2021, 14, x FOR PEER REVIEW 20 of 33 m). For this purpose, we conducted a thorough preliminary investigation of the resistance coefficient (CRC = 1.0 to 10 in 0.5 increments), based on which we adopted 2.5 for the coefficient (CRC), with a cosine-shaped distribution function in the swept area, as we have already shown in our previous paper [10]. Figure 20 shows the computational domain, computational grid, coordinate system, etc. employed in this study. Of particular note is the fact that the spatial grid resolution in each direction was set at equal intervals of 2 m. As a result, the total number of grid points was about 41 million. In this study, parallel simulations were performed using the supercomputer "SX-Aurora TSUBASA", owned by the RIAM, Kyushu University.
In this study, four patterns of simulations were performed as shown in Table 2. As already explained, the CFD PD wake model was applied to the swept area in all cases. Case 1 considered the effects of nacelle, tower and inflow turbulence. Case 2 considered the effects of nacelle and inflow turbulence. In Case 3, only the effect of inflow turbulence was considered. Case 4 did not consider the effects of nacelle, tower and inflow turbulence. By comparing the results of these simulations, the effect on the airflow characteristics in the wind turbine near-wake region (x/D = 3) was examined in detail.     Figure 21 shows a side view of the spatial distribution of the instantaneous U-velocity in Cases 1 to 4. Note first the simulation result shown in Figure 21a, which considers all the inflow turbulence, and the separation flow from the nacelle and tower. In the vicinity of the wind turbine (x/D < 1), the reduction (dimensionless U-velocity < 1.0) of the wind speed by the CFD PD wake model applied to the swept area and the separation flows from the wind turbine nacelle and tower are clearly observed. These flow fields interfere with each other, resulting in the formation of complex wake flows on the downstream side of the wind turbine (1 ≤ x/D). On the other hand, we can see that the inflow turbulence set in the inflow cross section approaches the wind turbine with almost no dissipation of its fluctuations. We can also see that inflow turbulence directly affects the flow field formed near the wind turbine. Additionally, wind speed fluctuations are clearly present even in the far wake region at x = 10 D. Next, note the numerical results of Cases 2 and 3 shown in Figure 21b,c, respectively. In the vicinity of the wind turbine (x/D < 1), there is a slight difference in the flow field formed behind the wind turbine due to the difference in the reproduction of the separated flow from the nacelle and tower. However, the overall behavior of the wake flows at 1 ≤ x/D is not significantly different from the numerical results in Figure 21a. The spatial variation of the wake flows seen in Figure 21a-c cannot be confirmed in Case 4, shown in Figure 21d, in which only the CFD PD wake model is applied to the wind turbine swept area without considering the inflow turbulence, and the separation flows from the nacelle and the tower. However, the overall behavior of wake flows including the wake width in 1 ≤ x/D is similar to that of Figure 21a-c. Figure 22 shows the rear view at x = 3D (near-wake region) corresponding to Figure 21. Even though there are differences in the flow patterns around the wind turbine with and without the consideration of inflow turbulence etc., no significant difference can be confirmed in the wind speed distribution (dimensionless U-velocity < 1.0) within the swept area in all the cases. Figure 23 shows a side view of the spatial distribution of the time-averaged U-velocity in Cases 1 to 4. Figure 24 shows a rear view at x = 3D (near-wake region) corresponding to Figure 23. In this study, the dimensionless time was defined using the radius of the wind turbine blade and the wind speed at the wind turbine hub height. We have taken averages over a sufficiently long time of 100 (dimensionless time). Even though there are differences in the simulation conditions in each case, we found that the qualitative behavior of the wake flows formed behind the wind turbine is extremely similar. In the rear view at x = 3D shown in Figure 24, almost concentric wake regions can be seen to have formed behind the wind turbine in all cases.
Energies 2021, 14, x FOR PEER REVIEW 21 of 3 Figure 21 shows a side view of the spatial distribution of the instantaneous U-velocity in Cases 1 to 4. Note first the simulation result shown in Figure 21a, which considers al the inflow turbulence, and the separation flow from the nacelle and tower. In the vicinity of the wind turbine (x/D < 1), the reduction (dimensionless U-velocity < 1.0) of the wind speed by the CFD PD wake model applied to the swept area and the separation flow from the wind turbine nacelle and tower are clearly observed. These flow fields interfer with each other, resulting in the formation of complex wake flows on the downstream side of the wind turbine (1 ≤ x/D). On the other hand, we can see that the inflow turbulenc set in the inflow cross section approaches the wind turbine with almost no dissipation o its fluctuations. We can also see that inflow turbulence directly affects the flow field formed near the wind turbine. Additionally, wind speed fluctuations are clearly presen even in the far wake region at x = 10 D. Next, note the numerical results of Cases 2 and 3 shown in Figure 21b,c, respectively. In the vicinity of the wind turbine (x/D < 1), there is a slight difference in the flow field formed behind the wind turbine due to the difference in the reproduction of the separated flow from the nacelle and tower. However, the overal behavior of the wake flows at 1 ≤ x/D is not significantly different from the numerica results in Figure 21a. The spatial variation of the wake flows seen in Figure 21a-c canno be confirmed in Case 4, shown in Figure 21d, in which only the CFD PD wake model i applied to the wind turbine swept area without considering the inflow turbulence, and the separation flows from the nacelle and the tower. However, the overall behavior o wake flows including the wake width in 1 ≤ x/D is similar to that of Figure 21a-c. Figure 22 shows the rear view at x = 3D (near-wake region) corresponding to Figur 21. Even though there are differences in the flow patterns around the wind turbine with and without the consideration of inflow turbulence etc., no significant difference can b confirmed in the wind speed distribution (dimensionless U-velocity < 1.0) within th swept area in all the cases. Figure 23 shows a side view of the spatial distribution of the time-averaged U-veloc ity in Cases 1 to 4. Figure 24 shows a rear view at x = 3D (near-wake region) corresponding to Figure 23. In this study, the dimensionless time was defined using the radius of th wind turbine blade and the wind speed at the wind turbine hub height. We have taken averages over a sufficiently long time of 100 (dimensionless time). Even though there ar differences in the simulation conditions in each case, we found that the qualitative behav ior of the wake flows formed behind the wind turbine is extremely similar. In the rea view at x = 3D shown in Figure 24, almost concentric wake regions can be seen to hav formed behind the wind turbine in all cases.
(a)   Figure 25 shows a comparison of the vertical distribution of the time-averaged Uvelocities at x = 3D (near-wake region). As already mentioned, the numerical results shown in Figure 25 were created based on the numerical data extracted from the time-averaged flow fields with a sufficiently long dimensionless time of 100. The figure also shows the range of the wind turbine swept area. The purple symbol is the observation result obtained by the Doppler lidar at a wind speed of 8 m/s (wind turbine hub height velocity), which was the target of the numerical simulations using the CFD PD wake model. The inflow profile set in the present simulations is shown by a solid black line for reference. Note first the numerical results for Case 1 shown by the solid green line-those considering all the inflow turbulence, and the separation flows from the nacelle and tower, in addition to the CFD PD wake model applied in the swept area. It is confirmed that the observation data are underestimated at an altitude of 35 m from the lower end of the wind turbine blade to 78 m at the wind turbine hub height. This is thought to be due to the effect of the separation flow from the tower in the numerical simulation. Next, note the numerical results for Case 2 shown by the solid blue line. In Case 2, the inflow turbulence and the separation flows from the nacelle are reproduced in addition to the CFD PD wake model without considering the separation flows from the tower. Therefore, the discrepancy between the numerical result and the observation data near the wind turbine hub height at 78 m is due to the fact that the separation flows from the nacelle were well reproduced in the numerical results. The results for Cases 3 and 4 are notable; the numerical results for both reproduce the observational data obtained using Doppler lidar well over the entire wind turbine swept area. This suggests the following: in practice, when evaluating the mean wind velocity deficit due to wind turbine wakes, the method of applying the CFD PD wake model in the wind turbine swept area is very effective. Figure 26 shows a comparison of the time-averaged U-velocity space distribution at the wind turbine hub height (78 m) under a wind speed of 8 m/s. For reference, the numerical results of the Park model etc. are also shown in this figure. Please note the position x/D = 3 downstream of the wind turbine No.2, which was the focus of this study. It can also be seen that the Park and Ishihara models give a larger velocity loss than the observation data obtained using the Doppler lidar. As mentioned earlier, it can be seen that both Cases 3 and 4 show good agreement with the observational data. However, the CFD PD wake model used in this study cannot reproduce the mean wind velocity deficit due to wind turbine wakes in the entire area downstream of the wind turbine with one parameter, as stated in our published paper [10]. Therefore, it is necessary to set a resistance coefficient (C RC ) suitable for the distance from the wind turbine that wind energy/wind farm developers want to evaluate. For example, if wind energy/wind farm developers want to know the exact time-averaged velocity deficit at x = 5D, they should set C RC to 9.0 or 10.0, or if they want to know the exact time-averaged velocity deficit at x = 10D, they should set C RC to 5.0 [10].
Finally, when verifying the behavior of wind turbine wake flows, it may be important to understand its location area. To deal with such a situation, we are developing a method for projecting the vertical flow visualization results of the simulation on Google Earth; the results are shown in Appendix B.
In addition, as explained at the beginning of this paper, the purpose of this study was to accurately predict the mutual interference of wind turbine wakes, especially in a largescale wind farm above the ocean consisting of multiple wind turbine groups. Therefore, in order to examine whether the CFD PD wake model can reproduce the mutual interference of the wind turbine wakes, we assumed a large-scale offshore wind farm consisting of multiple wind turbine groups in the nearshore area; the results are shown in Appendix C.

Conclusions
In the present study, the validity of the CFD PD wake model was verified by the wake flows of a 2 MW-class downwind turbine. We first clarified the airflow characteristics of the wind turbine wake flows using a vertically profiling remote sensing wind measurement device (lidar) for the wind turbines installed at the point where the sea is just in front of the wind turbines. A ground-based continuous-wave (CW), conically scanning wind lidar system was used. We paid particular attention to the behavior of the wind

Conclusions
In the present study, the validity of the CFD PD wake model was verified by the wake flows of a 2 MW-class downwind turbine. We first clarified the airflow characteristics of the wind turbine wake flows using a vertically profiling remote sensing wind measurement device (lidar) for the wind turbines installed at the point where the sea is just in front of the wind turbines. A ground-based continuous-wave (CW), conically scanning wind lidar system was used. We paid particular attention to the behavior of the wind

Conclusions
In the present study, the validity of the CFD PD wake model was verified by the wake flows of a 2 MW-class downwind turbine. We first clarified the airflow characteristics of the wind turbine wake flows using a vertically profiling remote sensing wind measurement device (lidar) for the wind turbines installed at the point where the sea is just in front of the wind turbines. A ground-based continuous-wave (CW), conically scanning wind lidar system was used. We paid particular attention to the behavior of the wind turbine near-wakes, which are about three times the diameter of the wind turbine's blade. It was apparent that the influence of the wind turbine wake, that is, the amount of wake loss (wind velocity deficit) is extremely large in the wind speed range of 6.0-10.0 m/s. Furthermore, it was shown that the wake loss (wind velocity deficit) was almost the same at wind speeds of 6.0-10.0 m/s. It was also shown, again, that these results correspond to the distribution of the thrust coefficient of the wind turbine. Currently, many bottom-mounted offshore wind farms are planned for the coastal areas of Japan, where wind speeds of 6.0-10.0 m/s are extremely frequent; in other words, wind speeds of 6.0-10.0 m/s are highly relevant for the realization of bottom-mounted offshore wind farms. Therefore, in the evaluation of the feasibility of bottom-mounted offshore wind farms, strict evaluation at a wind speed of 6.0-10.0 m/s is important; this range was selected in this study for reproducing the behavior of the wind turbine wakes with the CFD PD wake model. Targeting the vertical wind speed distribution in the near-wake region acquired in the lidar, we tuned the model parameters of the CFD PD wake model and considered the numerical results, including those from the conventional engineering wake models. A large parallel computation was performed using a vector supercomputer, showing that, in practice, when evaluating the mean wind velocity deficit due to wind turbine wakes, the method of applying the CFD PD wake model in the wind turbine swept area is very effective.

Data Availability Statement:
The data used to support the findings of this study are available from the corresponding author upon request.

Acknowledgments:
This research was carried out as a joint effort between Toshiba Energy Systems Co., Ltd. and Hitachi Zosen Corporation (Principal Investigator: Takanori Uchida), ongoing since April 2018. We would like to express our gratitude to the people concerned. Company, product and service names appearing in this publication include those that are trademarks or registered trademarks of their respective companies.

Conflicts of Interest:
The authors declare no conflict of interest. Figure A1 shows a schematic of the wake deficit and expansion by the Park model, which is formulated as described below based on that figure. Here, each variable in the equation is as follows:

Appendix A. Arrangement of Park Model Formulation
Uin is the (hub height) inflow wind speed; Uwake is the hub height wind speed downstream of the wind turbine at a distance x; Ct is the thrust coefficient; D is the rotor diameter; k is the wake decay constant.