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()]);
107 scalar cbrtAlphap(
pow(1.0 - alphac, 1.0/3.0));
114 *
sqr(1.0 -
sqr(cbrtAlphap))
120 /
sqr(1.0 -
sqr(cbrtAlphap));
126 *(A*(1.0 - alphac)/alphac + B*
Re)*muc/(alphac*
sqr(p.d()))
A list of keyword definitions, which are a keyword followed by any number of values (e...
dimensionedSymmTensor sqr(const dimensionedVector &dv)
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.
Abstract base class for particle forces.
Helper container for force Su and Sp terms.
PlessisMasliyahDragForce(CloudType &owner, const fvMesh &mesh, const dictionary &dict)
Construct from mesh.
virtual ~PlessisMasliyahDragForce()
Destructor.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
PlessisMasliyahDragForce drag model for solid spheres.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar pow3(const dimensionedScalar &ds)
Mesh data needed to do the Finite Volume discretisation.
Templated base class for dsmc cloud.
scalarField Re(const UList< complex > &cf)