High-order exponential integration for seismic wave modeling
CoRR(2024)
摘要
Seismic imaging is a major challenge in geophysics with broad applications.
It involves solving wave propagation equations with absorbing boundary
conditions (ABC) multiple times. This drives the need for accurate and
efficient numerical methods. This study examines a collection of exponential
integration methods, known for their good numerical properties on wave
representation, to investigate their efficacy in solving the wave equation with
ABC. The purpose of this research is to assess the performance of these
methods. We compare a recently proposed Exponential Integration based on Faber
polynomials with well-established Krylov exponential methods alongside a
high-order Runge-Kutta scheme and low-order classical methods. Through our
analysis, we found that the exponential integrator based on the Krylov subspace
exhibits the best convergence results among the high-order methods. We also
discovered that high-order methods can achieve computational efficiency similar
to lower-order methods while allowing for considerably larger time steps. Most
importantly, the possibility of undertaking large time steps could be used for
important memory savings in full waveform inversion imaging problems.
更多查看译文
AI 理解论文
溯源树
样例
生成溯源树,研究论文发展脉络
Chat Paper
正在生成论文摘要