CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > CFX

shaft power versus RMP curve in wind turbine simulation

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   August 25, 2012, 14:04
Default shaft power versus RMP curve in wind turbine simulation
  #1
New Member
 
michael
Join Date: Aug 2012
Posts: 13
Rep Power: 4
zhenduan306 is on a distinguished road
Hi, Dear friends,
I am simulating HAWT now. My purpose is to obtain the curve between shaft power and RMP(or wind speed , TSR).

I used the way ghorrocks mentioned in
Wind Turbine - Rotation of the Baldes

"Far easier, but a little less elegant is to sweep through a range of rotation speeds in a series of steady state runs. This will give you a torque/power versus speed curve. Then interpolate on this curve your intended torque/power and there you have it. If you like then run this point as well to improve the accuracy of the interpolation"

The weird thing is, when I fixed wind velocity, I increased rotational speed and torque increase too. The result is Cp is more than 0.59 (betz limit).
For example, wind speed=8m/s, angular velocity are 20 rad/s, 40 rad/s and 50 rad/s, respectively. The corresponding torque are 8.5, 22.35, 34.8 N.m. I think it is wrong. Because when wind speed is constant, the total wind power in the rot domain of the blades is fixed. Higher rotating speed means lower torque.



(my analysis type:steady)

Hence, I want to add CEL to my cfx-pre. This link gives me some clues:
Difficulty in calculating angular velocity of Savonius turbine simulation
from his CEL: he defined wind power as an expression. But I don't understand where he used it.

in my mind, I think the wind power and Betz limit should be as the constraint conditions. Otherwise, the results will be wrong.

These are my CFX pictures:



Any help would be greatly appreciated.
zhenduan306 is offline   Reply With Quote

Old   August 25, 2012, 20:22
Default
  #2
Member
 
Paulo
Join Date: Jun 2011
Posts: 31
Rep Power: 5
strobel is on a distinguished road
I'm having this problem aswell. I'm using this expression to cp:
P_dot= torque_z()@blades*2*pi*Angular velocity [rev s^-1]
P= 0.5*ave(Density)@inlet*ave(Velocity)@inlet^3*area( )@in_rotor
cp=P_dot/P

But it results in cp>1.
strobel is offline   Reply With Quote

Old   August 25, 2012, 20:56
Default
  #3
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,469
Rep Power: 76
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Before worrying about getting the speed/torque correct, have you done all the basic checks to ensure your simulation is accurate?

http://www.cfd-online.com/Wiki/Ansys..._inaccurate.3F
ghorrocks is offline   Reply With Quote

Old   August 25, 2012, 21:38
Default
  #4
New Member
 
michael
Join Date: Aug 2012
Posts: 13
Rep Power: 4
zhenduan306 is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
Before worrying about getting the speed/torque correct, have you done all the basic checks to ensure your simulation is accurate?

http://www.cfd-online.com/Wiki/Ansys..._inaccurate.3F
Hi, ghorrocks, I think I can't include all relevant physics conditions in my simulation. I wonder I should set wind power in the rot domain as a constrain condition. But I don't know where can I add it and how to use it.

By the way, my blades' TSR=7.2. According to the rated wind speed, the corresponding rotating speed can be obtained. In that case, I can obtain the reasonable shaft power and Cp.
zhenduan306 is offline   Reply With Quote

Old   August 26, 2012, 07:11
Default
  #5
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,469
Rep Power: 76
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Quote:
I wonder I should set wind power in the rot domain as a constrain condition.
I have no idea what you mean by that.

The method I recommend is quite simple - do a series of steady state (well, frozen rotor or other simple RFR model) at various rotation speeds and get the torque versus rotation speed line. You can then match that to the load point (or line) to get the operating point.

I am not a wind turbine person so I have no idea what TSR=7.2 means. But you seem to imply that you already know the operating speed of this turbine.
ghorrocks is offline   Reply With Quote

Old   August 26, 2012, 08:27
Default
  #6
New Member
 
michael
Join Date: Aug 2012
Posts: 13
Rep Power: 4
zhenduan306 is on a distinguished road
Hi, ghorrocks,Thank you for your reply. Actually, I know all the information of the blades in my wind turbine. So, I just use it to verify my simulation.

In fact, I found some guys have the same problem:
how to find wind turbine torque from CFD ?

Problem with torque output from a turbine simulation.

Anyway, I will try to figure it out. After that, I will respond my method.
zhenduan306 is offline   Reply With Quote

Old   August 26, 2012, 08:29
Default
  #7
New Member
 
michael
Join Date: Aug 2012
Posts: 13
Rep Power: 4
zhenduan306 is on a distinguished road
Quote:
Originally Posted by strobel View Post
I'm having this problem aswell. I'm using this expression to cp:
P_dot= torque_z()@blades*2*pi*Angular velocity [rev s^-1]
P= 0.5*ave(Density)@inlet*ave(Velocity)@inlet^3*area( )@in_rotor
cp=P_dot/P

But it results in cp>1.
Hi, strobel, did you resolve it?
zhenduan306 is offline   Reply With Quote

Old   August 26, 2012, 20:17
Default
  #8
Member
 
Paulo
Join Date: Jun 2011
Posts: 31
Rep Power: 5
strobel is on a distinguished road
Quote:
Originally Posted by zhenduan306 View Post
Hi, strobel, did you resolve it?
No, my problem is quite strange, because when i set a lower rotation i get cp<0.59, but if i increase the rotation, cp increase too and i get cp >1. I gonna try that way Glenn said, if it fix i post here the results.
strobel is offline   Reply With Quote

Old   August 26, 2012, 21:00
Default
  #9
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,469
Rep Power: 76
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
If you already know the operating point, then you not simulate that point, and check the torque it produces is correct? No need to generate an operating curve if you do not need it.
ghorrocks is offline   Reply With Quote

Old   August 26, 2012, 21:24
Default
  #10
New Member
 
michael
Join Date: Aug 2012
Posts: 13
Rep Power: 4
zhenduan306 is on a distinguished road
Quote:
Originally Posted by strobel View Post
No, my problem is quite strange, because when i set a lower rotation i get cp<0.59, but if i increase the rotation, cp increase too and i get cp >1. I gonna try that way Glenn said, if it fix i post here the results.
Good luck to you!

Last edited by zhenduan306; August 26, 2012 at 22:19.
zhenduan306 is offline   Reply With Quote

Old   August 26, 2012, 21:27
Default
  #11
New Member
 
michael
Join Date: Aug 2012
Posts: 13
Rep Power: 4
zhenduan306 is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
If you already know the operating point, then you not simulate that point, and check the torque it produces is correct? No need to generate an operating curve if you do not need it.
I just know the information of the blades. I am designing a new wind turbine now. Hence, I need this curve to choose a generator.
zhenduan306 is offline   Reply With Quote

Old   August 27, 2012, 06:24
Default
  #12
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,469
Rep Power: 76
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Quote:
Hence, I need this curve to choose a generator.
Then you need to generate the turbine operating curve.
ghorrocks is offline   Reply With Quote

Old   August 27, 2012, 07:18
Default
  #13
D.B
Member
 
DB
Join Date: Apr 2011
Posts: 60
Rep Power: 5
D.B is on a distinguished road
Try areaAve or VolAve and see if it works out
__________________
-D.B
D.B is offline   Reply With Quote

Old   August 27, 2012, 08:02
Default
  #14
New Member
 
michael
Join Date: Aug 2012
Posts: 13
Rep Power: 4
zhenduan306 is on a distinguished road
Quote:
Originally Posted by D.B View Post
Try areaAve or VolAve and see if it works out
Hi DB, can you tell me more details?

Thank you very much!
zhenduan306 is offline   Reply With Quote

Old   August 28, 2012, 01:23
Default
  #15
Member
 
Paulo
Join Date: Jun 2011
Posts: 31
Rep Power: 5
strobel is on a distinguished road
Zhenduan, maybe when you increase too much your rotation your turbine don't work freely in this rotation and doesn't work anymore as Turbine. So, in these conditions, it is no longer more a turbine and become a pump.
strobel is offline   Reply With Quote

Old   August 28, 2012, 07:25
Default
  #16
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,469
Rep Power: 76
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
I presume you are looking for the steady state operating condition - in that case, if the net torque is negative it simply means the rotor will decelerate, so the steady state speed is slower than the current speed.
ghorrocks is offline   Reply With Quote

Old   August 28, 2012, 09:44
Default
  #17
New Member
 
michael
Join Date: Aug 2012
Posts: 13
Rep Power: 4
zhenduan306 is on a distinguished road
Quote:
Originally Posted by strobel View Post
Zhenduan, maybe when you increase too much your rotation your turbine don't work freely in this rotation and doesn't work anymore as Turbine. So, in these conditions, it is no longer more a turbine and become a pump.
Yes, you are right! CFX is just a tool. it can't identify the minimum or maximum condition the wind turbine is able to run.
zhenduan306 is offline   Reply With Quote

Old   August 28, 2012, 09:48
Default
  #18
New Member
 
michael
Join Date: Aug 2012
Posts: 13
Rep Power: 4
zhenduan306 is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
I presume you are looking for the steady state operating condition - in that case, if the net torque is negative it simply means the rotor will decelerate, so the steady state speed is slower than the current speed.
OK, I will check it. Thank you very much!
zhenduan306 is offline   Reply With Quote

Old   October 31, 2012, 16:33
Default
  #19
Member
 
Paulo
Join Date: Jun 2011
Posts: 31
Rep Power: 5
strobel is on a distinguished road
Quote:
Originally Posted by zhenduan306 View Post
OK, I will check it. Thank you very much!
did you solve your problem?
i still getting high values of coefficient of shaft power(Cp). Have been strange, because in low rotation, i get Cp near of experimental values, but when i increase rotation Cp increase so fas and i get Cp higher than Beltz limit.
strobel is offline   Reply With Quote

Old   November 1, 2012, 15:24
Default
  #20
New Member
 
michael
Join Date: Aug 2012
Posts: 13
Rep Power: 4
zhenduan306 is on a distinguished road
Quote:
Originally Posted by strobel View Post
did you solve your problem?
i still getting high values of coefficient of shaft power(Cp). Have been strange, because in low rotation, i get Cp near of experimental values, but when i increase rotation Cp increase so fas and i get Cp higher than Beltz limit.
Unfortunately, I don't have time to continue it. If you just want to get Cp and Tip Speed Ratio curve, I suggest you use Qblade. It is a very good software. You can download from websit. Just google it.
zhenduan306 is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
FSI - Wind Turbine AUN CFX 13 August 29, 2012 16:44
Vertical axis wind turbine simulation vincentwks ANSYS 0 April 10, 2012 04:27
CFD analysis on wind turbine rotor Ken (Wind Turbine CFD Super Rookie) Main CFD Forum 42 July 18, 2011 22:11
1D analysis and simulation of gas turbine power plant sridhar.d009 Main CFD Forum 1 September 28, 2010 07:38
Pressure Drop/Porous Media within Shrouded Wind Turbine Simulation joshdsouza FLUENT 0 March 3, 2010 18:02


All times are GMT -4. The time now is 04:23.