Evaluation of Turbulence Models for the Prediction of Wind Turbine Aerodynamics

Author(s):  
Sarun Benjanirat ◽  
Lakshmi N. Sankar ◽  
Guanpeng Xu

The performance of the NREL Phase VI horizontal axis wind turbine has been studied with a 3-D unsteady Navier-Stokes solver. This solver is third order accurate in space and second order accurate in time, and uses an implicit time marching scheme. Calculations were done for a range of wind conditions from 7 m/s to 25 m/s where the flow conditions ranged from attached flow to massively separated flow. A variety of turbulence models were studied: Baldwin-Lomax Model, Spalart-Allmaras one-equation model, and k-ε two equations model with and without wall corrections. It was found all the models predicted the normal forces and associated bending moments well, but most of them had difficulties in modeling the chord wise forces, power generation, and pitching moments. It was found that the k-ε model with near wall corrections did the best job of predicting most the quantities with acceptable levels of accuracy. Additional studies aimed at transition model development, and grid sensitivity studies in the tip region are deemed necessary to improve the correlation with experiments.

2021 ◽  
pp. 1-13
Author(s):  
Khaoula Qaissi ◽  
Omer A Elsayed ◽  
Mustapha Faqir ◽  
Elhachmi Essadiqi

Abstract A wind turbine blade has the particularity of containing twisted and tapered thick airfoils. The challenge with this configuration is the highly separated flow in the region of high twist. This research presents a numerical investigation of the effectiveness of a Vortex Trapping Cavity (VTC) on the aerodynamics of the National renewable Energy laboratory (NREL) Phase VI wind turbine. First, simulations are conducted on the S809 profile to study the fluid flow compared to the airfoil with the redesigned VTC. Secondly, the blade is simulated with and without VTC to assess its effect on the torque and the flow patterns. The results show that for high angles of incidence at Rec=106, the lift coefficient increases by 10% and the wake region appears smaller for the case with VTC. For wind speeds larger than 10 m/s, the VTC improves the torque by 3.9%. This is due to the separation that takes place in the vicinity of the VTC and leads to trapping early separation eddies inside the cell. These eddies roll up forming a coherent laminar vortex structure, which in turn sheds periodically out of the cell. This phenomenon favourably reshapes excessive flow separation, reenergizes the boundary layer and globally improves blade torque.


Author(s):  
Earl P. N. Duque ◽  
Michael D. Burklund ◽  
Wayne Johnson

A vortex lattice code, CAMRAD II, and a Reynolds-Averaged Navier-Stoke code, OVERFLOW-D2, were used to predict the aerodynamic performance of a two-bladed horizontal axis wind turbine. All computations were compared with experimental data that was collected at the NASA Ames Research Center 80-by 120-Foot Wind Tunnel. Computations were performed for both axial as well as yawed operating conditions. Various stall delay models and dynamics stall models were used by the CAMRAD II code. Comparisons between the experimental data and computed aerodynamic loads show that the OVERFLOW-D2 code can accurately predict the power and spanwise loading of a wind turbine rotor.


2005 ◽  
Vol 127 (4) ◽  
pp. 475-487 ◽  
Author(s):  
Sven Schmitz ◽  
Jean-Jacques Chattot

A Navier-Stokes solver, CFX V5.6, is coupled with an in-house developed Vortex-Panel method for the numerical analysis of wind turbines. The Navier-Stokes zone is confined to the near-field around one wind turbine blade, the Vortex-Panel method models the entire vortex sheet of a two-bladed rotor and accounts for the far-field. This coupling methodology reduces both numerical diffusion and computational cost. The parallelized coupled solver (PCS) is applied to the NREL Phase VI rotor configuration under no-yaw conditions. Fully turbulent flow is assumed using the k - ϵ and k - ω turbulence models. Results obtained are very encouraging for fully attached flow. For separated and partially stalled flow, results are in good agreement with experimental data. Discrepancies observed between the turbulence models are attributed to different prediction of the onset of separation. This is revealed by two-dimensional (2D) results of the S809 airfoil.


2008 ◽  
Vol 130 (3) ◽  
Author(s):  
Alvaro Gonzalez ◽  
Xabier Munduate

This work undertakes an aerodynamic analysis over the parked and the rotating NREL Phase VI wind turbine blade. The experimental sequences from NASA Ames wind tunnel selected for this study respond to the parked blade and the rotating configuration, both for the upwind, two-bladed wind turbine operating at nonyawed conditions. The objective is to bring some light into the nature of the flow field and especially the type of stall behavior observed when 2D aerofoil steady measurements are compared to the parked blade and the latter to the rotating one. From averaged pressure coefficients together with their standard deviation values, trailing and leading edge separated flow regions have been found, with the limitations of the repeatability of the flow encountered on the blade. Results for the parked blade show the progressive delay from tip to root of the trailing edge separation process, with respect to the 2D profile, and also reveal a local region of leading edge separated flow or bubble at the inner, 30% and 47% of the blade. For the rotating blade, results at inboard 30% and 47% stations show a dramatic suppression of the trailing edge separation, and the development of a leading edge separation structure connected with the extra lift.


2002 ◽  
Vol 124 (2) ◽  
pp. 413-423 ◽  
Author(s):  
L. S. Hedges ◽  
A. K. Travin ◽  
P. R. Spalart

The flow around a generic airliner landing-gear truck is calculated using the methods of Detached-Eddy Simulation, and of Unsteady Reynolds-Averaged Navier-Stokes Equations, with the Spalart-Allmaras one-equation model. The two simulations have identical numerics, using a multi-block structured grid with about 2.5 million points. The Reynolds number is 6×105. Comparison to the experiment of Lazos shows that the simulations predict the pressure on the wheels accurately for such a massively separated flow with strong interference. DES performs somewhat better than URANS. Drag and lift are not predicted as well. The time-averaged and instantaneous flow fields are studied, particularly to determine their suitability for the physics-based prediction of noise. The two time-averaged flow fields are similar, though the DES shows more turbulence intensity overall. The instantaneous flow fields are very dissimilar. DES develops a much wider range of unsteady scales of motion and appears promising for noise prediction, up to some frequency limit.


2018 ◽  
Vol 45 (1) ◽  
pp. 53-65 ◽  
Author(s):  
Jelena Svorcan ◽  
Ognjen Pekovic ◽  
Toni Ivanov

Although much employed, wind energy systems still present an open, contemporary topic of many research studies. Special attention is given to precise aerodynamic modeling performed in the beginning since overall wind turbine performances directly depend on blade aerodynamic performances. Several models different in complexity and computational requirements are still widely used. Most common numerical approaches include: i) momentum balance models, ii) potential flow methods and iii) full computational fluid dynamics solutions. Short explanations, reviews and comparison of the existing computational concepts are presented in the paper. Simpler models are described and implemented while numerous numerical investigations of isolated horizontal-axis wind turbine rotor consisting of three blades have also been performed in ANSYS FLUENT 16.2. Flow field is modeled by Reynolds Averaged Navier-Stokes (RANS) equations closed by two different turbulence models. Results including global parameters such as thrust and power coefficients as well as local distributions along the blade obtained by different models are compared to available experimental data. Presented results include fluid flow visualizations in the form of velocity contours, sectional pressure distributions and values of power and thrust force coefficients for a range of operational regimes. Although obtained numerical results vary in accuracy, all presented numerical settings seem to slightly under- or over-estimate the global wind turbine parameters (power and thrust force coefficients). Turbulence can greatly affect the wind turbine aerodynamics and should be modeled with care.


Author(s):  
M. Sergio Campobasso ◽  
Mohammad H. Baba-Ahmadi

This paper presents the numerical models underlying the implementation of a novel harmonic balance compressible Navier-Stokes solver with low-speed preconditioning for wind turbine unsteady aerodynamics. The numerical integration of the harmonic balance equations is based on a multigrid iteration, and, for the first time, a numerical instability associated with the use of such an explicit approach in this context is discussed and resolved. The harmonic balance solver with low-speed preconditioning is well suited for the analyses of several unsteady periodic low-speed flows, such as those encountered in horizontal axis wind turbines. The computational performance and the accuracy of the technology being developed are assessed by computing the flow field past two sections of a wind turbine blade in yawed wind with both the time- and frequency-domain solvers. Results highlight that the harmonic balance solver can compute these periodic flows more than 10 times faster than its time-domain counterpart, and with an accuracy comparable to that of the time-domain solver.


2012 ◽  
Vol 134 (6) ◽  
Author(s):  
M. Sergio Campobasso ◽  
Mohammad H. Baba-Ahmadi

This paper presents the numerical models underlying the implementation of a novel harmonic balance compressible Navier-Stokes solver with low-speed preconditioning for wind turbine unsteady aerodynamics. The numerical integration of the harmonic balance equations is based on a multigrid iteration, and, for the first time, a numerical instability associated with the use of such an explicit approach in this context is discussed and resolved. The harmonic balance solver with low-speed preconditioning is well suited for the analyses of several unsteady periodic low-speed flows, such as those encountered in horizontal axis wind turbines. The computational performance and the accuracy of the technology being developed are assessed by computing the flow field past two sections of a wind turbine blade in yawed wind with both the time-and frequency-domain solvers. Results highlight that the harmonic balance solver can compute these periodic flows more than 10 times faster than its time-domain counterpart, and with an accuracy comparable to that of the time-domain solver.


Author(s):  
Nazia Binte Munir ◽  
Kyoungsoo Lee ◽  
Ziaul Huque ◽  
Raghava R. Kommalapati

The main purpose of the paper is to use Computational Fluid Dynamics (CFD) in 3-D analysis of aerodynamic forces of a Horizontal Axis Wind Turbine (HAWT) blade and compare the 3-D results with the 2-D experimental results. The National Renewable Energy Laboratory (NREL) Phase VI wind blade profile is used as a model for the analysis. The results are compared with the experimental data obtained by NREL at NASA Ames Research Center for the NREL Phase VI wind turbine blade. The aerodynamic forces are evaluated using 3-D Computational Fluid Dynamics (CFD) simulation. The commercial ANSYS CFX and parameterized 3-D CAD model of NREL Phase VI are used for the analysis. The Shear Stress Transport (SST) Gamma-Theta turbulence model and 0-degree yaw angle condition are adopted for CFD analysis. For the case study seven varying wind speeds (5 m/s, 7 m/s, 10 m/s, 13 m/s, 15 m/s, 20 m/s, 25 m/s) with constant blade rotational speed (72 rpm) are considered. To evaluate the 3-D aerodynamic effect sectional pressure coefficient (Cp) and integrated forces about primary axis such as normal, tangential, thrust and torque are evaluated for each of the seven wind speed cases and compared with the NREL experimental values. The numerical difference of values on wind blade surface between this study and 3-D results of NREL wind tunnel test are found negligible. The paper represents an important comparison between the 3-D lift & drag coefficient with the NREL 2-D experimental data. The results shows that though the current study is in good agreement with NREL 3-D experimental values there is large deviation between the NREL 2-D experimental data and current 3-D study which suggests that in case of 3-D analysis of aerodynamic force of blade surface it is better to use NREL 3-D values instead of 2-D experimental values.


2018 ◽  
Vol 42 (3) ◽  
pp. 243-251 ◽  
Author(s):  
Edison H Caicedo ◽  
Muhammad S Virk

This article describes a multiphase computational fluid dynamics–based numerical study of the aeroacoustics response of symmetric and asymmetric wind turbine blade profiles in both normal and icing conditions. Three different turbulence models (Reynolds-averaged Navier–Stokes, detached eddy simulation, and large eddy simulation) have been used to make a comparison of numerical results with the experimental data, where a good agreement is found between numerical and experimental results. Detached eddy simulation turbulence model is found suitable for this study. Later, an extended computational fluid dynamics–based aeroacoustics parametric study is carried out for both normal (clean) and iced airfoils, where the results indicate a significant change in sound levels for iced profiles as compared to clean.


Sign in / Sign up

Export Citation Format

Share Document