Poly-Grid Spectral Element Modeling for Wave Propagation in Complex Elastic Media

Abstract

Modeling elastic waves in complex media, with varying physical properties, require very accurate algorithms and a great computational effort to avoid nonphysical effects. Among the numerical methods the spectral elements (SEM) have a high precision and ease in modeling such problems and the physical domains can be discretized using very coarse meshes with elements of constant properties. In many cases, models with very complex geometries and small heterogeneities, shorter than the minimum wavelength, require grid resolution down to the thinnest scales, resulting in an extremely large problem size and greatly reducing accuracy and computational efficiency. In this paper, a poly-grid method (PG-CSEM) is presented that can overcome this limitation. To accurately deal with continuous variations or even small-scale fluctuations in elastic properties, temporary auxiliary grids are introduced that prevent the need to use large meshes, while at the macroscopic level wave propagation is solved maintaining the SEM accuracy and computational efficiency as confirmed by the numerical results.