Accelerated pseudo-transient method for elastic, viscoelastic, and coupled hydromechanical problems with applications
<p>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, wher...
Saved in:
| Main Authors: | , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Copernicus Publications
2025-01-01
|
| Series: | Geoscientific Model Development |
| Online Access: | https://gmd.copernicus.org/articles/18/563/2025/gmd-18-563-2025.pdf |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1832575982724186112 |
|---|---|
| author | Y. Alkhimenkov Y. Y. Podladchikov |
| author_facet | Y. Alkhimenkov Y. Y. Podladchikov |
| author_sort | Y. Alkhimenkov |
| collection | DOAJ |
| description | <p>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 study presents a comprehensive analysis of the APT method, focusing on its application to quasi-static elastic, viscoelastic, and coupled hydromechanical problems, specifically those governed by quasi-static Biot 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, nonlinearities, 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 – <span class="inline-formula">10 000<sup>2</sup></span> in 2D and <span class="inline-formula">512<sup>3</sup></span> 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 code, symbolic Maple scripts, and CUDA C codes made available in a permanent repository.</p> |
| format | Article |
| id | doaj-art-2b72c6f1bfda49fe8be09610b4c0fd17 |
| institution | Kabale University |
| issn | 1991-959X 1991-9603 |
| language | English |
| publishDate | 2025-01-01 |
| publisher | Copernicus Publications |
| record_format | Article |
| series | Geoscientific Model Development |
| spelling | doaj-art-2b72c6f1bfda49fe8be09610b4c0fd172025-01-31T14:01:22ZengCopernicus PublicationsGeoscientific Model Development1991-959X1991-96032025-01-011856358310.5194/gmd-18-563-2025Accelerated pseudo-transient method for elastic, viscoelastic, and coupled hydromechanical problems with applicationsY. Alkhimenkov0Y. Y. Podladchikov1Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USAInstitute of Earth Sciences, University of Lausanne, Lausanne, Switzerland<p>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 study presents a comprehensive analysis of the APT method, focusing on its application to quasi-static elastic, viscoelastic, and coupled hydromechanical problems, specifically those governed by quasi-static Biot 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, nonlinearities, 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 – <span class="inline-formula">10 000<sup>2</sup></span> in 2D and <span class="inline-formula">512<sup>3</sup></span> 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 code, symbolic Maple scripts, and CUDA C codes made available in a permanent repository.</p>https://gmd.copernicus.org/articles/18/563/2025/gmd-18-563-2025.pdf |
| spellingShingle | Y. Alkhimenkov Y. Y. Podladchikov Accelerated pseudo-transient method for elastic, viscoelastic, and coupled hydromechanical problems with applications Geoscientific Model Development |
| title | Accelerated pseudo-transient method for elastic, viscoelastic, and coupled hydromechanical problems with applications |
| title_full | Accelerated pseudo-transient method for elastic, viscoelastic, and coupled hydromechanical problems with applications |
| title_fullStr | Accelerated pseudo-transient method for elastic, viscoelastic, and coupled hydromechanical problems with applications |
| title_full_unstemmed | Accelerated pseudo-transient method for elastic, viscoelastic, and coupled hydromechanical problems with applications |
| title_short | Accelerated pseudo-transient method for elastic, viscoelastic, and coupled hydromechanical problems with applications |
| title_sort | accelerated pseudo transient method for elastic viscoelastic and coupled hydromechanical problems with applications |
| url | https://gmd.copernicus.org/articles/18/563/2025/gmd-18-563-2025.pdf |
| work_keys_str_mv | AT yalkhimenkov acceleratedpseudotransientmethodforelasticviscoelasticandcoupledhydromechanicalproblemswithapplications AT yypodladchikov acceleratedpseudotransientmethodforelasticviscoelasticandcoupledhydromechanicalproblemswithapplications |