65 if (
dgdt[celli] > 0.0)
70 else if (
dgdt[celli] < 0.0)
#define forAll(list, i)
Loop across all elements in list.
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
const dimensionSet & dimensions() const
Return dimensions.
static tmp< DimensionedField< Type, GeoMesh > > New(const word &name, const Mesh &mesh, const dimensionSet &, const Field< Type > &)
Return a temporary field constructed from name, mesh,.
const word & name() const
Return name.
A special matrix type and solver, designed for finite volume solutions of scalar equations....
tmp< volScalarField::Internal > Su() const
Return the explicit source.
tmp< volScalarField::Internal > Sp() const
Return the implicit source.
Foam::fvModels & fvModels() const
Return the fvModels that are created on demand.
const fvMesh & mesh
Region mesh.
virtual void alphaSuSp(tmp< volScalarField::Internal > &Su, tmp< volScalarField::Internal > &Sp)
Calculate the alpha equation sources.
volScalarField::Internal dgdt
Compressibility source.
volScalarField & alpha1
Reference to the phase1-fraction.
A class for managing temporary objects.
Calculate the divergence of the given field.
tmp< VolField< Type > > Su(const VolField< Type > &su, const VolField< Type > &vf)
tmp< VolField< Type > > Sp(const volScalarField &sp, const VolField< Type > &vf)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)