Wake-Skewed Angle Model

Discuss the theory and modeling of rotor aerodynamics.

Moderator: Bonnie.Jonkman

Albert.Meseguer
Posts: 10
Joined: Tue Feb 23, 2016 9:51 am
Organization: self
Location: Spain

Wake-Skewed Angle Model

Postby Albert.Meseguer » Tue Aug 23, 2016 4:02 am

Hi everyone,

I am running some simulations where I compare the effect of the two available skewed-wake correction model in Aerodyn v15.03 (Uncoupled and Pitt and Peters model) . My case presents a wind turbine with no tilt neither blade pre-cone. I use an uniform wind speed without shear effect.

Equation 30 in the report " Development and Validation of a New Blade Element Momentum Skewed-Wake Model within Aerodyn" ( http://www.nrel.gov/docs/fy15osti/63217.pdf ) , models the Pitt and Peters correction as function of the induction factor for each section:

equation_document.PNG
equation_document.PNG (24.36 KiB) Viewed 8755 times


Being X0 the angle between the vector normal to the rotor plane, in my case the yaw angle. I expect to obtain same results when using Uncoupled and Pitt and Peters models for yaw angle 0 ( the tangent term becomes 0). However, the models differ and I really do not understand why. I.e. I plot thrust of the turbine for both models:

Error_Pitt_uncoupled.png
Error_Pitt_uncoupled.png (30.26 KiB) Viewed 8755 times


Any idea why? Am I missing something?

Thank you,
Albert Meseguer Urbán

Jason.Jonkman
Posts: 5723
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

Re: Wake-Skewed Angle Model

Postby Jason.Jonkman » Tue Aug 23, 2016 6:45 am

Dear Albert,

I have not tried to reproduce your results, but I quickly skimmed the AeroDyn BEMT algorithm and it is not obvious to me how changing input SkewMod between uncoupled (SkewMod=1) and Pitt and Peters (SkewMod=2) would have any influence on the results if the inflow-skew angle were in fact zero. Have you output the inflow-skew angle (AeroDyn output RtSkew) to verify that it is in fact zero? Have you also noticed the inflow angle (output BαNβPhi), angle of attack (output BαNβAlpha), or the axial and/or tangential induction factors (outputs BαNβAxInd and/or BαNβTnInd) changing for certain blade nodes between these cases?

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Albert.Meseguer
Posts: 10
Joined: Tue Feb 23, 2016 9:51 am
Organization: self
Location: Spain

Re: Wake-Skewed Angle Model

Postby Albert.Meseguer » Wed Aug 24, 2016 2:40 am

Dear Jason,

Thank you for your fast reply.

In both simulations the inflow-skew angle is 0 in all the simulation (as expected) . I took a look into the aerodynamic parameters you mentioned in your post , for nodes close to the root both induction factors (B1N1AxInd , B1N1TnInd) are the same. Therefore , the inflow and attack angle (B1N1Alpha , B1N1Phi) present a really good agreement.

However, If I output other blade nodes , i.e. node 9 , the tangential induction factor value (B1N9TnInd) for the uncoupled option remains -1 all the simulation. On the other hand, the Pitt and Peters model is 0 the whole simulation. The axial induction factor (B1N9AxInd) has also discrepancies, the Uncoupled model predicts a value of +1 the complete simulation while the Pitt and Peters is 0 . This has an impact in inflow and attack angle since I'm considering the tangential induction factor in the calculation (TanInd = TRUE). This results in this discrepancies over time :

Uncoupled_Pitt_aero.png
Uncoupled_Pitt_aero.png (32.09 KiB) Viewed 8744 times


I do not understand why the induction factors are different if the module is only supposed to correct the yaw misalignment. I'm missing something probably.

Any idea?

P.S. : The first node "B1N1..." corresponds to the first node of the blade and "B1N9..." is the last node of the blade, defined in BlOutNd as 1,3,5,9,11,13,15,17,19.
Albert Meseguer Urbán

Jason.Jonkman
Posts: 5723
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

Re: Wake-Skewed Angle Model

Postby Jason.Jonkman » Wed Aug 24, 2016 12:45 pm

Dear Albert,

OK, thanks for looking into this more.

I looked in the source code (specifically SUBROUTINE ApplySkewedWakeCorrection() in source file BEMTUncoupled.f90) and noticed that if the axial and tangential induction factors (a and ap) are set to 1 or -1, respectively, both are reset to zero. This was added to avoid numerical problems in the calculation of the inflow-angle (phi) later in the code. However, earlier in the code, a=1 and ap=-1 is set at the root and tip when the root and tip loss models are enabled to ensure that no loads are computed at this nodes. So, in essence, the skewed-wake correction is over-writing the tip and hub loss effects at the root and tip in your case. Obviously, there is a problem in the logic here that we must fix. I've noted this as an issue that will be fixed in a future version.

Thanks for your help in debugging.

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Saptarshi.Sarkar
Posts: 19
Joined: Fri Mar 24, 2017 10:09 am
Organization: Trinity College Dublin
Location: Ireland

Re: Wake-Skewed Angle Model

Postby Saptarshi.Sarkar » Wed Jun 20, 2018 8:25 am

Dear Dr. Jonkman,

For the Pitt and Peter's skewed wake correction is the angle between the vector normal to the rotor plane and the wind vector calculated at every blade section? What I mean is, under turbulent wind every blade section will have a different wind vector while the normal vector to the rotor plane remains same. In this situation should I calculate the angle for every section separately and use that angle for the skewed wake correction?

I would really appreciate if you could help clear this confusion.

Thanks a million,
Saptarshi

Jason.Jonkman
Posts: 5723
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

Re: Wake-Skewed Angle Model

Postby Jason.Jonkman » Wed Jun 20, 2018 1:53 pm

Dear Saptarshi,

The skew angle in the Pitt and Peters skewed-wake correction is defined as the angle between the vector normal to the rotor plane and the vector of the wind averaged across the rotor. This angle is the same for all blades.

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Saptarshi.Sarkar
Posts: 19
Joined: Fri Mar 24, 2017 10:09 am
Organization: Trinity College Dublin
Location: Ireland

Re: Wake-Skewed Angle Model

Postby Saptarshi.Sarkar » Thu Jun 21, 2018 3:58 am

Dear Dr. Jonkman,

Thanks a lot for clarifying this.

Saptarshi

Lorenzo.Cottura
Posts: 56
Joined: Wed Jul 15, 2020 9:04 am
Organization: Politecnico di Torino
Location: Italy

Re: Wake-Skewed Angle Model

Postby Lorenzo.Cottura » Thu Sep 03, 2020 3:46 am

Dear Jason,
I introduce myself to this topic because my question is about skew correction.
In particular, I was wondering the formulation of Xo (angle between the vector normal to the rotor plane) in case of non-constant wind and all 6 active platform DOFs.
Could you bring me the equation or a document in which it is reported?
Thanks for your help,
Lorenzo.

Jason.Jonkman
Posts: 5723
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

Re: Wake-Skewed Angle Model

Postby Jason.Jonkman » Thu Sep 03, 2020 8:54 am

Dear Lorenzo,

The equation for the skew angle is calculated as described, i.e. the angle between the vector normal to the rotor plane and the vector of the wind averaged across the rotor. In the source code, it is calculated as:

skew = ACOS( V_DiskAvg dot xhat_Disk / 2NORM ( V_DiskAvg ) )

where V_DiskAvg is a the relative velocity (wind minus structure) averaged over the rotor disk and xhat_Disk is a unit vector normal to the disk.

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Lorenzo.Cottura
Posts: 56
Joined: Wed Jul 15, 2020 9:04 am
Organization: Politecnico di Torino
Location: Italy

Re: Wake-Skewed Angle Model

Postby Lorenzo.Cottura » Mon Nov 09, 2020 11:02 am

Dear Jason,
I'm not sure what you mean by "wind averaged across the rotor". Tell me if I understood correctly:
By defining V_DiskAvg = [V0_x, V0_y, V0_z] the vector of the wind speeds relative to the structure (wind minus structure) in the node under analysis, to be clear the same used in the BEM theory to obtain Vx and Vy according to the formula reported in "Development and Validation of a New Blade Element Momentum Skewed-Wake Model within AeroDyn" article where Vinf=V0_x.
Schermata 2020-11-09 alle 18.51.53.png
Schermata 2020-11-09 alle 18.51.53.png (31.45 KiB) Viewed 1072 times

The formula becomes:

skew = acos (V0_x / (sqrt (V0_x ^ 2 + V0_y ^ 2 + V0_z ^ 2)))

Confirm?
Thanks for the help,
Lorenzo.

Jason.Jonkman
Posts: 5723
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

Re: Wake-Skewed Angle Model

Postby Jason.Jonkman » Mon Nov 09, 2020 11:31 am

Dear Lorenzo,

Equation (28) from Ning et al's paper is not implemented in FAST / AeroDyn, as this formula is only valid for a rigid wind turbine and uniform flow.

In AeroDyn, the V_DiskAvg is computed by spatially averaging the relative wind velocity at each aerodynamic analysis node across the rotor, i.e.:

V_DiskAvg = SUM( SUM( V_wind(j,k) - V_structural(j,k), j = 1, 2, ..., NumBlNds ), k = 1, 2, ..., NumBlades )/( NumBlades*NumBlNds )

I hope that helps.

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Lorenzo.Cottura
Posts: 56
Joined: Wed Jul 15, 2020 9:04 am
Organization: Politecnico di Torino
Location: Italy

Re: Wake-Skewed Angle Model

Postby Lorenzo.Cottura » Mon Nov 09, 2020 12:03 pm

Dear Jason,
So V_DiskAvg becomes a constant vector for each time step like: V_DiskAvg = [V_DiskAvg_x, V_DiskAvg_y, V_DiskAvg_z]
and it will be
skew = acos (V_DiskAvg / (sqrt (V_DiskAvg_x^2 + V_DiskAvg_y^2 + V_DiskAvg_z^2)))
It is right?
Best regards,
Lorenzo

Jason.Jonkman
Posts: 5723
Joined: Thu Nov 03, 2005 4:38 pm
Location: Boulder, CO
Contact:

Re: Wake-Skewed Angle Model

Postby Jason.Jonkman » Mon Nov 09, 2020 12:22 pm

Dear Lorenzo,

If xhat_Disk is directed along the x axis, then:

skew = acos (V_DiskAvg_x / (sqrt (V_DiskAvg_x^2 + V_DiskAvg_y^2 + V_DiskAvg_z^2)))

Best regards,
Jason Jonkman, Ph.D.
Senior Engineer | National Wind Technology Center (NWTC)

National Renewable Energy Laboratory (NREL)
15013 Denver West Parkway | Golden, CO 80401
+1 (303) 384 – 7026 | Fax: +1 (303) 384 – 6901
nwtc.nrel.gov

Lorenzo.Cottura
Posts: 56
Joined: Wed Jul 15, 2020 9:04 am
Organization: Politecnico di Torino
Location: Italy

Re: Wake-Skewed Angle Model

Postby Lorenzo.Cottura » Mon Nov 09, 2020 12:50 pm

Dear Jason,
Yes, I meant that, sorry for the typo.
Best regards,
Lorenzo.


Return to “Rotor Aerodynamics”

Who is online

Users browsing this forum: No registered users and 1 guest