On the calibration of rotational augmentation models for wind turbine load estimation by means of CFD simulations

In this work the improved version of an engineering model which accounts for rotational augmentation effects by means of computational fluid dynamics (CFD) calibration is explored and discussed. Based on an analysis of the NREL Phase VI wind turbine, the novel modeling is presented, which uses as base line the formulation proposed by Chaviaropoulos and Hansen. The model is calibrated based on CFD simulations using OpenFOAM. The corresponding correction of the two dimensional polars is straightforward implemented within MoWiT, an in-house software for load calculation. The novel formulation results in improved lift and drag coefficients prediction in all considered cases, reducing the deviation with respect to the rotating CFD cases down to few percent. The optimal configuration including the correction for tip effects of Shen shows better agreements at the very tip of the blade. Furthermore the range of applicability for large wind turbine rotor blades based on a virtual 10MW\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10\,\hbox {MW}$$\end{document} rotor model is discussed.


Introduction
The Paris climate agreement in 2015 [1] aimed to reduce the greenhouse gas (GHG) and emissions in the biosphere, realizing the transition to a low-carbon economy. Wind energy sector has an important contribution to the fulfillment of these goals because of the targeted 28.5% of Europe's electricity demand by 2030 [2]. Added to the increasing demand for repowering of machines having reached the end of their lifetime [3], a number of challenges arose, among which enlarging the rotor blades of the wind turbines, while keeping the uncertainty about the related loads at a minimum. Robust tools such as aero-servo-elastic tools based on blade element momentum (BEM) theory are able to account for aerodynamic and structure interactions [4] with acceptable level of accuracy for most of the operating conditions. Within the last years, more sophisticated modeling based on high-fidelity aerodynamic approaches by the means of computational fluid dynamics (CFD) have emerged in parallel with a general development of aerodynamic modeling approaches [5]. The simplifica-B Elia Daniele elia.daniele@iwes.fraunhofer.de 1 Fraunhofer IWES, Am Seedeich 45, 27572 Bremerhaven, Germany 2 Flensburg University of Applied Sciences, Wind Energy Technology Institute, Nordstraße, 224943 Flensburg, Germany tions made in BEM limit the accuracy of the obtained results. Therefore semi-empirical engineering models are incorporated in it for improving its accuracy. The empirical nature of some of the available engineering models, which have been developed based on rather small-size experiments, puts in doubt their general validity [4]. Nevertheless, the improvements demonstrated by CFD have served as a base-line for some researchers to define their models [6]. Thus, the use of CFD as a bridge to improve modeling in wind energy, seems to be the key point for deriving more sophisticated models which can be implemented in BEM. In this work the improved version of an engineering model which accounts for rotational augmentation (RAM) effects by means of CFD calibration is explored and discussed. The goal is to employ CFD simulations for the axial inflow of the NREL Phase VI wind turbine [6] to calibrate an alternative rotational augmentation model. This leads to a reduction in BEM deviations and it has been made ready for the IWES in-house load simulation library "Modelica for Wind Turbines" (MoWiT). Three-dimensional and rotational effects are the origin of the rotational augmentation phenomenon which is mainly characterized by the increment of the lift coefficient as compared with the two-dimensional case (infinite airfoil) due to the delay in the flow separation at higher angles of attack. The lack of modeling of the rotational augmentation phenomenon causes over or under prediction of aerodynamic loads on wind turbines [7]. Therefore, there is an imperative need of research in this field of study for improving load and power production prediction. Over the last years many authors have described several approaches to explain the occurrence of rotational augmentation on blade turbines. For Madsen and Rasmussen [8], the combination of centrifugal and Coriolis forces plus radial pressure gradients are the cause of rotational augmentation. Lindenburg [9] assures that Coriolis effects are dominant over centrifugal forces, relieving the adverse pressure gradients on the boundary layer and consequently delaying the separation [10]. For other researches [11,12] the occurrence of the rotational augmentation is based on the blend of Coriolis, centrifugal forces and pressure gradients for the appearance of standing vortices on the blade suction surface. More recently, in Ref. [13], referring to large rotor size, it was observed that the root vortex system consist of counter-rotating trailing vortices which become stronger for the higher wind speed cases. In these same studies the ratio of the Coriolis to centrifugal forces weighted by the Rossby number was identified as capable of yielding information about the shear stress and boundary layer structure. In general, the effect of rotation is to increase the lift force in the separated flow region at inboard blade sections and decrease the lift force near the blade tip, but this part is still under research due to the inaccurate approach for tip corrections factors. However, the effect of rotation over the drag on the separated areas is still relative unknown. For instance, the work of Du and Selig [14] suggested that drag force decreases due to the separation moving towards the trailing edge, thus, shrinking the blade wake that reduces the sectional drag force on the rotating blade. Conversely, the work of other researchers, as Dowler et al. [10], indicated that the drag increases in a rotating system. The above mentioned different physical analysis led to alternative proposal for a semi-empirical modeling, to be later included within a load estimation tool. The investigations carried out by Snel [15] emphasized the analysis on boundary layer-equations for two-dimensional lift coefficient corrections, focusing on the ratio of local chord length to local radial location. Du and Selig [14] studied the integral threedimensional boundary layer including a tip speed ratio term and drag coefficient empirical correction. Chaviaropoulos and Hansen [16] developed a quasi-three-dimensional model based on the simplifications of the three-dimensional incompressible Navier-Stokes equations along the blade span-wise direction. The method described by Corrigan and Schillings [17] is characterized by the simplification of the boundary layer equations and was later evaluated by Tangler and Selig [18]. Bak et al. [19] proposed a model based on the analysis of the pressure distributions in a non-rotating and rotating blade. Dumitrescu et al. [20] introduced a correction model resulting from the simplification of the momentum integral equations (three-dimensional) for a general rotor blade. The model from Lindenburg [9] considered the centrifugal forces and the radial flow, the modified tip speed ratio and the local blade chord ratio as key parameter for the modification of the two-dimensional polars. Recently studies on large wind turbines blades show some interesting findings, as for instance, the work of Bangga et al. [21] indicating that the radial flow development at the inboard part of the blade creates a smaller size wake, inducing a reduction of the pressure level on the suction side, which in turn increases the lift coefficient. Until now, the mathematical "corrections" found in literature have been developed from experimental data, empirically or based on theoretical analyses. Some of these models correct both lift and drag, while others only the lift coefficient. The general formulation presented for the application of these corrections obeys to the following equations where C L represents the difference between C L 2D (airfoil lift) and the inviscid lift coefficient, C D is the difference between C D 2D (airfoil drag) and the drag coefficient at zero angle of attack, f L and f D are correction factors which depend on the model, and C L 3D and C D 3D are the corrected lift and drag coefficients for three-dimensional effects.
For the sake of brevity hereby a thorough description of all considered models, from Snel [15] to Lindenburg [9], from Du and Selig [14] to Corrigan and Schillings [17], is avoided. Instead the focus on the formulation proposed by Chaviaropoulos and Hansen [16] is reported. The latter represents a quasi-three-dimensional approach based on viscous-inviscid interactions applied for laminar and turbulent flows. The model considers that three dimensional and rotational effects are strongly related to geometrical features of the blade such as the local chord to radial position ratio c/r and the local twist angle θ of the blade. The maximum Mach number for this model is assumed to be 0.2 so then the flow can be studied as incompressible. The model has been derived from the quasi-three-dimensional incompressible Navier-Stokes equations in cylindrical coordinates taking the mean values and also making some assumptions for the mean values of the radial derivatives. A pressure correction algorithm is applied and the resulting pressure distribution around the airfoil is calculated, based on which lift and drag coefficients are deduced. Three-dimensional lift and drag coefficients are compared with the lift and drag from wind tunnel experiment, and finally three dimensional empirical relations have been deduced for the three dimensional correction. The corresponding correction factors for lift and drag coefficient are reported as where a ch = 2.2, h ch = 1 and n ch = 4 represent the constant, baseline parameters of this model [16]. According with this study, c/r and θ does not influence the aerodynamic performance in the attached flow regime, but they play a very important role at higher angles of attack where the flow is separated. In this case the Coriolis force sucks mass from the recirculation region and redirects it to the radial direction reducing its volume. This results in a pressure drop along the suction side, and an increases blade loading [16]. Clearly, deducing the value of the aerodynamic forces from both experiments or simulations required the estimation of the angle of attack (AoA). In a rotating blade, this aspect is directly connected with the determination of the local inflow velocity in the rotor plane which is free of bound circulation (local undisturbed velocity) [22]. There are several methods for calculating the angle of attack at the rotor, thus, the analyses of this work are restricted to three proposals: average azimuthal technique (AAT) [23], three-point (3P) [24] and zero gamma (ZG) [22]. Since there is no common agreement within the research community on which method should represent the most suitable option, for the remaining of this work the indication of the AoA from threedimensional CFD simulation refers always to an averaged value of the results obtained by applying the three aforementioned methods. The possible application and validity of the proposed calibrated model when extended to larger size turbines is discussed by means of a reference to literature results concerning the AVATAR 10 MW wind turbine virtual model based on the work described in Refs. [13,25,26]. The structure of the remaining part of the manuscript is as follows: a description of the CFD setup and its validation, together with the baseline modeling results compared to CFD, eevelopment of a model calibration based on the formulation of Chaviaropoulos and Hansen and its verification and validation, a concluding summary with outlooks. Figure 1 shows a flow-chart of the presented research. Standstill and rotating CFD simulations would be used to assess a hierarchy among the various modeling proposals available in the literature. By means of a comparison against selected experimental results the alternative model parameter are then derived and calibrated. Later this alternative is included into MoWiT.

Methodology
The NREL Phase VI model wind turbine was used for conducting the numerical simulations of this work, since a wide range of experimental data was available for this wind small scale turbine model, which made it the perfect candidate for validating different type of simulations, in particular standstill and rotating ones. Differently from many other experimental wind turbine models where as well more complex inflow scenarios have been analyzed, e.g. Ref. [27], in the case of the NREL Phase VI experiment detailed pressure measurements at several span-wise sections allowed for aerodynamic force integration and validation. In this section general information about the specification of the NREL Phase VI is presented. Later on, a compact description of the CFD library-OpenFOAM is given, along with the procedure and details regarding to the set up for the performed simulations.

Conducted simulations
The NREL Phase VI wind turbine and the S809 airfoil were simulated using OpenFOAM library version 4. 1 In order to reduce the number of simulations needed and allow for a broader verification of the results, both standstill and rotating cases were analyzed, in order to be able to compare multiple sections along the blade span-wise length by selecting different operating conditions for the two above mentioned cases. This circumstance lead to a fair agreement of the Reynolds number for the considered sections. For computing the twodimensional lift and drag sectional polars the blade has been pitched across a range of 78 • , both in stall and feather direction, with 2 • spacing. Firstly, steady simulations at several inflow velocities are performed with standstill configuration. Table 1 contains the Reynolds numbers (Re) along five spanwise positions corresponding to those for which experimental value are also available.
Standstill simulations provide information about the blade condition during the absence of rotation, although the finiteness effects at root and tip part of the blade does not allow a good agreement with two-dimensional airfoil simulations, available at Re equal to 0.65 × 10 6 , 0.75 × 10 6 , 0.85 × 10 6 and 0.95 × 10 6 . The simulation campaign for the rotating cases followed the same structure of that for the standstill configuration, except for the introduction of rotational effect by means of a source term in the momentum equation, as typically used in the advantageous case of single blade analysis. The related Reynolds number distribution along the spanwise direction is presented in Table 2.

Simulation setup
The OpenFOAM CFD library is used for fluid dynamic simulations, while for grid generation, pre-and post-processing were done by means of IWES in house software [28]. The turbulence model employed is the k-ω-SST [29] for its good stability in terms of convergence for the simulated cases. The initial values of boundary conditions are calculated based on what is reported in [30]. For all simulations the Reynolds-averaged Navier-Stokes equations (RANS) are  Table 1 Reynolds numbers along five span-wise sections on the blade for the standstill simulations   solved by means of the steady solver simpleFoam based on the SIMPLE algorithm [31], with a number of nonorthogonal correctors equal to 3. The convective terms are discretized with a second order scheme, while for the diffusive terms a fully orthogonal scheme was used. Regarding the linear solvers, for the pressure equation a generalized geometric-algebraic multi-grid solver (GAMG) was used.
For the equation of the velocity, the turbulent kinematic energy and its diffusion a smooth Gauss-Seidel solver was used. For the baseline two-dimensional simulations of the S809 airfoil a modified parameter set of the baseline k-ω-SST model was used to overcome the shortage in accuracy in predicting stall region. The HPC infrastructure of the University of Oldenburg [32] has been employed for all simulations.

Validation of the simulations
The CFD simulations were compared with the experimental results of NREL Phase VI [33] and S809 Airfoil [34]. In particular, for standstill cases, the sequence L3000ST0 is used, while for rotating case the sequence S0700000 is regarded. Although being RANS not the most suitable simulation approach for a configuration with extensive separated flow [35,36], it serves as a first step to measure the level of improvement a CFD driven parametrization may add to the baseline modeling. Another approximation resides in the choice of the turbulence model: Differently from the experiments, the presented simulation does not consider the presence of a laminar portion of the flow on blade and airfoil. This may impair the results for drag coefficient at low angle of attack as well as lift coefficient around the stall onset. In order to provide a qualitative assessment of the level of approxi-  Despite the possible validation for force and pressure coefficient, the proposal of an alternative model for rotational augmentation necessitates the capture of aerodynamic coefficient, thus, the AoA, or alternatively, the determination of the local velocity in the rotor plane without bound circulation [22], i.e. axial and angular induction factors. Several methods have been proposed for calculating the AoA, although studies of Guntur et al [37] reported the lack of a unique answer under different operating conditions. . Results correspond to the mean value of three methods selected from literature (AAT, 3P, ZG, respectively after [22][23][24]). In the same figure with the solid black line is indicated the value of the stall angle of attack for the S809 airfoil. For an inflow wind velocity value equal to 7 m s −1 are observed AoA larger than the two-dimensional stall value from root up to 50% of the span, indicating the need for correction. Figure 5b shows the comparison between two dimensional, stand still and rotating lift along the blade length, following the same schema derived from [21]. Results indicate the presence of augmentation effects at the inboard part of the blade, due to the increment of lift coefficients. Lift coefficient reduces as the radial position increases. The enhancement of C L is observed until r /R = 60%. The calculated values of two dimensional C L at middle span-wise positions have good agreement with the rotating case. However, standstill polar coefficients show lower values than those for the two-dimensional case, probably due to large finiteness effects at root affecting the solution up to the mid-span. In areas close to the tip, the presence

Results and discussion
In this section the focus lies on the upper row of the flowchart reported in Fig. 1, i.e. applying the standard RAM and compare them against the rotating CFD cases in order to determine the development of the novel model parametrization.

RAM verification
In this section, two-dimensional polars are corrected for tip and root effects by means of Prandtl [38] and Shen [39] formulations. Afterwards, the polars are corrected for rotational augmentation effects using the available RAM described in the introduction. A schema of this procedure is presented in Fig. 6.
The two-dimensional polars corrected by finiteness and rotational augmentation effects are compared to those resulting from CFD for the rotating configuration. At the root region of the blade, two span-wise sections are analyzed at  r /R = 33% and 39% for lift and drag. The two-dimensional polars are corrected by root effects using the formulation of Prandtl, whose modified proposal from Shen is used at tip. Subsequently, these finiteness-corrected polars are further manipulated to account for rotational augmentation effects using the following models: Du and Selig [14] (DS), Chaviaropoulos and Hansen [16] (CH), Snel [15] and Lindenburg [9] (SL). Considering the case with inflow wind velocity equal to 7 m s −1 and rotational speed of the rotor   Based on the results discussed quantitatively above non of the considered RAM alternatives is able to capture the features of the rotating polars along the entire blade span. Nevertheless, the CH model showed the smallest deviations among all other models.

Derivation of novel parametrization for RAM
Having identified a the good agreement between two-and three-dimensional polars by means of the adoption of the CH model in the previous section, a novel parametrization for RAM based on the CH formulation is proposed, aiming to extend its validity along the entire length of the blade, independently of the operating condition, being the latter for the NREL Phase VI related to the inflow wind velocity. , tip-CH-corrected two-dimensional case (··). Tip correction type is indicated in the sub-caption As it was described in the introduction, the CH model provides corrections for lift and drag coefficients based on the value of the variables c/r and φ, and the constant parameters a ch , h ch , n ch , see Eqs. (3) and (4).
The two-dimensional polars are corrected to fit to the CFD rotating polars by means of the baseline CH and nonlinear least squares method. While h ch and n ch are kept equal to the baseline literature values, an alternative value for a n is derived based on the minimization of the deviation between two and three-dimensional polars. Since RAM applies to rotating blade whose finiteness is not included in the two-dimensional polars, the same procedure is applied after correcting the two-dimensional polars for tip and root effects as well. Figure 10 shows an example of the fitting process for determining the optimal a n value. In this specific case the two-dimensional polars without tip and root corrections are reported, with h ch = 1, n ch = 4, and 15 m s −1 inflow velocity at a 33% span-wise location.
The value calculated for a n is 1.2 (Fig. 10a) and 2.2 ( Fig. 10b) for, respectively the case of lift and drag coef-  ficient, for which the calculated, optimal a n equals the value from the literature [16]. An analogous investigation is applied for all the conditions indicated in Table 2, resulting in a distribution of optimal values for a n depending on the inflow velocity and the spanwise position, as showed in Fig. 11.
Regarding the lift coefficient, a n values that minimize the deviation between two-dimensional and three-dimensional polars vary strongly with the span-wise position. Three different regions along the blade length are identified: 25-40%, 40-80%, and from 80% to the tip. After 10 m s −1 , i.e. as the turbine starts to be stall-regulated, the behavior of a n changes abruptly, and it is mainly due to the fact that the flow becomes extremely unstable and separated. Accordingly, after rated wind speed one single tendency is observed throughout the entire blade length.
The value of a n presents also a dependency on the inflow velocity value even when the flow is mostly non-separated along the blade, namely between 6 and 10 m s −1 . Within the inboard region from 25 to 40% r /R, no significant influence of the inflow velocity is noticed. Above 40% r /R the behavior of a n changes, reaching its inversion point towards the tip, as the inflow velocity increases. The last portion of the  blade also presents some dependency on inflow velocity as it is shown in Fig. 11a. It is more challenging to identify a pattern for the optimal values of a n for the drag coefficient. However, here also the a n values seem to depend on the span-wise position and the inflow velocity. Two regions are identified: from 25% until 70% and from 70% until tip. Regarding the inflow velocity there is a dependency of the calculation of a n when the flow is non-separated, and above the rated velocity, one tendency for a n emerges, see Fig. 11b.
In order to show the effects of tip and root corrections on the behavior of the optimal a n , the results of the fitting procedure based on the lift coefficient after the application of Prandtl and Shen models are reported, respectively, in Fig. 12a, b. The effect of Prandtl correction is more evident at the tip, where the optimal a n behaves completely different from what observed in Fig. 11a. The dependency of a n on the inflow velocity is kept similar as it was before the application of Prandtl model. The model of Shen provides variations at the tip region at both low and high inflow velocity values, leading to larger values for the optimal a n compared to those seen in Fig. 11a.

Alternative parameters' formulation
In this subsection is presented the formulation for an optimal a n for the lift coefficient obtained in conjunction with the application of the Shen tip correction model, that depends on the span-wise position and the inflow velocity value. Being the NREL Phase VI turbine model a stall-regulated machine, the distinction based on the operating condition is solely referred to the inflow velocity value. Based on the results of the previous section two types of parameters' formulations are proposed for the lift coefficient below and above an inflow velocity equal to 10 m s −1 . For low inflow velocities, Eqs. (5), (6) and (7) represent the alternative parameters' formulation of a n for the three distinct regions for the lift coefficient. The value of the fitting param- 3L and e 2L , e 3L are reported in Table 3. These are obtained by applying the fitting procedure that minimize the deviation of the polars using Eqs. (5), (6) and (7). Coefficients used for the optimal a n alternative parametric formulation for lift at low inflow velocities (6-10 m s −1 ), in conjunction with the application of the Shen correction model at the tip .78 e 3L = 91.91 Table 4 Coefficients used for the optimal a n alternative parametric formulation for lift at high inflow velocities (10-20 m s −1 ), in conjunction with the application of the Shen correction model at the tip The analysis performed for inflow velocities above 10 m s −1 led to Eq. (8). Table 4 shows the coefficients obtained by means of the fitting procedure that minimize the deviation of the polars using the above mentioned equation.

Assessment
The identified optimal configuration is assessed by comparing it against the results of the rotating CFD simulation cases used as reference. The two-dimensional polars are corrected by both the baseline and the alternative parametric CH formulation proposed by the authors. The accuracy of the alternative parametric formulation is investigated by comparing the optimal configurations with (models a n Prandtl and a n Shen ) and without (model a n ) tip and root corrections. Three sections along the blade are studied: r /R = 33%, 69% and 95%. Two inflow velocity values are considered, e.g. U ∞ = 7 m s −1 and U ∞ = 20 m s −1 , in order to assess the model performance at both low and high inflow velocities.
Regarding to lift coefficient, at r /R = 33% the alternative parametric model outperforms the original one, see Fig. 13).
For example, at 15 • AoA the baseline formulation of CH shows deviations of 15% and 40% at U ∞ = 7 m s −1 U ∞ = 20 m s −1 , respectively. Conversely, the alternative parametric formulation shows improvements, following closely the rotating polars with higher accuracy along the  (Fig. 13b), the model based on Prandtl correction for root effects further improves the results of about 5% below 15 • AoA, while for larger AoA the model with a n seems better. At lower inflow velocity values (Fig. 13a) the a n Prandtl and a n variants present similar results.
In the tip region, the formulation based on the Shen corrections for tip effect (model a n Shen ) is also presented. At low inflow speeds, the model a Prandtl exhibits excellent agreements with the rotating CFD polars (Fig. 14a). As the inflow velocity increases the a n Shen and a n variants show better results, and about 20% deviation is seen in both cases (Fig. 14b).
The alternative parametric formulation of the CH model results in improved lift and drag coefficients prediction in all cases. The optimal configuration including the correction for tip effects of Shen shows better agreements at the very tip of the blade. Slightly improvements were seen at the root when the configuration is accounting for the model of Prandtl.

Extension to large wind turbine rotor-blade scenarios
In this section the extension of the proposed modification is discussed based on the studies reported in Refs. [13,25,26] for the case of a 10 MW wind turbine rotor virtual model developed within the framework of the AVATAR EU project [40].  In Fig. 15 are reported the experimental and numerical results for a characteristics inboard section of the AVATAR rotor, where rotational augmentation effects are relevant. The related airfoil type is the DU-00-W-401, developed at TU Delft in The Netherlands, with the aim of having a reduced sensitivity to the roughness height for moderate to large thickness and large Reynolds number. Analogous airfoil families have been developed for low Reynolds number as well, e.g. in Ref. [41]. The impact of inflow turbulence and/or surface roughness is well represented by means of the CFD results available in same figure as computed in Ref. [25]. The high value for the airfoil maximum relative thickness introduce a pronounced difference in the lift coefficient polar between the free and forced transition scenarios. In particular, for the free transition case an early stall takes place between 2 and 4 degrees, while for the forced transition case a standard polar is observed.
The auxiliary input data for the application of the rotational augmentation model are taken from Ref. [43]. By means of a procedure similar to that described in the previous subsections, the a parameter for the CH RAM with the best fitting for the CFD unsteady RANS results can be determined. Depending on the choice for the input two-dimensional lift polar curve different scenarios can be investigated. In Fig. 16 the application of the fit procedure for the proposed alternative parametric formulation is shown. The fitting procedure results in a lam ≈ 0.7067 and a tur ≈ −1.3578 for, respectively, the laminar and turbulent input lift polar case. Moreover, a better approximation is possible when the laminar experimental lift polar curve is used as input for the fitted CH RAM. The effect of the Prandtl root correction for this  Fig. 16 Application of baseline and a-parameter based fitted CH RAM for the lift coefficient of the airfoil DU-00-W-401. Experimental results from Ref. [42] for a Reynolds number equal to 3 × 10 6 . CFD results from Ref. [25] for the same airfoil located at the 25% of the blade length for the AVATAR rotor. The CH RAM is applied to both laminar and fully turbulent experimental two-dimensional lift polar curve. The fitting procedure results in a lam ≈ 0.7067 and a tur ≈ −1.3578 for, respectively, the laminar and turbulent input lift polar case. Legend: experimental result with free transition (+), experimental result for forced transition at 2% suction side and 10% pressure side (•), baseline CH RAM applied on free transition experimental polar (-), baseline CH RAM applied on forced transition experimental polar (-), fitted CH RAM applied on free transition experimental polar (·), fitted CH RAM applied on forced transition experimental polar (·), CFD unsteady RANS from Ref. [25] for a comparable airfoil section having the maximum relative thickness equal to 40% and located at the 25% of the blade length for the AVATAR rotor (--) section, equal to 0.9768, does play a minor role in this case, leading to a value for the fitted a parameter equal to 0.7216.
Moreover, the correct capture of the lift polar curve identified by means of CFD simulation at the lowest inflow velocity speed value remains unattainable for the fitted CH RAM based on the laminar experimental lift polar curve. In fact, the lift polar point located at an AoA of about −2 • lies closer to the turbulent lift polar curve. For this reason, further threedimensional numerical simulation in absence of rotation may serve to determine the lift polar basis on which the RAM is to be applied, thus, separating the different levels related to the two-dimensional data referring to infinite body from the two-dimensional data referring to finite body, and both the latter from the two-dimensional data referring to finite rotating body.

Conclusions and outlook
In this work the rotational augmentation effects are studied for developing an alternative parametric formulation based on the baseline version of Chaviaropoulos and Hansen [16]. The novel formulation is calibrated based on the results of CFD simulations.
Starting from validated simulation setups, the finiteness and rotational augmentation effects are estimated by com-paring rotating CFD results against two-dimensional polars corrected for tip and root effects prior to the application of RAM selected from the literature.
The best agreements for both lift and drag coefficients were found for the model of Chaviaropoulos and Hansen [16], mainly at inboard regions where the rotational augmentation effects are stronger.
In order to extend the good performance of this model to a broader range of operating conditions and blade locations, a an alternative parametric formulation was derived and presented, together with its interaction in conjunction with the application of tip and root corrections. Subsequently, an optimization procedure based on curve fitting led to the determination of the optimal values for the a n parameter, resulting in a piece-wise modeling dependent on inflow wind velocity, i.e. operating tip speed ratio, and span-wise position. The novel formulation results in improved lift and drag coefficients prediction in all considered cases. The optimal configuration including the correction for tip effects of Shen shows better agreements at the very tip of the blade. One key point is to verify the proposed alternative parametric formulation in terms of independence with respect to the choice of the wind turbine model. The authors have reported some preliminary observations concerning the applicability of the proposed procedure for larger wind turbine rotor blades for the virtual model of the AVATAR rotor. With this regard, further issues should be addressed in future works, namely those concerning the determination of the initial lift (and drag) polar curve to be used as input for the application of the RAM. In fact, the usage of thicker aerodynamic shapes in the inboard part of the blade introduces large deviations from the standard thin airfoil section lift polar characteristics at low angle of attacks. This circumstance may be explored by means of additional DES simulations in order to capture with better approximation the highly separated region [44], possibly by means of laminar-turbulent transition modeling, in order to determine the non-rotating flow condition at the inboard region of the blade. Once these open points are clarified the application of the alternative parametric RAM formulation to these thick inboard sections may be conducted with reduced uncertainty. intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.