33 Foam::kineticTheoryModel::kineticTheoryModel
52 kineticTheoryProperties_
56 "kineticTheoryProperties",
57 Ua_.time().constant(),
63 kineticTheory_(kineticTheoryProperties_.lookup(
"kineticTheory")),
64 equilibrium_(kineticTheoryProperties_.lookup(
"equilibrium")),
70 kineticTheoryProperties_
77 kineticTheoryProperties_
84 kineticTheoryProperties_
87 granularPressureModel_
91 kineticTheoryProperties_
94 frictionalStressModel_
98 kineticTheoryProperties_
101 e_(kineticTheoryProperties_.lookup(
"e")),
102 alphaMax_(kineticTheoryProperties_.lookup(
"alphaMax")),
103 alphaMinFriction_(kineticTheoryProperties_.lookup(
"alphaMinFriction")),
104 Fr_(kineticTheoryProperties_.lookup(
"Fr")),
105 eta_(kineticTheoryProperties_.lookup(
"eta")),
106 p_(kineticTheoryProperties_.lookup(
"p")),
219 gs0_ = radialModel_->g0
221 min(
max(alpha_, scalar(1
e-6)), alphaMax_ - 0.01),
226 volScalarField PsCoeff = granularPressureModel_->granularPressureCoeff
235 kappa_ = conductivityModel_->kappa(alpha_, Theta_, gs0_, rhoa_, da_, e_);
238 mua_ = viscosityModel_->mua(alpha_, Theta_, gs0_, rhoa_, da_, e_);
252 12.0*(1.0 -
sqr(e_))*
sqr(alpha_)*rhoa_*gs0_*(1.0/da_)*ThetaSqrt/sqrtPi;
257 0.25*
sqr(betaPrim)*da_*
sqr(Ur)
258 /(
max(alpha_, scalar(1
e-6))*rhoa_*sqrtPi*(ThetaSqrt + TsmallSqrt));
261 lambda_ = (4.0/3.0)*
sqr(alpha_)*rhoa_*da_*gs0_*(1.0+e_)*ThetaSqrt/sqrtPi;
275 +
fvm::div(phi, Theta_,
"div(phi,Theta)")
282 +
fvm::Sp(J2/(Theta_ + Tsmall), Theta_)
295 (sqrtPi/(3.0*(3.0-e_)))
296 *(1.0 + 0.4*(1.0 + e_)*(3.0*e_ - 1.0)*alpha_*gs0_)
297 +1.6*alpha_*gs0_*(1.0 + e_)/sqrtPi
301 4.0*da_*rhoa_*(1.0 + e_)*alpha_*gs0_/(3.0*sqrtPi) - 2.0*K3/3.0;
314 Theta_ =
sqr((l1 +
sqrt(l2 + l3))/(2.0*(alpha_ + 1.0
e-4)*K4));
330 PsCoeff += pf/(Theta_+Tsmall);
332 PsCoeff.
min(1.0
e+10);
333 PsCoeff.
max(-1.0
e+10);
336 pa_ = PsCoeff*Theta_;
353 Info<<
"kinTheory: max(Theta) = " <<
max(Theta_).value() <<
endl;
357 Info<<
"kinTheory: min(nua) = " <<
min(ktn).value()
358 <<
", max(nua) = " <<
max(ktn).value() <<
endl;
360 Info<<
"kinTheory: min(pa) = " <<
min(pa_).value()
361 <<
", max(pa) = " <<
max(pa_).value() <<
endl;