34 template<
class CloudType>
43 positionsFile_(this->coeffDict().lookup(
"positionsFile")),
49 owner.db().time().constant(),
55 diameters_(positions_.size()),
56 injectorCoordinates_(positions_.size(),
barycentric::uniform(NaN)),
57 injectorCells_(positions_.size(), -1),
58 injectorTetFaces_(positions_.size(), -1),
59 injectorTetPts_(positions_.size(), -1),
60 massTotal_(this->readMassTotal(
dict, owner)),
61 U0_(this->coeffDict().lookup(
"U0")),
67 this->coeffDict().subDict(
"sizeDistribution"),
74 this->coeffDict().lookupOrDefault(
"ignoreOutOfBounds", false)
82 diameters_[i] = sizeDistribution_->sample();
87 template<
class CloudType>
94 positionsFile_(im.positionsFile_),
95 positions_(im.positions_),
96 diameters_(im.diameters_),
97 injectorCoordinates_(im.injectorCoordinates_),
98 injectorCells_(im.injectorCells_),
99 injectorTetFaces_(im.injectorTetFaces_),
100 injectorTetPts_(im.injectorTetPts_),
101 massTotal_(im.massTotal_),
103 sizeDistribution_(im.sizeDistribution_, false),
104 ignoreOutOfBounds_(im.ignoreOutOfBounds_)
110 template<
class CloudType>
117 template<
class CloudType>
128 !this->findCellAtPosition
131 injectorCoordinates_[pI],
133 injectorTetFaces_[pI],
153 Info<<
" " << nRejected
154 <<
" particles ignored, out of bounds" <<
endl;
159 template<
class CloudType>
168 template<
class CloudType>
176 if (0 >= time0 && 0 < time1)
178 return positions_.size();
187 template<
class CloudType>
195 if (0 >= time0 && 0 < time1)
206 template<
class CloudType>
219 coordinates = injectorCoordinates_[parcelI];
220 celli = injectorCells_[parcelI];
221 tetFacei = injectorTetFaces_[parcelI];
222 tetPti = injectorTetPts_[parcelI];
226 template<
class CloudType>
232 typename CloudType::parcelType::trackingData& td,
240 parcel.d() = diameters_[parcelI];
244 template<
class CloudType>
#define forAll(list, i)
Loop across all elements in list.
Templated base class for dsmc cloud.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
Templated injection model class.
virtual void topoChange()
Set injector locations when mesh is updated.
virtual ~ManualInjection()
Destructor.
virtual void setProperties(const label parcelI, const label nParcels, const scalar time, typename CloudType::parcelType::trackingData &td, typename CloudType::parcelType &parcel)
Set the parcel properties.
ManualInjection(const dictionary &dict, CloudType &owner, const word &modelName)
Construct from dictionary.
virtual void setPositionAndCell(const label parcelI, const label nParcels, const scalar time, barycentric &coordinates, label &celli, label &tetFacei, label &tetPti, label &facei)
Set the injection position and owner cell, tetFace and tetPt.
virtual label nParcelsToInject(const scalar time0, const scalar time1)
Number of parcels to introduce relative to SOI.
virtual bool fullyDescribed() const
Flag to identify whether model fully describes the parcel.
virtual scalar timeEnd() const
Return the end-of-injection time.
virtual scalar massToInject(const scalar time0, const scalar time1)
Parcel mass to introduce relative to SOI.
A list of keyword definitions, which are a keyword followed by any number of values (e....
Base class for statistical distributions.
A class for handling words, derived from string.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Ostream & endl(Ostream &os)
Add newline and flush stream.
void inplaceSubset(const UList< T > &select, const T &value, ListType &)
Inplace extract elements of List when select is a certain value.
const dimensionSet dimLength
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions)
randomGenerator rndGen(653213)