Linearised Calderón problem: Reconstruction of unbounded perturbations in 3D

1. Problématique et Contexte

L'article s'intéresse au problème de Calderón linéarisé en trois dimensions. Ce problème inverse consiste à reconstruire une perturbation de conductivité $\eta$ à l'intérieur d'un domaine (ici la boule unité $B \subset \mathbb{R}^3$) à partir de mesures de données de frontière.

• Modèle physique : Le problème est basé sur l'équation de conductivité $-\nabla \cdot (\gamma \nabla u) = 0$ avec des conditions de Neumann sur la frontière. Le paramètre de conductivité est $\gamma = 1 + \eta$.
• Linéarisation : Au lieu de reconstruire la conductivité complète $\gamma$ (problème non linéaire), les auteurs travaillent sur la dérivée de Fréchet de l'application Neumann-Dirichlet (ND) $\Lambda$ autour de la conductivité de fond $\gamma \equiv 1$. L'opérateur linéarisé $F = D\Lambda(1; \eta)$ relie la perturbation $\eta$ aux mesures de frontière.
• Défi majeur : L'objectif est de reconstruire des perturbations non bornées, c'est-à-dire des fonctions $\eta$ appartenant à l'espace $L^3(B)$ (et non seulement $L^\infty(B)$). Cela est rendu possible par le fait que l'opérateur linéarisé s'étend continûment aux perturbations dans $L^d(\Omega)$ pour un domaine de dimension $d$.

2. Méthodologie

La méthode proposée est une extension tridimensionnelle d'un algorithme précédent développé pour le cas bidimensionnel (2D). Elle repose sur les éléments suivants :

• Base de Zernike 3D : La perturbation $\eta$ est développée dans une base orthonormée de $L^2(B)$ constituée de fonctions de Zernike tridimensionnelles :
  $$\psi_{\ell}^{k,m}(r, \theta, \phi) = R_{\ell}^k(r) Y_{\ell}^m(\theta, \phi)$$
  où $Y_{\ell}^m$ sont les harmoniques sphériques et $R_{\ell}^k$ sont les polynômes radiaux de Zernike en 3D.
• Structure triangulaire (Substitution directe) : L'algorithme exploite une structure triangulaire dans les équations reliant les coefficients de la perturbation aux données de mesure. Cela permet une reconstruction exacte et directe par substitution directe (forward substitution), évitant ainsi la nécessité de résoudre un système linéaire global ou d'utiliser des méthodes itératives coûteuses.
• Données de mesure minimales : Contrairement à une base complète de densités de courant, la méthode nécessite uniquement des mesures correspondant à des harmoniques sphériques spécifiques de la forme $Y_{k+1}^0$ (indépendantes de l'indice $m$). Cela réduit considérablement la redondance des données par rapport au cas 2D où toutes les bases de Fourier sont utilisées.
• Traitement des coefficients : Les coefficients $c_{\ell}^{k,m}$ de la perturbation sont reconstruits de manière inductive :
  1. Les coefficients de rang $k=0$ (fonctions harmoniques) sont reconstruits directement à partir des données.
  2. Les coefficients de rang $k$ sont reconstruits en utilisant les données de rang $k$ et les coefficients déjà reconstruits pour les rangs $0, \dots, k-1$.

3. Résultats Principaux

Le résultat central est le Théorème 1.1, qui fournit une formule explicite pour le calcul des coefficients $c_{\ell}^{k,m}$ :

$$c_{\ell}^{k,m} = (Q_{\ell,0}^{k,m,k})^{-1} \left( \langle (F\eta)Y_{k+1}^0, Y_{\ell+k+1}^m \rangle_{L^2(\partial B)} - \sum_{q=0}^{k-1} \sum_{s=0}^{k-q} Q_{\ell,s}^{k,m,q} c_{\ell+2s}^{q,m} \right)$$

• Coefficients $Q$ : Ces scalaires sont indépendants de la perturbation $\eta$ et dépendent uniquement des indices $(k, \ell, m)$ et de coefficients de Gaunt (intégrales de produits de trois harmoniques sphériques).
• Validité pour $L^3$ : Le théorème garantit la reconstruction exacte pour toute perturbation $\eta \in L^3(B)$, élargissant ainsi le cadre d'application par rapport aux méthodes précédentes limitées à $L^2$ ou $L^\infty$.
• Stabilité numérique : Bien que la méthode en 3D soit potentiellement moins stable numériquement que son équivalent 2D (en raison d'une somme supplémentaire impliquant des coefficients avec des indices $\ell$ plus élevés), les simulations montrent que la méthode reste robuste.
• Exemples numériques : Les auteurs présentent des simulations montrant que :
  • Avec des données exactes, la reconstruction est quasi parfaite.
  • Avec des données bruitées (simulées via une méthode des éléments finis grossière), la méthode permet de reconstruire les projections orthogonales de la perturbation avec une bonne précision, même avec un troncature des indices.

4. Contributions Clés

1. Extension à la dimension 3 : Passage d'une méthode de reconstruction 2D à une méthode 3D utilisant une base de Zernike adaptée.
2. Reconstruction de perturbations non bornées : Preuve et algorithme permettant de traiter des perturbations dans $L^3(B)$, un espace plus large que $L^\infty(B)$, ce qui est crucial pour modéliser des défauts réalistes ou des singularités.
3. Efficacité algorithmique : La nature de "substitution directe" rend l'algorithme très efficace pour l'implémentation numérique, évitant l'inversion de grandes matrices.
4. Réduction des données : Démonstration qu'il n'est pas nécessaire d'utiliser toutes les harmoniques sphériques comme courants d'entrée pour une reconstruction exacte, contrairement au cas 2D.

5. Signification et Impact

Cet article représente une avancée significative dans le domaine de l'imagerie par tomographie de conductivité électrique (ECT) et du problème inverse de Calderón.

• Théorique : Il établit des résultats d'unicité et de reconstruction constructive pour des perturbations moins régulières (dans $L^3$), comblant un vide théorique entre les résultats classiques et les cas réalistes de singularités.
• Pratique : La méthode proposée est directement implémentable et efficace. Sa capacité à fonctionner avec des données partielles (sous-ensemble de mesures) et sa robustesse face au bruit la rendent prometteuse pour des applications réelles où les mesures complètes sont coûteuses ou impossibles à obtenir.
• Limitations et Perspectives : Les auteurs notent que la stabilité numérique en 3D est légèrement inférieure à celle en 2D en raison de la structure des sommes, mais que des techniques de régularisation (comme la troncature des indices $\ell$) permettent d'obtenir de bons résultats pratiques.

En résumé, Garde et Hirvensalo ont réussi à généraliser une méthode de reconstruction exacte et rapide au cas 3D, en élargissant simultanément la classe de fonctions admissibles pour la perturbation, offrant ainsi un outil puissant pour l'analyse inverse en trois dimensions.