diff --git a/applications/solvers/cfdemSolverPimple/UEqn.H b/applications/solvers/cfdemSolverPimple/UEqn.H
index 022516ac..d72c8c23 100644
--- a/applications/solvers/cfdemSolverPimple/UEqn.H
+++ b/applications/solvers/cfdemSolverPimple/UEqn.H
@@ -4,6 +4,7 @@ tmp<fvVectorMatrix> tUEqn
 (
     fvm::ddt(voidfraction,U) - fvm::Sp(fvc::ddt(voidfraction),U)
   + fvm::div(phi,U) - fvm::Sp(fvc::div(phi),U)
+  + particleCloud.divVoidfractionTau(U, voidfraction)
   - fOther/rho
   ==
     fvOptions(U)