diff --git a/applications/solvers/lagrangian/trackedReactingParcelFoam/pEqn.H b/applications/solvers/lagrangian/trackedReactingParcelFoam/pEqn.H
index b26a6594aa..b6207ce0c4 100644
--- a/applications/solvers/lagrangian/trackedReactingParcelFoam/pEqn.H
+++ b/applications/solvers/lagrangian/trackedReactingParcelFoam/pEqn.H
@@ -41,6 +41,7 @@ if (transonic)
           - lapTerm()
          ==
             reactingParcels.Srho()
+          + pointMassSources.Su()
         );
 
         pEqn.solve();
@@ -80,6 +81,7 @@ else
           - lapTerm()
          ==
             reactingParcels.Srho()
+          + pointMassSources.Su()
         );
 
         pEqn.solve();
diff --git a/applications/solvers/lagrangian/trackedReactingParcelFoam/rhoEqn.H b/applications/solvers/lagrangian/trackedReactingParcelFoam/rhoEqn.H
index 03eed2dc82..9517d2d912 100644
--- a/applications/solvers/lagrangian/trackedReactingParcelFoam/rhoEqn.H
+++ b/applications/solvers/lagrangian/trackedReactingParcelFoam/rhoEqn.H
@@ -37,6 +37,7 @@ Description
       + fvc::div(phi)
       ==
         reactingParcels.Srho()
+      + pointMassSources.Su()
     );
 }