void Foam::SurfaceFilmModel<CloudType>::setParcelProperties
(
parcelType& p,
const label filmFaceI
) const
{
// Set parcel properties
scalar vol = mathematical::pi/6.0*pow3(diameterParcelPatch_[filmFaceI]);
p.d() = diameterParcelPatch_[filmFaceI];
p.U() = UFilmPatch_[filmFaceI];
p.rho() = rhoFilmPatch_[filmFaceI];
p.nParticle() = massParcelPatch_[filmFaceI]/p.rho()/vol;
if (ejectedParcelType_ >= 0)
{
p.typeId() = ejectedParcelType_;
}
}
void Foam::CellAverageParticleVelocity<CloudType>::postMove
(
const parcelType& p,
const label cellI,
const scalar dt,
const point& position0,
bool&
)
{
volVectorField& U = UpPtr_();
volScalarField& pVol = pVolPtr_();
scalar pv = p.mass()/p.rho();
U[cellI] += dt*p.nParticle()*p.U()*pv;
pVol[cellI] += pv;
}