aspect icon indicating copy to clipboard operation
aspect copied to clipboard

Limit the growth of the time step size

Open bangerth opened this issue 3 years ago • 1 comments

@tjhei dug up a masters thesis at https://dc.uwm.edu/cgi/viewcontent.cgi?article=1405&context=etd that is about time step choices for BDF2. Among the conclusions there is that one should not let the time step size grow by too much from one time step to another. They call this grow the "amplification factor", and reference a number of papers that apparently show that it should not grow by more than either a factor of 1+sqrt(2) or 1.86 or 1.91 (pages 15 and 16), depending on who you trust. We should probably incorporate a limit in our code base to that end as well.

bangerth avatar Feb 18 '22 15:02 bangerth

We can basically just change the default of https://github.com/geodynamics/aspect/blob/edad8df4dacd0d43afbd8d57e0914d683630cea4/source/simulator/parameters.cc#L170-L179

tjhei avatar Feb 18 '22 16:02 tjhei