[cheng1988sjtu]

Hi All,

I've looked into the kEpsilon model in the twoPhaseEulerFoam, it seems that the equations are problematic:
fvScalarMatrix kEqn
(
fvm::ddt(beta, k)
+ fvm::div(phib, k)
- fvm::laplacian
(
alphak*nuEffb, k,
"laplacian(DkEff,k)"
)
==
beta*G
- fvm::Sp(beta*epsilon/k, k)
);

in the code, volScalarField G = 2*nutb*(tgradUb() && dev(symm(tgradUb())));
tgradUb = fvc::grad(Ub);

I think the terms in red should include the volumetric fraction "beta", any comments? Thank you,