Commit 215a3267 authored by andy's avatar andy
Browse files

ENH: Updated construction of he thermo boundary conditions from T field

parent 2914433a
......@@ -35,6 +35,39 @@ License
// * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
template<class BasicThermo, class MixtureType>
Foam::wordList Foam::heThermo<BasicThermo, MixtureType>::heBoundaryBaseTypes()
{
const volScalarField::GeometricBoundaryField& tbf =
this->T_.boundaryField();
wordList hbt(tbf.size(), word::null);
forAll(tbf, patchi)
{
if (isA<fixedJumpFvPatchScalarField>(tbf[patchi]))
{
const fixedJumpFvPatchScalarField& pf =
dynamic_cast<const fixedJumpFvPatchScalarField&>(tbf[patchi]);
hbt[patchi] = pf.interfaceFieldType();
}
else if (isA<fixedJumpAMIFvPatchScalarField>(tbf[patchi]))
{
const fixedJumpAMIFvPatchScalarField& pf =
dynamic_cast<const fixedJumpAMIFvPatchScalarField&>
(
tbf[patchi]
);
hbt[patchi] = pf.interfaceFieldType();
}
}
return hbt;
}
template<class BasicThermo, class MixtureType>
Foam::wordList Foam::heThermo<BasicThermo, MixtureType>::heBoundaryTypes()
{
......@@ -149,7 +182,8 @@ Foam::heThermo<BasicThermo, MixtureType>::heThermo
),
mesh,
dimEnergy/dimMass,
this->heBoundaryTypes()
this->heBoundaryTypes(),
this->heBoundaryBaseTypes()
)
{
init();
......@@ -179,7 +213,8 @@ Foam::heThermo<BasicThermo, MixtureType>::heThermo
),
mesh,
dimEnergy/dimMass,
this->heBoundaryTypes()
this->heBoundaryTypes(),
this->heBoundaryBaseTypes()
)
{
init();
......
......@@ -68,6 +68,10 @@ protected:
// by interrogating the temperature field boundary types
wordList heBoundaryTypes();
//- Return the enthalpy/internal energy field boundary base types
// by interrogating the temperature field boundary types
wordList heBoundaryBaseTypes();
//- Correct the enthalpy/internal energy field boundaries
void heBoundaryCorrection(volScalarField& he);
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment