increase stability of cfdemSolverRhoPimpleChem

analogous to 16560cfe4f
VSMALL is typically DBL_MIN (1e-37 or smaller) which is actually too
small here
This commit is contained in:
danielque
2022-03-01 10:13:52 +01:00
parent c96e0cc6fd
commit aa42bb9a28

View File

@ -60,13 +60,13 @@ tmp<fv::convectionScheme<scalar> > mvConvection
if (propagateInertSpecie) if (propagateInertSpecie)
{ {
if (inertIndex!=-1) Yt /= (1-Y[inertIndex] + VSMALL); if (inertIndex!=-1) Yt /= (1-Y[inertIndex] + ROOTVSMALL);
forAll(Y,i) forAll(Y,i)
{ {
if (i!=inertIndex) if (i!=inertIndex)
{ {
volScalarField& Yi = Y[i]; volScalarField& Yi = Y[i];
Yi = Yi/(Yt+VSMALL); Yi = Yi/(Yt+ROOTVSMALL);
} }
} }
} }