Enable turbulent walls for coupled patches

src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/regionCouple/regionCouplePolyPatch.C

@@ -172,13 +172,15 @@ Foam::regionCouplePolyPatch::regionCouplePolyPatch
     const polyBoundaryMesh& bm,
     const word& shadowRegionName,
     const word& shadowPatchName,
-    const bool attached
+    const bool attached,
+    const bool attachedWalls
 )
 :
     coupledPolyPatch(name, size, start, index, bm),
     shadowRegionName_(shadowRegionName),
     shadowPatchName_(shadowPatchName),
     attached_(attached),
+    attachedWalls_(attachedWalls),
     shadowIndex_(-1),
     patchToPatchPtr_(NULL),
     reconFaceCellCentresPtr_(NULL)
@@ -198,6 +200,7 @@ Foam::regionCouplePolyPatch::regionCouplePolyPatch
     shadowRegionName_(dict.lookup("shadowRegion")),
     shadowPatchName_(dict.lookup("shadowPatch")),
     attached_(dict.lookup("attached")),
+    attachedWalls_(dict.lookup("attachedWalls")),
     shadowIndex_(-1),
     patchToPatchPtr_(NULL),
     reconFaceCellCentresPtr_(NULL)
@@ -215,6 +218,7 @@ Foam::regionCouplePolyPatch::regionCouplePolyPatch
     shadowRegionName_(pp.shadowRegionName_),
     shadowPatchName_(pp.shadowPatchName_),
     attached_(pp.attached_),
+    attachedWalls_(pp.attachedWalls_),
     shadowIndex_(-1),
     patchToPatchPtr_(NULL),
     reconFaceCellCentresPtr_(NULL)
@@ -235,6 +239,7 @@ Foam::regionCouplePolyPatch::regionCouplePolyPatch
     shadowRegionName_(pp.shadowRegionName_),
     shadowPatchName_(pp.shadowPatchName_),
     attached_(pp.attached_),
+    attachedWalls_(pp.attachedWalls_),
     shadowIndex_(-1),
     patchToPatchPtr_(NULL),
     reconFaceCellCentresPtr_(NULL)

src/OpenFOAM/meshes/polyMesh/polyPatches/constraint/regionCouple/regionCouplePolyPatch.H

@@ -69,6 +69,9 @@ class regionCouplePolyPatch
         //- Are the regions attached
         mutable Switch attached_;
 
+        //- Are the region attached walls
+        mutable Switch attachedWalls_;
+
         //- Shadow patch index.  Delayed evaluation for construction
         mutable label shadowIndex_;
 
@@ -85,7 +88,7 @@ class regionCouplePolyPatch
         //- Calculate interpolation
         void calcInterpolation() const;
 
-
+        
 protected:
 
     // Protected Member functions
@@ -139,7 +142,8 @@ public:
             const polyBoundaryMesh& bm,
             const word& shadowRegionName,
             const word& shadowPatchName,
-            const bool attached
+            const bool attached,
+            const bool attachedWall
         );
 
         //- Construct from dictionary
@@ -233,6 +237,12 @@ public:
             return attached_;
         }
 
+        //AJ: read from dictionary if coupled patch is also a wall
+        bool isWall() const
+        {
+            return attachedWalls_;
+        }
+
         //- Attach regions
         void attach() const;
 

src/OpenFOAM/meshes/polyMesh/polyPatches/derived/wall/wallPolyPatch.H

@@ -114,6 +114,11 @@ public:
             );
         }
 
+        // Virtual function for wall handling of derived class
+        virtual bool isWall() const
+        {
+	    return true;
+        }
 
     // Destructor
 

src/OpenFOAM/meshes/polyMesh/polyPatches/polyPatch/polyPatch.H

@@ -368,6 +368,12 @@ public:
             return false;
         }
 
+        // Virtual function for wall handling of all derived calsses
+        virtual bool isWall() const
+        {
+            return false;
+        }       
+
         //- Return true if the given type is a constraint type
         static bool constraintType(const word& pt);
 

src/finiteVolume/fvMesh/fvPatches/fvPatch/fvPatch.H

@@ -170,6 +170,12 @@ public:
                 return polyPatch_.coupled();
             }
 
+            //- Return true if this patch is wall
+            bool isWall() const
+            {
+                return polyPatch_.isWall();
+            }
+
             //- Return true if the given type is a constraint type
             static bool constraintType(const word& pt);
 
@@ -200,10 +206,7 @@ public:
             }
 
             template<class T>
-            const typename Field<T>::subField patchSlice
+            const typename Field<T>::subField patchSlice(const Field<T>& l) const
-            (
-                const Field<T>& l
-            ) const
             {
                 return typename Field<T>::subField
                 (

src/finiteVolume/fvMesh/wallDist/nearWallDist.C

@@ -36,8 +36,19 @@ void Foam::nearWallDist::doAll()
 {
     cellDistFuncs wallUtils(mesh_);
 
+    // AJ: make sure to pick up all patches that are specified as a wall
+    const polyBoundaryMesh& bMesh = wallUtils.mesh().boundaryMesh();
+    labelHashSet wallPatchIDs(bMesh.size());
+    forAll(bMesh, patchI)
+    {
+        if (bMesh[patchI].isWall())
+        {
+            wallPatchIDs.insert(patchI);
+        }
+    }
+
     // Get patch ids of walls
-    labelHashSet wallPatchIDs(wallUtils.getPatchIDs<wallPolyPatch>());
+    // labelHashSet wallPatchIDs(wallUtils.getPatchIDs<wallPolyPatch>());
 
     // Size neighbours array for maximum possible
 
@@ -54,7 +65,9 @@ void Foam::nearWallDist::doAll()
 
         const fvPatch& patch = mesh_.boundary()[patchI];
 
-        if (isA<wallFvPatch>(patch))
+        // AJ: Allow other patch types to be seen as a wall type
+        // if (isA<wallFvPatch>(patch))
+        if (patch.isWall())
         {
             const polyPatch& pPatch = patch.patch();
 

src/finiteVolume/fvMesh/wallDist/wallDist.C

@@ -67,13 +67,24 @@ Foam::wallDist::~wallDist()
 // future (if only small topology change)
 void Foam::wallDist::correct()
 {
+    // AJ: make sure to pick up all patches that are specified as a wall
+    const polyBoundaryMesh& bMesh = cellDistFuncs::mesh().boundaryMesh();
+    labelHashSet wallPatchIDs(bMesh.size());
+    forAll(bMesh, patchI)
+    {
+        if (bMesh[patchI].isWall())
+        {
+            wallPatchIDs.insert(patchI);
+        }
+    }
+
     // Get patchids of walls
-    labelHashSet wallPatchIDs(getPatchIDs<wallPolyPatch>());
+    // labelHashSet wallPatchIDs(getPatchIDs<wallPolyPatch>());
 
     // Calculate distance starting from wallPatch faces.
     patchWave wave(cellDistFuncs::mesh(), wallPatchIDs, correctWalls_);
 
-    // Transfer cell values from wave into *this
+    // Transfer cell values from wave into *this 
     transfer(wave.distance());
 
     // Transfer values on patches into boundaryField of *this

src/finiteVolume/fvMesh/wallDist/wallDistData.C

@@ -73,14 +73,25 @@ Foam::wallDistData<TransferType>::~wallDistData()
 template<class TransferType>
 void Foam::wallDistData<TransferType>::correct()
 {
+  Info<< "wallDistData.correct() called" << endl;
     const polyMesh& mesh = cellDistFuncs::mesh();
 
     //
     // Fill data on wall patches with initial values
     //
 
+    const polyBoundaryMesh& bMesh = mesh.boundaryMesh();
+    labelHashSet wallPatchIDs(bMesh.size());
+    forAll(bMesh, patchI)
+    {
+        if (bMesh[patchI].isWall())
+        {
+            wallPatchIDs.insert(patchI);
+        }
+    }
+
     // Get patchids of walls
-    labelHashSet wallPatchIDs(getPatchIDs<wallPolyPatch>());
+    // labelHashSet wallPatchIDs(getPatchIDs<wallPolyPatch>());
 
     // Collect pointers to data on patches
     UPtrList<Field<Type> > patchData(mesh.boundaryMesh().size());

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

@@ -44,7 +44,9 @@ namespace RASModels
 
 void epsilonWallFunctionFvPatchScalarField::checkType()
 {
-    if (!isA<wallFvPatch>(patch()))
+    // AJ: Allow other patch types to be seen as wall type
+    // if (!isA<wallFvPatch>(patch()))
+    if (!this->patch().isWall())
     {
         FatalErrorIn("epsilonWallFunctionFvPatchScalarField::checkType()")
             << "Invalid wall function specification" << nl

src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/kqRWallFunctions/kqRWallFunction/kqRWallFunctionFvPatchField.C

@@ -28,6 +28,8 @@ License
 #include "fvPatchFieldMapper.H"
 #include "addToRunTimeSelectionTable.H"
 #include "wallFvPatch.H"
+#include "regionCoupleFvPatch.H"
+#include "movingWallVelocityFvPatchVectorField.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -43,7 +45,9 @@ namespace RASModels
 template<class Type>
 void kqRWallFunctionFvPatchField<Type>::checkType()
 {
-    if (!isA<wallFvPatch>(this->patch()))
+    // AJ: Allow other patch types to be seemovingWalln as wall type
+    // if (!isA<wallFvPatch>(this->patch())) 
+    if (!this->patch().isWall()) 
     {
         FatalErrorIn("kqRWallFunctionFvPatchField::checkType()")
             << "Invalid wall function specification" << nl