compressibleMultiphaseInterFoam: Corrected update of dgdt for multiple phases

Resolves bug-report https://bugs.openfoam.org/view.php?id=2677
2017-08-29 14:49:44 +01:00
parent 3fe7aa77e6
commit adac0dbb1d
1 changed files with 2 additions and 0 deletions
--- 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();