An FC-based spectral solver for elastodynamic problems in general three-dimensional domains

详细信息    查看全文

• 作者：Faisal Amlani ; Oscar P. Bruno ; ^{obruno@caltech.edu" class="auth_mail" title="E-mail the corresponding author}
• 关键词：Elastic wave equation ; Spectral solver ; General domain ; FC solver ; Seismic simulations ; Non-destructive testing
• 刊名：Journal of Computational Physics
• 出版年：2016
• 出版时间：15 February 2016
• 年：2016
• 卷：307
• 期：Complete
• 页码：333-354
• 全文大小：2357 K

文摘

This paper presents a spectral numerical algorithm for the solution of elastodynamics problems in general three-dimensional domains. Based on a recently introduced “Fourier continuation” (FC) methodology for accurate Fourier expansion of non-periodic functions, the proposed approach possesses a number of appealing properties: it yields results that are essentially free of dispersion errors, it entails mild CFL constraints, it runs at a cost that scales linearly with the discretization sizes, and it lends itself easily to efficient parallelization in distributed-memory computing clusters. The proposed algorithm is demonstrated in this paper by means of a number of applications to problems of isotropic elastodynamics that arise in the fields of materials science and seismology. These examples suggest that the new approach can yield solutions within a prescribed error tolerance by means of significantly smaller discretizations and shorter computing times than those required by other methods.