Dynamic Stall Theory and Implementation

This forum is dedicated to discussions of the redevelopment of AeroDyn.

Moderators: Bonnie.Jonkman, Jason.Jonkman

Marshall.Buhl
Posts: 437
Joined: Fri Oct 21, 2005 10:22 am
Organization: NREL/NWTC
Location: Boulder, CO
Location: Boulder, CO
Contact:

Dynamic Stall Theory and Implementation

Postby Marshall.Buhl » Fri Feb 08, 2008 11:36 am

This topic is dedicated to the discussion of dynamic-stall theory and implementation.
Mr. Marshall L. Buhl Jr.
NWTC-3811
National Renewable Energy Laboratory
Golden, CO 80401 USA
Marshall.Buhl@nrel.gov
Voice: +1 (303) 384-6914
Cell: +1 (303) 915-6623
Fax: +1 (303) 384-7079

Arash.Mohammady
Posts: 4
Joined: Thu May 30, 2013 3:40 am
Organization: DTU
Location: Denmark

Re: Dynamic Stall Theory and Implementation

Postby Arash.Mohammady » Tue Feb 18, 2014 4:50 am

Hi,

I am trying to add Stig Øye Dynamic stall model into AeroDyn and I would like to use stall angle of the airfoils in my code. Unfortunately, I did not manage to find where/how you read it in the AeroDyn source code. I could see that you read "Beddoes stall parameters" such as "zero lift angle of attack, AOL", "CNA", and etc. in READFL subroutine (AeroSubs.f90) but there is no stall angle of attack there.

I would really appreciate it if you could help me in this matter.

Thank you very much for you time.


Best regards,
Arash
Attachments
Capture.PNG
READFL subroutine
Capture.PNG (48.05 KiB) Viewed 9740 times

Marshall.Buhl
Posts: 437
Joined: Fri Oct 21, 2005 10:22 am
Organization: NREL/NWTC
Location: Boulder, CO
Location: Boulder, CO
Contact:

Re: Dynamic Stall Theory and Implementation

Postby Marshall.Buhl » Tue Feb 18, 2014 10:12 am

Arash,

I believe that input parameter is a relic of the old days when AeroDyn included the Gormont (sp?) dynamic-stall model. If you look at Page 21 of the "AeroDyn User's Guide," you will see that the (unnamed) parameter is "reserved for future use and backward compatibility."

I am currently redoing the input files and those relics will no longer be in them.


Marshall
Mr. Marshall L. Buhl Jr.
NWTC-3811
National Renewable Energy Laboratory
Golden, CO 80401 USA
Marshall.Buhl@nrel.gov
Voice: +1 (303) 384-6914
Cell: +1 (303) 915-6623
Fax: +1 (303) 384-7079

Arash.Mohammady
Posts: 4
Joined: Thu May 30, 2013 3:40 am
Organization: DTU
Location: Denmark

Re: Dynamic Stall Theory and Implementation

Postby Arash.Mohammady » Wed Feb 19, 2014 5:19 am

Hi Marshall,

Thanks for your help. These are the lines that you mentioned they are no longer used, right?

[attachment=0]Capture.PNG[/attachment]

And I guess i need to read the 5th line of AeroDyn input myself in order to have Stall angle of attack.


Best regards,
Arash
Attachments
Capture.PNG
Capture.PNG (18.04 KiB) Viewed 9724 times

Marshall.Buhl
Posts: 437
Joined: Fri Oct 21, 2005 10:22 am
Organization: NREL/NWTC
Location: Boulder, CO
Location: Boulder, CO
Contact:

Re: Dynamic Stall Theory and Implementation

Postby Marshall.Buhl » Wed Feb 19, 2014 10:04 am

Dear Arash,

Those lines look like the ones for the unused parameters. As you are modifying the code, you can do whatever you like. You can modify one to read in the stall angle (it doesn't have to be the 5th line) and delete the others. You would have to delete the lines in your input files if you did that. It's really up to you.

Please realize that you will have to make a similar effort if you want to use future versions of AeroDyn. The next version will not have any of those lines.


Marshall
Mr. Marshall L. Buhl Jr.
NWTC-3811
National Renewable Energy Laboratory
Golden, CO 80401 USA
Marshall.Buhl@nrel.gov
Voice: +1 (303) 384-6914
Cell: +1 (303) 915-6623
Fax: +1 (303) 384-7079

Arash.Mohammady
Posts: 4
Joined: Thu May 30, 2013 3:40 am
Organization: DTU
Location: Denmark

Re: Dynamic Stall Theory and Implementation

Postby Arash.Mohammady » Thu Feb 20, 2014 4:52 am

Thanks a lot for your help Marshall.

Best regards,
Arash

Srinivasa.Sudharsan
Posts: 39
Joined: Mon Jun 08, 2015 10:59 pm
Organization: thiagarajar college of engineering
Location: India

Re: Dynamic Stall Theory and Implementation

Postby Srinivasa.Sudharsan » Tue Mar 22, 2016 10:02 pm

Hi,
I am trying to model the dynamics of the extreme gust of around 60m/sec.I am also using BEDDOES. But when i try to run the simulation i get the error at some instant of simulation as:
"FAST_Solution:CalcOutputs_And_SolveForInputs:SolveOption2:AD_CalcOutput/Inflow:Inflow:infdist:MATINV:gaussj: Singular matrix encountered." and warning as:
""FAST_Solution:CalcOutputs_And_SolveForInputs:SolveOption2:AD_CalcOutput:ELEMFRC:BeddoesModel:ATTACH: Supersonic condition has
subsided with Blade #1 element #1.
AD_CalcOutput:ELEMFRC:BeddoesModel:ATTACH: Blade #1 element #5 is supersonic! Other elements are likely supersonic as well.
Supersonic mach nos. will be set to 0.7 to attempt continuation.
AD_CalcOutput:ELEMFRC:BeddoesModel:ATTACH: Supersonic condition has subsided with Blade #2 element #1.
AD_CalcOutput:ELEMFRC:BeddoesModel:ATTACH: Blade #2 element #10 is supersonic! Other elements are likely supersonic as well.
Supersonic mach nos. will be set to 0.7 to attempt continuation.
AD_CalcOutput:ELEMFRC:BeddoesModel:ATTACH: Supersonic condition has subsided with Blade #3 element #1.
AD_CalcOutput:ELEMFRC:BeddoesModel:ATTACH: Blade #3 element #9 is supersonic! Other elements are likely supersonic as well.
Supersonic mach nos. will be set to 0.7 to attempt continuation.""

can you please help with the solution for this problem so that i could complete my simulation to analyse the effect of extreme gust!!! what is the parameter to be changed to rectify the problem???? kindly help!!

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

Re: Dynamic Stall Theory and Implementation

Postby Jason.Jonkman » Wed Mar 23, 2016 8:37 am

Dear Srinivasa,

My guess is that your model is going numerically unstable for some reason. Warnings regarding a “small angle approximation violation,” warnings regarding “supersonic blades,” a simulation crash, or very large deflections that occur in the time series are good signs of a numerical instability. Without knowing more about your simulation settings, it's hard to know that what the problem is. Please use "Search..." in the upper right of the forum to search for "supersonic" to see how we've solved similar problems in the past.

If you have a gust around 60 m/s, I'm assuming you are modeling a situation with a parked/idling rotor with the blades feathered to 90 deg? What is a typical angle of attack? Please note that the BEDDOES model is applicable for angles of attack in the linear and stall regions of the airfoil. The model may not be accurate for much larger angles of attack.

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: Dynamic Stall Theory and Implementation

Postby Albert.Meseguer » Thu Aug 25, 2016 5:55 am

Hi everyone,

I have been reading the theory behind the different dynamic stall parameters and I'm a little bit confused about the definiton of alpha0 for the NACA64_A17.dat of the 5MW reference turbine. The value defined in the provided airfoil is -4.432. However It does not match with the provided definition. Using linear interpolation I obtain an angle of -3.8381 .

Is the given angle correct? If so, what is the theory behind it?

The rest of airfoils seems to match correctly,

Thank you,

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

Re: Dynamic Stall Theory and Implementation

Postby Jason.Jonkman » Tue Sep 06, 2016 5:09 pm

Dear Albert,

I see where you are getting the value of -3.8381, by using the angles of attack (AoAs) on either side of the 0-lift curve. The value of -4.432 came from using a best-fit line of the Cn versus AoA curve over a larger range of AoA, from 0 to 8 degrees, using the AirfoilPrep preprocessor: https://nwtc.nrel.gov/AirFoilPrep. The values are bit different due to a small amount of nonlinearity in the Cn versus AoA curve. The range of AoA from 0 to 8 degrees was somewhat arbitrary, and playing around with the AoA range in AirfoilPrep a bit, I can get values closer to yours. Have you noticed much difference in the results when changing this value?

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: Dynamic Stall Theory and Implementation

Postby Saptarshi.Sarkar » Wed Apr 05, 2017 9:58 am

Hi everyone,

Can somebody point me towards the theory of Beddoes-Leishman dynamic stall model which is used in FAST v15. Also, I have read in various places that the profile dependent parameters required in any dynamic stall model are generally obtained through experimental results. In that case, how is FAST assuming the parameters?

Thank you,

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

Re: Dynamic Stall Theory and Implementation

Postby Jason.Jonkman » Wed Apr 05, 2017 12:02 pm

Dear Saptarshi,

The implementation of the Beddoes-Leishman unsteady airfoil aerodynamics (UA) model of AeroDyn v15 is documented in the following draft report: https://wind.nrel.gov/nwtc/docs/UnsteadyAero_Manual.pdf (this is a draft report, and updates are in progress).

There are many UA model parameters specified in the airfoil data input files for AeroDyn v15, with DEFAULT settings for most parameters. See section 4.4 of the draft AeroDyn User's Guide and Theory Manual for more information: https://wind.nrel.gov/nwtc/docs/AeroDyn_Manual.pdf.

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: Dynamic Stall Theory and Implementation

Postby Saptarshi.Sarkar » Wed Apr 05, 2017 12:07 pm

Thanks a lot. I really appreciate it.

Marc.Winter
Posts: 0
Joined: Wed Oct 30, 2019 4:12 am
Organization: self
Location: Germany

Re: Dynamic Stall Theory and Implementation

Postby Marc.Winter » Wed Oct 30, 2019 8:04 am

In "The Unsteady Aerodynamics Module for FAST" (28th. February 2017), the onset of dyanmic stall is said to occur when the effective normal force coefficient C'_n exceeds a specific value C_n1 at leading edge separation. Since this is supposed to be used in HAWT applications, I assume this method is optimized for low Mach numbers (M < 0.3).
In "A Modified Dynamic Stall Model for Low Mach Numbers" by W. Sheng et al., which is also mentioned in the FAST 8 documentation, a different approach is chosen by delaying the effective AOA by a difference obtained from a linear dynamic_stall_onset_angle vs reduced_pitch_rate-fit from airfoil moment break data. This produces good results on NACA 0012, which is one of the airfoils I am interested in.

Do both approaches perform similarly for low mach numbers or would one be preferable for a certain set of airfoils?

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

Re: Dynamic Stall Theory and Implementation

Postby Jason.Jonkman » Wed Oct 30, 2019 4:29 pm

Dear Marc,

Sorry, but I don't know enough about these two approaches to comment. And unfortunately, Rick Damiani--who led that report--is no longer at NREL. Hopefully someone more knowledgeable about these approaches can respond.

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


Return to “AeroDyn Development”

Who is online

Users browsing this forum: No registered users and 0 guests