From bc4824dcd960767a76ebe4319695f93888f0d38c Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Mon, 11 Jun 2012 14:27:23 +0100
Subject: [PATCH] ENH: Thermo film - added vapour recoil pressure effect

---
 .../surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C    | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C b/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C
index ecedee0485..1cb392eb6a 100644
--- a/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C
+++ b/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C
@@ -255,6 +255,9 @@ void thermoSingleLayer::updateSubmodels()
     // Update source fields
     hsSp_ += primaryEnergyPCTrans_/magSf()/time().deltaT();
     rhoSp_ += primaryMassPCTrans_/magSf()/time().deltaT();
+
+    // Vapour recoil pressure
+    pSp_ -= sqr(primaryMassPCTrans_/magSf()/time_.deltaT())/2.0/rhoPrimary_;
 }