CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

How to find out k and epsilon

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

Like Tree5Likes
  • 2 Post By fredo490
  • 1 Post By fredo490
  • 1 Post By fredo490
  • 1 Post By fredo490

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 31, 2013, 07:51
Default How to find out k and epsilon
  #1
Senior Member
 
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 12
sam.ho is on a distinguished road
Hi Foamers,

I am simulating flow through valve in which i know only pressure at inlet and outlet of a valve. In this scenario how to find out k and epsilon ?
Or how to give boundary condition to inlet , outlet and walls ?
Inlet pressure = 5 bar
outlet pressure = 1 bar ( atmosphere )
Density = 800 kg/m3
dynamic viscosity = 6.68 X 10^-3 Pa s
sam.ho is offline   Reply With Quote

Old   December 31, 2013, 10:39
Default
  #2
Senior Member
 
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 249
Rep Power: 16
fredo490 is on a distinguished road
You have a tool on this website to help you: http://www.cfd-online.com/Tools/turbulence.php

After, it depends of your upstream geometry and properties (if you consider it fully turbulent ?).
Antimony and blttkgl like this.
fredo490 is offline   Reply With Quote

Old   January 2, 2014, 05:45
Default
  #3
Senior Member
 
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 12
sam.ho is on a distinguished road
Hi Fredo,

I have cross checks all the parameters and all are correct. As explained below

Boundary Conditions :
Inlet pressure = 5 Bar
Outlet Pressure = 1 Bar ( atmospheric pressure)
Density of fluid = 800 kg/m3
Dynamic Viscosity = 6.68 X 10^-3 Pa s

As i did not knew how to calculate k and epsilon values from pressure so i got nominal flow rate from ANSYS 14.5 results and calculated k and epsilon as follows
k = 0.00036 m2/s2
epsilon = 0.00014 m2/s3
and nuT = 1.04 X 10^-7 m2/s Calculated as muT= C_mu X k^2/epsilon
sam.ho is offline   Reply With Quote

Old   January 2, 2014, 11:06
Default
  #4
Senior Member
 
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 249
Rep Power: 16
fredo490 is on a distinguished road
There is no direct relation between the pressure and k/epsilon. You should seriously learn more about the turbulence theory ! There is no defaut value that always works... Even Fluent asks for some parameters to estimate k and epsilon (at the bottom of the inlet condition panel).

Usually people use the turbulent intensity to estimate k and epsilon. The turbulent intensity at the inlet highly depends of your case/geometry before your domain. If your inlet is a long pipe, it can be considered as fully turbulent and you can use the hydrolic diameter; if the inlet is "laminar" you need to estimate the level of turbulence intensity (it is never 0 in real cases)...
arashjkh likes this.
fredo490 is offline   Reply With Quote

Old   January 3, 2014, 02:21
Post
  #5
Senior Member
 
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 12
sam.ho is on a distinguished road
Hi Fredo,

I calculated turbulent intensity as 3% and then proceeded to calculate k, epsilon and nuT.
And I got the values as mentioned above.
Yesterday I changed the fvSolutions and fvSchemes by taking reference of motorBike tutorial and ran this case with k-omega SST model. It ran well but the velocity and pressure and tremendously high.
As I gave inlet pressure as 625 m2/s2 and outlet pressure as 125 m2/s2 which i obtained by dividing pressures with density.

May any body know why this behaviour of the simulation ?
sam.ho is offline   Reply With Quote

Old   January 4, 2014, 20:37
Default
  #6
Senior Member
 
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 21
jherb is on a distinguished road
A few thoughts:

Have you tried to use the results of your simulation with the k-epsilon turbulence model (especially the pressure and velocities) as initial conditions for your simulation with k-omega(SST)?

There are special boundary conditions for the turbulence related fields to use turbulence intensity and turbulent mixing length:
https://github.com/OpenFOAM/OpenFOAM...hScalarField.H
https://github.com/OpenFOAM/OpenFOAM...hScalarField.H
https://github.com/OpenFOAM/OpenFOAM...ld.H?source=cc
jherb is offline   Reply With Quote

Old   January 6, 2014, 05:08
Default
  #7
Senior Member
 
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 12
sam.ho is on a distinguished road
Hi jehrb,

Could you please explain how to calculate Turbulent Mixing length. ?
Because mine is a Ball Valve simulation .

I have used turbulent viscosity ratio for the calculation of k, epsilon and omega .
sam.ho is offline   Reply With Quote

Old   January 6, 2014, 06:09
Default
  #8
Senior Member
 
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 249
Rep Power: 16
fredo490 is on a distinguished road
Please read about the turbulence theory... You missunderstand many things.

In your case, your flow comes from a "long" pipe and it is often a good assumption to say that the flow is fully turbulent.

Therefor you can use the turbulence length scale and the turbulence intensity to get all your variables (k, epsilon, omega).
- There is an empirical relation between the turbulence length scale and the pipe diameter.
-There is an empirical relation between the turbulence intensity and the hydraulic diameter

At least read some documents like this:
http://jullio.pe.kr/fluent6.1/help/html/ug/node178.htm

And look in google: "turbulence in tube" ...
sam.ho likes this.
fredo490 is offline   Reply With Quote

Old   January 6, 2014, 06:29
Default
  #9
Senior Member
 
Join Date: Sep 2013
Location: Bangalore India
Posts: 134
Rep Power: 12
sam.ho is on a distinguished road
Hi ,

Thank you for the explanation.
But my case is not like that.. I know only pressures at inlet and outlet as boundary condition so that leading me to a lot of confusion.
Time being i took the flow rate from ANSYS results and calculated k, epsilon and omega.
And now i will calculate the values according to your suggestion and test my case.
Thank a ton
sam.ho is offline   Reply With Quote

Old   January 6, 2014, 17:35
Default
  #10
Senior Member
 
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 249
Rep Power: 16
fredo490 is on a distinguished road
jherb gave you the suitable boundary conditions. Indeed you don't have the velocity at your inlet but you can run a first simulation with I (turbulence intensity) at 5% and correct the value once you get a first convergence.

k is function of the Reynolds (and then the velocity), but epsilon/omega are only function of the diameter of your pipe (turbulent length scale).

Just use a first rough estimation of k and correct it in a second simulation that will be more accurate.

Quote:
Originally Posted by jherb View Post
There are special boundary conditions for the turbulence related fields to use turbulence intensity and turbulent mixing length:
https://github.com/OpenFOAM/OpenFOAM...hScalarField.H
https://github.com/OpenFOAM/OpenFOAM...hScalarField.H
https://github.com/OpenFOAM/OpenFOAM...ld.H?source=cc
sam.ho likes this.
fredo490 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
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 16 March 4, 2017 09:30
calculation stops after few time steps sivakumar OpenFOAM Running, Solving & CFD 7 March 17, 2013 07:37
Calculation of k and epsilon freezes Nigirim OpenFOAM Running, Solving & CFD 1 November 14, 2012 08:52
epsilon and K blowing up. sivakumar OpenFOAM Running, Solving & CFD 1 October 25, 2012 05:50
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 1 November 25, 2008 21:21


All times are GMT -4. The time now is 17:59.