Commit e8985cbd authored by andy's avatar andy
Browse files

BUG: Corrected commit 6a1b5673

parent 668cd823
......@@ -287,16 +287,13 @@ bool Foam::KinematicParcel<ParcelType>::move
const label cellI = p.cell();
const scalar magU = mag(U_);
if (p.active())
if (p.active() && moving && (magU > ROOTVSMALL))
{
const scalar d = dt*magU;
const scalar dCorr = min(d, maxCo*cbrt(V[cellI]));
dt *= dCorr/d;
if (moving && (magU > ROOTVSMALL))
{
dt *= p.trackToFace(p.position() + dCorr*U_/magU, td);
}
dt *=
dCorr/d
*p.trackToFace(p.position() + dCorr*U_/magU, td);
}
tEnd -= dt;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment