zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: corrected lagrangian evaporation models after changes to he thermo

Browse Source
This commit is contained in:
andy
2012-07-16 10:24:11 +01:00
parent 734f853ca9
commit 658485a60f
2 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C
View File

@ -216,7 +216,7 @@ Foam::scalar Foam::LiquidEvaporation<CloudType>::dh
} }
case (parent::etEnthalpyDifference): 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); scalar hp = liquids_.properties()[idl].h(p, T);
dh = hc - hp; dh = hc - hp;

6
src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporationBoil/LiquidEvaporationBoil.C
View File

@ -171,8 +171,8 @@ void Foam::LiquidEvaporationBoil<CloudType>::calculate
forAll(this->owner().thermo().carrier().Y(), i) forAll(this->owner().thermo().carrier().Y(), i)
{ {
scalar Yc = this->owner().thermo().carrier().Y()[i][cellI]; scalar Yc = this->owner().thermo().carrier().Y()[i][cellI];
Hc += Yc*this->owner().thermo().carrier().Hs(i, pc, Tc); Hc += Yc*this->owner().thermo().carrier().Ha(i, pc, Tc);
Hsc += Yc*this->owner().thermo().carrier().Hs(i, ps, Ts); Hsc += Yc*this->owner().thermo().carrier().Ha(i, ps, Ts);
Cpc += Yc*this->owner().thermo().carrier().Cp(i, ps, Ts); Cpc += Yc*this->owner().thermo().carrier().Cp(i, ps, Ts);
kappac += Yc*this->owner().thermo().carrier().kappa(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): 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); scalar hp = liquids_.properties()[idl].h(p, TDash);
dh = hc - hp; dh = hc - hp;