A Nonlocal $p$-Laplacian Interface Model with Sharp Interface

Resumo Técnico: Um Modelo de Interface $p$-Laplaciano Não Local com Interface Nítida

Enunciado do Problema
O artigo aborda o desafio de formular problemas de interface dentro de um arcabouço não local. Problemas de interface clássicos envolvem equações diferenciais parciais (EDPs) definidas em um domínio $\Omega$ dividido por uma interface $\Gamma$, onde a solução ou suas derivadas exibem descontinuidades. Especificamente, os autores consideram a extensão do problema de interface de Poisson e do problema de interface $p$-laplaciano com condições de membrana para configurações não locais. Uma dificuldade primária em formulações não locais é que as condições de interface locais clássicas (como a continuidade da solução e saltos prescritos no fluxo normal) não podem ser aplicadas diretamente. As abordagens existentes frequentemente aproximam a interface nítida $\Gamma$ com uma camada volumétrica de espessura $\delta$, o que pode complicar a modelagem da difusão independente nos subdomínios e o tratamento de dados de contorno.

Metodologia
Os autores propõem uma nova aproximação não local baseada em energia que mantém uma interface nítida $\Gamma$. Em vez de modificar a geometria do domínio, o modelo utiliza funcionais de energia específicos onde as condições de interface são codificadas implicitamente.

1. Formulação de Energia:

   • Para o caso linear ($p=2$), o funcional de energia local $E(u)$ é aproximado por um funcional não local $E_\delta(u)$. O termo de difusão não local é construído usando um núcleo $\Lambda(x,y)$ que restringe as interações para dentro dos subdomínios $\Omega_1$ e $\Omega_2$, garantindo a difusão não local independente em cada região.
   • A condição de interface de Dirichlet (salto no fluxo) é tratada via um termo de penalidade envolvendo médias não locais ($K_{\delta,i} * u$) ao longo da interface $\Gamma$.
   • A condição de contorno de Neumann é incorporada naturalmente através da estrutura variacional.
   • O arcabouço é estendido para o caso não linear do $p$-laplaciano com condições de membrana, onde o salto da solução através da interface é penalizado. O funcional não local resultante $F_\delta(u)$ compartilha a mesma forma estrutural de $E_\delta(u)$.

2. Análise Teórica:

   • Existência e Unicidade: Os autores estabelecem a existência e a unicidade de minimizadores para os funcionais não locais $E_\delta$ e $F_\delta$ usando o método direto do cálculo de variações. Isso se baseia em uma desigualdade de Poincaré não local recém-derivada que considera a penalidade de salto da interface.
   • $\Gamma$-Convergência: A principal contribuição teórica é a prova de $\Gamma$-convergência. Os autores demonstram que, à medida que o horizonte não local $\delta$ desaparece ($\delta \to 0$), os funcionais não locais $\Gamma$-convergem para seus contrapartes locais ($E$ e $F$).
   • Compacidade: Combinada com resultados de compacidade para operadores $p$-laplacianos não locais, a $\Gamma$-convergência implica que os minimizadores dos problemas não locais convergem para os minimizadores dos problemas de interface locais na norma $L^p$.

3. Implementação Numérica:

   • Os autores empregam um método de elementos finitos eficiente baseado na estrutura de produto tensorial de núcleos gaussianos e bases polinomiais multicúbicas.
   • Esta abordagem decompõe integrais de alta dimensão em produtos de integrais unidimensionais, reduzindo significamente o custo computacional.
   • Regras de quadratura especiais (Gauss-Legendre) são adaptadas para lidar com as camadas de contorno resultantes das funções de base cúbica próximas às interfaces das células.

Principais Resultados

• Preservação da Interface Nítida: Ao contrário das abordagens anteriores de restrição volumétrica, o modelo proposto mantém uma interface nítida, permitindo a difusão não local independente em subdomínios disjuntos sem estender os dados da interface para o interior (bulk).
• Prova de Convergência: Provas rigorosas confirmam que os minimizadores dos funcionais não locais convergem para os minimizadores locais quando $\delta \to 0$. A taxa de convergência observada é de primeira ordem.
• Validação Numérica: Experimentos numéricos em uma e duas dimensões verificam a convergência teórica. O erro $L^2$ entre a solução não local e a solução local de referência decai linearmente com $\delta$, e os resultados são robustos em diferentes resoluções de malha.

Significância e Alegações
O artigo alega que o modelo de interface nítida proposto oferece um arcabouço mais natural e flexível para problemas de interface em comparação com as aproximações volumétricas existentes. Ao formular o problema via funcionais de energia, as condições de interface de fluxo normal, difíceis de manipular, são satisfeitas implicitamente. Os autores destacam que esta abordagem:

1. Facilita a extensão para problemas de interface mais complexos, como demonstrado pela aplicação bem-sucedida ao $p$-laplaciano não linear com condições de membrana.
2. Evita a necessidade de estender os dados da interface (por exemplo, $g$) de $L^2(\Gamma)$ para $H^1(\Omega)$.
3. Fornece um elo matematicamente rigoroso entre modelos de interface não locais e locais através da $\Gamma$-convergência, garantindo que o modelo não local seja uma aproximação válida do problema físico clássico conforme o horizonte desaparece.

O trabalho não propõe novas aplicações físicas, mas sim fornece uma ferramenta matemática e numérica robusta para modelar fenômenos de interface existentes em campos como dinâmica de fluidos e ciência dos materiais, onde efeitos não locais são relevantes.