Abstract
The solution in terms of spherical harmonics to acoustic and elastic wave propagation problems in spherically stratified media is well established. However, except for very low ka values, a direct numerical implementation is unstable. Propagator matrix approaches have been applied for moderate ka values, whereas environment transformations have been developed to represent the spherical stratification by a plane stratification at higher ka values, in particular in relation to long‐range seismic propagation in the solid earth. Here, it is demonstrated that, by introducing a proper normalization of the spherical Bessel functions, unconditional stability can be obtained by using the Global Matrix approach [H. Schmidt and F. B. Jensen, J. Acoust. Soc. Am. 77, 813–825 (1985)]. This SAFARI code has therefore been modified to treat propagation in spherically stratified elastic media. The algorithm is stable at least up to ka = 103, allowing modeling of stratifications alternating between high‐ and low‐speed layer...
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call