Positive Jacobian constraints for elliptic boundary value problems with piecewise-regular coefficients arising from multi-wave inverse problems

Yves Capdeboscq,Tianrui Dai

DOI: https://doi.org/10.48550/arXiv.2301.01574

2023-01-04

Abstract:Multi-wave inverse problems are indirect imaging methods using the interaction of two different imaging modalities. One brings spatial accuracy, and the other contrast sensitivity. The inversion method typically involve two steps. The first step is devoted to accessing internal datum of quantities related to the unknown parameters being observed. The second step involves recovering the parameters themselves from the internal data. To perform that inversion, a typical requirement is that the Jacobian of fields involved does not vanish. A number of authors have considered this problem in the past two decades, and a variety of methods have been developed. Existing techniques require H{ö}lder continuity of the parameters to be reconstructed. In practical applications, the medium may present embedded elements, with distinct physical properties, leading to discontinuous coefficients. In this article we explain how a Jacobian constraint can imposed in the piecewise regular case, when the physical model is a divergence form second order linear elliptic boundary value problem.

Analysis of PDEs

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to impose a non - zero Jacobian determinant constraint in elliptic boundary value problems with piecewise regular coefficients in multi - wave inverse problems. Specifically, the paper focuses on how to ensure the stability and accuracy in the parameter reconstruction process through appropriate mathematical methods and techniques when there may be embedded elements in the medium, resulting in discontinuous coefficients. This problem is very important in practical applications because many physical systems, such as geological strata or human organs, may have different physical properties in different regions, resulting in discontinuities in the parameter field. The method proposed in the paper is not only applicable to two - dimensional cases, but also extended to higher dimensions, and is of great significance for dealing with complex multi - physical coupling imaging problems.