Open AccessExact extrapolation and immersive modelling with finite-difference injection
Author(s) -
DirkJan van Manen,
Xun Li,
Marlies Vasmel,
Filippo Broggini,
Johan O. A. Robertsson
Publication year - 2020
Publication title -
geophysical journal international
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0956-540X
DOI - 10.1093/gji/ggaa317
Subject(s) - extrapolation , helmholtz equation , finite difference , helmholtz free energy , representation (politics) , mathematical analysis , mathematics , range (aeronautics) , boundary (topology) , exact solutions in general relativity , wave equation , computer science , boundary value problem , physics , materials science , quantum mechanics , politics , political science , law , composite material
SUMMARY In numerical modelling of wave propagation, the finite-difference (FD) injection method enables the re-introduction of simulated wavefields in model subdomains with machine precision, enabling the efficient calculation of waveforms after localized model alterations. By rewriting the FD-injection method in terms of sets of equivalent sources, we show how the same principles can be applied to achieve on-the-fly wavefield extrapolation using Kirchhoff–Helmholtz (KH)-like integrals. The resulting extrapolation methods are numerically exact when used in conjunction with FD-computed Green’s functions. Since FD injection only relies on the linearity of the wave equation and compactness of FD stencils in space, the methods can be applied to both staggered and non-staggered discretizations with arbitrary-order spatial operators. Examples for both types of discretizations show how these extrapolators can be used to truncate models with exact absorbing or immersive boundary conditions. Such immersive modelling involves the evaluation of KH-type extrapolation and representation integrals in the same simulation, which include the long-range interactions missing from conventional FD injection.