8 dimensionedScalar k0("k0", k.dimensions(), SMALL);
9 dimensionedScalar epsilon0("epsilon0", epsilon.dimensions(), SMALL);
11 volScalarField divU = fvc::div(phi/fvc::interpolate(rho));
13 tmp<volTensorField> tgradU = fvc::grad(U);
14 volScalarField G = 2*mut*(tgradU() && dev(symm(tgradU())));
17 volScalarField Gcoef =
18 alphak*Cmu*k*(g & fvc::grad(rho))/(epsilon + epsilon0);
20 # include "wallFunctions.H"
22 // Dissipation equation
25 fvm::ddt(rho, epsilon)
26 + fvm::div(phi, epsilon)
29 alphaEps*mut + mul, epsilon,
30 "laplacian(DepsilonEff,epsilon)"
34 - fvm::SuSp(C1*(1.0 - C3)*Gcoef + (2.0/3.0*C1)*rho*divU, epsilon)
35 - fvm::Sp(C2*rho*epsilon/k, epsilon)
38 # include "wallDissipation.H"
43 bound(epsilon, epsilon0);
46 // Turbulent kinetic energy equation
59 - fvm::SuSp(Gcoef + 2.0/3.0*rho*divU, k)
60 - fvm::Sp(rho*epsilon/k, k)
69 //- Re-calculate viscosity
70 mut = rho*Cmu*sqr(k)/epsilon;
72 # include "wallViscosity.H"