Go to the documentation of this file.
34 template<
class CloudType>
47 template<
class CloudType>
63 template<
class CloudType>
72 template<
class CloudType>
79 gradkPtr_ =
fvc::grad(*this->kPtr_).ptr();
94 template<
class CloudType>
107 const scalar cps = 0.16432;
109 const scalar
k = this->kPtr_->primitiveField()[celli];
111 this->epsilonPtr_->primitiveField()[celli] + ROOTVSMALL;
112 const vector& gradk = this->gradkPtr_->primitiveField()[celli];
114 const scalar UrelMag =
mag(
U - Uc - UTurb);
116 const scalar tTurbLoc =
125 if (tTurb > tTurbLoc)
130 const vector dir = -gradk/(
mag(gradk) + SMALL);
virtual vector update(const scalar dt, const label celli, const vector &U, const vector &Uc, vector &UTurb, scalar &tTurb)
Update (disperse particles)
virtual ~GradientDispersionRAS()
Destructor.
virtual void cacheFields(const bool store)
Cache carrier fields.
static constexpr const zero Zero
Global zero (0)
Template functions to aid in the implementation of demand driven data.
bool ownGradK_
Take ownership of the grad(k)
tmp< GeometricField< typename outerProduct< vector, Type >::type, faPatchField, areaMesh >> grad(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
Random & rndGen()
Return reference to the random object.
void deleteDemandDrivenData(DataPtr &dataPtr)
const fvMesh & mesh() const
Return reference to the mesh.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Templated base class for dsmc cloud.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
const volVectorField * gradkPtr_
Gradient of k.
Calculate the second temporal derivative.
label nSolutionD() const
Return the number of valid solved-for dimensions in the mesh.
dimensionedScalar sqrt(const dimensionedScalar &ds)
Calculate the gradient of the given field.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
label k
Boltzmann constant.
The velocity is perturbed in the direction of -grad(k), with a Gaussian random number distribution wi...
dimensionedScalar sigma("sigma", dimMass/sqr(dimTime), transportProperties)
GradientDispersionRAS(const dictionary &dict, CloudType &owner)
Construct from components.
Base class for particle dispersion models based on RAS turbulence.