Optimally accurate second-order time-domain finite-difference scheme for acoustic, electromagnetic, and elastodynamic wave modeling: One-dimensional case

Numerical methods are extremely useful in solving real-life problems with complex materials and geometries. However, numerical methods in the time domain suffer from artificial numerical dispersion. Standard numerical techniques which are second-order in space and time, like the conventional rectangle nite-difference three point (FD3) method [5], rectangle nite-difference time-domain (FDTD) method [4], and rectangle nite integration technique (FIT) [2], provide estimates of the error of discretized numerical operators rather than the error of the numerical solutions computed using these operators [3]. Here, optimally accurate time-domain (TD) rectangle nite-difference (FD) operators which are second-order in time as well as in space are derived [1]. Optimal accuracy means the greatest attainable accuracy for a particular type of scheme, e.g., second-order FD, for some particular grid spacing. The modified FD scheme - FD modified: FDM - presented here attains reduction of numerical dispersion almost by a factor of 40 compared to the FD3, FDTD, and FIT. The CPU time for the FDM scheme is twice of that required by FD3 method. The modified operators lead to an implicit scheme, which is approximated by a predictor-corrector scheme yielding a two step explicit scheme. The possibility of extending this method to staggered grid approach is also presented. Finally the comparison between analytical solution, FDTD/FIT method, FD3 method and FDM scheme with simulation results are depicted. Further examples are given in the presentation.