title: Two-dimensional method for unconditionally stable elastic wave simulations.
abstract: The stability of the conventional staggered-grid finite-difference time-domain method for elastic wave simulations is limited by the Courant condition and material heterogeneity. Its computational efficiency is significantly hampered when the mesh size is much smaller than a wavelength (for geometric modeling accuracy) and/or with a high impedance contrast. An unconditionally stable alternating direction implicit method is proposed to overcome this issue. It is based on additive operator splitting that renders tri-diagonal matrices for the implicit update of selected field variables. Theoretical analysis of its stability and grid dispersion error are provided. Numerical examples further demonstrate its accuracy and efficiency.
authors.affiliation: Department of Electrical and Computer Engineering, Auburn University, Auburn, Alabama 36849.
authors.forename: Yu
authors.initials: Y
authors.lastname: Shao
authors.affiliation: Department of Electrical and Computer Engineering, Auburn University, Auburn, Alabama 36849.
authors.forename: Myoung
authors.initials: M
authors.lastname: An
authors.affiliation: Department of Electrical and Computer Engineering, Auburn University, Auburn, Alabama 36849.
authors.forename: Shumin
authors.initials: S
authors.lastname: Wang