33template<
class ParcelType>
34template<
class TrackCloudType>
37 const TrackCloudType&
cloud,
174template<
class ParcelType>
182template<
class ParcelType>
190template<
class ParcelType>
198template<
class ParcelType>
206template<
class ParcelType>
214template<
class ParcelType>
221template<
class ParcelType>
229template<
class ParcelType>
236template<
class ParcelType>
243template<
class ParcelType>
250template<
class ParcelType>
258template<
class ParcelType>
266template<
class ParcelType>
267template<
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);
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);
const uniformDimensionedVectorField & g
Base class for lagrangian averaging methods.
virtual void add(const barycentric &coordinates, const tetIndices &tetIs, const Type &value)=0
Member Functions.
virtual void average()
Calculate the average.
dimensioned< Type > average() const
Calculate and return arithmetic average.
Defines the attributes of an object for which implicit objectRegistry management is supported,...
const word & name() const noexcept
Return the object name.
const objectRegistry & db() const noexcept
Return the local objectRegistry.
const vector & Uc() const
Return the continuous phase velocity.
scalar muc() const
Return the continuous phase viscosity.
const interpolation< vector > & UInterp() const
scalar rhoc() const
Return the continuous phase density.
const interpolation< scalar > & muInterp() const
const interpolation< scalar > & rhoInterp() const
trackPart part() const
Return the part of the tracking operation taking place.
void updateAverages(const TrackCloudType &cloud)
Update the MPPIC averages.
static word timeName(const scalar t, const int precision=precision_)
Pointer management similar to std::unique_ptr, with some additional methods and type checking.
A cloud is a registry collection of lagrangian particles.
Class used to pass data into container.
Abstract base class for volume field interpolation.
const Time & time() const noexcept
Return time registry.
Selector class for relaxation factors, solver type and solution.
Storage and named access for the indices of a tet which is part of the decomposition of a cell.
const volScalarField & mu
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar pow3(const dimensionedScalar &ds)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
static constexpr const zero Zero
Global zero (0)
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh > > &tdf1, const word &name, const dimensionSet &dimensions)
Global function forwards to reuseTmpDimensionedField::New.
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for INVALID.
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)