From f1f2c0a9f123f487607b65fd55cd3e5d46f4c4ed Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Thu, 2 Jul 2015 22:50:44 +0100
Subject: [PATCH] reactingTwoPhaseEulerFoam/pU/pEqn: Added missing dmdt terms

---
 .../solvers/multiphase/reactingTwoPhaseEulerFoam/pU/pEqn.H    | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/pU/pEqn.H b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/pU/pEqn.H
index e352bbe91..19f2d473b 100644
--- a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/pU/pEqn.H
+++ b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/pU/pEqn.H
@@ -278,14 +278,14 @@ while (pimple.correct())
     {
         pEqnComp1 =
             (
-                phase1.continuityError()
+                phase1.continuityError() - fluid.dmdt()
               - fvc::Sp(fvc::ddt(alpha1) + fvc::div(alphaPhi1), rho1)
             )/rho1
           + (alpha1*psi1/rho1)*correction(fvm::ddt(p_rgh));
 
         pEqnComp2 =
             (
-                phase2.continuityError()
+                phase2.continuityError() + fluid.dmdt()
               - fvc::Sp(fvc::ddt(alpha2) + fvc::div(alphaPhi2), rho2)
             )/rho2
           + (alpha2*psi2/rho2)*correction(fvm::ddt(p_rgh));