Hello everybody,
For my master thesis I've developed a Matlab model of an onshore wind turbine, now I'm doing comparison between my model and FAST V8.16. I found a error in sidetoside tower top force (YawBrFyp) see figure.
I've made a comparison using Test 18 with the following changes:
all DOF set to "FALSE", except drivetrain and generator DOF;
0 deg shaft tilt;
0 deg yaw error;
0 wind power law exponent;
unsteady wind with: Umean component of 6.16 m/s, Vmean component of 0 m/s, Wmean component of 0 m/s.
In my model I have considered only the drag force to the nacelle:
Fyd= 0.5*Cd*D*L*V^2*rho
Where:
Cd = 1 (drag coefficient);
D = 5m (nacelle diameter);
L = 10 m (nacelle length);
V (Vcomponent of the wind)
rho = 1.2 kg/m^3 (air density).
I used BEM theory to model blade loads. Since I have three blades, I supposed that the Fy contribtion of each blade go to elide each other.
Are there phisical phenomena that I overlooked that justify the error between my model and FAST?
Thank you,
Riccardo.
sidetoside tower top force
Moderator: Bonnie.Jonkman

 Posts: 6
 Joined: Mon May 04, 2020 3:31 am
 Organization: Politecnico di Torino
 Location: Italia
sidetoside tower top force
 Attachments

 Fy comparison.jpg (76.17 KiB) Viewed 122 times

 Posts: 4777
 Joined: Thu Nov 03, 2005 4:38 pm
 Location: Boulder, CO
 Contact:
Re: sidetoside tower top force
Dear Riccardo,
It may help to plot the statistical distribution (histogram) of these time series, as well as the powerspectral density (PSD). From what I can tell, both FAST and your model have zero mean, but the standard deviation and tails of the distribution are larger from FAST. It is hard for me see from the time series how different the frequency content is between the models, but I would expect from FAST to see strong peaks at 3P, 6P, 9P etc., where "P" is per revolution; hopefully your model also captures these peaks?
I would guess the differences between FAST and your model are related to differences in how the aerodynamic loads and/or spatial variation of turbulence across the rotor are treated, but this is just a guess knowing nothing about your model.
Best regards,
It may help to plot the statistical distribution (histogram) of these time series, as well as the powerspectral density (PSD). From what I can tell, both FAST and your model have zero mean, but the standard deviation and tails of the distribution are larger from FAST. It is hard for me see from the time series how different the frequency content is between the models, but I would expect from FAST to see strong peaks at 3P, 6P, 9P etc., where "P" is per revolution; hopefully your model also captures these peaks?
I would guess the differences between FAST and your model are related to differences in how the aerodynamic loads and/or spatial variation of turbulence across the rotor are treated, but this is just a guess knowing nothing about your model.
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
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

 Posts: 6
 Joined: Mon May 04, 2020 3:31 am
 Organization: Politecnico di Torino
 Location: Italia
Re: sidetoside tower top force
Dear Jason,
Thanks for your prompt reply. As you suggested I plotted the power spectral density shown in figure.
As you can see from the figure, Fast presents the peaks you were talking about, while my model doesn't present peaks.
I guess that the reason of these differences is that FAST take into account some forces in Y direction derived from the 3bladed rotor revolution, while in my model I supposed nothing about this force for the reasons I wrote to you earlier.
If so, what effect I have neglected with my supposition?
Thank you for your help, best regards.
Riccardo
Thanks for your prompt reply. As you suggested I plotted the power spectral density shown in figure.
As you can see from the figure, Fast presents the peaks you were talking about, while my model doesn't present peaks.
I guess that the reason of these differences is that FAST take into account some forces in Y direction derived from the 3bladed rotor revolution, while in my model I supposed nothing about this force for the reasons I wrote to you earlier.
If so, what effect I have neglected with my supposition?
Thank you for your help, best regards.
Riccardo
 Attachments

 psd.jpg (57.64 KiB) Viewed 110 times

 Posts: 4777
 Joined: Thu Nov 03, 2005 4:38 pm
 Location: Boulder, CO
 Contact:
Re: sidetoside tower top force
Dear Riccardo,
I don't know any details about your model, so, I can't really comment on what is missing.
But the 3P, 6P, 9P, etc. peaks from OpenFAST are the result of the aerodynamic loads being calculated locally to each node along each blade as the blade passes through spatially varying turbulence, i.e., the socalled rotational sampling of turbulence. Does your model account for that?
Best regards,
I don't know any details about your model, so, I can't really comment on what is missing.
But the 3P, 6P, 9P, etc. peaks from OpenFAST are the result of the aerodynamic loads being calculated locally to each node along each blade as the blade passes through spatially varying turbulence, i.e., the socalled rotational sampling of turbulence. Does your model account for that?
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
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

 Posts: 6
 Joined: Mon May 04, 2020 3:31 am
 Organization: Politecnico di Torino
 Location: Italia
Re: sidetoside tower top force
Dear Jonkman,
Thank you for your reply. In my model I don't take into account any sidetoside Force when I calculate the load at each node of each blade, beacause I neglected that the wind varies spatially due the turbolence, which entails an imbalance of the loads between the blades. In fact I take into account anly the force in X direction and the tangential loads of each blade to calculate the Rotor torque. Now I understand the reason of the difference between FAST and my model.
Thank you for your help,
Best regards.
Riccardo.
Thank you for your reply. In my model I don't take into account any sidetoside Force when I calculate the load at each node of each blade, beacause I neglected that the wind varies spatially due the turbolence, which entails an imbalance of the loads between the blades. In fact I take into account anly the force in X direction and the tangential loads of each blade to calculate the Rotor torque. Now I understand the reason of the difference between FAST and my model.
Thank you for your help,
Best regards.
Riccardo.
Return to “Rotor Aerodynamics”
Who is online
Users browsing this forum: No registered users and 1 guest