diff --git a/applications/solvers/combustion/XiFoam/PDRFoam/setDeltaT.H b/applications/solvers/combustion/XiFoam/PDRFoam/setDeltaT.H index 98682b6af6..862cf015a8 100644 --- a/applications/solvers/combustion/XiFoam/PDRFoam/setDeltaT.H +++ b/applications/solvers/combustion/XiFoam/PDRFoam/setDeltaT.H @@ -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; diff --git a/applications/solvers/compressible/rhoCentralFoam/setDeltaT.H b/applications/solvers/compressible/rhoCentralFoam/setDeltaT.H index 235313dda6..43705f1be0 100644 --- a/applications/solvers/compressible/rhoCentralFoam/setDeltaT.H +++ b/applications/solvers/compressible/rhoCentralFoam/setDeltaT.H @@ -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; diff --git a/applications/solvers/multiphase/twoLiquidMixingFoam/setDeltaT.H b/applications/solvers/multiphase/twoLiquidMixingFoam/setDeltaT.H index 7b683fa736..7441fee79f 100644 --- a/applications/solvers/multiphase/twoLiquidMixingFoam/setDeltaT.H +++ b/applications/solvers/multiphase/twoLiquidMixingFoam/setDeltaT.H @@ -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; diff --git a/src/finiteVolume/cfdTools/general/include/setDeltaT.H b/src/finiteVolume/cfdTools/general/include/setDeltaT.H index e17ce9ccef..289400edd0 100644 --- a/src/finiteVolume/cfdTools/general/include/setDeltaT.H +++ b/src/finiteVolume/cfdTools/general/include/setDeltaT.H @@ -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; diff --git a/src/twoPhaseModels/twoPhaseMixture/VoF/setDeltaT.H b/src/twoPhaseModels/twoPhaseMixture/VoF/setDeltaT.H index 1e943bff9d..81741d7a8a 100644 --- a/src/twoPhaseModels/twoPhaseMixture/VoF/setDeltaT.H +++ b/src/twoPhaseModels/twoPhaseMixture/VoF/setDeltaT.H @@ -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;