Experimental and Numerical Analysis of the Effect of Vortex Generator Height on Vortex Characteristics and Airfoil Aerodynamic Performance

Xinkai Li 1,2,3,* , Ke Yang 1,2,3,4 and Xiaodong Wang 5 1 Institute of Engineering Thermophysics, Chinese Academy of Sciences, Beijing 100190, China; yangke@iet.cn 2 Key Laboratory of Wind Energy Utilization, Chinese Academy of Sciences, Beijing 100190, China 3 Dalian National Laboratory for Clean Energy, CAS, Beijing 100190, China 4 Graduate School, University of Chinese Academy of Sciences, Beijing 100049, China 5 College of Energy, Power and Mechanical Engineering, North China Electric Power University, Beijing 102206, China; wangxd@ncepu.edu.cn * Correspondence: lixinkai@iet.cn; Tel.: +86-010-8254-3038


Introduction
According to the latest data from the Global Wind Energy Council [1], as of 2017, the total installed capacity of the global wind power market has reached 539.6 GW, and wind power has become the third largest power source in the world.However, during the operation of wind turbines, the blade root is prone to flow separation, which decreases the power of wind turbines [2][3][4][5][6].Many novel methods have been proposed for improving the performance of wind turbine blades in recent years, such as vortex generators (VGs) [7], microflaps [8], microtabs [4], blowing and suction [9], synthetic jets [10], flexible wall [11], and plasma actuators [12].A comparison of these methods was performed by Barlas [13] and Johnson [14] to illustrate their efficiency.As reported by Lin [7] and Wang [15], vortex generators are among the most effective devices for improving the aerodynamic performance of blades.
Taylor [16] introduced vortex generators in 1947, and they were originally implemented in the field of flow control for aircraft wings.The principle underlying VGs' control of flow separation is that when a fluid flows over VGs, a concentrated vorticity is generated.Under the action of the concentrated vortices, the high-energy fluid outside the boundary layer will be drawn into the bottom of the boundary layer, and energy will be injected into the low-energy fluid in the boundary layer, thus delaying the separation of the boundary layer [17,18].Abdollah [19] used a numerical simulation method to verify whether installing VGs on wind turbine blades increases the efficiency of the wind turbine.Because of their simple structure and convenient installation, VGs have been widely utilized in the field of flow control [20][21][22] and heat transfer [23,24].Lin [7] studied the mechanism by which VGs control the flow separation in the boundary layer on a flat plate by using an experimental method.Godard [25][26][27] studied the effect of a triangular VG arrangement on the control of boundary-layer separation with an inverse pressure gradient by an experimental method.Urkiola et al. [22] described the vortex characteristics of a single VG on the basis of a Computational Fluid Dynamics (CFD) method.Velte et al. [28] measured the flow field of a rectangular VG on a flat plate; in their experiment, the VG height was the thickness of the local flat plate boundary layer.The relationships between the installation angle of the VGs and the radius of the vortex core, vortex circulation, and the induced velocity of the vortex core were studied by varying the installation angle of the VGs in the range of 20 • ≤ β ≤ 40 • .The results demonstrated that these quantities change linearly with the change in the installation angle of the VGs.It was also found that, although the VG is located in the turbulent boundary layer, the induced vortices are stable.Angele and Grewe [29] studied the control method of turbulent boundary-layer separation, and the transient eddy behavior of the turbulent boundary layer was analyzed in detail.To research airfoil flow separation control, Sørensen et al. [30] described a CFD procedure that can handle a geometrical resolved VG on an airfoil section.The airfoils of FFA-W3-301 and FFA-W3-360 were taken as the research objects.The results of calculation were compared with those of wind tunnel experiments.This method cannot accurately capture the details of a VG's flow field, but it can be used to qualitatively compare the effects of different VG parameters on the flow field.Mueller-Vahl et al. [31] studied the effect of chordwise position, spanwise spacing, and VG size on the aerodynamic characteristics of an airfoil by using force balance measurements for a 2D airfoil section.Wang et al. [15] investigated the influence of double-row VGs on the aerodynamic characteristics of the airfoil, and the development law of vortices along the flow direction was analyzed.The results showed that double-row VGs can effectively improve the lift coefficient of the airfoil and control airflow separation, which has a definite guiding significance for the application of VGs in wind turbine blades.The DU97-W-300 airfoil was adopted by Gao et al. [18] to study the effects of different combinations of VGs on the aerodynamic characteristics of an airfoil using a numerical method.Fouatih et al. [32] used an experimental method to study the effect of vortex generators on flow separation of the NACA 4415 airfoil, and they optimized the vortex generator's geometric parameters, including the installation position and installation angle.Velte and Hansen [33] carried out a three-dimensional particle image velocimetry for the DU91-w2-250 airfoil in a wind tunnel, and the impact of the VGs on the flow near stall was studied.The flow structures induced by the VGs and the separation-controlling effects on the airfoil were investigated in detail.
Hansen et al. [34] presented a VG aerodynamic configuration that can provide airfoils with a higher lift-drag ratio.The maximum lift of an airfoil reached 1.82, and the maximum lift-drag ratio increased from 110 to 114.5 by using this configuration of VGs.Kerho et al. [35] experimentally studied the effects of several VG structures on the aerodynamic performance of airfoils.The results demonstrated that for an airfoil with leading separation bubbles, the separation can be effectively controlled by installing VGs in the positions before the separation bubbles and after the peak pressure of the suction surface of the airfoil.Delft University also conducted relevant experimental studies on DU series airfoil sections [36].The vortex generator was installed at 0.2c from the leading edge.The results showed that the maximum lift coefficient increased from 1.55 to 1.97, and the drag coefficient increased before the stall angle and decreased after the stall angle.In addition to the above reports, the application of VGs in wind turbines is further discussed.Gyatt [37] analyzed and experimentally studied a small two-blade horizontal-axis wind turbine (diameter of 9.75 m) equipped with reverse-placement VGs.The arrangement of VGs can be divided into three types: inside half-blade spanwise, outside half-blade spanwise, and whole blade spanwise.Field measurements showed that the power of wind turbine blades with VGs increased by 20% when the wind speed was greater than 10 m/s, and when the wind speed was low, the reduction in the power was less than 4%.The results also showed that VGs should be installed in the outside half-blade spanwise.Sullivan [38] carried out wind farm measurements on the power and structural dynamic loads of a 2.5 MW horizontal-axis wind turbine blade with VGs.The results showed that the annual power could be increased by 20% when the VGs were installed at positions ranging from 20 to 100% of the blade spanwise, and the annual power could be increased by about 10% by arranging VGs in the 20-70% blade spanwise range.To solve the problem of the drag increasing as the lift increases in a rectangular arrangement of plate VGs, Kogaki et al. [28] adopted the airfoil VG, which not only ensures the function of VGs but also does not reduce the lift-drag ratio.However, this result can only be applied to small wind turbine blades with a low Reynolds number.Afjeh et al. [39] used lift line theory to simulate the flow around the blades, and the mirror eddy method was applied to study the effect of flow vortices generated by VGs on the flow around the blades.A fast calculation method for the aerodynamic characteristics of a horizontal-axis wind turbine with VGs was established and compared with the measured wind power of MOD-2, and the results showed that the method is feasible.
Although many scholars have studied the flow control mechanism of VGs, there are still many scientific problems related to the application of VGs for the flow control of wind turbine blades.For example, the optimum VG height to realize the best flow control effect has not been determined.Vortex generator height is often related to the local boundary layer thickness; the influence of the ratio of VG height to boundary layer thickness on the vorticity characteristics, as well as the influence of VG height on the lift and drag characteristics of airfoils, has not been reported yet.In this study, numerical simulation methods were used to study the influence of the VG height to boundary layer thickness ratio on the vortices' characteristics on a flat plate, and the influence of VG height on the aerodynamic characteristics of wind turbine airfoils was investigated in a wind tunnel.

Numerical Conditions
Fluent software was adopted for steady numerical calculations.The finite volume method was used for the discretization of the Navier-Stokes (N-S) equation with the second-order accuracy central difference scheme for spatial discretization, and the pressure-velocity coupling algorithm SIMPLE was adopted.The full-turbulence Spalart-Allmaras (S-A) model was selected as the turbulence model.The mesh of the computational domain was divided using ICEM software: 80 mesh nodes were arranged in the height and chord direction of the VG, and 200 mesh nodes were arranged downstream of the VG.The total number of mesh nodes was about 5,000,000.The height of the first layer near the VG was set to 1.2 × 10 −5 mm, with the wall y+ less than 1 over the VG surface.The mesh distribution on the surface and bottom of the VG is shown in Figure 1.
Energies 2019, 12, x FOR PEER REVIEW 3 of 19 studied a small two-blade horizontal-axis wind turbine (diameter of 9.75 m) equipped with reverseplacement VGs.The arrangement of VGs can be divided into three types: inside half-blade spanwise, outside half-blade spanwise, and whole blade spanwise.Field measurements showed that the power of wind turbine blades with VGs increased by 20% when the wind speed was greater than 10 m/s, and when the wind speed was low, the reduction in the power was less than 4%.The results also showed that VGs should be installed in the outside half-blade spanwise.Sullivan [38] carried out wind farm measurements on the power and structural dynamic loads of a 2.5 MW horizontal-axis wind turbine blade with VGs.The results showed that the annual power could be increased by 20% when the VGs were installed at positions ranging from 20 to 100% of the blade spanwise, and the annual power could be increased by about 10% by arranging VGs in the 20-70% blade spanwise range.To solve the problem of the drag increasing as the lift increases in a rectangular arrangement of plate VGs, Kogaki et al. [39] adopted the airfoil VG, which not only ensures the function of VGs but also does not reduce the lift-drag ratio.However, this result can only be applied to small wind turbine blades with a low Reynolds number.Afjeh et al. [40] used lift line theory to simulate the flow around the blades, and the mirror eddy method was applied to study the effect of flow vortices generated by VGs on the flow around the blades.A fast calculation method for the aerodynamic characteristics of a horizontal-axis wind turbine with VGs was established and compared with the measured wind power of MOD-2, and the results showed that the method is feasible.Although many scholars have studied the flow control mechanism of VGs, there are still many scientific problems related to the application of VGs for the flow control of wind turbine blades.For example, the optimum VG height to realize the best flow control effect has not been determined.Vortex generator height is often related to the local boundary layer thickness; the influence of the ratio of VG height to boundary layer thickness on the vorticity characteristics, as well as the influence of VG height on the lift and drag characteristics of airfoils, has not been reported yet.In this study, numerical simulation methods were used to study the influence of the VG height to boundary layer thickness ratio on the vortices' characteristics on a flat plate, and the influence of VG height on the aerodynamic characteristics of wind turbine airfoils was investigated in a wind tunnel.

Numerical Conditions
Fluent software was adopted for steady numerical calculations.The finite volume method was used for the discretization of the Navier-Stokes (N-S) equation with the second-order accuracy central difference scheme for spatial discretization, and the pressure-velocity coupling algorithm SIMPLE was adopted.The full-turbulence Spalart-Allmaras (S-A) model was selected as the turbulence model.The mesh of the computational domain was divided using ICEM software: 80 mesh nodes were arranged in the height and chord direction of the VG, and 200 mesh nodes were arranged downstream of the VG.The total number of mesh nodes was about 5,000,000.The height of the first layer near the VG was set to 1.2 × 10 −5 mm, with the wall y+ less than 1 over the VG surface.The mesh distribution on the surface and bottom of the VG is shown in Figure 1.In order to observe the motion of vortices developed downstream of the VG, the computational domain was set to 144H (where H is the height of the VG), and the VG was set to 9H away from the inlet.The width and height of the computational domain were 20H.The boundary conditions included the velocity inlet, pressure outlet, symmetrical boundary on both sides of the computational domain, free outflow boundary condition on the upper boundary, and non-sliding solid wall on the bottom and VG.The computational domain and boundary conditions are shown in Figure 2.
Energies 2019, 12, x FOR PEER REVIEW 4 of 19 In order to observe the motion of vortices developed downstream of the VG, the computational domain was set to 144H (where H is the height of the VG), and the VG was set to 9H away from the inlet.The width and height of the computational domain were 20H.The boundary conditions included the velocity inlet, pressure outlet, symmetrical boundary on both sides of the computational domain, free outflow boundary condition on the upper boundary, and non-sliding solid wall on the bottom and VG.The computational domain and boundary conditions are shown in Figure 2. The boundary layer often presents a turbulent form for VGs; therefore, this paper mainly studies the vortices' characteristics for VGs with different turbulent boundary layer heights.Changing the height of the VG (H) means altering the height of h/δ (as shown in Figure 3).Therefore, the purpose of studying VG height was achieved by changing the boundary layer height of the velocity type of the turbulence at the inlet.The inlet turbulence velocity type is given as follows: where u is the local velocity, δ is the thickness of the boundary layer, U is the main velocity, and y is the normal height of the wall.Table 1 shows the geometric parameters of the VG.The VG was set as a wall without thickness.Inflow velocity was U = 82 m/s, and the Reynolds number based on VG height was about 3 × 10 4 .Table 2 shows the ratio of VG Height to Different Boundary Layer (BL) Thickness.The boundary layer often presents a turbulent form for VGs; therefore, this paper mainly studies the vortices' characteristics for VGs with different turbulent boundary layer heights.Changing the height of the VG (H) means altering the height of h/δ (as shown in Figure 3).Therefore, the purpose of studying VG height was achieved by changing the boundary layer height of the velocity type of the turbulence at the inlet.The inlet turbulence velocity type is given as follows: where u is the local velocity, δ is the thickness of the boundary layer, U is the main velocity, and y is the normal height of the wall.In order to observe the motion of vortices developed downstream of the VG, the computational domain was set to 144H (where H is the height of the VG), and the VG was set to 9H away from the inlet.The width and height of the computational domain were 20H.The boundary conditions included the velocity inlet, pressure outlet, symmetrical boundary on both sides of the computational domain, free outflow boundary condition on the upper boundary, and non-sliding solid wall on the bottom and VG.The computational domain and boundary conditions are shown in Figure 2. The boundary layer often presents a turbulent form for VGs; therefore, this paper mainly studies the vortices' characteristics for VGs with different turbulent boundary layer heights.Changing the height of the VG (H) means altering the height of h/δ (as shown in Figure 3).Therefore, the purpose of studying VG height was achieved by changing the boundary layer height of the velocity type of the turbulence at the inlet.The inlet turbulence velocity type is given as follows: where u is the local velocity, δ is the thickness of the boundary layer, U is the main velocity, and y is the normal height of the wall.Table 1 shows the geometric parameters of the VG.The VG was set as a wall without thickness.Inflow velocity was U = 82 m/s, and the Reynolds number based on VG height was about 3 × 10 4 .Table 2 shows the ratio of VG Height to Different Boundary Layer (BL) Thickness.Table 1 shows the geometric parameters of the VG.The VG was set as a wall without thickness.Inflow velocity was U = 82 m/s, and the Reynolds number based on VG height was about 3 × 10 4 .Table 2 shows the ratio of VG Height to Different Boundary Layer (BL) Thickness.Equations ( 2) and ( 3) were used to define the lift and drag coefficients of the VG, respectively: where F l and F d are the total lift and drag forces of VG; ρ is the air density, and ρ = 1.225 kg/m 3 ; V ∞ is the free stream velocity of the VG, V ∞ = 82 m/s; L is the bottom chord length of the VG, L = 17 mm.In order to prove the adequacy of the grid resolution, grid-independent computation was carried out.Three sets of grids with different densities were calculated: 2,560,000 (coarse), 3,645,000 (medium), and 500,000 (fine).In addition, the Richardson extrapolation method [40] was used to verify the validity of the grid resolution.Table 3 shows the result of Richardson extrapolation.The calculation case is H = 1.0δ.In the table, C l , C d , and C l /C d are the results for the different grid densities, respectively.RE is the value obtained by Richardson extrapolation, R is the ratio of errors, and p is the order of accuracy.
The convergence conditions are 0 < R < 1 (monotonic convergence); R < 0 (oscillation convergence); and 1 < R (divergence).Since all the calculated results for R are less than 1, the mesh computation is monotonically convergent.Therefore, the fine grid was used for the presented computations.

Experiment Setup
The experiment was conducted in the wind tunnel of North China Electric Power University (NCEPU), and the wind tunnel is a large closed return wind tunnel.The airfoil was tested in a 2-D measurement section with dimensions of 1.5 × 3 × 4.5 m 3 .The maximum wind speed was 62 m/s.The measurement section was equipped with an electronic scanning valve of 256 channels and 64 channels, a PXI data acquisition system, etc. Upstream of the airfoil fitted with a Pitot tube, the total pressure of the inlet pipe was measured according to the Pitot tube principle.The side wall of the Pitot tube pipe had small holes, which were used to measure the static pressure of the upstream flow.Installed in the middle section of the airfoil surface were 96 pressure taps, which were used to measure the pressure coefficient distribution on the airfoil surface.The measuring accuracy is five ten thousandths.The Cp curve is the average result of time; the time interval is 0.1 second, and 10 seconds was used for sample time-averaging.Downstream of the airfoil, there were 111 total pressure measuring holes and 5 static pressure measuring holes (called the wake rake) that are used to calculate the airfoil drag coefficient.The wind tunnel measurement section and test airfoil are shown in Figure 4.The DU93-W-210 airfoil was adopted for the test; the relative thickness of the airfoil is 21%; the chord length (C) is 0.8 m; the span of the airfoil is 1.5 m.The DU93-W-210 airfoil is a special airfoil for wind turbines and was developed and designed by Delft University.It is usually arranged at 70% of the blade's outward span.Depending on the blade length, the actual flow conditions of the airfoil are different, but the optimum operating angle of attack is generally between 4° and 6°.The VGs were located at 0.2C from the leading edge.The test airfoil with VGs is illustrated in Figure 5.The arrangement and height of the VGs are presented in Figure 6, and the geometric and arrangement parameters of the VGs are shown in Table 4.The Reynolds number based on the chord of the airfoil is 1 × 10 6 .The tunnel speed was 18.2 m/s, and the angle of attack of the airfoil experiment ranged from −10° to 25°.The DU93-W-210 airfoil was adopted for the test; the relative thickness of the airfoil is 21%; the chord length (C) is 0.8 m; the span of the airfoil is 1.5 m.The DU93-W-210 airfoil is a special airfoil for wind turbines and was developed and designed by Delft University.It is usually arranged at 70% of the blade's outward span.Depending on the blade length, the actual flow conditions of the airfoil are different, but the optimum operating angle of attack is generally between 4 • and 6 • .The VGs were located at 0.2C from the leading edge.The test airfoil with VGs is illustrated in Figure 5.The arrangement and height of the VGs are presented in Figure 6, and the geometric and arrangement parameters of the VGs are shown in Table 4.The Reynolds number based on the chord of the airfoil is 1 × 10 6 .The tunnel speed was 18.2 m/s, and the angle of attack of the airfoil experiment ranged from −10 • to 25  The DU93-W-210 airfoil was adopted for the test; the relative thickness of the airfoil is 21%; the chord length (C) is 0.8 m; the span of the airfoil is 1.5 m.The DU93-W-210 airfoil is a special airfoil for wind turbines and was developed and designed by Delft University.It is usually arranged at 70% of the blade's outward span.Depending on the blade length, the actual flow conditions of the airfoil are different, but the optimum operating angle of attack is generally between 4° and 6°.The VGs were located at 0.2C from the leading edge.The test airfoil with VGs is illustrated in Figure 5.The arrangement and height of the VGs are presented in Figure 6, and the geometric and arrangement parameters of the VGs are shown in Table 4.The Reynolds number based on the chord of the airfoil is 1 × 10 6 .The tunnel speed was 18.2 m/s, and the angle of attack of the airfoil experiment ranged from −10° to 25°.       Figure 7b shows that the drag coefficient of the airfoil is low before the stall angle of attack (8 • ), and the drag coefficient measured in this work is similar to that measured by NPU and Delft.The drag coefficient measured after stall in this work is greater than that measured by NPU but similar to that measured by Delft.Because the aerodynamic performance of the airfoil is unstable after stall, the lift and drag coefficients have to be revised, and the mechanism of each correction method is different; however, our focus is on the aerodynamic performance before stall.Figure 7 shows compares this work with that by North China Electric Power University (NCEPU), Northwest Polytechnic University (NPU), and Delft University (Delft).It can be seen that in the linear section of the lift coefficient, the experimental data achieved in this work are similar to that obtained by NPU and Delft.After airfoil stall, the lift coefficient measured in this work is the highest compared with the results of the other three works.The lift coefficient first decreases and then increases slightly.Figure 7b shows that the drag coefficient of the airfoil is low before the stall angle of attack (8°), and the drag coefficient measured in this work is similar to that measured by NPU and Delft.The drag coefficient measured after stall in this work is greater than that measured by NPU but similar to that measured by Delft.Because the aerodynamic performance of the airfoil is unstable after stall, the lift and drag coefficients have to be revised, and the mechanism of each correction method is different; however, our focus is on the aerodynamic performance before stall.

Analysis of CFD Calculation Results
Figure 8 shows the static pressure iso-contours on the leeward side of the VG.It can be seen that the low-pressure area on the leeward side of the VG appears at the middle position of the chord, and it is not generated at the leading edge of the tip.The main reason for this phenomenon is that the kinetic energy of the bottom fluid in the boundary layer is very weak because of the viscous wall.Instead of forming leading-edge separation vortices, concentrated vortices are induced when the fluid passes over the swept leading edge of the VG from the pressure surface.However, there are still low-pressure zones on the leeward side of the VG that are caused by the concentrated vortices on that side.With the increase in the height of the VG (H) relative to the thickness of the boundary layer (δ), the pressure value on the leeward side of the VG is lower.The result shows that the larger the height of the VG, the greater the fluid kinetic energy in the boundary layer; thus, the concentrated

Analysis of CFD Calculation Results
Figure 8 shows the static pressure iso-contours on the leeward side of the VG.It can be seen that the low-pressure area on the leeward side of the VG appears at the middle position of the chord, and it is not generated at the leading edge of the tip.The main reason for this phenomenon is that the kinetic energy of the bottom fluid in the boundary layer is very weak because of the viscous wall.Instead of forming leading-edge separation vortices, concentrated vortices are induced when the fluid passes over the swept leading edge of the VG from the pressure surface.However, there are still low-pressure zones on the leeward side of the VG that are caused by the concentrated vortices on that side.With the increase in the height of the VG (H) relative to the thickness of the boundary layer (δ), the pressure value on the leeward side of the VG is lower.The result shows that the larger the height of the VG, the greater the fluid kinetic energy in the boundary layer; thus, the concentrated vortices generated by the VG will show a stronger trend, which leads to a lower pressure in the low-pressure zone on the leeward side of the VG.  Figure 9 shows the distribution of the static pressure coefficient on the leeward surface of the VG; in the figure, S is the local height, l is the bottom length of the VG.Firstly, the maximum amplitude of the low-pressure coefficients occurs at about x = 0.4l, and the strongest vortex intensity can be found at 0.4l.With the increase in the flow distance, the strength of the vortex core decreases gradually, and the low-pressure amplitude of the static pressure coefficient on the leeward side of the VG also presents a decreasing trend.This result shows that when the VG is placed on the wall, Figure 9 shows the distribution of the static pressure coefficient on the leeward surface of the VG; in the figure, S is the local height, l is the bottom length of the VG.Firstly, the maximum amplitude of the low-pressure coefficients occurs at about x = 0.4l, and the strongest vortex intensity can be found at 0.4l.With the increase in the flow distance, the strength of the vortex core decreases gradually, and the low-pressure amplitude of the static pressure coefficient on the leeward side of the VG also presents a decreasing trend.This result shows that when the VG is placed on the wall, concentrated vortices are generated at the upper end of the VG, and the intensity of the vortices reaches the maximum value near the position 0.4l.This is because the vortex gradually dissipates and is deflected away from the VG.Secondly, at any chord position, the magnitude of the static pressure coefficient increases with the increase in VG height, and the magnitudes of the static pressure coefficients are close to each other when H = 1.0δ, 1.5δ, and 2.0δ.When the value of H is 0.1δ and 0.2δ, the gentle pressure amplitude of the static pressure coefficient is lower, and when H = 0.5δ, the static pressure coefficient amplitude is in the middle position.Figure 9 shows the distribution of the static pressure coefficient on the leeward surface of the VG; in the figure, S is the local height, l is the bottom length of the VG.Firstly, the maximum amplitude of the low-pressure coefficients occurs at about x = 0.4l, and the strongest vortex intensity can be found at 0.4l.With the increase in the flow distance, the strength of the vortex core decreases gradually, and the low-pressure amplitude of the static pressure coefficient on the leeward side of the VG also presents a decreasing trend.This result shows that when the VG is placed on the wall, concentrated vortices are generated at the upper end of the VG, and the intensity of the vortices reaches the maximum value near the position 0.4l.This is because the vortex gradually dissipates and is deflected away from the VG.Secondly, at any chord position, the magnitude of the static pressure coefficient increases with the increase in VG height, and the magnitudes of the static pressure coefficients are close to each other when H = 1.0δ, 1.5δ, and 2.0δ.When the value of H is 0.1δ and 0.2δ, the gentle pressure amplitude of the static pressure coefficient is lower, and when H = 0.5δ, the static pressure coefficient amplitude is in the middle position.Figure 10 shows the distribution of static pressure coefficients at the center of the vortex core.It can be seen that the static pressure coefficient at the center of the vortex core first increases and then decreases in the direction of the flow and reaches the maximum at the position of 0.5l.The higher the H, the lower the static pressure at the center of the vortex core.
Energies 2019, 12, x FOR PEER REVIEW 9 of 19 Figure 10 shows the distribution of static pressure coefficients at the center of the vortex core.It can be seen that the static pressure coefficient at the center of the vortex core first increases and then decreases in the direction of the flow and reaches the maximum at the position of 0.5l.The higher the H, the lower the static pressure at the center of the vortex core.Figure 11 shows the vortex core area and vortex circulation varying with the flow direction in the range of the VG surface.The area of the vortex core is almost the same at different VG heights, but it is slightly different at the trailing edge of the VG surface.The area of the vortex core is approximately exponentially distributed as the flow direction distance increases.With the increase in VG height, the vortex circulation increases gradually, and the difference becomes progressively Figure 11 shows the vortex core area and vortex circulation varying with the flow direction in the range of the VG surface.The area of the vortex core is almost the same at different VG heights, but it is slightly different at the trailing edge of the VG surface.The area of the vortex core is approximately exponentially distributed as the flow direction distance increases.With the increase in VG height, the vortex circulation increases gradually, and the difference becomes progressively larger as the flow becomes closer to the trailing edge of the VG. Figure 11 shows the vortex core area and vortex circulation varying with the flow direction in the range of the VG surface.The area of the vortex core is almost the same at different VG heights, but it is slightly different at the trailing edge of the VG surface.The area of the vortex core is approximately exponentially distributed as the flow direction distance increases.With the increase in VG height, the vortex circulation increases gradually, and the difference becomes progressively larger as the flow becomes closer to the trailing edge of the VG. Figure 12 shows the variation in the vortex core area and vortex circulation with the flow direction distance downstream of the VG.The area of the vortex core downstream of the VG varies linearly with the distance of the flow direction, which is slightly different from that on the VG surface.The higher the VG height, the larger the area of the vortex core, and the vortex core area difference downstream of the VG is more pronounced.However, the variation in the core circulation in the downstream direction along the VG approximately exponentially dissipates, which is very different from that on the VG surface.It can be seen that the maximum position of the VG vorticity appears at the trailing edge of the VG. Figure 12 shows the variation in the vortex core area and vortex circulation with the flow direction distance downstream of the VG.The area of the vortex core downstream of the VG varies linearly with the distance of the flow direction, which is slightly different from that on the VG surface.The higher the VG height, the larger the area of the vortex core, and the vortex core area difference downstream of the VG is more pronounced.However, the variation in the core circulation in the downstream direction along the VG approximately exponentially dissipates, which is very different from that on the VG surface.It can be seen that the maximum position of the VG vorticity appears at the trailing edge of the VG.  Figure 13 shows the variation in the vortex core area and vortex circulation as VG height increases at its trailing edge.It is shown that the vortex core area is the largest when H/δ = 1.0.The area of the vortex core decreases slightly when H/δ > 1.0 and increases rapidly when H/δ < 1.0.Vortex circulation is proportional to VG height: with the increase in VG height, the vortex circulation increases gradually, but it does not possess a linear relationship with VG height as it increases.When H/δ < 1.0, the slope of vortex circulation increases when VG height is larger, but when H/δ > 1.0, the slope of vortex circulation increases as VG height decreases.Figure 13 shows the variation in the vortex core area and vortex circulation as VG height increases at its trailing edge.It is shown that the vortex core area is the largest when H/δ = 1.0.The area of the vortex core decreases slightly when H/δ > 1.0 and increases rapidly when H/δ < 1.0.Vortex circulation is proportional to VG height: with the increase in VG height, the vortex circulation increases gradually, but it does not possess a linear relationship with VG height as it increases.When H/δ < 1.0, the slope of vortex circulation increases when VG height is larger, but when H/δ > 1.0, the slope of vortex circulation increases as VG height decreases.
increases at its trailing edge.It is shown that the vortex core area is the largest when H/δ = 1.0.The area of the vortex core decreases slightly when H/δ > 1.0 and increases rapidly when H/δ < 1.0.Vortex circulation is proportional to VG height: with the increase in VG height, the vortex circulation increases gradually, but it does not possess a linear relationship with VG height as it increases.When H/δ < 1.0, the slope of vortex circulation increases when VG height is larger, but when H/δ > 1.0, the slope of vortex circulation increases as VG height decreases.Figure 14a shows the velocity profiles in the boundary layer of the inflow, and Figure 14b shows the average velocity of the fluid in the range of VG height.From Figure 14a, when the VG is in boundary layers with different thicknesses and velocity profiles of the inflow, the kinetic energy of the fluid in the boundary layer becomes progressively lower within the VG height range as the height of the boundary layer increases.It can also be seen from Figure 14b that with the increase in VG height, i.e., the decrease in the thickness of the boundary layer, the average velocity of the fluid in the VG height range becomes increasingly higher.The increase in average velocity is exponentially distributed with the height of the VG.When H/δ < 1, the slope of average velocity in the boundary layer increases as the height of the VG presents a larger trend; after H/δ > 1, the slope of average velocity in the boundary layer decreases with the increase in VG height.Figure 14a shows the velocity profiles in the boundary layer of the inflow, and Figure 14b shows the average velocity the fluid in the range of VG height.From Figure 14a, when the VG is in boundary layers with different thicknesses and velocity profiles of the inflow, the kinetic energy of the fluid in the boundary layer becomes progressively lower within the VG height range as the height of the boundary layer increases.It can also be seen from Figure 14b that with the increase in VG height, i.e., the decrease in the thickness of the boundary layer, the average velocity of the fluid in the VG height range becomes increasingly higher.The increase in average velocity is exponentially distributed with the height of the VG.When H/δ < 1, the slope of average velocity in the boundary layer increases as the height of the VG presents a larger trend; after H/δ > 1, the slope of average velocity in the boundary layer decreases with the increase in VG height.Figure 15 shows the variation in the lift, drag coefficient, and lift-drag ratio of the VG.With the increase in H/δ, the lift coefficient of the VG increases gradually.When H/δ = 2.0, the maximum lift coefficient of the VG is about 0.65.The drag coefficient of the VG increases with the increase in H/δ, and the maximum drag coefficient is 0.25 when H/δ = 2.0.The lift of the VG comes from the vortexinduced lift, and the drag is primarily pressure drag, that is, the profile drag of the VG.The variation in the lift coefficient and drag coefficient with the increase in H/δ basically shows the same trend.The dynamic energy of the fluid in VG height, vortex circulation, lift coefficient, and drag coefficient are essentially consistent with the variation in H/δ.The result shows that the vortex strength and drag of VG vortices are closely related to the H/δ value, and the greater the height of the VG, the less obvious the change effect.The lift-drag ratio of the VG reaches the maximum at H/δ = 0.5.It would be advisable to have a vorticity with a certain intensity to control the flow, but the profile drag is not too large.However, whether the ratio of H/δ is optimal should be discussed in relation to the kinetic energy of the fluid in the boundary layer.Figure 15 shows the variation in the lift, drag coefficient, and lift-drag ratio of the VG.With the increase in H/δ, the lift coefficient of the VG increases gradually.When H/δ = 2.0, the maximum lift coefficient of the VG is about 0.65.The drag coefficient of the VG increases with the increase in H/δ, and the maximum drag coefficient is 0.25 when H/δ = 2.0.The lift of the VG comes from the vortex-induced lift, and the drag is primarily pressure drag, that is, the profile drag of the VG.The variation in the lift coefficient and drag coefficient with the increase in H/δ basically shows the same trend.The dynamic energy of the fluid in VG height, vortex circulation, lift coefficient, and drag coefficient are essentially consistent with the variation in H/δ.The result shows that the vortex strength and drag of VG vortices are closely related to the H/δ value, and the greater the height of the VG, the less obvious the change effect.The lift-drag ratio of the VG reaches the maximum at H/δ = 0.5.It would be advisable to have a vorticity with a certain intensity to control the flow, but the profile drag is not too large.However, whether the ratio of H/δ is optimal should be discussed in relation to the kinetic energy of the fluid in the boundary layer.
induced lift, and the drag is primarily pressure drag, that is, the profile drag of the VG.The variation in the lift coefficient and drag coefficient with the increase in H/δ basically shows the same trend.The dynamic energy of the fluid in VG height, vortex circulation, lift coefficient, and drag coefficient are essentially consistent with the variation in H/δ.The result shows that the vortex strength and drag of VG vortices are closely related to the H/δ value, and the greater the height of the VG, the less obvious the change effect.The lift-drag ratio of the VG reaches the maximum at H/δ = 0.5.It would be advisable to have a vorticity with a certain intensity to control the flow, but the profile drag is not too large.However, whether the ratio of H/δ is optimal should be discussed in relation to the kinetic energy of the fluid in the boundary layer.Figure 16 shows the axial velocity contours in different flow directions downstream of the VG.The kinetic energy of the fluid in the vortex core height is ineffective when the value of H/δ is small.The high-energy fluid outside the vortex core is contained in the bottom layer of the boundary layer from the rotation of concentrated vortices when VG is used for flow control.From the axial velocity contours, it can be seen that with the rotation of the vortex core, the fluid kinetic energy outside the vortex core decreases, while the fluid's streamwise kinetic energy inside the vortex core increases.Figure 16 shows the axial velocity contours in different flow directions downstream of the VG.The kinetic energy of the fluid in the vortex core height is ineffective when the value of H/δ is small.The high-energy fluid outside the vortex core is contained in the bottom layer of the boundary layer from the rotation of concentrated vortices when VG is used for flow control.From the axial velocity contours, it can be seen that with the rotation of the vortex core, the fluid kinetic energy outside the vortex core decreases, while the fluid's streamwise kinetic energy inside the vortex core increases.Figure 17 shows the axial velocity distribution for the normal height at the vortex core position downstream of the VG.The horizontal line is placed at the center of the vortex core.With the increase in the flow distance, the area of the vortex core varies at different VG heights.The height of the vortex core marked in the figure is the height of the vortex core when H = 1.0δ.Figure 17 shows the axial velocity distribution for the normal height at the vortex core position downstream of the VG.The horizontal line is placed at the center of the vortex core.With the increase in the flow distance, the area of the vortex core varies at different VG heights.The height of the vortex core marked in the figure is the height of the vortex core when H = 1.0δ.H=1.5δ H=2.0δ Figure 17.Velocity profiles at the vortex core downstream of the VG.
Table 5 shows the relative increase in near-wall (y/H < 1) momentum for different VGs.Firstly, it can be observed that the center height of the vortex core increases with the increase in the flow direction distance.Secondly, the velocity profiles in the boundary layer are basically distributed in an S-line along the center of the vortex core, and the increase in fluid kinetic energy inside the vortex core corresponds to the decrease in fluid kinetic energy outside.From Table 5, comparing the kinetic energy of the bottom fluid of the boundary layer at different H/δ values, it can be seen that the kinetic energy for the bottom fluid of the boundary layer is the lowest when H = 0.1δ and 0.2δ, and there is little difference between them.When H = 0.5δ, the value is intermediate; when H = 1.0δ, 1.5δ, and 2.0δ, the values are almost the same, and the curves basically coincide.This result shows that the kinetic energy of the fluid in the boundary layer remains largely unchanged with the increase in VG height.However, the profile drag increases with the increase in VG height, as shown in Figure 15.For flow control, it is hoped that VG can generate concentrated vortices with sufficient intensity to increase the fluid kinetic energy in the boundary layer, and this will not increase the drag of the boundary layer by much.In terms of the increasing fluid kinetic energy in the boundary layer, the value of H = 1.0δ is the best choice.

Analysis of Experimental Results
The results of CFD calculation show that when the VG height is about 1.0δ, the fluid kinetic energy in the near-wall basically reaches the maximum, and as the VG height continues to increase, the kinetic energy remains virtually unchanged.The experimental design was guided by CFD calculation results.The thickness of the BL of the airfoil is different at different angles of attack.To quantify the relationship between VG height and BL thickness, the airfoil was computed using Xfoil; the thickness of the BL at the VG installation position is about 5.9 mm when the maximum lift coefficient angle of attack (8 • ) of the airfoil is reached.Considering the machining accuracy of VGs, three VGs were designed with different heights: 4 mm, 6 mm, and 8 mm.The ratio of VG height to BL thickness is 0.66δ, 1.0δ, and 1.33δ, respectively.
Figure 18 shows the lift coefficient, drag coefficient, and lift-drag ratio of the airfoil with and without the VG airfoils.Table 5 shows the relative variation in C l , C d , and lift-drag ratio.Before the stall angle of attack (8 • ), the lift coefficients for the airfoil with VGs and without VGs (clean) are similar to each other in the linear section; after the 8 • angle of attack, the lift coefficient of the airfoil without VGs remains basically unchanged, while the lift coefficient of the airfoil with VGs continues to increase.The stall angle of the clean airfoil is 8 • , while that of the airfoil with VGs is 18 • , so there is a 10 • increase in the stall angle of the airfoil.When the angle of attack reaches 23 • , the separation position of the clean airfoil is located upstream of 0.2C (VG installation position), and the VGs cannot control the flow separation, so the lift coefficient of the airfoil with VGs is no different from that of the clean airfoil.The comparison of the differences among the three VG heights shows that the maximum lift coefficient of the airfoil with VGs increases by 42.8% compared with that of the clean airfoil when H = 0.66δ; when H = 1.0δ and 1.33δ, the maximum lift coefficient increases by 48.7% and 48.6%, respectively.From the analysis of the maximum lift coefficient, it can be found that the VG flow control effect is the best when H = 1.0δ.
From Figure 18b and Table 6, it can be seen that the drag coefficient for H = 1.33δ is larger before the 8 • angle of attack.However, beyond the 8 • angle of attack, the drag coefficient of the clean airfoil increases rapidly.For this VG height, the drag coefficient of the airfoil with VGs is much smaller than that of the clean airfoil.With respect to the difference among the three VG heights, before the angle of attack of 18 • , the drag coefficient of the airfoil with H = 1.33δ is the largest, while the drag coefficients for the other two VG Heights do not much differ.The drag coefficient of the airfoil with H = 0.66δ is 84.9% lower than that of the clean airfoil; when H = 1.0δ, the drag coefficient of the airfoil decreases by 83.2%; when H = 1.33δ, the drag coefficient of the airfoil decreases by 76.8%.So, among these three heights, the drag coefficient of the airfoil with H = 0.66δ is reduced the most because the VG height is lower, and the profile drag of the VGs is minimum.
For the lift-drag ratio of the airfoil, the maximum lift-drag ratio of the clean airfoil is the highest, while the maximum lift-drag ratio of the airfoil with VGs is decreased.Therefore, the VGs cannot increase the maximum lift-drag ratio of the airfoil, but raising the angle of attack produces the best lift-drag ratio for the airfoil.The lift-drag ratio of the airfoil with VGs is larger than that of the clean airfoil after the angle of attack of 10 • .Compared with the maximum lift-drag ratio of the airfoil with and without VGs, the maximum lift-drag ratio of the airfoil with H = 0.66δ is 19.7% lower than that of the clean airfoil; the maximum lift-drag ratio of the airfoil with H = 1.0δ and 1.33δ is decreased by 19.3% and 51.2%, respectively.However, when the angle of attack is 18 • (stall angle of the airfoil with VGs), the lift-drag ratio of the airfoil with H = 0.66δ increases by 880.3% compared with that of the clean airfoil; the lift-drag ratio increases by 821.8% when H = 1.0δ; and the lift-drag ratio increases by 564.8% when H = 1.33δ.So, in summary, if considering increases in the lift coefficient and the lift-drag ratio, H = 1.0δ has a good effect on controlling airfoil flow separation.If considering the optimal lift-drag ratio, H = 0.66δ can meet the requirements.
19.3% and 51.2%, respectively.However, when the angle of attack is 18° (stall angle of the airfoil with VGs), the lift-drag ratio of the airfoil with H = 0.66δ increases by 880.3% compared with that of the clean airfoil; the lift-drag ratio increases by 821.8% when H = 1.0δ; and the lift-drag ratio increases by 564.8% when H = 1.33δ.So, in summary, if considering increases in the lift coefficient and the liftdrag ratio, H = 1.0δ has a good effect on controlling airfoil flow separation.If considering the optimal lift-drag ratio, H = 0.66δ can meet the requirements.Figure 19 shows the static pressure coefficient (Cp) distribution on the airfoil surface at different angles of attack.The pressure taps are located on the pressure side of the two VGs, so the Cp curve will be prominently downward near the installation locations of the VGs.The Cp curves of the airfoil with and without VGs are similar to each other at 0° and 5° angles of attack, and there are no differences among the three VG heights.A pressure plateau appears on the Cp curve of the clean airfoil about 88H away from the installation position of the VGs when the angle of attack is 10°.The reverse pressure gradient of the clean airfoil is about −0.77 kpa/m; with the action of the VGs, the  Figure 19 shows the static pressure coefficient (C p ) distribution on the airfoil surface at different angles of attack.The pressure taps are located on the pressure side of the two VGs, so the C p curve will be prominently downward near the installation locations of the VGs.The C p curves of the airfoil with and without VGs are similar to each other at 0 • and 5 • angles of attack, and there are no differences among the three VG heights.A pressure plateau appears on the C p curve of the clean airfoil about 88H away from the installation position of the VGs when the angle of attack is 10 • .The reverse pressure gradient of the clean airfoil is about −0.77 kpa/m; with the action of the VGs, the pressure plateau on the C p curve for the airfoil with VGs disappears, and the C p curves do not differ among the three VGs heights.With the increase in the angle of attack, the separation position of the clean airfoil is closer to the installation position of the VGs, and this configuration is beneficial to the control flow separation by the VGs when the angle of attack is 15 • , 18 • , and 20 • .However, with the increase in the angle of attack, the reverse pressure gradient for the clean airfoil becomes progressively larger.Therefore, with the increase in the angle of attack, the VGs cannot completely inhibit the flow separation, although they can postpone it.
clean airfoil is closer to the installation position of the VGs, and this configuration is beneficial to the control flow separation by the VGs when the angle of attack is 15°, 18°, and 20°.However, with the increase in the angle of attack, the reverse pressure gradient for the clean airfoil becomes progressively larger.Therefore, with the increase in the angle of attack, the VGs cannot completely inhibit the flow separation, although they can postpone it.The separation position of the airfoil is very close to the installation positions of the VGs, and the reverse pressure gradient is further increased where a large pressure plateau appears on the Cp curve for the airfoil with VGs when the angle of attack is 22°.In this scenario (22° angle of attack), the amplitude of the Cp curve is slightly lower at H = 1.33δ, and the position of the pressure plateau appears earlier than for the other two VG heights.When the angle of attack is 25°, because the flow separation of the clean airfoil has been occurring upstream of the VGs, the Cp curve of the airfoil with The separation position of the airfoil is very close to the installation positions of the VGs, and the reverse pressure gradient is further increased where a large pressure plateau appears on the C p curve for the airfoil with VGs when the angle of attack is 22 • .In this scenario (22 • angle of attack), the amplitude of the C p curve is slightly lower at H = 1.33δ, and the position of the pressure plateau appears earlier than for the other two VG heights.When the angle of attack is 25 • , because the flow separation of the clean airfoil has been occurring upstream of the VGs, the C p curve of the airfoil with VGs coincides with that of the clean airfoil.
Figure 20 shows the distribution of the total pressure coefficient of the wake rake.It can be seen that the total pressure loss of the airfoil with VGs is higher than that of the clean airfoil at 0 • , 5 • , and 10 • angles of attack.The total pressure loss for H = 1.33δ is higher than that for the other two VG heights at the 5 • and 10 • angles of attack.The total pressure loss of the airfoil with VGs is obviously smaller than that of the clean airfoil at 15 • and 18 • angles of attack; in this case, there is little difference among the three VG heights.The total pressure loss of the airfoil with VGs is still less than that of the clean airfoil at 20 • and 22 • angles of attack, for which the total pressure loss is the smallest when H = 0.66δ, and there is little difference between the other two VG heights.The total pressure coefficient curves of the airfoils with VGs coincide with that of the clean airfoil when the angle of attack is 25  The separation position of the airfoil is very close to the installation positions of the VGs, and the reverse pressure gradient is further increased where a large pressure plateau appears on the Cp curve for the airfoil with VGs when the angle of attack is 22°.In this scenario (22° angle of attack), the amplitude of the Cp curve is slightly lower at H = 1.33δ, and the position of the pressure plateau appears earlier than for the other two VG heights.When the angle of attack is 25°, because the flow separation of the clean airfoil has been occurring upstream of the VGs, the Cp curve of the airfoil with VGs coincides with that of the clean airfoil.
Figure 20 shows the distribution of the total pressure coefficient of the wake rake.It can be seen that the total pressure loss of the airfoil with VGs is higher than that of the clean airfoil at 0°, 5°, and 10° angles of attack.The total pressure loss for H = 1.33δ is higher than that for the other two VG heights at the 5° and 10° angles of attack.The total pressure loss of the airfoil with VGs is obviously smaller than that of the clean airfoil at 15° and 18° angles of attack; in this case, there is little difference among the three VG heights.The total pressure loss of the airfoil with VGs is still less than that of the clean airfoil at 20° and 22° angles of attack, for which the total pressure loss is the smallest when H = 0.66δ, and there is little difference between the other two VG heights.The total pressure coefficient curves of the airfoils with VGs coincide with that of the clean airfoil when the angle of attack is 25°.

Conclusions
When VGs are used for flow control, they are often put on the surface of an object.The characteristics of the vortices generated by VGs are quite different from those of aeronautical delta wings.The VGs are located on the surface of the object, and the concentrated vortices are not generated at the leading edge of the tip; instead, the fluid flows from the windward side of the VGs to the leeward surface, forming a concentrated vortex on the leeward side, and the intensity of the vortex reaches the maximum at the trailing edge of the VGs.There is a logarithmic relationship between the vortex circulation and VG height, and the average kinetic energy of the fluid in the boundary layer also has a logarithmic relationship with the VG height, so the vortex intensity of the VGs is proportional to the average kinetic energy of the fluid in the boundary layer.The greater the height of the VGs, the stronger the intensity of the concentrated vortices, and the stronger the flow control ability, but the greater the profile drag produced by the VGs.The lift and drag of the VG

Conclusions
When VGs are used for flow control, they are often put on the surface of an object.The characteristics of the vortices generated by VGs are quite different from those of aeronautical delta wings.The VGs are located on the surface of the object, and the concentrated vortices are not generated at the leading edge of the tip; instead, the fluid flows from the windward side of the VGs to the leeward surface, forming a concentrated vortex on the leeward side, and the intensity of the vortex reaches the maximum at the trailing edge of the VGs.There is a logarithmic relationship between the

Figure 2 .
Figure 2. Computational domain and boundary conditions.

Figure 2 .
Figure 2. Computational domain and boundary conditions.

Figure 2 .
Figure 2. Computational domain and boundary conditions.

Figure 4 .
Figure 4. Wind tunnel measurement section size and test airfoil.

Figure 4 .
Figure 4. Wind tunnel measurement section size and test airfoil.

Figure 4 .
Figure 4. Wind tunnel measurement section size and test airfoil.

Figure 7
Figure 7 shows compares this work with that by North China Electric Power University (NCEPU), Northwest Polytechnic University (NPU), and Delft University (Delft).It can be seen that in the linear section of the lift coefficient, the experimental data achieved in this work are similar to that obtained by NPU and Delft.After airfoil stall, the lift coefficient measured in this work is the highest compared with the results of the other three works.The lift coefficient first decreases and

Figure 7
Figure7shows compares this work with that by North China Electric Power University (NCEPU), Northwest Polytechnic University (NPU), and Delft University (Delft).It can be seen that in the linear section of the lift coefficient, the experimental data achieved in this work are similar to that obtained by NPU and Delft.After airfoil stall, the lift coefficient measured in this work is the highest compared with the results of the other three works.The lift coefficient first decreases and then increases slightly.Figure7bshows that the drag coefficient of the airfoil is low before the stall angle of attack (8 • ), and the drag coefficient measured in this work is similar to that measured by NPU and Delft.The drag coefficient measured after stall in this work is greater than that measured by NPU but similar to that measured by Delft.Because the aerodynamic performance of the airfoil is unstable after stall, the lift and drag coefficients have to be revised, and the mechanism of each correction method is different; however, our focus is on the aerodynamic performance before stall.

Energies 2019 ,
12, x FOR PEER REVIEW 8 of 19vortices generated by the VG will show a stronger trend, which leads to a lower pressure in the lowpressure zone on the leeward side of the VG.

Figure 10 .
Figure 10.Distribution of static pressure coefficient at vortex core center.

Figure 10 .
Figure 10.Distribution of static pressure coefficient at vortex core center.

Figure 11 .
Figure 11.Chordwise distribution of vortex core area and vortex circulation.

Figure 12 .
Figure 12.Distribution of vortex core area and vortex circulation with flow direction.

1 Figure 12 .
Figure 12.Distribution of vortex core area and vortex circulation with flow direction.

Figure 13 .
Figure 13.Change in vortex core area and vortex circulation at the VG trailing edge.

Figure 13 .
Figure 13.Change in vortex core area and vortex circulation at the VG trailing edge.

Figure 14 .
Figure 14.Inflow velocity profiles and average velocity in the VG height range.

Figure 14 .
Figure 14.Inflow velocity profiles and average velocity in the VG height range.

Figure 15 .
Figure 15.Lift and drag coefficient and lift-drag ratio with VG.

Figure 15 .
Figure 15.Lift and drag coefficient and lift-drag ratio with VG.

Figure 17 showsFigure 16 .
Figure 17 shows the axial velocity distribution for the normal height at the vortex core position downstream of the VG.The horizontal line is placed at the center of the vortex core.With the increase in the flow distance, the area of the vortex core varies at different VG heights.The height of the vortex core marked in the figure is the height of the vortex core when H = 1.0δ.

Figure 16 .
Figure 16.Axial velocity contours downstream of the VG.

Figure 17 .
Figure 17.Velocity profiles at the vortex core downstream of the VG.

Table 2 .
Ratio of VG Height to Different BL Thickness.
Height Ratio Case 1 Case 2 Case 3 Case 4 Case 5 Case 6

Table 2 .
Ratio of VG Height to Different BL Thickness.

Table 3 .
Results of the grid independence study.

Table 6 .
Relative variation of Cl, Cd, and lift-drag ratio.

Table 6 .
Relative variation of C l , C d , and lift-drag ratio.