Aerodynamic Performance of Vertical-Axis Wind Turbines

: The nonstationary separated incompressible ﬂows around Darrieus and Savonius rotors of vertical-axis wind turbines were investigated through computational simulation using the Reynolds averaged Navier–Stokes equations and Spalart–Allmaras turbulence model. The implicit ﬁnite-volume algorithm, the basis of which was artiﬁcial compressibility method, was chosen to obtain the numerical solution. The series of computational and physical experiments for Darrieus rotors with varied numbers and shapes of blades were performed. The detailed visualization of the ﬂow was presented. The turbulent ﬂows surrounding the Darrieus and Savonius rotors were studied, and as a part of these investigations, the major phases of vortex progress were identiﬁed. For this purpose, three series of computer tests on the aerodynamic and power properties of Savonius rotors with two and three buckets were performed, and their results are also presented. The inﬂuence of tip-speed ratio, solidity, and Reynolds numbers on the power coefﬁcients of the Darrieus and Savonius rotors was investigated. It has been demonstrated that increasing Reynolds number from 10 4 to 10 6 causes a rise in Darrieus rotors power coefﬁcient from 0.15 up to 0.5. The maximum values of power coefﬁcient are moved away from higher values of tip-speed ratio from 2 to 5 as a result of a decrease in Darrieus rotor solidity from 1.0 to 0.33. The greatest power coefﬁcient for a Savonius rotor with two blades is 0.23 and for a Savonius rotor with three blades is 0.19.


Introduction
A source of the cheapest electric power produced for commerce and industry in the near future could be wind energy, but the amount of the wind energy generation is still insufficient.Despite the obvious environmental benefits and constant cost reduction, generation of electricity from wind flows in the world is only 5%, which is several times lower than coal-and natural gas-based generation (38%) and (23%), respectively [1].
The biggest share of wind power generation (~85%) is contributed by horizontal-axis wind turbines (HAWT).HAWT has high power coefficients, but requires the use of complexto-manufacture blades with a variable profile and expensive electromechanical energy conversion devices.The need to reduce the per-unit costs of electricity production has led to the construction of large HAWT with a capacity of several megawatts [2].Placement of such wind turbines in offshore areas results in their distancing from the final energy consumer.Furthermore, the large size of the HAWT rotors blades (50 m and longer) leads to another problem-reduced service life of the blades due to constant lateral vibrations.
Obviously, aerodynamics is essential to a wind turbines operation.The wind turbine's ultimate efficiency is dependent on the rotor form optimization and aerodynamic properties of the spinning surfaces.
The wind turbine efficiency is defined by the coefficient of performance of a wind power.Results of experimental research showed that high-speed propeller type HAWT and also Darrieus [3] rotors with the fixed and variable pitch blades have high values of performance.
As reviewed by Eriksson et al. [4], vertical-axis wind turbines (VAWT) have a significant disadvantage compared to HAWTs, which is that the power coefficients achievable with this type of turbines are lower, with the power coefficients obtainable with HAWTs to be between 0.4 and 0.5 and the ones obtainable with VAWTs around 0.4.Despite this disadvantage, VAWTs provide numerous advantages compared to HAWTs.Therefore, VAWT could be adopted instead of HAWT.The Darrieus rotor blades have a squared shape and are cheaper to manufacture than the blades of the propeller type HAWT.Since VAWT are not influenced by the wind direction, electromechanical energy conversion devices become simpler, which reduce the operating costs.Because of the VAWT's ease of use, dependability, and capacity to operate in conditions of moderate wind, they can be deployed in remote locations of crucial social significance, see Paraschivoiu [5].
As mentioned previously, VAWTs perform less efficiently in terms of aerodynamics than HAWT.Such a situation is largely caused by a lack of understanding of the VAWT's aerodynamic processes.Another factor is the simplified methodology used to build industrial installations of high-power wind turbines, which ignores the latest achievements of the aerodynamics of the blades in the dynamic stall conditions and is based on empiric hypotheses about the flow pattern around the rotors of VAWTs.
Improving the VAWT designs can be achieved through the use of the Reynolds averaged Navier-Stokes equations (RANS), closed by differential turbulence models, together with the equations of the VAWT rotor dynamic and oscillations.Numerical simulation must be carried out within a wide range of key geometric, mass, and velocity parameters.Despite significant progress in using the methods of computational fluid dynamics (CFD) in solving industrial problems, systematic detailed studies of the flow around the VAWT rotors have actually not been carried out, as stated by Longhuan et al. [6].
Modern trends in design of complex machinery are related with application of complete mathematical models of mechanics of liquid and gas, based on general physical laws (conservation of mass, impulse, energy, rheological relations, dynamics of turbulent whirlwinds).From the mathematical point of view, such models are complicated systems of nonlinear differential equations, which require powerful calculation complexes.For creation of such models, many sections of theoretical and applied mathematics are utilized, namely, analytical and differential geometry, mathematical and, in particular, tensor analysis, mathematical physics, phenomenological principles of modeling, numerical analysis, and computational methodology.Solving of such systems creates a new level of learning when numeral calculations may replace natural experiments, substantially cutting down expenses on tests in wind-channels, which are present in a standard design procedure.
There are four levels of mathematical models of aerodynamics of wind turbines: I. Engineering, empirical technique; II.Inviscid gas model-The potential and Euler equations; III.Viscous gas model with averaging turbulent characteristics (RANS equations); IV.Viscous gas models with unsteady turbulent characteristics (DES, LES, DNS approaches).
Nowadays, the first two levels are crucial to wind power design.The specification of the aerodynamic, dynamic, and power properties of wind turbine rotors is commonly conducted using methods that are based on experimental data, such as in Brochier et al. [7] and Hu et al. [8], impulse and vortex theories, as in Strickland [9], and the computational solution of the potential and Euler equations.These methods rely among others on quasi-steady flow and ignorance viscous-inviscid interaction assumptions.The impact of dynamic stall on the computation of nonstationary stream around rotors of wind turbine is the main challenge, see Biadgo et al. [10].Until now, simplified models have provided aerodynamic estimates with insufficient accuracy, particularly in the situation of dynamic stall.
Passing to the third level of simulation of aerodynamic flows, it is possible to state that calculations on the basis of RANS equations are becoming widely used in industrial practice.Sun et al. [11] studied operating variations of a Darrieus VAWT by means of RANS equations, and Gonzalez Madina et al. [12] studied that of a Savonius [13] VAWT.
The fourth level is research in computational fluid dynamics.It is rather expensive even for industrialized countries, and its introduction in the industry is predicted for second half of the 21st century.Nevertheless, some authors have used the approaches mentioned above for VAWT analysis.For example, Abkar and Dabiri [14] used LES for flow characteristics prediction of the wake behind a VAWT, and Nguyen et al. [15] used DNS methods.Sheidani et al. [16] used both LES and RANS models for predicting the wake behind a VAWT, showing that RANS methods provide accurate predictions of the wake region length and width, and therefore, they are suitable for obtaining an overall view of the wake region.However, LES is able to predict more flow structures than RANS, especially in the internal wake region.
In recent works, Crooks et al. [17] conducted a computational design analysis study of a counter-rotating Darrieus turbine system, Karimian and Abdolahifar [18] proposed a new Darrieus design to benefit from straight-blade and helical-blade VAWTs, and Dadamoussa et al. [19] added Vortex Generators on the blades of a Darrieus VAWT, improving its performance.Al-Gburi et al. [20] compared the influence of several parameters on the performance of a Savonius rotor, and Altan and Gultekin [21] proposed additional design processes for Savonius rotor performance enhancement.
The aim of this paper is to study the structure and determine the main energy parameters of the unsteady flow around vertical-axis wind turbines with Darrieus and Savonius rotors based on the developed specialized CFD code.
The novelty and main contributions of this paper: 1.
The regularities of the development of the unsteady flow around Darrieus and Savonius rotors with different number of blades have been revealed.

2.
Qualitative comparisons of the flow pattern around a two-and three-blade Darrieus rotor with experimental data were carried out.

3.
The dependences of the power coefficient of the VAWT rotors on the Reynolds number, tip-speed ratio, and solidity are determined.

4.
The influence of the Reynolds number on the dynamic stall from the blades and the Darrieus rotor power coefficient is determined.

5.
The influence of the solidity and the tip-speed ratio on the power characteristics of the Darrieus rotor is determined.6.
The influence of the blades number and Reynolds number on the power characteristics of the Savonius rotor is determined.
The following manuscript is divided as follows: Section 2 explains the tested cases and numerical methodology, Section 3 details the numerical algorithm used for the simulations and describes the developed CFD code, Section 4 shows the obtained results, and Section 5 summarizes the main findings of this study.

Methodology
The orthogonal Darrieus and Savonius rotors are investigated in this study (Figure 1).The normal blades are substantially longer than the chord, which enables one to ignore end effects on the blades and to use the plane-parallel structure of flow hypothesis.As a result, the study assumes a two-dimensional statement in a plane that is perpendicular to the rotor's rotational axis.The rotors of Darrieus and Savonius are considered to be completely rigid.Incompressible RANS equations closed by the Spalart-Allmaras (SA) eddy viscosity turbulence model are used to describe the aerodynamic processes of wind turbines.These equations are straightforward, trustworthy, and suitable for a variety of external flows.For example, Shukla and Kaviti [22] demonstrated the reliability of this model for a VAWT with different airfoils.The normal blades are substantially longer than the chord, which enables one to ignore end effects on the blades and to use the plane-parallel structure of flow hypothesis.As a result, the study assumes a two-dimensional statement in a plane that is perpendicular to the rotor's rotational axis.The rotors of Darrieus and Savonius are considered to be completely rigid.Incompressible RANS equations closed by the Spalart-Allmaras (SA) eddy viscosity turbulence model are used to describe the aerodynamic processes of wind turbines.These equations are straightforward, trustworthy, and suitable for a variety of external flows.For example, Shukla and Kaviti [22] demonstrated the reliability of this model for a VAWT with different airfoils.

Mathematical Model
The low-speed air movement at low Mach numbers (M < 0.3) is examined in this work.Since the effects of air compressibility may be ignored in this situation, the simulation of the turbulent flow under investigation has been conducted using the incompressible RANS equations for a viscous flow (in Cartesian tensor form) [23] where M = u/u∞ is Mach number; u∞ is the free stream velocity, which used the definition of the Reynolds number too; ui are the Cartesian components of the velocity vector; xi, i = 1, 2 are the Cartesian coordinates, m; t is time, s; ρ is density, kg/m 3 ; m/s; p is pressure, Pa; ν and νt are the kinematic coefficients of molecular and turbulent viscosity, m 2 /s, correspondingly.

Turbulence Model
The RANS equations are closed by using the differential single-parametric SA turbulence model [18].The goal of the conventional SA model is to calculate the turbulent viscosity's dimensional kinematic coefficient where fv1 is a damping function of kinematic viscosities χ , t ṽ is the working variable, . In the SA model the equation to determine t ṽ takes the form [24]

Mathematical Model
The low-speed air movement at low Mach numbers (M < 0.3) is examined in this work.Since the effects of air compressibility may be ignored in this situation, the simulation of the turbulent flow under investigation has been conducted using the incompressible RANS equations for a viscous flow (in Cartesian tensor form) [23] where M = u/u ∞ is Mach number; u ∞ is the free stream velocity, which used the definition of the Reynolds number too; u i are the Cartesian components of the velocity vector; x i , i = 1, 2 are the Cartesian coordinates, m; t is time, s; ρ is density, kg/m 3 ; m/s; p is pressure, Pa; ν and ν t are the kinematic coefficients of molecular and turbulent viscosity, m 2 /s, correspondingly.

Turbulence Model
The RANS equations are closed by using the differential single-parametric SA turbulence model [18].The goal of the conventional SA model is to calculate the turbulent viscosity's dimensional kinematic coefficient where f v1 is a damping function of kinematic viscosities χ, v t is the working variable, m 2 /s, c v1 = 7.1.In the SA model the equation to determine v t takes the form [24] The source term for turbulence generation is the first term on the equation's right-hand side (4) where W ij is the vorticity tensor, 1/s, k = 0.41.
The function f v2 is determined from the relation The turbulence dissipation is caused by the second and third terms in the right-hand section of Equation ( 4).The fourth one has a function and is the dissipation of turbulence nearby a solid wall.
where d is distance from the nearest wall, m.The values for other non-dimensional constants are [24] The traditional terms of SA model, containing f t1 , f t2 and responsible for near-wall treatment for adequate reproduction of laminar-turbulent transition behavior, are not presented in the applied model version due to the fact that the laminar flow zones are negligible in comparison with the prevailing turbulent flow mode.

Boundary Conditions
No-slip boundary conditions are set on the body surfaces.By defining a zero-pressure gradient normal to the wall, the pressure is obtained.To determine the dependent variables at the outer boundary of the computational domain, non-reflecting boundary conditions based on the method of characteristics [25,26] are used.

Numerical Algorithm
The coordinate system for the governing equations was arbitrarily curvilinear.To calculate nonstationary flows, the artificial compressibility approach was used to agree to the pressure and velocity fields.The basic equations were converted into discrete forms using the structured grids.The multi-block computation methods were applied in places that were not just directly connected (Figure 2).Each block in an intersecting grid and may have a variable size.This strategy made it possible to provide a standardized system for estimating viscous fluid flows around objects with complex geometries.
The calculation grids were built by the method of many surfaces [27].The outer boundary is a circle with a radius of 40.All grids (domains) were rigidly connected to each other and rotated synchronously with the rotor blades relative to the center.
The dimension of the grids (the number of nodes in each direction) is not explicitly related to each other.The transfer of information from one computational domain to another is performed through their mutual overlap by means of interpolation relations linking dependent variables in the nodes of neighboring grids.
The computational grid around the Savonius rotor was a single-block O-type mesh with a total number of nodes of about one million (1050 × 980) (Figure 2a).A multi-block mesh with a total number of nodes of about 600 thousand was used around the Darrieus rotor (Figure 2b).A separate O-type mesh was built around each element of the Darrieus rotor, which intersected with neighboring meshes by at least two nodes.The first grid spacing near the surface was 10 −6 , which ensured the fulfillment of the condition y + < 1 for all computational cases.The computational grid around the Savonius rotor was a single-block O-type mesh with a total number of nodes of about one million (1050 × 980) (Figure 2a).A multi-block mesh with a total number of nodes of about 600 thousand was used around the Darrieus rotor (Figure 2b).A separate O-type mesh was built around each element of the Darrieus rotor, which intersected with neighboring meshes by at least two nodes.The first grid spacing near the surface was 10 −6 , which ensured the fulfillment of the condition y + < 1 for all computational cases.
In order to preserve the third order of approximation of convective terms in the transport equations, each node should have at least two neighboring nodes on each side when using a single-block grid.The order of approximation of convective terms can be decreased to a second one only for nodes located near the body surface or outer boundary of the computational domain.Based on these considerations, the minimum overlapping area of the computational grids should contain at least four nodes.Bicubic interpolation is used to determine the dependent variables at the inner boundaries of the computational blocks.
Using the finite volume method, the initial equations were numerically integrated.For convective terms, the third-order upwind Rogers-Kwak technique [28] was utilized.Convective terms in the turbulence model were approximated using the TVD scheme and a third-order ISNAS (Interpolation Scheme which is Non-oscillatory for Advected Scalars) flow limiter.A second-order central difference technique was used to approximate the derivatives in terms of viscosity.The algorithm for solving the equations is based on a three-level second-order time-accurate implicit system with sub-iterations in pseudo-time.The generalized minimal residual (GMRES) approach with ILU(k) pre-conditioning was employed to solve the resultant block-matrix system of linear algebraic equations.In order to preserve the third order of approximation of convective terms in the transport equations, each node should have at least two neighboring nodes on each side when using a single-block grid.The order of approximation of convective terms can be decreased to a second one only for nodes located near the body surface or outer boundary of the computational domain.Based on these considerations, the minimum overlapping area of the computational grids should contain at least four nodes.Bicubic interpolation is used to determine the dependent variables at the inner boundaries of the computational blocks.
Using the finite volume method, the initial equations were numerically integrated.For convective terms, the third-order upwind Rogers-Kwak technique [28] was utilized.Convective terms in the turbulence model were approximated using the TVD scheme and a third-order ISNAS (Interpolation Scheme which is Non-oscillatory for Advected Scalars) flow limiter.A second-order central difference technique was used to approximate the derivatives in terms of viscosity.The algorithm for solving the equations is based on a three-level second-order time-accurate implicit system with sub-iterations in pseudo-time.The generalized minimal residual (GMRES) approach with ILU(k) pre-conditioning was employed to solve the resultant block-matrix system of linear algebraic equations.

Developed CFD Code
The author's team used the specific CFD program created by Redchyts [29] in the context of this work.This software makes it possible to obtain the required level of compromise between the required processing resources and the quality of the results.The established computational methodology's foundational framework offers a comprehensive CFD approach based on the Navier-Stokes equations, including the RANS technique, a number of differential turbulence models, and a multi-block method for the flows in the multiple intersected domains.
On the one hand, a complete computational fluid dynamics approach based on the Navier-Stokes equations was implemented.On the other hand, the analysis of the physics of flow around orthogonal Darius and Savonius rotors allowed us to reduce the problem statement to a two-dimensional plane (Figure 1).This made it possible to reduce computational costs without compromising the reliability of the numerical experiment and to bring the calculation of unsteady flow to 50-100 h on a standard personal computer.
The developed specialized CFD package is written in Fortran 90 and allows solving two-dimensional unsteady laminar and turbulent flows of liquid and gas in arbitrary geometry.The package generates meshes using the multi-surface method [27].It is also possible to import grids from external mesh generators.
The development of CFD tools made it possible to handle problems in dynamics and aerodynamics as well as electrochemistry, multiphase fluids, combustion, and plasma kinetics.The fresh technical concepts developed, and new knowledge about the physics of flow separation and how to control it, enable to replicate the actual CFD code structure to many flow types over a wide speed range.
The pre-processor, the computing core, and the postprocessor are the three standard primary components of the generated specialized CFD code, just as in other comparable commercial and opensource software.The development of the basic geometry, the production of discrete space, and the movement of specific computational grid blocks are all intended uses for the pre-processor.In the computing core, the equations of dynamics of a viscous incompressible medium and a solid are integrated, the turbulence parameters are calculated, the boundary conditions are met, and the integral characteristics of the VAWT rotor are determined as functions in time.The post-processor is used to process integral flow parameters over time (such as the calculation of aerodynamic characteristics), prepare data for external visualizers, and perform built-in visualization of the findings (scalar and vector values).Benefits of house-developed package: the main task of the specialized package is to simulate the flow field numerically with maximum physical veracity; -specialization of the package and the efficiency of solving this class of problems; -no need to buy or rent expensive commercial packages.
This package has been developed since 2005 and was initially focused on solving the problems of unsteady flow around bodies with moving boundaries.The pre-and postprocessing tools, including the solver, are purposefully designed for maximum convenience and efficiency while being inexpensive for priority application to these flows.In addition, the presence of a package of our own design allows us to develop it in the future in any desired direction, depending on the tasks that will be relevant to the solution.This approach, for all its cost in terms of writing and debugging software code, makes it independent of the general trends in the CFD software market.
Modern verification methods and careful testing of the developed programs are absolutely necessary to verify the obtained numerical results and the reliability of the CFD methods.To verify the developed numerical methods, the test calculations of a series of problems such as laminar and turbulent flow around a stationary and rotating cylinder [29], subcritical and supercritical flow around stable and oscillating airfoils NACA 4412 and NACA 0015, aerodynamics of turbulent flow around a 30P30N three-element airfoil in a crash and landing configuration [30], and turbulent flow around a cargo truck-trailer [31] have been performed.

Darrieus Rotor
The results of numerical simulation of the flow around a rotating single blade of the Darrieus rotor are presented (Figure 3).The initial data for the numerical simulation were acquired from an experimental work [32] for the case of flow around a single Darrieus rotor blade in the presence of dynamic stall.At the tip-speed ratio λ = 2.5 and Reynolds number Re = 67,000, dynamic separation of the flow from the blade of the Darrieus rotor is observed over most of the trajectory.It is characterized by the flow separation from the leading edge of the blade and the formation of large vortex structures that are carried along the chord of the blade.The absolute velocity change along the circumference of the rotation blade leads to a larger area of dynamic stall in that part of the trajectory where the blade and flow velocities are in the same direction.In this zone, the incident stream carries the vortices in the direction of the blade movement.On the second half of the trajectory, the flow carries the vortices in the direction opposite to the movement of the blade.In this case, the duration of the dynamic stall is less than in the previous one.
leading edge of the blade and the formation of large vortex structures that are carried along the chord of the blade.The absolute velocity change along the circumference of the rotation blade leads to a larger area of dynamic stall in that part of the trajectory where the blade and flow velocities are in the same direction.In this zone, the incident stream carries the vortices in the direction of the blade movement.On the second half of the trajectory, the flow carries the vortices in the direction opposite to the movement of the blade.In this case, the duration of the dynamic stall is less than in the previous one.The value of the angle θ = 0° corresponds to the rotor position when the leading edge of the blade is directed towards the flow.On the windward part of the blade trajectory, the blade angle increases, the pressure coefficient on the inner surface of the blade decreases.This leads to an increase in the absolute values of the lift coefficient, which makes the main contribution to the torque (Figure 4).The value of the angle θ = 0 • corresponds to the rotor position when the leading edge of the blade is directed towards the flow.On the windward part of the blade trajectory, the blade angle increases, the pressure coefficient on the inner surface of the blade decreases.This leads to an increase in the absolute values of the lift coefficient, which makes the main contribution to the torque (Figure 4).The boundary layer separated near the trailing edge of the blade, and the position of the separation point is shifted towards the leading edge of the blade.This leads to the separation of the vortex from the leading edge of the blade and subsequent movement along the chord towards the trailing edge of the blade.A drop in the lift and drag coefficients is observed (Figure 4).On the leeward part of the blade trajectory, the streamline pattern and aerodynamic coefficient behavior is somewhat different from the windward part.The boundary layer separated near the trailing edge of the blade, and the position of the separation point is shifted towards the leading edge of the blade.This leads to the separation of the vortex from the leading edge of the blade and subsequent movement along the chord towards the trailing edge of the blade.A drop in the lift and drag coefficients is observed (Figure 4).On the leeward part of the blade trajectory, the streamline pattern and aerodynamic coefficient behavior is somewhat different from the windward part.
The dependence of non-averaged coefficients of lift force C L , drag force C D , and torque C Q on the angular position of the rotor are shown in Figure 4.At the tip-speed ratio λ = 2.5, the maximum of torque coefficient C Q = 0.4 is observed at rotor angle θ = 80 • .After reaching the maximum of torque coefficient, it drops to values of C Q = −0.2,which corresponds to the rotor angle θ = 145 • .In this range of angular positions, the Darrieus rotor is subject to dynamic stall.Subsequently, when flow attachment θ = 190 • occurs, the value of torque coefficient has a positive value up to angle θ = 250 • .Further increasing of angle leads to dynamic stall of flow from the outside surface of the blade, and torque coefficient becomes negative (Figure 4).
The values of the torque and power coefficients of a single blade of the Darrieus rotor averaged over one revolution have the following values C Q = 0.11 and C P = 0.044, respectively.
Only symmetric airfoils from NACA00XX family are considered in the majority of experimental and theoretical studies [7][8][9], analyzing the flow structures of the Darrieus rotor and its power coefficients.The overall structure solidity may vary if one or several blades are added.In present paper, the VAWT with NACA 0018 airfoil blades [10,11] at zero angle of attack is studied.
For a two-blade Darrieus rotor, flow visualization and reconstruction have been carried out using physical [7] and computational experiments (Figure 5).Only the vortices with the highest intensity are shown for better clarity.At various points along the trajectory of the blade, the stages of vortex creation, development, and stall are presented.The flow around the blade has an associated character at the start of a leeward side of a trajectory and before rotor angle of 234 degrees.At the rotor angle (240 • ) specified, a separated flow is formed.Vortices split from the blade leading edge and start moving along a surface when there is dynamic stall on the inner surface [32][33][34].
The calculations were carried out in two stages.In the first stage, an undisturbed flow was set in the computational domain.In the second stage, the flow structure was studied at a fixed tip-speed ratio, and the main parameters of the unsteady flow of the Darrieus rotor were determined.
The three-blade Darrieus rotor's flow pattern is presented by the physical investigation [35] (Figure 6a).Numerical simulation was carried out with tip-speed ratio of λ = 3 (Figure 6b).Dynamic separation of the flow from the blade of the Darrieus rotor is observed on most of the trajectory.A series of large, counter-rotating vortices make up the instantaneous flow structure.The main reason for the dominance of dynamic vortex shedding over most of the blade trajectory is the low Reynolds number Re = 1470, which corresponds to the initial stage of the transition from laminar to turbulent flow.Different areas of the blade trajectory exhibit asymmetries.The majority of shed vortices are more intense than the vortices traveling with the flow.This is as a result of the relative flow velocity being higher at this point in the blade trajectory than it is at the opposite point.A cluster of vortices that were created by the first blade is present in the front portion of the Darrieus rotor.There is an accumulation of vortices created by earlier blades downstream.
corresponds to the initial stage of the transition from laminar to turbulent flow.Different areas of the blade trajectory exhibit asymmetries.The majority of shed vortices are more intense than the vortices traveling with the flow.This is as a result of the relative flow velocity being higher at this point in the blade trajectory than it is at the opposite point.A cluster of vortices that were created by the first blade is present in the front portion of the Darrieus rotor.There is an accumulation of vortices created by earlier blades downstream.For the parameters under consideration, the stream pattern is distinguished by fundamentally unstable effects, such as stalling, the creation of intricate vortex ensembles, and interactions between vortices of various intensities, sizes, and speeds and the rotor surface.The impact of tip-speed ratio, solidity, and Reynolds numbers on the power coefficients of the Darrieus rotors was investigated (Figure 7).The main parameters of computational experiments: For the parameters under consideration, the stream pattern is distinguished by fundamentally unstable effects, such as stalling, the creation of intricate vortex ensembles, and interactions between vortices of various intensities, sizes, and speeds and the rotor surface.The impact of tip-speed ratio, solidity, and Reynolds numbers on the power coefficients of the Darrieus rotors was investigated (Figure 7).The main parameters of computational experiments:

Savonius Rotor
This paper presents the findings from three sets of simulations on the aerodynam and power properties of Savonius rotors with two and three blades.In [35], an expe mental study of a two-and three-bladed Savonius rotor in a low-speed wind tunnel w carried out in order to determine the aerodynamic and energy characteristics.It has been demonstrated that increasing Reynolds number from 10 4 to 10 6 causes a rise in power coefficient from 0.15 up to 0.5.The maximum values of power coefficient are moved away from higher values of tip-speed ratio from 2 to 5 as a result of a decrease in Darrieus rotor solidity from 1.0 to 0.33 [32,33].
The overestimated values of the Darrieus rotor power factor by 10-20% are the result of a somewhat idealized mathematical model.The developed model does not take into account the three-dimensional effects, end effects, large-scale turbulence, the effects of the laminar-turbulent transition, and rotor power losses due to the electro-mechanical part of the rotor.Accounting for these effects will be the focus of our further research.

Savonius Rotor
This paper presents the findings from three sets of simulations on the aerodynamic and power properties of Savonius rotors with two and three blades.In [36], an experimental study of a two-and three-bladed Savonius rotor in a low-speed wind tunnel was carried out in order to determine the aerodynamic and energy characteristics.The undisturbed flow velocity was U ∞ = 7 m/s and U ∞ = 14 m/s, which corresponds to the Reynolds numbers calculated from the rotor diameter Re = 4.32 × 10 5 and Re = 8.64 × 10 5 , respectively.
The geometry of a two-blade Savonius rotor.The rotor blades have semicircles with radius r = 0.25 m in their cross section [36].The rotation radius was measured from the rotation axis to the outer edge of the blades.The computational grid contained 800 × 600 nodes.
The geometry of a three-blade Savonius rotor.The arc of the blades in the Savonius threebladed rotor was reduced by 20 at the outer edge [36].The calculated grid contained 1050 × 980 ≈ 10 6 nodes.A detailed grid is necessary to adequately represent the complex geometry of the Savonius rotor and to qualitatively account for the processes in the turbulent boundary layer.
The duration of one computational experiment on a personal computer was about 100 h.
A stable Savonius rotor that was fixed at different angles of attack was the subject of the first set of computational experiments.The time-averaged torque coefficients for the two-blade Savonius rotor are positive for most angular positions, but they are negative for angles between 55 and 80 degrees.Torque coefficients for three-blade designs are always positive.
The Savonius rotor's fixed tip-speed ratio was the subject of the second set of simulations (Figure 8).The maximum torque coefficient for Savonius rotors with two blades is 0.4, while for those with three blades, it is 0.35.(Figure 9).They match a tip-speed ratio of 0.40.The greatest power coefficient for a Savonius rotor with two blades is 0.23 and for a Savonius rotor with three blades is 0.19.The obtained results and the experimental data are adequately matched [32,33].A stable Savonius rotor that was fixed at different angles of attack was the subject of the first set of computational experiments.The time-averaged torque coefficients for the two-blade Savonius rotor are positive for most angular positions, but they are negative for angles between 55 and 80 degrees.Torque coefficients for three-blade designs are always positive.
The Savonius rotor's fixed tip-speed ratio was the subject of the second set of simulations (Figure 8).The maximum torque coefficient for Savonius rotors with two blades is 0.4, while for those with three blades, it is 0.35.(Figure 9).They match a tip-speed ratio of 0.40.The greatest power coefficient for a Savonius rotor with two blades is 0.23 and for a Savonius rotor with three blades is 0.19.The obtained results and the experimental data are adequately matched [32,33].The solution to the joined problem of the aerodynamics and dynamics of the threeblade Savonius rotor is represented by the third series of computational experiments.Three phases of calculations were performed.The first stage's goal was to identify a periodic stream with a Karman vortex street-like shape (Figure 10a).The coupled equations governing the rotation and aerodynamics of Savonius rotors were solved in the second and third stages (Figure 10a,b).When the rotor is unlocked, the incoming wind flow causes it to start rotating, intensifying the vortex formation.Incoming stream speed, distinctive size, and VAWT rotation velocity all affect how often vortices were detached.The Savonius rotor was subjected to load torque in the third stage (Figure 10c), which stabilized the rotor's angular speed and caused close to cyclic fluctuations in the drag, lift, and torque coefficients [32,33].The solution to the joined problem of the aerodynamics and dynamics of the three-blade Savonius rotor is represented by the third series of computational experiments.Three phases of calculations were performed.The first stage's goal was to identify a periodic stream with a Karman vortex street-like shape (Figure 10a).The coupled equations governing the rotation and aerodynamics of Savonius rotors were solved in the    ) leads to an increase of the Darrieus rotor power coefficient (from 0.1 to 0.5).It was found that a decrease of the Darrieus rotor solidity (from 1 to 0.33), the power coefficient becomes less sensitive to changes of the tip-speed ratio.It is shown that the optimal values of the tip-speed ratio (from 2 to 5) significantly depend on the Reynolds number and the Darrieus rotor solidity.On the basis of computational and natural experiments, it is shown that the maximum power coefficient of a two-blade Savonius rotor (0.23) is higher than that of a three-blade rotor (0.19).
The developed specialized CFD code can be used in the design and development of new wind turbines and selection of a rational rotor shape and wind turbine operating modes.The obtained parametric dependencies can be used to improve existing engineering and semi-empirical methods.
The direction of further research will be the application of the γ-Reθ transition model for this class of flows together with the SA turbulence model and the comparison of the obtained results with experimental and computational data of other authors.It is also planned to develop a specialized package for solving three-dimensional incompressible Navier-Stokes equations for the calculation of arbitrarily shaped bodies.

Figure 3 .
Figure 3. Reconstruction of the flow field around one-blade Darrieus rotor.

Figure 3 .
Figure 3. Reconstruction of the flow field around one-blade Darrieus rotor.

Figure 4 .
Figure 4. Change of lift CL, drag CD, and torque CQ coefficients of the blade depending on the rotor angular position θ.

Figure 4 .
Figure 4. Change of lift C L , drag C D , and torque C Q coefficients of the blade depending on the rotor angular position θ.

Figure 5 .
Figure 5. Visualization (a,c) and reconstruction (b,d) of the flow around Darrieus rotor of two blades based on experimental [7] (a,b) and computational (c,d) results.Figure 5. Visualization (a,c) and reconstruction (b,d) of the flow around Darrieus rotor of two blades based on experimental [7] (a,b) and computational (c,d) results.

Figure 5 .Figure 6 .
Figure 5. Visualization (a,c) and reconstruction (b,d) of the flow around Darrieus rotor of two blades based on experimental [7] (a,b) and computational (c,d) results.Figure 5. Visualization (a,c) and reconstruction (b,d) of the flow around Darrieus rotor of two blades based on experimental [7] (a,b) and computational (c,d) results.J. Mar.Sci.Eng.2023, 11, x FOR PEER REVIEW 11 of 17

Figure 6 .
Figure 6.Visualization of the flow field during operation of the three-blade Darrieus rotor for the speed coefficient λ = 3 based on physical (a) and computational (b) experiments [35].

Figure 8 .
Figure 8. Flow fields around the rotating Savonius rotor with two (a) and three (b) blades, respectively, based on computer experiments.

Figure 8 .
Figure 8. Flow fields around the rotating Savonius rotor with two (a) and three (b) blades, respectively, based on computer experiments.

Figure 10 .
Figure 10.Flow fields around the non-rotating (a) and rotating (b,c) three-blade Savonius roto 5. Conclusions Using the developed specialized CFD code based on the Reynolds averaged vier-Stokes equations with turbulence model, new solutions for the unsteady around the rotors of the VAWT with Darrieus and Savonius rotors are performed finite volume approach was used to numerically integrate the set of governing equat Convective terms were approximated to third-order accuracy using the Rogers-K counter-flow scheme.The turbulence model's convective component was simulate ing a TVD scheme with a third-order ISNAS flow limiter.Two different varieties of rieus and Savonius rotors, with different numbers of blades and geometries, have properties estimated and researched.The investigation of the growth of turbulen stable incompressible flows around the VAWT Darrieus and Savonius rotors was ducted.Qualitative comparisons of the flow pattern around a two-and three-blade rieus rotor with experimental data have been carried out.The regularities of unst flow development around Darrieus and Savonius rotors with different number of b have been revealed.The dependences of the power coefficient of the VAWT rotors o Reynolds number, tip-speed ratio, and solidity are determined.It is shown that th namic stall from the blades of the Darrieus rotor at low Reynolds numbers is observe most of the blade trajectory.It is established that an increase of the Reynolds nu
Using the developed specialized CFD code based on the Reynolds averaged Navier-Stokes equations with turbulence model, new solutions for the unsteady flow around the rotors of the VAWT with Darrieus and Savonius rotors are performed.The finite volume approach was used to numerically integrate the set of governing equations.Convective terms were approximated to third-order accuracy using the Rogers-Kwak counter-flow scheme.The turbulence model's convective component was simulated using a TVD scheme with a third-order ISNAS flow limiter.Two different varieties of Darrieus and Savonius rotors, with different numbers of blades and geometries, have their properties estimated and researched.The investigation of the growth of turbulent unstable incompressible flows around the VAWT Darrieus and Savonius rotors was conducted.Qualitative comparisons of the flow pattern around a two-and three-blade Darrieus rotor with experimental data have been carried out.The regularities of unsteady flow development around Darrieus and Savonius rotors with different number of blades have been revealed.The dependences of the power coefficient of the VAWT rotors on the Reynolds number, tip-speed ratio, and solidity are determined.It is shown that the dynamic stall from the blades of the Darrieus rotor at low Reynolds numbers is observed on most of the blade trajectory.It is established that an increase of the Reynolds number (from 10 4 to 10 6