Uniqueness of an inverse electromagnetic coefficient problem with partial boundary data and its numerical resolution through an iterated sensitivity equation

Jérémy Heleine

2024-10-03

Abstract:In this paper we study an inverse boundary value problem for Maxwell's equations. The goal is to reconstruct perturbations in the refractive index of the medium inside an object from the knowledge of the tangential trace of an electric field on a part of the boundary of the domain. We first provide a uniqueness result for this inverse problem. Then, we propose a complete procedure to reconstruct numerically the perturbations, based on the minimization of a cost functional involving an iterated sensitivity equation.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to reconstruct the perturbation of the refractive index inside the medium through partial boundary data. Specifically, the author studied an inverse problem, that is, to reconstruct the change of the refractive index inside the medium from the known tangential traces of the electric field in a partial area of the object boundary.

### Problem Background

1. **Physical Background**
- The research object is the inverse boundary value problem of Maxwell's equations.
- The goal is to reconstruct the change of the refractive index inside the medium from partial boundary measurement data.
- This technique has potential uses in medical imaging (such as stroke diagnosis) and industrial applications.

2. **Mathematical Model**
- Let \(\Omega\) be a bounded and simply - connected \(p\)-dimensional region with the boundary \(\Gamma=\partial\Omega\).
- The magnetic permeability is the vacuum magnetic permeability \(\mu_0\).
- The permittivity and conductivity are \(\varepsilon(x)\) and \(\sigma(x)\) respectively, and the complex refractive index \(\kappa(x)\) is defined as:
\[
\kappa(x)=\frac{\varepsilon(x)+i\frac{\sigma(x)}{\omega}}{\varepsilon_0}
\]
- The electric field \(E\) satisfies the time - harmonic Maxwell equation:
\[
\nabla\times(\nabla\times E)-k^2\kappa E = 0
\]
where \(k = \omega\sqrt{\mu_0\varepsilon_0}\) is the wave number.

3. **Inverse Problem Description**
- Reconstruct the refractive index in the entire domain \(\Omega\) from the measurement data on the partial boundary \(\Gamma_0\subset\Gamma\).
- The measurement data on the accessible partial boundary \(\Gamma_0\) is used to reconstruct the information of the inaccessible part \(\Gamma_1=\Gamma\setminus\Gamma_0\).

### Main Contributions of the Paper

1. **Uniqueness Results**
- Provide a uniqueness proof of the solution to the inverse problem in the case of partial data without imposing specific geometric or physical conditions on the inaccessible boundary.

2. **Numerical Method**
- Propose a new method based on the iterative sensitivity equation to numerically solve this inverse problem.
- Improve the quasi - reversible method by introducing an auxiliary unknown \(F = \nabla\times E\) to handle the error in the high - frequency case more accurately.

3. **Data Completion**
- Use the iterative quasi - reversible method to transfer the partial boundary data to the artificial boundary \(\Gamma_{\text{int}}\), thereby obtaining complete boundary data.

4. **Regularization Strategy**
- Introduce a linearized cost function and reduce the computational complexity through Taylor expansion approximation.

### Summary
This paper aims to solve the inverse problem of reconstructing the change of the refractive index inside the medium from partial boundary measurement data, provides a uniqueness proof of the solution and proposes a new numerical solution method. By introducing the iterative sensitivity equation and the improved quasi - reversible method, the paper effectively deals with data completion and regularization problems and improves the accuracy and efficiency of the algorithm.