[idefix]

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