diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/alphaEqn.H b/applications/solvers/multiphase/twoPhaseEulerFoam/alphaEqn.H
index f4b97ccd9c551677baa932fb7e8719680f1d4fb1..8ce175b85c4e82dbcc4590525a7a9bccfa6f8dc9 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/alphaEqn.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/alphaEqn.H
@@ -92,6 +92,7 @@
                 )
             );
 
+            // Ensure that the flux at inflow BCs is preserved
             alphaPhic1.boundaryField() = min
             (
                 phi1.boundaryField()*alpha1.boundaryField(),