diff --git a/applications/solvers/multiphase/bubbleFoam/UEqns.H b/applications/solvers/multiphase/bubbleFoam/UEqns.H
index 82d53d84b0093c70f0ed81385a35c6fffba17eba..3239df21cb3385c4c9a55eece0d8de26e67530bd 100644
--- a/applications/solvers/multiphase/bubbleFoam/UEqns.H
+++ b/applications/solvers/multiphase/bubbleFoam/UEqns.H
@@ -2,7 +2,7 @@ fvVectorMatrix UaEqn(Ua, Ua.dimensions()*dimVol/dimTime);
 fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
 
 {
-    volTensorField Rca(-nuEffa*(fvc::grad(Ua)().T()));
+    volTensorField Rca(-nuEffa*(T(fvc::grad(Ua))));
     Rca = Rca + (2.0/3.0)*sqr(Ct)*I*k - (2.0/3.0)*I*tr(Rca);
 
     surfaceScalarField phiRa
@@ -36,7 +36,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
     UaEqn.relax();
 
 
-    volTensorField Rcb(-nuEffb*fvc::grad(Ub)().T());
+    volTensorField Rcb(-nuEffb*T(fvc::grad(Ub)));
     Rcb = Rcb + (2.0/3.0)*I*k - (2.0/3.0)*I*tr(Rcb);
 
     surfaceScalarField phiRb
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H b/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
index b07f3e70dd396e95e9f38637a6b5ba4fb06c7fbb..0a590e6bf936ac844072a7d3ee9ebed78b6c1658 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
@@ -3,7 +3,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
 
 {
     {
-        volTensorField gradUaT(fvc::grad(Ua)().T());
+        volTensorField gradUaT(T(fvc::grad(Ua)));
 
         if (kineticTheory.on())
         {
@@ -58,7 +58,7 @@ fvVectorMatrix UbEqn(Ub, Ub.dimensions()*dimVol/dimTime);
     }
 
     {
-        volTensorField gradUbT(fvc::grad(Ub)().T());
+        volTensorField gradUbT(T(fvc::grad(Ub)));
         volTensorField Rcb
         (
             "Rcb",