added check for dt > ROOTVSMALL before allowing to calc (un)coupled

This commit is contained in:
andy
2008-08-19 18:33:39 +01:00
parent 53133e9c86
commit 58f43193da

View File

@ -234,13 +234,17 @@ bool Foam::KinematicParcel<ParcelType>::move
// Update cell based properties
p.updateCellQuantities(td, dt, celli);
if (td.cloud().coupled())
// Avoid problems with extremely small timesteps
if (dt > ROOTVSMALL)
{
p.calcCoupled(td, dt, celli);
}
else
{
p.calcUncoupled(td, dt, celli);
if (td.cloud().coupled())
{
p.calcCoupled(td, dt, celli);
}
else
{
p.calcUncoupled(td, dt, celli);
}
}
if (p.onBoundary() && td.keepParticle)