33template<
class CloudType>
46 this->coeffs().getWord(
"alphac")
52template<
class CloudType>
63 this->coeffs().getWord(
"alphac")
71template<
class CloudType>
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()))
static const Foam::dimensionedScalar A("", Foam::dimPressure, 611.21)
static const Foam::dimensionedScalar B("", Foam::dimless, 18.678)
Templated base class for dsmc cloud.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
Abstract base class for particle forces.
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.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
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 pow3(const dimensionedScalar &ds)
static constexpr const zero Zero
Global zero (0)
dimensionedScalar cbrt(const dimensionedScalar &ds)