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 |
