33template<
class CloudType>
42 inputFileName_(this->coeffDict().
lookup(
"inputFile")),
43 duration_(this->coeffDict().getScalar(
"duration")),
44 parcelsPerSecond_(this->coeffDict().getScalar(
"parcelsPerSecond")),
45 randomise_(this->coeffDict().getBool(
"randomise")),
57 injectorCells_(injectors_.size()),
58 injectorTetFaces_(injectors_.size()),
59 injectorTetPts_(injectors_.size())
69 this->
volumeTotal_ += injectors_[i].mDot()/injectors_[i].rho();
75template<
class CloudType>
82 inputFileName_(im.inputFileName_),
83 duration_(im.duration_),
84 parcelsPerSecond_(im.parcelsPerSecond_),
85 randomise_(im.randomise_),
86 injectors_(im.injectors_),
87 injectorCells_(im.injectorCells_),
88 injectorTetFaces_(im.injectorTetFaces_),
89 injectorTetPts_(im.injectorTetPts_)
95template<
class CloudType>
99 bitSet reject(injectors_.size());
106 !this->findCellAtPosition
109 injectorTetFaces_[i],
112 !this->ignoreOutOfBounds_
121 const label nRejected = reject.
count();
131 Info<<
" " << nRejected
132 <<
" positions rejected, out of bounds" <<
endl;
137template<
class CloudType>
140 return this->SOI_ + duration_;
144template<
class CloudType>
151 if ((time0 >= 0.0) && (time0 < duration_))
153 return floor(injectorCells_.size()*(time1 - time0)*parcelsPerSecond_);
160template<
class CloudType>
168 if ((time0 >= 0.0) && (time0 < duration_))
172 volume += injectors_[i].mDot()/injectors_[i].rho()*(time1 - time0);
180template<
class CloudType>
184 const label nParcels,
195 Random& rnd = this->owner().rndGen();
196 injectorI = rnd.
position<label>(0, injectorCells_.size() - 1);
200 injectorI = parcelI*injectorCells_.size()/nParcels;
203 position = injectors_[injectorI].
x();
204 cellOwner = injectorCells_[injectorI];
205 tetFacei = injectorTetFaces_[injectorI];
206 tetPti = injectorTetPts_[injectorI];
210template<
class CloudType>
214 const label nParcels,
219 label injectorI = parcelI*injectorCells_.size()/nParcels;
222 parcel.U() = injectors_[injectorI].U();
225 parcel.d() = injectors_[injectorI].d();
228 parcel.rho() = injectors_[injectorI].rho();
231 parcel.T() = injectors_[injectorI].T();
234 parcel.Cp() = injectors_[injectorI].Cp();
237 parcel.Y() = injectors_[injectorI].Y();
241template<
class CloudType>
248template<
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.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
const objectRegistry & db() const noexcept
Return the local objectRegistry.
Templated injection model class.
Type position(const Type &start, const Type &end)
Return a sample on the interval [start,end].
Particle injection sources read from look-up table. Each row corresponds to an injection site.
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)
virtual void setProperties(const label parcelI, const label nParcels, const scalar time, typename CloudType::parcelType &parcel)
Set the parcel properties.
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.
virtual scalar userTimeToTime(const scalar theta) const
Convert the user-time (e.g. CA deg) to real-time (s).
const Cmpt & x() const
Access to the vector x component.
A bitSet stores bits (elements with only two states) in packed internal format and supports a variety...
unsigned int count(const bool on=true) const
Count number of bits set.
void flip()
Invert all bits in the addressable region.
void set(const bitSet &bitset)
Set specified bits from another bitset.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
const Time & time() const noexcept
Return time registry.
constant condensation/saturation model.
Lookup type of boundary radiation properties.
A class for handling words, derived from Foam::string.
void inplaceSubset(const BoolListType &select, ListType &input, const bool invert=false)
Inplace extract elements of the input list when select is true.
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
#define forAll(list, i)
Loop across all elements in list.