From ad2c46f7d89403abf21019f0deb13a3df49ad8f0 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Fri, 9 Aug 2013 12:37:36 +0100
Subject: [PATCH] ENH: PhaseChangeModel - updated TMax function and added Tvap
 function

---
 .../PhaseChangeModel/PhaseChangeModel.C             | 13 ++++++++++++-
 .../PhaseChangeModel/PhaseChangeModel.H             |  5 ++++-
 2 files changed, 16 insertions(+), 2 deletions(-)

diff --git a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/PhaseChangeModel/PhaseChangeModel.C b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/PhaseChangeModel/PhaseChangeModel.C
index 3ba13f27e1..9ec9307278 100644
--- a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/PhaseChangeModel/PhaseChangeModel.C
+++ b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/PhaseChangeModel/PhaseChangeModel.C
@@ -179,12 +179,23 @@ Foam::scalar Foam::PhaseChangeModel<CloudType>::dh
 
 
 template<class CloudType>
-Foam::scalar Foam::PhaseChangeModel<CloudType>::TMax(const scalar) const
+Foam::scalar Foam::PhaseChangeModel<CloudType>::TMax
+(
+    const scalar,
+    const scalarField&
+) const
 {
     return GREAT;
 }
 
 
+template<class CloudType>
+Foam::scalar Foam::PhaseChangeModel<CloudType>::Tvap(const scalarField& Y) const
+{
+    return -GREAT;
+}
+
+
 template<class CloudType>
 void Foam::PhaseChangeModel<CloudType>::addToPhaseChangeMass(const scalar dMass)
 {
diff --git a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/PhaseChangeModel/PhaseChangeModel.H b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/PhaseChangeModel/PhaseChangeModel.H
index f61031c758..6f3a171708 100644
--- a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/PhaseChangeModel/PhaseChangeModel.H
+++ b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/PhaseChangeModel/PhaseChangeModel.H
@@ -184,8 +184,11 @@ public:
             const scalar T
         ) const;
 
+        //- Return vapourisation temperature
+        virtual scalar Tvap(const scalarField& Y) const;
+
         //- Return maximum/limiting temperature
-        virtual scalar TMax(const scalar pIn) const;
+        virtual scalar TMax(const scalar p, const scalarField& Y) const;
 
         //- Add to phase change mass
         void addToPhaseChangeMass(const scalar dMass);