35 namespace incompressible
54 RASModel(typeName, U, phi, lamTransportModel),
173 eta_(k_/epsilon_*
sqrt(2.0*
magSqr(0.5*(gradU_ + gradU_.T())))),
174 ksi_(k_/epsilon_*
sqrt(2.0*
magSqr(0.5*(gradU_ - gradU_.T())))),
175 Cmu_(2.0/(3.0*(A1_ + eta_ + alphaKsi_*ksi_))),
176 fEta_(A2_ +
pow(eta_, 3.0)),
182 magSqr(gradU_ + gradU_.T())
183 -
magSqr(gradU_ - gradU_.T())
206 pow(k_, 3.0)/
sqr(epsilon_)
211 + (gradU_ & gradU_)().
T()
213 + Ctau2_/fEta_*(gradU_ & gradU_.T())
214 + Ctau3_/fEta_*(gradU_.T() & gradU_)
217 - 20.0*
pow(k_, 4.0)/
pow(epsilon_, 3.0)
220 ((gradU_ & gradU_) & gradU_.T())
221 + ((gradU_ & gradU_.T()) & gradU_.T())
222 - ((gradU_.T() & gradU_) & gradU_)
223 - ((gradU_.T() & gradU_.T()) & gradU_)
228 nut_ = Cmu_*
sqr(k_)/(epsilon_ + epsilonSmall_) + C5viscosity_;
229 nut_.correctBoundaryConditions();
251 ((2.0/3.0)*
I)*k_ - nut_*
twoSymm(gradU_) + nonlinearStress_,
291 if (RASModel::read())
330 Cmu_*
sqr(k_)/epsilon_*S2 - (nonlinearStress_ && gradU_)
344 -
fvm::Sp(C2_*epsilon_/k_, epsilon_)
374 eta_ = k_/epsilon_*
sqrt(2.0*
magSqr(0.5*(gradU_ + gradU_.
T())));
375 ksi_ = k_/epsilon_*
sqrt(2.0*
magSqr(0.5*(gradU_ - gradU_.
T())));
376 Cmu_ = 2.0/(3.0*(A1_ + eta_ + alphaKsi_*ksi_));
377 fEta_ = A2_ +
pow(eta_, 3.0);
380 - 2.0*
pow(Cmu_, 3.0)*
pow(k_, 4.0)/
pow(epsilon_, 3.0)
383 nut_ = Cmu_*
sqr(k_)/epsilon_ + C5viscosity_;
386 nonlinearStress_ =
symm
389 pow(k_, 3.0)/
sqr(epsilon_)*
394 + (gradU_ & gradU_)().
T()
396 + Ctau2_/fEta_*(gradU_ & gradU_.T())
397 + Ctau3_/fEta_*(gradU_.T() & gradU_)
400 - 20.0*
pow(k_, 4.0)/
pow(epsilon_, 3.0)
403 ((gradU_ & gradU_) & gradU_.T())
404 + ((gradU_ & gradU_.T()) & gradU_.T())
405 - ((gradU_.T() & gradU_) & gradU_)
406 - ((gradU_.T() & gradU_.T()) & gradU_)