From be16ebf24f44873981b895c6b1ecd017f31fb4de 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 e352bbe914d..19f2d473b9b 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));
--
GitLab