Bugfix: use function form of transpose

applications/solvers/DNS/dnsFoam/globalProperties.H

@@ -9,7 +9,7 @@
              << (
                      0.5*nu*average
                      (
-                         magSqr(fvc::grad(U) + fvc::grad(U)().T())
+                         magSqr(fvc::grad(U) + T(fvc::grad(U)))
                      )
                 ).value() << endl;
 

applications/solvers/combustion/PDRFoam/PDRModels/turbulence/PDRkEpsilon/PDRkEpsilon.C

@@ -219,7 +219,7 @@ tmp<fvVectorMatrix> PDRkEpsilon::divDevRhoReff(volVectorField& U) const
 {
     return
     (
-      - fvm::laplacian(muEff(), U) - fvc::div(muEff()*dev2(fvc::grad(U)().T()))
+      - fvm::laplacian(muEff(), U) - fvc::div(muEff()*dev2(T(fvc::grad(U))))
     );
 }
 

applications/solvers/compressible/rhoCentralFoam/rhoCentralFoam.C

@@ -139,7 +139,7 @@ int main(int argc, char *argv[])
           + aphiv_neg*(rho_neg*(e_neg + 0.5*magSqr(U_neg)) + p_neg)
           + aSf*p_pos - aSf*p_neg;
 
-        volTensorField tauMC("tauMC", mu*dev2(fvc::grad(U)().T()));
+        volTensorField tauMC("tauMC", mu*dev2(T(fvc::grad(U))));
 
         // --- Solve density
         solve(fvm::ddt(rho) + fvc::div(phi));

applications/solvers/multiphase/compressibleInterDyMFoam/UEqn.H

@@ -11,7 +11,7 @@
       + fvm::div(rhoPhi, U)
       - fvm::laplacian(muEff, U)
       - (fvc::grad(U) & fvc::grad(muEff))
-    //- fvc::div(muf*(mesh.Sf() & fvc::interpolate(fvc::grad(U)().T())))
+      - fvc::div(muEff*(mesh.Sf() & fvc::interpolate(T(fvc::grad(U)))))
     );
 
     UEqn.relax();

applications/solvers/multiphase/multiphaseInterFoam/UEqn.H

@@ -6,5 +6,5 @@
       + fvm::div(mixture.rhoPhi(), U)
       - fvm::laplacian(muf, U)
       - (fvc::grad(U) & fvc::grad(muf))
-    //- fvc::div(muf*(mesh.Sf() & fvc::interpolate(fvc::grad(U)().T())))
+    //- fvc::div(muf*(mesh.Sf() & fvc::interpolate(T(fvc::grad(U)))))
     );