From 491ee7efcc40b3de596152b3999acd74aaefebb7 Mon Sep 17 00:00:00 2001
From: Will Bainbridge <http://cfd.direct>
Date: Mon, 12 Aug 2019 09:24:26 +0100
Subject: [PATCH] KocamustafaogullariIshiiDepartureFrequency: Stabilisation

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

Patch contributed by Juho Peltola, VTT.
---
 .../KocamustafaogullariIshiiDepartureFrequency.C                | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/applications/solvers/multiphase/reactingEulerFoam/derivedFvPatchFields/wallBoilingSubModels/departureFrequencyModels/KocamustafaogullariIshiiDepartureFrequency/KocamustafaogullariIshiiDepartureFrequency.C b/applications/solvers/multiphase/reactingEulerFoam/derivedFvPatchFields/wallBoilingSubModels/departureFrequencyModels/KocamustafaogullariIshiiDepartureFrequency/KocamustafaogullariIshiiDepartureFrequency.C
index b23a89ca73..a977897e19 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/derivedFvPatchFields/wallBoilingSubModels/departureFrequencyModels/KocamustafaogullariIshiiDepartureFrequency/KocamustafaogullariIshiiDepartureFrequency.C
+++ b/applications/solvers/multiphase/reactingEulerFoam/derivedFvPatchFields/wallBoilingSubModels/departureFrequencyModels/KocamustafaogullariIshiiDepartureFrequency/KocamustafaogullariIshiiDepartureFrequency.C
@@ -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
     (