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

BUG: Improving robustness of tracking on cells with detA zero

Browse Source 
1) Small modification to the tracking logic for detA zero.
2) Adding small vector displacement to locate function to avoid error where
particle is inserted at the cell centre.
This commit is contained in:
sergio
2021-02-15 12:20:02 -08:00
parent 8f4b080542
commit 48fdec46c7
1 changed files with 8 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
src/lagrangian/basic/particle/particle.C
View File

@ -432,7 +432,12 @@ void Foam::particle::locate
}
}
const vector displacement = position - mesh_.cellCentres()[celli_];
// Perturbing displacement to avoid zero in case position is the
// cellCentre
const vector displacement =
position
- mesh_.cellCentres()[celli_]
+ vector(VSMALL, VSMALL, VSMALL);
// Loop all cell tets to find the one containing the position. Track
// through each tet from the cell centre. If a tet contains the position
@ -748,7 +753,7 @@ Foam::scalar Foam::particle::trackToStationaryTri
// Calculate the hit fraction
label iH = -1;
scalar muH = std::isnormal(detA) && detA <= 0 ? VGREAT : 1/detA;
scalar muH = detA <= 0 ? VGREAT : 1/detA;
for (label i = 0; i < 4; ++ i)
{
if (Tx1[i] < - detA*SMALL)
@ -903,7 +908,7 @@ Foam::scalar Foam::particle::trackToMovingTri
// Calculate the hit fraction
label iH = -1;
scalar muH = std::isnormal(detA[0]) && detA[0] <= 0 ? VGREAT : 1/detA[0];
scalar muH = detA[0] <= 0 ? VGREAT : 1/detA[0];
for (label i = 0; i < 4; ++i)
{
const Roots<3> mu = hitEqn[i].roots();