From c9c3c4626953b37bca16bb44c96d05595d07401e Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Tue, 6 Aug 2013 11:08:50 +0100
Subject: [PATCH] BUG: SprayParcel - limit T to Tmax - mantis #942

---
 .../Templates/SprayParcel/SprayParcel.C       | 23 ++++---------------
 1 file changed, 5 insertions(+), 18 deletions(-)

diff --git a/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C b/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
index d0d3814bc7..876fa675a1 100644
--- a/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
+++ b/src/lagrangian/spray/parcels/Templates/SprayParcel/SprayParcel.C
@@ -76,6 +76,10 @@ void Foam::SprayParcel<ParcelType>::calc
         }
     }
 
+    // set the maximum temperature limit
+    const scalar TMax = td.cloud().composition().liquids().TMax(this->pc_);
+    td.cloud().constProps().TMax() = TMax;
+
     // store the parcel properties
     const scalarField& Y(this->Y());
     scalarField X(td.cloud().composition().liquids().X(Y));
@@ -301,15 +305,12 @@ Foam::scalar Foam::SprayParcel<ParcelType>::chi
 {
     // modifications to take account of the flash boiling on primary break-up
 
-    static label nIter = 200;
-
     typedef typename TrackData::cloudType::reactingCloudType reactingCloudType;
     const CompositionModel<reactingCloudType>& composition =
         td.cloud().composition();
 
     scalar chi = 0.0;
     scalar T0 = this->T();
-    scalar Tc0 = this->Tc();
     scalar p0 = this->pc();
     scalar pAmb = td.cloud().pAmbient();
 
@@ -322,21 +323,7 @@ Foam::scalar Foam::SprayParcel<ParcelType>::chi
             // liquid is boiling - calc boiling temperature
 
             const liquidProperties& liq = composition.liquids().properties()[i];
-            scalar TBoil = T0;
-
-            for (label k=0; k<nIter; k++)
-            {
-                scalar pBoil = liq.pv(p0, TBoil);
-
-                if (pBoil > p0)
-                {
-                    TBoil = TBoil - (T0 - Tc0)/nIter;
-                }
-                else
-                {
-                    break;
-                }
-            }
+            scalar TBoil = liq.pvInvert(p0);
 
             scalar hl = liq.hl(pAmb, TBoil);
             scalar iTp = liq.h(pAmb, T0) - liq.rho(pAmb, T0);