Validation of a Mid-Fidelity Numerical Approach for Wind Turbine Aerodynamics Characterization

: This work is aimed at investigating the capabilities and limits of the mid-ﬁdelity numerical solver DUST for the evaluation of wind turbines aerodynamic performance. In particular, this study was conducted by analysing the benchmarks NREL-5MW and Phase VI wind turbines, widely investigated in the literature via experimental and numerical activities. The work was started by simulating a simpler conﬁguration of the NREL-5MW turbine to progressively integrate complexities such as shaft tilt, cone effects and yawed inﬂow conditions, offering a detailed portrayal of their collective impact on turbine performance. A particular focus was then given to the evaluation of aerodynamic responses from the tower and nacelle, as well as aerodynamic behavior in yawed inﬂow condition, crucial for optimizing farm layouts. In the second phase, the work was focused on the NREL Phase VI turbine due to the availability of experimental data on this benchmark case. A comparison of DUST simulation results with both experimental data and high-ﬁdelity CFD tools shows the robustness and adaptability of this mid-ﬁdelity solver for these applications, thus opening a new scenario for the use of such mid-ﬁdelity tools for the preliminary design of novel wind turbine conﬁgurations or complex environments as wind farms, characterised by robust interactional aerodynamics.


Introduction
In recent years, the challenge of developing sustainable ways to produce energy has caused the amount of interest in wind energy to rapidly increase.This challenge has been accompanied by an evolution in the methods used to evaluate the aerodynamic performance of wind turbines, representing a key factor for the development of more effective solutions in this field.These methods vary widely in complexity and fidelity, each serving different needs within the industry.To evaluate the aerodynamic performance of wind turbines, previous research studies were mainly directed at experiments considering both wind tunnels [1] or wind field measurements.Although the aerodynamics of wind turbines are extremely complicated, accurate measurements are still essential in the field of wind energy research even if they require a substantial economical effort.Nevertheless, numerical methods have been shown to be a valid and cheaper alternative to accurately calculate wind turbine aerodynamic performance.Indeed, numerical methods offer a different degree of fidelity based on their mathematical formulation and the computational effort required.
Starting from lower order numerical methods widely used for the study of wind turbines, blade element momentum (BEM) theory [2] represents a two-dimensional method suitable to provide a fast estimation of aerodynamic forces on rotating blades.However, the BEM method does not consider span-wise flow features on the blade and cannot capture the fine details related to boundary layer modelling.Moreover, this method relies on empirical functions to reproduce dynamic stall effects or blade tip loss due to compressibility [3].
Thus, the BEM method, while computationally efficient, is limited in accuracy, particularly for the simulations of complex flow conditions characterised by interactional aerodynamics or when stall phenomena occurs.On the other hand, Computational Fluid Dynamics (CFD) methods based on solving Navier-Stokes equations provide more comprehensive and accurate results in wind turbines applications, thanks to an effective modelling of blade boundary layers and of a physical representation of their wakes [4].In particular, a great effort was spent on the use of high-fidelity commercial CFD software for wind turbine system applications, as shown for example by the works by Rezaeiha et al. [5,6].Moreover, several works in the literature use uncertainty quantification methods in order to evaluate error propagation in high-fidelity CFD numerical simulations, as shown for example in [7,8].The trade-off for the increasing accuracy is a significant increase in computational resources and time, making CFD less practical for early-stage design and more suited for final validation and detailed analysis of wind turbine applications.
Mid-fidelity numerical methods, widely investigated in recent years particularly for rotary wing applications, strike a balance between computational efficiency and accuracy.These methods, mainly based on the Vortex Particle Method (VPM) for wake modelling [9,10], showed to be promising to provide a more accurate representation of both aerodynamic loads and interactional flow physics of wind turbines at comparatively lower computational costs with respect to high-fidelity CFD.Thus, these tools can be considered quite interesting to be used for the preliminary design phase of wind turbine configurations or complex environments as wind farms, characterised by robust interactional aerodynamics.Indeed, several mid-fidelity aerodynamic codes were developed in recent years to support the design process of rotary wing applications.Just to cite few examples, DLR developed UPM [11], an unsteady panel and free-wake code originally intended for aeroacoustic simulations of helicopters but recently applied on arbitrary complex configurations as compound rotorcraft [12].Moreover, GENeral Unsteady Vortex Particle (GENUVP) software, based on a panel method coupled with a VPM solver, was developed at the National Technical University of Athens (NTUA) for both aerodynamic and aeroacoustic simulations of rotary wing applications, from helicopters to wind turbines [13].
Starting from a collaboration with Airbus, Politecnico di Milano developed the midfidelity aerodynamic solver DUST implementing a robust formulation of different aerodynamic elements from lifting lines to surface panels.This numerical solver, based on the use of VPM for wake modelling, was thoroughly validated against experiments and high-fidelity CFD over several rotorcraft configurations with increasing complexity [14][15][16][17] and now represents one of the most mature mid-fidelity solvers widely used both in academics and industry, particularly for the preliminary design of multi-rotors architectures in aeronautics.Thus, the present work aims to provide a thorough validation of the DUST solver for a different challenge, i.e., the aerodynamic characterisation of wind turbines.The validation of DUST within this work will consist of a meticulous comparison of aerodynamic load and response predictions with those from high-fidelity CFD codes and experimental standards.With this aim, two benchmark test cases widely investigated in wind energy field were selected.
The first benchmark is the NREL (National Renewable Energy Laboratory) offshore 5 MW baseline wind turbine (NREL-5 MW).It was developed as a theoretical framework for validating turbine models [18] and was extensively validated using the NREL BEM model FAST [19].Indeed, due to the fact that all details of the geometry and airfoil properties of the turbine are publicly available, this three-blade turbine, in an upwind, horizontal axis configuration, was widely investigated in the literature.URANS simulations of the NREL-5 MW were performed in [20][21][22].Then, particular interest was given to the study of wake prediction, as shown in [23][24][25].This wind turbine reference was also used as a good benchmark to validate Fluid Structure Interaction tool chains, as shown in [26][27][28][29][30].
The second benchmark used in this work is represented by NREL Unsteady Aerodynamics Experiment Phase VI ('NREL Phase VI'), performed on a 20 kW, 10 m diameter, twobladed rotor wind turbine in the NASA Ames Research Center wind tunnel [1].The main purpose of the tests was to acquire measurements on the aerodynamic and structural response of a horizontal axis wind turbine.Many researchers have performed simulations of the NREL Phase VI experiment for validation of their solvers.Sorensen et al. [31] and Langer Moller et al. [32] used a RANS simulation to predict loading and force coefficients, while in Hsu et al. [33], fluid-structure interaction simulations were performed with stabilized multiscale methods on this test case.
The specific goals of the present work were therefore to corroborate the accuracy of DUST against high-fidelity simulations for the NREL-5 MW wind turbine, with a particular emphasis on the impact of principal wind turbine parameters such as rotor shaft tilt, rotor cone angle, yawed inflow and blade-tower interaction.Moreover, DUST simulations were also performed for the NREL Phase VI small wind turbine to confirm the code's robustness by comparison with experimental data, focusing on the code's ability to accurately forecast performance metrics particularly focused on high yawed conditions.Generally, this research aims to establish the foundation for the broader application of mid-fidelity codes in wind energy, particularly during iterative design and analysis phases requiring a proper balance between computational swiftness and accuracy.

DUST Software
The mid-fidelity open-source software DUST www.dust.polimi.it(accessed on 1 March 2024) has been developed by Politecnico di Milano since 2017 for the simulation of the interactional aerodynamics of rotorcraft and unconventional aircraft configurations [34].The code is released under the open-source MIT license.The capabilities of the code have been extended in recent years and DUST has been also coupled to the open-source multibody solver MBDyn [35], also developed at Politecnico di Milano, enabling the performance of aeroelastic analysis of complete rotorcraft configurations.The mathematical formulation of DUST relies on an integral boundary element formulation of the aerodynamic problem and VPM of the wakes.This choice naturally fits the Helmholtz decomposition of the velocity field from a mathematical point of view and avoids the numerical instabilities occurring with connected models of the wake.A model can be composed of several components, connected to user-defined reference frames, whose position and motion can be defined in a hierarchical way.Different aerodynamic elements allow for different levels of fidelity in the model, ranging from lifting line elements to zero-thickness lifting surfaces and surface panels for thick solid bodies.
Since the VPM is a grid-free Lagrangian approach that models the free vorticity of wakes, a DUST model needs only surface meshes eliminating the need for a volume mesh of the flow surrounding the object being studied.This approach enables a robust representation of the interacting wakes generated by lifting surfaces and bodies, as is typically seen in rotary-wing vehicle applications.A DUST model can be made up of multiple components that use different aerodynamic elements, offering varying levels of accuracy in the model ranging from lifting line elements [36,37] to vortex lattices and surface panels [38].
The lifting line (LL) element used in DUST is a one-dimensional model that represents thin, slender lifting bodies.The element considers the sectional aerodynamic coefficients of lift, drag, and pitching moment calculated on the base on the local angle of attack, the local Reynolds number and the local Mach number.By taking into account the airfoil's camber and thickness, as well as the effects of viscosity and compressibility, the model is able to accurately calculate aerodynamic loads.However, since the problem is stated in explicit form, numerical instability may occur during aerodynamic or aeroelastic analysis.
The vortex lattice (VL) method offers a discrete representation of the average surface of thin lifting bodies by modeling them as a sheet of vortex rings with an intensity of Γ.This is equivalent to a piecewise-uniform surface doublet distribution.While the VL method accounts for compressibility effects by applying a Prandtl-Glauert correction to the loads, it cannot capture the non-linear behavior of aerodynamic loads.To ensure numerical stability, the VL element utilizes a fully implicit numerical scheme.The non-linear vortex lattice element (NL-VL) has been implemented in order combine a surface discretisation with the capability of LL to use an aerodynamic database to introduce viscous corrections.This hybrid element formulated according to an implicit method allows for a higher robustness in terms of accuracy and numerical stability.Further details about the formulations of this quite novel element can be found in [39].
The surface panels (SP) are formulated as a Morino-like problem [38] in implicit form.This allows the SP element to accurately describe the real shape of the component by applying non-penetration boundary conditions in the physical position.However, the SP element is not able to capture the non-linear effects on aerodynamic loads.
The reader is referred to [34] for a complete mathematical formulation of the code.

NREL-5 MW Configuration
The first test case considered is the NREL-5 MW offshore wind turbine.It is a threebladed horizontal axis utility scale turbine with a rotor diameter D = 126 m developed by Jonkman et al. [18].This is a widely used test case for aerodynamics and aeroelasticity research as all details of the turbine geometry and of the operating conditions are public, which is not the case for large-scale commercial turbines.The NREL-5 MW machine is a purely theoretical wind turbine and therefore no experimental data is available to validate simulations.However, the extensive simulations conducted on this test case lend substantial credibility to its applicability for this study.A summary of the turbine's essential features is presented in Table 1 and illustrated in Figure 1.features is presented in Table 1 and illustrated in Figure 1.The blade includes two non-lifting cylindrical segments and a combination of five DU airfoils and one NACA airfoil distributed over 18 spanwise sections.Specifications regarding the airfoil positioning, chord lengths and angles of twist used to build the DUST numerical model are detailed in Table 2.The hub of the NREL-5 MW wind turbine is positioned 5 m upwind from the tower's centerline and elevated 90 m above the ground in an undeflected state.This configuration aligns with the vertical distance from the tower top to the hub height, set at 2.4 m.Consequently, the yaw bearing's elevation above ground or mean sea level (MSL) is established at 87.6 m.
The wind turbine configuration and parameters used in DUST simulations for the NREL-5 MW wind turbine are listed in Table 3.A total of four rotor configurations are explored, each addressing different aspects of rotor and wind turbine performance.Configurations C1, C3, and C4 are assessed under direct rated inflow conditions, whereas configuration C2 is evaluated with a yaw inflow angle up to 30 • .Finally, to investigate the blade tower interaction, a full wind turbine configuration was considered, dealing with C4 case rotor modelling.
Thus, the DUST numerical model of this wind turbine includes blades, tower, hub and nacelle.The tower and the nacelle surface are modelled using surface panels (SP), specifically 2700 elements for the tower and 600 elements for the nacelle.Two different models of the blades are built using both lifting lines (LL) and non-linear vortex lattices (NL-VL).Both aerodynamic elements require input tables containing the aerodynamic coefficients of two-dimensional airfoils to introduce the viscous corrections on the load obtained with the potential solution.The airloads tables are obtained by XFoil [40] simulations performed before the stall region, i.e., for a range of angle of attack between −15 • and 15 • .Airfoil aerodynamic coefficients in the angle of attack range between −180 • and 180 • are obtained through the Viterna method [41].The selection of the time-step discretisation of the simulations over the rotor revolution and of the spatial spanwise discretisation of the blades was dictated by the will to balance computational efficiency with the fidelity of load pattern capture, thus serving the dual aims of computational economy and detailed aerodynamic characterization.Indeed, one of the main goals of the activity is to evaluate DUST proficiency in capturing the fundamental aerodynamic phenomena pertinent to wind turbine simulations and to ensure that the tool remains computationally economical.Indeed, the scope of the present work is to show the capability of the tool to be used for the preliminary design of novel wind turbine concepts or multiple rotors configurations as wind farms.As a matter of fact, the requirement of this work was to perform all the simulations using a commercial laptop.With this aim, the validity of the choice of parameters was verified by comparison with high-fidelity CFD results by Dose [26] for the C1 configuration. Figure 2 shows DUST simulations results of the NREL-5 MW wind turbine obtained with a discretisation of 75 elements along the span-wise direction and 72 time-steps per revolution for a total duration of the simulation of 50 s, corresponding to about 10 rotor revolutions.In particular, DUST simulations results are reported for both LL and NL-VL elements with the same span-wise discretisation but different chord-wise elements for the latter method.A convergence of integral loads can be observed for DUST simulations after 30 s corresponding to almost eight revolutions.A good correlation in terms of aerodynamic power between CFD and DUST results, considering both LL and NL-VL elements, is obtained.A reasonably good agreement was found in terms of aerodynamic thrust, particularly considering NL-VL elements, while LL elements seem to slightly under-predict the CFD thrust output.In particular, the sensitivity analysis to chord-wise discretisation shown in Figure 2 leads to the choice of five elements for all the simulations performed on the NREL-5 MW wind turbine in this work.Generally, the good accuracy observed with respect to CFD for the calculation of integral loads confirms the suitability of the time and spatial discretisation selected for DUST simulations of this wind turbine.
Figure 3 shows the mesh of the complete DUST numerical model of the NREL-5 MW wind turbine.

NREL Phase VI Configuration
The second test case considered is the NREL Unsteady Aerodynamics Experiment (UAE) Phase VI wind turbine [1] tested at NASA Ames wind tunnel.A summary of the turbine's essential features is presented in Table 4.
The rotor blades use the NREL S809 airfoil and the relative specifications regarding the chord and twist distributions are detailed in Table 5.The hub of the NREL Phase VI wind turbine is positioned 1.5 m upwind from the tower's centerline and elevated 11.5 m above the ground in an undeflected state.The base and tip diameter of the tower are 0.61 m and 0.41 m, respectively.The model of the NREL Phase VI wind turbine is shown in Figure 4.The DUST numerical model includes the rotor and the tower-nacelle assembly in order to be consistent for the comparison with experimental results.The tower and the nacelle are modelled by surface panels using 2700 and 600 elements, respectively.For the present test case, the DUST numerical models of the blades are built using NL-VL elements only.Following the same approach considered for the previous wind turbine, the selection of the number of time-steps and span-wise elements discretisation was validated by comparison with experimental results, now available for this wind turbine.DUST simulation results obtained by discretising each blade with 50 elements along the span-wise direction and varying the chord-wise elements distribution are compared in Figure 5 with experimental thrust and torque measured in [1] and CFD data [32].The representation of wind turbine loads evaluated by DUST, plotted as function of the number of revolutions, highlights the bumps related to the passage of two blades in front of the pylon and shows a converged behaviour after almost four rotor revolutions.The DUST numerical model presents reasonable accuracy with experimental data for thrust and slightly underestimates the torque, showing a weak sensitivity with respect to chord-wise element discretisation.Thus, as performed for the previous wind turbine, a discretisation of five chord elements was considered as a good trade-off and therefore was used for all the simulations performed on this wind turbine in this work.
The wind turbine configuration and parameters used in DUST simulations for the NREL Phase VI wind turbine are listed in Table 6.In particular, the simulations consider the full model with blade pitch at 3 • windward and different yaw angles.Table 6.Simulations parameters of the NREL Phase VI wind turbine.

NREL-5 MW Wind Turbine
The first step of the validation process involves simulating the simplest configuration of the wind turbine C1 at rated conditions, not taking into account neither the shaft tilt angle α T nor blade cone angle α C .DUST results are compared to the ones obtained with high-fidelity numerical simulation conducted by Dose et al. [26] using the RANS approach and by Leng et al. [28].
Figure 6 shows the comparison of the normal and tangential sectional forces acting along the normalised spanwise coordinate of the blade computed by the DUST NL-VL model and high-fidelity CFD simulation by Dose et al. [26].The normal force density comparison illustrated in Figure 6a exhibits a close alignment between DUST and highfidelity CFD data along almost the majority of the blade span.However, near the blade tip, DUST slightly underestimated the normal force.This behaviour is expected as wing tips are characterized by a more complex flow regime where three-dimensional phenomena and tip vortices generation typically leads to significant viscous effects and separated flow regions, which cannot be accurately captured by the DUST mid-fidelity approach.Tangential force density comparison depicted in Figure 6b shows an underestimation of DUST computation before the first 30% of blade span and a slight overestimation in the central spanwise region.The high-fidelity CFD model, however, shows a more pronounced decrease in tangential force towards the blade tip with respect to the DUST results.Generally speaking, the behaviour of the sectional force comparison confirms the capability of DUST to competently reproduce the general force distribution along the blade span.The performance of a wind turbine is significantly affected by the wind turbine wake.Therefore, a crucial step for DUST validation is the ability to correctly capture the wake development and its main features in terms of velocity field.Then, the deficit of the axial velocity in the wake calculated by DUST is compared with the high-fidelity simulation by [42], considering four streamwise positions located at different distances x downstream from the rotor plane.The comparison of velocity profiles is shown in Figure 7 as function of z, representing the distance from the shaft axis parallel to the rotor plane.
A reasonably good agreement was found between axial velocity profiles in the wake of the wind turbine, particularly near to the rotor disk, i.e., x = 0.125D and x = 0.25D.Some discrepancies were found at higher distances where CFD results show a slightly larger extension of the wake, probably due to numerical diffusion.Generally, Figure 8 illustrates the wake structure computed by DUST shown as an iso-surface of the vortex Q-criterion.Readers are referred to results reported in [42] to appreciate the capability of DUST to capture the behaviour of the wake structure computed by CFD.
The influence of blade cone angle α C and shaft tilt angle α T on the aerodynamic performance under rated inflow conditions is now explored by comparing results obtained with DUST and high-fidelity CFD by Dose et al. [26] over C1, C3 and C4 cases.Generally, the sensitivity of integral loads to cone and tilt angles variations is minimal.Indeed, by adding a blade cone angle, a slight increase in the rotor power is observed.On the other hand, introducing a shaft tilt slightly reduces rotor power and thrust.In particular, the integral power and thrust on the rotor are compared in Table 7.Both DUST models provide values quite similar to high-fidelity CFD findings.A difference below 2% is found for power evaluation with both LL and NL-VL elements, while thrust evaluation obtained by DUST with NL-VL is in quite good agreement with CFD.An underestimation of rotor thrust of about 8% is found for the LL DUST model with respect to CFD.Case C4 enabled us to provide an insight into the effects of varying freestream velocity, thanks to high-fidelity CFD results from Sørensen and Johansen [20], Lin et al. [21], Wu and Nguyen [22].Figure 9 compares the computed power and thrust as function of wind speeds.Generally, both DUST models capture the behaviour of power and thrust curves evaluated by CFD, thus confirming the capability to correctly model the effects of wind inflow on rotor performance.In particular, the LL element model tends to slightly overestimate power production below the rated wind speed, while the NL-VL element model slightly overestimates rotor power above the rated speed.Now considering thrust evaluation, the LL element model tends to underestimate CFD high-fidelity data along the entire velocity spectrum, while a comparatively better agreement is obtained with the NL-VL element model.Different yawed configurations are then investigated.At first, a yaw angle γ = 30 • is considered to assess the integral and the sectional parameters with respect to highfidelity CFD simulation by Dose et al. [26].Then, other two yaw angles configurations are considered, namely γ = 10 • , and 20 • , respectively, to be compared with the LES study by yaw [24].Table 8 compares the integral rotor power and thrust computed for C2 configuration with γ = 30 • where errors in terms of power and thrust with respect to CFD by Dose et al. [26] are below 10% for the LL model and below 5% for the NL-VL model.Compared to the axial flow condition case, power and thrust outputs decrease in accordance to CFD evaluation by Dose et al. [26].Results by Dose et al. [26] are based on the blade resolving CFD method, so dynamic stall effects occurring particularly in yawed conditions are implicitly included.On the other hand, no dynamic stall models are considered in DUST formulation.In order to quantify the local effect on blades subjected to yawed inflow, sectional normal and tangential forces are compared at three different spanwise coordinates in Figure 10 as a function of blade azimuthal angle ψ in the last computed revolution.DUST tends to calculate higher normal and tangential force densities up to 60% of the span coordinate, while a slightly lower normal force density is found near the blade tip compared to the high-fidelity CFD results.Nevertheless, generally the global behaviour of the sectional forces acting on blades are captured by DUST representation also in yawed condition.Another important feature to better understand the capability to capture the wind turbine wake behaviour for yawed configuration is the study of the evolution of the wake center-lines and the associated skew angles.The wake center-lines are imaginary lines that run through the middle of the wake, marking the trajectory of the region with the highest velocity deficit caused by the wind turbine and representing the path that the wake follows downstream of the turbine.This line can be calculated as the integral of the velocity deficit on a wake section perpendicular to the flow direction divided by the area of that section [43].This method gives an average value that accounts for the different velocities across the entire wake area, providing a singular measure of the wake's strength.The skew angles, instead, give a more quantitative analysis of the wake center trajectories.Skew angle χ is the angular measure between the position vector p extending from the rotor disk midpoint and the normal vector to the rotor disc, as shown in Figure 11 and defined in Equation (1). Figure 12 shows the skew angles behaviour calculated by DUST compared with LES results by Kim and Lee [25], for yawed conditions ranging from γ = 0 • to γ = 30 • .The comparison shows how DUST is capable to correctly evaluate the wake center line evolution and the corresponding skew angles that change with respect to yaw angle and show a peak around x/D = 2.A more general evaluation of DUST capabilities to capture wind turbine wake is provided by the visualization of the flow fields.Figures 13 and 14 show the contours of the averaged and instantaneous streamwise velocity computed by DUST (NL-VL) at different yaw angles.This flow field representation highlights the capability of DUST to evaluate the evolution and deflection of the wakes occurring at yawed conditions, as well as the highly unsteady behaviour of the flow fields in these conditions.If compared with LES results by Kim and Lee [25], the global behaviour of the flow field is captured but the speed of deformation of the sectional shape of the wakes computed by DUST is higher.Indeed, the coherence of the wake shape is lost at a lower distance from the rotor disk with respect to LES computation by Kim and Lee [25], as shown by the contours of the averaged streamwise velocity on transversal planes shown in Figure 15.The final configuration considered deals with the complete case equipped with tower and nacelle for C4 rotor condition.Figure 16 illustrates a visualisation of the wake structure computed by DUST for the complete wind turbine shown as an iso-surface of the vortex Q-criterion.As performed before, DUST simulation results using both LL and NL-VL models are compared with the high-fidelity CFD calculation by Yu and Kwon [30] and Lin et al. [21].The results in terms of integral loads comparison are presented in Table 9 for a single rotor and a full wind turbine.The introduction of nacelle and tower provides a reduction of the rotor performance for both power and thrust.This behaviour is captured by DUST simulations, thus showing the capability to reproduce interactional effects between rotor wake and tower, as performed by CFD.In particular, the NL-VL model shows integral values in quite good agreement with high-fidelity CFD [21] as evaluated differences are below than 5% for both power and thrust.

NREL Phase VI
The second part of the validation deals with the NREL Phase VI wind turbine and was focused in particular on the analysis of highly yawed conditions for 7 m/s freestream velocity.Firstly, an overview of rotor thrust and torque computed by DUST at zero yaw is given in Table 10 and compared with experiments by Hand et al. [1] and high-fidelity CFD data by Länger-Möller et al. [32].
Table 10.Rotor loads comparison for a full NREL Phase VI wind turbine at 7 m/s and zero yaw with high-fidelity CFD results from [32] and experimental results from [1].A good agreement of the integral load values are found for DUST simulations for the complete wind turbine model, as thrust is near to the mean experimental value, while torque is slightly underestimated similarly to CFD.A detailed insight on loads behaviour is provided by the time history of the single blade aerodynamic torque evaluated on the last computed revolution, compared with experiments and CFD [33] in Figure 17.Results distinctly illustrate the capability of DUST to reproduce time evolution of the blade torque observed by experiments and particularly the bump occurring as a blade passes in front of the tower around the azimuthal angle of 180 • .Now, the attention is focused on the evaluation of loads for yawed conditions.Thus, blade sectional normal and tangential force are compared in Figures 18 and 19, respectively, for different yaw angles.In particular, DUST results are compared with NREL measurements [1] and with numerical results obtained by RANS simulations [44] and by a non-linear vortex lattice method (NVLM) described in [45].
The trends of the normal sectional forces computed by DUST resemble the experimental data for the considered yawed angles, just showing a slight downward offset with respect to the mean experimental values at γ = 60 • , similarly to that found by both the RANS and NVLM approach used as a comparison.Concerning tangential force density, a reasonably good agreement with experiments and NVLM approach results [44] is obtained particularly at γ = 0 • .Some discrepancies are observed at higher yaw angles in the inner blade region, where also the RANS approach shows larger error with respect to experiments.On the other hand, a quite good agreement with both experimental and numerical data is also found for higher yaw angles in the outer half region of the blade.Generally, the sectional load values obtained by DUST resemble an accuracy comparable with the similar NVLM method used in [45] and somewhat higher than the one obtained by high-fidelity CFD, thus further confirming the promising capability of the DUST midfidelity approach for the evaluation of wind turbine aerodynamic performance also under deeply yawed conditions.[33] and experimental results from [1] at V ∞ = 7 m s −1 .The standard deviation of the experimental load is reported as a banded region.
. Averaged normal force coefficient densities for different yaw angles compared with the NREL measurement data [1] and numerical results using RANS [44] and non-linear vortex lattice method (NVLM) [45].
. Averaged tangential force coefficient densities for different yaw angles compared with the NREL measurement data [1] and numerical results using RANS [44] and non-linear vortex lattice method (NVLM) [45].
Finally, a representation of the wind turbine wakes under the different considered yawed conditions is shown in Figure 20 by means of iso-surfaces of Q-criterion coloured by non-dimensional axial velocity component contours.To impose the yaw rotation, the entire turbine was rotated as performed during the wind tunnel tests, assuming the was flow always directed from left to right.This comparison shows the capability of DUST to evaluate the coherent vortical structure of the wind turbine wake as well as their deflection related to the different yawed conditions.related to the different yawed conditions.

Conclusions
The present activity was aimed to validate the mid-fidelity numerical solver DUST for the evaluation of wind turbine aerodynamics.This study involved the numerical modelling of two wind turbines widely investigated in the literature for numerical tool validation purposes.A thorough comparison of DUST simulations results with both experimental and high-fidelity CFD data showed the limits and capabilities of the mid-fidelity solver to reproduce loads and flow field characteristics in different operating conditions including high yaw angles.
In particular, the study of the NREL-5 MW wind turbine indicates that DUST is capable to provide a reasonably good degree of accuracy with respect to high-fidelity CFD for the computation on integral loads, showing errors always below 10%.In particular, the nonlinear vortex lattice method implemented in DUST showed higher fidelity with respect to the lifting line method, particularly at higher inflow speeds.A similar accuracy in terms of integral loads is found in yawed conditions even if the local comparison of force density along blade span indicates some appreciable discrepancies with respect to CFD, particularly for the tangential component.Concerning the reproduction of the wind turbine wake, DUST shows a global good evaluation of wake evolution and trajectory also under yawed conditions.Nevertheless, quantitative comparison shows a better representation of the near field, while higher discrepancies are found further from the rotor disk where deformation of the sectional shape of the wake is higher with respect to LES representation, particularly under yawed conditions.
Thanks to the availability of experimental results, the study of the NREL Phase VI wind turbine provides a reasonably strong evaluation of the capabilities of DUST, particularly focused on the analysis of loads under yawed conditions.The integral loads computed with DUST are in reasonably good agreement with both experimental and CFD data.A more detailed insight on the blade force densities along the span showed a general agreement between DUST and numerical approaches characterised by different degrees of fidelity, particularly for the normal force component.Higher discrepancies with experiments are found for the tangential components, particularly on the inner region of the blade.
Generally, the outcomes of the present work encourage the use of thee DUST midfidelity approach for the preliminary aerodynamic design of wind turbine configurations requiring a huge number of simulations that could be performed with a somewhat lower computational cost with respect to high-fidelity CFD.In particular, the capability of this VPM-based solver to reproduce interactional flow physics mechanisms of rotary wing machines opens a new scenario for the use of mid-fidelity solvers in the design phase of wind farms.

Figure 1 .
Figure 1.Main geometrical quantities of the NREL-5 MW MW wind turbine.

Figure 2 .
Figure 2. Comparison of the aerodynamic power and thrust computed by DUST and CFD by Dose et al. [26] for the NREL-5 MW wind turbine in C1 configuration; V ∞ = 11.4 m s −1 .

Figure 3 .
Figure 3. Mesh of the NREL-5 MW wind turbine built for DUST simulations (blades modelled with NLVL with n chord = 5).

Figure 7 .
Figure 7.Comparison of the axial velocity distributions at different downstream locations computed by DUST and CFD by Fan et al. [42] for the NREL-5 MW wind turbine in C1 configuration; V ∞ = 11.4 m s −1 .

Figure 8 .
Figure 8. Instantaneous iso-surface for the vortex criterion Q coloured by the velocity magnitude calculated by DUST for the NREL-5 MW wind turbine in C1 configuration; V ∞ = 11.4 m s −1 .

Figure 9 .
Figure 9.Comparison of integral rotor power (a) and thrust (b) computed by DUST and high-fidelity CFD by Sørensen and Johansen[20], Lin et al.[21], Wu and Nguyen[22] as function of inflow wind speeds for the NREL-5 MW wind turbine in C4 configuration.

Figure 16 .
Figure 16.Instantaneous flow field computed by DUST for the full turbine configuration: iso-surface of Q-criterion colored by adimensional free-stream velocity.

Figure 17 .
Figure 17.Comparison of the time history of the single blade aerodynamic torque over a full revolution compared with high-fidelity results from[33] and experimental results from[1] at V ∞ = 7 m s −1 .The standard deviation of the experimental load is reported as a banded region.

Figure 20 .
Figure 20.Comparison of the iso-surfaces of Q-criterion coloured by non-dimensional axial velocity component contours for the NREL Phase VI wind turbine at 7 m/s at different yaw.

Table 1 .
Summary of the NREL-5 MW MW wind turbine properties.

Table 2 .
Summary of the main properties of the NREL-5 MW MW wind turbine blade.

Table 3 .
Simulation parameters of the NREL-5 MW wind turbine.

Table 5 .
Summary of the main properties of the NREL Phase VI wind turbine blade.

Table 9 .
[21]arison of integral rotor power and thrust computed by DUST and high-fidelity CFD at V ∞ = 11.4 m s −1 by Lin et al.[21]for configuration C4.