From cfdd46ef155dc8b78dc6d8c6e75cf0da98446977 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Mon, 17 Mar 2014 22:38:53 +0000
Subject: [PATCH] Separate convection and diffusion terms in the MULEScorr
 predictor to remove interaction of the diffusion term with the
 flux-correction.

---
 applications/solvers/multiphase/driftFluxFoam/alphaEqn.H | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/applications/solvers/multiphase/driftFluxFoam/alphaEqn.H b/applications/solvers/multiphase/driftFluxFoam/alphaEqn.H
index 817a0c9124..d56b371619 100644
--- a/applications/solvers/multiphase/driftFluxFoam/alphaEqn.H
+++ b/applications/solvers/multiphase/driftFluxFoam/alphaEqn.H
@@ -17,6 +17,11 @@
                 phi,
                 upwind<scalar>(mesh, phi)
             ).fvmDiv(phi, alpha1)
+        );
+
+        solve
+        (
+            alpha1Eqn
           - fv::gaussLaplacianScheme<scalar, scalar>
             (
                 mesh,
@@ -25,8 +30,6 @@
             ).fvmLaplacian(fvc::interpolate(mut/rho), alpha1)
         );
 
-        alpha1Eqn.solve();
-
         Info<< "Phase-1 volume fraction = "
             << alpha1.weightedAverage(mesh.Vsc()).value()
             << "  Min(alpha1) = " << min(alpha1).value()