diff --git a/applications/solvers/lagrangian/kinematicParcelFoam/UEqn.H b/applications/solvers/lagrangian/kinematicParcelFoam/UEqn.H
index 0e4177b8e1..c2ea5f8bc7 100644
--- a/applications/solvers/lagrangian/kinematicParcelFoam/UEqn.H
+++ b/applications/solvers/lagrangian/kinematicParcelFoam/UEqn.H
@@ -6,7 +6,7 @@
       + MRF.DDt(U)
       + turbulence->divDevReff(U)
      ==
-        parcels.SU(U, true)
+        invRhoInf*parcels.SU(U)
       + fvOptions(U)
     );
 
diff --git a/applications/solvers/lagrangian/kinematicParcelFoam/createFields.H b/applications/solvers/lagrangian/kinematicParcelFoam/createFields.H
index 0182144d86..67eda1592a 100644
--- a/applications/solvers/lagrangian/kinematicParcelFoam/createFields.H
+++ b/applications/solvers/lagrangian/kinematicParcelFoam/createFields.H
@@ -39,6 +39,11 @@ dimensionedScalar rhoInfValue
     laminarTransport
 );
 
+const dimensionedScalar invRhoInf
+(
+    dimless/dimDensity, scalar(1)/rhoInfValue.value()
+);
+
 volScalarField rhoInf
 (
     IOobject