Controlling time step

Hello,
I’m rather new to ASPECT. I’m trying to run a simple model that simulates the development of A type Olivine undergoing simple shear, starting from zero strain and gradually increasing to a strain of 5. I attached the relevant prm file. The velocity vector (in subsection Velocity function, which is located within the subsection Prescribed Stokes solution) is defined is x10e-5;0;0; so the shear strain rate is 0.5*(10e-5) = 5*10e-6. The end time is 2e6 (seconds).
Since strain increment = strain rate * time step, I need to be able to control the time step size. If possible, I prefer to keep CFL as is. I read the documentation that there is a Time Stepping section although I haven’t been able to find it (and it also uses CFL) and I feel a bit lost in the pile of subsections and parameters. Any help would be appreciated.
Thanks!
Neta
drex.prm (5.1 KB)

Hey Neta,

Thanks for posting on the forum! Navigating the manual can certainly be overwhelming for a code as big as ASPECT :sweat_smile:. It sounds like the parameters that you’re interested in are actually Global parameters, specifically the parameters Maximum time step size and perhaps also Maximum first time step size. The CFL number parameter can also be changed here, though you shouldn’t need to adjust it if you tweak these two parameters.

Cheers,
Daniel

Hi Daniel,
Thank you very much for your quick and kind response. I feel privileged to receive a reply from someone like you! This is exactly what I’ve been looking for.
Thanks!!!
Neta

‫בתאריך יום א׳, 24 בנוב׳ 2024 ב-17:54 מאת ‪Daniel Douglas via Computational Infrastructure for Geodynamics‬‏ <‪notifications@geodynamics.discoursemail.com‬‏>:‬