From 9e79d990f3fcc66465a6b45c2ea4b9b0ef125e1f Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Tue, 20 Aug 2013 12:30:44 +0100
Subject: [PATCH] multiphaseEulerFoam: Ensure that the flux at inflow BCs is
 preserved

---
 .../multiphaseSystem/multiphaseSystem.C                    | 7 +++++++
 1 file changed, 7 insertions(+)

diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C
index 2a472a943f..f347de82a0 100644
--- a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C
+++ b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C
@@ -125,6 +125,13 @@ void Foam::multiphaseSystem::solveAlphas()
             );
         }
 
+        // Ensure that the flux at inflow BCs is preserved
+        phiAlphaCorr.boundaryField() = min
+        (
+            phase1.phi().boundaryField()*alpha1.boundaryField(),
+            phiAlphaCorr.boundaryField()
+        );
+
         MULES::limit
         (
             geometricOneField(),