36template<
class CloudType>
46 phiName_(this->coeffDict().template getOrDefault<
word>(
"phi",
"phi")),
47 rhoName_(this->coeffDict().template getOrDefault<
word>(
"rho",
"rho")),
48 duration_(this->coeffDict().getScalar(
"duration")),
60 this->coeffDict().getScalar(
"parcelConcentration")
66 this->coeffDict().subDict(
"sizeDistribution"),
74 concentration_->userTimeToTime(time);
85template<
class CloudType>
93 phiName_(im.phiName_),
94 rhoName_(im.rhoName_),
95 duration_(im.duration_),
96 concentration_(im.concentration_.clone()),
97 parcelConcentration_(im.parcelConcentration_),
98 sizeDistribution_(im.sizeDistribution_.clone())
104template<
class CloudType>
111template<
class CloudType>
118template<
class CloudType>
121 return this->SOI_ + duration_;
125template<
class CloudType>
135 scalar flowRateIn = 0.0;
138 flowRateIn =
max(0.0, -
sum(phip));
146 flowRateIn =
max(0.0, -
sum(phip/rhop));
155template<
class CloudType>
162 if ((time0 >= 0.0) && (time0 < duration_))
164 scalar dt = time1 - time0;
166 scalar c = concentration_->
value(0.5*(time0 + time1));
168 scalar nParcels = parcelConcentration_*c*flowRate()*dt;
170 Random& rnd = this->owner().rndGen();
172 label nParcelsToInject = floor(nParcels);
180 nParcels - scalar(nParcelsToInject)
188 return nParcelsToInject;
195template<
class CloudType>
204 if ((time0 >= 0.0) && (time0 < duration_))
206 scalar c = concentration_->
value(0.5*(time0 + time1));
208 volume = c*(time1 - time0)*flowRate();
211 this->volumeTotal_ = volume;
212 this->massTotal_ = volume*this->owner().constProps().rho0();
218template<
class CloudType>
232 this->owner().
mesh(),
242template<
class CloudType>
252 parcel.U() = this->owner().U()[parcel.cell()];
255 parcel.d() = sizeDistribution_->sample();
259template<
class CloudType>
266template<
class CloudType>
const CloudType & owner() const
Return const access to the owner cloud.
Templated base class for dsmc cloud.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
const fvMesh & mesh() const
Return reference to the mesh.
const dimensionSet & dimensions() const
Return dimensions.
Top level data entry class for use in dictionaries. Provides a mechanism to specify a variable as a c...
const Boundary & boundaryField() const
Return const-reference to the boundary field.
const objectRegistry & db() const noexcept
Return the local objectRegistry.
Templated injection model class.
scalar massTotal_
Total mass to inject [kg].
Patch injection, by using patch flow rate to determine concentration and velocity.
virtual scalar volumeToInject(const scalar time0, const scalar time1)
Volume of parcels to introduce relative to SOI.
virtual label parcelsToInject(const scalar time0, const scalar time1)
Number of parcels to introduce relative to SOI.
virtual void setPositionAndCell(const label parcelI, const label nParcels, const scalar time, vector &position, label &cellOwner, label &tetFacei, label &tetPti)
Set the injection position and owner cell, tetFace and tetPt.
virtual bool validInjection(const label parcelI)
Return flag to identify whether or not injection of parcelI is.
virtual ~PatchFlowRateInjection()
Destructor.
virtual void setProperties(const label parcelI, const label nParcels, const scalar time, typename CloudType::parcelType &parcel)
Set the parcel properties.
virtual scalar flowRate() const
Return the total volumetric flow rate across the patch [m3/s].
virtual void updateMesh()
Set injector locations when mesh is updated.
virtual bool fullyDescribed() const
Flag to identify whether model fully describes the parcel.
scalar timeEnd() const
Return the end-of-injection time.
Type globalPosition(const Type &start, const Type &end)
Return a sample on the interval [start,end].
virtual scalar userTimeToTime(const scalar theta) const
Convert the user-time (e.g. CA deg) to real-time (s).
Class to control time during OpenFOAM simulations that is also the top-level objectRegistry.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
const Type & value() const
Return const reference to value.
A library of runtime-selectable doubly-truncated probability distribution models. Returns random samp...
void updateMesh()
Update for new mesh topology.
const Time & time() const noexcept
Return time registry.
const Type & lookupObject(const word &name, const bool recursive=false) const
label setPositionAndCell(const fvMesh &mesh, const scalar fraction01, Random &rnd, vector &position, label &cellOwner, label &tetFacei, label &tetPti)
Set the injection position and owner cell, tetFace and tetPt.
Mesh consisting of general polyhedral cells.
A class for handling words, derived from Foam::string.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensioned< Type > sum(const DimensionedField< Type, GeoMesh > &df)
const dimensionSet dimArea(sqr(dimLength))
const dimensionSet dimVelocity
void reduce(const List< UPstream::commsStruct > &comms, T &value, const BinaryOp &bop, const int tag, const label comm)
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh > > &tdf1, const word &name, const dimensionSet &dimensions)
Global function forwards to reuseTmpDimensionedField::New.