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

1. Probleemstelling

Het artikel behandelt het lineaire Calderón-probleem in drie dimensies. Dit is een inverse probleem waarbij men een onbekende geleidbaarheidsverstorende term (perturbatie) $\eta$ in een eenheidsbol $B \subset \mathbb{R}^3$ wil reconstrueren op basis van lineaire data.

• Context: Het klassieke Calderón-probleem probeert de geleidbaarheid $\gamma$ te reconstrueren uit de Neumann-naar-Dirichlet (ND) kaart $\Lambda(\gamma)$. In de lineaire benadering wordt $\gamma$ geschreven als $1 + \eta$, waarbij $\eta$ een kleine verstoring is.
• De Operator: De lineaire operator $F = D\Lambda(1; \eta)$ beschrijft de relatie tussen de verstoring $\eta$ en de verandering in de ND-kaart. Voor harmonische functies $u_f$ en $u_g$ met Neumann-randen $f$ en $g$ geldt:
  $$\langle (F\eta)f, g \rangle_{L^2(\partial B)} = -\int_B \eta \nabla u_f \cdot \nabla u_g \, dx$$
• Unieke Uitdaging: Het doel is om $\eta$ te reconstrueren uit $F\eta$. Een belangrijk aspect van dit werk is dat de methode geldig is voor verstoringen in de ruimte $L^3(B)$, wat ongebonden (unbounded) functies toelaat, in tegenstelling tot eerdere methoden die vaak beperkt waren tot $L^2$ of $L^\infty$.

2. Methodologie

De auteurs ontwikkelen een exacte directe reconstructiemethode gebaseerd op een 3D Zernike-basis. De aanpak is een uitbreiding van een eerder 2D-algoritme (Garde & Hyvönen, 2024) naar drie dimensies.

• Basisfuncties: De verstoring $\eta$ wordt ontwikkeld in een orthonormale basis bestaande uit producten van 3D radiale Zernike-polynomen $R^k_\ell(r)$ en sferische harmonischen $Y^m_\ell(\theta, \phi)$:
  $$\psi^{m}_{k,\ell}(r, \theta, \phi) = R^k_\ell(r) Y^m_\ell(\theta, \phi)$$
  Hierbij is $k$ gerelateerd aan het radiale gedrag en $\ell, m$ aan de hoekvariabelen.
• Forward Substitution (Voorwaartse Substitutie): De kern van de methode is dat de reconstructie van de coëfficiënten $c^{m}_{k,\ell}$ van de verstoring $\eta$ een triangulaire structuur heeft.
  • De coëfficiënten met een lage index $k$ kunnen direct worden berekend uit de meetdata.
  • Coëfficiënten met een hogere $k$ worden berekend door de reeds gevonden coëfficiënten met lagere $k$ te gebruiken (iteratief).
• Meetdata: Een opvallend kenmerk is dat voor de exacte reconstructie slechts een relatief klein subset van randmetingen nodig is. Specifiek worden alleen Neumann-randvoorwaarden gebruikt die corresponderen met sferische harmonischen van de vorm $Y^0_{k+1}$ (d.w.z. met $m=0$). Dit is in contrast met een volledige basis van stroomdichtheden; in 3D is er minder redundantie nodig dan in 2D.

3. Belangrijkste Bijdragen en Resultaten

Hoofdstelling (Theorem 1.1)

De auteurs bewijzen dat voor elke $\eta \in L^3(B)$, de coëfficiënten $c^{m}_{k,\ell}$ in de Zernike-ontwikkeling exact kunnen worden gereconstrueerd via een expliciete formule:
$$c^{m}_{k,\ell} = (Q^{m,k}_{\ell,0})^{-1} \left( \langle (F\eta)Y^0_{k+1}, Y^m_{\ell+k+1} \rangle_{L^2(\partial B)} - \sum_{q=0}^{k-1} \sum_{s=0}^{k-q} Q^{m,q}_{\ell,s} c^{m}_{\ell+2s} \right)$$

• De termen $Q$ zijn scalair en onafhankelijk van $\eta$, gedefinieerd via Gaunt-coëfficiënten (integralen van producten van drie sferische harmonischen) en Pochhammer-symbolen.
• De methode werkt inductief: eerst wordt $\eta_0$ (harmonische componenten) gereconstrueerd, dan $\eta_1$, enzovoort.

Uitbreiding naar $L^3$

Het artikel toont aan dat de lineaire operator $F$ continu uitbreidbaar is naar perturbaties in $L^3(B)$ (voor $d=3$). Dit is significant omdat het toelaat om ongebonden verstoringen te reconstrueren, wat in de praktijk relevant is voor materialen met sterke lokale variaties.

Numerieke Stabiliteit en Efficiëntie

• Efficiëntie: Omdat het een "forward substitution" is, is de methode computatieel zeer efficiënt en eenvoudig te implementeren.
• Data-efficiëntie: In 3D zijn minder metingen nodig dan in 2D. In 2D moeten alle Fourier-basisfuncties worden gebruikt, terwijl in 3D alleen de $m=0$ componenten van de sferische harmonischen nodig zijn voor de invoer.
• Stabiliteit: De methode is numeriek stabiel, hoewel de 3D-versie iets complexer is dan de 2D-versie vanwege extra sommaties die ontstaan door het product van twee sferische harmonischen (wat leidt tot een eindige lineaire combinatie van harmonischen van verschillende graden).

4. Numerieke Voorbeelden

De auteurs demonstreren de methode met een getransleerde, zeer gelokaliseerde Gaussische verstoring:

• Perfecte Data: Met nauwkeurige, gesimuleerde data (16 cijfers nauwkeurig) levert de methode bijna perfecte reconstructies op, zelfs met afgeknotte indices.
• Ruwe Data: Zelfs bij gebruik van ruwe data gegenereerd door een eindige-elementenmethode (FEM) met ruis, produceert de methode redelijke reconstructies. De triangulaire structuur zorgt ervoor dat fouten in hogere indices de berekening van lagere indices niet beïnvloeden (in tegenstelling tot veel andere inverse methoden).

5. Significatie

Dit werk is een belangrijke stap in de theoretische en numerieke oplossing van het lineaire Calderón-probleem in drie dimensies:

1. Generalisatie: Het generaliseert succesvolle 2D-methoden naar 3D met behoud van de exacte directe reconstructie-eigenschap.
2. Ruimtelijke Robuustheid: Het toont aan dat reconstructie mogelijk is voor ongebonden verstoringen ($L^3$), wat de toepasbaarheid vergroot voor realistische, onregelmatige materialen.
3. Praktische Haalbaarheid: De methode vereist minder meetdata dan theoretisch mogelijk zou zijn (gebruikmakend van symmetrie) en is computatieel goedkoop door de triangulaire structuur.
4. Toepassing: De resultaten zijn relevant voor elektrische impedantie-tomografie (EIT) en andere toepassingen waar inwendige geleidbaarheidsvariaties moeten worden gekarakteriseerd aan de hand van oppervlaktmetingen.

Samenvattend biedt dit artikel een wiskundig onderbouwde, efficiënte en robuuste algoritme voor het direct reconstrueren van complexe 3D-geleidbaarheidsprofielen uit lineaire randmetingen.