A direct imaging method for inverse scattering problem of biharmonic wave with phased and phaseless data

Resumo Técnico: Um Método de Imageamento Direto para Problemas de Espalhamento Inverso Biharmônico

Formulação do Problema
Este artigo aborda o problema de espalhamento inverso biharmônico, que visa determinar a forma e a localização de um obstáculo estendido $D$ em $\mathbb{R}^2$ utilizando dados de medição. Ao contrário dos problemas de espalhamento inverso mais comuns associados a operadores diferenciais de segunda ordem (acústico, elástico, eletromagnético), este trabalho foca na equação de onda biharmônica, uma equação diferencial parcial de quarta ordem. O estudo considera tanto dados de fase (onde a informação de fase está disponível) quanto dados sem fase (onde apenas o módulo do campo é medido). O obstáculo está sujeito a uma de quatro condições de contorno enraizadas na teoria das placas: engastada, simplesmente apoiada, apoiada sobre roletes ou livre.

Metodologia
Os autores propõem um método de imageamento direto baseado na técnica de Migração de Tempo Reverso (RTM). Ao contrário dos métodos iterativos, esta abordagem não requer suposições geométricas prévias sobre o espalhador, nem necessita da resolução numérica do problema de espalhamento direto em cada iteração.

O cerne da metodologia envolve a construção de funções de imageamento específicas ($I_j(z)$) para os pontos de amostragem $z$ dentro de um domínio $\Omega$. Estas funções são projetadas para exibir características distintas dependendo de o ponto de amostragem estar próximo ou distante da fronteira do obstáculo.

1. Dados de Fase (IP-1 a IP-10): Os autores definem dez funções de imageamento utilizando várias combinações de dados de medição, incluindo:

   • Campo espalhado ($u^{sc}$) e sua derivada normal ($\partial_n u^{sc}$).
   • Momento fletor ($M u^{sc}$) e força transversal ($N u^{sc}$).
   • Padrões de campo distante ($u^\infty$).
   • Dados gerados por fontes pontuais e ondas planas.

   As funções de imageamento são construídas utilizando a solução fundamental da equação de Helmholtz ($\Phi_\kappa$) em vez da solução fundamental da equação de onda biharmônica. Esta formulação destina-se a reduzir a complexidade computacional. As funções envolvem integrais sobre as fronteiras de medição ($\Gamma_r, \Gamma_s$) ou o círculo unitário ($S$), combinando os dados medidos com a solução fundamental de Helmholtz e suas derivadas.

2. Dados Sem Fase (IP-11): Para o caso em que apenas o módulo do campo total $|u|$ está disponível, uma função de imageamento específica $I_{11}(z)$ é construída. Esta função baseia-se na diferença entre o quadrado do módulo do campo total e o quadrado do módulo do campo incidente, normalizado pelo campo incidente.

Análise Teórica e Contribuições Principais
O artigo fornece uma análise matemática rigorosa das funções de imageamento propostas:

• Análise de Resolução: Os autores provam que, para qualquer ponto de amostragem $z$, a função de imageamento $I_j(z)$ pode ser decomposta em um termo principal e um termo de ruído/erro. O termo principal é proporcional à integral do quadrado do módulo do padrão de campo distante de uma solução auxiliar $\psi(\cdot, z)$, que satisfaz a equação de onda biharmônica com dados de fronteira derivados da parte imaginária da solução fundamental de Helmholtz.
• Mecanismo de Contraste: A análise demonstra que as funções de imageamento exibem um alto contraste quando o ponto de amostragem $z$ está próximo da fronteira do obstáculo e decaem à medida que $z$ se afasta. Este comportamento é atribuído às propriedades da parte imaginária da solução fundamental da equação de Helmholtz (uma função de Bessel), que atinge o pico quando os pontos de fonte e observação coincidem.
• Novidade na Formulação: Uma contribuição significativa é a derivação de funções de imageamento que dependem apenas da solução fundamental da equação de Helmholtz (ou de uma onda plana acústica) em vez da solução fundamental biharmônica. Os autores observam que estudos anteriores de RTM para ondas acústicas, elásticas e eletromagnéticas não podem ser estendidos diretamente ao caso biharmônico devido à ordem superior do operador diferencial.
• Robustez às Condições de Contorno: O arcabouço teórico é estabelecido como aplicável a todos os quatro tipos de condições de contorno (engastada, simplesmente apoiada, apoiada sobre roletes e livre), com comportamentos assintóticos específicos derivados para cada uma.

Resultados Numéricos
O artigo apresenta experimentos numéricos para validar as descobertas teóricas:

• Reconstrução de Obstáculo Único: O algoritmo reconstrui com sucesso um círculo unitário utilizando todas as 11 funções de imageamento (10 de fase, 1 sem fase) sem ruído.
• Robustez ao Ruído: Experimentos com um obstáculo em forma de pipa demonstram a robustez do algoritmo contra níveis de ruído de 5% e 10%. As formas reconstruídas permanecem precisas e distintas.
• Múltiplos Obstáculos: O método é mostrado como eficaz para localizar e dar forma a dois obstáculos distintos (um círculo e uma pipa) simultaneamente, mesmo na presença de ruído.

Significado e Alegações
O artigo alega que o método de imageamento direto proposto oferece duas vantagens principais sobre os métodos iterativos:

1. Elimina a necessidade de suposições geométricas prévias sobre o espalhador.
2. Evita o ônus computacional de resolver repetidamente o problema de espalhamento direto.

Os autores enfatizam que sua abordagem é a primeira a aplicar imageamento direto baseado em RTM a problemas de espalhamento inverso biharmônico usando um único tipo de dados de medição (variando de campos espalhados a campos totais sem fase), enquanto confia na solução fundamental mais simples de Helmholtz. A análise de resolução confirma que o método fornece um contraste claro para identificar a fronteira do obstáculo, e os experimentos numéricos verificam sua eficiência na reconstrução de geometrias complexas e sua estabilidade em ambientes ruidosos. O trabalho estende a aplicabilidade dos métodos de imageamento direto de equações de onda de segunda ordem para equações de onda biharmônicas de quarta ordem.