diff --git a/applications/modules/incompressibleDenseParticleFluid/momentumPredictor.C b/applications/modules/incompressibleDenseParticleFluid/momentumPredictor.C
index edb602fbbe..84c4f27491 100644
--- a/applications/modules/incompressibleDenseParticleFluid/momentumPredictor.C
+++ b/applications/modules/incompressibleDenseParticleFluid/momentumPredictor.C
@@ -42,7 +42,7 @@ void Foam::solvers::incompressibleDenseParticleFluid::momentumPredictor()
       - fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc)
       + momentumTransport->divDevSigma(Uc)
      ==
-        fvModels().source(Uc)
+        fvModels().source(alphac, Uc)
     );
     fvVectorMatrix& UcEqn = tUcEqn.ref();