33 template<
class CloudType>
43 inputFileName_(this->coeffDict().
lookup(
"inputFile")),
44 duration_(this->coeffDict().getScalar(
"duration")),
45 parcelsPerSecond_(this->coeffDict().getScalar(
"parcelsPerSecond")),
46 randomise_(this->coeffDict().getBool(
"randomise")),
52 owner.db().time().constant(),
58 injectorCells_(injectors_.size()),
59 injectorTetFaces_(injectors_.size()),
60 injectorTetPts_(injectors_.size())
70 this->
volumeTotal_ += injectors_[i].mDot()/injectors_[i].rho();
76 template<
class CloudType>
84 inputFileName_(im.inputFileName_),
85 duration_(im.duration_),
86 parcelsPerSecond_(im.parcelsPerSecond_),
87 randomise_(im.randomise_),
88 injectors_(im.injectors_),
89 injectorCells_(im.injectorCells_),
90 injectorTetFaces_(im.injectorTetFaces_),
91 injectorTetPts_(im.injectorTetPts_)
97 template<
class CloudType>
101 bitSet reject(injectors_.size());
108 !this->findCellAtPosition
111 injectorTetFaces_[i],
114 !this->ignoreOutOfBounds_
123 label nRejected = reject.count();
133 Info<<
" " << nRejected
134 <<
" positions rejected, out of bounds" <<
endl;
139 template<
class CloudType>
143 return this->SOI_ + duration_;
147 template<
class CloudType>
155 if ((time0 >= 0.0) && (time0 < duration_))
157 return floor(injectorCells_.size()*(time1 - time0)*parcelsPerSecond_);
164 template<
class CloudType>
173 if ((time0 >= 0.0) && (time0 < duration_))
177 volume += injectors_[i].mDot()/injectors_[i].rho()*(time1 - time0);
185 template<
class CloudType>
189 const label nParcels,
200 Random& rnd = this->owner().rndGen();
201 injectorI = rnd.
position<label>(0, injectorCells_.size() - 1);
205 injectorI = parcelI*injectorCells_.size()/nParcels;
208 position = injectors_[injectorI].x();
209 cellOwner = injectorCells_[injectorI];
210 tetFacei = injectorTetFaces_[injectorI];
211 tetPti = injectorTetPts_[injectorI];
215 template<
class CloudType>
219 const label nParcels,
224 label injectorI = parcelI*injectorCells_.size()/nParcels;
227 parcel.U() = injectors_[injectorI].U();
230 parcel.d() = injectors_[injectorI].d();
233 parcel.rho() = injectors_[injectorI].rho();
236 parcel.T() = injectors_[injectorI].T();
239 parcel.Cp() = injectors_[injectorI].Cp();
242 parcel.Y() = injectors_[injectorI].Y();
245 parcel.YGas() = injectors_[injectorI].YGas();
248 parcel.YLiquid() = injectors_[injectorI].YLiquid();
251 parcel.YSolid() = injectors_[injectorI].YSolid();
255 template<
class CloudType>
263 template<
class CloudType>
scalar timeEnd() const
Return the end-of-injection time.
A list of keyword definitions, which are a keyword followed by a number of values (eg...
Templated injection model class.
Ostream & endl(Ostream &os)
Add newline and flush stream.
virtual void setProperties(const label parcelI, const label nParcels, const scalar time, typename CloudType::parcelType &parcel)
Set the parcel properties.
Lookup type of boundary radiation properties.
void inplaceSubset(const BoolListType &select, ListType &input, const bool invert=false)
Inplace extract elements of the input list when select is true.
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.
const CloudType & owner() const
Return const access to the owner cloud.
#define forAll(list, i)
Loop across all elements in list.
virtual scalar userTimeToTime(const scalar theta) const
Convert the user-time (e.g. CA deg) to real-time (s).
A class for handling words, derived from Foam::string.
const Time & time() const noexcept
Return time registry.
virtual label parcelsToInject(const scalar time0, const scalar time1)
Number of parcels to introduce relative to SOI.
virtual bool validInjection(const label parcelI)
Return flag to identify whether or not injection of parcelI is permitted.
const objectRegistry & db() const noexcept
Return the local objectRegistry.
ReactingMultiphaseLookupTableInjection(const dictionary &dict, CloudType &owner, const word &modelName)
Construct from dictionary.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
const wordList volume
Standard volume field types (scalar, vector, tensor, etc)
A Vector of values with scalar precision, where scalar is float/double depending on the compilation f...
Particle injection sources read from look-up table. Each row corresponds to an injection site...
virtual void updateMesh()
Set injector locations when mesh is updated.
A bitSet stores bits (elements with only two states) in packed internal format and supports a variety...
messageStream Info
Information stream (stdout output on master, null elsewhere)
scalar volumeTotal_
Total volume of particles introduced by this injector [m^3] Note: scaled to ensure massTotal is achie...
virtual bool fullyDescribed() const
Flag to identify whether model fully describes the parcel.
virtual scalar volumeToInject(const scalar time0, const scalar time1)
Volume of parcels to introduce relative to SOI.
Defines the attributes of an object for which implicit objectRegistry management is supported...
Templated base class for dsmc cloud.
Type position(const Type &start, const Type &end)
Return a sample on the interval [start,end].