Articles | Volume 15, issue 14
https://doi.org/10.5194/gmd-15-5757-2022
© Author(s) 2022. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/gmd-15-5757-2022
© Author(s) 2022. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Assessing the robustness and scalability of the accelerated pseudo-transient method
Laboratory of Hydraulics, Hydrology and Glaciology (VAW), ETH Zurich, Zurich, Switzerland
Swiss Federal Institute for Forest, Snow and Landscape Research (WSL), Birmensdorf, Switzerland
Ivan Utkin
Laboratory of Hydraulics, Hydrology and Glaciology (VAW), ETH Zurich, Zurich, Switzerland
Swiss Federal Institute for Forest, Snow and Landscape Research (WSL), Birmensdorf, Switzerland
Faculty of Mechanics and Mathematics, Lomonosov Moscow State University, Moscow, Russia
Thibault Duretz
Institut für Geowissenschaften, Goethe‐Universität Frankfurt, Frankfurt, Germany
Univ. Rennes, CNRS, Géosciences Rennes UMR 6118, 35000 Rennes, France
Samuel Omlin
Swiss National Supercomputing Centre (CSCS), ETH Zurich, Lugano, Switzerland
Yuri Y. Podladchikov
Faculty of Mechanics and Mathematics, Lomonosov Moscow State University, Moscow, Russia
Institute of Earth Sciences, University of Lausanne, Lausanne, Switzerland
Swiss Geocomputing Centre, University of Lausanne, Lausanne, Switzerland
Related authors
Arne Spang, Marcel Thielmann, Casper Pranger, Albert de Montserrat, and Ludovic Räss
EGUsphere, https://doi.org/10.5194/egusphere-2025-2417, https://doi.org/10.5194/egusphere-2025-2417, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
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. 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.
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.
Arne Spang, Marcel Thielmann, Casper Pranger, Albert de Montserrat, and Ludovic Räss
EGUsphere, https://doi.org/10.5194/egusphere-2025-2417, https://doi.org/10.5194/egusphere-2025-2417, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
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. 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.
Thomas Geffroy, Philippe Yamato, Philippe Steer, Benjamin Guillaume, and Thibault Duretz
EGUsphere, https://doi.org/10.5194/egusphere-2025-1962, https://doi.org/10.5194/egusphere-2025-1962, 2025
Short summary
Short summary
While erosion's role in mountain building is well known, deformation from valley incision in inactive regions is less understood. Using our numerical models, we show that incision alone can cause significant crustal deformation and drive lower crust exhumation. This is favored in areas with thick crust, weak lower crust, and high plateaux. Our results show surface processes can reshape Earth's surface over time.
Yury Alkhimenkov and Yury Y. Podladchikov
Geosci. Model Dev., 18, 563–583, https://doi.org/10.5194/gmd-18-563-2025, https://doi.org/10.5194/gmd-18-563-2025, 2025
Short summary
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.
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.
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.
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.
Lorenzo G. Candioti, Thibault Duretz, Evangelos Moulas, and Stefan M. Schmalholz
Solid Earth, 12, 1749–1775, https://doi.org/10.5194/se-12-1749-2021, https://doi.org/10.5194/se-12-1749-2021, 2021
Short summary
Short summary
We quantify the relative importance of forces driving the dynamics of mountain building using two-dimensional computer simulations of long-term coupled lithosphere–upper-mantle deformation. Buoyancy forces can be as high as shear forces induced by far-field plate motion and should be considered when studying the formation of mountain ranges. The strength of rocks flooring the oceans and the density structure of the crust control deep rock cycling and the topographic elevation of orogens.
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.
Lorenzo G. Candioti, Stefan M. Schmalholz, and Thibault Duretz
Solid Earth, 11, 2327–2357, https://doi.org/10.5194/se-11-2327-2020, https://doi.org/10.5194/se-11-2327-2020, 2020
Short summary
Short summary
With computer simulations, we study the interplay between thermo-mechanical processes in the lithosphere and the underlying upper mantle during a long-term (> 100 Myr) tectonic cycle of extension–cooling–convergence. The intensity of mantle convection is important for (i) subduction initiation, (ii) the development of single- or double-slab subduction zones, and (iii) the forces necessary to initiate subduction. Our models are applicable to the opening and closure of the western Alpine Tethys.
Andrey Afanasyev and Ivan Utkin
Adv. Geosci., 54, 89–98, https://doi.org/10.5194/adgeo-54-89-2020, https://doi.org/10.5194/adgeo-54-89-2020, 2020
Short summary
Short summary
Understanding and predicting of subsurface geophysical and geotechnical processes, such as hydrothermal flows and oil and gas production, is limited by the fact that such processes cannot be directly observed. Limited monitoring can often be done only through surface observations including gravimetry and ground uplift or subsidence. The article describes the software that allows to constrain subsurface parameters by the observations at surface.
Cited articles
Alamatian, J.: A new formulation for fictitious mass of the Dynamic Relaxation
method with kinetic damping, Comput. Struct., 90–91, 42–54,
https://doi.org/10.1016/j.compstruc.2011.10.010, 2012. a
Alkhimenkov, Y., Khakimova, L., and Podladchikov, Y.: Stability of discrete
schemes of Biot's poroelastic equations,
Geophys. J. Int.,
225, 354–377, https://doi.org/10.1093/gji/ggaa584, 2021a. a, b
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, 7, https://doi.org/10.1029/2020JB021175, 2021b. a
Bakhvalov, N. S.: On the convergence of a relaxation method with natural
constraints on the elliptic operator,
USSR Comp. Math. Math.+, 6, 101–135, https://doi.org/10.1016/0041-5553(66)90118-2, 1966. a
Barnes, M. R.: Form Finding and Analysis of Tension Structures by Dynamic
Relaxation, International Journal of Space Structures, 14, 89–104,
https://doi.org/10.1260/0266351991494722, 1999. 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
Besard, T., Foket, C., and De Sutter, B.: Effective Extensible Programming:
Unleashing Julia on GPUs,
IEEE T. Parall. Distr., 30, 827–841, https://doi.org/10.1109/TPDS.2018.2872064, 2018. a
Besard, T., Churavy, V., Edelman, A., and De Sutter, B.: Rapid software
prototyping for heterogeneous and distributed platforms,
Adv. Eng. Softw., 132, 29–46, 2019. 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, b
Brandt, A.: Multi-level adaptive solutions to boundary-value problems,
Math. Comput., 31, 333–390,
https://doi.org/10.1090/S0025-5718-1977-0431719-X, 1977. a, b
Byrne, S., Wilcox, L. C., and Churavy, V.: MPI.jl: Julia bindings for the
Message Passing Interface, Proceedings of the JuliaCon Conferences, 1, 68,
https://doi.org/10.21105/jcon.00068, 2021. a
Cassell, A. C. and Hobbs, R. E.: Numerical stability of dynamic relaxation
analysis of non-linear structures,
Int. J. Numer. Meth. Eng., 10, 1407–1410, https://doi.org/10.1002/nme.1620100620, 1976. a
Chester, M.: Second sound in solids, Phys. Rev., 131, 2013–2015,
https://doi.org/10.1103/PhysRev.131.2013, 1963. a
Costa, A.: Permeability-porosity relationship: A reexamination of the
Kozeny-Carman equation based on a fractal pore-space geometry assumption,
Geophys. Res. Lett., 33, L02318, https://doi.org/10.1029/2005GL025134, 2006. a
Cox, S. and Zuazua, E.: The rate at which energy decays in a damped string,
Commun. Part. Diff. Eq., 19, 213–243,
https://doi.org/10.1080/03605309408821015, 1994. a
Cundall, A. P.: Explicit finite differnce method in geomechanics, Second Int.
Conf. Numerical Methods in Geomechanics, Blacksburg, 1976, 1, 132–150, 1976. a
Eager, D. L., Zahorjan, J., and Lazowska, E. D.: Speedup versus efficiency in
parallel systems, IEEE T. Comput., 38, 408–423,
https://doi.org/10.1109/12.21127, 1989. a
Fuhrer, O., Chadha, T., Hoefler, T., Kwasniewski, G., Lapillonne, X., Leutwyler, D., Lüthi, D., Osuna, C., Schär, C., Schulthess, T. C., and Vogt, H.: Near-global climate simulation at 1 km resolution: establishing a performance baseline on 4888 GPUs with COSMO 5.0, Geosci. Model Dev., 11, 1665–1681, https://doi.org/10.5194/gmd-11-1665-2018, 2018. a
Gaitonde, A. L.: A dual-time method for two-dimensional unsteady incompressible
flow calculations, Int. J. Numer. Meth. Eng., 41, 1153–1166,
https://doi.org/10.1002/(SICI)1097-0207(19980330)41:6<1153::AID-NME334>3.0.CO;2-9,
1998. a
Gustafson, J. L.: Reevaluating Amdahl's law, Commun. ACM, 31,
532–533, https://doi.org/10.1145/42411.42415, 1988. a
Hackbusch, W.: Multi-Grid Methods and Applications, vol. 4, Springer
Series in Computational Mathematics, Springer, Berlin,
Heidelberg, https://doi.org/10.1007/978-3-662-02427-0, 1985. a
Hoemmen, M. F.: Communication-avoiding Krylov subspace methods, PhD thesis,
EECS Department, University of California, Berkeley,
http://www2.eecs.berkeley.edu/Pubs/TechRpts/2010/EECS-2010-37.html (last access: 16 May 2022),
2010. a
Jordan, P. M. and Puri, A.: Digital signal propagation in dispersive media,
J. Appl. Phys., 85, 1273–1282, https://doi.org/10.1063/1.369258, 1999. a
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, Jülich, Germany, 11–12 February 2016,
edited by: Binder, K., Müller, M., Kremer, A., and Schnurpfeil, A.,
48, 299–307, Forschungszentrum Jülich, Jülich, 2016. a
Kilic, B. and Madenci, E.: Structural stability and failure analysis using
peridynamic theory, Int. J. NonLin. Mech., 44,
845–854, https://doi.org/10.1016/j.ijnonlinmec.2009.05.007, 2009. a
Kumar, V., Grama, A., Gupta, A., and Karypis, G.: Introduction to parallel
computing, Benjamin/Cummings Redwood City, CA, 110, ISBN 0805331700, 1994. a
Mandal, J., Sonawane, C., Iyer, A., and GosaviInamdar, S.: Incompressible flow
computations over moving boundary using a novel upwind method, Comput.
Fluid., 46, 348–352, 2011. a
Maxwell, J. C.: IV. On the dynamical theory of gases,
Philos. T. R. Soc. Lond., 157, 49–88,
https://doi.org/10.1098/rstl.1867.0004, 1867. 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
McCalpin, J. D.: Memory bandwidth and machine balance in current high
performance computers, IEEE computer society technical committee on computer
architecture (TCCA) newsletter, 2, 19–25, 1995. a
McKee, S., Tomé, M., Ferreira, V., Cuminato, J., Castelo, A., Sousa, F.,
and Mangiavacchi, N.: The MAC method, Comput. Fluid., 37, 907–930,
https://doi.org/10.1016/j.compfluid.2007.10.006, 2008. a
Omlin, S.: Development of massively parallel near peak performance solvers for
three-dimensional geodynamic modelling, PhD thesis, University of Lausanne, https://doi.org/10.48550/ARXIV.2207.08716,
2017. a
Omlin, S. and Räss, L.: ImplicitGlobalGrid.jl,
https://github.com/eth-cscs/ImplicitGlobalGrid.jl (last access: 16 May 2022),
2021a. a
Omlin, S., Räss, L., and Podladchikov, Y.: HPC.m – the MATLAB HPC
Compiler and its Use for Solving 3D Poromechanics on Supercomputers, in:
Platform for Advanced Scientific Computing Conference (PASC15), 1–3 June 2015, Zurich, Switzerland, ETH Zurich,
Zurich, Switzerland, 2015a. a
Omlin, S., Räss, L., and Podladchikov, Y.: From MATLAB to
Supercomputing: the MATLAB HPC-Compiler, in: 2nd Frontiers in Computational
Physics: Energy Sciences Conference 2015, 3–5 June 2015, Zurich, Switzerland, ETH Zurich, Zurich, Switzerland,
2015b. a
Otter, J.: Computations for prestressed concrete reactor pressure vessels
using dynamic relaxation, Nuclear Structural Engineering, 1, 61–75,
https://doi.org/10.1016/0369-5816(65)90097-9, 1965. a
Otter, J. R. H., Cassell, A. C., and Hobbs, R. E.: Dynamic relaxation,
35,
633–656,
https://doi.org/10.1680/IICEP.1966.8604, 1966. a
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
Pascal, H.: Pressure wave propagation in a fluid flowing through a porous
medium and problems related to interpretation of Stoneley's wave attenuation
in acoustical well logging, Int. J. Eng. Sci., 24,
1553–1570, https://doi.org/10.1016/0020-7225(86)90163-1, 1986. a
Patankar, S.: Numerical heat transfer and fluid flow, CRC press, https://doi.org/10.1201/9781482234213, 1980. a
Podladtchikov, I. and Podladchikov, Y. Y.: Memory Bound Wave Propagation at
Hardware Limit,
http://on-demand.gputechconf.com/gtc/2013/video/S3176-Memory-Bound-Wave-Propagation-Hardware-Limit.flv (last access: 16 May 2022),
2013. a
Poliakov, A. N., 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, Flow and creep in the solar system: observations, modeling and
theory, NATO ASI Series, vol. 391, Springer, Dordrecht, 175–195, https://doi.org/10.1007/978-94-015-8206-3_12, 1993. a
Poliakov, A. N. B., Herrmann, H. J., Podladchikov, Y. Y., and Roux, S.:
Fractal plastic shear bands, Fractals, 2, 567–581,
https://doi.org/10.1142/s0218348x9400079x, 1994. a
Ramesh, G. and Krishnamoorthy, C. S.: Post‐buckling analysis of structures
by dynamic relaxation, Int. J. Numer. Meth. Eng., 36, 1339–1364, https://doi.org/10.1002/nme.1620360806, 1993. a
Räss, L. and Omlin, S.: ParallelRandomFields.jl,
https://github.com/luraess/ParallelRandomFields.jl (last access: 16 May 2022), 2021. a
Räss, L. and Utkin, I.: PTsolvers/PseudoTransientDiffusion.jl:
PseudoTransientDiffusion.jl 1.0.0, Zenodo [code], https://doi.org/10.5281/zenodo.6553699,
2022a. a, b
Räss, L. and Utkin, I.: PTsolvers/PseudoTransientStokes.jl:
PseudoTransientStokes.jl 1.0.0, Zenodo [code], https://doi.org/10.5281/zenodo.6553714,
2022b. a, b
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
Räss, L., Kolyukhin, D., and Minakov, A.: Efficient parallel random field
generator for large 3-D geophysical problems, Comput. Geosci.,
131, 158–169, https://doi.org/10.1016/j.cageo.2019.06.007, 2019b. a
Räss, L., Omlin, S., and Podladchikov, Y. Y.: Resolving Spontaneous
Nonlinear Multi-Physics Flow Localization in 3-D: Tackling Hardware Limit,
https://developer.nvidia.com/gtc/2019/video/S9368 (last access: 16 May 2022), 2019c. 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, b, c, d, e, f
Reuber, G. S., Holbach, L., and Räss, L.: Adjoint-based inversion for
porosity in shallow reservoirs using pseudo-transient solvers for non-linear
hydro-mechanical processes, J. Comput. Phys., 423, 109797,
https://doi.org/10.1016/j.jcp.2020.109797, 2020. a
Rezaiee-Pajand, M., Kadkhodayan, M., Alamatian, J., and Zhang, L. C.: A new
method of fictitious viscous damping determination for the dynamic relaxation
method, Comput. Struct., 89, 783–794,
https://doi.org/10.1016/J.COMPSTRUC.2011.02.002, 2011. a, b
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, https://doi.org/10.1098/rsta.1911.0009, 1911. a
Riley, J. D.: Iteration Procedures for the Dirichlet Difference
Problem, Mathematical Tables and Other Aids to Computation, 8, 125–131,
https://doi.org/10.2307/2001924, 1954. a, b
Saad, Y.: Iterative methods for linear systems of equations: A brief
historical journey, in: Contemporary Mathematics, edited by: Brenner, S.,
Shparlinski, I., Shu, C.-W., and Szyld, D., American
Mathematical Society, Providence, Rhode Island, 754, 197–215, https://doi.org/10.1090/conm/754/15141,
2020. a, b
Shin, D. and Strikwerda, J. C.: Inf-sup conditions for finite-difference
approximations of the Stokes equations, J. Aust.
Math. Soc. B, 39, 121–134, 1997. a
Swirydowicz, K., Langou, J., Ananthan, S., Yang, U., and Thomas, S.: Low
synchronization Gram-Schmidt and generalized minimal residual algorithms,
Numer. Linear Algebr., 28, e2343, https://doi.org/10.1002/nla.2343, 2020.
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
Trottenberg, U., Oosterlee, C. W. C. W., and Schüller, A.: Multigrid,
Academic Press, hardcover ISBN 9780127010700,
eBook ISBN 9780080479569, 2001. a
Virieux, J.: P-SV wave propagation in heterogeneous media: Velocity‐stress
finite‐difference method, Geophysics, 51, 889–901,
https://doi.org/10.1190/1.1442147, 1986. a
Wang, L. H., Yarushina, V., 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, 2021. a, b, c
Wulf, W. A. and McKee, S. A.: Hitting the memory wall: implications of the
obvious, ACM Comp. Ar., 23, 20–24,
https://doi.org/10.1145/216585.216588, 1995. a
Young, D. M.: Second-degree iterative methods for the solution of large linear
systems, J. Approx. Theory, 5, 137–148,
https://doi.org/10.1016/0021-9045(72)90036-6, 1972. a
Zheng, L., Gerya, T., Knepley, M., Yuen, D. A., Zhang, H., and Shi, Y.: GPU
implementation of multigrid solver for stokes equation with strongly variable
viscosity, in: Lecture Notes in Earth System Sciences, 0, 321–333,
Springer, Berlin, Heidelberg, https://doi.org/10.1007/978-3-642-16405-7_21, 2013. a
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.
Continuum mechanics-based modelling of physical processes at large scale requires huge...