International Journal of Mathematics and Mathematical Sciences
Volume 2008 (2008), Article ID 291968, 31 pages
doi:10.1155/2008/291968

Higher-order splitting method for elastic wave propagation

Jürgen Geiser

Institut für Mathematik, Humboldt Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany

Abstract

Motivated by seismological problems, we have studied a fourth-order split scheme for the elastic wave equation. We split in the spatial directions and obtain locally one-dimensional systems to be solved. We have analyzed the new scheme and obtained results showing consistency and stability. We have used the split scheme to solve problems in two and three dimensions. We have also looked at the influence of singular forcing terms on the convergence properties of the scheme.