ENH: Turbulence wal functions - apply zero-gradient condition on start-up

src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C

@@ -311,6 +311,9 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
     cornerWeights_()
 {
     checkType();
+
+    // apply zero-gradient condition on start-up
+    this->operator==(patchInternalField());
 }
 
 

src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C

@@ -324,6 +324,9 @@ omegaWallFunctionFvPatchScalarField::omegaWallFunctionFvPatchScalarField
     cornerWeights_()
 {
     checkType();
+
+    // apply zero-gradient condition on start-up
+    this->operator==(patchInternalField());
 }
 
 

src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/epsilonWallFunctions/epsilonWallFunction/epsilonWallFunctionFvPatchScalarField.C

@@ -311,6 +311,9 @@ epsilonWallFunctionFvPatchScalarField::epsilonWallFunctionFvPatchScalarField
     cornerWeights_()
 {
     checkType();
+
+    // apply zero-gradient condition on start-up
+    this->operator==(patchInternalField());
 }
 
 

src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/omegaWallFunctions/omegaWallFunction/omegaWallFunctionFvPatchScalarField.C

@@ -322,6 +322,9 @@ omegaWallFunctionFvPatchScalarField::omegaWallFunctionFvPatchScalarField
     cornerWeights_()
 {
     checkType();
+
+    // apply zero-gradient condition on start-up
+    this->operator==(patchInternalField());
 }