Based on the formula for stability of finite-difference methods with second-order in time and general-order in space for the scalar wave equation, I obtain a stability formula for Lax-Wendroff methods with fourth-order in time and general-order in space. Unlike the formula for methods with second-order in time, this formula depends on two parameters: one parameter is related to the weights for approximations of second spatial derivatives; the other parameter is related to the weights for approximations of fourth spatial derivatives. When discretizing the mixed derivatives properly, the formula can be generalized to the case where the spacings in different directions are different. This formula can be useful in high-accuracy seismic modeling using the scalar wave equation on rectangular grids, which involves both high-order spatial discretizations and high-order temporal approximations. I also prove the instability of methods obtained by applying high-order finite-difference approximations directly to the second temporal derivative, and this result solves the “Bording’s conjecture.”
A step-by-step approach in the time-domain BEM formulation for the scalar wave equation
