The Aerodynamic Performance of Horizontal Axis Wind Turbines under Rotation Condition

: The near-surface unsteady incoming ﬂow in the atmospheric boundary layer has a great inﬂuence on the aerodynamic performance of horizontal axis wind turbines. To consider the effect of the rotation of the blade on the aerodynamic state of a wind turbine near the ground, the ﬂuid-structure interaction (FSI) method based on the shear stress transfer (SST) turbulence model is applied to analyze the unsteady aerodynamic interaction characteristics including solving the velocity ﬁeld, pressure ﬁeld, structural response state, variation of deformation


Introduction
Wind energy is one of the most important clean energy to support sustainable energy transformation [1][2][3][4]. As a type of large-scale rotary power generation equipment, the wind turbine frequently works in an erratic airflow environment during the energy conversion [5], such as wind shear, tower shadow effect, rotation effect, atmospheric turbulence, tilt, and yaw imbalance, which may lead to an unstable load of wind turbines [6][7][8][9][10][11]. This directly contributes to wind turbine failure, shorter service lives, lower output power, and more expensive operation and maintenance [12]. At the same time, complex operating conditions will cause separation, stall, and other phenomena in the process of wind turbine blade rotation. This causes instability of the flow field in the three-dimensional space and is accompanied by changes in aerodynamic load. Furthermore, it complicates the operation state of the wind turbine. Therefore, the aerodynamic modeling and the accurate acquisition of unsteady load under the rotating state of the wind turbine are faced with some difficulties [9,13]. To solve this problem, the fluid-structure coupling method is used to extend the computational fluid dynamics characteristics to the structural surface part to study the important method for the operation of the whole wind turbine [14][15][16][17][18]. It provides an effective numerical simulation method to accurately obtain the aerodynamic and structural response state during wind turbine rotation, and evaluate the flow field and aerodynamic performance of wind turbines. Therefore, it is of great significance to study the aerodynamic characteristics of wind turbines under unsteady inflow conditions [5,10,19].
Currently, the numerical methods mainly contain the Blade Element Momentum method (BEM) [20][21][22], Vortex Wake method (VWM) [23,24], and Computational Fluid Dynamics (CFD) [25,26]. These are the three main numerical research techniques used to study the aerodynamic characteristics of wind turbines. Many researchers have studied the aerodynamic performance of wind turbines with different airfoil shapes and sizes by using the blade element momentum method and modification method [21,[27][28][29]. Sang et al. [30] confirmed the feasibility of the new technique through numerical simulation analysis and comparison. They proposed a traditional boundary element correction method for dynamic wind vector changes of wind turbine blades. However, due to the simplification of basic assumptions, the BEM has some limitations to accurately obtain the three-dimensional aerodynamic effects during the operation of wind turbines. Based on potential flow theory, the VWM is a kind of method between BEM and CFD to predict the aerodynamic performance of rotating machinery by calculating the eddy current trajectory. For this method, Dong et al. [31] and Shaler et al. [32] have improved the algorithm on the computational efficiency of the free vortex circulation model. The increase in the size of the wind turbine makes the blades more flexible and increases the interaction deformation, which is followed by unstable load fluctuations. These fluctuations damage the structural stability and aerodynamic performance of the blades and affect the stable operation of the wind turbine. The CFD method does not simplify the actual flow field too much but obtains the complete flow field information by solving the fluid control equation numerically. It shows its powerful ability to calculate the aerodynamic performance of wind turbines and analyze the flow mechanism. Therefore, it is an important means for the aerodynamic research of wind turbines. The computational fluid mechanics method is more and more favored by scholars.
For the aerodynamic performance research of wind turbines with unsteady influence, the current research mainly focuses on the influence of the wind turbine blade structure part. The research on the influence of the overall structure of wind turbines is more neglected. The influence of the external field environment is mostly considered under normal incoming wind conditions, so the physical characteristics of the full-size wind turbine flow field and the fluid-structure interaction mechanism cannot be accurately analyzed.
In this paper, the unsteady conditions under the joint action of the uniform incoming flow and three-dimensional rotational effects are considered. ANSYS platform fluidstructure interaction analysis is combined. The results of the steady-state calculation of the flow field through the structural field are received. The interaction surface pressure, velocity, and other data are transferred to study the aerodynamic interaction characteristics of the solid structure of a wind turbine under the action of the fluid. Section 2 introduces the research method and mathematical model. It summarizes the superiority of CFD computational fluid dynamics simulations. Also, it solves the CFD-CSD equations using weakly coupled separation. Section 3 describes the characteristics and material properties of the research object. It sets up the wind turbine geometry models for different working conditions. Moreover, it establishes the coupled computational models for the internal and external flow field environments and the full-size wind turbine geometry. Section 4 presents all the analysis results and explains the related physical phenomena and reasons.

Governing Equations
The weakly coupled separated solution method is used in conjunction with the numerical simulation software solution module. To transfer the computational data of the fluid-domain fluid dynamics (CFD) equations and the solid-domain structural dynamics (CSD) equations in a predetermined order, realizing the fluid-structure interaction at each Sustainability 2023, 15, 12553 3 of 15 time step. The air in this study can be thought of as being incompressible because the fluid velocity is low and far below Mach 0.3 [33]. The gas density is taken to be 1.1691 kg/m 3 at 25 • C, one standard atmosphere. The assumed constant viscosity is 1.84 × 10 −5 kg/ms −1 . The whole flow area is solved for the incompressible, unsteady Reynolds-averaged Navier-Stokes (URANS) equation. To increase accuracy, a solver using a linked pressure basis and a second-order implicit formulation is selected. The solution variables are solved using a second-order windward discretization scheme.

Fluid Governing Equations
Mass conservation equation: Momentum conservation equation: where ρ f is the fluid density, v is the fluid velocity vector, f f is the body force vector, τ f is the shear stress tensor, the expression is: where p is fluid pressure, µ is dynamic viscosity, I is the identity matrix, and e is the velocity stress tensor, the expression is:

Structure Governing Equations
The conservation equation of the solid part system is derived from Newton's second law as follows: where ρ s is structure density, .. d s is a local structure domain acceleration vector, σ s is the Cauchy stress tensor, and f s is the volume force vector.

Basic Conservation Equation of Fluid-Structure Interaction
The conservation of stress and displacement at the interaction junction is calculated. The following equations are satisfied: where τ is stress, d is displacement, n is the direction cosine, the subscript f represents the fluid, and the subscript s stands for solid.

Turbulence Modeling
The investigation of the aerodynamic connection of horizontal axis wind turbines employs a hybrid k-ω/k-ε shear stress transfer (SST) turbulence model [34]. The k-ω SST model is a two-equation hybrid model applicable to the k-ω/k-ε turbulence model interconverted to forecast far-field flow. It can take turbulent shear transport into account and successfully predict flow transport under various pressure gradient situations [35]. For exterior areas and free shear flow, the k-ε model is employed [36]. This model has been widely utilized by academics to examine the aerodynamic functionality of wind turbines with good convergence and accuracy [5,21,37,38]. The accuracy of the convergence of the residuals for flow continuity, velocity, turbulent kinetic energy, and specific dissipation rate was adjusted to 1 × 10 −6 and 1 × 10 −3 . Indicator residual monitoring was utilized to determine the convergence of the computational findings. Pressure, momentum, turbulent kinetic energy, specific dissipation rate, and specific energy were all roughly discretized spatially to second order. Here are the details of the SST model: where k is turbulent kinetic energy, u i is the directional speed, x i is the directional displacement, µ t is eddy viscosity,P k is turbulence kinetic energy transmitted forward, ω is the specific dissipation rate, S is the invariant measure of strain rate, σ k , β * , σ ω , σ ω2 are model-related constants. The mixed function F 1 is defined as: where CD kω = max 2ρσ ω2 , 10 −10 is the cross-diffusion term, y is the distance to the closest wall, v is the kinematic viscosity, F 1 is one inside the boundary layer (k-ω model), and both values are zero away from the surface.

Geometric Model
The NACA44XX airfoil coordinates are selected in the Profili program. The leaf element coordinates of the blade section are added to the airfoil coordinates. To create the fundamental framework of the blade, the leaf element curves were scaled and rotated around the angle of attack of the airfoil. The laminate surface was used to simulate the leaf element coordinates. The model of the 1.5 MW NREL horizontal axis wind turbine is referenced in the geometric structure. The engine room is a cube structure measuring 0.9 m × 0.7 m × 0.525 m thanks to a simplified design, as shown in Figure 1. The diameter of the wind wheel is 6.5 m, the height of the tower barrel is 8.75 m, the top section diameter of the tower barrel is 0.25 m, and the bottom section diameter is 0.45 m.

Computational Domain Configuration
To replicate the wind turbine in its rotating state. The moving cylindrical rotating inner domain enclosing the entire rotor of the wind turbine is joined with the stationary rectangular outer domain containing the engine room and tower barrel. The diameter and thickness of the inner cylinder are 6.8 m and 0.6 m, respectively. By using a Boolean operation, the blade and inner rotation domains are produced. The middle tower and inner

Computational Domain Configuration
To replicate the wind turbine in its rotating state. The moving cylindrical rotating inner domain enclosing the entire rotor of the wind turbine is joined with the stationary rectangular outer domain containing the engine room and tower barrel. The diameter and thickness of the inner cylinder are 6.8 m and 0.6 m, respectively. By using a Boolean operation, the blade and inner rotation domains are produced. The middle tower and inner basin were removed. In order to simulate the unsteady influence of uniform wind speed and reduce the influence of the boundary environment, the method of digging out the tower cylinder in the outer basin and the inner basin is adopted. The entire calculation area is 40 m × 30 m × 21 m, as shown in Figure 2. The two regions moving relative to each other share a sliding interface that acts as a coupling surface for data transmission.
To replicate the wind turbine in its rotating state. The moving cylindrical rotating inner domain enclosing the entire rotor of the wind turbine is joined with the stationary rectangular outer domain containing the engine room and tower barrel. The diameter and thickness of the inner cylinder are 6.8 m and 0.6 m, respectively. By using a Boolean operation, the blade and inner rotation domains are produced. The middle tower and inner basin were removed. In order to simulate the unsteady influence of uniform wind speed and reduce the influence of the boundary environment, the method of digging out the tower cylinder in the outer basin and the inner basin is adopted. The entire calculation area is 40 m × 30 m × 21 m, as shown in Figure 2. The two regions moving relative to each other share a sliding interface that acts as a coupling surface for data transmission.
Before CFD calculation analysis, grid independence verification was carried out. Three different grid sizes were used in the grid test. Grid details and results are shown in Table 1. Verify that the grid is divided into fine, medium, and coarse, and the number of grid units increases in equal proportion, and the overall grid growth rate between adjacent grids is 25%. The three grids are simulated and calculated under standard working conditions, and the final output power results are compared. It can be seen that even on the coarse grid, the power difference between the coarse grid and the medium grid is 4.32%, and the power difference between the medium grid and the fine grid is only 1.42%. Therefore, the accuracy of the medium grid is sufficient.  Before CFD calculation analysis, grid independence verification was carried out. Three different grid sizes were used in the grid test. Grid details and results are shown in Table 1. Verify that the grid is divided into fine, medium, and coarse, and the number of grid units increases in equal proportion, and the overall grid growth rate between adjacent grids is 25%. The three grids are simulated and calculated under standard working conditions, and the final output power results are compared. It can be seen that even on the coarse grid, the power difference between the coarse grid and the medium grid is 4.32%, and the power difference between the medium grid and the fine grid is only 1.42%. Therefore, the accuracy of the medium grid is sufficient.

CFD Mesh and Boundary Conditions
The flow field setting is realized by the FLUENT module in ANSYS Workbench. The flow field domain calculation inhibits the structural field part of the overall calculation domain. The blade surface is a rotating wall without slip shear and is also the interface of the fluid domain and the structural domain that transfer data between the two solvers. The area is discretized with unstructured meshes. Finer meshes and additional expansion layers are added to the blade surface to improve resolution. In meshing, the size of the symmetrical wall and ground elements was set at 2.719 m. The size of the volume mesh was set at 0.1 m in the internal rotation domain. The curvature capture division method was adopted, and the normal angle of curvature was defined as 18 • . The mesh division of The flow field setting is realized by the FLUENT module in ANSYS Workbench. The flow field domain calculation inhibits the structural field part of the overall calculation domain. The blade surface is a rotating wall without slip shear and is also the interface of the fluid domain and the structural domain that transfer data between the two solvers. The area is discretized with unstructured meshes. Finer meshes and additional expansion layers are added to the blade surface to improve resolution. In meshing, the size of the symmetrical wall and ground elements was set at 2.719 m. The size of the volume mesh was set at 0.1 m in the internal rotation domain. The curvature capture division method was adopted, and the normal angle of curvature was defined as 18°. The mesh division of the flow field is shown in Figure 3. The details of the gridʹs local magnification are marked in red. The upstream velocity inlet boundary param was set at free flow wind speed at a distance of 10 m from the entire machine. The initial standard flow rate is 12 m/s, turbulence kinetic energy intensity at 0.05, and turbulence viscosity ratio at 10. A standard atmosphere of pressure is present at the downstream pressure outlet border, which is 30 m from the machine. The fixed symmetric wall surface, ground surface, and blade surface of the inner rotating domain in the outer domain were given non-slip bounds. Frame motion was applied to create the rotation domain to rotate around the Z-axis relative to the stationary domain. As shown in Figure 4, the connection between the wind turbine and the engine room is the rotation axis. The center point of the connection surface is the initial coordinate of the rotation axis. The initial coordinates of the rotation axis were (−4.9567 × 10 −3 , 2.541, 0.7535), and the rotation speed was defined. Since the blade and the rotation domain were in relative motion during the simulation process, a wall motion mode similar to that in the rotation domain was adopted for the blade interface. The moving wall was selected, the movement type was rotation, and the relative speed was 0. The upstream velocity inlet boundary param was set at free flow wind speed at a distance of 10 m from the entire machine. The initial standard flow rate is 12 m/s, turbulence kinetic energy intensity at 0.05, and turbulence viscosity ratio at 10. A standard atmosphere of pressure is present at the downstream pressure outlet border, which is 30 m from the machine. The fixed symmetric wall surface, ground surface, and blade surface of the inner rotating domain in the outer domain were given non-slip bounds. Frame motion was applied to create the rotation domain to rotate around the Z-axis relative to the stationary domain. As shown in Figure 4, the connection between the wind turbine and the engine room is the rotation axis. The center point of the connection surface is the initial coordinate of the rotation axis. The initial coordinates of the rotation axis were (−4.9567 × 10 −3 , 2.541, 0.7535), and the rotation speed was defined. Since the blade and the rotation domain were in relative motion during the simulation process, a wall motion mode similar to that in the rotation domain was adopted for the blade interface. The moving wall was selected, the movement type was rotation, and the relative speed was 0.

Structural Model Setup
The calculation and setting of the structural field consider the pressure data from the numerical calculation of the flow field and ignore the volume part of the flow field. Before mesh division, material properties are assigned to the wind wheel and tower barrel, respectively, as shown in Table 2.

Structural Model Setup
The calculation and setting of the structural field consider the pressure data from the numerical calculation of the flow field and ignore the volume part of the flow field. Before mesh division, material properties are assigned to the wind wheel and tower barrel, respectively, as shown in Table 2. During blade mesh division, to avoid element distortion at the front and rear flange Angle tips during mesh division, surface mesh division units were embedded, and blade and tower barrel element sizes were set at 0.02 m and 0.05 m, respectively, as shown in Figure 4a. Since the blade surface is a coupled boundary, the aerodynamic force and structural displacement are exchanged, and the large deflection assumption is adopted. The centrifugal force takes into account the rotational speed of the blade structure and defines the rotation mode of the wind wheel. To solve the setting, add the rotation axis and amplitude value of the wind turbine, select the bottom of the wind turbine support, and apply a fixed constraint. In the solution of fluid-structure interaction and structure, the blade surface was taken as the interaction interface, and the flow field pressure data was imported. The maximum pressure value was 2448 Pa, and the minimum pressure value was −3559.1 Pa, as shown in Figure 4b.

Results and Discussion
To evaluate the convergence of CFD analysis, the residual, velocity, turbulent kinetic energy, and specific dissipation rate were monitored, 2000 iteration steps were set, the standard initialization method was adopted, and the initial value was calculated from the inlet boundary. The convergence was proved by determining that the residual value was less than 1 × 10 −6 or the net mass percentage was less than 0.1% [39,40]. The result shows that the net export mass flow difference is much less than 0.1% of the total import flow. Figure 5 shows the velocity distributions at the wind turbine of z = 1 m plane and the blade, respectively. As plotted in Figure 5a, the global velocity has a range from 0 m/s to 89.25 m/s, and there is a velocity expansion along the blade, in which the velocity increases rapidly from the blade root to the blade tip, and the rotating annular flow pattern can be observed. Considering the blocking effect of the blade structure on the airflow, the area where the blade is rotated forms wind eddy currents, leading to aerodynamic changes. Such eddy currents are generated immediately and disappear quickly with the rotation of the rotor. The vortex alternately and periodically occurs and falls off, dissipating a large amount of energy in the mainstream, so the wind speed in the flow field area where the blade is rotated is lower than that in other areas, which is the reason for the trailing phenomenon of the blade. Figure 5b shows that the hub region is in a relatively static state, and the speed starts to show a linear growth trend at the blade root and gradually develops to the tip at the mid-span. Such eddy currents are generated immediately and disappear quickly with the rotation of the rotor. The vortex alternately and periodically occurs and falls off, dissipating a large amount of energy in the mainstream, so the wind speed in the flow field area where the blade is rotated is lower than that in other areas, which is the reason for the trailing phenomenon of the blade. Figure 5b shows that the hub region is in a relatively static state, and the speed starts to show a linear growth trend at the blade root and gradually develops to the tip at the mid-span.  Figure 6 shows the wind turbine velocity distribution contours at x = 0 m and y = 2.5 m planes. From Figure 6a, it can be seen that the rotation domain is blocked under the incoming wind speed. The velocity variation in the near wind turbine region is concentrated, resulting in wake drag in the downstream region of the wind turbine, and the range of velocity variation expands. Without the barrier of the structure, the airflow exchanges to a steady state. The wind speed of the far wind area along the inflow direction gradually decreases and finally reaches the initial speed. Figure 6b shows that as the wind speed decreases in the flow field created by the wind turbine, the influence range of the airflow disturbance grows steadily bigger. The velocity gradient gradually decreases along the wake region and eventually dissipates in the atmosphere to return to normal wind levels. Due to changes in the adjacent flow field brought on by the tower barrel effect, the distribution position of the tower barrel and engine room, along with the geometric shape of the blades, induce velocity loss in the near wake area, lowering wind speed behind the tower.   Figure 6 shows the wind turbine velocity distribution contours at x = 0 m and y = 2.5 m planes. From Figure 6a, it can be seen that the rotation domain is blocked under the incoming wind speed. The velocity variation in the near wind turbine region is concentrated, resulting in wake drag in the downstream region of the wind turbine, and the range of velocity variation expands. Without the barrier of the structure, the airflow exchanges to a steady state. The wind speed of the far wind area along the inflow direction gradually decreases and finally reaches the initial speed. Figure 6b shows that as the wind speed decreases in the flow field created by the wind turbine, the influence range of the airflow disturbance grows steadily bigger. The velocity gradient gradually decreases along the wake region and eventually dissipates in the atmosphere to return to normal wind levels. Due to changes in the adjacent flow field brought on by the tower barrel effect, the distribution position of the tower barrel and engine room, along with the geometric shape of the blades, induce velocity loss in the near wake area, lowering wind speed behind the tower.

Aerodynamics Distribution
Such eddy currents are generated immediately and disappear quickly with the rotation of the rotor. The vortex alternately and periodically occurs and falls off, dissipating a large amount of energy in the mainstream, so the wind speed in the flow field area where the blade is rotated is lower than that in other areas, which is the reason for the trailing phenomenon of the blade. Figure 5b shows that the hub region is in a relatively static state, and the speed starts to show a linear growth trend at the blade root and gradually develops to the tip at the mid-span.  Figure 6 shows the wind turbine velocity distribution contours at x = 0 m and y = 2.5 m planes. From Figure 6a, it can be seen that the rotation domain is blocked under the incoming wind speed. The velocity variation in the near wind turbine region is concentrated, resulting in wake drag in the downstream region of the wind turbine, and the range of velocity variation expands. Without the barrier of the structure, the airflow exchanges to a steady state. The wind speed of the far wind area along the inflow direction gradually decreases and finally reaches the initial speed. Figure 6b shows that as the wind speed decreases in the flow field created by the wind turbine, the influence range of the airflow disturbance grows steadily bigger. The velocity gradient gradually decreases along the wake region and eventually dissipates in the atmosphere to return to normal wind levels. Due to changes in the adjacent flow field brought on by the tower barrel effect, the distribution position of the tower barrel and engine room, along with the geometric shape of the blades, induce velocity loss in the near wake area, lowering wind speed behind the tower.   Figure 7 shows the pressure distribution contours of the wind turbine rotation domain at the profile z = 1 m and the blade. The pressure of the rotation domain at the profile z = 1 m is shown in Figure 7a. As the wind wheel is shearing the air, negative pressure starts to develop on the lower surface of the blade along the direction of the rotation. The pressure across the middle to the tip of the blade begins to increase. Considering that at the same angular velocity, the linear velocity increases as it stretches along the chord length of the blade. The pressure difference between the upper and lower surfaces of the blade leads to an increase in the airflow supplement exchange speed. It can be observed that the airflow disturbance around the rotation region weakens and tends to be stable with the increase of the radial distance.
The pressure of the blade is shown in Figure 7b. The negative pressure on the upwind surface of the blade ranges from −3160 Pa to −217 Pa. In the middle of the blade, clear negative pressure zones start to form, and negative pressure is concentrated on the blade surface. Additionally, as the blades are propelled by the normal resultant force, the pressure is evenly distributed on the windward. The leeward sides are no longer completely attached to the blade surface due to the action of centrifugal force but instead, eventually dissipate through energy conversion. It complies with the law of energy conservation and also proves the rationality of the aerodynamic characteristics of the wind turbine.
blade leads to an increase in the airflow supplement exchange speed. It can be observed that the airflow disturbance around the rotation region weakens and tends to be stable with the increase of the radial distance.
The pressure of the blade is shown in Figure 7b. The negative pressure on the upwind surface of the blade ranges from −3160 Pa to −217 Pa. In the middle of the blade, clear negative pressure zones start to form, and negative pressure is concentrated on the blade surface. Additionally, as the blades are propelled by the normal resultant force, the pressure is evenly distributed on the windward. The leeward sides are no longer completely attached to the blade surface due to the action of centrifugal force but instead, eventually dissipate through energy conversion. It complies with the law of energy conservation and also proves the rationality of the aerodynamic characteristics of the wind turbine.  The wind shear mainly changes the stress gradient on the blade surface. Under the action of clockwise torque, the surface stress on the windward side is greater than that on the leeward side. The minimum stress on the leeward side is 70.058 Pa. The stress concentration area of the wind wheel is located in the area where the blade root is connected to the hub. The stress concentration area is mainly on the mid-span surface of the blade caused by the forward flow shear, and the stress peak value is 5.9239 MPa, which is far less than the designed strength of the blade material.

Structural Response
The strain contours are shown in Figure 8b. The analysis reveals that stress and strain exhibit a linear association because the distribution of the stress-strain contours is essentially the same. The large strain concentration area is located in the middle and upper part of the blade and the middle and lower part of the tower, and the maximum strain is located in the middle and upper part of the blade and reaches 8.6116 × 10 −5 m/m, so the design strength of the middle and upper part of the blade should be improved to prevent   Figure 9 shows the deformation curve of wind turbine blades during shutdown and rotation. As shown in the lower right corner of Figure 9a, we numbered each blade. The starting point of the displacement curve does not reflect the zero displacement, but the displacement of the blade at the blade root.

Structural Response
The starting points of the three curves represent the displacement of the blade root and are not zero because the displacement on the blade is extracted on the basis of the deformation of the tower barrel, so the initial value corresponds to the displacement value at the top of the tower barrel. The endpoints of the three curves represent the displacement on the tip of the blade. It can be seen from the figure as the position with the largest deformation on the entire blade. The comparative analysis shows that the total displacement curves of the windward side of blade 1 and blade 2 both increase in the positive exponential form under these two conditions. The blade tip deformation is the largest in the shutdown state, and the transverse displacement of the blades in the windward area is basically the same. The same rule applies to rotation conditions. The wind shear mainly changes the stress gradient on the blade surface. Under the action of clockwise torque, the surface stress on the windward side is greater than that on the leeward side. The minimum stress on the leeward side is 70.058 Pa. The stress concentration area of the wind wheel is located in the area where the blade root is connected to the hub. The stress concentration area is mainly on the mid-span surface of the blade caused by the forward flow shear, and the stress peak value is 5.9239 MPa, which is far less than the designed strength of the blade material.
The strain contours are shown in Figure 8b. The analysis reveals that stress and strain exhibit a linear association because the distribution of the stress-strain contours is essentially the same. The large strain concentration area is located in the middle and upper part of the blade and the middle and lower part of the tower, and the maximum strain is located in the middle and upper part of the blade and reaches 8.6116 × 10 −5 m/m, so the design strength of the middle and upper part of the blade should be improved to prevent fracture. Therefore, the design strength of the upper part of the blade should be improved to prevent breaking. Figure 9 shows the deformation curve of wind turbine blades during shutdown and rotation. As shown in the lower right corner of Figure 9a, we numbered each blade. The starting point of the displacement curve does not reflect the zero displacement, but the displacement of the blade at the blade root. The total deformation change curve of blade 3 in the downwind area of the wind turbine under different speeds is shown in Figure 10. The tower shadow effect affects the structural deformation at standard operating conditions (200 rpm), and the influence of the tower shadow effect on air interference is weakened at lower rotational speeds, resulting in the deformation mutation of the blades in the downwind area being advanced in the structural response. The change of blade deformation at lower rotation speeds (110 rpm, 140 rpm, 170 rpm) is also shown in the figure, and a comparative analysis is made with that at standard speed. The deformation of the blades in the downwind region also has a similar mutation as that under the standard speed, resulting in an advance in the structural response of the blades in the downwind region. The most intuitive representation in the figure is the advance of the deformation mutation point in the upward diffusion process of the blade. At lower rotational speeds, the influence of the tower shadow effect on air interference is enhanced, making it more difficult to get rid of the influence of tower The starting points of the three curves represent the displacement of the blade root and are not zero because the displacement on the blade is extracted on the basis of the deformation of the tower barrel, so the initial value corresponds to the displacement value at the top of the tower barrel. The endpoints of the three curves represent the displacement on the tip of the blade. It can be seen from the figure as the position with the largest deformation on the entire blade. The comparative analysis shows that the total displacement curves of the windward side of blade 1 and blade 2 both increase in the positive exponential form under these two conditions. The blade tip deformation is the largest in the shutdown state, and the transverse displacement of the blades in the windward area is basically the same. The same rule applies to rotation conditions. Under rotating conditions, the response frequency of the blade structure becomes larger, and the overall displacement increases significantly compared with that when the blade is stopped, and the peak increment in the upwind area reaches 2.1437 mm and 0.8674 mm, respectively. It shows that the aerodynamic force and centrifugal force of the wind turbine significantly affect the deformation of the blades during operation. Similarly, the displacement of blade 3 in the downwind area decreases at first and then increases.
Unlike when the blade is stopped, the displacement mutation appears in the second half span and then recovers to the initial value. This is because when the blade overlaps with the tower, the airflow disturbance in the area around the structure is enhanced, and mutual interference leads to the deformation and fluctuation of the blade. When the blade skimmed over the tower, the airflow interference in the interference region was weakened, and the response frequency of the structure began to increase. At this time, the structural response state of blade 3 recovers rapidly until the lower half span begins to show a changing trend similar to the blade deformation in the windward region. It is proved that the azimuth angle has a great influence on the dynamic response of wind turbine structures.
The total deformation change curve of blade 3 in the downwind area of the wind turbine under different speeds is shown in Figure 10. The tower shadow effect affects the structural deformation at standard operating conditions (200 rpm), and the influence of the tower shadow effect on air interference is weakened at lower rotational speeds, resulting in the deformation mutation of the blades in the downwind area being advanced in the structural response. The change of blade deformation at lower rotation speeds (110 rpm, 140 rpm, 170 rpm) is also shown in the figure, and a comparative analysis is made with that at standard speed. The deformation of the blades in the downwind region also has a similar mutation as that under the standard speed, resulting in an advance in the structural response of the blades in the downwind region. The most intuitive representation in the figure is the advance of the deformation mutation point in the upward diffusion process of the blade. At lower rotational speeds, the influence of the tower shadow effect on air interference is enhanced, making it more difficult to get rid of the influence of tower shadow downwind. When the rotational speed increases, the deformation of the blades in the downwind region develops towards the blades in the upwind region, resulting in the mutation point moving toward the blade tip.

Torque and Output Power
To study the output power changes of wind turbines with different wind speeds under uniform inflow conditions, the inflow conditions were set at 5 m/s, 8 m/s, 10 m/s, 12 m/s, 15 m/s, 18 m/s, and 20 m/s. The above inflow conditions were simulated and calculated, respectively, and the wind turbine torque value M under different working conditions was obtained through data processing. The wind turbine output power P was obtained by using the following formula, and then the output data was analyzed. 15 m/s, 18 m/s, and 20 m/s. The above inflow conditions were simulated and calculated, respectively, and the wind turbine torque value M under different working conditions was obtained through data processing. The wind turbine output power P was obtained by using the following formula, and then the output data was analyzed.
where, ω is the angular velocity, ω = RPM·π/30, rad/s. Figure 11 shows the variation curve of wind turbine output power under different incoming flow conditions.

Torque and Output Power
To study the output power changes of wind turbines with different wind speeds under uniform inflow conditions, the inflow conditions were set at 5 m/s, 8 m/s, 10 m/s, 12 m/s, 15 m/s, 18 m/s, and 20 m/s. The above inflow conditions were simulated and calculated, respectively, and the wind turbine torque value M under different working conditions was obtained through data processing. The wind turbine output power P was obtained by using the following formula, and then the output data was analyzed.

P Mω =
where, ω is the angular velocity, , rad/s. Figure 11 shows the variation curve of wind turbine output power under different incoming flow conditions. It can be concluded from the figure above that the output power of the wind turbine presents approximately linear changes under different uniform inlet flows of the wind field. When the inlet wind speed increases from 5 m/s to 20 m/s, the output power of the wind turbine increases from 68.33 kW to 84.33 kW, respectively. In the process of increasing output power, considering that the torque of wind turbines increases with the increase It can be concluded from the figure above that the output power of the wind turbine presents approximately linear changes under different uniform inlet flows of the wind field. When the inlet wind speed increases from 5 m/s to 20 m/s, the output power of the wind turbine increases from 68.33 kW to 84.33 kW, respectively. In the process of increasing output power, considering that the torque of wind turbines increases with the increase of windward speed. The aerodynamic pressure captured by wind turbine blades under the larger inlet wind speed is also higher. Due to the structural characteristics of the blades, under the action of the flow field, blade surface load changes blade shape param, and the lift-drag ratio of wind turbines increases. This results in better performance of the blades and an increase in output power.
The wind energy captured in the operation of wind turbines is provided to the rotor to generate torque. The centrifugal force diffuses the flow field data to the surrounding area, enhances the interaction of wind speed, and increases the central torque of wind turbines. Such changes in radial flow velocity and pressure are more significant, and ultimately, the output power is generated in essence. The mechanical properties of wind turbine blade speed and pressure are different under different wind conditions, which provides topic guidance for the study of wind turbine overall aerodynamic performance or aerodynamics.
The output power increases linearly with the incoming flow speed, and there is a strong correlation between the two. It is conceivable that this variation may represent the general characteristics of wind turbines at real-time wind speeds. In addition, Table 3 and Figure 12 show three previous studies [26,38,41] for power monitoring under changes in wind speed under different model sizes and turbulence models. The parameter values and results of similar CFD models under different studies are shown in Table 3, and the variation trend of the standardized power of the literature model with the standardization speed is shown in Figure 12. The velocity curves of the standardized calculation results of the three similar studies are also similar to those in this paper.

Conclusions
The finite element software ANSYS is used to study the downstream solid interaction aerodynamic characteristics of wind turbines under the action of the flow field, and the following conclusions are drawn.
(1) The fluid-structure interaction wind field model of the wind turbine is established and verified in light of the complexity of the aerodynamic environment. Under the influence of inflow wind and rotation in the flow field, the airflow disturbance closest to the wind turbine is the most important. Under the interference of the tower shadow effect, the wind speed in the near wake zone and the far wake zone decreases to varying degrees along the incoming wind direction. (2) While the centrifugal force of the normal force that drives the pressure around the wind turbine has a greater impact and the pressure around the wind turbine changes primarily along the radial direction with the weakening of the airflow disturbance, the windward side of the blade first experiences the action of airflow, and the pressure changes in a wide range. The maximum stress in the structural reaction of the tower is eight times greater at its top than it is at its bottom. The peak blade stress reaches 5.92 MPa.  Wang et al. [26] and Liu et al. [41] are consistent with the variation range of incoming wind speed in this paper, but the CFD turbulence model adopted by them is different. The standardized power variation of the two models has a good correlation with the standardized speed, and the fitting correlation is relatively consistent with the CFD calculation results of the k-ω SST model in this paper. While the wind speed range of Mo et al. [38] is large, although the R 2 value is small, the results still confirm that the power has good growth at lower wind speeds. There is also a basically consistent trend of change. Compared with the CFD results of Wang et al. [26], Mo et al. [38], and Liu et al. [41], it can be seen that this study has a good similar variation rule with the literature data. It confirms the reliability of the CFD results and model.

Conclusions
The finite element software ANSYS is used to study the downstream solid interaction aerodynamic characteristics of wind turbines under the action of the flow field, and the following conclusions are drawn.
(1) The fluid-structure interaction wind field model of the wind turbine is established and verified in light of the complexity of the aerodynamic environment. Under the influence of inflow wind and rotation in the flow field, the airflow disturbance closest to the wind turbine is the most important. Under the interference of the tower shadow effect, the wind speed in the near wake zone and the far wake zone decreases to varying degrees along the incoming wind direction. (2) While the centrifugal force of the normal force that drives the pressure around the wind turbine has a greater impact and the pressure around the wind turbine changes primarily along the radial direction with the weakening of the airflow disturbance, the windward side of the blade first experiences the action of airflow, and the pressure changes in a wide range. The maximum stress in the structural reaction of the tower is eight times greater at its top than it is at its bottom. The peak blade stress reaches 5.92 MPa. (3) The change in wind turbine blade displacement in the upwind direction under standard (12 m/s and 200 rpm) and shutdown conditions is the same and grows exponentially. Under the action of rotation, the transverse deformation of blade 3 occurred in 73%, and the transverse deformation of downwind blade 3 was almost unchanged at 20%. With the increase of wind turbine speed, the abrupt deformation point of the blade in the downwind area moves towards the blade tip. (4) The aerodynamic impact of the wind turbine is improved by increasing the entering wind speed, catching more wind energy, producing a considerable bending moment and power, and increasing output power by a net 20 kW.