33 template<
class CloudType>
46 this->coeffs().getWord(
"alphac")
52 template<
class CloudType>
63 this->coeffs().getWord(
"alphac")
71 template<
class CloudType>
75 const typename CloudType::parcelType::trackingData& td,
82 const scalar alphac = alphac_[
p.cell()];
84 const scalar cbrtAlphap =
cbrt(1.0 - alphac);
91 *
sqr(1.0 -
sqr(cbrtAlphap))
96 const scalar
B =
sqr(alphac)/
sqr(1.0 -
sqr(cbrtAlphap));
102 *(
A*(1.0 - alphac)/alphac +
B*
Re)*muc/(alphac*
sqr(
p.d()))
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.
dimensionedScalar cbrt(const dimensionedScalar &ds)
PlessisMasliyahDragForce(CloudType &owner, const fvMesh &mesh, const dictionary &dict)
Construct from mesh.
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.
scalarField Re(const UList< complex > &cmplx)
Extract real component.
dimensionedScalar pow3(const dimensionedScalar &ds)
Mesh data needed to do the Finite Volume discretisation.
static const Foam::dimensionedScalar A("", Foam::dimPressure, 611.21)
Templated base class for dsmc cloud.
static const Foam::dimensionedScalar B("", Foam::dimless, 18.678)
static constexpr const zero Zero
Global zero (0)