33template<
class CloudType>
38 if (alphacRe < 1000.0)
40 return 24.0*(1.0 + 0.15*
pow(alphacRe, 0.687));
49template<
class CloudType>
62 this->coeffs().getWord(
"alphac")
68template<
class CloudType>
79 this->coeffs().getWord(
"alphac")
87template<
class CloudType>
98 const scalar alphac = alphac_[
p.cell()];
106 *(150.0*(1.0 - alphac)/alphac + 1.75*
Re)*muc/(alphac*
sqr(
p.d()))
115 *0.75*CdRe(alphac*
Re)*muc*
pow(alphac, -2.65)/(alphac*
sqr(
p.d()))
Templated base class for dsmc cloud.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
Particle-drag model wherein drag forces (per unit carrier-fluid velocity) are dynamically computed ba...
virtual forceSuSp calcCoupled(const typename CloudType::parcelType &p, const typename CloudType::parcelType::trackingData &td, const scalar dt, const scalar mass, const scalar Re, const scalar muc) const
Calculate the coupled force.
Abstract base class for particle forces.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
virtual tmp< volScalarField > CdRe() const
Drag coefficient.
Class used to pass data into container.
Helper container for force Su and Sp terms.
Mesh data needed to do the Finite Volume discretisation.
scalarField Re(const UList< complex > &cf)
Extract real component.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
static constexpr const zero Zero
Global zero (0)