Preprints
https://doi.org/10.5194/egusphere-2025-2417
https://doi.org/10.5194/egusphere-2025-2417
12 Aug 2025
 | 12 Aug 2025
Status: this preprint is open for discussion and under review for Geoscientific Model Development (GMD).

Overcoming the numerical challenges owing to rapid ductile localization

Arne Spang, Marcel Thielmann, Casper Pranger, Albert de Montserrat, and Ludovic Räss

Abstract. Strain localization is among the most challenging mechanical phenomena for computational Earth sciences. Accurately capturing it is difficult because strain localization initiates spontaneously, is self-accelerating, and its characteristic length and time scales are typically significantly smaller than the spatial and temporal resolutions of the model. This results in an undesirable dependence of the model behavior on numerical parameters and a large computational cost. Strain localization is most commonly associated with brittle failure, but ductile processes such as thermal runaway can also result in rapid ductile localization. Here, we present a numerical model to investigate thermal runaway, and further propose strategies to overcome the challenges associated with resolving rapid localization: (i) adaptive time stepping; (ii) adaptive rescaling; and (iii) two types of regularization. We demonstrate the effect of these strategies in one- and two-dimensional models. We rely on the accelerated pseudo-transient method to solve the governing equations and use graphics processing units to accelerate two-dimensional computations. Our adaptive time stepping strategy allows us to accurately capture spontaneous and rapid stress release during thermal runaway while reducing time steps by more than ten orders of magnitude. Adaptive rescaling further reduces rounding errors and the number of required iterations by two orders of magnitude. Viscosity regularization and gradient regularization enable us to mitigate resolution dependencies but may differently impact the physical response of the model.

Competing interests: At least one of the (co-)authors is a member of the editorial board of Geoscientific Model Development.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Share
Arne Spang, Marcel Thielmann, Casper Pranger, Albert de Montserrat, and Ludovic Räss

Status: open (until 08 Oct 2025)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
Arne Spang, Marcel Thielmann, Casper Pranger, Albert de Montserrat, and Ludovic Räss
Arne Spang, Marcel Thielmann, Casper Pranger, Albert de Montserrat, and Ludovic Räss

Viewed

Total article views: 462 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
450 8 4 462 13 8
  • HTML: 450
  • PDF: 8
  • XML: 4
  • Total: 462
  • BibTeX: 13
  • EndNote: 8
Views and downloads (calculated since 12 Aug 2025)
Cumulative views and downloads (calculated since 12 Aug 2025)

Viewed (geographical distribution)

Total article views: 447 (including HTML, PDF, and XML) Thereof 447 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
1
 
 
 
 
Latest update: 03 Sep 2025
Download
Short summary
Concentration of deformation is difficult to capture accurately in computer simulations. We present a number of challenges associated with concentrated viscous deformation and demonstrate strategies to overcome them. These strategies include automatic selection of appropriate time steps to react to rapid changes in model behavior, automatic rescaling to avoid rounding errors, and two methods to prevent model instability. This way, we are able to accurately capture very fast viscous deformation.
Share