diff --git a/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H b/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H
index 5dc9058d9e..3c9cfdeade 100644
--- a/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H
+++ b/applications/solvers/lagrangian/LTSReactingParcelFoam/pEqn.H
@@ -31,7 +31,7 @@
         parcels.Srho()
     );
 
-    sources.apply(pDDtEqn);
+    sources.apply(pDDtEqn, rho.name());
 
     while (pimple.correctNonOrthogonal())
     {
diff --git a/applications/solvers/lagrangian/coalChemistryFoam/pEqn.H b/applications/solvers/lagrangian/coalChemistryFoam/pEqn.H
index 5ef8472210..b8bb134522 100644
--- a/applications/solvers/lagrangian/coalChemistryFoam/pEqn.H
+++ b/applications/solvers/lagrangian/coalChemistryFoam/pEqn.H
@@ -26,7 +26,7 @@ if (pimple.transonic())
             coalParcels.Srho()
         );
 
-        sources.apply(pEqn);
+        sources.apply(pEqn, rho.name());
 
         pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
 
@@ -56,7 +56,7 @@ else
             coalParcels.Srho()
         );
 
-        sources.apply(pEqn);
+        sources.apply(pEqn, rho.name());
 
         pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
 
diff --git a/applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/pEqn.H b/applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/pEqn.H
index 581c03e104..705a4bae92 100644
--- a/applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/pEqn.H
+++ b/applications/solvers/lagrangian/porousExplicitSourceReactingParcelFoam/pEqn.H
@@ -31,7 +31,7 @@
         parcels.Srho()
     );
 
-    sources.apply(pDDtEqn);
+    sources.apply(pDDtEqn, rho.name());
 
     while (pimple.correctNonOrthogonal())
     {