exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Bugfix: Area tolerances needed in thin face cut

Browse source
This commit is contained in:
Hrvoje Jasak 2018-09-06 18:32:35 +01:00
parent 7508fede0b
commit 58421f0713
1 changed files with 6 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
src/immersedBoundary/immersedBoundary/immersedPoly/ImmersedFace.C
View file

@ -274,7 +274,8 @@ void Foam::ImmersedFace<Distance>::createSubfaces
if
(
drySubface_.mag(facePointsAndIntersections_)
< immersedPoly::tolerance_()
< immersedPoly::tolerance_()*
localFace.mag(facePointsAndIntersections_)
)
{
// The face is practically wet
@ -297,10 +298,13 @@ void Foam::ImmersedFace<Distance>::createSubfaces
wetSubface_.setSize(nWet);
// Check area: if it is very small, reset the face
// Note: must use relative tolerance, ie divide with original face
// area magnitude. HJ, 30/Aug/2018
if
(
wetSubface_.mag(facePointsAndIntersections_)
< immersedPoly::tolerance_()
< immersedPoly::tolerance_()*
localFace.mag(facePointsAndIntersections_)
)
{
// The face is practically dry