Inverse problems for hyperbolic partial differential equations
Lead Research Organisation:
University College London
Department Name: Mathematics
Abstract
A typical inverse boundary value problem arises when an object is probed by sending waves to it in order to determine some properties of the object from the measured responses. The proposed research focuses on the case where the measurements are modelled using the acoustic wave equation. The acoustic case can be seen as a prototype case for several geophysical imaging problems. These problems can be local, for example, when the goal is to detect geological deposits, or global when the goal is to chart the internal structure of the Earth.
The aim of the proposed research is to significantly extend the mathematical theory of inverse boundary value problems and to bridge the gap between this theory and practical applications in geophysical imaging. The main theoretical objective is to show that the inverse boundary value problem for the wave equation has a unique solution in the fully coordinate invariant context.
The main practical objective is to develop a new computational method for inverse boundary value problems. The presently used computational methods in geophysical imaging are not based on the mathematical theory of inverse boundary value problems. Bridging the gap between the theory and applications has a potential for a high impact on the latter, and the resulting physical understanding will also guide further mathematical study.
The research brings together ideas from analysis of partial differential equations, differential geometry, and numerical analysis.
The aim of the proposed research is to significantly extend the mathematical theory of inverse boundary value problems and to bridge the gap between this theory and practical applications in geophysical imaging. The main theoretical objective is to show that the inverse boundary value problem for the wave equation has a unique solution in the fully coordinate invariant context.
The main practical objective is to develop a new computational method for inverse boundary value problems. The presently used computational methods in geophysical imaging are not based on the mathematical theory of inverse boundary value problems. Bridging the gap between the theory and applications has a potential for a high impact on the latter, and the resulting physical understanding will also guide further mathematical study.
The research brings together ideas from analysis of partial differential equations, differential geometry, and numerical analysis.
Planned Impact
The proposed research has transformative potential since solving the inverse boundary value problem for the wave equation in a coordinate invariant manner will guide solving other inverse boundary value problems of geometric nature, since it explores a promising new paradigm to solve unique continuation problems, and since it builds bridges between mathematics and geophysical imaging. Moreover, the methods to be developed are based on geometric techniques that have applications beyond inverse boundary value problems, for example, in cosmology.
The fellow has experience in communicating with the researchers in geometry and analysis. To disseminate the results to numerical analysis and geophysical imaging communities the connections of the collaborators will be used. Drawing from the fellow's experience in inter-disciplinary discussions with the Photoacoustic Imaging Group at University College London, in the context of his EPSRC 1st grant project (EP/L026473/1), a plan to communicate the results to geophysicists has been based on construction of computational examples.
The fellow has experience in communicating with the researchers in geometry and analysis. To disseminate the results to numerical analysis and geophysical imaging communities the connections of the collaborators will be used. Drawing from the fellow's experience in inter-disciplinary discussions with the Photoacoustic Imaging Group at University College London, in the context of his EPSRC 1st grant project (EP/L026473/1), a plan to communicate the results to geophysicists has been based on construction of computational examples.
People |
ORCID iD |
Lauri Oksanen (Principal Investigator / Fellow) |
Publications
Alexakis S
(2022)
Lorentzian Calderón problem under curvature bounds
in Inventiones mathematicae
Alexakis S
(2020)
Lorentzian Calderón problem under curvature bounds
Burman E
(2020)
Stability estimate for scalar image velocimetry
Burman E
(2019)
A fully discrete numerical control method for the wave equation
Burman E
(2018)
Primal-Dual Mixed Finite Element Methods for the Elliptic Cauchy Problem
in SIAM Journal on Numerical Analysis
Burman E
(2018)
Fully discrete finite element data assimilation method for the heat equation
in ESAIM: Mathematical Modelling and Numerical Analysis
Burman E
(2020)
A finite element data assimilation method for the wave equation
in Mathematics of Computation
Description | We have developed a method to solve inverse boundary value problems modelling geophysical imaging problems. Our method is guaranteed to converge to the unique, global solution of the inverse problem, but it is unstable and a regularization is needed. At the moment we use a rather naive Tikhonov regularization for certain linear control problems that appear as steps in the method solving the full non-linear problem. The regularization can be seen as a way to trade resolution in the solution to stability in the computation, and a major objective of the ongoing research is to get a better exchange rate by using stabilization techniques from the analysis of finite element methods. The control problems are dual problems to unique continuation problems, and we started approaching the regularization of the hyperbolic control problems by studying first some less complicated unique continuation problems. We developed stabilized finite element theory for unique continuation problems for the heat equation, the Helmholtz equation, the time-domain acoustic wave equation, and convection-diffusion equations. We developed also finite element theory for the time-domain acoustic wave equation (arXiv:1903.02320). In a more theoretical line of research, we have obtained new uniqueness results related to determination of time-dependent vector and scalar potentials in time-domain wave equations in ultrastatic geomeries. We proved also uniqueness result for scalar potentials in the case of stationary geometries (arXiv:1911.04834). In the past year we obtained a major breakthrough result (arXiv:2008.07508) allowing for non-stationary geomeries without any real analytic features. In particular, we can recover scalar potentials in a set of geometries that is open with respect to smooth non-stationary perturbations. |
Exploitation Route | The stabilized finite element methods we have developed can very likely to be adapted to other ill-posed problems. Our focus is on geophysical imaging, but such problems arise in many other applied fields, for instance in material sciences. |
Sectors | Energy Healthcare |
Description | UCL-Inria associate team IMFIBIO |
Organisation | The National Institute for Research in Computer Science and Control (INRIA) |
Country | France |
Sector | Public |
PI Contribution | Collaboration between our group at UCL Maths and COMMEDIA project-team at Inria Paris focusing on Innovative Methods for Forward and Inverse problems in BIO-medical applications. The stabilized finite element methods developed by our team, in the context of this award, will be used in the collaboration. |
Collaborator Contribution | Inria Paris team has developed computational models for fluid-structure-contact interaction and inverse problems related to cardiac electrophysiology and data driven blood flow simulations. |
Impact | None yet. The associate team has been just formally established. |
Start Year | 2020 |
Description | Applied inverse problems conference |
Form Of Engagement Activity | A talk or presentation |
Part Of Official Scheme? | No |
Geographic Reach | International |
Primary Audience | Other audiences |
Results and Impact | Talk in a scientific conference |
Year(s) Of Engagement Activity | 2017 |
Description | Lectures in a summer school at Max Planck Institute, Leipzig |
Form Of Engagement Activity | Participation in an activity, workshop or similar |
Part Of Official Scheme? | No |
Geographic Reach | International |
Primary Audience | Postgraduate students |
Results and Impact | Lectures on inverse problem for hyperbolic partial differential equations |
Year(s) Of Engagement Activity | 2018 |
URL | https://www.mis.mpg.de/calendar/conferences/2018/operators2018.html |
Description | Math + X Symposium on Data Science and Inverse Problems in Geophysics |
Form Of Engagement Activity | A talk or presentation |
Part Of Official Scheme? | No |
Geographic Reach | International |
Primary Audience | Other audiences |
Results and Impact | Talk at a scientific conference |
Year(s) Of Engagement Activity | 2018 |
Description | Recent Advances in Seismic Modeling and Inversion: From Analysis to Applications workshop |
Form Of Engagement Activity | A talk or presentation |
Part Of Official Scheme? | No |
Geographic Reach | International |
Primary Audience | Other audiences |
Results and Impact | Talk at a scientific conference |
Year(s) Of Engagement Activity | 2017 |