From aa0b101cad5a60a2465ce7dbe7d6bcc224e21b47 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Fri, 31 Mar 2023 11:54:34 +0100
Subject: [PATCH] MomentumTransferPhaseSystem: Removed unnecessary caching of D

---
 .../MomentumTransferPhaseSystem.C                           | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

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