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(),