On the required number of electrodes for uniqueness and convex reformulation in an inverse coefficient problem

Andrej Brojatsch,Bastian Harrach
2024-11-01
Abstract:We introduce a computer-assisted proof for uniqueness and global reconstruction for the inverse Robin transmission problem, where the corrosion function on the boundary of an interior object is to be determined from current-voltage measurements on the boundary of an outer domain. We consider the shunt electrode model where, in contrast to the standard Neumann boundary condition, the applied electrical current is only partially known. The aim is to determine the corrosion coefficient with a finite number of measurements. In this paper, we present a numerically verifiable criterion that ensures unique solvability of the inverse problem, given a desired resolution. This allows us to explicitly determine the required number and position of the required electrodes. Furthermore, we will present an error estimate for noisy data. By rewriting the problem as a convex optimisation problem, our aim is to develop a globally convergent reconstruction algorithm.
Optimization and Control,Analysis of PDEs,Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: in the inverse Robin transmission problem, determine the number and position of electrodes required to ensure uniqueness and global reconstruction. Specifically, the author studies how to determine the corrosion function on the boundary of an internal object from current - voltage measurements on the boundary of an external domain. The paper mainly focuses on the shunt electrode model, in which part of the applied current is known. ### Specific description of the problem 1. **Inverse Robin transmission problem**: - The goal is to reconstruct the corrosion coefficient on the boundary of an internal object from current - voltage measurements on the boundary of an external domain. - Traditional methods usually use regularized data - fitting methods, but due to the non - convexity of the residual function, these methods are difficult to achieve global convergence. 2. **Shunt electrode model**: - Different from the standard Neumann boundary condition, the shunt electrode model assumes that part of the applied current is known. - This model is closer to the actual situation, but it increases the complexity of the problem. 3. **Uniqueness and global reconstruction**: - The paper proposes a numerically verifiable criterion to ensure that the inverse problem has a unique solution at a given resolution. - By reformulating the problem as a convex optimization problem, a globally convergent reconstruction algorithm is developed. ### Main contributions - **Explicit calculation of the required number of electrodes**: For the first time, a method is proposed that can explicitly calculate the required number of electrodes to ensure a unique solution to the problem of reconstructing non - linear spatially - dependent coefficients. - **Error estimation**: Provides error estimates for noisy data. - **Theoretical generalization**: Extends existing research results from idealized models to the shunt electrode model and establishes a more accurate uniqueness criterion. ### Mathematical formula representation To ensure the correctness and readability of the formulas, the following are some of the key formulas involved in the paper: - **Forward mapping**: \[ F(\gamma)I = U, \quad \text{where} \quad u(I)_\gamma \text{ satisfies equation (1)} \] where \( F: \mathbb{R}^n_+ \to S^m_+ \), and \( S^m_+ \) represents the space of symmetric positive - definite matrices. - **Derivative relationship**: \[ I^T (F'(\gamma) \delta) J = -\sum_{i = 1}^n \delta_i \int_{\Gamma_i} u(I)_\gamma u(J)_\gamma \, dx \] - **Uniqueness and stability conditions**: \[ F'(\gamma)d \not\succeq 0 \quad \text{for all} \quad d \in \mathbb{R}^n\setminus\{0\} \] where \( \succeq \) represents the Loewner order relation. Through these formulas and methods, the paper provides a systematic method to determine the number and position of electrodes required, thereby ensuring the uniqueness and global reconstruction of the inverse problem.