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.

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.
 
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