diff --git a/src/TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyz/maxDeltaxyz.C b/src/TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyz/maxDeltaxyz.C
index 351d0fb11f35f966a5bc3be5ca6965e71bd0c7a0..284af6ecf2c7eeb8c4976300356182c556a1c55a 100644
--- a/src/TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyz/maxDeltaxyz.C
+++ b/src/TurbulenceModels/turbulenceModels/LES/LESdeltas/maxDeltaxyz/maxDeltaxyz.C
@@ -64,14 +64,14 @@ void Foam::LESModels::maxDeltaxyz::calcDelta()
             const point& fc = faceC[facei];
             const vector& n = faceN[facei];
 
-            scalar tmp = magSqr(n*(n & (fc - cc)));
+            scalar tmp = mag(n & (fc - cc));
             if (tmp > deltaMaxTmp)
             {
                 deltaMaxTmp = tmp;
             }
         }
 
-        hmax[celli] = deltaCoeff_*Foam::sqrt(deltaMaxTmp);
+        hmax[celli] = deltaCoeff_*deltaMaxTmp;
     }
 
     if (nD == 3)
@@ -113,7 +113,7 @@ Foam::LESModels::maxDeltaxyz::maxDeltaxyz
         dict.optionalSubDict(type() + "Coeffs").lookupOrDefault<scalar>
         (
             "deltaCoeff",
-            1
+            2
         )
     )
 {