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.…
Verfasserangaben: | Seoleun Shin, Gert ZöllerORCiDGND, Matthias HolschneiderORCiDGND, Sebastian ReichORCiDGND |
---|---|
DOI: | https://doi.org/10.1016/j.cageo.2010.11.011 |
ISSN: | 0098-3004 |
Titel des übergeordneten Werks (Englisch): | 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 |
Verlag: | Elsevier |
Verlagsort: | Oxford |
Publikationstyp: | Wissenschaftlicher Artikel |
Sprache: | Englisch |
Jahr der Erstveröffentlichung: | 2011 |
Erscheinungsjahr: | 2011 |
Datum der Freischaltung: | 26.03.2017 |
Freies Schlagwort / Tag: | Multigrid; Multiple time stepping; Strike-slip fault model |
Band: | 37 |
Ausgabe: | 8 |
Seitenanzahl: | 8 |
Erste Seite: | 1075 |
Letzte Seite: | 1082 |
Fördernde Institution: | German Research Society (DFG), Collaborative research centre [SFB555]; DFG [HA4789/2-1] |
Organisationseinheiten: | Mathematisch-Naturwissenschaftliche Fakultät / Institut für Mathematik |
Peer Review: | Referiert |