We propose compact implicit finite-difference (FD) schemes in time-space domain based on second-order FD approximation for accurate solution of the acoustic wave equation in 1D, 2D, and 3D. Our method is based on weighted linear combination of the second-order FD operators with different spatial orientations to mitigate numerical error anisotropy and weighted averaging of the mass acceleration term over the grid-points of the second-order FD stencil to reduce the overall numerical dispersion error. We present derivation of the schemes for 1D, 2D, and 3D cases and obtain their corresponding dispersion equations, then we find optimum weights by optimization of the time-space domain dispersion function and finally tabulate the optimized weights for each case. We analyze the numerical dispersion, stability and convergence rates of the proposed schemes and compare their numerical dispersion characteristics with the standard high-order ones. We also discuss efficient solution of the system of equations associated with the proposed implicit schemes using conjugate gradient method. The comparison of dispersion curves and the numerical solutions with the analytical and the pseudo-spectral solutions reveals that the proposed schemes have better performance than the standard spatial high-order schemes and remain stable for relatively large time-steps.
Scalar Wave Equation Modeling with Time-Space Domain Dispersion-Relation-Based Staggered-Grid Finite-Difference Schemes
