Hadamard integrators for wave equations in time and frequency domain: Eulerian formulations via butterfly algorithms (2401.01423v2)
Abstract: Starting from the Kirchhoff-Huygens representation and Duhamel's principle of time-domain wave equations, we propose novel butterfly-compressed Hadamard integrators for self-adjoint wave equations in both time and frequency domain in an inhomogeneous medium. First, we incorporate the leading term of Hadamard's ansatz into the Kirchhoff-Huygens representation to develop a short-time valid propagator. Second, using the Fourier transform in time, we derive the corresponding Eulerian short-time propagator in frequency domain; on top of this propagator, we further develop a time-frequency-time (TFT) method for the Cauchy problem of time-domain wave equations. Third, we further propose the time-frequency-time-frequency (TFTF) method for the corresponding point-source Helmholtz equation, which provides Green's functions of the Helmholtz equation for all angular frequencies within a given frequency band. Fourth, to implement TFT and TFTF methods efficiently, we introduce butterfly algorithms to compress oscillatory integral kernels at different frequencies. As a result, the proposed methods can construct wave field beyond caustics implicitly and advance spatially overturning waves in time naturally with quasi-optimal computational complexity and memory usage. Furthermore, once constructed the Hadamard integrators can be employed to solve both time-domain wave equations with various initial conditions and frequency-domain wave equations with different point sources. Numerical examples for two-dimensional wave equations illustrate the accuracy and efficiency of the proposed methods.
- Dover Publications, Inc., New York. (1965)
- Babich, V.M.: The short wave asymptotic form of the solution for the problem of a point source in an inhomogeneous medium. USSR Computational Mathematics and Mathematical Physics 5(5), 247–251 (1965)
- SIAM Review 42, 451–484 (2000)
- J. Comput. Phys. 261, 36–64 (2014)
- J. Comput. Phys. 59, 396–404 (1985)
- SIAM Multiscale Model. Simul. 7, 1727–1750 (2009)
- John Wiley-Sons (1962)
- J. Comput. Phys. 228, 6440–6455 (2009)
- Academic Press, New York and London (1964)
- Hadamard, J.: Lectures on Cauchy’s Problem in Linear Partial Differential Equations. Yale University Press; (reprinted Dover Publications, New York 1952) (1923)
- J. Sci. Comput. 88, 54 (2021)
- SIAM J. Sci. Comput. 21, 2126–2143 (2000)
- J. Comput. Phys. 196, 367–391 (2004)
- J. Comput. Phys. 348, 108–138 (2017)
- J. Comput. Phys. 228, 2951–2977 (2009)
- J. Comput. Phys. 229, 8888–8917 (2010)
- Geophysics 72, SM61–SM76 (2007)
- Methods Appl. Analy. 21, 031–066 (2014)
- Comm. in Comp. Phys. 8, 758–796 (2010)
- SIAM Journal on Scientific Computing 39(2), A503–A531 (2017)
- Multiscale Modeling & Simulation 13(2), 714–732 (2015)
- J. Comput. Phys. 115, 200–212 (1994)
- IEEE Antennas and Wireless Propagation Letters 16, 1179–1183 (2016)
- SIAM Multiscale Model. Simul. 21, 269–308 (2023)
- SIAM J. Sci. Comput. 43, A883–A907 (2021)
- SIAM J. Multiscale Model. Simul. 14(3), 1089–1122 (2016)
- J. Comput. Phys. 313, 478–510 (2016)
- SIAM J. Multiscale Model. Simul. 16, 727–751 (2018)
- J. Comput. Phys. 270, 378–401 (2014)
- SIAM J. Numer. Analy. 52, 23–44 (2014)
- IEEE Transactions on Antennas and Propagation 44(8), 1086–1093 (1996)
- SIAM J. Numer. Analy. 28, 907–922 (1991)
- SIAM J. Multiscale Modeling and Simulation 16, 595–636 (2016)
- SIAM J. Numer. Analy. 59, 2536–2570 (2021)
- SIAM Multiscale Model. Simul. 19, 46–86 (2021)
- Minimax Theory and its Applications 8, 171–212 (2023)
- Geophysics 67, 167–176 (2002)
- J. Comput. Phys. 229, 7848–7873 (2010)
- SIAM J. Multiscale Modeling and Simulation 8, 1803–1837 (2010)
- Journal of Scientific Computing 67, 883–908 (2016)
- J. Sci. Comp. 19, 501–526 (2003)
- J. Comput. Phys. 228, 8856–8871 (2009)
- SIAM J. Multiscale Modeling and Simulation 6, 688–709 (2007)
- Journal of Computational Physics 497, 112637 (2024)
- Geophys. J. Internat. 162, 1–8 (2005)
- Zhao, H.K.: Fast sweeping method for eikonal equations. Math. Comp. 74, 603–627 (2005)
Collections
Sign up for free to add this paper to one or more collections.
Paper Prompts
Sign up for free to create and run prompts on this paper using GPT-5.