A multigrid solver for modeling complex interseismic stress fields
- We develop a multigrid, multiple time stepping scheme to reduce computational efforts for calculating complex stress interactions in a strike-slip 2D planar fault for the simulation of seismicity. The key elements of the multilevel solver are separation of length scale, grid-coarsening, and hierarchy. In this study the complex stress interactions are split into two parts: the first with a small contribution is computed on a coarse level, and the rest for strong interactions is on a fine level. This partition leads to a significant reduction of the number of computations. The reduction of complexity is even enhanced by combining the multigrid with multiple time stepping. Computational efficiency is enhanced by a factor of 10 while retaining a reasonable accuracy, compared to the original full matrix-vortex multiplication. The accuracy of solution and computational efficiency depend on a given cut-off radius that splits multiplications into the two parts. The multigrid scheme is constructed in such a way that it conserves stress in theWe develop a multigrid, multiple time stepping scheme to reduce computational efforts for calculating complex stress interactions in a strike-slip 2D planar fault for the simulation of seismicity. The key elements of the multilevel solver are separation of length scale, grid-coarsening, and hierarchy. In this study the complex stress interactions are split into two parts: the first with a small contribution is computed on a coarse level, and the rest for strong interactions is on a fine level. This partition leads to a significant reduction of the number of computations. The reduction of complexity is even enhanced by combining the multigrid with multiple time stepping. Computational efficiency is enhanced by a factor of 10 while retaining a reasonable accuracy, compared to the original full matrix-vortex multiplication. The accuracy of solution and computational efficiency depend on a given cut-off radius that splits multiplications into the two parts. The multigrid scheme is constructed in such a way that it conserves stress in the entire half-space.…
Author details: | Seoleun Shin, Gert ZöllerORCiDGND, Matthias HolschneiderORCiDGND, Sebastian ReichORCiDGND |
---|---|
DOI: | https://doi.org/10.1016/j.cageo.2010.11.011 |
ISSN: | 0098-3004 |
Title of parent work (English): | Computers & geosciences : an international journal devoted to the publication of papers on all aspects of geocomputation and to the distribution of computer programs and test data sets ; an official journal of the International Association for Mathematical Geology |
Publisher: | Elsevier |
Place of publishing: | Oxford |
Publication type: | Article |
Language: | English |
Year of first publication: | 2011 |
Publication year: | 2011 |
Release date: | 2017/03/26 |
Tag: | Multigrid; Multiple time stepping; Strike-slip fault model |
Volume: | 37 |
Issue: | 8 |
Number of pages: | 8 |
First page: | 1075 |
Last Page: | 1082 |
Funding institution: | German Research Society (DFG), Collaborative research centre [SFB555]; DFG [HA4789/2-1] |
Organizational units: | Mathematisch-Naturwissenschaftliche Fakultät / Institut für Mathematik |
Peer review: | Referiert |