From 658485a60f6763f3a41119ae8ed95bbac3e97b7e Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Mon, 16 Jul 2012 10:24:11 +0100
Subject: [PATCH] BUG: corrected lagrangian evaporation models after changes to
he thermo
---
.../PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C | 2 +-
.../LiquidEvaporationBoil/LiquidEvaporationBoil.C | 6 +++---
2 files changed, 4 insertions(+), 4 deletions(-)
diff --git a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C
index 475430bd83e..97b84a6049f 100644
--- a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C
+++ b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C
@@ -216,7 +216,7 @@ Foam::scalar Foam::LiquidEvaporation<CloudType>::dh
}
case (parent::etEnthalpyDifference):
{
- scalar hc = this->owner().composition().carrier().Hs(idc, p, T);
+ scalar hc = this->owner().composition().carrier().Ha(idc, p, T);
scalar hp = liquids_.properties()[idl].h(p, T);
dh = hc - hp;
diff --git a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C
index 32821f89662..803b849fa43 100644
--- a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C
+++ b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C
@@ -171,8 +171,8 @@ void Foam::LiquidEvaporationBoil<CloudType>::calculate
forAll(this->owner().thermo().carrier().Y(), i)
{
scalar Yc = this->owner().thermo().carrier().Y()[i][cellI];
- Hc += Yc*this->owner().thermo().carrier().Hs(i, pc, Tc);
- Hsc += Yc*this->owner().thermo().carrier().Hs(i, ps, Ts);
+ Hc += Yc*this->owner().thermo().carrier().Ha(i, pc, Tc);
+ Hsc += Yc*this->owner().thermo().carrier().Ha(i, ps, Ts);
Cpc += Yc*this->owner().thermo().carrier().Cp(i, ps, Ts);
kappac += Yc*this->owner().thermo().carrier().kappa(i, ps, Ts);
}
@@ -315,7 +315,7 @@ Foam::scalar Foam::LiquidEvaporationBoil<CloudType>::dh
}
case (parent::etEnthalpyDifference):
{
- scalar hc = this->owner().composition().carrier().Hs(idc, p, TDash);
+ scalar hc = this->owner().composition().carrier().Ha(idc, p, TDash);
scalar hp = liquids_.properties()[idl].h(p, TDash);
dh = hc - hp;
--
GitLab