30 template<
class CloudType>
42 return 24.0*(1.0 + (1.0/6.0)*
pow(Re, 2.0/3.0));
49 template<
class CloudType>
61 template<
class CloudType>
73 template<
class CloudType>
80 template<
class CloudType>
93 scalar
y =
min(
max(p.y(), 0), 1);
95 value.
Sp() = mass*0.75*muc*CdRe(Re)*(1 + 2.632*
y)/(p.rho()*
sqr(p.d()));
A list of keyword definitions, which are a keyword followed by any number of values (e...
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
DistortedSphereDragForce(CloudType &owner, const fvMesh &mesh, const dictionary &dict)
Construct from mesh.
Abstract base class for particle forces.
Helper container for force Su and Sp terms.
Drag model based on assumption of distorted spheres according to:
virtual ~DistortedSphereDragForce()
Destructor.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Mesh data needed to do the Finite Volume discretisation.
Templated base class for dsmc cloud.
scalarField Re(const UList< complex > &cf)
scalar Sp() const
Return const access to the implicit coefficient [kg/s].
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.