3-D elastic simulation of wave propagation using cylindrical coordinates

The paper presents a method for the simulation of wave propagation on a three-dimensional cylindrical grid using a multi-domain approach, where the numerical grid is divided into a number of subdomains consisting of concentric cylindrical shells. Discretization of the wavefield is carried out in r, θ and z, where r is the radial distance from the center, θ is the polar angle and z is depth. Spatial derivative approximations are calculated using the Chebychev expansion along the radial direction, and using the Fourier expansion along the angular and the vertical directions. The equations of conservation of momentum and the stress-strain relation for an isotropic elastic medium are the basis for the numerical technique. -from Authors