30 template<
class CloudType>
43 template<
class CloudType>
50 template<
class CloudType>
57 template<
class CloudType>
68 const scalar
cps = 0.16432;
73 const scalar UrelMag =
mag(U - Uc - UTurb);
75 const scalar tTurbLoc =
min
77 k[celli]/epsilon[celli],
78 cps*
pow(k[celli], 1.5)/epsilon[celli]/(UrelMag + SMALL)
90 scalar sigma =
sqrt(2.0*k[celli]/3.0);
91 vector dir = 2.0*this->owner().rndGen().vector01() - vector::one;
92 dir /=
mag(dir) + SMALL;
98 while ((rsq > 1.0) || (rsq == 0.0))
100 x1 = 2.0*this->owner().rndGen().scalar01() - 1.0;
101 x2 = 2.0*this->owner().rndGen().scalar01() - 1.0;
105 scalar fac =
sqrt(-2.0*
log(rsq)/rsq);
109 UTurb = sigma*fac*dir;
116 UTurb = vector::zero;