43 namespace diameterModels
61 return phase()*kappai_;
69 const dictionary& diameterProperties,
70 const phaseModel& phase
73 diameterModel(diameterProperties, phase),
78 IOobject::groupName(
"kappai", phase.
name()),
86 dMax_(
"dMax",
dimLength, diameterProperties),
87 dMin_(
"dMin",
dimLength, diameterProperties),
88 residualAlpha_(
"residualAlpha",
dimless, diameterProperties),
90 sources_(diameterProperties.
lookup(
"sources"), IATEsource::iNew(*this))
106 return max(6/
max(kappai_, 6/dMax_), dMin_);
144 R += sources_[j].R(alphaAv, kappai_);
181 IATEsource::iNew(*
this)
182 ).transfer(sources_);
fvMatrix< scalar > fvScalarMatrix
#define forAll(list, i)
Loop across all elements in list.
tmp< fvMatrix< Type > > SuSp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
tmp< fvMatrix< Type > > Sp(const volScalarField::Internal &, const GeometricField< Type, fvPatchField, volMesh > &)
Fundamental dimensioned constants.
Macros for easy insertion into run-time selection tables.
Area-weighted average a surfaceField creating a volField.
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
virtual tmp< volScalarField > calcA() const
Get the surface area per unit volume field.
Calculate the first temporal derivative.
virtual bool read(const dictionary &phaseProperties)
Read phaseProperties dictionary.
stressControl lookup("compactNormalStress") >> compactNormalStress
const phaseModel & phase() const
Return the phase.
virtual bool read(const dictionary &phaseProperties)
Read phaseProperties dictionary.
virtual ~IATE()
Destructor.
tmp< fvMatrix< Type > > ddt(const GeometricField< Type, fvPatchField, volMesh > &vf)
Calculate the matrix for the first temporal derivative.
surfaceScalarField alphaPhi(phi.name()+alpha1.name(), fvc::flux(phi, alpha1, alphaScheme))
IATE(const dictionary &diameterProperties, const phaseModel &phase)
Construct from components.
Calculate the divergence of the given field.
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
defineTypeNameAndDebug(combustionModel, 0)
word name(const complex &)
Return a string representation of a complex.
Calculate the matrix for the divergence of the given field and flux.
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
#define R(A, B, C, D, E, F, K, M)
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
tmp< GeometricField< Type, fvPatchField, volMesh > > average(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Area-weighted average a surfaceField creating a volField.
virtual void correct()
Correct the diameter field.
virtual tmp< volScalarField > calcD() const
Get the diameter field.
A class for managing temporary objects.
static options & New(const fvMesh &mesh)
Construct fvOptions and register to datbase if not present.
const dictionary & diameterProperties() const
Return the phase diameter properties dictionary.
Calculate the matrix for implicit and explicit sources.
ITstream & lookup(const word &, bool recursive=false, bool patternMatch=true) const
Find and return an entry data stream.