chpjz0391 |
May 25, 2016 05:21 |
What does rhoEff mean in KinematicCloud.H?
Hi all,
What does rhoEff mean in KinematicCloud.H?
Code:
inline const Foam::tmp<Foam::volScalarField>
Foam::KinematicCloud<ParcelType>::rhoEff() const
{
tmp<volScalarField> trhoEff
(
new volScalarField
(
IOobject
(
this->name() + "RhoEff",
this->db().time().timeName(),
this->db(),
IOobject::NO_READ,
IOobject::NO_WRITE,
false
),
mesh_,
dimensionedScalar("zero", dimDensity, 0.0)
)
);
scalarField& rhoEff = trhoEff().internalField();
forAllConstIter(typename KinematicCloud<ParcelType>, *this, iter)
{
const ParcelType& p = iter();
const label cellI = p.cell();
rhoEff[cellI] += p.nParticle()*p.mass();
}
rhoEff /= mesh().V();
return trhoEff;
}
Does this mean the particles concentration per cell?
If it does, how do I output this parameter?
|