Next Article in Journal
Study of a Bi-Vertical Axis Turbines Farm Using the Actuator Cylinder Method
Next Article in Special Issue
Structural Modeling and Failure Assessment of Spar-Type Substructure for 5 MW Floating Offshore Wind Turbine under Extreme Conditions in the East Sea
Previous Article in Journal
Energy Efficiency and Decarbonization in the Context of Macroeconomic Stabilization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Accuracy of uRANS and LES-Based CFD Modeling Approaches for Rotor and Wake Aerodynamics of the (New) MEXICO Wind Turbine Rotor Phase-III

by
Shantanu Purohit
,
Ijaz Fazil Syed Ahmed Kabir
and
E. Y. K. Ng
*
School of Mechanical and Aerospace Engineering, Nanyang Technological University, 50 Nanyang Avenue, Block N3, Singapore 639798, Singapore
*
Author to whom correspondence should be addressed.
Energies 2021, 14(16), 5198; https://doi.org/10.3390/en14165198
Submission received: 7 June 2021 / Revised: 12 August 2021 / Accepted: 17 August 2021 / Published: 23 August 2021
(This article belongs to the Special Issue Advancement in Wind Turbine Technology)

Abstract

:
This work presents a comparison study of the CFD modeling with two different turbulence modeling approaches viz. unsteady RANS and LES, on a full-scale model of the (New) MEXICO rotor wind turbine. The main emphasis of the paper is on the rotor and wake aerodynamics. Simulations are carried out for the three wind speeds considered in the MEXICO experiment (10, 15, and 24 ms−1). The results of uRANS and LES are compared against the (New) MEXICO experimental measurements of pressure distributions, axial, radial, and azimuth traverse of three velocity components. The near wake characteristics and vorticity are also analyzed. The pressure distribution results show that the LES can predict the onset of flow separation more accurately than uRANS when the turbine operates in the stall condition. The LES can compute the flow structures in wake significantly better than the uRANS for the stall condition of the blade. For the design condition, the mean absolute error in axial and radial velocity components along radial traverse is less than 10% for both the modeling approaches, whereas tangential component error is less than 2% from the LES approach. The results also reveal that wake recovers faster in the uRANS approach, requiring further research of the far wake region using both CFD modeling approaches.

1. Introduction

Wind energy is becoming an important source of renewable energy. Due to low power density, wind turbines are clustered in wind farms to maximize energy output. In fact, an important limitation of such arrangement is the wake shed by the upstream turbines on the downstream turbines, which leads to lower power production and increased loads on the downstream wind turbines [1,2]. A wind turbine wake has flow scales ranging from millimeters to several hundred meters. Therefore, it necessitates studying the dynamics of wind turbine wakes to design next-gen wind turbines and wind farms.
To design efficient and reliable wind turbines, it is vital to have a proper understanding of rotor and wake aerodynamics. The two most popular experimental studies conducted to examine these effects are NREL Phase VI wind turbine experiment [3] and MEXICO rotor experiments [4,5]. Several authors have used these studies to validate their computational models. Accordingly, a multitude of numerical studies has been done on the MEXICO rotor. Bechmann et al. [6] employed EllipSys3D CFD code, a RANS-based flow solver with k-ω SST turbulence model, to validate the MEXICO rotor measurements and extracted 3-D airfoil coefficients from the numerical simulations. The same flow solver has been applied to study the near wake of the MEXICO rotor [7] and validation of the Phase-II MEXICO results [8]. The actuator line method integrated into EllipSys3D is also used to study the near wake of the MEXICO rotor [9,10]. Nathan et al. [11] compared the results of the Phase-II MEXICO rotor with actuator line modeling using OpenFOAM and EllipSys3D and found that both of these flow solvers predict lower speed wake measurements well but fail when the flow becomes separated. Other researchers have used the MEXICO rotor to explain the rotational effects at the inboard region which leads to stall delay [12,13]. Carrión et al. [14] proposed an all-Mach scheme for flow around wind turbines and used the MEXICO rotor for validation. The wake velocity results obtained from this scheme proved to be consistent with the experimental measurements. Recently, Li et al. [15] used the MEXICO rotor and employed a Lattice Boltzmann/ LES method for the investigation of wake characteristics. The evolution of velocity and wake characteristics agreed well with the experiments.
Due to the computational cost associated with the LES, it is commonly used in indirect rotor modeling with the actuator line and actuator disk methods [16]. Several other CFD studies have been done on the MEXICO rotor [17,18,19,20,21,22,23,24]. The RANS-based flow modeling approaches have been widely studied in the literature. Regodeseves et al. [17] conducted unsteady RANS simulations of the MEXICO rotor including the effects of nacelle and tower and found good agreement of numerical computations against the experiments. More recently, Zhang et al. [18] compared RANS and DDES simulations for the stall condition of the MEXICO rotor and found improved load prediction in the outboard region of the blade compared to those obtained from the RANS. Tsalicoglou et al. [22] studied uniform and yawed inflow conditions and Rodrigues and Lengsfeld [23] investigated the axial conditions of the MEXICO rotor using the RANS approach. Qian et al. [24] compared the RANS and DES numerical computation approaches for the yaw condition of the MEXICO rotor. Johansen et al. [25] used DES to study parked conditions and sinusoidal pitching of the NREL Phase VI wind turbine about its axis. Amiri et al. [19] assessed the performance of different Reynolds Stress Transport (RST) models to predict wind turbine aerodynamics at the stalled condition of the NREL Phase VI wind turbine. The computational results suggest that out of several RST models, the quadratic pressure-strain and elliptic bending RST models perform better for the entire wind speed spectrum considered for the analysis. Hsu et al. [20] used FEM to study the aerodynamics of the NREL Phase VI wind turbine using ALE-VMS formulation and applied the same discrete formulation for all the wind speed cases. The proposed formulation predicted the flow solution accurately, thus, circumventing the need to adjust the RANS-based turbulence model for different flow regimes. Kabir and Ng [26] applied the uRANS method and modeled the rotating effect with a sliding mesh approach to improve the BEM method.
Hybrid RANS and LES techniques bridge the gap between computationally expensive LES and less accurate RANS modeling. It solves the attached boundary layers using the RANS approach and separated flows using the LES approach. This is extensively used in modeling flow through the horizontal axis wind turbines (HAWTs) [27], but little or no improvement over unsteady RANS is observed for the prediction of blade characteristics [19,25,28,29,30].
As seen from the above literature review, a host of numerical studies on HAWTs are based on the RANS approach. To the best of the authors’ knowledge, LES is not used widely in full rotor modeling. Mo et al. [21] used large-eddy simulations to study the wake characteristics of the NREL Phase VI wind turbine. It is noted that no study employing LES has been conducted for the MEXICO rotor that gives a detailed comparison of both rotor and wake aerodynamics. This paper critically assesses the aerodynamic performance of the (New) MEXICO rotor using uRANS and LES. It is a well-known fact that the spatio-temporal wake structure of the turbulence can be better captured by the LES method [31]. This is investigated here for the (New) MEXICO rotor. The pressure distribution is compared at the five radial locations of r R = 0.25, 0.35, 0.60, 0.82, and 0.92. The results are further complemented with the wake velocity comparison downstream of the blade. Here, axial, radial, and azimuthal traverses of velocity components are also compared against the experimental measurements. A discussion on wake velocity and vorticity contours from the two approaches is also presented.

2. Numerical Framework

This section describes the CFD analysis setup of a wind turbine for predicting its aerodynamic performance. This section begins with the description of the MEXICO experiment and the computational geometry used for the analysis, followed by the mesh details and boundary conditions applied to the domain. The governing equations used in uRANS based k-ω SST model and WALE-based LES are highlighted next. Finally, a short discussion on the number of revolutions required in LES for the results to reach a steady-state is included.

2.1. MEXICO Experiment

The Mexnext Phase-III Model Experiments In Controlled Conditions (MEXICO) project was an international collaboration led by the Energy Research Centre of the Netherlands (ECN). The goal of the project was to generate big datasets that can be used to improve and develop the new aerodynamic models of wind turbines [3]. In this study, results of the (New) MEXICO experiment—which was made available in 2018—are used, which is a follow-up of Phase 1 and Phase 2 MEXICO experiments and resolves discrepancies of the results obtained for loads and velocities in Phase-1. The diameter of the MEXICO rotor blade is 4.5 m and consists of 3 airfoils viz. the DU91-W2-250 airfoil, RISØ-A2-21 airfoil, and NACA 64-418 airfoil from 20 to 45.6% span, 54.4% to 65.6% span, and from 74.4% till the tip of the blade, respectively, as shown in Figure 1. The MEXICO rotor is chosen here for the comparison of uRANS and LES turbulence modeling because along with the rotor measurements, it also includes detailed measurement data of the flow-field, unlike the NREL Phase VI experiment which only includes rotor measurements. The pressure measurement is obtained using 148 Kulite pressure sensors, distributed over 5 radial locations of r R = 0.25 ,   0.35 ,   0.60 ,   0.82 ,   0.92 , and velocity components are measured using stereo PIV measurement technique. The dotted lines in Figure 1 indicate the locations of pressure measurements. The blade is twisted and tapered throughout its length, and the global pitch angle is kept at −2.3°. The blade is rotated at 425.1 rpm, and the experiments were performed in a wind tunnel for the wind speeds of 10, 15, and 24 ms−1 which corresponds to turbulent wake state, design condition, and stall condition of the blade, respectively. In this study, the flow is perpendicular to the plane of rotation, and the effect of blade yaw is not considered.

2.2. Computational Domain, Boundary Conditions, and Solution Methodology

The computational domain used for the CFD analysis is divided into 3 domains, as shown in Figure 2 and Figure 3. The blade is enclosed in a small cylindrical domain (called the rotating domain) which is given a rotation of 425.1 rpm. A fine mesh domain is created in the region adjacent to the rotating domain, called stationary domain (fine), to capture near wake characteristics more accurately. These two domains are enclosed in a stationary domain, as depicted in Figure 2. The rotating domain is placed at the center of the stationary domain. The other dimensions of the computational domain are shown in Figure 2. The effect of the tower is not considered in the analysis. To model the relative motion between the rotating domain and the stationary domain, a sliding mesh approach is used. This approach, however, is more computationally demanding than the other commonly used technique, called the Multiple Reference Frame (MRF) approach, but is more accurate and describes the transient phenomenon accurately [27].
The inlet is given a uniform velocity and 0 Pa gauge pressure is specified as pressure outlet boundary condition. The far-field surfaces are modeled with a symmetrical boundary condition. Turbulence is generally not high in the wind tunnels and the experiment does not provide any turbulence intensity value, keeping this in mind, a low value of 0.1% is chosen as turbulence intensity value and viscosity ratio ( μ t μ ) is set as 10. Residual for continuity was kept as 1 × 10−3 and all other variables were kept as 1 × 10−5.
The unsteady RANS simulations were initially run at higher time steps to allow the flow to reach the outlet. The converged solution was then run at a low value of time-step corresponding to 1° rotation of blade per time-step is used (3.92 × 10−4 s) [17]. Running LES from scratch is very computationally demanding, therefore, to save computational time, the LES uses the initial converged solution of uRANS at a larger time-step. The RANS model is then changed to the WALE-LES model and a lower time-step value (3.92 × 10−4 s) is used for 4 revolutions of the blade.
The ANSYS Fluent software (version 20.1) is used to perform the uRANS and LES simulations. Double-precision-based segregated solvers were applied. Pressure-velocity coupling is obtained with the SIMPLE algorithm. For better accuracy, the momentum and turbulence equations are spatially discretized using the second-order upwind scheme in uRANS. The momentum equations in the LES are discretized using a bounded central differencing approach. The Least Squares Cell-Based method is used to compute the gradients in both approaches. A first-order temporal discretization is adopted for uRANS, whereas a second-order bounded implicit temporal discretization is used for the LES approach [21].

2.3. Mesh Details

A good quality unstructured mesh composed of tetrahedral and prism elements was generated using Ansys ICEM-CFD software. The mesh of the blade with a zoomed view of the mesh at the tip and root region is shown in Figure 4. A sufficiently fine mesh was created on the blade with 15 inflation layers with the initial height of 1 × 10−5 m which gives y+ < 1, to resolve the boundary layer effects. The growth rate of 1.2 is used to extrude the inflation layer at the blade surface. Figure 5 shows zoomed view of the inflation layer mesh at the root region of the blade. The rotating and stationary domains were meshed independently and connected via mesh interfaces in Ansys FLUENT. The mesh consists of a total of 38.7 million elements, of which the rotating domain consisted of nearly 19 million elements and the rest of the elements were in the stationary domains.

2.4. Governing Equations

The flow fields in the large part of wind turbines are usually incompressible. Therefore, continuity and momentum equations can be written, respectively, as follows:
x i ( u i ) = 0
ρ t ( u i ) + x j ( ρ u i u j + p δ i j τ i j ) = ρ f i
In the above equations, ρ is the air density, u i and u j are the ith and jth direction velocity component (ms−1), p is the pressure (Pa), f i is the body force per unit mass (N/kg), δ i j is the Kronecker delta, and τ i j are the stresses acting on the fluid particle (N/m2). Here, i, and j take values of 1, 2, 3 corresponding to x, y, and z-direction, respectively.
Turbulence models in CFD are methods for modeling the evolution of turbulence in fluid flows. Turbulent flows are commonplace in nature, and they occur in the majority of engineering applications; therefore, a turbulence model is required for the majority of simulations. Turbulence models of various sorts, such as RANS, DES, and LES, all have strengths and limitations owing to the nature of the modeling assumptions. The turbulence modeling chosen for this study is uRANS and LES, which are based on transient analysis and are widely considered to be more accurate than steady analysis. There are several types of closure models for RANS, with k-ω SST being the most commonly used for wind turbine CFD study. Similarly, among the various LES models, the Wall-Adaptive Local Eddy-Viscosity (WALE) model was chosen for comparison. The equations involved in both of these closure models are described next.

2.4.1. uRANS k-ω SST Model

In this study, the closure model for uRANS equations is chosen as the k-ω SST turbulence model. It is a two-equation eddy-viscosity model that makes use of k-ω formulation in the boundary layer and switches to k-ε in the free-stream. Among several eddy-viscosity turbulence models, the previous studies for flow through HAWTs [7,8,12,17,22,23,32,33] have shown that the k-ω SST model performs better for the flow simulations through a wind turbine. This model is especially reliable in adverse pressure gradients and separating flows—the two common phenomena encountered in wind turbines. The equations associated with the k-ω SST turbulence model can be described as follows:
t ( ρ k ) + x i ( ρ k u i ) = x j ( Γ k k x j ) + G ˜ k Y k + S k
t ( ρ ω ) + x i ( ρ ω u i ) = x j ( Γ ω ω x j ) + G ω Y ω + D ω + S ω
where k is the turbulence kinetic energy (m2/s2); ω is the specific dissipation rate (s−1), Y k and Y ω are the dissipation rate of k and ω , respectively, Γ k and Γ ω are the effective diffusivity of k and ω , respectively; S k and S ω are the user-defined source terms of k and ω , respectively, D ω is the cross-diffusion term, G ω is the specific dissipation rate generation.
The reader is requested to refer to References [34,35] for the detailed description of the model and the equations involved.

2.4.2. WALE LES Model

The turbulent flows are characterized by both large and small eddies spanning a wide spectrum of time and length scales. The main principle behind LES is that it resolves larger eddies and models smaller eddies using Sub Grid-Scale (SGS) modeling. It does so by using a spatial and temporal filtering operation on the Navier-Stokes Equations (1) and (2). Filter width and grid size are the two important parameters that influence the filtering process. If the scale of the eddies is smaller than either of these two parameters, then such eddies get filtered out. In this way, it lies between RANS and DNS approaches in terms of the fraction of the resolved eddies.
The filtered Navier-Stokes equations can be written as:
ρ t + x i ( ρ u i ¯ ) = 0
t ( ρ u i ¯ ) + x j ( ρ u i ¯ u j ¯ ) = x j ( σ i j ) p ¯ x i τ i j x j
The stress tensor due to molecular viscosity, σij, can be defined as:
σ i j   [ μ ( u i ¯ x j + u j ¯ x i ) ] 2 3 μ u l ¯ x l δ i j
And the subgrid-scale stress, τij, is defined as:
τ i j = ρ u i u j ¯ ρ u i ¯ u j ¯
The above subgrid-scale stresses brought by the application of filtering operation are not known a priori and therefore must be modeled. The subgrid-scale turbulence stress can be obtained from:
τ i j 1 3 τ k k δ i j = 2 μ t S i j ¯
The WALE model is chosen to obtain eddy-viscosity, µt, because it does not require any wall or damping functions and its algebraic nature performs simulations faster.
μ t = ρ L s 2 ( S i j d S i j d ) 3 / 2 ( S ¯ i j S ¯ i j ) 5 / 2 + ( S i j d S i j d ) 5 / 4
where, S i j and S ¯ i j are square of velocity gradient tensor and shear stress tensor, respectively. The WALE models, Ls and S i j d , are defined, respectively, as:
L s = m i n ( κ d ,   C w V 1 / 3 )
S i j d = 1 2 ( g ¯ i j 2 + g ¯ j i 2 ) 1 3 δ i j g ¯ k k 2
In the above equations, κ is the von Karman constant, d is the distance to the closest wall (m), V is the computational cell volume (m3), g ¯ is the velocity gradient tensor The default value of WALE constant, C w , is 0.325.

2.5. Number of Rotor Revolutions for Steady-State Solution in LES

It is a well-known fact that to achieve stability in the LES solutions, it needs to be run for sufficiently longer flow times than the RANS simulations. Therefore, a sensitivity study is presented here to test the number of blade revolutions needed to achieve steady-state solutions. Figure 6 shows a comparison of axial, radial, and tangential velocity plots with experiments, at different revolutions of the turbine blade, for the wind speed of 24 ms−1. The mean absolute deviation for all three components of velocity is less than 1% for the three different revolutions of the blade considered in the analysis. Keeping the computational cost associated with LES in mind, the simulations here are, thus, run for four revolutions of the turbine blade.

3. Results and Discussions

3.1. Comparison of Pressure Distribution at 5 Radial Locations

The pressure coefficient (Cp) distributions at the 5 radial locations for the wind speeds of 10, 15, and 24 ms−1 which corresponds to the Tip speed ratio (λ) = Ω R U = 10, 6.67, and 4.17, respectively are shown in Figure 7, Figure 8 and Figure 9. The Cp equation can be given as:
C p = P P 0.5 ρ ( U 2 + ( ω r ) 2 )
The Cp values predicted from the two modeling approaches, both on the suction and pressure surfaces, overlap each other for the tip speed ratio (λ) of 10, as shown in Figure 7. The Cp distributions obtained from the CFD computations are in poor agreement in the inboard region of the blade, especially at the radial location of 0.25R, due to insufficient resolution of the pressure sensors to resolve the flow physics [3]. Similar Cp distribution at the radial location of 0.25R was also obtained in previous works [6,8,13,17,22,32,33]. The agreement between CFD predictions and experimental measurements at the mid-span and outboard region is, however, excellent, as evident from Figure 7. The Cp distribution predicted for the wind speed case of 15 ms−1 by both models, as presented in Figure 8, is in good agreement with the experimental values, with slight differences on the suction side. The LES model underpredicts pressure in the leading-edge region at the suction side in the inboard and mid-span region. On the other head, both models give nearly the same pressure prediction at the outboard region of the blade when the wind turbine operates in the design condition ( U = 15 ms−1, λ = 6.67). The differences in pressure distribution obtained from the two modeling approaches can be seen clearly for the wind speed of 24 ms−1, as evident from Figure 9. At the wind speed of 24 ms−1 (λ = 4.167), the turbine operates in stalled flow conditions with massive flow separation. At such high-speed flows, enhancement in suction peak is observed and the slope of pressure gradient reduces. The separation is said to occur when the pressure gradient slope ( d P d x = 0 ) does not change on the blade surface. This separation point is better predicted by the LES model as compared to the uRANS model, which predicts delay in the flow separation, as can be seen from Figure 9. The LES model overpredicts pressure distribution at the leading-edge region located at outboard stations (at spans of 0.82R and 0.92R). Overall, pressure distributions obtained from both models at different blade spans match well with the experimental measurements.

3.2. Velocity Distributions

3.2.1. Axial Traverse

Figure 10 depicts the axial traverse conventions used in plotting the three components of the velocity. The axial, radial and tangential components of velocities are plotted at a distance of z = −1.5 m when the rotor is in the position shown in Figure 10a. The axial traverse covers the measurement range from x = −4.5 m (upstream) to x = 5.9 m (downstream). The corresponding velocity components plots for the wind speed of 15 and 24 ms−1 are included in Figure 11. The vertical line at x/D =0 represents the rotor plane.
The three components of velocity i.e., axial, radial, tangential, give information about velocity deficit, wake expansion, and wake rotation, respectively. All the three components of velocity are predicted well by both modeling approaches, as can be seen in Figure 11a–f. The velocity components upstream of the blade are predicted accurately by both modeling approaches indicating the outer domain size chosen in the analysis is appropriate to exclude the tunnel effects. The oscillations in the radial and tangential velocity for the wind speed case of 15 ms−1 are not predicted accurately by both modeling approaches.
The differences in velocity components predicted by LES and uRANS can be seen much more clearly at the wind speed of 24 ms−1 (Figure 11b,d,f). The undulation in the axial velocity due to flow separation on the blades (Figure 11b) is satisfactorily predicted by the LES approach. Due to unsteady wake effects, the oscillations in the radial and tangential velocity components (Figure 11d,f) behind the rotor plane are better captured by the LES approach. The uRANS approach, on the other hand, averages out such oscillatory patterns.

3.2.2. Radial Traverse

This section compares the radial traverse of velocity components obtained using the LES and uRANS modeling approaches, which are plotted downstream of the turbine blade at a distance of 0.3 m from the rotor plane as shown in Figure 12. The velocity components are plotted when the blade orientation is as indicated in Figure 12a, along a solid, vertical line in Figure 12b.
The axial component of velocity for the wind speed of 10 and 15 m/s (Figure 13a,d), are predicted slightly closer by LES than uRANS against the experimental measurements. For the wind speed of 24 m/s (Figure 13g), the LES model follows the experimental trend of velocity deficit and recovery from 0.2 < r R < 0.4 but shows faster wake recovery from 0.4 < r R < 1, than seen in the experiments. The uRANS, on the other hand, is unable to capture the wake deficit and recovery in the inboard region (0.2 < r R < 0.4) but predicts much closer wake velocity at the mid-span and outboard region than the LES model, as evident from Figure 13g. The tip vortices, as can be seen from the sudden jump in axial velocity (Figure 13a,d) at r/R = 1, are predicted well by both modeling approaches for the wind speed of 10 and 15 ms−1. Whereas, such tip vortices are not predicted well when the rotor operates at the wind speed of 24 ms−1.
When the blade operates in stall condition (λ = 4.17), the separated flow in the mid-span and outboard region is accurately predicted by both of these modeling approaches (Figure 13h,i). The rotation in the wake (tangential component of velocity) is predicted satisfactorily by both of these modeling approaches as can be seen from Figure 13c,f,i. The vortex shed by the blade tip is predicted more accurately for the wind speeds of 10 and 15 ms−1 than for the wind speed of 24 ms−1. Overall, the predicted trends in the radial traverse of velocity components are closer to the experimental measurements.

3.2.3. Azimuthal Traverse

Figure 14 shows the conventions and definitions of azimuthal traverse of velocity components. Figure 15, Figure 16 and Figure 17 show a comparison of the azimuthal traverse of the three components of velocity (axial, radial, and tangential) with the experimental measurements for the λ of 10, 6.67, and 4.17, respectively. The figures are plotted at a distance of 0.3 m downstream from a rotor plane (Figure 14b) at five radial locations r R = 0.25, 0.35, 0.60, 0.82, and 0.92 (Figure 14a). Radial velocity can be defined as the component of wind velocity directed away from the rotor center, perpendicular to the axial direction, whereas tangential velocity is the wind component aligned along the rotor rotational direction [36]. Comparing the axial velocity plot (Figure 15a), it can be observed that the prediction is closer to the measurements at the inboard and mid-section of the blade. Slight deviation can be seen in the outboard region (r = 0.82R and 0.92R), where the prediction of the axial velocity component is closer to the measured value than the uRANS approach. Overall, the trend in axial velocity components of both the upwash and downwash of the blade is well captured by uRANS and LES. The axial velocity is, however, overpredicted by both uRANS and LES at the outboard region. The vortex shed by the inboard region (r = 0.25R and 0.35R) is not captured well in the downwash of the blade. The radial and the tangential components of velocity during upwash of the blade are not captured well by both of these approaches. Due to the separated blade wake and the root vortices in the inboard region for the lower wind speed case (10 ms−1), this region is difficult to predict. At the outboard region, radial component prediction obtained using uRANS is closer to the measurements, whereas LES overpredicts the radial velocity component. The tangential velocity predicted by both of these approaches nearly overlaps each other for all the radial positions considered (Figure 15c).
An excellent agreement between numerical and experimental measurements of the three velocity components of the azimuthal traverse is observed when the blade operates at the design condition (i.e., λ = 6.67), as can be seen in Figure 16. Compared to the turbulent wake state azimuthal plots, the tip and root vortices from the blade are better captured at this condition. The only appreciable difference in the axial velocity is observed at the outermost radial location of the blade (r = 0.92R), where the LES computations are closer to the experimental measurements than the uRANS predictions. The mean absolute error (MAE) in axial velocity computations and experimental measurements at r = 0.92R for all azimuthal positions is 18.40% and 12.04% from the uRANS and LES approaches, respectively. The MAE can be defined as:
M A E = | X C F D X e x p X e x p | × 100 %
Similarly, the radial velocity component (Figure 16b) is also approximated well by both modeling approaches with marginal deviations from the measurements. The tangential velocity computations from the LES approach are closer to the experimental measurements than the uRANS approach, with the average error being 18.5% and 10% from uRANS and LES, respectively. Moreover, as evident from Figure 16c, the tangential velocity during the downwash of the blade is better predicted by the LES approach than the uRANS approach.
Finally, a comparison of azimuthal traverse for the lowest tip speed ratio (λ = 4.17) corresponding to the wind speed of 24 ms−1 is shown in Figure 17. The axial velocity component (left side of Figure 17) in the inboard region is predicted well by the LES approach, whereas agreement is better between measurements and uRANS computations at the outboard region. A good agreement is also observed for the radial component of velocity as well, as evident from the plots in Figure 17b. The jump in the tangential velocity at 70° rotor azimuth due to induction effects, is captured well by the LES method. Overall, both approaches follow the general trend of the experimental observations.

3.3. Discussion on Near Wake Velocity at Several Downstream Locations

Insight into near wake is vital to improve wind farm efficiency as average wind farm losses due to wake could reach up to 20% of the total output of a large wind farm [37]. The near wake region can be defined as the downstream distance up to which blade influence can be experienced and after that distance, the wake profile is Gaussian in shape. Figure 18 represents velocity at various downstream locations obtained from uRANS and LES. It should be noted that wake radius increases with an increase in downstream distance. Therefore, to show the complete wake profile, r/R exceeds 1. Here, ‘r’ represents the local wake radius while ‘R’ is the blade radius. As mentioned, the overall wake radius might be greater than the rotor radius, therefore, we have included r/R greater than 1. The point where the velocity ratio (Uwake/U) reaches 1 is the free-stream region and the corresponding ‘r’ is the wake radius. The wake recovery for the wind speed of 10 ms−1 is found to be similar from both uRANS and LES approaches. The wake is recovered faster by the uRANS approach than LES for the wind speeds of 15 and 24 ms−1, as evident from Figure 18c,d. The curves are symmetric around the rotor axis because the effect of the tower is not considered in the analysis. Due to intense wake mixing at higher wind speeds, velocity deficit is lower at higher TSR, and vice-versa. The wake profile is Gaussian for all the tip-speed ratios after 3 rotor diameters downstream. Similar trends as shown in Figure 18d were also observed at farther downstream locations.

3.4. Discussion on Vorticity

Figure 19 shows the non-dimensional vorticity contours, given by the expression:
ζ n = ζ D u r e f
where ζ is the local vorticity magnitude (in s−1), D is the rotor diameter, and uref is the hub height velocity. From Figure 19, it can be observed that near wake vortices are resolved much better by LES than uRANS, as evident from the small vortical turbulence structures resolved by the LES approach. Such structures are averaged out in the uRANS approach.
Vortex is identified by using the λ 2 criterion developed by Jeong and Hussain [38]. To discern the presence of the vortex, they make use of pressure minimum definition by eliminating the effects of unsteady straining and viscous effects. Mathematically, pressure minimum created by vortex motion is defined as S 2 + Ω 2 , where S and Ω are the symmetric and anti-symmetric parts of the velocity gradient tensor. According to this definition, a vortex is defined when the eigenvalue, λ 2 , is less than zero. Figure 20 depicts the vortex shed by the blades, on a vertical plane, for the λ of 6.67 and 4.17. For the λ of 4.17, the tangle of vortices shed by the nacelle and the blade is captured better by LES than uRANS. From Figure 20, it can be seen that the vortex diffusion is reduced by applying large eddy simulation. Moreover, small scaler turbulent structures are also better resolved by the LES approach.

3.5. Comparison of the Computational Resources Required by uRANS and LES

The computational cost associated with the uRANS and LES modeling techniques is discussed in this section. The computational cost comparison is presented here for the λ = 6.67. Since the same higher time-step converged solutions were used for both uRANS and LES, the computational cost is compared when the blade completes four revolutions at low time-step (3.92 × 10−4 s) with both uRANS and LES techniques. Simulations were performed in a high-processing computing facility with 32 cores per host. The bar plot in Figure 21 compares the wall time (h) obtained using both the modeling techniques. Here, wall time means the number of hours it takes to complete four revolutions of the wind turbine blade. As expected, LES takes longer than uRANS as it requires higher computational resources. It can be commented that by expending computational cost, which is not drastically higher than the uRANS, better flow characteristics of near wake can be obtained using the LES approach.

4. Conclusions

A critical investigation is carried out to study the aerodynamic performance of the LES and uRANS flow modeling approaches. The Phase-III Mexnext (New) MEXICO rotor is chosen to carry out this numerical study, as it gives exhaustive information about the rotor measurements and flow field. The numerical results obtained using the two modeling methods are compared with the experimental results of pressure distributions, along with the flow variables downstream of the blade for axial, radial, and azimuthal traverses.
From this study, the following conclusions can be drawn:
  • The inception of flow separation for the stall condition of the blade is predicted accurately by the LES approach than by the uRANS approach;
  • When λ = 10, the azimuthal traverse of the radial and tangential components of velocity in the inboard region are not predicted more accurately by both of the modeling approaches;
  • The fluctuations in the wake when the blade operates at λ = 4.17 are properly captured by the LES approach;
  • The wake gets recovered faster in the uRANS approach than in the LES method;
  • LES can capture the small-scale vortices shed by the blade much more accurately.
The results presented in this research have given important insights into the near and far wake characteristics that will be useful in wind farm optimization. In our future works, presented results will be valuable in modeling wind turbines using actuator line modeling. Other important conditions such as different atmospheric boundary layer conditions and turbulent intensities will also be investigated in the future using machine learning techniques.

Author Contributions

S.P.: conceptualization, data curation, formal analysis, visualization, investigation, validation, software, writing—original draft. I.F.S.A.K.: data curation, conceptualization. E.Y.K.N.: supervision, writing—review & editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The data used have been supplied by the consortium which carried out the EU FP5 project Mexico: ‘Model rotor EXperiments In COntrolled conditions’. The consortium received additional support to perform the New Mexico measurements from the EU projects ESWIRP and INNWIND.EU.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

CpPressure coefficient
DRotor diameter (m)
f i Body force per unit mass (N/kg)
kTurbulence kinetic energy (m2/s2)
LsMixing length of sub-grid scales (m)
pPressure (Pa)
PPressure on the blade (Pa)
PAtmospheric pressure (Pa)
p ¯ Mean component of pressure term (Pa)
uVelocity (m/s)
u i ¯ Mean velocity along x-axis (m/s)
u j ¯ Mean velocity along y-axis (m/s)
urefFree-stream air velocity (m/s)
UaxAxial component of velocity (ms−1)
UradRadial component of velocity (ms−1)
UtanTangential component of velocity (ms−1)
VVolume (m3)
xiDistance from the origin along the positive x-axis (m)
yNon-dimensional distance normal to the wall
δ i j Kronecker delta
τ i j Stress tensor (N/m2)
ωSpecific dissipation rate (1/s)
εTurbulent dissipation rate (m2/s3)
μDynamic viscosity (Ns/m2)
ΨAzimuthal angle of the blade (°)
ρDensity (kg/m3)
ζnNon-dimensional vorticity
ζLocal vorticity (s−1)
λTip-speed ratio

Abbreviations

CFDComputational Fluid Dynamics
uRANSunsteady Reynolds-Averaged Navier Stokes
LESLarge Eddy Simulation
MEXICOModel Experiments In Controlled Conditions
NRELNational Renewable Energy Laboratory
SSTShear Stress Transport
DESDetached-Eddy Simulation
DDESDelayed Detached-Eddy Simulation
RSTReynolds Stress Transport
FEMFinite Element Method
ALE-VMSArbitrary Lagrangian-Eulerian Variational Multiscale
HAWTsHorizontal Axis Wind Turbines
WALEWall-Adaptive Local Eddy Viscosity
SIMPLESemi-Implicit Method for Pressure Linked Equations
MRFMultiple Rotating Frame
SGSSub-Grid Scale
DNSDirect Numerical Simulation
MAEMean Absolute Error
Subscripts
i, j, k1,2,3 or x, y, z-directions
axAxial
radradial
tantangential

References

  1. Kabir, I.F.S.A.; Ng, E.Y.K. Effect of different atmospheric boundary layers on the wake characteristics of NREL phase VI wind turbine. Renew. Energy 2019, 130, 1185–1197. [Google Scholar] [CrossRef]
  2. Kabir, I.F.S.A.; Safiyullah, F.; Ng, E.Y.K.; Tam, V.W.Y. New analytical wake models based on artificial intelligence and rivalling the benchmark full-rotor CFD predictions under both uniform and ABL inflows. Energy 2020, 193, 116761. [Google Scholar] [CrossRef]
  3. Hand, M.M.; Simms, D.A.; Fingersh, L.J.; Jager, D.W.; Cotrell, J.R.; Schreck, S.; Larwood, S.M. Unsteady Aerodynamics Experiment Phase VI: Wind Tunnel Test Configurations and Available Data Campaigns; Technical Report NREL/TP-500-29955; National Renewable Energy Laboratory: Golden, CO, USA, 2001. [Google Scholar]
  4. Schepers, G.J.; Boorsma, K.; Cho, T.; Gomez-Iradi, S.; Schaffarczyk, P.; Jeromin, A.; Shen, W.Z.; Lutz, T.; Meister, K.; Stoevesandt, B.; et al. Analysis of Mexico Wind Tunnel Measurements: Final Report of IEA Task 29, Mexnext (Phase 1); ECN-E, No. 12-004; Energy Research Center of the Netherlands (ECN): Petten, The Netherlands, 2012. [Google Scholar]
  5. Boorsma, K.; Schepers, G.J.; Gomez-Iradi, S.; Herraez, I.; Lutz, T.; Weihing, P.; Oggiano, L.; Pirrung, G.; Madsen, H.A.; Shen, W.Z.; et al. Final Report of IEA Wind Task 29 Mexnext (Phase 3); ECN-E-18-003; Energy Research Center of the Netherlands (ECN): Petten, The Netherlands, 2018. [Google Scholar]
  6. Bechmann, A.; Sørensen, N.N.; Zahle, F. CFD simulations of the MEXICO rotor. Wind Energy 2011, 14, 677–689. [Google Scholar] [CrossRef]
  7. Sørensen, N.N.; Bechmann, A.; Réthoré, P.-E.; Zahle, F. Near wake Reynolds-averaged Navier-Stokes predictions of the wake behind the MEXICO rotor in axial and yawed flow conditions. Wind Energy 2014, 17, 75–86. [Google Scholar] [CrossRef] [Green Version]
  8. Sørensen, N.N.; Zahle, F.; Boorsma, K.; Schepers, G. CFD computations of the second round of MEXICO rotor measurements. J. Phys. Conf. Ser. 2016, 753, 022054. [Google Scholar] [CrossRef]
  9. Shen, W.Z.; Zhu, W.J.; Sørensen, J.N. Actuator line/Navier-Stokes computations for the MEXICO rotor: Comparison with detailed measurements. Wind Energy 2012, 15, 811–825. [Google Scholar] [CrossRef]
  10. Nilsson, K.; Shen, W.Z.; Sørensen, J.N.; Breton, S.-P.; Ivanell, S. Validation of the actuator line method using near wake measurements of the MEXICO rotor. Wind Energy 2015, 18, 499–514. [Google Scholar] [CrossRef]
  11. Nathan, J.; Meyer Forsting, A.R.; Troldberg, N.; Masson, C. Comparison of OpenFOAM and EllipSys3D actuator line methods with (NEW) MEXICO results. J. Phys. Conf. Ser. 2017, 854, 12033. [Google Scholar] [CrossRef] [Green Version]
  12. Guntur, S.; Sørensen, N.N. A study on rotational augmentation using CFD analysis of flow in the inboard region of the MEXICO rotor blades. Wind Energy 2015, 18, 745–756. [Google Scholar] [CrossRef]
  13. Herráez, I.; Stoevesandt, B.; Peinke, J. Insight into Rotational Effects on a Wind Turbine Blade Using Navier–Stokes Computations. Energies 2014, 7, 6798–6822. [Google Scholar] [CrossRef] [Green Version]
  14. Carrión, M.; Steijl, R.; Woodgate, M.; Barakos, G.; Munduate, X.; Gomez-Iradi, S. Computational fluid dynamics analysis of the wake behind the MEXICO rotor in axial flow conditions. Wind Energy 2015, 18, 1023–1045. [Google Scholar] [CrossRef]
  15. Li, L.; Xu, C.; Shi, C.; Han, X.; Shen, W. Investigation of wake characteristics of the MEXICO wind turbine using lattice Boltzmann method. Wind Energy 2021, 24, 116–132. [Google Scholar] [CrossRef]
  16. Martínez-Tossas, L.A.; Churchfield, M.J.; Leonardi, S. Large eddy simulations of the flow past wind turbines: Actuator line and disk modeling. Wind Energy 2015, 18, 1047–1060. [Google Scholar] [CrossRef]
  17. Regodeseves, P.G.; Morros, C.S. Unsteady numerical investigation of the full geometry of a horizontal axis wind turbine: Flow through the rotor and wake. Energy 2020, 202, 117674. [Google Scholar] [CrossRef]
  18. Zhang, Y.; Deng, S.; Wang, X. RANS and DDES simulations of a horizontal-axis wind turbine under stalled flow condition using OpenFOAM. Energy 2019, 167, 1155–1163. [Google Scholar] [CrossRef]
  19. Amiri, M.M.; Shadman, M.; Estefen, S.F. URANS simulations of a horizontal axis wind turbine under stall condition using Reynolds stress turbulence models. Energy 2020, 213, 118766. [Google Scholar] [CrossRef]
  20. Hsu, M.-C.; Akkerman, I.; Bazilevs, Y. Finite element simulation of wind turbine aerodynamics: Validation study using NREL Phase VI experiment. Wind Energy 2014, 17, 461–481. [Google Scholar] [CrossRef]
  21. Mo, J.-O.; Choudhry, A.; Arjomandi, M.; Lee, Y.-H. Large eddy simulation of the wind turbine wake characteristics in the numerical wind tunnel model. J. Wind Eng. Ind. Aerodyn. 2013, 112, 11–24. [Google Scholar] [CrossRef]
  22. Tsalicoglou, C.; Jafari, S.; Chokani, N.; Abhari, R.S. RANS Computations of MEXICO Rotor in Uniform and Yawed Inflow. J. Eng. Gas Turbines Power 2014, 136, 011202. [Google Scholar] [CrossRef]
  23. Rodrigues, R.V.; Lengsfeld, C. Development of a Computational System to Improve Wind Farm Layout, Part 1: Model Validation and Near Wake Analysis. Energies 2019, 12, 940. [Google Scholar] [CrossRef] [Green Version]
  24. Qian, Y.; Zhang, Z.; Wang, T. Comparative Study of the Aerodynamic Performance of the New MEXICO Rotor under Yaw Conditions. Energies 2018, 11, 833. [Google Scholar] [CrossRef] [Green Version]
  25. Johansen, J.; Sørensen, N.N.; Michelsen, J.A.; Schreck, S. Detached-eddy simulation of flow around the NREL Phase VI blade. Wind Energy 2002, 5, 185–197. [Google Scholar] [CrossRef]
  26. Kabir, I.F.S.A.; Ng, E.Y.K. Insight into stall delay and computation of 3D sectional aerofoil characteristics of NREL phase VI wind turbine using inverse BEM and improvement in BEM analysis accounting for stall delay effect. Energy 2017, 120, 518–536. [Google Scholar] [CrossRef]
  27. Thé, J.; Yu, H. A critical review on the simulations of wind turbine aerodynamics focusing on hybrid RANS-LES methods. Energy 2017, 138, 257–289. [Google Scholar] [CrossRef]
  28. Li, Y.; Paik, K.-J.; Xing, T.; Carrica, P.M. Dynamic overset CFD simulations of wind turbine aerodynamics. Renew. Energy 2012, 37, 285–298. [Google Scholar] [CrossRef]
  29. Stone, C.; Lynch, C.; Smith, M. Hybrid RANS/LES simulations of a horizontal axis wind turbine. In Proceedings of the 48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Orlando, FL, USA, 4–7 January 2010. [Google Scholar]
  30. Lynch, C.; Smith, M. Unstructured overset incompressible computational fluid dynamics for unsteady wind turbine simulations. Wind Energy 2013, 16, 1033–1048. [Google Scholar] [CrossRef]
  31. Wu, Y.-T.; Porté-Agel, F. Atmospheric Turbulence Effects on Wind-Turbine Wakes: An LES Study. Energies 2012, 12, 5340–5362. [Google Scholar] [CrossRef]
  32. Plaza, B.; Bardera, R.; Visiedo, S. Comparison of BEM and CFD results for MEXICO rotor aerodynamics. J. Wind Eng. Ind. Aerodyn. 2015, 145, 115–122. [Google Scholar] [CrossRef]
  33. Lienard, C.; Boisard, R. Investigation of the MEXICO rotor aerodynamics in axial flow, including boundary layer transition effects. In Proceedings of the Wind Energy Symposium, AIAA SciTech Forum, Kissimmee, FL, USA, 8–12 January 2018. [Google Scholar]
  34. Menter, F.R. Zonal two equation k-ω turbulence models for aerodynamic flows. AIAA J. 1993, 93, 2906. [Google Scholar]
  35. Menter, F.R. Two-equation eddy-viscosity turbulence models for engineering applications. AIAA J. 1994, 32, 1598–1605. [Google Scholar] [CrossRef] [Green Version]
  36. Simley, E.; Angelou, N.; Mikkelsen, T.; Sjöholm, M.; Mann, J.; Pao, L.Y. Characterization of wind velocities in the upstream induction zone of a wind turbine using scanning continuous-wave lidars. J. Renew. Sustain. Energy 2016, 8, 013301. [Google Scholar] [CrossRef] [Green Version]
  37. Barthelmie, R.J.; Hansen, K.; Frandsen, S.T.; Rathmann, O.; Schepers, J.G.; Schlez, W.; Phillips, J.; Rados, K.; Zervos, A.; Politis, E.S.; et al. Modelling and measuring flow and wind turbine wakes in large wind farms offshore. Wind Energy 2009, 12, 431–444. [Google Scholar] [CrossRef]
  38. Jeong, J.; Hussain, F. On the identification of vortex. J. Fluid Mech. 1995, 285, 69–94. [Google Scholar] [CrossRef]
Figure 1. MEXICO rotor blade with airfoil distribution at different sections of the blade. The dotted vertical line represents the pressure measurement locations.
Figure 1. MEXICO rotor blade with airfoil distribution at different sections of the blade. The dotted vertical line represents the pressure measurement locations.
Energies 14 05198 g001
Figure 2. Entire computational domain used for the CFD analysis.
Figure 2. Entire computational domain used for the CFD analysis.
Energies 14 05198 g002
Figure 3. Zoomed view of the rotating domain and stationary domain (fine): (left) rotating domain inside a stationary domain of fine mesh (right) blade inside a rotating domain.
Figure 3. Zoomed view of the rotating domain and stationary domain (fine): (left) rotating domain inside a stationary domain of fine mesh (right) blade inside a rotating domain.
Energies 14 05198 g003
Figure 4. Mesh on the blade with a zoomed view of the mesh near the tip and hub region.
Figure 4. Mesh on the blade with a zoomed view of the mesh near the tip and hub region.
Energies 14 05198 g004
Figure 5. Zoomed view of the inflation layer in the root region of the blade.
Figure 5. Zoomed view of the inflation layer in the root region of the blade.
Energies 14 05198 g005
Figure 6. Investigation of the number of revolutions needed to achieve state solution in LES; (a) axial; (b) radial; (c) tangential velocity at z = −1.5 m for 24 ms−1 wind speed.
Figure 6. Investigation of the number of revolutions needed to achieve state solution in LES; (a) axial; (b) radial; (c) tangential velocity at z = −1.5 m for 24 ms−1 wind speed.
Energies 14 05198 g006
Figure 7. Pressure coefficient distribution at the 5 radial locations: r R = (a) 0.25; (b) 0.35; (c) 0.60; (d) 0.82, and (e) 0.92; for the freestream wind velocity (U) of 10 ms−1 (λ = 10).
Figure 7. Pressure coefficient distribution at the 5 radial locations: r R = (a) 0.25; (b) 0.35; (c) 0.60; (d) 0.82, and (e) 0.92; for the freestream wind velocity (U) of 10 ms−1 (λ = 10).
Energies 14 05198 g007aEnergies 14 05198 g007b
Figure 8. Pressure coefficient distribution at the 5 radial locations: r R = (a) 0.25; (b) 0.35; (c) 0.60; (d) 0.82, and (e) 0.92; for the freestream wind velocity (U) of 15 ms−1 (λ = 6.67).
Figure 8. Pressure coefficient distribution at the 5 radial locations: r R = (a) 0.25; (b) 0.35; (c) 0.60; (d) 0.82, and (e) 0.92; for the freestream wind velocity (U) of 15 ms−1 (λ = 6.67).
Energies 14 05198 g008aEnergies 14 05198 g008b
Figure 9. Pressure coefficient distribution at the 5 radial locations: r R = (a) 0.25; (b) 0.35; (c) 0.60; (d) 0.82, and (e) 0.92; for the freestream wind velocity (U) of 24 ms−1 (λ = 4.17).
Figure 9. Pressure coefficient distribution at the 5 radial locations: r R = (a) 0.25; (b) 0.35; (c) 0.60; (d) 0.82, and (e) 0.92; for the freestream wind velocity (U) of 24 ms−1 (λ = 4.17).
Energies 14 05198 g009aEnergies 14 05198 g009b
Figure 10. Axial traverse conventions and its definitions in the MEXICO rotor: (a) Front view; (b) Top view.
Figure 10. Axial traverse conventions and its definitions in the MEXICO rotor: (a) Front view; (b) Top view.
Energies 14 05198 g010
Figure 11. Axial traverse of 3 components of velocity at the location z = 1.5 m for the free-stream wind speed of 15 m/s and 24 m/s: (a,b) axial velocity; (c,d) radial velocity; (e,f) tangential velocity.
Figure 11. Axial traverse of 3 components of velocity at the location z = 1.5 m for the free-stream wind speed of 15 m/s and 24 m/s: (a,b) axial velocity; (c,d) radial velocity; (e,f) tangential velocity.
Energies 14 05198 g011
Figure 12. Radial traverse conventions and its definitions in the MEXICO rotor: (a) front view; (b) top view.
Figure 12. Radial traverse conventions and its definitions in the MEXICO rotor: (a) front view; (b) top view.
Energies 14 05198 g012
Figure 13. Radial traverse of 3 components of velocity at the location x = 0.3m for the free-stream wind speed of 10, 15, and 24 m/s: (a,d,g) axial velocity; (b,e,h) radial velocity; (c,f,i) tangential velocity.
Figure 13. Radial traverse of 3 components of velocity at the location x = 0.3m for the free-stream wind speed of 10, 15, and 24 m/s: (a,d,g) axial velocity; (b,e,h) radial velocity; (c,f,i) tangential velocity.
Energies 14 05198 g013aEnergies 14 05198 g013b
Figure 14. Azimuthal traverse conventions and its definitions in the MEXICO rotor.
Figure 14. Azimuthal traverse conventions and its definitions in the MEXICO rotor.
Energies 14 05198 g014
Figure 15. Azimuthal traverse of (a) axial; (b) radial; (c) tangential velocity for the wind speed of 10 ms−1. Notations: * Experimental, – – – uRANS, –––– LES.
Figure 15. Azimuthal traverse of (a) axial; (b) radial; (c) tangential velocity for the wind speed of 10 ms−1. Notations: * Experimental, – – – uRANS, –––– LES.
Energies 14 05198 g015
Figure 16. Azimuthal traverse of (a) axial; (b) radial; (c) tangential velocity for the wind speed of 15 ms−1. Notations: * Experimental, – – – uRANS, –––– LES.
Figure 16. Azimuthal traverse of (a) axial; (b) radial; (c) tangential velocity for the wind speed of 15 ms−1. Notations: * Experimental, – – – uRANS, –––– LES.
Energies 14 05198 g016
Figure 17. Azimuthal traverse of (a) axial; (b) radial; (c) tangential velocity for the wind speed of 15 ms−1. Notations: * Experimental, – – – uRANS, –––– LES.
Figure 17. Azimuthal traverse of (a) axial; (b) radial; (c) tangential velocity for the wind speed of 15 ms−1. Notations: * Experimental, – – – uRANS, –––– LES.
Energies 14 05198 g017
Figure 18. Wake velocities at various downstream locations of the wind turbine.
Figure 18. Wake velocities at various downstream locations of the wind turbine.
Energies 14 05198 g018
Figure 19. Vorticity contours for λ = 6.67 at two downstream locations of x = 1D (top) and x = 2D (bottom).
Figure 19. Vorticity contours for λ = 6.67 at two downstream locations of x = 1D (top) and x = 2D (bottom).
Energies 14 05198 g019
Figure 20. λ 2 criterion plot: (a) TSR = 6.67; (b) TSR = 4.17.
Figure 20. λ 2 criterion plot: (a) TSR = 6.67; (b) TSR = 4.17.
Energies 14 05198 g020
Figure 21. Wall time comparison of uRANS and LES for the wind speed case of 15 ms−1 (λ = 6.67).
Figure 21. Wall time comparison of uRANS and LES for the wind speed case of 15 ms−1 (λ = 6.67).
Energies 14 05198 g021
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Purohit, S.; Kabir, I.F.S.A.; Ng, E.Y.K. On the Accuracy of uRANS and LES-Based CFD Modeling Approaches for Rotor and Wake Aerodynamics of the (New) MEXICO Wind Turbine Rotor Phase-III. Energies 2021, 14, 5198. https://doi.org/10.3390/en14165198

AMA Style

Purohit S, Kabir IFSA, Ng EYK. On the Accuracy of uRANS and LES-Based CFD Modeling Approaches for Rotor and Wake Aerodynamics of the (New) MEXICO Wind Turbine Rotor Phase-III. Energies. 2021; 14(16):5198. https://doi.org/10.3390/en14165198

Chicago/Turabian Style

Purohit, Shantanu, Ijaz Fazil Syed Ahmed Kabir, and E. Y. K. Ng. 2021. "On the Accuracy of uRANS and LES-Based CFD Modeling Approaches for Rotor and Wake Aerodynamics of the (New) MEXICO Wind Turbine Rotor Phase-III" Energies 14, no. 16: 5198. https://doi.org/10.3390/en14165198

APA Style

Purohit, S., Kabir, I. F. S. A., & Ng, E. Y. K. (2021). On the Accuracy of uRANS and LES-Based CFD Modeling Approaches for Rotor and Wake Aerodynamics of the (New) MEXICO Wind Turbine Rotor Phase-III. Energies, 14(16), 5198. https://doi.org/10.3390/en14165198

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop