Preprints
https://doi.org/10.5194/gmd-2024-160
https://doi.org/10.5194/gmd-2024-160
Submitted as: development and technical paper
 | 
10 Sep 2024
Submitted as: development and technical paper |  | 10 Sep 2024
Status: this preprint is currently under review for the journal GMD.

Accelerated pseudo-transient method for elastic, viscoelastic, and coupled hydro-mechanical problems with applications

Yury Alkhimenkov and Yury Y. Podladchikov

Abstract. The Accelerated Pseudo-Transient (APT) method is a matrix-free approach used to solve partial differential equations (PDEs), characterized by its reliance on local operations, which makes it highly suitable for parallelization. With the advent of the memory-wall phenomenon around 2005, where memory access speed overtook floating-point operations as the bottleneck in high-performance computing, the APT method has gained prominence as a powerful tool for tackling various PDEs in geosciences. Recent advancements have demonstrated the APT method's computational efficiency, particularly when applied to quasi-static nonlinear problems using Graphical Processing Units (GPUs). This manuscript presents a comprehensive analysis of the APT method, focusing on its application to quasi-static elastic, viscoelastic, and coupled hydro-mechanical problems, specifically those governed by quasi-static Biot's poroelastic equations, across 1D, 2D, and 3D domains. We systematically investigate the optimal numerical parameters required to achieve rapid convergence, offering valuable insights into the method's applicability and efficiency for a range of physical models. Our findings are validated against analytical solutions, underscoring the robustness and accuracy of the APT method in both homogeneous and heterogeneous media. We explore the influence of boundary conditions, non-linearities, and coupling on the optimal convergence parameters, highlighting the method's adaptability in addressing complex and realistic scenarios. To demonstrate the flexibility of the APT method, we apply it to the nonlinear mechanical problem of strain localization using a poro-elasto-viscoplastic rheological model, achieving extremely high resolutions – 10,0002 voxels in 2D and 5123 voxels in 3D – that, to our knowledge, have not been previously explored for such models. Our study contributes significantly to the field by providing a robust framework for the effective implementation of the APT method in solving challenging geophysical problems. Importantly, the results presented in this paper are fully reproducible, with Matlab, symbolic Maple scripts, and CUDA C codes made available in a permanent repository.

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 preprint. The responsibility to include appropriate place names lies with the authors.
Yury Alkhimenkov and Yury Y. Podladchikov

Status: open (until 05 Nov 2024)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
Yury Alkhimenkov and Yury Y. Podladchikov
Yury Alkhimenkov and Yury Y. Podladchikov

Viewed

Total article views: 98 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
66 9 23 98 1 1
  • HTML: 66
  • PDF: 9
  • XML: 23
  • Total: 98
  • BibTeX: 1
  • EndNote: 1
Views and downloads (calculated since 10 Sep 2024)
Cumulative views and downloads (calculated since 10 Sep 2024)

Viewed (geographical distribution)

Total article views: 89 (including HTML, PDF, and XML) Thereof 89 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
1
 
 
 
 
Latest update: 04 Oct 2024
Download
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 hydro-mechanical 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 can handle complex problems in high-resolution models.