BUG: mantis #1408: added small value protection to division

This commit is contained in:
william
2014-10-15 11:50:55 +01:00
committed by Andrew Heather
parent 02d1f5d6f0
commit 2dfb664e1f

View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2012 OpenFOAM Foundation \\ / A nd | Copyright (C) 2012-2014 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -54,7 +54,7 @@ Foam::scalar Foam::SaffmanMeiLiftForce<CloudType>::SaffmanMeiLiftForce::Cl
Cld = 6.46*0.0524*sqrt(beta*Re); 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;
} }