From 8e552bd2b4e9561a8bc9bbd7e85bd12ee1b2695f Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Mon, 14 Feb 2022 16:02:19 +0000
Subject: [PATCH] ENH: snappyHexMesh: parallel consistency. Fixes #2331.
Coupled boundary faces should behave as internal faces
---
.../snappyHexMesh/meshRefinement/meshRefinementBaffles.C | 9 ++++++++-
1 file changed, 8 insertions(+), 1 deletion(-)
diff --git a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C
index 8bc18ac8846..8638d0386d6 100644
--- a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C
+++ b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C
@@ -737,6 +737,13 @@ void Foam::meshRefinement::getZoneFaces
nBaffles.setSize(zoneIDs.size());
nBaffles = Zero;
+
+ //- Get per face whether it is internal or coupled
+ const bitSet isInternalOrCoupled
+ (
+ syncTools::getInternalOrCoupledFaces(mesh_)
+ );
+
forAll(zoneIDs, j)
{
label zoneI = zoneIDs[j];
@@ -756,7 +763,7 @@ void Foam::meshRefinement::getZoneFaces
forAll(fz, i)
{
label faceI = fz[i];
- if (mesh_.isInternalFace(faceI))
+ if (isInternalOrCoupled[faceI])
{
if (fz.flipMap()[i])
{
--
GitLab