From 92ef068675f205a00e7367f26d81c5c376575fcf Mon Sep 17 00:00:00 2001
From: Henry Weller <http://openfoam.org>
Date: Tue, 29 Aug 2017 14:48:32 +0100
Subject: [PATCH] compressibleMultiphaseInterFoam: Corrected update of dgdt for
 multiple phases

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

diff --git a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/pEqn.H b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/pEqn.H
index 649063da7..f0d7e05ed 100644
--- a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/pEqn.H
+++ b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/pEqn.H
@@ -106,6 +106,8 @@
                 phase().dgdt() =
                     pos0(phase())
                   *(p_rghEqnComps[phasei] & p_rgh)/phase().thermo().rho();
+
+                phasei++;
             }
 
             phi = phiHbyA + p_rghEqnIncomp.flux();