Added rothalpy support in boundary conditions

This commit is contained in:
Hrvoje Jasak 2014-08-14 17:03:32 +01:00
parent c5a2b9d772
commit ffb866bb9a
3 changed files with 15 additions and 3 deletions

View file

@ -103,7 +103,11 @@ void Foam::fixedEnthalpyFvPatchScalarField::updateCoeffs()
const_cast<fvPatchScalarField&>(thermo.T().boundaryField()[patchi]);
Tw.evaluate();
if (dimensionedInternalField().name() == db().mangleFileName("h"))
if
(
dimensionedInternalField().name() == db().mangleFileName("h")
|| dimensionedInternalField().name() == db().mangleFileName("i")
)
{
operator==(thermo.h(Tw, patchi));
}

View file

@ -104,7 +104,11 @@ void Foam::gradientEnthalpyFvPatchScalarField::updateCoeffs()
Tw.evaluate();
if (dimensionedInternalField().name() == db().mangleFileName("h"))
if
(
dimensionedInternalField().name() == db().mangleFileName("h")
|| dimensionedInternalField().name() == db().mangleFileName("i")
)
{
gradient() = thermo.Cp(Tw, patchi)*Tw.snGrad()
+ patch().deltaCoeffs()*

View file

@ -112,7 +112,11 @@ void Foam::mixedEnthalpyFvPatchScalarField::updateCoeffs()
valueFraction() = Tw.valueFraction();
if (dimensionedInternalField().name() == db().mangleFileName("h"))
if
(
dimensionedInternalField().name() == db().mangleFileName("h")
|| dimensionedInternalField().name() == db().mangleFileName("i")
)
{
refValue() = thermo.h(Tw.refValue(), patchi);
refGrad() = thermo.Cp(Tw, patchi)*Tw.refGrad()