From f989d8e181f3bef197a04cf1dcaf2f28f14d643b 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

---
 .../parcels/Templates/KinematicParcel/KinematicParcel.C    | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.C
index cccb87fea1..1be59d9d06 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
@@ -342,7 +342,10 @@ bool Foam::KinematicParcel<ParcelType>::move
 
                 if (p.active())
                 {
-                    dt *= p.trackToFace(p.position() + dt*U_, td);
+                    scalar d = mag(dt*U_);
+                    const vector dir = U_/mag(U_);
+                    d = min(d, dir & mesh.bounds().span());
+                    dt *= p.trackToFace(p.position() + d*dir, td);
                 }
 
                 tEnd -= dt;