From 110968a5da63b0b21c2cec9630257d32a054a1de Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Tue, 18 Mar 2014 23:19:41 +0000
Subject: [PATCH] Improved calculation of div(U) using the volumetric flux
---
applications/solvers/multiphase/driftFluxFoam/kEpsilon.H | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/applications/solvers/multiphase/driftFluxFoam/kEpsilon.H b/applications/solvers/multiphase/driftFluxFoam/kEpsilon.H
index 929311548e6..66216b87ce3 100644
--- 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()))));
--
GitLab