From e2b6e3238b9336fae48fb1551d083bc1968838de Mon Sep 17 00:00:00 2001
From: andy <a.heather@opencfd.co.uk>
Date: Thu, 13 Jan 2011 11:29:36 +0000
Subject: [PATCH] ENH: Limit tracking distance for steady calcs

---
 .../Templates/KinematicParcel/KinematicParcel.C        | 10 +++++++---
 1 file changed, 7 insertions(+), 3 deletions(-)

diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C
index cccb87fea1..81ba4c1023 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) 1991-2010 OpenCFD Ltd.
+    \\  /    A nd           | Copyright (C) 1991-2011 OpenCFD Ltd.
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -340,9 +340,13 @@ bool Foam::KinematicParcel<ParcelType>::move
                 // will change if a face is hit
                 label cellI = p.cell();
 
-                if (p.active())
+                const scalar magU = mag(U_);
+                if (p.active() && magU > ROOTVSMALL)
                 {
-                    dt *= p.trackToFace(p.position() + dt*U_, td);
+                    const scalar d = mag(dt*U_) + ROOTVSMALL;
+                    const vector n = U_/magU;
+                    const scalar dCorr = min(d, mag(n & mesh.bounds().span()));
+                    dt *= dCorr/d*p.trackToFace(p.position() + dCorr*n, td);
                 }
 
                 tEnd -= dt;