Improved calculation of div(U) using the volumetric flux

110968a5 · Henry · Andrew Heather · acb63f51 · 110968a5
Commit 110968a5 authored 11 years ago by Henry Committed by Andrew Heather 11 years ago
--- a/applications/solvers/multiphase/driftFluxFoam/kEpsilon.H
+++ b/applications/solvers/multiphase/driftFluxFoam/kEpsilon.H
@@ -10,7 +10,7 @@ if (turbulence)
    dimensionedScalar epsilon0("epsilon0", epsilon.dimensions(), 0);
    dimensionedScalar epsilonMin("epsilonMin", epsilon.dimensions(), SMALL);

-    volScalarField divU(fvc::div(rhoPhi/fvc::interpolate(rho)));
+    volScalarField divU(fvc::div(phi));

    tmp<volTensorField> tgradU = fvc::grad(U);
    volScalarField G(mut*(tgradU() && dev(twoSymm(tgradU()))));