30 template<
class CloudType>
40 inputFileName_(this->coeffDict().lookup(
"inputFile")),
52 owner.db().time().constant(),
62 duration_ = owner.db().time().userTimeToTime(duration_);
65 injectorCells_.setSize(injectors_.size());
66 injectorTetFaces_.setSize(injectors_.size());
67 injectorTetPts_.setSize(injectors_.size());
72 this->volumeTotal_ = 0.0;
75 this->volumeTotal_ += injectors_[i].mDot()/injectors_[i].rho();
77 this->volumeTotal_ *= duration_;
81 template<
class CloudType>
89 inputFileName_(im.inputFileName_),
90 duration_(im.duration_),
91 parcelsPerSecond_(im.parcelsPerSecond_),
92 randomise_(im.randomise_),
93 injectors_(im.injectors_),
94 injectorCells_(im.injectorCells_),
95 injectorTetFaces_(im.injectorTetFaces_),
96 injectorTetPts_(im.injectorTetPts_)
102 template<
class CloudType>
110 template<
class CloudType>
116 this->findCellAtPosition
119 injectorTetFaces_[i],
127 template<
class CloudType>
131 return this->SOI_ + duration_;
135 template<
class CloudType>
143 if ((time0 >= 0.0) && (time0 < duration_))
145 return floor(injectorCells_.size()*(time1 - time0)*parcelsPerSecond_);
154 template<
class CloudType>
163 if ((time0 >= 0.0) && (time0 < duration_))
167 volume += injectors_[i].mDot()/injectors_[i].rho()*(time1 - time0);
175 template<
class CloudType>
179 const label nParcels,
191 injectorI = rnd.
position<
label>(0, injectorCells_.size() - 1);
195 injectorI = parcelI*injectorCells_.size()/nParcels;
198 position = injectors_[injectorI].
x();
199 cellOwner = injectorCells_[injectorI];
200 tetFacei = injectorTetFaces_[injectorI];
201 tetPtI = injectorTetPts_[injectorI];
205 template<
class CloudType>
209 const label nParcels,
214 label injectorI = parcelI*injectorCells_.size()/nParcels;
217 parcel.U() = injectors_[injectorI].U();
220 parcel.d() = injectors_[injectorI].d();
223 parcel.rho() = injectors_[injectorI].rho();
226 parcel.T() = injectors_[injectorI].T();
229 parcel.Cp() = injectors_[injectorI].Cp();
232 parcel.Y() = injectors_[injectorI].Y();
235 parcel.YGas() = injectors_[injectorI].YGas();
238 parcel.YLiquid() = injectors_[injectorI].YLiquid();
241 parcel.YSolid() = injectors_[injectorI].YSolid();
245 template<
class CloudType>
253 template<
class CloudType>
#define forAll(list, i)
Loop across all elements in list.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
A list of keyword definitions, which are a keyword followed by any number of values (e...
Templated injection model class.
virtual void setProperties(const label parcelI, const label nParcels, const scalar time, typename CloudType::parcelType &parcel)
Set the parcel properties.
stressControl lookup("compactNormalStress") >> compactNormalStress
A class for handling words, derived from string.
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.
virtual ~ReactingMultiphaseLookupTableInjection()
Destructor.
scalar timeEnd() const
Return the end-of-injection time.
virtual bool fullyDescribed() const
Flag to identify whether model fully describes the parcel.
ReactingMultiphaseLookupTableInjection(const dictionary &dict, CloudType &owner, const word &modelName)
Construct from dictionary.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
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.
virtual scalar volumeToInject(const scalar time0, const scalar time1)
Volume of parcels to introduce relative to SOI.
Type position(const Type &start, const Type &end)
Return a sample between start and end.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
Templated base class for dsmc cloud.
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.