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

Wind Turbine (VAWT) - Negative Moment in Fluent

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

Like Tree4Likes
  • 1 Post By dalecooper
  • 1 Post By blackmask
  • 1 Post By dalecooper
  • 1 Post By JUDASHU

Reply
 
LinkBack Thread Tools Display Modes
Old   April 30, 2013, 08:04
Default Wind Turbine (VAWT) - Negative Moment in Fluent
  #1
New Member
 
dalecooper
Join Date: Jul 2012
Posts: 10
Rep Power: 4
dalecooper is on a distinguished road
Hello!

I am doing some research about Vertical Axis Wind Turbine (VAWT) for my master thesis. I have been doing some simulations with Fluent with a 4 bladed turbine. The radius of the turbine is 3 meters.

The main target is to obtain a profile of the Coefficient of Perfomance (Cp) versus the Tip Speed Ratio (TSR). In Fluent I'm using a sliding mesh technique with 1 interface to separate the rotating zone from the stationary zone.

The problem I found is that whatever value of TSR I choose I never get a positive moment coefficient. I have tried to refine the mesh and reduce the time step but that didn't solve the problem.

Any idea or solution would be well received. Thanks in advance!
Ansys2015 likes this.
dalecooper is offline   Reply With Quote

Old   April 30, 2013, 08:58
Default
  #2
Senior Member
 
Join Date: Aug 2011
Posts: 315
Rep Power: 11
blackmask will become famous soon enough
You should do the steady simulation with SRF/MRF first. This way you need not warry about the time step size and yet the converged result can be used as a starting point for unsteady sliding-mesh based simulation. If this is a realy turbine, you could ask the designer about the optimal TSR. Do you happen to rotate in the wrong direction?
Ansys2015 likes this.
blackmask is offline   Reply With Quote

Old   April 30, 2013, 10:19
Default
  #3
New Member
 
dalecooper
Join Date: Jul 2012
Posts: 10
Rep Power: 4
dalecooper is on a distinguished road
Quote:
Originally Posted by blackmask View Post
You should do the steady simulation with SRF/MRF first. This way you need not warry about the time step size and yet the converged result can be used as a starting point for unsteady sliding-mesh based simulation. If this is a realy turbine, you could ask the designer about the optimal TSR. Do you happen to rotate in the wrong direction?
Hello blackmask.

I have already used this approach of beginning with a steady simulation with an angular velocity equal to 0 rad/s as my thesis director suggested to me. After that and when the results have converged I move to transient simulation using a constant inlet velocity of 10 m/s and changing the value of the rotational speed in order to simulate different values of the TSR.

I'm using 'Transition SST (4 eqn)' with the default options as a viscous model and I'm not sure if that could be the problem. I have also tried with different solution methods but apart from the covergence velocity no significant changes were detected in Cm.
dalecooper is offline   Reply With Quote

Old   April 30, 2013, 10:32
Default
  #4
Senior Member
 
Join Date: Aug 2011
Posts: 315
Rep Power: 11
blackmask will become famous soon enough
No. You should start your simulation with 10m/s inlet velocity and an appropriate rotational speed using the MRF (moving reference frame) technique. You should check the fluent user's guide if you are not familiar with MRF. The default parameters of specifying b.c. for TR-SST could results in an unresonable high turbulent viscosity at the inlet but I do not think it is the reason that you get negtive torque.
blackmask is offline   Reply With Quote

Old   May 4, 2013, 01:57
Default
  #5
New Member
 
Reza
Join Date: Sep 2010
Posts: 9
Rep Power: 6
CFD_LES is on a distinguished road
Hello dalecooper

I work on a 4-Bladed VAWT with semi-circle blade. I obtain ((Negative)) moment for this VAWT.
You obtained negative moment too. Did you understand the reason?
Could you solve your problem?
Also I have to say, My problem is steady.
Thanks so much
CFD_LES is offline   Reply With Quote

Old   July 8, 2013, 06:58
Default
  #6
New Member
 
dalecooper
Join Date: Jul 2012
Posts: 10
Rep Power: 4
dalecooper is on a distinguished road
Quote:
Originally Posted by CFD_LES View Post
Hello dalecooper

I work on a 4-Bladed VAWT with semi-circle blade. I obtain ((Negative)) moment for this VAWT.
You obtained negative moment too. Did you understand the reason?
Could you solve your problem?
Also I have to say, My problem is steady.
Thanks so much
Hello CFD_LES,

Finally I was able to solve my problem of negative Cm values in my simulations. All I did was to restart again with the geometry and meshing. So I think there was some problem with the interfaces between the stationary and rotating domains, but I'm not sure if that was the reason. All I can tell you is that after that I began to obtain positive Cm values.

Nowadays I'm doing some simulations in order to obtain the coefficient of performance of the turbine as a function of the tip speed ratio (TSR). The results look pretty good except, for a unique point (TSR=1.5) where the Cm has an abnormal value. I was expecting to obtain an increasing Cm value up to a certain point where the turbine is most efficient, and after that the Cm should begin to decrease.

I have tried with different inlet velocities, but the behavior is exactly the same, and also the critical point is located at TSR=1.5

Best regards
israa likes this.
dalecooper is offline   Reply With Quote

Old   October 8, 2013, 03:53
Default
  #7
New Member
 
Michael Bausas
Join Date: May 2013
Posts: 10
Rep Power: 4
mikebausas is on a distinguished road
hello dale cooper,

I am working with a 2d model VAWT which i simulated using fluent. And now, I generated some results of Cm from the iterations. My question is, what would be my next step in order for me to conclude that the results are acceptable or accurate for the specific turbine I am testing? Where shall I compare my results with? And how did you compute for the other parameters, such as Cp, Torqu, etc? Does Fluent also gives these data? Thanks

Regards,
Mike
mikebausas is offline   Reply With Quote

Old   October 8, 2013, 05:06
Default
  #8
New Member
 
dalecooper
Join Date: Jul 2012
Posts: 10
Rep Power: 4
dalecooper is on a distinguished road
Quote:
Originally Posted by mikebausas View Post
hello dale cooper,

I am working with a 2d model VAWT which i simulated using fluent. And now, I generated some results of Cm from the iterations. My question is, what would be my next step in order for me to conclude that the results are acceptable or accurate for the specific turbine I am testing? Where shall I compare my results with? And how did you compute for the other parameters, such as Cp, Torqu, etc? Does Fluent also gives these data? Thanks

Regards,
Mike
Hello mikebausas,

In turbines, one of the most important graphic to obtain is the one that compares the Cp vs TSR (Tip Speed Ratio). From Fluent you can obtain directly the Cm (Coefficient of Moment). Once you get the Cm you must average the values obtained through a number of periods assuring that you have reached a steady solution and use the following equation to obtain the coefficient of performance (Cp):

Cp=Cm*TSR

What I have been doing is simulate for different values of TSR. To do that you have to change the value of the rotational velocity of the sliding mesh zone in which you locate the airfoils for every TSR value.

Regards.
dalecooper is offline   Reply With Quote

Old   April 3, 2014, 04:42
Default
  #9
New Member
 
Ali
Join Date: Dec 2012
Location: Turkey
Posts: 10
Rep Power: 4
mct90 is on a distinguished road
Quote:
Originally Posted by dalecooper View Post
Hello mikebausas,

In turbines, one of the most important graphic to obtain is the one that compares the Cp vs TSR (Tip Speed Ratio). From Fluent you can obtain directly the Cm (Coefficient of Moment). Once you get the Cm you must average the values obtained through a number of periods assuring that you have reached a steady solution and use the following equation to obtain the coefficient of performance (Cp):

Cp=Cm*TSR

What I have been doing is simulate for different values of TSR. To do that you have to change the value of the rotational velocity of the sliding mesh zone in which you locate the airfoils for every TSR value.

Regards.
Hi Dalecooper. i am also studying similar project.
I want to know, to obtain a CP-TSR graph like the image, does the cm decrease at a value of rotational velocity(w).

For ex; w=1,tsr=1 cm=0.12 ---w=2,tsr=2 cm=0.20---- w=3 trsr=3 cm=0.18 ---.
cp.JPG
mct90 is offline   Reply With Quote

Old   April 6, 2014, 22:56
Default
  #10
New Member
 
behzad shahizare
Join Date: Feb 2014
Location: malaysia
Posts: 5
Rep Power: 3
shahizare is on a distinguished road
Send a message via Skype™ to shahizare
Quote:
Originally Posted by mct90 View Post
Hi Dalecooper. i am also studying similar project.
I want to know, to obtain a CP-TSR graph like the image, does the cm decrease at a value of rotational velocity(w).

For ex; w=1,tsr=1 cm=0.12 ---w=2,tsr=2 cm=0.20---- w=3 trsr=3 cm=0.18 ---.
Attachment 29809

Hi Dalecooper
kindly please help me . i simulate vawt in fluent in 2d and 3d model . i have problem with a reference numbers in fluent , what do you think how much i should put for area and length . and after getting results (cm) how can i convert it to Cp?
shahizare is offline   Reply With Quote

Old   May 5, 2014, 14:03
Default
  #11
New Member
 
shafqat
Join Date: Apr 2014
Posts: 8
Rep Power: 3
shafqat is on a distinguished road
hi i just want to know why my Cm keeps on increasing at high tsr..i mean it has to decrease after optimal tsr
shafqat is offline   Reply With Quote

Old   May 14, 2014, 06:32
Default
  #12
New Member
 
Ali
Join Date: Dec 2012
Location: Turkey
Posts: 10
Rep Power: 4
mct90 is on a distinguished road
Quote:
Originally Posted by shahizare View Post
Hi Dalecooper
kindly please help me . i simulate vawt in fluent in 2d and 3d model . i have problem with a reference numbers in fluent , what do you think how much i should put for area and length . and after getting results (cm) how can i convert it to Cp?
Hi,

I am working 2d vawt. In reference values i took length=R (rotor radius) and Area= 2*R*h .(h=1 because of 2D). references are below;

http://www.iaa.ncku.edu.tw/~aeromems...%20control.pdf (page 341 )

Impacts of solidity and hybrid system in small wind turbines
performance M.H. Mohamed (page 497)
mct90 is offline   Reply With Quote

Old   May 15, 2014, 02:46
Default
  #13
New Member
 
behzad shahizare
Join Date: Feb 2014
Location: malaysia
Posts: 5
Rep Power: 3
shahizare is on a distinguished road
Send a message via Skype™ to shahizare
Quote:
Originally Posted by mct90 View Post
Hi,

I am working 2d vawt. In reference values i took length=R (rotor radius) and Area= 2*R*h .(h=1 because of 2D). references are below;

http://www.iaa.ncku.edu.tw/~aeromems...%20control.pdf (page 341 )

Impacts of solidity and hybrid system in small wind turbines
performance M.H. Mohamed (page 497)


abi nasilsin


this is link of my work https://www.dropbox.com/sh/c2hl306cm...Tm6mI8i8xX2T8a

actually in what reason i dont know why i can not get good validation with my reference . is it possible for you to check my case file?

tasakor lar
shahizare is offline   Reply With Quote

Old   May 30, 2014, 00:05
Default
  #14
New Member
 
Aalisha
Join Date: Jan 2014
Posts: 18
Rep Power: 3
aalisha is on a distinguished road
Hello CFD Mates!!

I am a Masters student working on a vertical axis marine current turbine. I am using Fluent for the required simulations. To calculate the Cp, I have monitored the value of Cm, Cl and Cd.

However, the value obtained for Cm is very high ( 85, 70 etc). If anyone can guide me what is the reason for this. It will be a great help. Thanks.
aalisha is offline   Reply With Quote

Old   August 7, 2014, 10:19
Default
  #15
New Member
 
amin veysi
Join Date: May 2013
Posts: 5
Rep Power: 4
amin_veysi is on a distinguished road
Quote:
Originally Posted by aalisha View Post
Hello CFD Mates!!

I am a Masters student working on a vertical axis marine current turbine. I am using Fluent for the required simulations. To calculate the Cp, I have monitored the value of Cm, Cl and Cd.

However, the value obtained for Cm is very high ( 85, 70 etc). If anyone can guide me what is the reason for this. It will be a great help. Thanks.
hi aalisha change the reference value,
reference value panel, change the velocity, area and length
amin_veysi is offline   Reply With Quote

Old   August 11, 2014, 22:56
Default
  #16
New Member
 
Aalisha
Join Date: Jan 2014
Posts: 18
Rep Power: 3
aalisha is on a distinguished road
Hey Amin,

Thanks for your kind reply. I tried it and its working now.
aalisha is offline   Reply With Quote

Old   January 20, 2015, 22:30
Default
  #17
New Member
 
WenBin.JU
Join Date: Jul 2014
Posts: 1
Rep Power: 0
JUDASHU is on a distinguished road
Quote:
Originally Posted by mct90 View Post
Hi,

I am working 2d vawt. In reference values i took length=R (rotor radius) and Area= 2*R*h .(h=1 because of 2D). references are below;

http://www.iaa.ncku.edu.tw/~aeromems...%20control.pdf (page 341 )

Impacts of solidity and hybrid system in small wind turbines
performance M.H. Mohamed (page 497)
hi
I am doing a 2D VAWT CFD simulation with Fluent, and I have my reference valus as you suggested, I still wonder what the "compute from" and "reference area" should be, and I have got a CM curve in my attachment, it seems fluctuates drastic, so is this a normal case or did I have something wrong in my case?
Ansys2015 likes this.
JUDASHU 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
Wind Turbine on FLUENT federvo.mala ANSYS 21 July 3, 2013 02:05
non-orthogonal faces and incorrect orientation? nennbs OpenFOAM Native Meshers: blockMesh 7 April 17, 2013 05:42
Moving reference frame in wind turbine kongl1986 FLUENT 0 March 30, 2013 11:50
2D Simulation of Savonius Wind Turbine ravindersingh FLUENT 4 December 9, 2011 14:00
Boundary Condition for ducted wind turbine using FLUENT Styward FLOW-3D 2 October 30, 2009 01:53


All times are GMT -4. The time now is 22:25.