Generalized Dynamic Wake 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:

Generalized Dynamic Wake Theory and Implementation

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

This topic is dedicated to the discussion of GDW 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

Scott.Larwood
Posts: 98
Joined: Thu May 17, 2007 9:24 am
Location: Stockton, CA

Turning of dynamic wake at 9 m/s

Postby Scott.Larwood » Tue Sep 23, 2008 10:42 am

Hi,

I've noticed that I often get unstable results for average wind speeds at 9 m/s, and I have to re-run these cases with equilibrium inflow. Perhaps the cutoff for the switch to equilibrium should be raised from 8 m/s to 9 m/s until the problem with the dynamic wake at low speeds gets fixed.

Regards
--
Scott Larwood
Assistant Professor of Mechanical Engineering
University of the Pacific
3601 Pacific Avenue
Stockton, CA 95211

Email: slarwood@pacific.edu
Phone: 209.946.3080
Website: http://web.pacific.edu/x19551.xml

Pat.Moriarty
Posts: 40
Joined: Mon Nov 07, 2005 4:59 am
Location: NWTC/NREL

GDW velocity cutoff

Postby Pat.Moriarty » Tue Sep 23, 2008 11:38 am

Hi Scott,

You're welcome to do this yourself in line 5442 in AeroSubs.f90

IF (DYNINFL .AND. UBar / ConvFact <= 8.0) THEN

by changing 8.0 to 9.0 or any number you prefer. This essentially forces it to run in BEM mode, which you are doing manually now. So it may be easier to just do it manually if you don't want to recompile.

The next version of AeroDyn will have a more elegant correction (vortex ring correction) to this problem. We're currently testing the stability of the algorithms. Until then we're stuck with this implementation.

cheers,
Pat
Big whirls have little whirls, That feed on their velocity; And little whirls have lesser whirls, And so on to viscosity - Lewis Richardson.

Scott.Larwood
Posts: 98
Joined: Thu May 17, 2007 9:24 am
Location: Stockton, CA

Postby Scott.Larwood » Tue Sep 23, 2008 11:56 am

Thank you for the reply Pat. I probably stick to manual re-runs unless it becomes a burden. Good luck with the new algorithms; I hope that it makes for interesting work.
--
Scott Larwood
Assistant Professor of Mechanical Engineering
University of the Pacific
3601 Pacific Avenue
Stockton, CA 95211

Email: slarwood@pacific.edu
Phone: 209.946.3080
Website: http://web.pacific.edu/x19551.xml

Scott.Larwood
Posts: 98
Joined: Thu May 17, 2007 9:24 am
Location: Stockton, CA

Postby Scott.Larwood » Mon Sep 29, 2008 6:29 pm

Hi Pat,

Just writing to let you that I couldn't take it any longer so I made the code change you suggested and it works fine. No more wasted runs, no more manual re-runs.

Regards,
--
Scott Larwood
Assistant Professor of Mechanical Engineering
University of the Pacific
3601 Pacific Avenue
Stockton, CA 95211

Email: slarwood@pacific.edu
Phone: 209.946.3080
Website: http://web.pacific.edu/x19551.xml

Marco.Capuzzi
Posts: 23
Joined: Sun May 09, 2010 6:44 am

Re: Generalized Dynamic Wake Theory and Implementation

Postby Marco.Capuzzi » Fri May 11, 2012 9:28 am

Hi everyone,

I have a question about the GDW model, because I'm implementing such a model.

In the "AeroDyn Theory Manual", at pg.18 there is expression 88.

Now, my question is: do the signs of the quantities in this expression refer to the figure 12?
So for a WT operating in a standard idealized condition, being the induced flow opposite to the incoming flow direction,
the total inflow velocity is equal to the difference between the inflow due to free-stream wind and the induced flow?

Is this correct?

Thanks

Marco

Rick.Damiani
Posts: 5
Joined: Wed Nov 02, 2011 4:18 pm
Organization: NWTC/NREL
Location: CO

Re: Generalized Dynamic Wake Theory and Implementation

Postby Rick.Damiani » Mon May 14, 2012 11:23 am

Hi Marco,

We do have some inconsistencies in the signs for GDW, i.e. sometimes Aerodyn follows the helicopter conventions, and sometimes the wind turbine conventions, but in the end via double minus signs things work out just fine. We have in the plans to take care of that and make nomenclature uniform, in fact in Peters's model the sign convention can take care of all situations from helicopter to wind turbine states.
Figure 12 is confusing, and we are planning on replacing that as well, I have produced another figure in a work I did for the lab in the past that is meant to replace it.
Aerodyn assumes lamda_m<0 and lamda_f>0.

The more generic Peters's model assumes those quantities positive when pointed towards incoming flow (upwdind).
So in turbine state, lambda_f<0 and lambda_m>0 (both assumed positive when pointed towards incoming flow).

Hope that helps,
Rick

Marco.Capuzzi
Posts: 23
Joined: Sun May 09, 2010 6:44 am

Re: Generalized Dynamic Wake Theory and Implementation

Postby Marco.Capuzzi » Tue May 15, 2012 8:41 am

Hi Rick,

yes your answer helps and I think I understood Aerodyn's convention.

thanks

Marco

Marco.Capuzzi
Posts: 23
Joined: Sun May 09, 2010 6:44 am

Re: Generalized Dynamic Wake Theory and Implementation

Postby Marco.Capuzzi » Fri May 18, 2012 5:58 am

Hi everyone,

I have a question about the GDW (Generalized Dynamic Wake) model implemented within Aerodyn.

In the manual It is stated that the initial values of induced flow are obtained from the BEM routine. Then, the software switches on the GDW model, by using
the BEM solution as initial value for the differential equations, which define the GDW modelling.
Now, clearly you need to project the BEM solution on to the expansion shape functions and the azimuthal modes in order to calculate the initial values in the
appropriate reference.

My question is: when projecting the BEM solution, do you consider the inflow value upon a certain blade element as distributed over the respective annular surface
or is it considered as distributed only over the part of the rotor occupied by the blades (as it is done for the thrust, when calculating the pressure coefficients for
the GDW model)?

Thanks,

Marco

Khanh.Nguyen
Posts: 8
Joined: Wed Mar 14, 2012 2:14 pm
Organization: NWTC
Location: Boulder, CO, USA

Re: Generalized Dynamic Wake Theory and Implementation

Postby Khanh.Nguyen » Wed Jun 20, 2012 9:13 am

Marco,

The initial value for the GDW is the average of the BEM solution, and the average is over the blade radius.

Khanh

Marco.Capuzzi
Posts: 23
Joined: Sun May 09, 2010 6:44 am

Re: Generalized Dynamic Wake Theory and Implementation

Postby Marco.Capuzzi » Thu Jun 21, 2012 9:43 am

Hi Khanh,

thanks for the reply; I understand that you consider the average over the radius, but this value is considered as
constant also over the azimuth? i.e. do you consider the initial value of induced flow to be different from zero only on the first harmonic (m=0)?

I have some questions about the GDW as described within the Aerodyn Theory Manual. Sorry if I'm picky or sometime not precise, but I'm try to
implement the same model and I need to use the correct expressions.

1) In expressions 92,93 (pg.19), the azimuthal coordinate is left as an independent variable. Shouldn't it refer to the azimuthal coordinate of the considered blade?
What I mean is, the sine and cosine should not be calculate as sin(m*ψ_q) and cos(m*ψ_q)? Where "q" refers to the blade identification number.

2) In the same expressions you are considering a summation to approximate the integral of the lift distribution multiplied by the shape function over the dimensionless blade's length.
Why the term "∆r", approximating the infinitesimal variation of the radius "dr", is not present in the expression? Are you assuming an uniform aerodynamic mesh? Because, if the mesh is not
uniform, even though the radial coordinate is dimensionless and varies between zero and one, I thought you had to consider this term in the approximated numeric integration.

3) By comparing expressions 92-93 and the same terms as expressed in Peters (1991,Correlation of measured induced velocities with a finite-state wake model), I found a difference, which can be important.
In Peters, they consider the thrust per unit span (in N/m) and then they non-dimensionalize it. In the manual, the resultant thrust (in N) of a certain blade element is considered and then non-dimensionalized. If you compare
these two non-dimensionalizations, you have that in the manual the thrust per unit span is equal to the elements' thrusts divided by the blade's length (indeed, in your expression the R term at the denominator
has an exponent 4 instead of 3, as in Peters). Clearly, this is dimensionally correct, but it yields different results, because Peters thrust per unit span is equal to the element's thrust divided by its radial width.
Can you try to clarify this? Probably I misunderstood something, but I do not see what it is.

Thanks in advance,

Marco

Marco.Capuzzi
Posts: 23
Joined: Sun May 09, 2010 6:44 am

Re: Generalized Dynamic Wake Theory and Implementation

Postby Marco.Capuzzi » Fri Jun 22, 2012 9:17 am

Hi Khanh,

with regard to questions 2 and 3 in my last post, you can ignore them. I redid all of the calculations again and I found why
I was obtaining different results. Sorry!!

Question 1 still stands.

Regards

Marco

Khanh.Nguyen
Posts: 8
Joined: Wed Mar 14, 2012 2:14 pm
Organization: NWTC
Location: Boulder, CO, USA

Re: Generalized Dynamic Wake Theory and Implementation

Postby Khanh.Nguyen » Fri Jun 29, 2012 1:54 pm

Marco,

The initial condition of the GDW model requires the inflow value (calculated with BEM in this case) at a rotor azimuth (or at one instance in time). Thus, there is no need for the azimuth-averaged value.

Regarding question (1), you are correct. The variable ψ in Eqs. [92] and [93] in the Theory Manual should be corrected to ψ_q to be consistent with the outer summation of q, from 1 to B, B being the blade number.

Thanks,

Khanh

James.Nichols
Posts: 2
Joined: Tue Feb 03, 2015 2:39 am
Organization: Lloyd's Register EMEA
Location: UK

Re: Generalized Dynamic Wake Theory and Implementation

Postby James.Nichols » Thu Feb 05, 2015 8:54 am

Hi Everyone,

I have been doing some work to compare CFD models of rotor dynamics with dynamic wake models. When using the GDW model, I found that the radial distribution of inflow factor was not fine enough to make an accurate comparison. I was wondering how to include more radial mode shapes in the calculation.
I found the following parameter arrays which are used in the creation of the gamma matrix:
MRvector = (/ 0, 0, 1, 1, 2, 3 /)
NJVector = (/ 1, 3, 2, 4, 3, 4 /)
However, I'm not sure of the logic for assigning these arrays; could anyone give any tips for this?
Also, is the reason for the current choice of number of mode shapes governed by the speed of calculation, stability, or both?

Thanks

James


Return to “AeroDyn Development”

Who is online

Users browsing this forum: No registered users and 1 guest