|
[Sponsors] |
October 31, 2020, 02:39 |
Are sudden changes to the time step bad?
|
#1 |
Member
Join Date: Sep 2018
Posts: 53
Rep Power: 7 |
Hi all,
I heard somewhere that in transient simulations, the time step should not change too quickly. Is this true? If yes, why? For example, if a constant time step of 0.1 seconds produces results of acceptable accuracy, why would reducing the time step to 10^-6 seconds mid-run be a bad thing? Neglecting the huge increase in computation time of course. |
|
October 31, 2020, 04:40 |
|
#2 |
Senior Member
|
Assuming we are talking about methods that can actually handle that in an accurate way, I would make a distinction between methods only requiring previous the time step solution and those requiring more time levels.
I don't think the former (e.g. runge kutta) have such issues. The latter do, in the same way a sudden change in grid spacing is bad. But with time integration it is all over your grid. |
|
November 1, 2020, 06:07 |
|
#3 |
Member
Join Date: Sep 2018
Posts: 53
Rep Power: 7 |
It isn't obvious to me why multi time level schemes would suffer from a sudden change in the time step. An increase in the truncation error maybe? Could you explain a little more?
|
|
November 2, 2020, 00:23 |
|
#4 |
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,272
Rep Power: 34 |
The fluxes have time step dependency and when you suddenly change the timestep the change in flux is huge, most often you do not have number of inner iterations to handle this.
If you wish to make sudden change then at least for few time steps you would have to have very large number of inner iterations. |
|
November 2, 2020, 07:56 |
|
#5 | |
Senior Member
|
Quote:
As I said, it is in the exact same way it happens for a grid, non-uniformity introduces further errors which were not present for constant steps. Problem is that the nature of the additional terms is typically different from that of the errors also present for the uniform steps and their magnitude an order smaller. Also, the error sign depends from the change of the step (increase vs decrease). Not sure what to add, Ferziger explains the matter pretty clearly. |
||
November 2, 2020, 11:32 |
|
#6 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71 |
Quote:
Would you use a mesh having a cell of size 0.1 with an adjacent cell of size 10^-6? |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[solidMechanics] Support thread for "Solid Mechanics Solvers added to OpenFOAM Extend" | bigphil | OpenFOAM CC Toolkits for Fluid-Structure Interaction | 686 | December 22, 2022 09:10 |
pimpleDyMFoam computation randomly stops | babapeti | OpenFOAM Running, Solving & CFD | 5 | January 24, 2018 05:28 |
pressure in incompressible solvers e.g. simpleFoam | chrizzl | OpenFOAM Running, Solving & CFD | 13 | March 28, 2017 05:49 |
High Courant Number @ icoFoam | Artex85 | OpenFOAM Running, Solving & CFD | 11 | February 16, 2017 13:40 |
pisoFoam with k-epsilon turb blows up - Some questions | Heroic | OpenFOAM Running, Solving & CFD | 26 | December 17, 2012 03:34 |