33 template<
class ParcelType>
34 template<
class TrackCloudType>
37 const TrackCloudType&
cloud,
41 ParcelType::trackingData(
cloud),
52 cloud.solution().dict(),
66 cloud.solution().dict(),
80 cloud.solution().dict(),
94 cloud.solution().dict(),
108 cloud.solution().dict(),
122 cloud.solution().dict(),
136 cloud.solution().dict(),
144 template<
class ParcelType>
145 template<
class TrackCloudType>
148 const TrackCloudType&
cloud
152 volumeAverage_() = 0;
153 radiusAverage_() = 0;
157 frequencyAverage_() = 0;
171 cloud.solution().dict(),
178 for (
const typename TrackCloudType::parcelType&
p :
cloud)
182 const scalar m =
p.nParticle()*
p.mass();
184 volumeAverage_->add(
p.coordinates(), tetIs,
p.nParticle()*
p.volume());
185 rhoAverage_->add(
p.coordinates(), tetIs, m*
p.rho());
186 uAverage_->add(
p.coordinates(), tetIs, m*
p.U());
187 massAverage_->add(
p.coordinates(), tetIs, m);
189 volumeAverage_->average();
190 massAverage_->average();
191 rhoAverage_->average(*massAverage_);
192 uAverage_->average(*massAverage_);
195 for (
const typename TrackCloudType::parcelType&
p :
cloud)
199 const vector u = uAverage_->interpolate(
p.coordinates(), tetIs);
208 uSqrAverage_->average(*massAverage_);
211 radiusAverage_() = volumeAverage_();
213 for (
const typename TrackCloudType::parcelType&
p :
cloud)
221 p.nParticle()*
pow(
p.volume(), 2.0/3.0)
225 radiusAverage_->average(weightAverage);
229 for (
const typename TrackCloudType::parcelType&
p :
cloud)
233 const scalar a = volumeAverage_->interpolate(
p.coordinates(), tetIs);
234 const scalar r = radiusAverage_->interpolate(
p.coordinates(), tetIs);
235 const vector u = uAverage_->interpolate(
p.coordinates(), tetIs);
237 const scalar
f = 0.75*a/
pow3(r)*
sqr(0.5*
p.d() + r)*
mag(
p.U() - u);
239 frequencyAverage_->add(
p.coordinates(), tetIs,
p.nParticle()*
f*
f);
241 weightAverage.
add(
p.coordinates(), tetIs,
p.nParticle()*
f);
243 frequencyAverage_->average(weightAverage);
247 template<
class ParcelType>
255 template<
class ParcelType>