Hi, all.
I am working on the controller for offshore floating wind turbines in Region 2 (between cutin and rated velocity), and using FAST in Simulink.
I would like to have an estimator for obtaining the effective/relative wind velocity (considering platform's moving...) at hub height.
Therefore, I chose the FAST’s output “RotTorq”, “RotSpeed” and the “Cqlameda curve” to solve the equation as follows by iteration.
Tr=0.5*pho*pi*R^3*Cq*Vrel^2
However, the result was questionable.
Then, I tried it on an onshore wind turbine, and prepared a wind file with a step from 8m/s to 6m/s. The result was still slow, and the estimator needed a transient process to move to 6m/s, about 20s.
Later, I tried it again on the onshore wind turbine using the same wind file above mentioned. The only difference is that I CLOSE the GenDOF. The result seemed good, and the estimator would suddenly move to 6m/s. But I still needed the GenDOF.
Could somebody tell me the reason? In my opinion, I use the torque and speed of the rotor, not the generator, so I cannot understand the result.
And could somebody tell me the method for obtaining the effective/relative wind speed? It is important to me.
Thanks in advance.
Best regards,
Lei DUAN
estimator of effective/relative wind velocity
Moderator: Bonnie.Jonkman

 Posts: 34
 Joined: Tue Apr 02, 2013 9:27 pm
 Organization: NREL
 Location: Colorado
Re: estimator of effective/relative wind velocity
Hi,
When you say you closed the generator DOF, do you mean you switched it to False? This would mean then the rotor speed is now fixed. If the rotor speed is not fixed, then changes in wind speed will cause transients while the rotor speed accelerates/decelerates to a new state, which I think is the transient you're seeing when you don't fix rotor speed.
That said, I do think it should be possible to develop a better estimator even with turbine dynamics, but I think it will be important to include turbine acceleration. See for example the discussion of the issue it the introduction of this paper:
http://iopscience.iop.org/17426596/75/1/012082
Paul
When you say you closed the generator DOF, do you mean you switched it to False? This would mean then the rotor speed is now fixed. If the rotor speed is not fixed, then changes in wind speed will cause transients while the rotor speed accelerates/decelerates to a new state, which I think is the transient you're seeing when you don't fix rotor speed.
That said, I do think it should be possible to develop a better estimator even with turbine dynamics, but I think it will be important to include turbine acceleration. See for example the discussion of the issue it the introduction of this paper:
http://iopscience.iop.org/17426596/75/1/012082
Paul

 Posts: 5
 Joined: Mon Oct 24, 2016 7:13 am
 Organization: Delft University of Technology
 Location: Delft
Re: estimator of effective/relative wind velocity
Hi to all,
for the purposes of my research project I am trying to develop an effective wind speed estimator in Simulink for the 5 MW Fast Wind Turbine, in line with the work http://iopscience.iop.org/17426596/75/1/012082 appearing in the above post. My progress so far, is summarized as follows:
 Assuming rigid drivetrain (DrTrDOF=False in ElastoDyn ), the Aerodynamic torque Ta (Qa in the paper above) is computed in Simulink through a Matlab function block by the following relation:
Ta=Jr*omegar_dot+LSShftTq,
where Jr denotes the rotor inertia (I took Jr=115.926 kg*m^2, hub inertia about low speed shaft found in "Definition of a 5MW Reference Wind Turbine for Offshore System Development"), omegar_dot is the acceleration of the rotor (I take RotSpeed as output from FAST convert it to rad/sec by multiplying with pi/30 and then taking the derivative using a derivative block in Simulink), and LSShftTq is given, again, as output from FAST.
 Then from
2*Ta/rho*pi*R^5*omegar^2=cp(beta,TSR)/TSR^3=F(beta,TSR)
if one computes Ta and omegar (aerodynamic torque and rotor speed), the leftside term of the quation above is known and F (beta,TSR) is computed.
The table F vs TSR can be computed offline, using the knowledge of Cp(beta,TSR), thus for a given F, the tipspped ratio TSR can be computed.To compute F=Cp/TSR^3 i have used the tables found in "https://wind.nrel.gov/forum/wind/viewtopic.php?f=2&t=582&p=6299&hilit=torque+coefficients#p6299". My results for Cp and F, for beta=0 when the wind turbine operates below rated power, are attached.
 Then, with TSR known the effective wind speed is given from V=omegar*R/TSR. I have developed all the previous steps in Simulink using a stepwise changing simulated wind speed, with wind speed below the rated wind speed. My concern is that, although the estimated wind speed converges to the simulated one there is a significant delay for the actual convergence which I can not explain. I have attached the comparison between the simulated wind speed and the estimated one. The actual response of the estimated wind speed resembles a lot the actual behavior of a lowpass filter could dynamic inflow explain this behavior? Do you see any mistake i have done in my modelling design? I have also attached the AeroDyn file I have used for the simulations.
I am looking forward to your responses, thank you in advance for your comments.
Bart
for the purposes of my research project I am trying to develop an effective wind speed estimator in Simulink for the 5 MW Fast Wind Turbine, in line with the work http://iopscience.iop.org/17426596/75/1/012082 appearing in the above post. My progress so far, is summarized as follows:
 Assuming rigid drivetrain (DrTrDOF=False in ElastoDyn ), the Aerodynamic torque Ta (Qa in the paper above) is computed in Simulink through a Matlab function block by the following relation:
Ta=Jr*omegar_dot+LSShftTq,
where Jr denotes the rotor inertia (I took Jr=115.926 kg*m^2, hub inertia about low speed shaft found in "Definition of a 5MW Reference Wind Turbine for Offshore System Development"), omegar_dot is the acceleration of the rotor (I take RotSpeed as output from FAST convert it to rad/sec by multiplying with pi/30 and then taking the derivative using a derivative block in Simulink), and LSShftTq is given, again, as output from FAST.
 Then from
2*Ta/rho*pi*R^5*omegar^2=cp(beta,TSR)/TSR^3=F(beta,TSR)
if one computes Ta and omegar (aerodynamic torque and rotor speed), the leftside term of the quation above is known and F (beta,TSR) is computed.
The table F vs TSR can be computed offline, using the knowledge of Cp(beta,TSR), thus for a given F, the tipspped ratio TSR can be computed.To compute F=Cp/TSR^3 i have used the tables found in "https://wind.nrel.gov/forum/wind/viewtopic.php?f=2&t=582&p=6299&hilit=torque+coefficients#p6299". My results for Cp and F, for beta=0 when the wind turbine operates below rated power, are attached.
 Then, with TSR known the effective wind speed is given from V=omegar*R/TSR. I have developed all the previous steps in Simulink using a stepwise changing simulated wind speed, with wind speed below the rated wind speed. My concern is that, although the estimated wind speed converges to the simulated one there is a significant delay for the actual convergence which I can not explain. I have attached the comparison between the simulated wind speed and the estimated one. The actual response of the estimated wind speed resembles a lot the actual behavior of a lowpass filter could dynamic inflow explain this behavior? Do you see any mistake i have done in my modelling design? I have also attached the AeroDyn file I have used for the simulations.
I am looking forward to your responses, thank you in advance for your comments.
Bart
 Attachments

 Aerodyn.JPG (139.57 KiB) Viewed 1872 times

 Cp and F.JPG (25.79 KiB) Viewed 1872 times

 Simulated wind speed vs estimated wind speed.jpg (51.41 KiB) Viewed 1872 times

 Posts: 3857
 Joined: Thu Nov 03, 2005 4:38 pm
 Location: Boulder, CO
 Contact:
Re: estimator of effective/relative wind velocity
Dear Bert,
Regarding your calculation of Ta:
Best regards,
Regarding your calculation of Ta:
 The hub inertia of the NREL 5MW turbine is 115926 kg*m^2, not 115.926 kg*m^2. However, the equation you are using requires the rotor inertia instead of the hub inertia, which for the NREL 5MW turbine is 38677040.613 kg*m^2.
 When using AeroDyn v15, you can get the aerodynamic applied torque directly as an output i.e. RotAeroMxh.
 Your AeroDyn model does not have any dynamics enabled, so, the delay is not coming from that.
 Instead, the delay is likely caused by structural dynamics i.e. the structural degreesoffreedom that are enabled.
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
Who is online
Users browsing this forum: No registered users and 1 guest