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

calculation of k, epsilon and omega

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

Like Tree52Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 28, 2010, 10:50
Default calculation of k, epsilon and omega
  #1
New Member
 
Philipp Bachmann
Join Date: May 2010
Location: Esslingen, Germany
Posts: 7
Rep Power: 16
kroetenechse is on a distinguished road
Hello, i´m a new OpenFOAM user.

Could somebody give me a advice, how I can calculate k (turbulent kinematic energy) and epsilon (dissipation rate) on every boundary for my k-epsilon-model.

And the calculation of epsilon, k and omega for my k-w-SST-Model. Do I need also nuTilda for the k-w-SST-Model?

If someone can help me, I would be very happy.

Thank you in advance.
phil
Gang Wang, labyrinth01 and erinsam like this.
kroetenechse is offline   Reply With Quote

Old   October 28, 2010, 12:02
Default
  #2
Member
 
Simon Lapointe
Join Date: May 2009
Location: Québec, Qc, Canada
Posts: 33
Rep Power: 17
Simon Lapointe is on a distinguished road
Hi,

At the inlet, the values of k and epsilon can be calculated using the level of turbulence intensity and viscosity ratio (or turbulent length scale) you want to achieve. The relations found on this page can be useful:
http://www.cfd-online.com/Wiki/Turbu...ary_conditions

On wall boundary conditions, you use kqRWallFunction and epsilonWallFunction since the k-epsilon model makes use of wall functions.

For the k-w SST model, you need to specify BCs for k and omega, not epsilon nor nuTilda.

Hope this helps.
Simon Lapointe is offline   Reply With Quote

Old   November 4, 2010, 14:50
Default
  #3
New Member
 
Philipp Bachmann
Join Date: May 2010
Location: Esslingen, Germany
Posts: 7
Rep Power: 16
kroetenechse is on a distinguished road
Hi Simon,

thank you very much! Now I could calculate epsilon, k and omega. This was one of my main problems. You don´t know, how you helped me with that information :-). Also for these three variables I found some functions:
turbulentIntensityKineticEnergyInlet for k in the inlet (under that you make the description of the intensity and the value)
turbulentMixingLengthFrequencyInlet for omega in the inlet (then description mixing length, k and value)
turbulentMixingLengthDissipationRateInlet for epsilon in the inlet (then description mixing length and value)
...like that for example k

inlet
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.02;
value uniform 0.0006;
}

But relating to the wallfunctions I need a little bit more support, please. I´m really sorry.

- For the k-epsilon-model I can use kqRWallFunction for k and for epsilon I can use
epsilonWallFunction?
like that?
[k]
wall
{
type kqRWallFunction;
value uniform ...;
}

- Are there more wallfunctions for the k-omega-SST model?

- what can I use for the outlet?

- what means type zerogradient? The variabel will not change in that boundary?

Thank you in Advance
phil
kroetenechse is offline   Reply With Quote

Old   November 9, 2010, 09:33
Default
  #4
New Member
 
Join Date: Nov 2010
Posts: 6
Rep Power: 15
Gerard is on a distinguished road
Hi Philipp,

as far as I know you can use kqRWallFunction for k, epsilonWallFunction for epsilon and omegaWallFunction for omega. Furthermore there are Wall Functions for nut (implying that you want to do incompressible calculations), such as nutWallFunction or nutSpalartAllmarasWallFunction. Have a look at the code in your Openfoam directory /src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions for details and further wall functions.
If you want to use kOmegaSST it is clear that you need to specify omega instead of epsilon.
I hope this helps a little bit. I am afraid, but I am new both to CFD in general and to OpenFoam in particular.


Hi everybody,

as I just mentioned, I'm new to CFD and OpenFoam. So I hope it is not to annoying to ask you a simple question:
the value specified for a wall function, should it be the value at the wall (k=0,...) or at the first grid point let's say at 40 y+?

Thanks.
Gerard
mm.abdollahzadeh likes this.
Gerard is offline   Reply With Quote

Old   November 9, 2010, 10:10
Default
  #5
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20
vkrastev is on a distinguished road
Quote:
Originally Posted by Gerard View Post
Hi Philipp,

as far as I know you can use kqRWallFunction for k, epsilonWallFunction for epsilon and omegaWallFunction for omega. Furthermore there are Wall Functions for nut (implying that you want to do incompressible calculations), such as nutWallFunction or nutSpalartAllmarasWallFunction. Have a look at the code in your Openfoam directory /src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions for details and further wall functions.
If you want to use kOmegaSST it is clear that you need to specify omega instead of epsilon.
I hope this helps a little bit. I am afraid, but I am new both to CFD in general and to OpenFoam in particular.


Hi everybody,

as I just mentioned, I'm new to CFD and OpenFoam. So I hope it is not to annoying to ask you a simple question:
the value specified for a wall function, should it be the value at the wall (k=0,...) or at the first grid point let's say at 40 y+?

Thanks.
Gerard
The value you specify for a wall function is only an initial guess (as it is the initial value for the internal field), so theoretically speaking it could be quite arbitrary. However, you can (of course) "accomodate" the convergence of the solver by choosing a more "reasonable" value: one common practice (you can see it also in some of the tutorials, for instance the pitzDaily incompressible steady-state case) is to estimate the initial turbulent quantities for the internal field and then use the same values as initial guesses for the wall function entry.
Hope this helps

V.
amin_jalalian, serles, aow and 6 others like this.
vkrastev is offline   Reply With Quote

Old   November 9, 2010, 10:41
Default
  #6
New Member
 
Join Date: Nov 2010
Posts: 6
Rep Power: 15
Gerard is on a distinguished road
Hi Vesselin,

thank you for your reply!
Do I get you right, that the values specified for the wall functions are initial values and not "classical boundary conditions", that need to be fulfilled within every time step?
So theoretically it would be desirable to set it to the estimated value of the first grid point next to the wall? At the moment I just define uniform values for the internal field
Thanks again.

Gerard
Gerard is offline   Reply With Quote

Old   November 9, 2010, 10:53
Default
  #7
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20
vkrastev is on a distinguished road
Quote:
Originally Posted by Gerard View Post
Do I get you right, that the values specified for the wall functions are initial values and not "classical boundary conditions", that need to be fulfilled within every time step?
Yes. If you have a look in the source code, you will see that the wall function for k simply impose a zeroGradient condition at the wall, while the while function for epsilon calculate (for each timestep) the first grid point value by using an algebraic expression derived from the classical logarithmic law-of the wall approach.

Quote:
Originally Posted by Gerard View Post
So theoretically it would be desirable to set it to the estimated value of the first grid point next to the wall?
Yes, but if you use a reasonable assumption for the internal field initial value (for instance, estimating a turbulence intensity and a turbulent lenght scale for the bulk or free-stream flow), then it should do a good job as well.

Best regards

V.
takahisa, aow and SAJESH P like this.
vkrastev is offline   Reply With Quote

Old   November 9, 2010, 11:23
Default
  #8
New Member
 
Join Date: Nov 2010
Posts: 6
Rep Power: 15
Gerard is on a distinguished road
Thank you very much! That really helped me.

Best regards
Gerard
Gerard is offline   Reply With Quote

Old   November 25, 2010, 11:02
Default
  #9
New Member
 
Philipp Bachmann
Join Date: May 2010
Location: Esslingen, Germany
Posts: 7
Rep Power: 16
kroetenechse is on a distinguished road
Thanks to all,

my calculation did run now. You were a great help.

Thanks a lot.

phil
kroetenechse is offline   Reply With Quote

Old   January 17, 2011, 17:51
Default
  #10
Member
 
Vishal Jambhekar
Join Date: Mar 2009
Location: University Stuttgart, Stuttgart Germany
Posts: 90
Blog Entries: 1
Rep Power: 17
vishal is on a distinguished road
Hi,

I have some same problem please have a look at following post...



http://www.cfd-online.com/Forums/ope...urbulence.html
__________________
Cheers,

Vishal Jambhekar...
"Simulate the way ahead......!!!"
vishal is offline   Reply With Quote

Old   March 9, 2011, 12:52
Default
  #11
jms
Member
 
José
Join Date: Jan 2011
Posts: 73
Rep Power: 15
jms is on a distinguished road
Dear all,

while solving one case of flow over an airfoil with the turbulence model k-omega SST I would like to compute the values of omega at the wall. In OpenFOAM it seems that you can only call it as a "fixedValue", not as "calculated". I tried to call the wall as a "patch" instead of a "wall" and it does not work either. I get the following error:

gradientInternalCoeffs cannot be called for a calculatedFvPatchField
on patch blabla of field omega in file "/.../omega"
You are probably trying to solve for a field with a default boundary condition.

The reason why I want to set it as "calculated" is because I have run the same case in ANSYS CFX and I saw that it gave different values for this omega at the wall. Thus, by fixing it, maybe it is affecting to the results obtained.

Thank you very much for your help,

Regards,

José
jms is offline   Reply With Quote

Old   March 9, 2011, 17:19
Default
  #12
New Member
 
Join Date: Nov 2010
Posts: 6
Rep Power: 15
Gerard is on a distinguished road
Hi José,

I'm not sure because I have just used OF for a few months now, but maybe you use an older version of OF which automatically applies wallfunctions to omega with SST? Did you check that?

Gerard
Gerard is offline   Reply With Quote

Old   October 19, 2011, 03:10
Default
  #13
New Member
 
Muhammad Umer Ijaz Chaudrey
Join Date: Aug 2011
Location: Eindhoven, The Netherlands
Posts: 26
Rep Power: 15
umer.chaudrey is on a distinguished road
Send a message via Skype™ to umer.chaudrey
Hey everyone,

I have a query regarding K-epsilon values too.

I am using a transient compressible rhoPimpleFoam solver, RAS turbulence model but without any wall functions for k and epsilon.

Whatever I read in various threads on this forum, I understand that initial values of internal fields of k & epsilon are independent of the solution.

But however in my case, if I use different values my solution behaviour and values change. Can anyone please comment on that?

I had written on another thread too, but no one replied regarding specific k & epsilon intial values for internal field. I am not using a velocity boundary condition, rather pressure difference.

Thanks, I will really appreciate some advice.
umer.chaudrey is offline   Reply With Quote

Old   October 19, 2011, 05:20
Default
  #14
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Eindhoven, NL
Posts: 374
Rep Power: 21
romant is on a distinguished road
Quote:
Originally Posted by umer.chaudrey View Post
Hey everyone,

I have a query regarding K-epsilon values too.

I am using a transient compressible rhoPimpleFoam solver, RAS turbulence model but without any wall functions for k and epsilon.

Whatever I read in various threads on this forum, I understand that initial values of internal fields of k & epsilon are independent of the solution.

But however in my case, if I use different values my solution behaviour and values change. Can anyone please comment on that?

I had written on another thread too, but no one replied regarding specific k & epsilon intial values for internal field. I am not using a velocity boundary condition, rather pressure difference.

Thanks, I will really appreciate some advice.

If you use wall function less RANS models (RAS in OF) then you need to check with the turbulence model which boundary conditions it suggests for k and epsilon. those models usually require also that the first cell is within y+ < 1.
__________________
~roman
romant is offline   Reply With Quote

Old   October 19, 2011, 05:23
Default
  #15
New Member
 
Muhammad Umer Ijaz Chaudrey
Join Date: Aug 2011
Location: Eindhoven, The Netherlands
Posts: 26
Rep Power: 15
umer.chaudrey is on a distinguished road
Send a message via Skype™ to umer.chaudrey
First of all thanks for your reply Roman.

But how do I check y+ values in OF when I am not using any wall functions.

According to my knowledge, the yPlusRAS utility in OF only gives y+ values when wall functions are being used. Not without them.

Correct me if I am wrong, as I am only new to this.
umer.chaudrey is offline   Reply With Quote

Old   October 19, 2011, 05:30
Default
  #16
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Eindhoven, NL
Posts: 374
Rep Power: 21
romant is on a distinguished road
you can check with this thread

http://www.cfd-online.com/Forums/ope...-testcase.html

where somebody has posted a tool in order to calculate the y+ values for low reynolds number models which do not use wall functions.
atg likes this.
__________________
~roman
romant is offline   Reply With Quote

Old   October 19, 2011, 05:33
Default
  #17
New Member
 
Muhammad Umer Ijaz Chaudrey
Join Date: Aug 2011
Location: Eindhoven, The Netherlands
Posts: 26
Rep Power: 15
umer.chaudrey is on a distinguished road
Send a message via Skype™ to umer.chaudrey
Roman

Thanks, for the link. I will just have a look at it.
umer.chaudrey is offline   Reply With Quote

Old   October 19, 2011, 05:37
Default One Last Thing
  #18
New Member
 
Muhammad Umer Ijaz Chaudrey
Join Date: Aug 2011
Location: Eindhoven, The Netherlands
Posts: 26
Rep Power: 15
umer.chaudrey is on a distinguished road
Send a message via Skype™ to umer.chaudrey
And Roman, please comment on this last thing before I check the y+ code.

The thing is the documentations I have read about OpenFOAM and in the threads here, for k-epsilon model to calculate the initial values of k and epsilon, I need mean flow velocity and a turbulent length scale.

k= 3/2*(UI)^2 or k = 1/2 U'.U'

I can determine the turbulent length scale based on my geometry. But I do not know the mean flow velocity, because I am generating flow based on pressure difference. The initial velocity everywhere in the model is zero at time zero. The initial conditions I mean.

So I was wondering just for that initial values in my 0 folder, how do I determine the k and epsilon value based on these relations given for my internal field and outlet, even before I run my solution.

Thanks.
umer.chaudrey is offline   Reply With Quote

Old   October 19, 2011, 05:39
Default
  #19
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Eindhoven, NL
Posts: 374
Rep Power: 21
romant is on a distinguished road
I would use an estimate of the mean flow velocity and then run the simulation once. you then will see what your velocity becomes. you can then also see your k values and epsilon values, which can then be used as your inlet conditions.
elham usefi likes this.
__________________
~roman
romant is offline   Reply With Quote

Old   October 19, 2011, 05:42
Default
  #20
New Member
 
Muhammad Umer Ijaz Chaudrey
Join Date: Aug 2011
Location: Eindhoven, The Netherlands
Posts: 26
Rep Power: 15
umer.chaudrey is on a distinguished road
Send a message via Skype™ to umer.chaudrey
Okay thanks again for the help, appreciate it. I'll just do that.
umer.chaudrey is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Bounding epsilon or bounding omega Stylianos OpenFOAM 8 February 23, 2018 13:41
boundary conditions for simpleFoam calculation foam_noob OpenFOAM Running, Solving & CFD 8 July 1, 2015 08:07
Help with k epsilon values of turbulence vw.cfd OpenFOAM 17 December 12, 2013 05:57
Use of k omega turbulence model john_w OpenFOAM Running, Solving & CFD 2 September 22, 2009 05:15
Values of k & epsilon in the moving wall ghlee Main CFD Forum 1 August 16, 1999 12:42


All times are GMT -4. The time now is 00:01.