diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C index 22e01b9a98ced336dbbfb5e0745ed8cc93034a69..eee3b7e8ecba6536f809a8c11fe64ee56d84c0ed 100644 --- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C +++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -28,6 +28,12 @@ License #include "IntegrationScheme.H" #include "meshTools.H" +// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * // + +template<class ParcelType> +Foam::label Foam::KinematicParcel<ParcelType>::maxTrackAttempts = 1; + + // * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * * // template<class ParcelType> @@ -271,7 +277,8 @@ bool Foam::KinematicParcel<ParcelType>::move scalar tEnd = (1.0 - p.stepFraction())*trackTime; const scalar dtMax = tEnd; - bool moving = true; + bool tracking = true; + label nTrackingStalled = 0; while (td.keepParticle && !td.switchProcessor && tEnd > ROOTVSMALL) { @@ -287,7 +294,7 @@ bool Foam::KinematicParcel<ParcelType>::move const label cellI = p.cell(); const scalar magU = mag(U_); - if (p.active() && moving && (magU > ROOTVSMALL)) + if (p.active() && tracking && (magU > ROOTVSMALL)) { const scalar d = dt*magU; const scalar dCorr = min(d, maxCo*cellLengthScale[cellI]); @@ -300,13 +307,25 @@ bool Foam::KinematicParcel<ParcelType>::move scalar newStepFraction = 1.0 - tEnd/trackTime; - if - ( - mag(p.stepFraction() - newStepFraction) - < particle::minStepFractionTol - ) + if (tracking) { - moving = false; + if + ( + mag(p.stepFraction() - newStepFraction) + < particle::minStepFractionTol + ) + { + nTrackingStalled++; + + if (nTrackingStalled > maxTrackAttempts) + { + tracking = false; + } + } + else + { + nTrackingStalled = 0; + } } p.stepFraction() = newStepFraction; diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H index 63d15118452296ea86d92c22af0b05bc1c831499..2f42d4c02e840a5538316de1cca05255636c0d5e 100644 --- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H +++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -218,6 +218,10 @@ public: }; + //- Number of particle tracking attempts before we assume that it stalls + static label maxTrackAttempts; + + protected: // Protected data