×

Discretization of multipole sources in a finite difference setting for wave propagation problems. (English) Zbl 1452.65194

Summary: Seismic sources are commonly idealized as point-sources due to their small spatial extent relative to seismic wavelengths. The acoustic isotropic point-radiator is inadequate as a model of seismic wave generation for seismic sources that are known to exhibit directivity. Therefore, accurate modeling of seismic wavefields must include source representations generating anisotropic radiation patterns. Such seismic sources can be modeled as multipoles, that is, a time-dependent linear combination of spatial derivatives of the spatial delta function. Since the solutions of linear hyperbolic systems with point-source right hand sides are necessarily singular, standard results on convergence of grid-based numerical methods (finite difference or finite element) do not imply convergence of numerical solutions. We present a method for discretizing multipole sources in a finite difference setting, an extension of the moment matching conditions developed for the Dirac delta function in other applications, along with numerical evidence demonstrating the accuracy of these approximations. Using this analysis, we develop a weak convergence theory for the discretization of a family of symmetric hyperbolic systems of first-order partial differential equations, with singular source terms, solved via staggered-grid finite difference methods: we show that grid-independent space-time averages of the numerical solutions converge to the same averages of the continuum solution, and provide an estimate for the error in terms of moment matching and truncation error conditions. Numerical experiments confirm this result, but also suggest a stronger one: optimal convergence rates appear to be achieved point-wise in space away from the source.

MSC:

65M12 Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
35L05 Wave equation
35A21 Singularity in context of PDEs
PDFBibTeX XMLCite
Full Text: DOI Link

References:

[1] Blazek, K.; Stolk, C. C.; Symes, W., A mathematical framework for inverse wave problems in heterogeneous media, Inverse Probl., 29, 065001, 1-34 (2013) · Zbl 1273.35309
[2] Brenner, P.; Thomée, V.; Wahlbin, L. B., Besov Spaces and Applications to Difference Methods for Initial Value Problems, vol. 434 (2006), Springer
[3] Cohen, G. C., Higher Order Numerical Methods for Transient Wave Equations (2002), Springer: Springer New York · Zbl 0985.65096
[4] Courant, R.; Hilbert, D., Methods of Mathematical Physics, vol. II (1962), Wiley-Interscience: Wiley-Interscience New York · Zbl 0099.29504
[5] Epanomeritakis, I.; Akcelik, V.; Ghattas, O.; Bielak, J., A Newton-CG method for large-scale three-dimensional elastic full-waveform seismic inversion, Inverse Probl., 24, Article 034015 pp. (2009) · Zbl 1142.65052
[6] Gustafsson, B.; Mossberg, E., Time compact high order difference methods for wave propagation, SIAM J. Sci. Comput., 26, 259-271 (2004) · Zbl 1075.65112
[7] Gustafsson, B.; Wahlund, P., Time compact difference methods for wave propagation in discontinuous media, SIAM J. Sci. Comput., 26, 272-293 (2004) · Zbl 1077.65092
[8] Gustafsson, B.; Wahlund, P., Time compact high order difference methods for wave propagation, 2D, J. Sci. Comput., 25, 195-211 (2004) · Zbl 1203.65136
[9] Hörmander, L., Linear Partial Differential Operators (1969), Springer-Verlag: Springer-Verlag New York · Zbl 0177.36401
[10] Hosseini, B.; Nigam, N.; Stockie, J. M., On regularizations of the Dirac delta distribution, J. Comput. Phys., 305, 423-447 (2016) · Zbl 1349.35008
[11] Hu, W.; Abubakar, A.; Habashy, T., Application of the nearly perfectly matched layer in acoustic wave modeling, Geophysics, 72, SM169-SM176 (2007)
[12] Komatisch, D.; Barnes, C.; Tromp, J., Simulation of anisotropic wave propagation based upon a spectral element method, Geophysics, 65, 1251-1260 (2000)
[13] Königsberger, K., Analysis, vol. 2 (2013), Springer-Verlag
[14] Lax, P. D., Hyperbolic Partial Differential Equations, Courant Lecture Notes (2006), American Mathematical Society: American Mathematical Society Providence · Zbl 1113.35002
[15] Layton, W. J., Simplified \(l^\infty\) estimates for difference schemes for partial differential equations, Proc. Am. Math. Soc., 86, 3, 491-495 (1982) · Zbl 0497.65049
[16] Leveque, R., Finite Difference Methods for Ordinary and Partial Differential Equations (2007), SIAM: SIAM Philadelphia · Zbl 1127.65080
[17] Moczo, P.; Robertsson, J. O.A.; Eisner, L., The finite-difference time-domain method for modeling of seismic wave propagation, Adv. Geophys., 48, 421-516 (2006)
[18] Padula, A. D.; Symes, W.; Scott, S. D., A software framework for the abstract expression of coordinate-free linear algebra and optimization algorithms, ACM Trans. Math. Softw., 36, 8:1-8:36 (2009) · Zbl 1364.65130
[19] Peskin, C. S., The immersed boundary method, Acta Numer., 11, 479-517 (2002) · Zbl 1123.74309
[20] Petersson, N. A.; Sjögreen, B., Stable grid refinement and singular source discretization for seismic wave simulations, Commun. Comput. Phys., 8, 1074-1110 (2010) · Zbl 1364.86010
[21] Petersson, N. A.; O’Reilly, O.; Sjögreen, B.; Bydlon, S., Discretizing singular point sources in hyperbolic wave propagation problems, J. Comput. Phys., 321, 532-555 (2016) · Zbl 1349.65329
[22] Richtmyer, R.; Morton, K., Difference Methods for Initial-Value Problems (1967), Interscience Publishers, a division of John Wiley & Sons: Interscience Publishers, a division of John Wiley & Sons New York · Zbl 0155.47502
[23] Shearer, P. M., Introduction to Seismology (2009), Cambridge University Press
[24] Symes, W.; Sun, D.; Enriquez, M., From modelling to inversion: designing a well-adapted simulator, Geophys. Prospect., 59, 814-833 (2011)
[25] Tornberg, A.; Engquist, B., Regularization techniques for numerical approximation of pdes with singularities, J. Sci. Comput., 19, 527-552 (2003) · Zbl 1035.65085
[26] Tornberg, A.; Engquist, B., Numerical approximation of singular source terms in differential equations, J. Comput. Phys., 200, 462-488 (2004) · Zbl 1115.76392
[27] Virieux, J., SH-wave propagation in heterogeneous media: velocity stress finite-difference method, Geophysics, 49, 1933-1957 (1984)
[28] Waldén, J., On the approximation of singular source terms in differential equations, Numer. Methods Partial Differ. Equ., 15, 4, 503-520 (1999) · Zbl 0938.65112
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.