32 template<
class CloudType>
36 { contaminationType::PURE,
"pure" },
37 { contaminationType::SLIGHT,
"slight" },
38 { contaminationType::FULL,
"full" },
44 template<
class CloudType>
47 const scalar
f = 1 + 0.15*
pow(
Re, 0.687);
49 switch (contaminationType_)
51 case contaminationType::PURE:
56 case contaminationType::SLIGHT:
61 case contaminationType::FULL:
76 template<
class CloudType>
85 sigma_(this->coeffs().getScalar(
"sigma")),
88 contaminationTypeNames.
get(
"contamination", this->coeffs())
93 template<
class CloudType>
101 contaminationType_(df.contaminationType_)
107 template<
class CloudType>
111 const typename CloudType::parcelType::trackingData&
td,
118 const scalar Eo =
p.Eo(
td, sigma_);
119 const scalar CdRe =
max(this->CdRe(
Re),
Re*8/3*Eo/(Eo + 4));
List< ReturnType > get(const UPtrList< T > &list, const AccessOp &aop)
List of values generated by applying the access operation to each list item.
A list of keyword definitions, which are a keyword followed by a number of values (eg...
wallPoints::trackData td(isBlockedFace, regionToBlockSize)
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Abstract base class for particle forces.
TomiyamaDragForce(CloudType &owner, const fvMesh &mesh, const dictionary &dict)
Construct from mesh.
Helper container for force Su and Sp terms.
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.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
scalarField Re(const UList< complex > &cmplx)
Extract real component.
Particle-drag model wherein drag forces (per unit carrier-fluid velocity) are dynamically computed us...
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Enum is a wrapper around a list of names/values that represent particular enumeration (or int) values...
Mesh data needed to do the Finite Volume discretisation.
Templated base class for dsmc cloud.
static constexpr const zero Zero
Global zero (0)