wallFunctions: Handle moving meshes

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

@@ -96,11 +96,6 @@ void epsilonWallFunctionFvPatchScalarField::setMaster()
 
 void epsilonWallFunctionFvPatchScalarField::createAveragingWeights()
 {
-    if (initialised_)
-    {
-        return;
-    }
-
     const volScalarField& epsilon =
         static_cast<const volScalarField&>(this->dimensionedInternalField());
 
@@ -108,6 +103,11 @@ void epsilonWallFunctionFvPatchScalarField::createAveragingWeights()
 
     const fvMesh& mesh = epsilon.mesh();
 
+    if (initialised_ && !mesh.changing())
+    {
+        return;
+    }
+
     volScalarField weights
     (
         IOobject

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

@@ -98,11 +98,6 @@ void omegaWallFunctionFvPatchScalarField::setMaster()
 
 void omegaWallFunctionFvPatchScalarField::createAveragingWeights()
 {
-    if (initialised_)
-    {
-        return;
-    }
-
     const volScalarField& omega =
         static_cast<const volScalarField&>(this->dimensionedInternalField());
 
@@ -110,6 +105,11 @@ void omegaWallFunctionFvPatchScalarField::createAveragingWeights()
 
     const fvMesh& mesh = omega.mesh();
 
+    if (initialised_ && !mesh.changing())
+    {
+        return;
+    }
+
     volScalarField weights
     (
         IOobject

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

@@ -98,11 +98,6 @@ void epsilonWallFunctionFvPatchScalarField::setMaster()
 
 void epsilonWallFunctionFvPatchScalarField::createAveragingWeights()
 {
-    if (initialised_)
-    {
-        return;
-    }
-
     const volScalarField& epsilon =
         static_cast<const volScalarField&>(this->dimensionedInternalField());
 
@@ -110,6 +105,11 @@ void epsilonWallFunctionFvPatchScalarField::createAveragingWeights()
 
     const fvMesh& mesh = epsilon.mesh();
 
+    if (initialised_ && !mesh.changing())
+    {
+        return;
+    }
+
     volScalarField weights
     (
         IOobject

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

@@ -100,11 +100,6 @@ void omegaWallFunctionFvPatchScalarField::setMaster()
 
 void omegaWallFunctionFvPatchScalarField::createAveragingWeights()
 {
-    if (initialised_)
-    {
-        return;
-    }
-
     const volScalarField& omega =
         static_cast<const volScalarField&>(this->dimensionedInternalField());
 
@@ -112,6 +107,11 @@ void omegaWallFunctionFvPatchScalarField::createAveragingWeights()
 
     const fvMesh& mesh = omega.mesh();
 
+    if (initialised_ && !mesh.changing())
+    {
+        return;
+    }
+
     volScalarField weights
     (
         IOobject

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

@@ -98,11 +98,6 @@ void epsilonWallFunctionFvPatchScalarField::setMaster()
 
 void epsilonWallFunctionFvPatchScalarField::createAveragingWeights()
 {
-    if (initialised_)
-    {
-        return;
-    }
-
     const volScalarField& epsilon =
         static_cast<const volScalarField&>(this->dimensionedInternalField());
 
@@ -110,6 +105,11 @@ void epsilonWallFunctionFvPatchScalarField::createAveragingWeights()
 
     const fvMesh& mesh = epsilon.mesh();
 
+    if (initialised_ && !mesh.changing())
+    {
+        return;
+    }
+
     volScalarField weights
     (
         IOobject

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

@@ -100,11 +100,6 @@ void omegaWallFunctionFvPatchScalarField::setMaster()
 
 void omegaWallFunctionFvPatchScalarField::createAveragingWeights()
 {
-    if (initialised_)
-    {
-        return;
-    }
-
     const volScalarField& omega =
         static_cast<const volScalarField&>(this->dimensionedInternalField());
 
@@ -112,6 +107,11 @@ void omegaWallFunctionFvPatchScalarField::createAveragingWeights()
 
     const fvMesh& mesh = omega.mesh();
 
+    if (initialised_ && !mesh.changing())
+    {
+        return;
+    }
+
     volScalarField weights
     (
         IOobject