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

(incompressible) interFoam - where is the density in the k, epsilon and continuity eq

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 22, 2015, 15:40
Default (incompressible) interFoam - where is the density in the k, epsilon and continuity eq
  #1
Member
 
Join Date: Aug 2011
Posts: 82
Rep Power: 6
idefix is on a distinguished road
Hello,

the following equation for epsilon is used to calculate this value when using the k-epsilon turbulence model

fvm::ddt(epsilon_)
+ fvm::div(phi_, epsilon_)
-fvm::Sp(fvc::div(phi_), epsilon_)
- fvm::laplacian(DepsilonEff(), epsilon_)
==
C1_*G*epsilon_/k_
- fvm::Sp(C2_*epsilon_/k_, epsilon_)



Now I want to use the k-epsilon turbulence model together with the interFoam-solver, which is incompressible.
Do I have to include the densitiy for example in
fvm::ddt(epsilon_)
--> fvm::ddt(epsilon_ * density) ?

I am asking because the solver is described in the UserGuide as solver for two incompressible fluids. That is clear for pure gas and pure liquid phase.
But what happens for the area where the interface is?

Is it necessary there to add the density in the epsilon-equation?

I foundn in Rusche´s thesis (eq. 3.25) an explanation that the velocity is weighted with the liquid volume fraction. Is this valid for k and epsilon too?
If this is the case where does it happen in the code?


What happens to the contiunity equation?
Is it:
d rho / dt + d (rho * u_i)/ dx_i ?

Thanks a lot for your help
idefix 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
simpleFoam error - "Floating point exception" mbcx4jc2 OpenFOAM Running, Solving & CFD 12 August 4, 2015 02:20
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 13 November 4, 2013 15:13
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
dynamic Mesh is faster than MRF???? sharonyue OpenFOAM Running, Solving & CFD 14 August 26, 2013 07:47
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 06:41.