zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

lagrangian/spray/submodels/BreakupModel/SHF: Prevent divide-by-zero

Browse Source 
Resolves bug report https://bugs.openfoam.org/view.php?id=3495
This commit is contained in:
Will Bainbridge
2020-05-22 09:06:22 +01:00
parent 8604c1eb5f
commit 9b23ae3b1b
1 changed files with 15 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
src/lagrangian/spray/submodels/BreakupModel/SHF/SHF.C
View File

@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2011-2019 OpenFOAM Foundation
\\ / A nd | Copyright (C) 2011-2020 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@ -156,9 +156,21 @@ bool Foam::SHF<CloudType>::update
scalar weGasCorr = weGas/(1.0 + weCorrCoeff_*ohnesorge);
// droplet deformation characteristic time
// update the droplet characteristic time
tc += dt;
scalar tChar = d/Urmag*sqrt(rho/rhoc);
// droplet deformation characteristic rate
scalar rChar = Urmag/d*sqrt(rhoc/rho);
// return if the characteristic deformation rate is too low for the
// following modelling to be calculable
if (tc*rChar < small)
{
return false;
}
// droplet deformation characteristic time
scalar tChar = 1/rChar;
scalar tFirst = cInit_*tChar;
@ -171,9 +183,6 @@ bool Foam::SHF<CloudType>::update
bool success = false;
// update the droplet characteristic time
tc += dt;
if (weGas > weConst_)
{
if (weGas < weCrit1_)