Skip to content
Snippets Groups Projects
Commit 57d7158f authored by Andrew Heather's avatar Andrew Heather
Browse files

added calcInterval - defaults to 1 if not present to recapture earlier behaviour

parent 105e99da
Branches
Tags
No related merge requests found
...@@ -115,6 +115,10 @@ vanDriestDelta::vanDriestDelta ...@@ -115,6 +115,10 @@ vanDriestDelta::vanDriestDelta
Cdelta_ Cdelta_
( (
dd.subDict(type() + "Coeffs").lookupOrDefault<scalar>("Cdelta", 0.158) dd.subDict(type() + "Coeffs").lookupOrDefault<scalar>("Cdelta", 0.158)
),
calcInterval_
(
dd.subDict(type() + "Coeffs").lookupOrDefault<label>("calcInterval", 1)
) )
{ {
delta_ = geometricDelta_(); delta_ = geometricDelta_();
...@@ -131,14 +135,18 @@ void vanDriestDelta::read(const dictionary& d) ...@@ -131,14 +135,18 @@ void vanDriestDelta::read(const dictionary& d)
d.readIfPresent<scalar>("kappa", kappa_); d.readIfPresent<scalar>("kappa", kappa_);
dd.readIfPresent<scalar>("Aplus", Aplus_); dd.readIfPresent<scalar>("Aplus", Aplus_);
dd.readIfPresent<scalar>("Cdelta", Cdelta_); dd.readIfPresent<scalar>("Cdelta", Cdelta_);
dd.readIfPresent<label>("calcInterval", calcInterval_);
calcDelta(); calcDelta();
} }
void vanDriestDelta::correct() void vanDriestDelta::correct()
{ {
geometricDelta_().correct(); if (mesh().time().timeIndex() % calcInterval_ == 0)
calcDelta(); {
geometricDelta_().correct();
calcDelta();
}
} }
......
...@@ -61,6 +61,7 @@ class vanDriestDelta ...@@ -61,6 +61,7 @@ class vanDriestDelta
scalar kappa_; scalar kappa_;
scalar Aplus_; scalar Aplus_;
scalar Cdelta_; scalar Cdelta_;
label calcInterval_;
// Private Member Functions // Private Member Functions
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment