From 2dfb664e1f6c3871704b6a39bd331c0523aabc61 Mon Sep 17 00:00:00 2001
From: william <william>
Date: Wed, 15 Oct 2014 11:50:55 +0100
Subject: [PATCH] BUG: mantis #1408: added small value protection to division

---
 .../ParticleForces/Lift/SaffmanMeiLift/SaffmanMeiLiftForce.C  | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/Lift/SaffmanMeiLift/SaffmanMeiLiftForce.C b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/Lift/SaffmanMeiLift/SaffmanMeiLiftForce.C
index d6d8dda976..26da78568f 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/Lift/SaffmanMeiLift/SaffmanMeiLiftForce.C
+++ b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/Lift/SaffmanMeiLift/SaffmanMeiLiftForce.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2012 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2012-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -54,7 +54,7 @@ Foam::scalar Foam::SaffmanMeiLiftForce<CloudType>::SaffmanMeiLiftForce::Cl
         Cld = 6.46*0.0524*sqrt(beta*Re);
     }
 
-    return 3.0/(mathematical::twoPi*sqrt(Rew))*Cld;
+    return 3.0/(mathematical::twoPi*sqrt(Rew + ROOTVSMALL))*Cld;
 }