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

Tow-stage Immersed face tolerance consistency. QUick reject using shortest face in edge and detailed analysis edge-by-edge. Inno Gatin

Browse source
This commit is contained in:
Hrvoje Jasak 2019-05-15 13:02:15 +01:00
parent 56eb4231ed
commit b6c9f0c577
1 changed files with 10 additions and 9 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

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

@ -257,10 +257,11 @@ void Foam::ImmersedFace<Distance>::createSubfaces
// than 3 points // than 3 points
if (nDry < 3) if (nDry < 3)
{ {
FatalErrorInFunction // The face is wet
<< "There are fewer than three points" isAllWet_ = true;
<< " on the wet part of the face." isAllDry_ = false;
<< abort(FatalError);
drySubface_.clear();
} }
else else
{ {
@ -272,10 +273,11 @@ void Foam::ImmersedFace<Distance>::createSubfaces
if (nWet < 3) if (nWet < 3)
{ {
FatalErrorInFunction // The face is dry
<< "There are fewer than three points" isAllWet_ = false;
<< " on the dry part of the face." isAllDry_ = true;
<< abort(FatalError);
wetSubface_.clear();
} }
else else
{ {
@ -342,7 +344,6 @@ Foam::ImmersedFace<Distance>::ImmersedFace
absTol = minEdgeLength*immersedPoly::tolerance_(); absTol = minEdgeLength*immersedPoly::tolerance_();
} }
// Check if all points are wet or dry, using absolute tolerance // Check if all points are wet or dry, using absolute tolerance
if (max(depth) < absTol) if (max(depth) < absTol)
{ {