diff --git a/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C b/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C
index 8f1c2f7c..809ba052 100644
--- a/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C
+++ b/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimple.C
@@ -121,6 +121,7 @@ int main(int argc, char *argv[])
         {
             #include "rhoEqn.H"
         }
+        rhoeps = rho*voidfraction;
 
         // --- Pressure-velocity PIMPLE corrector loop
         while (pimple.loop())