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

Added wet surface check for dry face

Browse source
This commit is contained in:
Hrvoje Jasak 2019-03-15 10:47:56 +00:00
parent 85e0d4f65e
commit 7fc1abd4b7
1 changed files with 26 additions and 26 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

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

@ -272,19 +272,19 @@ void Foam::ImmersedFace<Distance>::createSubfaces
{
drySubface_.setSize(nDry);
// // Check area: if it is very small, reset the face
// if
// (
// drySubface_.mag(facePointsAndIntersections_)
// < immersedPoly::liveFactor_()*
// localFace.mag(facePointsAndIntersections_)
// )
// {
// // The face is practically wet
// isAllWet_ = true;
// Check area: if it is very small, reset the face
if
(
drySubface_.mag(facePointsAndIntersections_)
< immersedPoly::liveFactor_()*
localFace.mag(facePointsAndIntersections_)
)
{
// The face is practically wet
isAllWet_ = true;
// drySubface_.clear();
// }
drySubface_.clear();
}
}
if (nWet < 3)
@ -299,21 +299,21 @@ 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::liveFactor_()*
// localFace.mag(facePointsAndIntersections_)
// )
// {
// // The face is practically dry
// isAllDry_ = true;
// 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::liveFactor_()*
localFace.mag(facePointsAndIntersections_)
)
{
// The face is practically dry
isAllDry_ = true;
// wetSubface_.clear();
// }
wetSubface_.clear();
}
}
}
}