33 template<
class CloudType>
52 template<
class CloudType>
63 template<
class CloudType>
90 const scalar We = rhoc*
sqr(Urmag)*d/
sigma;
96 const scalar Wec = 12.0*(1.0 + 1.077*
pow(Oh, 1.6));
101 scalar taubBar = 5.5;
108 taubBar = 0.766*
pow(We - 12.0, 0.25);
113 taubBar = 14.1*
pow(We - 12.0, -0.25);
118 taubBar = 2.45*
pow(We - 12.0, 0.25);
123 taubBar = 6.0*
pow(We - 12.0, -0.25);
132 const scalar rho12 =
sqrt(rhoc/
rho);
135 const scalar Vd = Urmag*rho12*(B1_*taubBar + B2_*
sqr(taubBar));
138 const scalar Vd1 =
max(
sqr(1.0 - Vd/Urmag), SMALL);
139 const scalar dStable = Wec*
sigma/(Vd1*rhoc*
sqr(Urmag));
149 const scalar semiMass = nParticle*
pow3(d);
152 const scalar taub = taubBar*d/(Urmag*rho12);
155 const scalar frac = dt/taub;
156 d = (d + frac*dStable)/(1.0 + frac);
159 nParticle = semiMass/
pow3(d);
DSMCCloud< dsmcParcel > CloudType
A list of keyword definitions, which are a keyword followed by a number of values (eg...
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar sqrt(const dimensionedScalar &ds)
virtual bool update(const scalar dt, const vector &g, scalar &d, scalar &tc, scalar &ms, scalar &nParticle, scalar &KHindex, scalar &y, scalar &yDot, const scalar d0, const scalar rho, const scalar mu, const scalar sigma, const vector &U, const scalar rhoc, const scalar muc, const vector &Urel, const scalar Urmag, const scalar tMom, scalar &dChild, scalar &massChild)
Update the parcel properties.
dimensionedScalar sigma("sigma", dimMass/sqr(dimTime), transportProperties)
bool readEntry(const word &keyword, T &val, enum keyType::option matchOpt=keyType::REGEX, IOobjectOption::readOption readOpt=IOobjectOption::MUST_READ) const
Find entry and assign to T val. FatalIOError if it is found and the number of tokens is incorrect...
const dictionary & coeffDict() const
Return const access to the coefficients dictionary.
const uniformDimensionedVectorField & g
A Vector of values with scalar precision, where scalar is float/double depending on the compilation f...
virtual bool defaultCoeffs(const bool printMsg) const
Returns true if defaultCoeffs is true and outputs on printMsg.
const dimensionedScalar mu
Atomic mass unit.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar pow3(const dimensionedScalar &ds)
Templated break-up model class.
Templated base class for dsmc cloud.
PilchErdman(const dictionary &, CloudType &)
Construct from dictionary.