zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

fvModels: phaseChange: Corrections

Browse Source 
The form of the sensible energy/enthalpy transfer and the sign of the
latent heat contribution have both been corrected
This commit is contained in:
Will Bainbridge
2024-06-21 16:40:13 +01:00
parent adfc80c412
commit 03bfaaa4fa
2 changed files with 34 additions and 28 deletions
Download Patch File Download Diff File Expand all files Collapse all files

30
src/fvModels/derived/phaseChange/multicomponentPhaseChange.C
View File

@ -126,7 +126,19 @@ void Foam::fv::multicomponentPhaseChange::addSup
const label speciei =
specieThermos()[i].species()[species()[mDoti]];
const volScalarField::Internal mDotIn(s*mDot(mDoti));
const volScalarField::Internal mDot(this->mDot(mDoti));
// Direct transfer of energy due to mass transfer
tmp<volScalarField::Internal> hs =
vfToVif(specieThermos()[i].hsi(speciei, p, Tchange));
if (energySemiImplicit_)
{
eqn += -fvm::SuSp(-s*mDot, heOrYi) + s*mDot*(hs - heOrYi);
}
else
{
eqn += s*mDot*hs;
}
// Absolute enthalpies at the interface
Pair<tmp<volScalarField::Internal>> has;
@ -135,20 +147,10 @@ void Foam::fv::multicomponentPhaseChange::addSup
has[j] = vfToVif(specieThermos()[j].hai(speciei, p, Tchange));
}
// Direct transfer of energy due to mass transfer
if (energySemiImplicit_)
{
eqn += -fvm::SuSp(-mDotIn, heOrYi) + mDotIn*(has[i]() - heOrYi);
}
else
{
eqn += mDotIn*has[i]();
}
// Latent heat of phase change
eqn +=
(i == 0 ? Lfraction() - 1 : Lfraction())
*mDotIn*(has.second() - has.first());
eqn -=
(i == 0 ? 1 - Lfraction() : Lfraction())
*mDot*(has.second() - has.first());
}
}
// Mass fraction equation

32
src/fvModels/derived/phaseChange/singleComponentPhaseChange.C
View File

@ -120,7 +120,21 @@ void Foam::fv::singleComponentPhaseChange::addSup
const volScalarField& p = this->p();
const volScalarField Tchange(vifToVf(this->Tchange()));
const volScalarField::Internal mDotIn(s*mDot());
const volScalarField::Internal mDot(this->mDot());
// Direct transfer of energy due to mass transfer
tmp<volScalarField::Internal> hs =
specieThermos().valid()[i]
? vfToVif(specieThermos()[i].hsi(specieis_[i], p, Tchange))
: vfToVif(thermos()[i].hs(p, Tchange));
if (energySemiImplicit_)
{
eqn += -fvm::SuSp(-s*mDot, heOrYi) + s*mDot*(hs - heOrYi);
}
else
{
eqn += s*mDot*hs;
}
// Absolute enthalpies at the interface
Pair<tmp<volScalarField::Internal>> has;
@ -132,20 +146,10 @@ void Foam::fv::singleComponentPhaseChange::addSup
: vfToVif(thermos()[j].ha(p, Tchange));
}
// Direct transfer of energy due to mass transfer
if (energySemiImplicit_)
{
eqn += -fvm::SuSp(-mDotIn, heOrYi) + mDotIn*(has[i]() - heOrYi);
}
else
{
eqn += mDotIn*has[i]();
}
// Latent heat of phase change
eqn +=
(i == 0 ? Lfraction() - 1 : Lfraction())
*mDotIn*(has.second() - has.first());
eqn -=
(i == 0 ? 1 - Lfraction() : Lfraction())
*mDot*(has.second() - has.first());
}
// Mass fraction equation
else if