Unsteady Aerodynamics

Discuss the theory and modeling of rotor aerodynamics.

Moderator: Bonnie.Jonkman

Mehran.Saeedi
Posts: 6
Joined: Wed May 13, 2020 2:04 pm
Organization: self
Location: Germany

Unsteady Aerodynamics

Postby Mehran.Saeedi » Mon May 18, 2020 12:17 am

Hi everyone,

I have a question on the unsteady aerodynamic (dynamic stall) model.

Currently, I am running stand-alone AeroDyn test cases to get to know the tool and how it works.
I've noticed that activating the UA model causes a drop in Cl and as a result a drop of the predictid AeroPwr.
The BEM part should be fine, becaues the converged AOA and induction coefficients are not affected by switching UA on.
My initial guess was that is should be due to worng dynamic stall parameters (which are obtained from AirfoilPrep, so I played with them and the results don't seem to be that sensitive to the inputs. I always get the remarkable drop in lift and power in the second time step, please see the attached pictures.

UA inactive:
UA_Off.png
UA_Off.png (47.25 KiB) Viewed 180 times

UA active:
UA_On.png
UA_On.png (32.34 KiB) Viewed 180 times

Any hint is highly appreciated.

BR,
Mehran

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

Re: Unsteady Aerodynamics

Postby Jason.Jonkman » Mon May 18, 2020 7:36 am

Dear Mehran,

I would normally expect the unsteady airfoil aerodynamics (UA) model of FAST / AeroDyn influence the variations of lift about the mean, but not to directly influence the mean lift itself. Your results show an unexpected impact of UA on the mean lift.

Which version of FAST / OpenFAST and which UA model of AeroDyn are you using? What UA parameters have been set?

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

Mehran.Saeedi
Posts: 6
Joined: Wed May 13, 2020 2:04 pm
Organization: self
Location: Germany

Re: Unsteady Aerodynamics

Postby Mehran.Saeedi » Tue May 19, 2020 1:03 am

Dear Jason,

thanks for your reply. I was not anticipating a change in mean values either.

My results are from AeroDyn_v15.04.02 (standalone mode).
I am using the Beddoes-Leishman unsteady model with Gonzalez's method (UAMod = 2)
The UA parametes are the following:
-6.0 alpha0
9.0 alpha1
-9.0 alpha2
0.9 eta_e
5.0994 C_nalpha
"Default" T_f0
"Default" T_V0
"Default" T_p
"Default" T_VL
"Default" b1
"Default" b2
"Default" b5
"Default" A1
"Default" A2
"Default" A5
0 S1
0 S2
0 S3
0 S4
1.5339 Cn1
-0.70 Cn2
"Default" St_sh
0.0180 Cd0
0 Cm0
0 k0
0 k1
0 k2
0 k3
0 k1_hat
"Default" x_cp_bar
"Default" UACutout
200 filtCutOff

They originate from the AirfoilPrep excel-tool.

Best regards,
Mehran

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

Re: Unsteady Aerodynamics

Postby Jason.Jonkman » Tue May 19, 2020 7:17 am

Dear Mehran,

Although filtCutOff = 200 Hz seems quite high, I don't see any obvious errors in your other input. Do you get improved results by setting filtCutOff = "DEFAULT" (20 Hz), or lower? Do you have similar issues when using other UAMod options, e.g. UAMod = 3?

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

Mehran.Saeedi
Posts: 6
Joined: Wed May 13, 2020 2:04 pm
Organization: self
Location: Germany

Re: Unsteady Aerodynamics

Postby Mehran.Saeedi » Wed May 20, 2020 6:23 am

Dear Jason,

thanks for your hints and for checking my input.
I've tried UAMod=3 and also reducing filtCutOff down to 2 Hz. But I am still getting the same results.

Best regards,
Mehran Saeedi

Mehran.Saeedi
Posts: 6
Joined: Wed May 13, 2020 2:04 pm
Organization: self
Location: Germany

Re: Unsteady Aerodynamics

Postby Mehran.Saeedi » Wed May 20, 2020 8:05 am

Dear Jason,

I was checking the AeroDyn 5MW benchmark, where I noticed that the same issue happens there as well.
Please have a look at the figure below. It compares the results with UA being switched On and Off.
They are generated from the default AeroDyne example: "NRELOffshrBsline5MW_Onshore_AeroDyn15_Test01.dat"
I've not made any changes in the input files rather than activating/deactivating UA.
UA_comparison.png
UA_comparison.png (114.32 KiB) Viewed 147 times


Activating UA causes a remarkable drop in RtAeroPwr and the lift coeffienct in two blade sections out of the plotted four.


Best regards,
Mehran

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

Re: Unsteady Aerodynamics

Postby Jason.Jonkman » Wed May 20, 2020 1:47 pm

Dear Mehran,

Perhaps the UA hysteresis loop is not symmetrically balanced about the static lift curve for this airfoil and UA settings, such that adding UA changes the mean lift? Have you tried plotting the UA hysteresis loop and the static lift curve versus angle of attack on the same plot?

It may help to use the UA driver to understand this behaviour. If you are unaware, there is a yet-to-be-documented standalone driver for running only the UA submodel of AeroDyn v15 in isolation that is included in the AeroDyn v15 and OpenFAST archives. You may want to play around with this to better understand the influence of UA.

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

Mehran.Saeedi
Posts: 6
Joined: Wed May 13, 2020 2:04 pm
Organization: self
Location: Germany

Re: Unsteady Aerodynamics

Postby Mehran.Saeedi » Fri May 22, 2020 1:38 am

Dear Jason,

it is an usymmetric hysteresis loop in deed. Here the comparison between static lift and UA lift for B1N4 and B1N5:
N4:
DU30_A17_10Hz.PNG
DU30_A17_10Hz.PNG (10.1 KiB) Viewed 131 times

N5:
DU25_A17_10Hz.PNG
DU25_A17_10Hz.PNG (9.41 KiB) Viewed 131 times


I am a newbie in UA and should definitely dig deeper into its theory. Do you think it makes sense (physically) for the unsteady lift coefficient to drift away from the static one? Or can it be due to not prefectly set UA parameters in the airfoil input file?

I'll try the standalone UA driver too, thanks for the hint.

Best regards,
Mehran

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

Re: Unsteady Aerodynamics

Postby Jason.Jonkman » Fri May 22, 2020 7:00 am

Dear Mehran,

Well, that explains why enabling UA shifts the mean lift. I would expect inconsistencies between the UA parameters and the static airfoil data to be the cause. Are you still using unmodified airfoil data from the AeroDyn v15.04.02 CertTest? For example, are you using unmodified NRELOffshrBsline5MW_Onshore_Test01.dvr, NRELOffshrBsline5MW_Onshore_AeroDyn15_Test01.dat, NRELOffshrBsline5MW_AeroDyn_blade.dat, and DU25_A17_10Hz.dat?

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

Mehran.Saeedi
Posts: 6
Joined: Wed May 13, 2020 2:04 pm
Organization: self
Location: Germany

Re: Unsteady Aerodynamics

Postby Mehran.Saeedi » Mon May 25, 2020 12:39 am

Dear Jason,

from the files you've mentioned, I'va just playes with 'NRELOffshrBsline5MW_Onshore_AeroDyn15_Test01.dat' to activate/deactivate UA. The rest are left untouched.


Best regards,
Mehran

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

Re: Unsteady Aerodynamics

Postby Jason.Jonkman » Tue May 26, 2020 2:17 pm

Dear Mehran.

OK, thanks for confirming. I've now posted this problem as an issue in the OpenFAST r-test repository: https://github.com/OpenFAST/r-test/issues/20. Hopefully this will start the process to get this bug in the airfoil data input file(s) fixed. Please feel free to add further comments there. And let us know if you correct the airfoil data input file(s) yourself.

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 “Rotor Aerodynamics”

Who is online

Users browsing this forum: No registered users and 1 guest