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.


Article metrics loading...

Loading full text...

Full text loading...


  1. Courant, R., Friedrichs, K., and Lewy, H.
    [1928] Über die partiellen Differenzengleichungen der mathematischen Physik. Mathematische Annalen (In German), 100(1), 32–74.
    [Google Scholar]
  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.
    [Google Scholar]
  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.
    [Google Scholar]
  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.
    [Google Scholar]
  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.
    [Google Scholar]
  6. Wang, M., and Xu, S.
    [2015] Finite-difference time dispersion transforms for wave propagation. Geophysics, 80(6), WD19–WD25.
    [Google Scholar]

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error