Foam::CompositionModel<CloudType>::CompositionModel
(
const dictionary& dict,
CloudType& owner,
const word& type
)
:
dict_(dict),
owner_(owner),
coeffDict_(dict.subDict(type + "Coeffs")),
mcCarrierThermo_(owner.mcCarrierThermo()),
liquids_
(
liquidMixture::New
(
owner.mesh().objectRegistry::lookupObject<dictionary>
(
owner.carrierThermo().name()
)
)
),
solids_
(
solidMixture::New
(
owner.mesh().objectRegistry::lookupObject<dictionary>
(
owner.carrierThermo().name()
)
)
),
phaseProps_
(
coeffDict_.lookup("phases"),
mcCarrierThermo_.species(),
liquids_().components(),
solids_().components()
)
{}
Foam::LiquidEvaporation<CloudType>::LiquidEvaporation
(
const dictionary& dict,
CloudType& owner
)
:
PhaseChangeModel<CloudType>(dict, owner, typeName),
liquids_
(
liquidMixture::New
(
owner.mesh().objectRegistry::template
lookupObject<dictionary>
(
owner.carrierThermo().name()
)
)
),
activeLiquids_(this->coeffDict().lookup("activeLiquids")),
liqToCarrierMap_(activeLiquids_.size(), -1),
liqToLiqMap_(activeLiquids_.size(), -1)
{
if (activeLiquids_.size() == 0)
{
WarningIn
(
"Foam::LiquidEvaporation<CloudType>::LiquidEvaporation"
"("
"const dictionary& dict, "
"CloudType& owner"
")"
) << "Evaporation model selected, but no active liquids defined"
<< nl << endl;
}
// Determine mapping between liquid and carrier phase species
forAll(activeLiquids_, i)
{
liqToCarrierMap_[i] =
owner.composition().globalCarrierId(activeLiquids_[i]);
}
