--- LiquidEvaporationBoil.C.org	2013-08-07 13:46:54.780906438 +0300
+++ LiquidEvaporationBoil.C	2013-08-07 13:48:06.294096459 +0300
@@ -151,18 +151,37 @@ void Foam::LiquidEvaporationBoil<CloudTy
     scalarField& dMassPC
 ) const
 {
-    // construct carrier phase species volume fractions for cell, cellI
-    const scalarField XcMix(calcXc(cellI));
-
     // liquid volume fraction
     const scalarField X(liquids_.X(Yl));
 
+    // immediately evaporate mass that has reached critical condition
+    if (mag(T - liquids_.Tc(X)) < SMALL)
+    {
+        if (debug)
+        {
+            WarningIn
+            (
+                "void Foam::LiquidEvaporationBoil<CloudType>::calculate"
+            )   << "Parcel reached critical conditions" << endl;
+        }
+
+        forAll(activeLiquids_, i)
+        {
+            const label lid = liqToLiqMap_[i];
+            dMassPC[lid] = GREAT;
+        }
+        return;
+    }
+
     // droplet surface pressure assumed to surface vapour pressure
     scalar ps = liquids_.pv(pc, Ts, X);
 
     // vapour density at droplet surface [kg/m3]
     scalar rhos = ps*liquids_.W(X)/(specie::RR*Ts);
 
+    // construct carrier phase species volume fractions for cell, cellI
+    const scalarField XcMix(calcXc(cellI));
+
     // carrier thermo properties
     scalar Hsc = 0.0;
     scalar Hc = 0.0;

From e215f19926038940914c1cf4308ae63494aa6ae9 Mon Sep 17 00:00:00 2001
From: Armin Wehrfritz <armin.wehrfritz@aalto.fi>
Date: Wed, 14 Aug 2013 12:21:16 +0300
Subject: [PATCH] Revise check for crit mass

---
 .../PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C   |    2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C
index d78d6b1..adde074 100644
--- a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C
+++ b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C
@@ -155,7 +155,7 @@ void Foam::LiquidEvaporationBoil<CloudType>::calculate
     const scalarField X(liquids_.X(Yl));
 
     // immediately evaporate mass that has reached critical condition
-    if (mag(T - liquids_.Tc(X)) < SMALL)
+    if ((liquids_.Tc(X) - T) < SMALL)
     {
         if (debug)
         {
-- 
1.7.10.4