From 852726c25fbe6b4bd4e947bec7ff32bc92bb103e Mon Sep 17 00:00:00 2001
From: Henry Weller <http://openfoam.org>
Date: Wed, 7 Feb 2018 15:36:51 +0000
Subject: [PATCH] compressibleInterDyMFoam: Improved conservation for morphing
 mesh cases

Resolves bug-report https://bugs.openfoam.org/view.php?id=2824
---
 .../compressibleInterFoam/compressibleInterDyMFoam/pEqn.H       | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/pEqn.H b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/pEqn.H
index dd949f5da0..063e2145f1 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/pEqn.H
+++ b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/pEqn.H
@@ -82,6 +82,7 @@
                   - (fvOptions(alpha1, mixture.thermo1().rho())&rho1)
                 )/rho1
               - fvc::ddt(alpha1) - fvc::div(alphaPhi1)
+              + fvc::div(mesh.phi())*alpha1
               + (alpha1*psi1/rho1)*correction(fvm::ddt(p_rgh))
             );
 
@@ -93,6 +94,7 @@
                  - (fvOptions(alpha2, mixture.thermo2().rho())&rho2)
                )/rho2
              - fvc::ddt(alpha2) - fvc::div(alphaPhi2)
+             + fvc::div(mesh.phi())*alpha2
              + (alpha2*psi2/rho2)*correction(fvm::ddt(p_rgh))
             );
     }