Go to the documentation of this file.
33 template<
class ParcelType>
34 template<
class TrackCloudType>
37 const TrackCloudType&
cloud,
41 ParcelType::trackingData(
cloud),
46 cloud.solution().interpolationSchemes(),
54 cloud.solution().interpolationSchemes(),
62 cloud.solution().interpolationSchemes(),
80 cloud.solution().dict(),
94 cloud.solution().dict(),
108 cloud.solution().dict(),
122 cloud.solution().dict(),
136 cloud.solution().dict(),
150 cloud.solution().dict(),
164 cloud.solution().dict(),
169 g_(
cloud.g().value()),
174 template<
class ParcelType>
182 template<
class ParcelType>
190 template<
class ParcelType>
198 template<
class ParcelType>
206 template<
class ParcelType>
214 template<
class ParcelType>
221 template<
class ParcelType>
229 template<
class ParcelType>
236 template<
class ParcelType>
243 template<
class ParcelType>
250 template<
class ParcelType>
258 template<
class ParcelType>
266 template<
class ParcelType>
267 template<
class TrackCloudType>
271 const TrackCloudType&
cloud
275 volumeAverage_() = 0;
276 radiusAverage_() = 0;
280 frequencyAverage_() = 0;
294 cloud.solution().dict(),
301 for (
const typename TrackCloudType::parcelType&
p :
cloud)
305 const scalar m =
p.nParticle()*
p.mass();
307 volumeAverage_->add(
p.coordinates(), tetIs,
p.nParticle()*
p.volume());
308 rhoAverage_->add(
p.coordinates(), tetIs, m*
p.rho());
309 uAverage_->add(
p.coordinates(), tetIs, m*
p.U());
310 massAverage_->add(
p.coordinates(), tetIs, m);
312 volumeAverage_->average();
313 massAverage_->average();
314 rhoAverage_->average(*massAverage_);
315 uAverage_->average(*massAverage_);
318 for (
const typename TrackCloudType::parcelType&
p :
cloud)
322 const vector u = uAverage_->interpolate(
p.coordinates(), tetIs);
331 uSqrAverage_->average(*massAverage_);
334 radiusAverage_() = volumeAverage_();
336 for (
const typename TrackCloudType::parcelType&
p :
cloud)
344 p.nParticle()*
pow(
p.volume(), 2.0/3.0)
348 radiusAverage_->average(weightAverage);
352 for (
const typename TrackCloudType::parcelType&
p :
cloud)
356 const scalar a = volumeAverage_->interpolate(
p.coordinates(), tetIs);
357 const scalar r = radiusAverage_->interpolate(
p.coordinates(), tetIs);
358 const vector u = uAverage_->interpolate(
p.coordinates(), tetIs);
360 const scalar
f = 0.75*a/
pow3(r)*
sqr(0.5*
p.d() + r)*
mag(
p.U() - u);
362 frequencyAverage_->add(
p.coordinates(), tetIs,
p.nParticle()*
f*
f);
364 weightAverage.
add(
p.coordinates(), tetIs,
p.nParticle()*
f);
366 frequencyAverage_->average(weightAverage);
Base class for lagrangian averaging methods.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
const vector & Uc() const
Return the continuous phase velocity.
void updateAverages(const TrackCloudType &cloud)
Update the MPPIC averages.
static constexpr const zero Zero
Global zero (0)
static word timeName(const scalar t, const int precision=precision_)
trackPart part() const
Return the part of the tracking operation taking place.
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
trackingData(const TrackCloudType &cloud, trackPart part=tpLinearTrack)
Construct from components.
dimensionedScalar pow3(const dimensionedScalar &ds)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
const interpolation< vector > & UInterp() const
virtual void add(const barycentric &coordinates, const tetIndices &tetIs, const Type &value)=0
Member Functions.
Abstract base class for interpolation.
const word & name() const noexcept
Return name.
Storage and named access for the indices of a tet which is part of the decomposition of a cell.
Pointer management similar to std::unique_ptr, with some additional methods and type checking.
A cloud is a registry collection of lagrangian particles.
const interpolation< scalar > & muInterp() const
dimensionedSymmTensor sqr(const dimensionedVector &dv)
scalar rhoc() const
Return the continuous phase density.
scalar muc() const
Return the continuous phase viscosity.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const interpolation< scalar > & rhoInterp() const
virtual void average()
Calculate the average.
const Time & time() const noexcept
Return time registry.
const objectRegistry & db() const noexcept
Return the local objectRegistry.