diff --git a/src/TurbulenceModels/turbulenceModels/RAS/kEpsilon/kEpsilon.C b/src/TurbulenceModels/turbulenceModels/RAS/kEpsilon/kEpsilon.C
index 5303fa1f58f98b776f822a2bffce8aaa3e6a8b75..e4a45b01fb8157e1def74f58fe306873e8caf3d8 100644
--- a/src/TurbulenceModels/turbulenceModels/RAS/kEpsilon/kEpsilon.C
+++ b/src/TurbulenceModels/turbulenceModels/RAS/kEpsilon/kEpsilon.C
@@ -236,7 +236,7 @@ void kEpsilon<BasicTurbulenceModel>::correct()
     volScalarField divU(fvc::div(fvc::absolute(this->phi(), U)));
 
     tmp<volTensorField> tgradU = fvc::grad(U);
-    volScalarField G(this->GName(), nut*(tgradU() && dev(twoSymm(tgradU()))));
+    volScalarField G(this->GName(), nut*(dev(twoSymm(tgradU())) && tgradU()));
     tgradU.clear();
 
     // Update epsilon and G at the wall
@@ -260,7 +260,6 @@ void kEpsilon<BasicTurbulenceModel>::correct()
     solve(epsEqn);
     bound(epsilon_, this->epsilonMin_);
 
-
     // Turbulent kinetic energy equation
     tmp<fvScalarMatrix> kEqn
     (