[remidemol]

Hello foamers,

I'm a bit confuse about the different formulations of mixtureKEpsilon models. It seems to me that in the version 3.0.1 mixtureKEpsilon model is written with an incompressible form, which is contradictory with the description of the twoPhaseEulerFoam I use.

Do anyone have an explanation ? Maybe it's written in an other way.

Thank you for your help !

In version 2.3 it's like this (I suppose compressible form)

tmp<fvScalarMatrix> epsEqn
(
fvm::ddt(rhom, epsilonm)
+ fvm::div(phim, epsilonm)
- fvm::Sp(fvc::ddt(rhom) + fvc::div(phim), epsilonm)
- fvm::laplacian(DepsilonEff(rhom*nutm), epsilonm)
==
C1_*rhom*Gm*epsilonm/km
- fvm::SuSp(((2.0/3.0)*C1_)*rhom*divUm, epsilonm)
- fvm::Sp(C2_*rhom*epsilonm/km, epsilonm)
+ epsilonSource()
);

Whereas in version 3.0.1 it is (I suppose incompressible form)

    tmp<fvScalarMatrix> epsEqn
    (
        fvm::ddt(epsilonm)
      + fvm::div(phim, epsilonm)
      - fvm::Sp(fvc::div(phim), epsilonm)
      - fvm::laplacian(DepsilonEff(nutm), epsilonm)
     ==
        C1_*Gm*epsilonm/km
      - fvm::SuSp(((2.0/3.0)*C1_)*divUm, epsilonm)
      - fvm::Sp(C2_*epsilonm/km, epsilonm)
      + epsilonSource()
    );