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

Resumo Técnico: Reconstrução de Perturbações Ilimitadas no Problema de Calderón Linearizado em 3D

1. O Problema

O artigo aborda o Problema de Calderón Linearizado em três dimensões espaciais, especificamente dentro de uma bola unitária $B \subset \mathbb{R}^3$.

• Contexto: O objetivo é reconstruir uma perturbação de condutividade $\eta$ (onde a condutividade total é $\gamma = 1 + \eta$) a partir de dados de fronteira linearizados.
• Dados: O dado disponível é o operador derivada de Fréchet $F = D\Lambda(1; \eta)$ do mapa Neumann-para-Dirichlet (ND) avaliado em $\gamma \equiv 1$. A relação é dada por:
  $$\langle (F\eta)f, g \rangle_{L^2(\partial B)} = -\int_B \eta \nabla u_f \cdot \nabla u_g \, dx$$
  onde $u_f$ e $u_g$ são funções harmônicas com traços de Neumann $f$ e $g$.
• Inovação Principal do Escopo: Diferente de trabalhos anteriores que restringiam $\eta$ ao espaço $L^2(B)$, este trabalho demonstra que é possível realizar a reconstrução exata para perturbações no espaço $L^3(B)$. Isso permite lidar com perturbações ilimitadas (unbounded perturbations), uma vez que $L^3(B) \supset L^\infty(B)$.

2. Metodologia

A abordagem proposta baseia-se em uma expansão em série utilizando uma base de Zernike 3D, adaptada para o domínio esférico.

• Base de Expansão:
  A perturbação $\eta$ é expandida na base ortonormal $\{\psi_{\ell}^{k,m}\}$, definida como o produto de polinômios radiais de Zernike 3D ($R_\ell^k$) e harmônicos esféricos ($Y_\ell^m$):
  $$\psi_{\ell}^{k,m}(r, \theta, \phi) = R_\ell^k(r) Y_\ell^m(\theta, \phi)$$
  onde $k, \ell \in \mathbb{N}_0$ e $m \in \{-\ell, \dots, \ell\}$.

• Estrutura de Substituição Direta (Forward Substitution):
  O método central é um algoritmo de substituição direta (triangular). Ao invés de resolver um sistema global inverso, os coeficientes da expansão de $\eta$ são calculados sequencialmente:

  1. Os coeficientes correspondentes a $k=0$ (subespaço de funções harmônicas) são recuperados diretamente de uma medição específica.
  2. Os coeficientes para $k=1$ são recuperados usando os dados de medição e os coeficientes já calculados de $k=0$.
  3. Geralmente, os coeficientes para um índice $k$ são recuperados usando os dados de fronteira e todos os coeficientes de índices menores ($0, \dots, k-1$).

• Redução de Dados de Fronteira:
  Uma característica notável é a eficiência na coleta de dados. Para reconstruir o termo de ordem $k$, o método requer apenas medições com densidades de corrente de Neumann específicas da forma $Y_{k+1}^0$ (harmônicos esféricos com ordem $m=0$).

  • Isso contrasta com o caso 2D, onde todas as bases de Fourier são necessárias.
  • Em 3D, não é necessário variar o índice $m$ nas medições, reduzindo significativamente a redundância de dados (já que o número de modos $m$ cresce como $2\ell+1$).

• Fórmula de Reconstrução (Teorema 1.1):
  Os coeficientes $c_{\ell}^{k,m}$ são dados explicitamente por:
  $$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 - \sum_{q=0}^{k-1} \sum_{s=0}^{k-q} Q_{\ell,s}^{k,m,q} c_{\ell+2s}^{q,m} \right)$$
  Onde os escalares $Q$ dependem de coeficientes de Gaunt e símbolos 3j de Wigner, mas são independentes de $\eta$.

3. Contribuições Chave

1. Extensão para $L^3(B)$: O trabalho prova que o mapa linearizado $F$ estende-se continuamente para perturbações em $L^3(B)$, permitindo a reconstrução de condutividades que podem ser ilimitadas (desde que em $L^3$).
2. Algoritmo Exato e Direto: Fornece uma fórmula analítica fechada para a reconstrução exata de qualquer $\eta \in L^3(B)$, sem necessidade de iteração ou otimização não linear.
3. Eficiência Computacional: O método é computacionalmente eficiente devido à sua natureza de substituição direta e à necessidade de um subconjunto relativamente pequeno de medições de fronteira (apenas harmônicos com $m=0$).
4. Generalização 3D: Estende o método de reconstrução baseado em Zernike 2D (anteriormente desenvolvido por Garde & Hyvönen) para três dimensões, lidando com a complexidade adicional da multiplicação de harmônicos esféricos (que gera combinações lineares de vários graus, ao contrário da ortogonalidade simples em Fourier 2D).

4. Resultados e Validação Numérica

• Exemplo Numérico: Os autores testaram o método com uma perturbação gaussiana localizada e deslocada.
  • Cenário Ideal: Com dados simulados de alta precisão (16 dígitos), a reconstrução aproximada seguiu o padrão exato da perturbação à medida que o número de termos $K$ na soma aumentava.
  • Cenário Realista: O método foi testado com dados gerados por um modelo de Elementos Finitos (FE) grosseiro e com ruído. Mesmo com dados imprecisos, a reconstrução dos primeiros termos da série ($k=0$ a $4$) foi "quase perfeita" e numericamente estável.
• Estabilidade: Embora o método em 3D seja potencialmente menos estável numericamente do que em 2D (devido à necessidade de somar coeficientes com índices $\ell$ maiores do que o atual), os experimentos mostraram que ele produz aproximações robustas mesmo com dados ruidosos.

5. Significado e Impacto

• Avanço Teórico: A prova de que a reconstrução exata é possível para $L^3$ em 3D expande os limites teóricos do problema inverso de condutividade, permitindo modelar materiais com descontinuidades mais severas.
• Aplicabilidade Prática: A redução na quantidade de medições de fronteira necessárias (apenas modos $m=0$) torna o método mais viável para implementações experimentais onde a aquisição de dados completos de todos os modos de harmônicos esféricos seria proibitivamente cara ou complexa.
• Ferramenta Numérica: O método oferece uma alternativa direta e rápida aos métodos iterativos tradicionais (como Newton-Raphson) para problemas de tomografia de impedância elétrica (EIT), especialmente útil como inicializador para métodos não lineares ou para aplicações onde a linearização é uma aproximação válida.

Em resumo, o artigo estabelece um marco na teoria inversa 3D ao fornecer um método de reconstrução direta, eficiente e matematicamente rigoroso para perturbações de condutividade em espaços funcionais mais amplos, superando limitações anteriores de $L^2$.