35 template<
class CloudType>
46 template<
class CloudType>
58 template<
class CloudType>
65 template<
class CloudType>
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 =
105 const scalar u = 2*rnd.
sample01<scalar>() - 1;
107 const scalar a =
sqrt(1 -
sqr(u));
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
A list of keyword definitions, which are a keyword followed by a number of values (eg...
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sqrt(const dimensionedScalar &ds)
label k
Boltzmann constant.
virtual ~StochasticDispersionRAS()
Destructor.
dimensionedScalar sigma("sigma", dimMass/sqr(dimTime), transportProperties)
Type sample01()
Return a sample whose components lie in the range [0,1].
dimensionedScalar cos(const dimensionedScalar &ds)
constexpr scalar twoPi(2 *M_PI)
Random & rndGen()
Return reference to the random object.
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
virtual vector update(const scalar dt, const label celli, const vector &U, const vector &Uc, vector &UTurb, scalar &tTurb)
Update (disperse particles)
dimensionedScalar sin(const dimensionedScalar &ds)
A Vector of values with scalar precision, where scalar is float/double depending on the compilation f...
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Type GaussNormal()
Return a sample whose components are normally distributed with zero mean and unity variance N(0...
Base class for particle dispersion models based on RAS turbulence.
Templated base class for dsmc cloud.
StochasticDispersionRAS(const dictionary &dict, CloudType &owner)
Construct from components.
The velocity is perturbed in random direction, with a Gaussian random number distribution with varian...
static constexpr const zero Zero
Global zero (0)