|December 22, 2015, 15:40||
(incompressible) interFoam - where is the density in the k, epsilon and continuity eq
Join Date: Aug 2011
Posts: 82Rep Power: 6
the following equation for epsilon is used to calculate this value when using the k-epsilon turbulence model
+ fvm::div(phi_, epsilon_)
- fvm::laplacian(DepsilonEff(), epsilon_)
- 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_ * 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?
d rho / dt + d (rho * u_i)/ dx_i ?
Thanks a lot for your help
|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|