zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: vanDriest: change calcInterval to updateInterval

Browse Source
This commit is contained in:
mattijs
2022-12-07 11:31:18 +00:00
parent 9fbe2ef0b3
commit 3e260b84db
2 changed files with 32 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

23
src/TurbulenceModels/turbulenceModels/LES/LESdeltas/vanDriestDelta/vanDriestDelta.C
View File

@ -6,7 +6,7 @@
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
Copyright (C) 2011-2017 OpenFOAM Foundation Copyright (C) 2011-2017 OpenFOAM Foundation
Copyright (C) 2020 OpenCFD Ltd. Copyright (C) 2020,2022 OpenCFD Ltd.
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
This file is part of OpenFOAM. This file is part of OpenFOAM.
@ -133,16 +133,14 @@ Foam::LESModels::vanDriestDelta::vanDriestDelta
"Cdelta", "Cdelta",
0.158 0.158
) )
),
calcInterval_
(
dict.optionalSubDict(type() + "Coeffs").getOrDefault<label>
(
"calcInterval",
1
)
) )
{ {
const dictionary& coeffsDict(dict.optionalSubDict(type() + "Coeffs"));
if (!coeffsDict.readIfPresent("calcInterval", calcInterval_))
{
coeffsDict.readIfPresent("updateInterval", calcInterval_);
}
delta_ = geometricDelta_(); delta_ = geometricDelta_();
} }
@ -157,7 +155,10 @@ void Foam::LESModels::vanDriestDelta::read(const dictionary& dict)
dict.readIfPresent<scalar>("kappa", kappa_); dict.readIfPresent<scalar>("kappa", kappa_);
coeffsDict.readIfPresent<scalar>("Aplus", Aplus_); coeffsDict.readIfPresent<scalar>("Aplus", Aplus_);
coeffsDict.readIfPresent<scalar>("Cdelta", Cdelta_); coeffsDict.readIfPresent<scalar>("Cdelta", Cdelta_);
coeffsDict.readIfPresent<label>("calcInterval", calcInterval_); if (!coeffsDict.readIfPresent<label>("calcInterval", calcInterval_))
{
coeffsDict.readIfPresent("updateInterval", calcInterval_);
}
calcDelta(); calcDelta();
} }
@ -165,7 +166,7 @@ void Foam::LESModels::vanDriestDelta::read(const dictionary& dict)
void Foam::LESModels::vanDriestDelta::correct() void Foam::LESModels::vanDriestDelta::correct()
{ {
if (turbulenceModel_.mesh().time().timeIndex() % calcInterval_ == 0) if ((turbulenceModel_.mesh().time().timeIndex() % calcInterval_) == 0)
{ {
geometricDelta_().correct(); geometricDelta_().correct();
calcDelta(); calcDelta();

21
src/TurbulenceModels/turbulenceModels/LES/LESdeltas/vanDriestDelta/vanDriestDelta.H
View File

@ -6,7 +6,7 @@
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
Copyright (C) 2011-2015 OpenFOAM Foundation Copyright (C) 2011-2015 OpenFOAM Foundation
Copyright (C) 2019 OpenCFD Ltd. Copyright (C) 2019,2022 OpenCFD Ltd.
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
This file is part of OpenFOAM. This file is part of OpenFOAM.
@ -30,6 +30,25 @@ Class
Description Description
Simple cube-root of cell volume delta used in incompressible LES models. Simple cube-root of cell volume delta used in incompressible LES models.
Example of the delta specification in turbulenceProperties:
\verbatim
delta vanDriest;
vanDriestCoeffs
{
// Underlying geometric delta
delta cubeRootVol;
// Coefficients for geometric delta
cubeRootVolCoeffs
{
deltaCoeff 2.0;
}
// Optional update interval
// Default is 1 (update every step)
updateInterval 1000;
}
\endverbatim
SourceFiles SourceFiles
vanDriestDelta.C vanDriestDelta.C