CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Post-Processing

calculate epsilon from RANS k-Omega simulation

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

Like Tree3Likes
  • 1 Post By CFD_Monkey
  • 2 Post By CFD_Monkey

LinkBack Thread Tools Display Modes
Old   June 20, 2013, 06:07
Default calculate epsilon from RANS k-Omega simulation
New Member
Join Date: Dec 2012
Posts: 7
Rep Power: 6
CFD_Monkey is on a distinguished road
Hi FOAMers,

in order to get the epsilon value from a RANS simulation (kOmega, kOmegaSST respectively) I was searching the source code. This is what I've found in kOmegaSST.H:

00250         //- Return the turbulence kinetic energy dissipation rate
 00251         virtual tmp<volScalarField> epsilon() const
 00252         {
 00253             return tmp<volScalarField>
 00254             (
 00255                 new volScalarField
 00256                 (
 00257                     IOobject
 00258                     (
 00259                         "epsilon",
 00260                         mesh_.time().timeName(),
 00261                         mesh_
 00262                     ),
 00263                     betaStar_*k_*omega_,
 00264                     omega_.boundaryField().types()
 00265                 )
 00266             );
 00267         }
My question: How can I run this section to get epsilon (for e.g. in case of the reynoldsstresstensor it's just the command R)

immortality likes this.
CFD_Monkey is offline   Reply With Quote

Old   June 27, 2013, 10:55
New Member
Join Date: Dec 2012
Posts: 7
Rep Power: 6
CFD_Monkey is on a distinguished road
I don't know if this is the best way to get the epsilon value but it works!

1. create the directory ~/Openfoam/$USER-2.1.x/utilities/postProcessing/turbulence/

2. copy the attached files into.

3. run the command "wmake" from the epsilon directory

This will create the folders platforms/.../bin in the ~/Openfoam/$USER-2.1.x/ directory and compiles the epsilon.C file. The executive "epsilon" is created there and can be called by "epsilon" from any case directory.

Note: This should work for kOmega as well as for kOmegaSST. While the difference is simply a change of constants:
kOmega: epsilon=Cmu*k*omega
kOmegaSST: epsilon=betaStar*k*omega
Which is normally the same value: betaStar=Cmu=0.09
Attached Files
File Type: zip (94.6 KB, 57 views)
McCharles and weizl08 like this.
CFD_Monkey is offline   Reply With Quote


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
How to calculate torque from the simulation result sivakumar OpenFOAM Post-Processing 18 January 6, 2017 04:00
Unsteady RANS simulation of Jet in crossflow harshad88 FLUENT 0 June 2, 2013 13:29
epsilon and K blowing up. sivakumar OpenFOAM Running, Solving & CFD 1 October 25, 2012 04:50
Do RANS simulation overpredict maximum Cp? Sergio Rossi Main CFD Forum 0 May 15, 2007 09:23
RANS simulation prasat Main CFD Forum 0 June 4, 2003 11:04

All times are GMT -4. The time now is 20:19.