33template<
class CloudType>
52template<
class CloudType>
63template<
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);
const uniformDimensionedVectorField & g
Templated break-up model class.
Templated base class for dsmc cloud.
Particle secondary breakup model based on Pilch-Erdman total droplet breakup model.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
bool readEntry(const word &keyword, T &val, enum keyType::option matchOpt=keyType::REGEX, bool mandatory=true) const
virtual bool update()
Update the mesh for both mesh motion and topology change.
const dictionary & coeffDict() const
Return const access to the coefficients dictionary.
virtual bool defaultCoeffs(const bool printMsg) const
Returns true if defaultCoeffs is true and outputs on printMsg.
const volScalarField & mu
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 pow3(const dimensionedScalar &ds)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar sqrt(const dimensionedScalar &ds)