diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H b/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
index 6c765ec66b9235efced56355b7b3047732ef7f45..abd3c8b6461699285aa81682b71291e30c8d4a87 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
@@ -31,6 +31,7 @@ volScalarField dragCoeff(fluid.dragCoeff());
         U1Eqn += fvm::Sp(dragCoeff, U1);
         fvOptions.constrain(U1Eqn);
         U1.correctBoundaryConditions();
+        fvOptions.correct(U1);
     }
 
     {
@@ -54,5 +55,6 @@ volScalarField dragCoeff(fluid.dragCoeff());
         U2Eqn += fvm::Sp(dragCoeff, U2);
         fvOptions.constrain(U2Eqn);
         U2.correctBoundaryConditions();
+        fvOptions.correct(U2);
     }
 }