Energy-conserving 3D elastic wave simulation with finite difference discretization on staggered grids with nonconforming interfaces

by   Longfei Gao, et al.

In this work, we describe an approach to stably simulate the 3D isotropic elastic wave propagation using finite difference discretization on staggered grids with nonconforming interfaces. Specifically, we consider simulation domains composed of layers of uniform grids with different grid spacings, separated by planar interfaces. This discretization setting is motivated by the observation that wave speeds of earth media tend to increase with depth due to sedimentation and consolidation processes. We demonstrate that the layer-wise finite difference discretization approach has the potential to significantly reduce the simulation cost, compared to its counterpart that uses holistically uniform grids. Such discretizations are enabled by summation-by-parts finite difference operators, which are standard finite difference operators with special adaptations near boundaries or interfaces, and simultaneous approximation terms, which are penalty terms appended to the discretized system to weakly impose boundary or interface conditions. Combined with specially designed interpolation operators, the discretized system is shown to preserve the energy-conserving property of the continuous elastic wave equation, and a fortiori ensure the stability of the simulation. Numerical examples are presented to corroborate these analytical developments.



There are no comments yet.


page 18

page 19


Explicit coupling of acoustic and elastic wave propagation in finite difference simulations

We present a mechanism to explicitly couple the finite-difference discre...

Energy conservative SBP discretizations of the acoustic wave equation in covariant form on staggered curvilinear grids

We develop a numerical method for solving the acoustic wave equation in ...

Non-stiff narrow-stencil finite difference approximations of the Laplacian on curvilinear multiblock grids

The Laplacian appears in several partial differential equations used to ...

A high order finite difference method for the elastic wave equations in bounded domains with nonconforming interfaces

We develop a stable finite difference method for the elastic wave equati...

A non-balanced staggered-grid finite-difference scheme for the first-order elastic wave-equation modeling

We introduce an efficient and accurate staggered-grid finite-difference ...

Geometry of martensite needles in shape memory alloys

We study the geometry of needle-shaped domains in shape-memory alloys. N...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.