research/eReactingFoam-4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

added testing infinity for pure condition

Browse source
This commit is contained in:
Yeongdo Park 2019-01-11 03:55:15 -05:00
parent 0da17683df
commit 93500dd389
1 changed files with 4 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

5
diffusivityModel/diffusivityModel/diffusivityModel.C
View file

@ -25,6 +25,8 @@ License
#include "diffusivityModel.H"
#include <cmath>
#include "IFstream.H"
#include "speciesTable.H"
#include "volFieldsFwd.H"
@ -74,7 +76,8 @@ Foam::diffusivityModel::mixAvgDi(UList<scalar>& Di, const scalarSymmetricSquareM
label pureSpecieI = -1;
forAll (X, i)
{
if (1. - X[i] < 1e-12)
if (1. - X[i] < 1e-12 || std::isinf(X[i]/(1.-Y[i])))
{
pureSpecieI = i;
pure = true;