diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/diameterModels/velocityGroup/sizeGroup/shapeModels/fractal/fractal.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/diameterModels/velocityGroup/sizeGroup/shapeModels/fractal/fractal.C
index eaa5ef2301..62fadc35fd 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/diameterModels/velocityGroup/sizeGroup/shapeModels/fractal/fractal.C
+++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/diameterModels/velocityGroup/sizeGroup/shapeModels/fractal/fractal.C
@@ -27,6 +27,7 @@ License
 #include "sinteringModel.H"
 #include "fvmDdt.H"
 #include "fvmDiv.H"
+#include "fvcDiv.H"
 #include "fvmSup.H"
 #include "fvcSup.H"
 #include "fvcDdt.H"
@@ -214,7 +215,15 @@ void Foam::diameterModels::shapeModels::fractal::correct()
         fvc::ddt(alpha, rho, fi)*kappa_.oldTime()
       + alpha*rho*fi*fvm::ddt(kappa_)
       + fvm::div(fAlphaRhoPhi, kappa_)
-      + fvm::SuSp(-phase.continuityError()*fi, kappa_)
+      + fvm::SuSp
+        (
+            fi
+           *(
+                fi.VelocityGroup().dmdt()
+              - (fvc::ddt(alpha, rho) + fvc::div(phase.alphaRhoPhi()))
+            ),
+            kappa_
+        )
       ==
       - sinteringModel_->R()
       + fvc::Su(Su_*rho, kappa_)
diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C
index dd176dafaa..c64d72eb54 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C
+++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C
@@ -1304,7 +1304,12 @@ void Foam::diameterModels::populationBalanceModel::solve()
                         phase.alphaRhoPhi(),
                         fi
                     )
-                  + fvm::SuSp(-phase.continuityError(), fi)
+                  + fvm::SuSp
+                    (
+                        fi.VelocityGroup().dmdt()
+                      - (fvc::ddt(alpha, rho) + fvc::div(phase.alphaRhoPhi())),
+                        fi
+                    )
                   ==
                     fvc::Su(Su_[i]*rho, fi)
                   - fvm::SuSp(SuSp_[i]*rho, fi)