33 namespace compressible
46 void Smagorinsky::updateSubGridScaleFields(
const volTensorField& gradU)
56 muSgs_ = ck_*
rho()*delta()*
sqrt(k_);
57 muSgs_.correctBoundaryConditions();
59 alphaSgs_ = muSgs_/Prt_;
60 alphaSgs_.correctBoundaryConditions();
66 Smagorinsky::Smagorinsky
74 LESModel(typeName, rho, U, phi, thermoPhysicalModel),
98 updateSubGridScaleFields(gradU());
102 bool Smagorinsky::read()
104 if (GenEddyVisc::read())
106 ck_.readIfPresent(coeffDict());