35template<
class CloudType>
46template<
class CloudType>
58template<
class CloudType>
65template<
class CloudType>
76 Random& rnd = this->owner().rndGen();
78 const scalar cps = 0.16432;
80 const scalar
k = this->kPtr_->primitiveField()[celli];
82 this->epsilonPtr_->primitiveField()[celli] + ROOTVSMALL;
84 const scalar UrelMag =
mag(
U - Uc - UTurb);
86 const scalar tTurbLoc =
99 const scalar sigma =
sqrt(2*
k/3.0);
105 const scalar u = 2*rnd.
sample01<scalar>() - 1;
107 const scalar a =
sqrt(1 -
sqr(u));
Templated base class for dsmc cloud.
Base class for particle dispersion models based on RAS turbulence.
Type sample01()
Return a sample whose components lie in the range [0,1].
The velocity is perturbed in random direction, with a Gaussian random number distribution with varian...
virtual ~StochasticDispersionRAS()
Destructor.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
virtual bool update()
Update the mesh for both mesh motion and topology change.
constexpr scalar twoPi(2 *M_PI)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
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)
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)
dimensionedScalar cos(const dimensionedScalar &ds)