diff --git a/src/functionObjects/forces/forces/forces.C b/src/functionObjects/forces/forces/forces.C
index e3dbc07f1f967eb828828a3555061f1f431a7a5f..70a34aba1877f5a17ef5b5e24cc1992ca24fb686 100644
--- a/src/functionObjects/forces/forces/forces.C
+++ b/src/functionObjects/forces/forces/forces.C
@@ -350,27 +350,27 @@ Foam::functionObjects::forces::devRhoReff() const
     typedef compressible::turbulenceModel cmpTurbModel;
     typedef incompressible::turbulenceModel icoTurbModel;
 
+    const auto& U = lookupObject<volVectorField>(UName_);
+
     if (foundObject<cmpTurbModel>(cmpTurbModel::propertiesName))
     {
         const cmpTurbModel& turb =
             lookupObject<cmpTurbModel>(cmpTurbModel::propertiesName);
 
-        return turb.devRhoReff();
+        return turb.devRhoReff(U);
     }
     else if (foundObject<icoTurbModel>(icoTurbModel::propertiesName))
     {
         const incompressible::turbulenceModel& turb =
             lookupObject<icoTurbModel>(icoTurbModel::propertiesName);
 
-        return rho()*turb.devReff();
+        return rho()*turb.devReff(U);
     }
     else if (foundObject<fluidThermo>(fluidThermo::dictName))
     {
         const fluidThermo& thermo =
             lookupObject<fluidThermo>(fluidThermo::dictName);
 
-        const volVectorField& U = lookupObject<volVectorField>(UName_);
-
         return -thermo.mu()*dev(twoSymm(fvc::grad(U)));
     }
     else if (foundObject<transportModel>("transportProperties"))
@@ -378,8 +378,6 @@ Foam::functionObjects::forces::devRhoReff() const
         const transportModel& laminarT =
             lookupObject<transportModel>("transportProperties");
 
-        const volVectorField& U = lookupObject<volVectorField>(UName_);
-
         return -rho()*laminarT.nu()*dev(twoSymm(fvc::grad(U)));
     }
     else if (foundObject<dictionary>("transportProperties"))
@@ -389,8 +387,6 @@ Foam::functionObjects::forces::devRhoReff() const
 
         dimensionedScalar nu("nu", dimViscosity, transportProperties);
 
-        const volVectorField& U = lookupObject<volVectorField>(UName_);
-
         return -rho()*nu*dev(twoSymm(fvc::grad(U)));
     }
     else