diff --git a/applications/solvers/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C b/applications/solvers/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
index 8a7d0c1868..2a2f797881 100644
--- a/applications/solvers/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
+++ b/applications/solvers/modules/multiphaseEuler/phaseSystems/PhaseSystems/MomentumTransferPhaseSystem/MomentumTransferPhaseSystem.C
@@ -596,8 +596,10 @@ Foam::MomentumTransferPhaseSystem<BasePhaseSystem>::Fs()
         const phaseInterface& interface =
             turbulentDispersionModelIter()->interface();
 
-        const volScalarField D(turbulentDispersionModelIter()->D());
-        const surfaceScalarField Df(fvc::interpolate(D));
+        const surfaceScalarField Df
+        (
+            fvc::interpolate(turbulentDispersionModelIter()->D())
+        );
 
         const volScalarField alpha12(interface.phase1() + interface.phase2());
         const surfaceScalarField snGradAlpha1By12