1887

### Abstract

Summary

The explicit finite-difference (FD) scheme is popular for solving the wave equation due to its simplicity in numerical implementation. However, its maximum time step is strictly restricted by the Courant-Friedrichs-Lewy (CFL) stability limit, which leads to a heavy computational burden in the presence of small-scale structures and high-velocity targets.

We extend the CFL stability limit of the explicit FD scheme using the eigenvalue perturbation. The eigenvalue perturbation can normalize the unstable eigenvalues and guarantee the stability of the update matrix by using an arbitrary time step. The update matrix can be preprocessed before the numerical simulation, thus retaining the computational efficiency well. We further incorporate the forward time-dispersion transform (FTDT) and the inverse time-dispersion transform (ITDT) to reduce the time dispersion error caused by using an unusually large time step. Our numerical experiments indicate that the combination of the eigenvalue perturbation, the FTDT method, and the ITDT method can simulate highly accurate waveforms when applying a time step beyond the CFL stability limit. The time step can be extended even toward the Nyquist limit. This means that we could save many iteration steps without suffering from time-dispersion error and stability problems.

/content/papers/10.3997/2214-4609.201901541
2019-06-03
2019-12-16

References

1. Courant, R., Friedrichs, K., and Lewy, H.
[1928] Über die partiellen Differenzengleichungen der mathematischen Physik. Mathematische Annalen (In German), 100(1), 32–74.
2. Gao. Y., Zhang, J., YaoZ.
[2018] Removing the stability limit of the explicit finite-difference scheme with eigenvalue perturbation. Geophysics, 83(6), A93–A98.
3. Koene, E. F. M., Robertsson, J. O. A., Broggini, F., and Andersson, F.
[2018] Liminating time dispersion from seismic wave modelling. Geophysical Journal International, 213, 169–180.
4. Li, X., Sarris, C. D., and Triverio, P.
[2014] Overcoming the FDTD stability limit via model order reduction and eigenvalue perturbation. Microwave Symposium (IMS), 2014 IEEE MTT-S International. IEEE. Paper read at Microwave Symposium. 1–3.
5. Liu, Y., and Sen, M. K.
[2009] A practical implicit finite-difference method: Examples from seismic modelling. Journal of Geophysics and Engineering, 6, 231–249.
6. Wang, M., and Xu, S.
[2015] Finite-difference time dispersion transforms for wave propagation. Geophysics, 80(6), WD19–WD25.