A multigrid solver for modeling complex interseismic stress fields

详细信息    查看全文

• 作者：Seoleun Shin ; ^a ; ^{seshin@uni-potsdam.de"" rel=""nofollow} ; Gert Zö ; ller^a ; Matthias Holschneider^a ; Sebastian Reich^a
• 关键词：Multigrid ; Multiple time stepping ; Strike-slip fault model
• 刊名：Computers and Geosciences
• 出版年：2011
• 出版时间：August 2011
• 年：2011
• 卷：37
• 期：8
• 页码：1075-1082
• 全文大小：728 K

文摘

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 the entire half-space.