CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   kEpsilon.H in twoPhaseEulerFoam (http://www.cfd-online.com/Forums/openfoam/93133-kepsilon-h-twophaseeulerfoam.html)

cheng1988sjtu October 5, 2011 21:28

kEpsilon.H in twoPhaseEulerFoam
 
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,:)


All times are GMT -4. The time now is 10:49.