zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: snappyHexMesh: more parallel consistency. See #2331

Browse Source
This commit is contained in:
mattijs
2024-04-24 17:10:27 +01:00
parent 31b74cd627
commit 690d02f9e8
1 changed files with 17 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

27
src/mesh/snappyHexMesh/meshRefinement/meshRefinementProblemCells.C
View File

@ -6,7 +6,7 @@
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
Copyright (C) 2011-2015 OpenFOAM Foundation Copyright (C) 2011-2015 OpenFOAM Foundation
Copyright (C) 2015-2022 OpenCFD Ltd. Copyright (C) 2015-2022,2024 OpenCFD Ltd.
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
This file is part of OpenFOAM. This file is part of OpenFOAM.
@ -473,13 +473,16 @@ void Foam::meshRefinement::markFacesOnProblemCells
// Baffle all faces of cells that need to be removed // Baffle all faces of cells that need to be removed
forAllConstIters(problemCells, iter) forAllConstIters(problemCells, iter)
{ {
const cell& cFaces = mesh_.cells()[iter.key()]; for (const label facei : mesh_.cells()[iter.key()])
forAll(cFaces, i)
{ {
label facei = cFaces[i]; const label patchi = patches.whichPatch(facei);
if (facePatch[facei] == -1 && mesh_.isInternalFace(facei)) if
(
facePatch[facei] == -1
// && mesh_.isInternalFace(facei)
&& (patchi == -1 || patches[patchi].coupled())
)
{ {
facePatch[facei] = nearestAdaptPatch[facei]; facePatch[facei] = nearestAdaptPatch[facei];
faceZone[facei] = nearestAdaptZone[facei]; faceZone[facei] = nearestAdaptZone[facei];
@ -719,14 +722,15 @@ void Foam::meshRefinement::markFacesOnProblemCells
else else
{ {
// Block all faces of cell // Block all faces of cell
forAll(cFaces, cf) for (const label facei : cFaces)
{ {
label facei = cFaces[cf]; const label patchi = patches.whichPatch(facei);
if if
( (
facePatch[facei] == -1 facePatch[facei] == -1
&& mesh_.isInternalFace(facei) // && mesh_.isInternalFace(facei)
&& (patchi == -1 || patches[patchi].coupled())
) )
{ {
facePatch[facei] = nearestAdaptPatch[facei]; facePatch[facei] = nearestAdaptPatch[facei];
@ -805,10 +809,13 @@ void Foam::meshRefinement::markFacesOnProblemCells
for (const label facei : cFaces) for (const label facei : cFaces)
{ {
const label patchi = patches.whichPatch(facei);
if if
( (
facePatch[facei] == -1 facePatch[facei] == -1
&& mesh_.isInternalFace(facei) // && mesh_.isInternalFace(facei)
&& (patchi == -1 || patches[patchi].coupled())
) )
{ {
facePatch[facei] = nearestAdaptPatch[facei]; facePatch[facei] = nearestAdaptPatch[facei];