Articles | Volume 19, issue 12
https://doi.org/10.5194/gmd-19-5343-2026
© Author(s) 2026. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Special issue:
https://doi.org/10.5194/gmd-19-5343-2026
© Author(s) 2026. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Automatic tuning of iterative pseudo-transient solvers for modeling the deformation of heterogeneous media
Thibault Duretz
CORRESPONDING AUTHOR
Institut für Geowissenschaften, Goethe‐Universität Frankfurt, Frankfurt, Germany
Albert de Montserrat
Geophysical Fluid Dynamics, Institute of Geophysics, ETH Zurich, Zurich, Switzerland
Rubén Sevilla
Zienkiewicz Centre for Computational Engineering, Faculty of Science and Engineering, Swansea University, Wales, UK
Ludovic Räss
Institute of Earth sciences, University of Lausanne, Lausanne, Switzerland
Ivan Utkin
Laboratory of Hydraulics, Hydrology and Glaciology (VAW), ETH Zurich, Zurich, Switzerland
Arne Spang
Bayerisches Geoinstitut, Universität Bayreuth, Bayreuth, Germany
Related authors
No articles found.
Arne Spang, Marcel Thielmann, Casper Pranger, Albert de Montserrat, and Ludovic Räss
Geosci. Model Dev., 19, 369–388, https://doi.org/10.5194/gmd-19-369-2026, https://doi.org/10.5194/gmd-19-369-2026, 2026
Short summary
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. The strategies include automatic selection of appropriate time steps to react to rapid changes in model behavior, automatic rescaling to avoid rounding errors, and three methods to prevent model instability. This way, we are able to accurately capture very fast viscous deformation.
Nicolas Riel, Boris J. P. Kaus, Albert de Montserrat, Evangelos Moulas, Eleanor C. R. Green, and Hugo Dominguez
Geosci. Model Dev., 18, 6951–6962, https://doi.org/10.5194/gmd-18-6951-2025, https://doi.org/10.5194/gmd-18-6951-2025, 2025
Short summary
Short summary
Our research focuses on improving the way we predict mineral assemblage. Current methods, while accurate, are slowed by complex calculations. We developed a new approach that simplifies these calculations and speeds them up significantly using a technique called the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm. This breakthrough reduces computation time by more than five times, potentially unlocking new horizons in modeling reactive magmatic systems.
Manuele Faccenda, Brandon P. VanderBeek, Albert de Montserrat, Jianfeng Yang, Francesco Rappisi, and Neil Ribe
Solid Earth, 15, 1241–1264, https://doi.org/10.5194/se-15-1241-2024, https://doi.org/10.5194/se-15-1241-2024, 2024
Short summary
Short summary
The Earth's internal dynamics and structure can be well understood by combining seismological and geodynamic modelling with mineral physics, an approach that has been poorly adopted in the past. To this end we have developed ECOMAN, an open-source software package that is intended to overcome the computationally intensive nature of this multidisciplinary methodology and the lack of a dedicated and comprehensive computational framework.
Anjali Sandip, Ludovic Räss, and Mathieu Morlighem
Geosci. Model Dev., 17, 899–909, https://doi.org/10.5194/gmd-17-899-2024, https://doi.org/10.5194/gmd-17-899-2024, 2024
Short summary
Short summary
We solve momentum balance for unstructured meshes to predict ice flow for real glaciers using a pseudo-transient method on graphics processing units (GPUs) and compare it to a standard central processing unit (CPU) implementation. We justify the GPU implementation by applying the price-to-performance metric for up to million-grid-point spatial resolutions. This study represents a first step toward leveraging GPU processing power, enabling more accurate polar ice discharge predictions.
Ludovic Räss, Ivan Utkin, Thibault Duretz, Samuel Omlin, and Yuri Y. Podladchikov
Geosci. Model Dev., 15, 5757–5786, https://doi.org/10.5194/gmd-15-5757-2022, https://doi.org/10.5194/gmd-15-5757-2022, 2022
Short summary
Short summary
Continuum mechanics-based modelling of physical processes at large scale requires huge computational resources provided by massively parallel hardware such as graphical processing units. We present a suite of numerical algorithms, implemented using the Julia language, that efficiently leverages the parallelism. We demonstrate that our implementation is efficient, scalable and robust and showcase applications to various geophysical problems.
Cited articles
Alkhimenkov, Y., Räss, L., Khakimova, L., Quintal, B., and Podladchikov, Y.: Resolving Wave Propagation in Anisotropic Poroelastic Media Using Graphical Processing Units (GPUs), J. Geophys. Res.-Sol. Ea., 126, https://doi.org/10.1029/2020JB021175, 2021. a
Anzt, H., Chow, E., Saak, J., and Dongarra, J.: Updating incomplete factorization preconditioners for model order reduction, Numer. Algorithms, 73, 611–630, https://doi.org/10.1007/s11075-016-0110-2, 2016. a
Balachandar, S. and Yuen, D. A.: Three-dimensional fully spectral numerical method for mantle convection with depth-dependent properties, J. Comput. Phys., 113, https://doi.org/10.1006/jcph.1994.1118, 1994. a
Baumgardner, J. R.: Three-dimensional treatment of convective flow in the earth's mantle, J. Stat. Phys., 39, 501–511, https://doi.org/10.1007/BF01008348, 1985. a
Bezanson, J., Edelman, A., Karpinski, S., and Shah, V. B.: Julia: A Fresh Approach to Numerical Computing, SIAM Rev., 59, 65–98, https://doi.org/10.1137/141000671, 2017. a
Burcet, M., Oliveira, B., Afonso, J. C., and Zlotnik, S.: A face-centred finite volume approach for coupled transport phenomena and fluid flow, Appl. Math. Model., 125, 293–312, https://doi.org/10.1016/j.apm.2023.08.031, 2024. a
Chow, E., Anzt, H., Scott, J., and Dongarra, J.: Using Jacobi iterations and blocking for solving sparse triangular systems in incomplete factorization preconditioning, J. Parallel Distr. Com., 119, 219–230, https://doi.org/10.1016/j.jpdc.2018.04.017, 2018. a
Cockburn, B. and Gopalakrishnan, J.: New hybridization techniques, GAMM-Mitteilungen, 28, 154–182, https://doi.org/10.1002/gamm.201490017, 2005. a
Cundall, P. and Board, M.: A microcomputer program for modelling large-strain plasticity problems, Proceedings of the Sixth International Conference on Numerical Methods in Geomechanics, 11–15 April 1988, Innsbruck, Austria, Vols. 1–3, ISBN 978-90-6191-80, 1988. a
Dabrowski, M., Krotkiewski, M., and Schmid, D. W.: MILAMIN: MATLAB-based finite element method solver for large problems, Geochem. Geophy. Geosy., 9, https://doi.org/10.1029/2007GC001719, 2008. a, b
de Borst, R. and Duretz, T.: On viscoplastic regularisation of strain-softening rocks and soils, Int. J. Numer. Anal. Met., 44, 890–903, https://doi.org/10.1002/nag.3046, 2020. a
de la Puente, J., Ferrer, M., Hanzich, M., Castillo, J. E., and Cela, J. M.: Mimetic seismic wave modeling including topography on deformed staggered grids, Geophysics, https://doi.org/10.1190/geo2013-0371.1, 2014. a
de Montserrat, A., Aellig, P. S., Schuler, C., Navarrete, I., Räss, L., Fuchs, L., Kaus, B. J. P., and Dominguez, H.: JustRelax.jl: A Julia package for geodynamic modeling with matrix-free solvers, Journal of Open Source Software, 11, 9365, https://doi.org/10.21105/joss.09365, 2026. a
Duretz, T., May, D. A., Gerya, T. V., and Tackley, P. J.: Discretization errors and free surface stabilization in the finite difference and marker-in-cell method for applied geodynamics: A numerical study, Geochem. Geophy. Geosy., 12, https://doi.org/10.1029/2011GC003567, 2011. a, b, c
Duretz, T., May, D. A., and Yamato, P.: A free surface capturing discretization for the staggered grid finite difference scheme, Geophys. J. Int., 204, 1518–1530, https://doi.org/10.1093/gji/ggv526, 2016. a
Duretz, T., Räss, L., Podladchikov, Y., and Schmalholz, S.: Resolving thermomechanical coupling in two and three dimensions: spontaneous strain localization owing to shear heating, Geophys. J. Int., 216, 365–379, https://doi.org/10.1093/gji/ggy434, 2019. a, b, c
Duretz, T., de Borst, R., Yamato, P., and Le Pourhiet, L.: Toward Robust and Predictive Geodynamic Modeling: The Way Forward in Frictional Plasticity, Geophys. Res. Lett., 47, e2019GL086027, https://doi.org/10.1029/2019GL086027, 2020. a
Feng, Y. T.: On the discrete dynamic nature of the conjugate gradient method, J. Comput. Phys., 211, 91–98, https://doi.org/10.1016/j.jcp.2005.05.006, 2006. a
Fu, Z., James Lewis, T., Kirby, R. M., and Whitaker, R. T.: Architecting the finite element method pipeline for the GPU, J. Comput. Appl. Math., 257, 195–211, https://doi.org/10.1016/j.cam.2013.09.001, 2014. a
Fuchs, L. and Becker, T. W.: On the Role of Rheological Memory for Convection-Driven Plate Reorganizations, Geophys. Res. Lett., 49, e2022GL099574, https://doi.org/10.1029/2022GL099574, 2022. a
Fullsack, P.: An arbitrary Lagrangian-Eulerian formulation for creeping flows and its application in tectonic models, Geophys. J. Int., 120, 1–23, https://doi.org/10.1111/j.1365-246X.1995.tb05908.x, 1995. a
Furuichi, M., May, D. A., and Tackley, P. J.: Development of a Stokes flow solver robust to large viscosity jumps using a Schur complement approach with mixed precision arithmetic, J. Comput. Phys., 230, 8835–8851, https://doi.org/10.1016/j.jcp.2011.09.007, 2011. a
Gerya, T.: Introduction to Numerical Geodynamic Modelling, Cambridge University Press, Cambridge, England, UK, ISBN 978-1-10714314-2, 2019. a
Gerya, T. V.: Initiation of transform faults at rifted continental margins: 3D petrological-thermomechanical modeling and comparison to the Woodlark Basin, Petrology, 21, 550–560, https://doi.org/10.1134/S0869591113060039, 2013. a, b
Gerya, T. V. and Yuen, D. A.: Characteristics-based marker-in-cell method with conservative finite-differences schemes for modeling geological flows with strongly variable transport properties, Phys. Earth Planet. In., 140, 293–318, https://doi.org/10.1016/j.pepi.2003.09.006, 2003. a, b
Giacomini, M. and Sevilla, R.: A second-order face-centred finite volume method on general meshes with automatic mesh adaptation, Int. J. Numer. Meth. Eng., 121, 5227–5255, https://doi.org/10.1002/nme.6428, 2020. a
Goldade, R., Wang, Y., Aanjaneya, M., and Batty, C.: An adaptive variational finite difference framework for efficient symmetric octree viscosity, ACM Trans. Graph., 38, 1–14, https://doi.org/10.1145/3306346.3322939, 2019. a
Gray, T. S., Tackley, P. J., and Gerya, T.: Variational Stokes method applied to free surface boundaries in numerical geodynamic models, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-6354, 2025. a
Halter, W. R., Macherel, E., and Schmalholz, S. M.: A simple computer program for calculating stress and strain rate in 2D viscous inclusion-matrix systems, J. Struct. Geol., 160, 104617, https://doi.org/10.1016/j.jsg.2022.104617, 2022. a, b, c
Harlow, F. H. and Welch, J. E.: Numerical Calculation of Time-Dependent Viscous Incompressible Flow of Fluid with Free Surface, Phys. Fluids, 8, 2182–2189, https://doi.org/10.1063/1.1761178, 1965. a
Hassani, R., Jongmans, D., and Chéry, J.: Study of plate deformation and stress in subduction processes using two-dimensional numerical models, J. Geophys. Res.-Sol. Ea., 102, 17951–17965, https://doi.org/10.1029/97JB01354, 1997. a
Hillebrand, B., Thieulot, C., Geenen, T., van den Berg, A. P., and Spakman, W.: Using the level set method in geodynamical modeling of multi-material flows and Earth's free surface, Solid Earth, 5, 1087–1098, https://doi.org/10.5194/se-5-1087-2014, 2014. a
Joldes, G. R., Wittek, A., and Miller, K.: An adaptive Dynamic Relaxation method for solving nonlinear finite element problems. Application to brain shift estimation, Int. J. Numer. Meth. Bio., 27, 173–185, https://doi.org/10.1002/cnm.1407, 2011. a, b
Kaus, B. J. P., Popov, A. A., Baumann, T. S., Püsök, A. E., Bauville, A., Fernandez, N., and Collignon, M.: Forward and Inverse Modelling of Lithospheric Deformation on Geological Timescales, in: NIC Symposium 2016, edited by: Binder, K., Müller, M., Kremer, A., and Schnurpfeil, A., Vol. 48, 299–307, Forschungszentrum Jülich, Jülich, ISBN 978-3-95806-109-5, 2016. a
Kiran, U., Sharma, D., and Gautam, S. S.: GPU-warp based finite element matrices generation and assembly using coloring method, J. Comput. Des. Eng., 6, 705–718, https://doi.org/10.1016/j.jcde.2018.11.001, 2019. a
Kronbichler, M., Heister, T., and Bangerth, W.: High Accuracy Mantle Convection Simulation through Modern Numerical Methods, Geophys. J. Int., 191, 12–29, https://doi.org/10.1111/j.1365-246X.2012.05609.x, 2012. a
Larionov, E., Batty, C., and Bridson, R.: Variational stokes: a unified pressure-viscosity solver for accurate viscous liquids, ACM T. Graphic., 36, 1–11, https://doi.org/10.1145/3072959.3073628, 2017. a
Luisier, C., Tajčmanová, L., Yamato, P., and Duretz, T.: Garnet microstructures suggest ultra-fast decompression of ultrahigh-pressure rocks, Nat. Commun., 14, 1–8, https://doi.org/10.1038/s41467-023-41310-w, 2023. a
Macherel, E., Räss, L., and Schmalholz, S. M.: 3D Stresses and Velocities Caused by Continental Plateaus: Scaling Analysis and Numerical Calculations With Application to the Tibetan Plateau, Geochem. Geophy. Geosy., 25, e2023GC011356, https://doi.org/10.1029/2023GC011356, 2024. a
May, D., Brown, J., and Le Pourhiet, L.: A scalable, matrix-free multigrid preconditioner for finite element discretizations of heterogeneous Stokes flow, Comput. Method. Appl. M., 290, 496–523, https://doi.org/10.1016/j.cma.2015.03.014, 2015. a, b
May, D. A. and Moresi, L.: Preconditioned iterative methods for Stokes flow problems arising in computational geodynamics, Phys. Earth Planet. In., 171, 33–47, https://doi.org/10.1016/j.pepi.2008.07.036, 2008. a
May, D. A., Brown, J., and Pourhiet, L. L.: pTatin3D: High-Performance Methods for Long-Term Lithospheric Dynamics, in: SC'14: Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis, 274–284, https://doi.org/10.1109/SC.2014.28, 2014. a
Omlin, S., Räss, L., and Utkin, I.: Distributed Parallelization of xPU Stencil Computations in Julia, Proceedings of the JuliaCon Conferences, 137 pp., https://doi.org/10.21105/jcon.00137, 2022. a, b
Papadrakakis, M.: A method for the automatic evaluation of the dynamic relaxation parameters, Comput. Method. Appl. M., 25, 35–48, https://doi.org/10.1016/0045-7825(81)90066-9, 1981. a, b, c
Poliakov, A. N. B., Cundall, P. A., Podladchikov, Y. Y., and Lyakhovsky, V. A.: An Explicit Inertial Method for the Simulation of Viscoelastic Flow: An Evaluation of Elastic Effects on Diapiric Flow in Two- and Three- Layers Models, in: Flow and Creep in the Solar System: Observations, Modeling and Theory, 175–195, Springer, Dordrecht, the Netherlands, ISBN 978-94-015-8206-3, https://doi.org/10.1007/978-94-015-8206-3_12, 1993. a
Popov, A. A. and Sobolev, S. V.: SLIM3D: A Tool for Three-dimensional Thermomechanical Modeling of Lithospheric Deformation with Elasto-visco-plastic Rheology, Phys. Earth Planet. In., 171, 55–75, https://doi.org/10.1016/j.pepi.2008.03.007, 2008. a
Räss, L. and Duretz, T.: PTsolvers/AutoTuningGMD: Release 1.0, Zenodo [code], https://doi.org/10.5281/zenodo.20269899, 2025. a
Räss, L., Duretz, T., Podladchikov, Y. Y., and Schmalholz, S. M.: M2Di: Concise and efficient MATLAB 2-D Stokes solvers using the Finite Difference Method, Geochem. Geophy. Geosy., 18, 755–768, https://doi.org/10.1002/2016GC006727, 2017. a, b, c
Räss, L., Duretz, T., and Podladchikov, Y. Y.: Resolving hydro-mechanical coupling in two and three dimensions: Spontaneous channelling of porous fluids owing to decompaction weakening, Geophys. J. Int., https://doi.org/10.1093/gji/ggz239, 2019. a
Rezaiee-Pajand, M. and Sarafrazi, S. R.: Nonlinear dynamic structural analysis using dynamic relaxation with zero damping, Comput. Struct., 89, 1274–1285, https://doi.org/10.1016/j.compstruc.2011.04.005, 2011. a
Rudi, J., Stadler, G., and Ghattas, O.: Weighted BFBT Preconditioner for Stokes Flow Problems with Highly Heterogeneous Viscosity, SIAM J. Sci. Comput., 39, S272–S297, https://doi.org/10.1137/16M108450X, 2017. a
Samuel, H. and Evonuk, M.: Modeling advection in geophysical flows with particle level sets, Geochem. Geophy. Geosy., 11, Q08020, https://doi.org/10.1029/2010GC003081, 2010. a
Sanan, P., May, D. A., Bollhöfer, M., and Schenk, O.: Pragmatic solvers for 3D Stokes and elasticity problems with heterogeneous coefficients: evaluating modern incomplete LDLT preconditioners, Solid Earth, 11, 2031–2045, https://doi.org/10.5194/se-11-2031-2020, 2020. a
Schmalholz, S. M. and Podladchikov, Y. Yu.: Strain and competence contrast estimation from fold shape, Tectonophysics, 340, 195–213, https://doi.org/10.1016/S0040-1951(01)00151-2, 2001. a
Schmeling, H., Babeyko, A. Y., Enns, A., Faccenna, C., Funiciello, F., Gerya, T., Golabek, G. J., Grigull, S., Kaus, B. J. P., Morra, G., Schmalholz, S. M., and van Hunen, J.: A benchmark comparison of spontaneous subduction models – Towards a free surface, Phys. Earth Planet. Int., 171, 198–223, https://doi.org/10.1016/j.pepi.2008.06.028, 2008. a, b
Sevilla, R. and Duretz, T.: A face-centred finite volume method for high-contrast Stokes interface problems, Int. J. Numer. Meth. Eng., 124, 3709–3732, https://doi.org/10.1002/nme.7294, 2023. a, b
Sevilla, R. and Duretz, T.: Face-centred finite volume methods for Stokes flows with variable viscosity, Int. J. Numer. Meth. Eng., 125, e7450, https://doi.org/10.1002/nme.7450, 2024. a
Sevilla, R., Giacomini, M., and Huerta, A.: A face-centred finite volume method for second-order elliptic problems, Int. J. Numer. Meth. Eng., 115, 986–1014, https://doi.org/10.1002/nme.5833, 2018. a, b, c
Shih, Y.-h., Stadler, G., and Wechsung, F.: Robust Multigrid Techniques for Augmented Lagrangian Preconditioning of Incompressible Stokes Equations with Extreme Viscosity Variations, SIAM J. Sci. Comput., https://doi.org/10.1137/21M1430698, 2022. a, b
Shin, D. and Strikwerda, J. C.: Inf-sup conditions for finite-difference approximations of the stokes equations, ANZIAM J., 39, 121–134, https://doi.org/10.1017/S0334270000009255, 1997. a
Spang, A., Thielmann, M., de Montserrat, A., and Duretz, T.: Transient Propagation of Ductile Ruptures by Thermal Runaway, J. Geophys. Res.-Sol. Ea., 130, e2025JB031240, https://doi.org/10.1029/2025JB031240, 2025. a
Tackley, P. J.: Effects of strongly variable viscosity on three-dimensional compressible convection in planetary mantles, J. Geophys. Res.-Sol. Ea., 101, 3311–3332, https://doi.org/10.1029/95JB03211, 1996. a, b
Tackley, P. J.: Modelling compressible mantle convection with large viscosity contrasts in a three-dimensional spherical shell using the yin-yang grid, Phys. Earth Planet. In., 171, 7–18, https://doi.org/10.1016/j.pepi.2008.08.005, 2008. a
Tackley, P. J.: Tectono-Convective Modes on Earth and Other Terrestrial Bodies, in: Dynamics of Plate Tectonics and Mantle Convection, 159–180, Elsevier, Walthm, MA, USA, ISBN 978-0-323-85733-8, https://doi.org/10.1016/B978-0-323-85733-8.00006-8, 2023. a
Vieira, L. M., Giacomini, M., Sevilla, R., and Huerta, A.: A second-order face-centred finite volume method for elliptic problems, Comput. Method. Appl. M., 358, 112655, https://doi.org/10.1016/j.cma.2019.112655, 2020. a
Yarushina, V. M. and Podladchikov, Y. Y.: (De)compaction of porous viscoelastoplastic media: Model formulation, J. Geophys. Res.-Sol. Ea., 120, 4146–4170, https://doi.org/10.1002/2014JB011258, 2015. a
Zhong, S., McNamara, A., Tan, E., Moresi, L., and Gurnis, M.: A benchmark study on mantle convection in a 3-D spherical shell using CitcomS, Geochem. Geophy. Geosy., 9, https://doi.org/10.1029/2008GC002048, 2008. a
Short summary
Geodynamic modeling helps scientists understand how the Earth deforms. New computer methods make these simulations faster and more efficient, especially on powerful computers. They automatically adjust settings for better performance and can handle complex materials and flow types. This approach makes it easier to study large, detailed models of Earth processes.
Geodynamic modeling helps scientists understand how the Earth deforms. New computer methods make...
Special issue