Variational Data Assimilation for the Wave Equation in Heterogeneous Media: Numerical Investigation of Stability

Abstract

In recent years, several numerical methods for solving the unique continuation problem for the wave equation in a homogeneous medium with given data on the lateral boundary of the space-time cylinder have been proposed. This problem enjoys Lipschitz stability if the geometric control condition (GCC) is fulfilled, which allows devising optimally convergent numerical methods. In this article, we investigate whether these results carry over to the case in which the medium exhibits a jump discontinuity. Our numerical experiments suggest a positive answer. However, we also observe that the presence of discontinuities in the medium renders the computations far more demanding than in the homogeneous case.

Keywords

Unique continuation / Data assimilation / Wave equation / Discontinuous coefficients / Geometric control condition (GCC) / 65M32 / 65M60 / 35L05

References

[1] Alnaes, M.S., Logg, A., Ølgaard, K.B., Rognes, M.E., Wells, G.N.: Unified form language: a domain-specific language for weak formulations of partial differential equations. ACM Trans. Math. Software 40 (2014). https://doi.org/10.1145/2566630

[2] Baratta, I.A., Dean, J.P., Dokken, J.S., Habera, M., Hale, J.S., Richardson, C.N., Rognes, M.E., Scroggs, M.W., Sime, N., Wells, G.N.: DOLFINx: the next generation FEniCS problem solving environment (2023). https://doi.org/10.5281/zenodo.10447666

[3] Bardos C, Lebeau G, Rauch J. Sharp sufficient conditions for the observation, control, and stabilization of waves from the boundary. SIAM J. Control. Optim., 1992, 30(5): 1024-1065

[4] Burman, E., Feizmohammadi, A., Münch, A., Oksanen, L.: Space time stabilized finite element methods for a unique continuation problem subject to the wave equation. ESAIM Math. Model. Numer. Anal. 55, 969–991 (2021). https://doi.org/10.1051/m2an/2020062

[5] Burman, E., Feizmohammadi, A., Münch, A., Oksanen, L.: Spacetime finite element methods for control problems subject to the wave equation. ESAIM Control Optim. Calc. Var. 29, 41 (2023). https://doi.org/10.1051/cocv/2023028

[6] Burman E, Delay G, Ern A. The unique continuation problem for the wave equation discretized with a high-order space-time nonconforming method. Numer. Math., 2025, 157: 1259-1284

[7] Burman E, Feizmohammadi A, Oksanen L. A finite element data assimilation method for the wave equation. Math. Comput., 2020, 89(324): 1681-1709

[8] Burman, E., Preuss, J.: Unique continuation for the wave equation based on a discontinuous Galerkin time discretization. IMA J. Numer. Anal. (2025). https://doi.org/10.1093/imanum/draf036

[9] Burman E, Preuss J. Unique continuation for an elliptic interface problem using unfitted isoparametric finite elements. SMAI J. Comput. Math., 2025, 11: 165-202

[10] Cîndea N, Münch A. Inverse problems for linear hyperbolic equations using mixed formulations. Inverse Prob., 2015, 31(7 075001

[11] Dahmen, W., Monsuur, H., Stevenson, R.: Least squares solvers for ill-posed PDEs that are conditionally stable. ESAIM Math. Model. Numer. Anal. 57(4), 2227–2255 (2023). https://doi.org/10.1051/m2an/2023050

[12] Filippas, S.: Whispering gallery modes for a transmission problem. ESAIM: COCV 29, 75 (2023). https://doi.org/10.1051/cocv/2023067

[13] Filippas, S.: Quantitative unique continuation for wave operators with a jump discontinuity across an interface and applications to approximate control. Bull. Soc. Math. France 152(2), 199–294 (2024). https://doi.org/10.24033/bsmf.2889

[14] Gagnon, L.: Sufficient conditions for the controllability of wave equations with a transmission condition at the interface. ESAIM: COCV 29, 51 (2023). https://doi.org/10.1051/cocv/2023034

[15] Lehrenfeld, C., Heimann, F., Preuß, J., Wahl, H.: ngsxfem: Add-on to NGSolve for geometrically unfitted finite element discretizations. J. Open Source Softw. 6(64), 3237 (2021). https://doi.org/10.21105/joss.03237

[16] Mattsson K, Ham F, Iaccarino G. Stable and accurate wave-propagation in discontinuous media. J. Comput. Phys., 2008, 22719): 8753-8767

[17] Mishra S, Molinaro R. Estimates on the generalization error of physics-informed neural networks for approximating a class of inverse problems for PDEs. IMA J. Numer. Anal., 2021, 42(2): 981-1022

[18] Schöberl J. NETGEN an advancing front 2D/3D-mesh generator based on abstract rules. Comput. Vis. Sci., 1997, 1(1): 41-52

[19] Schöberl, J.: C++11 implementation of finite elements in NGSolve. Technical report, ASC-2014-30. Institute for Analysis and Scientific Computing, TU Wien. http://hdl.handle.net/20.500.12708/28346

[20] Scroggs, M.W., Baratta, I.A., Richardson, C.N., Wells, G.N.: Basix: a runtime finite element basis evaluation library. J. Open Source Softw. 7(73), 3982 (2022). https://doi.org/10.21105/joss.03982

[21] Scroggs, M.W., Dokken, J.S., Richardson, C.N., Wells, G.N.: Construction of arbitrary order finite element degree-of-freedom maps on polygonal and polyhedral cell meshes. ACM Trans. Math. Software 48(2), 18 (2022). https://doi.org/10.1145/3524456

[22] Stolk, C.: On the Modeling and Inversion of Seismic Data. Ph.D. thesis. Utrecht University, Utrecht(2000). https://dspace.library.uu.nl/handle/1874/855

Funding

Georg-August-Universität Göttingen (1018)

RIGHTS & PERMISSIONS

The Author(s)

Just Accepted

This article has successfully passed peer review and final editorial review, and will soon enter typesetting, proofreading and other publishing processes. The currently displayed version is the accepted final manuscript. The officially published version will be updated with format, DOI and citation information upon launch. We recommend that you pay attention to subsequent journal notifications and preferentially cite the officially published version. Thank you for your support and cooperation.