Fast wave computation via Fourier integral operators
Name
Demanet_Fast Wave.pdf
Size
1.73 MB
Format
Adobe PDF
Checksum (MD5)
ea03acb4b9f5ec61d68a46b68e361082
Author(s)
Demanet, Laurent
Ying, Lexing
Date Issued
February 2012
Journal
Mathematics of Computation
Publisher
American Mathematical Society (AMS)
Citation
Demanet, Laurent, and Lexing Ying. “Fast Wave Computation via Fourier Integral Operators.” Mathematics of Computation 81.279 (2012): 1455–1486. Web. 11 Apr. 2012.
Version
Author's final manuscript
Abstract
This paper presents a numerical method for ``time upscaling'' wave equations, i.e., performing time steps not limited by the Courant-Friedrichs-Lewy (CFL) condition. The proposed method leverages recent work on fast algorithms for pseudodifferential and Fourier integral operators (FIO). This algorithmic approach is not asymptotic: it is shown how to construct an exact FIO propagator by 1) solving Hamilton-Jacobi equations for the phases, and 2) sampling rows and columns of low-rank matrices at random for the amplitudes. The setting of interest is that of scalar waves in two-dimensional smooth periodic media (of class C∞ over the torus), where the bandlimit $ N$ of the waves goes to infinity. In this setting, it is demonstrated that the algorithmic complexity for solving the wave equation to fixed time T ≃ 1 can be as low as O(N [superscript 2] log N) with controlled accuracy. Numerical experiments show that the time complexity can be lower than that of a spectral method in certain situations of physical interest.
MIT Department
Massachusetts Institute of Technology. Department of Mathematics
Terms of Use
Creative Commons Attribution-Noncommercial-Share Alike 3.0
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1090/S0025-5718-2012-02557-9
