40 template<
class BasePhaseModel>
43 const phaseSystem& fluid,
44 const word& phaseName,
45 const bool referencePhase,
49 BasePhaseModel(fluid, phaseName, referencePhase, index),
54 fluid.
mesh().solverDict(
"Yi")
57 PtrList<volScalarField>& Y = this->thermo_->composition().Y();
61 if (this->thermo_->composition().solve(i))
63 const label j = YActive_.size();
64 YActive_.resize(j + 1);
65 YActive_.set(j, &Y[i]);
73 template<
class BasePhaseModel>
80 template<
class BasePhaseModel>
83 this->thermo_->composition().normalise();
88 template<
class BasePhaseModel>
95 template<
class BasePhaseModel>
102 const tmp<surfaceScalarField> talphaRhoPhi(this->
alphaRhoPhi());
119 template<
class BasePhaseModel>
123 return this->thermo_->composition().Y();
127 template<
class BasePhaseModel>
131 return this->thermo_->composition().Y(name);
135 template<
class BasePhaseModel>
139 return this->thermo_->composition().Y();
143 template<
class BasePhaseModel>
151 template<
class BasePhaseModel>
virtual PtrList< volScalarField > & YRef()
Access the species mass fractions.
virtual void correctSpecies()
Correct the species fractions.
#define forAll(list, i)
Loop across all elements in list.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
fluidReactionThermo & thermo
virtual tmp< fvScalarMatrix > YiEqn(volScalarField &Yi)
Return the species fraction equation.
Calculate the matrix for the laplacian of the field.
const dimensionSet dimless
virtual tmp< fvScalarMatrix > divj(volScalarField &Yi) const
Return the source term for the given specie mass-fraction.
MultiComponentPhaseModel(const phaseSystem &fluid, const word &phaseName, const bool referencePhase, const label index)
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Calculate the first temporal derivative.
virtual tmp< surfaceScalarField > alphaRhoPhi() const
Return the mass flux of the phase.
virtual const PtrList< volScalarField > & Y() const
Return the species mass fractions.
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
Calculate the matrix for the first temporal derivative.
A templated 1D list of pointers to objects of type <T>, where the size of the array is known and used...
virtual const UPtrList< volScalarField > & YActive() const
Return the active species mass fractions.
Calculate the divergence of the given field.
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
virtual ~MultiComponentPhaseModel()
Destructor.
Calculate the matrix for the divergence of the given field and flux.
#define R(A, B, C, D, E, F, K, M)
A templated 1D list of pointers to objects of type <T>, where the size of the array is known and used...
A class for managing temporary objects.
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
virtual UPtrList< volScalarField > & YActiveRef()
Access the active species mass fractions.
virtual bool pure() const
Return whether the phase is pure (i.e., not multi-component)
Calculate the matrix for implicit and explicit sources.