diff --git a/applications/modules/twoPhaseSolver/alphaPredictor.C b/applications/modules/twoPhaseSolver/alphaPredictor.C
index 874c111364..96037c1d51 100644
--- a/applications/modules/twoPhaseSolver/alphaPredictor.C
+++ b/applications/modules/twoPhaseSolver/alphaPredictor.C
@@ -256,8 +256,8 @@ void Foam::solvers::twoPhaseSolver::alphaSolve
                     alpha1,
                     talphaPhi1Un(),
                     talphaPhi1Corr.ref(),
-                    Sp(),
-                    (-Sp()*alpha1)(),
+                    (Sp() + divU())(),
+                    (-(Sp() + divU())*alpha1)(),
                     oneField(),
                     zeroField()
                 );