This post is be about the validation and verification of the computational fluid dynamics analysis of a three blade vertical axis wind turbine. The turbine had a diameter of 2 m with each blade being 1 m tall. The blades had an NACA-0018 airfoil cross section.
The computational fluid dynamics analysis employed the κ-ε turbulence model with damping functions as the turbulence model, SIMPLE-R as the numerical algorithm. The spatial discretization schemes for the convective fluxes and diffusive terms used are the second order upwind and central approximations, respectively. An implicit first-order Euler scheme is employed to approximate the time derivatives.
The Cartesian computational mesh with immersed boundary method had a total of 769,357 cells. Among those 769,357 cells, 166,188 cells were around the turbine blades. Mesh controls were employed to refine the mesh near the turbine blades. A time step of 3e-3 was employed. The computational domain inlet was 1.5 D away from the turbine and the outlet was 3D away. The computational domain walls on the sides were 1D x 1.5D, where D represents the turbine diameter. The mesh and the computational domain are shown in Fig. 1. The vertical teal arrow represents the force of gravity, the curved teal arrow represents the direction of turbine rotation. The dark blue arrow represents the direction of free stream velocity.
The simulations ran at a tip-speed ratio of 1.87 at a wind speed of 4.03 m.s-1. The velocity distribution around the turbine after 4 revolutions is shown in Fig. 2. Validation of the numerical analysis was carried out using [1]. The results of power produced by the turbine were with in 4% of the experimental results [1]. An animation of the numerical analysis is also shown.
The computational fluid dynamics analysis employed the κ-ε turbulence model with damping functions as the turbulence model, SIMPLE-R as the numerical algorithm. The spatial discretization schemes for the convective fluxes and diffusive terms used are the second order upwind and central approximations, respectively. An implicit first-order Euler scheme is employed to approximate the time derivatives.
The Cartesian computational mesh with immersed boundary method had a total of 769,357 cells. Among those 769,357 cells, 166,188 cells were around the turbine blades. Mesh controls were employed to refine the mesh near the turbine blades. A time step of 3e-3 was employed. The computational domain inlet was 1.5 D away from the turbine and the outlet was 3D away. The computational domain walls on the sides were 1D x 1.5D, where D represents the turbine diameter. The mesh and the computational domain are shown in Fig. 1. The vertical teal arrow represents the force of gravity, the curved teal arrow represents the direction of turbine rotation. The dark blue arrow represents the direction of free stream velocity.
Fig. 1, Mesh and computational domain.
The simulations ran at a tip-speed ratio of 1.87 at a wind speed of 4.03 m.s-1. The velocity distribution around the turbine after 4 revolutions is shown in Fig. 2. Validation of the numerical analysis was carried out using [1]. The results of power produced by the turbine were with in 4% of the experimental results [1]. An animation of the numerical analysis is also shown.
Fig. 1, Flow field around the turbine.
Thank you for reading. If you would like to contribute to the research, both financially and scientifically, please feel free to reach out.
[1] Yi-Xin Peng, You-Lin Xu, Sheng Zhan and Kei-Man Shum, High-solidity straight-bladed vertical axis wind turbine: Aerodynamic force measurements, Journal of Wind Engineering and Industrial Aerodynamics, January 2019.
[1] Yi-Xin Peng, You-Lin Xu, Sheng Zhan and Kei-Man Shum, High-solidity straight-bladed vertical axis wind turbine: Aerodynamic force measurements, Journal of Wind Engineering and Industrial Aerodynamics, January 2019.
