34 namespace incompressible
53 RASModel(typeName, U, phi, lamTransportModel),
192 eta_(k_/epsilon_*
sqrt(2.0*
magSqr(0.5*(gradU_ + gradU_.T())))),
193 ksi_(k_/epsilon_*
sqrt(2.0*
magSqr(0.5*(gradU_ - gradU_.T())))),
194 Cmu_(2.0/(3.0*(A1_ + eta_ + alphaKsi_*ksi_))),
195 fEta_(A2_ +
pow(eta_, 3.0)),
197 nut_(
"nut", Cmu_*
sqr(k_)/(epsilon_ + epsilonSmall_)),
204 pow(k_, 3.0)/
sqr(epsilon_)
209 + (gradU_ & gradU_)().
T()
211 + Ctau2_/fEta_*(gradU_ & gradU_.T())
212 + Ctau3_/fEta_*(gradU_.T() & gradU_)
241 ((2.0/3.0)*
I)*k_ - nut_*
twoSymm(gradU_) + nonlinearStress_,
281 if (RASModel::read())
323 Cmu_*
sqr(k_)/epsilon_*S2
324 - (nonlinearStress_ && gradU_)
337 -
fvm::Sp(C2_*epsilon_/k_, epsilon_)
367 eta_ = k_/epsilon_*
sqrt(2.0*
magSqr(0.5*(gradU_ + gradU_.T())));
368 ksi_ = k_/epsilon_*
sqrt(2.0*
magSqr(0.5*(gradU_ - gradU_.T())));
369 Cmu_ = 2.0/(3.0*(A1_ + eta_ + alphaKsi_*ksi_));
370 fEta_ = A2_ +
pow(eta_, 3.0);
372 nut_ = Cmu_*
sqr(k_)/epsilon_;
376 nonlinearStress_ =
symm
378 pow(k_, 3.0)/
sqr(epsilon_)
383 + (gradU_ & gradU_)().
T()
385 + Ctau2_/fEta_*(gradU_ & gradU_.T())
386 + Ctau3_/fEta_*(gradU_.T() & gradU_)