36 namespace diameterModels
48 tmp<volScalarField> tInvDsm
62 const sizeGroup& fi = sizeGroups_[i];
64 invDsm += fi.a()*fi/fi.x();
73 tmp<volScalarField> tN
87 N += phase()*sizeGroups_[i]/sizeGroups_[i].x();
95 Foam::diameterModels::velocityGroup::fSum()
const 97 tmp<volScalarField> tsumSizeGroups
111 sumSizeGroups += sizeGroups_[i];
114 return tsumSizeGroups;
118 void Foam::diameterModels::velocityGroup::scale()
120 Info<<
"Scaling sizeGroups for velocityGroup " << phase().name() <<
endl;
124 sizeGroups_[i].max(0);
131 sizeGroups_[i] /= f_;
133 sizeGroups_[i].correctBoundaryConditions();
142 const dictionary& diameterProperties,
143 const phaseModel& phase
146 diameterModel(diameterProperties, phase),
147 popBalName_(diameterProperties.
lookup(
"populationBalance")),
166 diameterProperties.
lookup(
"sizeGroups"),
167 sizeGroup::iNew(phase, *this)
169 d_(IOobject::groupName(
"d", phase.
name()), dsm())
189 tmp<volScalarField> tA
205 a += fi.a()*fi/fi.x();
216 sizeGroups_[
i].correct();
221 phase().
mesh().solverDict(popBalName_).lookupOrDefault<Switch>
233 f_.correctBoundaryConditions();
235 Info<<
phase().
name() <<
" sizeGroups-sum volume fraction, min, max = " 236 << f_.weightedAverage(
phase().
mesh().V()).value()
237 <<
' ' <<
min(f_).value()
238 <<
' ' <<
max(f_).value()
243 Info<< this->
phase().
name() <<
" Sauter mean diameter, min, max = " 244 << d_.weightedAverage(d_.mesh().V()).value()
245 <<
' ' <<
min(d_).value()
246 <<
' ' <<
max(d_).value()
253 const dictionary& phaseProperties
#define forAll(list, i)
Loop across all elements in list.
virtual tmp< volScalarField > d() const
Get the diameter field.
dimensionedSphericalTensor inv(const dimensionedSphericalTensor &dt)
static tmp< GeometricField< scalar, fvPatchField, volMesh > > New(const word &name, const Internal &, const PtrList< fvPatchField< scalar >> &)
Return a temporary field constructed from name,.
Ostream & endl(Ostream &os)
Add newline and flush stream.
sizeGroup(const word &name, const dictionary &dict, const phaseModel &phase, const velocityGroup &velocityGroup, const fvMesh &mesh)
const word & name() const
const dimensionSet dimless
velocityGroup(const dictionary &diameterProperties, const phaseModel &phase)
Construct from dictionary and phase.
const phaseModel & phase() const
Return const-reference to the phase.
Macros for easy insertion into run-time selection tables.
virtual tmp< volScalarField > a() const
Get the surface area per unit volume field.
const dimensionSet dimLength
GeometricField< scalar, fvPatchField, volMesh > volScalarField
const tmp< volScalarField > a() const
Return representative surface area of the sizeGroup.
stressControl lookup("compactNormalStress") >> compactNormalStress
virtual bool read(const dictionary &phaseProperties)
Read phaseProperties dictionary.
virtual ~velocityGroup()
Destructor.
void min(const dimensioned< Type > &)
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
const Mesh & mesh() const
Return mesh.
defineTypeNameAndDebug(combustionModel, 0)
word name(const complex &)
Return a string representation of a complex.
Internal & ref()
Return a reference to the dimensioned internal field.
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
void max(const dimensioned< Type > &)
const dimensionSet dimVolume
const autoPtr< label > & i() const
Return label of the sizeGroup within the population balance.
virtual bool read(const dictionary &diameterProperties)
Read diameterProperties dictionary.
A class for managing temporary objects.
virtual void correct()
Correct the model.