diff --git a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUWallFunction/nutUWallFunctionFvPatchScalarField.C b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUWallFunction/nutUWallFunctionFvPatchScalarField.C
index af072fe580cda356498053548975895cbe0972f2..f64f65edcac371cda6b58bd780dade7dce81375a 100644
--- a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUWallFunction/nutUWallFunctionFvPatchScalarField.C
+++ b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUWallFunction/nutUWallFunctionFvPatchScalarField.C
@@ -6,7 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2011-2016 OpenFOAM Foundation
-    Copyright (C) 2019-2021 OpenCFD Ltd.
+    Copyright (C) 2019-2023 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -68,12 +68,11 @@ Foam::nutUWallFunctionFvPatchScalarField::calcNut() const
     forAll(yPlus, facei)
     {
         // Viscous sublayer contribution
-        const scalar nutVis = 0;
+        const scalar nutVis = nuw[facei];
 
         // Inertial sublayer contribution
         const scalar nutLog =
-            nuw[facei]
-           *(yPlus[facei]*kappa/log(max(E*yPlus[facei], 1 + 1e-4)) - 1.0);
+            nuw[facei]*yPlus[facei]*kappa/log(max(E*yPlus[facei], 1 + 1e-4));
 
         switch (blender_)
         {
@@ -132,6 +131,8 @@ Foam::nutUWallFunctionFvPatchScalarField::calcNut() const
                 break;
             }
         }
+
+        nutw[facei] -= nuw[facei];
     }
 
     return tnutw;
diff --git a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkWallFunction/nutkWallFunctionFvPatchScalarField.C b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkWallFunction/nutkWallFunctionFvPatchScalarField.C
index 8d9e40ea65dccc6dd24f33240b5636f152eacd17..d0d89ec71802002892788ceb9d4e7b93081c69b0 100644
--- a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkWallFunction/nutkWallFunctionFvPatchScalarField.C
+++ b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutkWallFunction/nutkWallFunctionFvPatchScalarField.C
@@ -6,7 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2011-2016, 2019 OpenFOAM Foundation
-    Copyright (C) 2019-2022 OpenCFD Ltd.
+    Copyright (C) 2019-2023 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -72,11 +72,11 @@ calcNut() const
         const scalar yPlus = Cmu25*y[facei]*sqrt(k[celli])/nuw[facei];
 
         // Viscous sublayer contribution
-        const scalar nutVis = 0;
+        const scalar nutVis = nuw[facei];
 
         // Inertial sublayer contribution
         const scalar nutLog =
-            nuw[facei]*(yPlus*kappa/log(max(E*yPlus, 1 + 1e-4)) - 1.0);
+            nuw[facei]*yPlus*kappa/log(max(E*yPlus, 1 + 1e-4));
 
         switch (blender_)
         {
@@ -134,6 +134,8 @@ calcNut() const
                 break;
             }
         }
+
+        nutw[facei] -= nuw[facei];
     }
 
     return tnutw;