36 namespace incompressible
55 RASModel(typeName, U, phi, lamTransportModel),
222 eta_(k_/epsilon_*
sqrt(2.0*
magSqr(0.5*(gradU_ + gradU_.T())))),
223 ksi_(k_/epsilon_*
sqrt(2.0*
magSqr(0.5*(gradU_ - gradU_.T())))),
224 Cmu_(2.0/(3.0*(A1_ + eta_ + alphaKsi_*ksi_))),
225 fEta_(A2_ +
pow(eta_, 3.0)),
229 -2.0*
pow(Cmu_, 3.0)*
pow(k_, 4.0)/
pow(epsilon_, 3.0)
230 *(
magSqr(gradU_ + gradU_.T()) -
magSqr(gradU_ - gradU_.T()))
233 yStar_(
sqrt(k_)*y_/
nu() + SMALL),
254 pow(k_, 3.0)/
sqr(epsilon_)
259 + (gradU_ & gradU_)().
T()
261 + Ctau2_/fEta_*(gradU_ & gradU_.T())
262 + Ctau3_/fEta_*(gradU_.T() & gradU_)
265 - 20.0*
pow(k_, 4.0)/
pow(epsilon_, 3.0)
268 ((gradU_ & gradU_) & gradU_.T())
269 + ((gradU_ & gradU_.T()) & gradU_.T())
270 - ((gradU_.T() & gradU_) & gradU_)
271 - ((gradU_.T() & gradU_.T()) & gradU_)
284 scalar(1) -
exp(-Am_*yStar_))
285 /(scalar(1) -
exp(-Aepsilon_*yStar_) + SMALL
287 *
sqr(k_)/(epsilon_ + epsilonSmall_)
295 nut_.correctBoundaryConditions();
317 ((2.0/3.0)*
I)*k_ - nut_*
twoSymm(gradU_) + nonlinearStress_,
357 if (RASModel::read())
403 yStar_ =
sqrt(k_)*y_/
nu() + SMALL;
407 (scalar(1) -
exp(-Am_*yStar_))
408 /(scalar(1) -
exp(-Aepsilon_*yStar_) + SMALL);
413 Cmu_*fMu*
sqr(k_)/epsilon_*S2 - (nonlinearStress_ && gradU_);
424 + C2_*f2*
pow(Cmu_, 0.75)*
pow(k_, scalar(0.5))
425 /(kappa_*y_*(scalar(1) -
exp(-Aepsilon_*yStar_)))
426 *
exp(-Amu_*
sqr(yStar_))*epsilon_
427 -
fvm::Sp(C2_*f2*epsilon_/k_, epsilon_)
458 eta_ = k_/epsilon_*
sqrt(2.0*
magSqr(0.5*(gradU_ + gradU_.
T())));
459 ksi_ = k_/epsilon_*
sqrt(2.0*
magSqr(0.5*(gradU_ - gradU_.
T())));
460 Cmu_ = 2.0/(3.0*(A1_ + eta_ + alphaKsi_*ksi_));
461 fEta_ = A2_ +
pow(eta_, 3.0);
464 - 2.0*
pow(Cmu_, 3.0)*
pow(k_, 4.0)/
pow(epsilon_, 3.0)
468 Cmu_*fMu*
sqr(k_)/epsilon_
476 nonlinearStress_ =
symm
479 pow(k_, 3.0)/
sqr(epsilon_)
484 + (gradU_ & gradU_)().
T()
486 + Ctau2_/fEta_*(gradU_ & gradU_.T())
487 + Ctau3_/fEta_*(gradU_.T() & gradU_)
490 - 20.0*
pow(k_, 4.0)/
pow(epsilon_, 3.0)
493 ((gradU_ & gradU_) & gradU_.T())
494 + ((gradU_ & gradU_.T()) & gradU_.T())
495 - ((gradU_.T() & gradU_) & gradU_)
496 - ((gradU_.T() & gradU_.T()) & gradU_)