void PrandtlDelta::calcDelta()
{
delta_ = min
(
static_cast<const volScalarField&>(geometricDelta_()),
(kappa_/Cdelta_)*wallDist(mesh_).y()
);
}
void Foam::LESModels::PrandtlDelta::calcDelta()
{
delta_ = min
(
static_cast<const volScalarField&>(geometricDelta_()),
(kappa_/Cdelta_)*wallDist::New(turbulenceModel_.mesh()).y()
);
}
void PrandtlDelta::correct()
{
geometricDelta_().correct();
if (mesh_.changing())
{
calcDelta();
}
}
void PrandtlDelta::read(const dictionary& d)
{
const dictionary& dd(d.subDict(type() + "Coeffs"));
geometricDelta_().read(dd);
d.readIfPresent<scalar>("kappa", kappa_);
dd.readIfPresent<scalar>("Cdelta", Cdelta_);
calcDelta();
}
void Foam::LESModels::PrandtlDelta::correct()
{
geometricDelta_().correct();
if (turbulenceModel_.mesh().changing())
{
calcDelta();
}
}
void Foam::LESModels::PrandtlDelta::read(const dictionary& dict)
{
const dictionary& coeffDict(dict.optionalSubDict(type() + "Coeffs"));
geometricDelta_().read(coeffDict);
dict.readIfPresent<scalar>("kappa", kappa_);
coeffDict.readIfPresent<scalar>("Cdelta", Cdelta_);
calcDelta();
}
void PrandtlDelta::read(const dictionary& d)
{
const dictionary& dd(d.subDict(type() + "Coeffs"));
geometricDelta_().read(dd);
kappa_ = dimensionedScalar(d.lookup("kappa")).value();
Cdelta_ = dimensionedScalar(dd.lookup("Cdelta")).value();
calcDelta();
}
