Dear,
I'm writing a Wind Generators FE model in Matlab. I'm implementing the BEM model you suggest in the AaeroDyn Theory Manual. I use the TipLoss Model, the HubLoss Model and the Buhl Correction for evaluate the axial induction factor (a) when CT (thrust coefficient) >0.96*F. I don't consider the Dynamic Stall Model and the Skewed Wake Correction (I assume the wind flow always perpendicular to the rotor plane).
In the iterative process to evaluate the axial induction factor (a) and the tangential induction factor (a'), i encounter these problems:
1) Often when the local angle of attack (α) is negative, CT results negative and for the relation [27] page 10 of the AeroDyn theory Manual, the axial induction factor (a) results >1 and the local flow angle (ϕ) results negative. So the TipLoss Factor and the HubLoss Factor results complex number! I have to evaluate the TipLoss Factor and the HubLoss Factor assuming the absolute value of the local flow angle, i.e. ϕ ?
Furthermore, the pressure normal to the rotor plane acting on this blade portion results of opposite direction respect to the wind flow, is physically correct ?
2) Sometimes, the iterative process doesn't converge to calculate the axial induction factor (a) ? So i force the convergence. Have you ever encounter this event ? How can i do ? Which could be the reason ?
3) If the wind flow i considered stochastic, do you assume U∞ (= mean wind speed) + the fluctuating contribution of the wind speed for the calculation of the axial induction factor ? In your relations doesn't appear the fluctuating part of the wind speed, why ?
I really thank you for your support, best regards
Lorenzo Montanari
BEM method Implementation
Moderator: Bonnie.Jonkman

 Posts: 15
 Joined: Mon Jul 04, 2011 4:06 am
 Organization: University of Parma
 Location: ITALY

 Posts: 437
 Joined: Fri Oct 21, 2005 10:22 am
 Organization: NREL/NWTC
 Location: Boulder, CO
 Location: Boulder, CO
 Contact:
Re: BEM method Implementation
Lorenzo,
I moved your topic to the Rotor Aerodynamics section because it really has nothing to do with testing. One of our interns (Dave Maniaci of Penn State) wrote a paper on convergence issues for WT_Perf. You should read it:
http://wind.nrel.gov/designcodes/papers/AIAA2011150703.pdf
I moved your topic to the Rotor Aerodynamics section because it really has nothing to do with testing. One of our interns (Dave Maniaci of Penn State) wrote a paper on convergence issues for WT_Perf. You should read it:
http://wind.nrel.gov/designcodes/papers/AIAA2011150703.pdf
Mr. Marshall L. Buhl Jr.
NWTC3811
National Renewable Energy Laboratory
Golden, CO 80401 USA
Marshall.Buhl@nrel.gov
Voice: +1 (303) 3846914
Cell: +1 (303) 9156623
Fax: +1 (303) 3847079
NWTC3811
National Renewable Energy Laboratory
Golden, CO 80401 USA
Marshall.Buhl@nrel.gov
Voice: +1 (303) 3846914
Cell: +1 (303) 9156623
Fax: +1 (303) 3847079

 Posts: 15
 Joined: Mon Jul 04, 2011 4:06 am
 Organization: University of Parma
 Location: ITALY
Re: BEM method Implementation
Dear Marshall Buhl,
I have read the paper you suggest to me; and I have insert in my Matlab program the sine curve singularity transition model but still it doesn't converge? The problem is with low wind speed. If i have correctly understood with that modification now the your program converge always ?
Is it necessary to write your new iterative process with a twodimensional NewtonRaphson iteration routine, iterating on both a and a‟ rather than separating them into two separate iteration loops ?
Does AeroDyn use the BEM model routine that is implemented in WT_Perf ?
In FullField wind condition, when you calculate the axial induction factor (a) and the tangential induction factor (a') do you use the local wind velocity in the current element (U_local = U_mean + U_fluctuating in the current blade element) ?
In reference to the AeroDyn Theory Manual formules (from page 4 to page 10), I don't understand exactly if in the term U∞ (= mean wind speed) AeroDyn considers also the fluctating component of the wind ?
Doesn't AeroDyn calculate a and a' factors also in function of the local fluctuating speed acting on the current blade element (U_fluctuating) ? Is the induced speed function only of the U_mean of the wind ?
And, in upwind rotor condition, do you add to Ulocal the term  Utower_shadow_effect ((Ulocal = Umean + Ufluctuating  Utower_shadow_effect in the current blade element) ?
a and a' in AeroDyn are function of the wind speed U∞ = Umean + Ufluctuating  Utower_shadow_effect in the current blade element ?
I really thank you for your fundamental help, best regards
Lorenzo Montanari
I have read the paper you suggest to me; and I have insert in my Matlab program the sine curve singularity transition model but still it doesn't converge? The problem is with low wind speed. If i have correctly understood with that modification now the your program converge always ?
Is it necessary to write your new iterative process with a twodimensional NewtonRaphson iteration routine, iterating on both a and a‟ rather than separating them into two separate iteration loops ?
Does AeroDyn use the BEM model routine that is implemented in WT_Perf ?
In FullField wind condition, when you calculate the axial induction factor (a) and the tangential induction factor (a') do you use the local wind velocity in the current element (U_local = U_mean + U_fluctuating in the current blade element) ?
In reference to the AeroDyn Theory Manual formules (from page 4 to page 10), I don't understand exactly if in the term U∞ (= mean wind speed) AeroDyn considers also the fluctating component of the wind ?
Doesn't AeroDyn calculate a and a' factors also in function of the local fluctuating speed acting on the current blade element (U_fluctuating) ? Is the induced speed function only of the U_mean of the wind ?
And, in upwind rotor condition, do you add to Ulocal the term  Utower_shadow_effect ((Ulocal = Umean + Ufluctuating  Utower_shadow_effect in the current blade element) ?
a and a' in AeroDyn are function of the wind speed U∞ = Umean + Ufluctuating  Utower_shadow_effect in the current blade element ?
I really thank you for your fundamental help, best regards
Lorenzo Montanari

 Posts: 437
 Joined: Fri Oct 21, 2005 10:22 am
 Organization: NREL/NWTC
 Location: Boulder, CO
 Location: Boulder, CO
 Contact:
Re: BEM method Implementation
Lorenzo,
First, there is not always a solution to the full BEM equations. If you can't get it to converge, eliminate the tangential induction.
You can iterate a and a' separately, but it will probably take longer to converge.
The AeroDyn implementation of BEM is slightly different and it incorrectly applies the tip loss factor to the tangential induction.
I don't know the answers to the rest of the questions, and, sadly, do not have time to investigate. Sorry. If I did, I would review the AeroDyn theory manual and also read the code, which is included in the archive.
Marshall
First, there is not always a solution to the full BEM equations. If you can't get it to converge, eliminate the tangential induction.
You can iterate a and a' separately, but it will probably take longer to converge.
The AeroDyn implementation of BEM is slightly different and it incorrectly applies the tip loss factor to the tangential induction.
I don't know the answers to the rest of the questions, and, sadly, do not have time to investigate. Sorry. If I did, I would review the AeroDyn theory manual and also read the code, which is included in the archive.
Marshall
Mr. Marshall L. Buhl Jr.
NWTC3811
National Renewable Energy Laboratory
Golden, CO 80401 USA
Marshall.Buhl@nrel.gov
Voice: +1 (303) 3846914
Cell: +1 (303) 9156623
Fax: +1 (303) 3847079
NWTC3811
National Renewable Energy Laboratory
Golden, CO 80401 USA
Marshall.Buhl@nrel.gov
Voice: +1 (303) 3846914
Cell: +1 (303) 9156623
Fax: +1 (303) 3847079
Return to “Rotor Aerodynamics”
Who is online
Users browsing this forum: No registered users and 1 guest