33template<
class CloudType>
45 return 24.0*(1.0 + (1.0/6.0)*
pow(Re, 2.0/3.0));
51template<
class CloudType>
63template<
class CloudType>
75template<
class CloudType>
87 const scalar
y =
min(
max(
p.y(), scalar(0)), scalar(1));
94 mass*0.75*muc*CdRe(
Re)*(1.0 + 2.632*
y)/(
p.rho()*
sqr(
p.d()))
Templated base class for dsmc cloud.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
Particle-drag model wherein drag forces (per unit carrier-fluid velocity) are dynamically computed by...
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.
Abstract base class for particle forces.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
virtual tmp< volScalarField > CdRe() const
Drag coefficient.
Class used to pass data into container.
Helper container for force Su and Sp terms.
Mesh data needed to do the Finite Volume discretisation.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
scalarField Re(const UList< complex > &cf)
Extract real component.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
static constexpr const zero Zero
Global zero (0)