From ba717a55103ae8e5d61c5651362f4f8bc8fea74b Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Mon, 12 Nov 2012 11:43:48 +0000
Subject: [PATCH] ENH: Calculate film boiling temperature instead of reading
from user input
---
.../standardPhaseChange/standardPhaseChange.C | 6 ++++--
.../standardPhaseChange/standardPhaseChange.H | 5 +----
2 files changed, 5 insertions(+), 6 deletions(-)
diff --git a/src/regionModels/surfaceFilmModels/submodels/thermo/phaseChangeModel/standardPhaseChange/standardPhaseChange.C b/src/regionModels/surfaceFilmModels/submodels/thermo/phaseChangeModel/standardPhaseChange/standardPhaseChange.C
index e632712ff48..4a08755ee3e 100644
--- a/src/regionModels/surfaceFilmModels/submodels/thermo/phaseChangeModel/standardPhaseChange/standardPhaseChange.C
+++ b/src/regionModels/surfaceFilmModels/submodels/thermo/phaseChangeModel/standardPhaseChange/standardPhaseChange.C
@@ -77,7 +77,6 @@ standardPhaseChange::standardPhaseChange
)
:
phaseChangeModel(typeName, owner, dict),
- Tb_(readScalar(coeffs_.lookup("Tb"))),
deltaMin_(readScalar(coeffs_.lookup("deltaMin"))),
L_(readScalar(coeffs_.lookup("L"))),
TbFactor_(coeffs_.lookupOrDefault<scalar>("TbFactor", 1.1))
@@ -130,8 +129,11 @@ void standardPhaseChange::correctModel
// cell pressure [Pa]
const scalar pc = pInf[cellI];
+ // calculate the boiling temperature
+ const scalar Tb = liq.pvInvert(pc);
+
// local temperature - impose lower limit of 200 K for stability
- const scalar Tloc = min(TbFactor_*Tb_, max(200.0, T[cellI]));
+ const scalar Tloc = min(TbFactor_*Tb, max(200.0, T[cellI]));
// saturation pressure [Pa]
const scalar pSat = liq.pv(pc, Tloc);
diff --git a/src/regionModels/surfaceFilmModels/submodels/thermo/phaseChangeModel/standardPhaseChange/standardPhaseChange.H b/src/regionModels/surfaceFilmModels/submodels/thermo/phaseChangeModel/standardPhaseChange/standardPhaseChange.H
index 1fff3384ff8..ee1b21d6d0a 100644
--- a/src/regionModels/surfaceFilmModels/submodels/thermo/phaseChangeModel/standardPhaseChange/standardPhaseChange.H
+++ b/src/regionModels/surfaceFilmModels/submodels/thermo/phaseChangeModel/standardPhaseChange/standardPhaseChange.H
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -69,9 +69,6 @@ protected:
// Protected data
- //- Boiling temperature / [K]
- const scalar Tb_;
-
//- Minimum film height for model to be active
const scalar deltaMin_;
--
GitLab