From 47c6b9e12e70e0fb73d37cbc634c321d7d849811 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Thu, 4 Aug 2016 22:25:27 +0100
Subject: [PATCH] epsilonWallFunction, omegaWallFunction: Set generation to
 zero in laminar sub-layer

Improves stability for complex flows
---
 .../epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C | 2 --
 .../omegaWallFunction/omegaWallFunctionFvPatchScalarField.C     | 2 --
 2 files changed, 4 deletions(-)

diff --git a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C
index 0334befa91..cda78252e6 100644
--- a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C
+++ b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C
@@ -255,8 +255,6 @@ void Foam::epsilonWallFunctionFvPatchScalarField::calculate
         else
         {
             epsilon0[celli] += w*2.0*k[celli]*nuw[facei]/sqr(y[facei]);
-
-            G0[celli] += w*G[celli];
         }
     }
 }
diff --git a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C
index 2c2af10947..b8633f62bf 100644
--- a/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C
+++ b/src/TurbulenceModels/turbulenceModels/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C
@@ -280,8 +280,6 @@ void omegaWallFunctionFvPatchScalarField::calculate
             {
                 omega0[celli] += w*omegaVis;
             }
-
-            G0[celli] += w*G[celli];
         }
     }
 }