CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Post-Processing

How to calculate power and thrust using an actuator disk model for turbines?

Register Blogs Community New Posts Updated Threads Search

Like Tree3Likes
  • 2 Post By tareqkh
  • 1 Post By pbachant

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 31, 2017, 01:10
Default How to calculate power and thrust using an actuator disk model for turbines?
  #1
New Member
 
Josh Bowman
Join Date: Oct 2015
Location: US
Posts: 3
Rep Power: 11
Bowman is on a distinguished road
Hi,
I am a student and I have been working with OpenFOAM for a few months. I have seen various posts about turbine simulations within OpenFOAM here on CFD-online. I am currently using the actuationDisksource within fvOptions in OpenFOAM for a horizontal axis hydrokinetic turbines. My plan of study is to test each type of model (actuationDisksource, radialActuationDisksource, rotorDisksource, possibly an actuator line model, and finally a fully resolved turbine). I am just starting and I have noticed that within OpenFOAM there are the force and forceCoeffs functions. I am of course, interested in the power and thrust on the Disk(the swept area of the turbine blade). What is considered the typical/standard process among users to calculate these performance parameters? Is this done through a matlab or python script reading in the force output or is there a calculation feature with OpenFOAM?

The versions of OF that I have access to at my facility are 2.3 and v1606+. All grid generation is done through Pointwise and post-processing is done through either ParaView or Ensight.

Thank you and any help would be greatly appreciated,
Josh
Bowman is offline   Reply With Quote

Old   January 31, 2017, 11:30
Default
  #2
Senior Member
 
tareqkh's Avatar
 
CFD
Join Date: Nov 2010
Location: United States
Posts: 243
Rep Power: 17
tareqkh is on a distinguished road
Hi Josh,

I have attached a file where you will be able to extract the volume forces. I am assuming that you are using the actuator disk model and you have a file named volume VolumeForce in the 0 folder! You need to copy the file to the system folder. Perform the simulation once it finishes type sample in the terminal. Finally, you need to adjust the coordinates of the sampleDict to the coordinates of your turbine in the mesh i.e. the radius of your turbine and the location.

Also, you can integrate the pressure over the whole disk to get the total thrust. You could extract the pressure over the disk then use Matlab (trapz) to perform the integration.

Best,
Attached Files
File Type: zip sampleDict.zip (556 Bytes, 116 views)
Bowman and Simurgh like this.
tareqkh is offline   Reply With Quote

Old   January 31, 2017, 17:22
Default
  #3
New Member
 
Josh Bowman
Join Date: Oct 2015
Location: US
Posts: 3
Rep Power: 11
Bowman is on a distinguished road
Thank you Khamlaj!
I do not have a file named VolumeForce in the /0 folder. That is something I'll need to do. Could you point me in the right direction to see how this is implemented? I use Pointwise for meshing. So the disk volume is set to cellZone. I am assuming that within the file VolumeForce, I will specify my cellZone for disk? I'll upload my case when I get back to the HPCC. However, I am essentially just modifying the turbineSiting tutorial in simpleFoam to fit my mesh.

Here is what I am trying to simulate. An actuator disk model with a hub/support structure. As of right now, I only have the actuator disk model without the hub/support structure but it does have a hole through it.

U_magnitude_Iso_surface_and_vector_clip.jpg

AD_mesh_with_hole.jpg

AD_mesh_support_model.jpg
Bowman is offline   Reply With Quote

Old   January 31, 2017, 19:45
Default
  #4
Senior Member
 
Pete Bachant
Join Date: Jun 2012
Location: Boston, MA
Posts: 173
Rep Power: 14
pbachant is on a distinguished road
I see two options here:

1. Use an instance of a probes functionObject at the actuator disk's upstreamPoint, which can then be used to calculate power with 1/2*rho*A*U**3.
2. Modify actuationDiskSourceTemplates.C to print the upstream velocity at each iteration. At line 57, you could add

Code:
Info<< "Upstream velocity for " << name_ << ": " << upU << endl;
Bowman likes this.
__________________
Home | Twitter | GitHub
pbachant is offline   Reply With Quote

Old   October 27, 2017, 04:11
Default
  #5
New Member
 
Maria
Join Date: Feb 2017
Posts: 25
Rep Power: 9
hoemmaria is on a distinguished road
Hey Josh!

Did you find a good way to calculate the power from an actuator disk? I am doing the same for my master thesis. However, I do not have a constant loaded disk, but more of a gaussian load distribution.

Kind regards,
Maria
hoemmaria is offline   Reply With Quote

Old   March 15, 2021, 11:45
Default
  #6
New Member
 
Kim Best
Join Date: Jan 2021
Posts: 20
Rep Power: 5
KimBest is on a distinguished road
Hi,


I am having to do a similar thing. Would I be right in saying that the purpose of upstreamPoint is to record a value of velocity at a specified coordinate, for the purpose of calculating power?


If so, does this mean pasting the code posted above:



Info<< "Upstream velocity for " << name_ << ": " << upU << endl;


in the aforementioned file will output this velocity to be used in the calculation?


Thanks,


Kim
KimBest is offline   Reply With Quote

Reply

Tags
actuator, calculate, forces and force coeff., performance, turbine


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 Off
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 12:53.