zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

setDeltaT: Corrected time-step increase limiter

Browse Source
This commit is contained in:
Henry Weller
2022-06-21 15:26:36 +01:00
parent eeccdceb26
commit 0f8a020c70
5 changed files with 25 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
applications/solvers/combustion/XiFoam/PDRFoam/setDeltaT.H
View File

@ -35,7 +35,11 @@ if (adjustTimeStep)
{
scalar deltaT = maxCo*runTime.deltaTValue()/(CoNum + StCoNum + small);
deltaT = min(deltaT, fvModels.maxDeltaT());
deltaT = min(min(deltaT, 1.0 + 0.1*deltaT), 1.2*runTime.deltaTValue());
deltaT = min
(
min(deltaT, runTime.deltaTValue() + 0.1*deltaT),
1.2*runTime.deltaTValue()
);
runTime.setDeltaT(min(deltaT, maxDeltaT));
Info<< "deltaT = " << runTime.deltaTValue() << endl;

6
applications/solvers/compressible/rhoCentralFoam/setDeltaT.H
View File

@ -34,7 +34,11 @@ Description
if (adjustTimeStep)
{
scalar deltaT = maxCo*runTime.deltaTValue()/(CoNum + small);
deltaT = min(min(deltaT, 1.0 + 0.1*deltaT), 1.2*runTime.deltaTValue());
deltaT = min
(
min(deltaT, runTime.deltaTValue() + 0.1*deltaT),
1.2*runTime.deltaTValue()
);
runTime.setDeltaT(min(deltaT, maxDeltaT));
Info<< "deltaT = " << runTime.deltaTValue() << endl;

6
applications/solvers/multiphase/twoLiquidMixingFoam/setDeltaT.H
View File

@ -36,7 +36,11 @@ if (adjustTimeStep)
scalar deltaT =
min(maxCo/(CoNum + small), maxAlphaCo/(alphaCoNum + small))
*runTime.deltaTValue();
deltaT = min(min(deltaT, 1.0 + 0.1*deltaT), 1.2*runTime.deltaTValue());
deltaT = min
(
min(deltaT, runTime.deltaTValue() + 0.1*deltaT),
1.2*runTime.deltaTValue()
);
runTime.setDeltaT(min(deltaT, maxDeltaT));
Info<< "deltaT = " << runTime.deltaTValue() << endl;

6
src/finiteVolume/cfdTools/general/include/setDeltaT.H
View File

@ -35,7 +35,11 @@ if (adjustTimeStep)
{
scalar deltaT = maxCo*runTime.deltaTValue()/(CoNum + small);
deltaT = min(deltaT, fvModels.maxDeltaT());
deltaT = min(min(deltaT, 1.0 + 0.1*deltaT), 1.2*runTime.deltaTValue());
deltaT = min
(
min(deltaT, runTime.deltaTValue() + 0.1*deltaT),
1.2*runTime.deltaTValue()
);
runTime.setDeltaT(min(deltaT, maxDeltaT));
Info<< "deltaT = " << runTime.deltaTValue() << endl;

6
src/twoPhaseModels/twoPhaseMixture/VoF/setDeltaT.H
View File

@ -37,7 +37,11 @@ if (adjustTimeStep)
min(maxCo/(CoNum + small), maxAlphaCo/(alphaCoNum + small))
*runTime.deltaTValue();
deltaT = min(deltaT, fvModels.maxDeltaT());
deltaT = min(min(deltaT, 1.0 + 0.1*deltaT), 1.2*runTime.deltaTValue());
deltaT = min
(
min(deltaT, runTime.deltaTValue() + 0.1*deltaT),
1.2*runTime.deltaTValue()
);
runTime.setDeltaT(min(deltaT, maxDeltaT));
Info<< "deltaT = " << runTime.deltaTValue() << endl;