33 template<
class CloudType>
37 if (alphacRe < 1000.0)
39 return 24.0*(1.0 + 0.15*
pow(alphacRe, 0.687));
48 template<
class CloudType>
61 this->coeffs().getWord(
"alphac")
67 template<
class CloudType>
78 this->coeffs().getWord(
"alphac")
86 template<
class CloudType>
90 const typename CloudType::parcelType::trackingData& td,
97 const scalar alphac = alphac_[
p.cell()];
103 *0.75*CdRe(alphac*
Re)*muc*
pow(alphac, -2.65)/(alphac*
sqr(
p.d()))
WenYuDragForce(CloudType &owner, const fvMesh &mesh, const dictionary &dict)
Construct from mesh.
Particle-drag model wherein drag forces (per unit carrier-fluid velocity) are dynamically computed ba...
A list of keyword definitions, which are a keyword followed by a number of values (eg...
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Abstract base class for particle forces.
Helper container for force Su and Sp terms.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
scalarField Re(const UList< complex > &cmplx)
Extract real component.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Mesh data needed to do the Finite Volume discretisation.
Templated base class for dsmc cloud.
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.
static constexpr const zero Zero
Global zero (0)