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

1. Il Problema: Il Problema di Calderón Linearizzato in 3D

Il lavoro si concentra sul problema di Calderón linearizzato in tre dimensioni spaziali, definito su una palla unitaria $B \subset \mathbb{R}^3$.

• Contesto: Si considera un problema di conducibilità elettrica dove il coefficiente di conducibilità $\gamma$ è una perturbazione $\eta$ attorno a un valore di riferimento costante (assunto $\gamma \equiv 1$).
• Modello Fisico: La potenziale elettrica interna $u$ soddisfa l'equazione di conducibilità $-\nabla \cdot (\gamma \nabla u) = 0$ con condizioni al contorno di Neumann (densità di corrente superficiale $f$).
• Mappa Dati: La mappa di interesse è l'operatore di Neumann-to-Dirichlet (ND) $\Lambda(\gamma)$, che mappa la corrente applicata $f$ al potenziale misurato sulla superficie $\partial B$.
• Linearizzazione: Si studia la derivata di Fréchet $F = D\Lambda(1; \eta)$ dell'operatore ND rispetto alla perturbazione $\eta$. L'obiettivo è ricostruire $\eta$ conoscendo l'operatore lineare $F\eta$.
• Sfida Principale: A differenza dei lavori precedenti che trattavano perturbazioni limitate ($L^\infty$ o $L^2$), questo studio affronta la ricostruzione di perturbazioni illimitate, specificamente nello spazio $L^3(B)$. È dimostrato che l'operatore $F$ si estende continuamente a perturbazioni in $L^3(B)$, permettendo di gestire singolarità più forti rispetto al caso bidimensionale classico ($L^2$).

2. Metodologia: La Base di Zernike 3D e Sostituzione Diretta

Il cuore della metodologia risiede nell'uso di una base di Zernike tridimensionale per espandere la perturbazione incognita $\eta$.

• Espansione della Perturbazione: La funzione $\eta \in L^3(B)$ viene espansa in una serie di funzioni di base ortonormali $\psi_{\ell}^{k,m}$, definite come il prodotto di polinomi radiali di Zernike 3D ($R_\ell^k$) e armoniche sferiche ($Y_\ell^m$):
  $$\eta = \sum_{k, \ell, m} c_{\ell}^{k,m} \psi_{\ell}^{k,m}$$
• Struttura Triangolare (Forward Substitution): Il risultato principale è una formula di sostituzione diretta (forward substitution) per calcolare i coefficienti $c_{\ell}^{k,m}$.
  • I coefficienti possono essere ricostruiti in modo induttivo partendo da $k=0$.
  • Per calcolare i coefficienti con indice radiale $k$, è necessario conoscere solo i coefficienti con indici $q < k$ (e specifici indici angolari $\ell$).
  • La formula esplicita (Teorema 1.1) lega i coefficienti incogniti ai dati misurati $(F\eta)$ proiettati su specifiche armoniche sferiche, sottraendo i contributi dei coefficienti già calcolati.
• Efficienza Computazionale: Il metodo è estremamente efficiente perché:
  1. Evita la risoluzione di sistemi lineari globali o l'uso di tecniche iterative complesse.
  2. Richiede un sottoinsieme relativamente piccolo di misurazioni al contorno. A differenza del caso 2D dove tutte le armoniche di Fourier sono necessarie, in 3D è sufficiente applicare correnti corrispondenti alle armoniche sferiche $Y_{k+1}^0$ (dove l'indice azimutale $m$ è zero) per ricostruire esattamente la proiezione ortogonale $\eta_k$.

3. Contributi Chiave

1. Estensione a $L^3(B)$: Il paper dimostra che il problema di ricostruzione lineare è ben posto per perturbazioni nello spazio $L^3(B)$, un passo avanti rispetto alle limitazioni $L^2$ o $L^\infty$ dei metodi precedenti. Questo è cruciale per modellare perturbazioni fisiche con singolarità.
2. Formula di Ricostruzione Esatta: Viene fornito un algoritmo analitico esatto (Teorema 1.1) basato su coefficienti scalari pre-calcolabili ($Q_{\ell,s}^{k,m}$) che coinvolgono coefficienti di Gaunt e simboli 3j di Wigner.
3. Ridondanza dei Dati: Si dimostra che in 3D la ricostruzione richiede meno dati rispetto alla base completa delle correnti. L'uso esclusivo di correnti con simmetria assiale ($m=0$) è sufficiente per ricostruire l'intera struttura della perturbazione grazie alle proprietà di ortogonalità e alle relazioni di ricorrenza delle armoniche sferiche.
4. Analisi della Stabilità Numerica: Viene discusso il compromesso tra efficienza e stabilità. Sebbene il metodo sia diretto, la presenza di una somma aggiuntiva (necessaria perché il prodotto di due armoniche sferiche è una combinazione lineare di armoniche di diversi gradi, a differenza del caso 2D) rende il metodo potenzialmente meno stabile numericamente rispetto alla controparte 2D, richiedendo una regolarizzazione attenta.

4. Risultati e Validazione Numerica

Gli autori hanno validato il metodo attraverso esempi numerici:

• Caso Ideale (Dati Precisi): Utilizzando dati simulati con alta precisione (16 cifre decimali), il metodo ha ricostruito perturbazioni gaussiane localizzate con estrema accuratezza, anche con indici di troncamento ridotti. La natura triangolare dell'algoritmo garantisce che i coefficienti calcolati non siano influenzati dal troncamento degli indici superiori.
• Caso Realistico (Rumore e Discretizzazione): Simulando misurazioni ottenute da una discretizzazione agli elementi finiti (FE) "ruvida" (con rumore numerico intrinseco), il metodo ha prodotto ricostruzioni approssimate ma fisicamente significative.
• Regolarizzazione: È stato mostrato che una troncatura dipendente dall'indice $k$ degli indici angolari $\ell$ agisce come una regolarizzazione efficace, permettendo di ottenere ricostruzioni stabili anche in presenza di dati imperfetti.

5. Significato e Implicazioni

Questo lavoro rappresenta un avanzamento significativo nella teoria inversa per l'elettroimpedenza tomografica (EIT):

• Generalizzazione Dimensionale: Estende con successo le tecniche di ricostruzione diretta basate su Zernike dal piano (2D) allo spazio (3D), adattandosi alle complessità matematiche aggiuntive (come la moltiplicazione delle armoniche sferiche).
• Gestione delle Singolarità: La capacità di trattare perturbazioni in $L^3$ apre la porta alla modellazione di difetti o inclusioni con gradienti di conducibilità molto elevati o singolari, che erano precedentemente fuori portata per i metodi lineari diretti.
• Efficienza Pratica: La necessità di un set di misurazioni ridotto (solo correnti con $m=0$) suggerisce che l'implementazione pratica di sensori EIT potrebbe essere semplificata, riducendo il numero di elettrodi o configurazioni di misura necessarie per una ricostruzione di base.
• Fondamento Teorico: Fornisce un quadro teorico solido per futuri sviluppi algoritmici, dimostrando che la struttura triangolare delle equazioni di ricostruzione può essere preservata anche in dimensioni superiori e per spazi di funzioni più ampi.

In sintesi, Garde e Hirvensalo offrono un metodo diretto, esatto ed efficiente per la ricostruzione 3D di perturbazioni di conducibilità, superando i limiti di regolarità delle perturbazioni e ottimizzando il numero di misurazioni necessarie, pur gestendo le sfide di stabilità numerica intrinseche alla dimensionalità superiore.