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: improve offsets handling in extrudeToRegionMesh (#1933)

Browse Source 
- non-uniform offsets are generated due to truncation errors,
  which can lead to problems later on (e.g. redistributePar).

  Detect if the offsets are close to being uniform.
This commit is contained in:
mattijs
2020-11-23 19:58:31 +00:00
committed by Mark Olesen
parent 89ea011585
commit ae49b51182
1 changed files with 49 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

61
applications/utilities/mesh/generation/extrude/extrudeToRegionMesh/extrudeToRegionMesh.C
View File

@ -1134,18 +1134,55 @@ void setCouplingInfo
if (isA<mappedWallPolyPatch>(pp))
{
newPatches[patchi] = new mappedWallPolyPatch
(
pp.name(),
pp.size(),
pp.start(),
patchi,
sampleRegion, // sampleRegion
mode, // sampleMode
pp.name(), // samplePatch
offsets[zoneI], // offset
patches
);
const boundBox bb(pp.points(), pp.meshPoints(), true);
const vector avgOffset = gAverage(offsets[zoneI]);
const scalar mergeSqrDist =
gMax(magSqr(offsets[zoneI]-avgOffset));
// Verify uniformity of offset
// (same check as blockMesh geom merge)
if (mergeSqrDist < magSqr(10*SMALL*bb.span()))
{
Info<< "Adding on " << mesh.name()
<< " coupling patch " << pp.name()
<< " with uniform offset " << avgOffset << endl;
// Uniform offset
newPatches[patchi] = new mappedWallPolyPatch
(
pp.name(),
pp.size(),
pp.start(),
patchi,
sampleRegion, // sampleRegion
mode, // sampleMode
pp.name(), // samplePatch
avgOffset, // uniform offset
patches
);
}
else
{
Info<< "Adding on " << mesh.name()
<< " coupling patch " << pp.name()
<< " with non-uniform offset" << endl;
// Uniform offset
newPatches[patchi] = new mappedWallPolyPatch
(
pp.name(),
pp.size(),
pp.start(),
patchi,
sampleRegion, // sampleRegion
mode, // sampleMode
pp.name(), // samplePatch
offsets[zoneI], // non-uniform offsets
patches
);
}
}
}
}