33 template<
class CloudType>
44 template<
class CloudType>
56 template<
class CloudType>
63 template<
class CloudType>
74 Random& rnd = this->owner().rndGen();
76 const scalar cps = 0.16432;
78 const scalar
k = this->kPtr_->primitiveField()[celli];
80 this->epsilonPtr_->primitiveField()[celli] + rootVSmall;
82 const scalar UrelMag =
mag(
U - Uc - UTurb);
84 const scalar tTurbLoc =
103 const scalar u = 2*rnd.
scalar01() - 1;
105 const scalar a =
sqrt(1 -
sqr(u));
Templated base class for dsmc cloud.
Base class for particle dispersion models based on RAS turbulence.
scalar scalar01()
Advance the state and return a scalar sample from a uniform.
scalar scalarNormal()
Advance the state and return a scalar sample from a normal.
The velocity is perturbed in random direction, with a Gaussian random number distribution with varian...
virtual ~StochasticDispersionRAS()
Destructor.
StochasticDispersionRAS(const dictionary &dict, CloudType &owner)
Construct from components.
virtual vector update(const scalar dt, const label celli, const vector &U, const vector &Uc, vector &UTurb, scalar &tTurb)
Update (disperse particles)
A list of keyword definitions, which are a keyword followed by any number of values (e....
const scalar twoPi(2 *pi)
const dimensionedScalar sigma
Stefan-Boltzmann constant: default SI units: [W/m^2/K^4].
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sin(const dimensionedScalar &ds)
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensioned< scalar > mag(const dimensioned< Type > &)
dimensionedScalar cos(const dimensionedScalar &ds)