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: trackedParticle: could get stuck on processor face

Browse Source
This commit is contained in:
mattijs
2013-12-20 15:57:56 +00:00
parent 6af4a60a53
commit 4d42057791
1 changed files with 20 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files

31
src/mesh/autoMesh/autoHexMesh/trackedParticle/trackedParticle.C
View File

@ -90,27 +90,36 @@ Foam::trackedParticle::trackedParticle
bool Foam::trackedParticle::move
(
trackingData& td,
const scalar trackedParticle
const scalar trackTime
)
{
td.switchProcessor = false;
td.keepParticle = true;
scalar tEnd = (1.0 - stepFraction())*trackedParticle;
scalar tEnd = (1.0 - stepFraction())*trackTime;
scalar dtMax = tEnd;
while (td.keepParticle && !td.switchProcessor && tEnd > SMALL)
if (tEnd <= SMALL)
{
// set the lagrangian time-step
scalar dt = min(dtMax, tEnd);
// Remove the particle
td.keepParticle = false;
}
else
{
td.keepParticle = true;
// mark visited cell with max level.
td.maxLevel()[cell()] = max(td.maxLevel()[cell()], level_);
while (td.keepParticle && !td.switchProcessor && tEnd > SMALL)
{
// set the lagrangian time-step
scalar dt = min(dtMax, tEnd);
dt *= trackToFace(end_, td);
// mark visited cell with max level.
td.maxLevel()[cell()] = max(td.maxLevel()[cell()], level_);
tEnd -= dt;
stepFraction() = 1.0 - tEnd/trackedParticle;
dt *= trackToFace(end_, td);
tEnd -= dt;
stepFraction() = 1.0 - tEnd/trackTime;
}
}
return td.keepParticle;