35template<
class CloudType>
43 Tref_(this->coeffDict().getScalar(
"Tref"))
49template<
class CloudType>
56template<
class CloudType>
63template<
class CloudType>
72 label typeIdP = pP.typeId();
73 label typeIdQ = pQ.typeId();
78 cloud.constProps(typeIdP).d()
79 +
cloud.constProps(typeIdQ).d()
85 cloud.constProps(typeIdP).omega()
86 +
cloud.constProps(typeIdQ).omega()
89 scalar cR =
mag(pP.U() - pQ.U());
96 scalar mP =
cloud.constProps(typeIdP).mass();
98 scalar mQ =
cloud.constProps(typeIdQ).mass();
100 scalar mR = mP*mQ/(mP + mQ);
105 *
pow(2.0*physicoChemical::k.value()*Tref_/(mR*cR*cR), omegaPQ - 0.5)
112template<
class CloudType>
121 label typeIdP = pP.typeId();
122 label typeIdQ = pQ.typeId();
128 scalar mP =
cloud.constProps(typeIdP).mass();
130 scalar mQ =
cloud.constProps(typeIdQ).mass();
132 vector Ucm = (mP*UP + mQ*UQ)/(mP + mQ);
134 scalar cR =
mag(UP - UQ);
138 scalar sinTheta =
sqrt(1.0 - cosTheta*cosTheta);
142 vector postCollisionRelU =
151 UP = Ucm + postCollisionRelU*mQ/(mP + mQ);
153 UQ = Ucm - postCollisionRelU*mP/(mP + mQ);
Templated DSMC particle collision class.
Templated base class for dsmc cloud.
ParcelType parcelType
Type of parcel the cloud was instantiated for.
Type sample01()
Return a sample whose components lie in the range [0,1].
Variable Hard Sphere BinaryCollision Model.
virtual scalar sigmaTcR(const typename CloudType::parcelType &pP, const typename CloudType::parcelType &pQ) const
Return the collision cross section * relative velocity product.
virtual ~VariableHardSphere()
Destructor.
virtual bool active() const
Flag to indicate whether model activates collision model.
A cloud is a registry collection of lagrangian particles.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
constexpr scalar pi(M_PI)
constexpr scalar twoPi(2 *M_PI)
dimensionedScalar exp(const dimensionedScalar &ds)
dimensionedScalar lgamma(const dimensionedScalar &ds)
dimensionedScalar sin(const dimensionedScalar &ds)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
dimensionedScalar cos(const dimensionedScalar &ds)