31 template<
class CloudType>
43 return 24.0*(1.0 + 0.15*
pow(Re, 0.687));
50 template<
class CloudType>
61 this->
mesh().template lookupObject<volScalarField>
63 this->coeffs().lookup(
"alphac")
69 template<
class CloudType>
78 this->
mesh().template lookupObject<volScalarField>
80 this->coeffs().lookup(
"alphac")
88 template<
class CloudType>
95 template<
class CloudType>
105 scalar alphac(alphac_[p.cell()]);
113 *(150.0*(1.0 - alphac)/alphac + 1.75*
Re)*muc/(alphac*
sqr(p.d()))
122 *0.75*CdRe(alphac*Re)*muc*
pow(alphac, -2.65)/(alphac*
sqr(p.d()))
ErgunWenYuDragForce(CloudType &owner, const fvMesh &mesh, const dictionary &dict)
Construct from mesh.
A list of keyword definitions, which are a keyword followed by any number of values (e...
Ergun-Wen-Yu drag model for solid spheres.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Abstract base class for particle forces.
virtual forceSuSp calcCoupled(const typename CloudType::parcelType &p, const scalar dt, const scalar mass, const scalar Re, const scalar muc) const
Calculate the coupled force.
Helper container for force Su and Sp terms.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Mesh data needed to do the Finite Volume discretisation.
virtual ~ErgunWenYuDragForce()
Destructor.
Templated base class for dsmc cloud.
scalarField Re(const UList< complex > &cf)