diff --git a/src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/qZeta/qZeta.C b/src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/qZeta/qZeta.C
index b9f99c117c..a12b61f1a6 100644
--- a/src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/qZeta/qZeta.C
+++ b/src/TurbulenceModels/incompressible/turbulentTransportModels/RAS/qZeta/qZeta.C
@@ -49,7 +49,7 @@ tmp<volScalarField> qZeta::fMu() const
 
     if (anisotropic_)
     {
-        return exp((-scalar(2.5) + Rt/20.0)/pow(scalar(1) + Rt/130.0, 3.0));
+        return exp((-scalar(2.5) + Rt/20.0)/pow3(scalar(1) + Rt/130.0));
     }
     else
     {