diff --git a/applications/solvers/lagrangian/denseParticleFoam/denseParticleFoam.C b/applications/solvers/lagrangian/denseParticleFoam/denseParticleFoam.C
index 00f796597e..1e7d6dec4f 100644
--- a/applications/solvers/lagrangian/denseParticleFoam/denseParticleFoam.C
+++ b/applications/solvers/lagrangian/denseParticleFoam/denseParticleFoam.C
@@ -164,13 +164,10 @@ int main(int argc, char *argv[])
         // Face continuous-dispersed phase drag coefficient
         const surfaceScalarField Dcf(fvc::interpolate(Dc));
 
-        // Face dispersed phase drag force
-        const surfaceScalarField Fdf(fvc::flux(Fd));
-
         // Effective flux of the dispersed phase
         const surfaceScalarField phid
         (
-            Fdf/(Dcf + dimensionedScalar(Dc.dimensions(), small))
+            fvc::flux(Fd)/(Dcf + dimensionedScalar(Dc.dimensions(), small))
         );
 
         // Face buoyancy force
diff --git a/applications/solvers/lagrangian/denseParticleFoam/pEqn.H b/applications/solvers/lagrangian/denseParticleFoam/pEqn.H
index 2e7dbf9c94..124bb89a4f 100644
--- a/applications/solvers/lagrangian/denseParticleFoam/pEqn.H
+++ b/applications/solvers/lagrangian/denseParticleFoam/pEqn.H
@@ -3,7 +3,7 @@
     const volScalarField r1ADUc(1/(1 + rAUc*Dc));
 
     const surfaceScalarField rAUcf(fvc::interpolate(rAUc));
-    const surfaceScalarField r1ADUcf(1/(1 + rAUcf*fvc::interpolate(Dc)));
+    const surfaceScalarField r1ADUcf(1/(1 + rAUcf*Dcf));
     const surfaceScalarField rADUcf("Dp", r1ADUcf*rAUcf);
 
     volVectorField HbyA(constrainHbyA(rAUc*UcEqn.H(), Uc, p));
@@ -27,7 +27,7 @@
         fvc::makeAbsolute(phiHbyAD, Uc);
     }
 
-    phiHbyAD += rADUcf*(Fgf + Fdf);
+    phiHbyAD += rADUcf*(Fgf + Dcf*phid);
 
     // Update the pressure BCs to ensure flux consistency
     constrainPressure(p, Uc, phiHbyAD, rADUcf);