INVERSE PROBLEMS

Keywords: inverse problems

## Full Waveform Inversion

Most optimization approaches for reconstruction of velocity models used in Full Waveform Inversion
(FWI) are based on smooth techniques such as the Tikhonov regularization. However, realistic velocity
profiles often present various discontinuities, sharp interfaces and high contrasts, as in the important
case of the presence of salt bodies. Such discontinuities arise from the fact that sound waves travel with
greater velocity inside salt bodies compared to the neighboring sediments.
In this work we propose to use nonsmooth optimization techniques for the reconstruction of sharp
interfaces. We consider a simplified setting where the velocity is piecewise constant, with known
constant but distinct values in the salt body and in the sediment region. In this way, the optimization
problem is recast as a shape optimization problem, where the
interface of the salt region becomes the unknown. The problem is formulated as the minimization of
a tracking-type cost functional with respect to the geometry of the interface, and the evolution of the
interface is performed using a level set method.

> Reconstruction of Sharp Interfaces in Time-Domain Full Waveform Inversion

> Reconstruction of Sharp Interfaces in Time-Domain Full Waveform Inversion

## The inverse source problem

The inverse source problem consists in reconstructing a mass distribution in a geometrical domain from boundary measurements of the associated potential and its normal derivative.
In this paper the inverse source problem is reformulated as a topology optimization problem, where the support of the mass distribution is the unknown variable.
The Kohn-Vogelius functional is minimized. It measures the misfit between the solutions of two auxiliary problems containing information about the boundary measurements.
The Newtonian potential is used to complement the unavailable information on the hidden boundary.
The resulting topology optimization algorithm is based on an analytic formula for the variation of the Kohn-Vogelius functional with respect to a class of mass distributions consisting of a finite number of ball-shaped trial anomalies.
Finally, in order to show the effectiveness of the devised reconstruction algorithm, some numerical experiments in two and three spatial dimensions are presented.

> A new reconstruction method for the inverse source problem from partial boundary measurements

> A new reconstruction method for the inverse potential problem

> A Non-Iterative Method for the Inverse Potential Problem Based on the Topological Derivative

> A new reconstruction method for the inverse source problem from partial boundary measurements

> A new reconstruction method for the inverse potential problem

> A Non-Iterative Method for the Inverse Potential Problem Based on the Topological Derivative

## Numerical algorithms for inverse problems

Two approaches are proposed for solving inverse problems in shape optimization. We are looking
for the unknown position of a small hole in a domain . First, the asymptotic analysis of the underlying
p.d.e. defined in a perturbed domain is performed and the so-called topological derivative is defined. Then,
in the first approach, the self-adjoint extensions of elliptic operators are used to model the solution of a
partial differential equation defined in the singularly perturbed domain. A least-square functional is then
minimized to identify the hole. In the second approach, neural networks are used to determine the inverse
of the mapping which associates a set of shape functionals to the position of the unknown hole. In both
approaches the topological derivatives are used to approximate the shape functionals.

> Numerical algorithms for an inverse problem in shape optimization

> Numerical algorithms for an inverse problem in shape optimization

## Self-adjoint extensions

Self-adjoint extensions of elliptic operators are used to model the solution of
a partial differential equation defined in a singularly perturbed domain. The asymptotic expansion
of the solution of a Laplacian with respect to a small parameter " is first performed
in a domain perturbed by the creation of a small hole. The resulting singular perturbation is
approximated by choosing an appropriate self-adjoint extension of the Laplacian, according
to the previous asymptotic analysis. The sensitivity with respect to the position of the
center of the small hole is then studied for a class of functionals depending on the domain.
A numerical application for solving an inverse problem is presented. Error estimates are
provided and a link to the notion of topological derivative is established.

> Using self-adjoint extensions in shape optimization

> Using self-adjoint extensions in shape optimization