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

KocamustafaogullariIshiiDepartureFrequency: Stabilisation

Browse Source 
Prevent floating point exception if vapor density is higher than liquid
density

Patch contributed by Juho Peltola, VTT.
This commit is contained in:
Will Bainbridge
2019-08-12 09:24:26 +01:00
parent af2baeb6d1
commit 491ee7efcc
1 changed files with 1 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
applications/solvers/multiphase/reactingEulerFoam/derivedFvPatchFields/wallBoilingSubModels/departureFrequencyModels/KocamustafaogullariIshiiDepartureFrequency/KocamustafaogullariIshiiDepartureFrequency.C
View File

@ -90,7 +90,7 @@ KocamustafaogullariIshiiDepartureFrequency::fDeparture
liquid.mesh().lookupObject<uniformDimensionedVectorField>("g");
const scalarField rhoLiquid(liquid.thermo().rho(patchi));
const scalarField rhoVapor(vapor.thermo().rho(patchi));
const scalarField rhoVapor(min(vapor.thermo().rho(patchi), rhoLiquid));
const tmp<volScalarField> tsigma
(