Articles | Volume 18, issue 2
https://doi.org/10.5194/gmd-18-563-2025
© Author(s) 2025. 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-18-563-2025
© Author(s) 2025. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Accelerated pseudo-transient method for elastic, viscoelastic, and coupled hydromechanical problems with applications
Yury Alkhimenkov
CORRESPONDING AUTHOR
Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
Yury Y. Podladchikov
Institute of Earth Sciences, University of Lausanne, Lausanne, Switzerland
Related authors
Yury Alkhimenkov, Lyudmila Khakimova, and Yury Podladchikov
EGUsphere, https://doi.org/10.5194/egusphere-2024-3237, https://doi.org/10.5194/egusphere-2024-3237, 2024
Short summary
Short summary
This study examines stress drops and earthquake nucleation in elasto-plastic media using 2D simulations, highlighting the importance of high temporal and spatial resolutions in capturing stress evolution and strain fields. Stress drops reflect fault rupture mechanics and emulate earthquake behavior. The non-Gaussian distribution of stress drop amplitudes resembles "solid turbulence." Elasto-plastic models simulate key earthquake processes and could improve seismic hazard assessment.
Yury Alkhimenkov and Yury Y. Podladchikov
EGUsphere, https://doi.org/10.5194/egusphere-2024-3238, https://doi.org/10.5194/egusphere-2024-3238, 2024
Short summary
Short summary
This paper presents a rigorous derivation of Gassmann's equations, grounded in thermodynamic principles and conservation laws, addressing gaps and potential inconsistencies in the original formulation. It also explores Biot's poroelastic equations, demonstrating that Gassmann's equations are a specific case within Biot’s framework. The study affirms the robustness of Gassmann's equations when assumptions are met, and symbolic Maple routines are provided to ensure reproducibility of the results.
Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
Geosci. Model Dev., 14, 7749–7774, https://doi.org/10.5194/gmd-14-7749-2021, https://doi.org/10.5194/gmd-14-7749-2021, 2021
Short summary
Short summary
We propose an implementation of the material point method using graphical processing units (GPUs) to solve elastoplastic problems in three-dimensional configurations, such as the granular collapse or the slumping mechanics, i.e., landslide. The computational power of GPUs promotes fast code executions, compared to a traditional implementation using central processing units (CPUs). This allows us to study complex three-dimensional problems tackling high spatial resolution.
Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
Geosci. Model Dev., 13, 6265–6284, https://doi.org/10.5194/gmd-13-6265-2020, https://doi.org/10.5194/gmd-13-6265-2020, 2020
Short summary
Short summary
In this work, we present an efficient and fast material point method (MPM) implementation in MATLAB. We first discuss the vectorization strategies to adapt this numerical method to a MATLAB implementation. We report excellent agreement of the solver compared with classical analysis among the MPM community, such as the cantilever beam problem. The solver achieves a performance gain of 28 compared with a classical iterative implementation.
Yury Alkhimenkov, Lyudmila Khakimova, and Yury Podladchikov
EGUsphere, https://doi.org/10.5194/egusphere-2024-3237, https://doi.org/10.5194/egusphere-2024-3237, 2024
Short summary
Short summary
This study examines stress drops and earthquake nucleation in elasto-plastic media using 2D simulations, highlighting the importance of high temporal and spatial resolutions in capturing stress evolution and strain fields. Stress drops reflect fault rupture mechanics and emulate earthquake behavior. The non-Gaussian distribution of stress drop amplitudes resembles "solid turbulence." Elasto-plastic models simulate key earthquake processes and could improve seismic hazard assessment.
Yury Alkhimenkov and Yury Y. Podladchikov
EGUsphere, https://doi.org/10.5194/egusphere-2024-3238, https://doi.org/10.5194/egusphere-2024-3238, 2024
Short summary
Short summary
This paper presents a rigorous derivation of Gassmann's equations, grounded in thermodynamic principles and conservation laws, addressing gaps and potential inconsistencies in the original formulation. It also explores Biot's poroelastic equations, demonstrating that Gassmann's equations are a specific case within Biot’s framework. The study affirms the robustness of Gassmann's equations when assumptions are met, and symbolic Maple routines are provided to ensure reproducibility of the results.
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.
Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
Geosci. Model Dev., 14, 7749–7774, https://doi.org/10.5194/gmd-14-7749-2021, https://doi.org/10.5194/gmd-14-7749-2021, 2021
Short summary
Short summary
We propose an implementation of the material point method using graphical processing units (GPUs) to solve elastoplastic problems in three-dimensional configurations, such as the granular collapse or the slumping mechanics, i.e., landslide. The computational power of GPUs promotes fast code executions, compared to a traditional implementation using central processing units (CPUs). This allows us to study complex three-dimensional problems tackling high spatial resolution.
Martin Franz, Michel Jaboyedoff, Ryan P. Mulligan, Yury Podladchikov, and W. Andy Take
Nat. Hazards Earth Syst. Sci., 21, 1229–1245, https://doi.org/10.5194/nhess-21-1229-2021, https://doi.org/10.5194/nhess-21-1229-2021, 2021
Short summary
Short summary
A landslide-generated tsunami is a complex phenomenon that involves landslide dynamics, wave dynamics and their interaction. This phenomenon threatens numerous lives and infrastructures around the world. To assess this natural hazard, we developed an efficient numerical model able to simulate the landslide, the momentum transfer and the wave all at once. The good agreement between the numerical simulations and physical experiments validates our model and its novel momentum transfer approach.
Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
Geosci. Model Dev., 13, 6265–6284, https://doi.org/10.5194/gmd-13-6265-2020, https://doi.org/10.5194/gmd-13-6265-2020, 2020
Short summary
Short summary
In this work, we present an efficient and fast material point method (MPM) implementation in MATLAB. We first discuss the vectorization strategies to adapt this numerical method to a MATLAB implementation. We report excellent agreement of the solver compared with classical analysis among the MPM community, such as the cantilever beam problem. The solver achieves a performance gain of 28 compared with a classical iterative implementation.
Cited articles
Alkhimenkov, Y.: Numerical validation of Gassmann’s equations, Geophysics, 88, A25–A29, 2023. a
Alkhimenkov, Y. and Podladchikov, Y.: APTsolver, Zenodo [code and data set], https://doi.org/10.5281/zenodo.14056939, 2024. a
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, e2020JB021175, https://doi.org/10.1029/2020JB021175, 2021b. a, b, c
Alkhimenkov, Y., Khakimova, L., and Podladchikov, Y.: Shear bands triggered by solitary porosity waves in deforming fluid-saturated porous media, Geophys. Res. Lett., 51, e2024GL108789, https://doi.org/10.1029/2024GL108789, 2024a. a, b
Alkhimenkov, Y., Khakimova, L., Utkin, I., and Podladchikov, Y.: Resolving Strain Localization in Frictional and Time-Dependent Plasticity: Two- and Three-Dimensional Numerical Modeling Study Using Graphical Processing Units (GPUs), J. Geophys. Res.-Sol. Ea., 129, e2023JB028566, https://doi.org/10.1029/2023JB028566, 2024b. a, b, c, d
Benyamin, M., Calder, J., Sundaramoorthi, G., and Yezzi, A.: Accelerated variational PDEs for efficient solution of regularized inversion problems, J. Math. Imaging Vis., 62, 10–36, 2020. a
Biot, M. A.: Mechanics of deformation and acoustic propagation in porous media, J. Appl. Phys., 33, 1482–1498, 1962. a
Calder, J. and Yezzi, A.: PDE acceleration: a convergence rate analysis and applications to obstacle problems, Res. Math. Sci., 6, 35, https://doi.org/10.1007/s40687-019-0197-x, 2019. a
Cox, S. and Zuazua, E.: The rate at which energy decays in a damped string, Commun. Part. Diff. Eq., 19, 213–243, 1994. a
Cundall, P.: Explicit finite-difference methods in geomechanics, Proc. 2nd Int. Cof. Num. Meth. Geomech., ASCE, New York, 132–150, https://cir.nii.ac.jp/crid/1571417124088087168 (last access: 30 January 2025), 1976. a
Duretz, T., de Borst, R., and Le Pourhiet, L.: Finite thickness of shear bands in frictional viscoplasticity and implications for lithosphere dynamics, Geochem. Geophy. Geosy., 20, 5598–5616, 2019. a
Gassmann, F.: Uber die elastizitat poroser medien, Vierteljahrsschrift der Naturforschenden Gesellschaft in Zurich, 96, 1–23, 1951. a
Hirsch, C.: Numerical computation of internal and external flows, Volume 1: Fundamentals of numerical discretization, John Wiley and Sons, ISBN 0471917621, 1988. a
Landau, L. D. and Lifshitz, E. M.: Course of Theoretical Physics Vol 7: Theory of Elasticity, Pergamon press, 1959. a
Li, S. and Wang, G.: Introduction to micromechanics and nanomechanics, World Scientific Publishing Company, https://doi.org/10.1142/6834, 2008. a
Nemat-Nasser, S. and Hori, M.: Micromechanics: overall properties of heterogeneous materials, Elsevier, https://doi.org/10.1115/1.2788912, 2013. a, b
Omlin, S., Malvoisin, B., and Podladchikov, Y. Y.: Pore fluid extraction by reactive solitary waves in 3-D, Geophys. Res. Lett., 44, 9267–9275, 2017. a
Otter, J.: Computations for prestressed concrete reactor pressure vessels using dynamic relaxation, Nuclear Structural Engineering, 1, 61–75, 1965. a
Otter, J. R. H., Cassell, A. C., Hobbs, R. E., and POISSON: Dynamic relaxation, P. I. Civil Eng., 35, 633–656, 1966. a
Poliakov, A., Cundall, P., Podladchikov, Y., and Lyakhovsky, V.: 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, Springer, 175–195, Print ISBN 978-90-481-4245-3, https://doi.org/10.1007/978-94-015-8206-3_12, 1993a. a, b
Poliakov, A., Podladchikov, Y., and Talbot, C.: Initiation of salt diapirs with frictional overburdens: numerical experiments, Tectonophysics, 228, 199–210, 1993b. a
Poliakov, A., Herrmann, H. J., Podladchikov, Y. Y., and Roux, S.: Fractal plastic shear bands, Fractals, 2, 567–581, 1994. a
Polyak, B. T.: Some methods of speeding up the convergence of iteration methods, USSR Comp. Math. Math.+, 4, 1–17, https://doi.org/10.1016/0041-5553(64)90137-5, 1964. a
Räss, L., Duretz, T., and Podladchikov, Y.: Resolving hydromechanical coupling in two and three dimensions: spontaneous channelling of porous fluids owing to decompaction weakening, Geophys. J. Int., 218, 1591–1616, 2019. a
Räss, L., Licul, A., Herman, F., Podladchikov, Y. Y., and Suckale, J.: Modelling thermomechanical ice deformation using an implicit pseudo-transient method (FastICE v1.0) based on graphical processing units (GPUs), Geosci. Model Dev., 13, 955–976, https://doi.org/10.5194/gmd-13-955-2020, 2020. a
Richardson, L. F.: IX. The approximate arithmetical solution by finite differences of physical problems involving differential equations, with an application to the stresses in a masonry dam, Philos. T. R. Soc. A, 210, 307–357, 1911. a
Riley, J. D.: Iteration procedures for the Dirichlet difference problem, Mathematical Tables and Other Aids to Computation, 8, 125–131, 1954. a
Wang, L. H., Yarushina, V. M., Alkhimenkov, Y., and Podladchikov, Y.: Physics-inspired pseudo-transient method and its application in modelling focused fluid flow with geological complexity, Geophys. J. Int., 229, 1–20, https://doi.org/10.1093/gji/ggab426, 2022. a, b
Young, D. M.: Second-degree iterative methods for the solution of large linear systems, J. Approx. Theory, 5, 137–148, 1972. a
Short summary
The accelerated pseudo-transient (APT) method is an efficient way to solve partial differential equations, particularly well-suited for parallel computing. This paper explores the APT method's effectiveness in solving elastic, viscoelastic, and hydromechanical problems, focusing on quasi-static conditions in 1D, 2D, and 3D. The study examines the best numerical settings for fast and accurate solutions. The paper shows how the APT method can handle complex problems in high-resolution models.
The accelerated pseudo-transient (APT) method is an efficient way to solve partial differential...
Special issue