Team apics

Overall Objectives
Scientific Foundations
Application Domains
New Results
Contracts and Grants with Industry
Other Grants and Activities

Section: Application Domains

Geometric inverse problems for elliptic partial differential equations

Participants : Laurent Baratchart, Yannick Fischer, José Grimm, Juliette Leblond, Ana-Maria Nicu, Jonathan R. Partington, Stéphane Rigat [ Univ. Aix-Marseille I ] , Emmanuel Russ [ Univ. Aix-Marseille III ] , Edward Saff, Meriem Zghal.

We are mainly concerned with classical inverse problems like the one of localizing defaults (as cracks, pointwise sources or occlusions) in a two or three dimensional domain from boundary data (which may correspond to thermal, electrical, or magnetic measurements), of a solution to Laplace or to some conductivity equation in the domain. These defaults can be expressed as a lack of analyticity of the solution of the associated Dirichlet-Neumann problem that may be approached, in balls, using techniques of best rational or meromorphic approximation on the boundary of the object (see section 3.1 ).

Indeed, it turns out that traces of the boundary data on 2-D cross sections (disks) coincide with analytic functions in the slicing plane, that has branched singularities inside the disk [5] . These singularities are related to the actual location of the sources (namely, they reach in turn a maximum in modulus when the plane contains one of the sources). Hence, we are back to the 2-D framework where approximately recovering these singularities can be performed using best rational approximation.

In this connection, the realistic case where data are available on part of the boundary only offers a typical opportunity to apply the analytic extension techniques (see section 3.1.1 ) to Cauchy type issues, a somewhat different kind of inverse problems in which the team is strongly interested.

The approach proposed here consists in recovering, from measured data on a subset K of the boundary $ \partial$D of a domain D of R2 or R3 , say the values FK on K of some function F , the subset Im9 ${\#947 \#8834 D}$ of its singularities (typically, a crack or a discrete set of pointwise sources), provided that F is an analytic function in Im10 ${D\#8726 \#947 }$ .

The two above steps are shown in [20] to provide a robust way of locating sources from incomplete boundary data in a 2-D situation with several annular layers. Numerical experiments have already yielded excellent results in 3-D situations and we are now on the way to process real experimental magneto-encephalographic data, see also sections 5.7 , 6.3.2 . The PhD theses of A.-M. Nicu and M. Zghal are concerned with these applications, in collaboration with the Odyssée team of Inria Sophia Antipolis, and with neuroscience teams in partner-hospitals (hosp. Timone, Marseille).

Such methods are currently being generalized to problems with variable conductivity governed by a 2-D Beltrami equation, see [46] , [64] , [65] . The application we have in mind is to plasma confinement for thermonuclear fusion in a Tokamak, more precisely with the extrapolation of magnetic data on the boundary of the chamber from the outer boundary of the plasma, which is a level curve for the poloidal flux solving the original div-grad equation. Solving this inverse problem of Bernoulli type is of importance to determine the appropriate boundary conditions to be applied to the chamber in order to shape the plasma [53] . These issues are the topics of the PhD theses of S. Chaabi and Y. Fischer, and of a joint collaboration with the CEA-IRFM (Cadarache), the Laboratoire J.-A. Dieudonné at the Univ. of Nice-SA, and the CMI-LATP at the Univ. of Marseille I (see section 6.3.3 ), see [64] , [65] .

Inverse potential problems are also naturally encountered in magnetization issues that arise in nondestructive control. A particular application, which the object of a joint NSF-supported project with Vanderbilt University and MIT, is to geophysics where the remanent magnetization a rock is to be analyzed using a squid-magnetometer in order to analyze the history of the object; specifically, the analysis of Martian rocks is conducted at MIT, for instance to understand if inversions of the magnetic field took place there. Mathematically speaking, the problem is to recover the (3-D valued) magnetization m from measurements of the vector potential:

Im11 ${\#8747 _\#937 \mfrac {\mtext div~m{(x^')}~dx^'}{{|x-}x^'{|}}~,}$

outside the volume $ \upper_omega$ of the object.


Logo Inria