Singularities in phase separation models: a spectral element approach for the nonlocal Cahn-Hilliard equation

Título: Singularidades em modelos de separação de fases: uma abordagem de elementos espectrais para a equação de Cahn–Hilliard não local

1. Problema e Contexto

O artigo aborda a aproximação numérica da equação de Cahn–Hilliard não local, um modelo fundamental para descrever a separação de fases em misturas binárias (como ligas metálicas ou misturas de fluidos). Diferente do modelo clássico local, a versão não local incorpora interações de longo alcance através de um núcleo de convolução singular (especificamente o potencial newtoniano em 2D, que possui uma singularidade logarítmica na origem).

O sistema estudado inclui:

• Mobilidade constante.
• Potencial logarítmico (singular nas fases puras $\pm 1$), que garante que a variável de fase permaneça fisicamente dentro do intervalo $[-1, 1]$.
• Condições de contorno de fluxo nulo.

Desafios Principais:

1. Singularidade do Núcleo: A presença de uma singularidade no núcleo de convolução $K$ (no caso, $K(x) \sim \log|x|$) torna a discretização direta difícil, pois operadores de convolução tornam-se densos e exigem tratamento especial para evitar instabilidades numéricas.
2. Potencial Singular: O potencial logarítmico $F(\rho)$ tende a $\pm \infty$ quando $\rho \to \pm 1$, exigindo que os esquemas numéricos respeitem estritamente a separação das fases puras.
3. Custo Computacional: A natureza não local gera operadores densos, o que pode ser computacionalmente proibitivo em malhas de alta resolução.

2. Metodologia

Os autores desenvolvem um esquema numérico eficiente e controlado por erro baseado em uma combinação de análise matemática rigorosa e métodos espectrais.

A. Fundamentação Teórica

• Estendem a teoria de existência e unicidade de soluções fracas para domínios limitados com fronteira Lipschitz (incluindo domínios poligonais com cantos), que são comuns em malhas computacionais.
• Provam a propriedade de separação estrita instantânea: para $t > 0$, a solução $\rho(x,t)$ permanece uniformemente afastada das fases puras $\pm 1$ (ou seja, $|\rho| \le 1 - \delta$).
• Estabelecem a convergência das soluções para um equilíbrio único com uma taxa de decaimento polinomial, utilizando desigualdades de Łojasiewicz–Simon.

B. Método Numérico: Elementos Espectrais com "Multishape"

Para lidar com a convolução singular $K \star \rho$, os autores propõem uma técnica de decomposição de corte (cutoff decomposition) combinada com o método de elementos espectrais (SEM):

1. Decomposição da Integral: A integral de convolução é dividida em duas partes:

   • Parte Regular ($I_1$): Integração sobre o domínio excluindo uma pequena vizinhança $\epsilon$ ao redor do ponto de avaliação. Esta parte é livre de singularidades.
   • Parte Singular ($I_2$): Integração sobre a pequena vizinhança $\epsilon$. Aqui, utiliza-se uma aproximação de Taylor de ordem zero (assumindo que $\rho$ é constante localmente).
   • O termo singular é calculado analiticamente como $\rho(x) \int_{\text{vizinhança}} K(u) du$, onde a integral do núcleo sobre a geometria local é pré-calculada em forma fechada.

2. Abordagem "Multishape" (MultiShape):

   • Para calcular a parte regular $I_1$ em um domínio quadrado unitário, utiliza-se o pacote de software MultiShape.
   • O método decompõe o domínio (com um "buraco" circular ou quadrado removido ao redor do ponto de colocalização) em uma malha de elementos quadriláteros.
   • Aplica-se um método pseudoespectral (baseado em pontos de Chebyshev) em cada elemento para calcular a ação do operador de convolução com alta precisão.
   • Isso permite pré-calculá uma matriz de convolução que é esparsa e eficiente, evitando a necessidade de calcular integrais densas a cada passo de tempo.

3. Integração Temporal: O sistema resultante de EDPs é resolvido usando um solver rígido (stiff solver) de passo variável (ode15s no MATLAB), tratando as condições de contorno e a não linearidade do potencial.

3. Contribuições Chave

• Novo Método de Aproximação: Desenvolvimento de um método de elementos espectrais específico para núcleos de convolução singulares, utilizando uma decomposição analítica da singularidade e uma malha adaptativa ("multishape") para a parte regular.
• Análise de Erro Rigorosa: Estabelecimento de estimativas de erro que mostram uma convergência de ordem $O(\epsilon^2 \log \epsilon)$ para a aproximação da parte singular, além de convergência espectral para a parte regular.
• Extensão Teórica: Prova de taxas de convergência para soluções fracas em 3D e em domínios com cantos (Lipschitz), preenchendo lacunas na literatura existente.
• Código de Código Aberto: Disponibilização de uma implementação completa e reprodutível no GitHub.

4. Resultados Numéricos

Os autores realizaram extensos testes numéricos validando o método:

• Validação da Aproximação: Testes com funções constantes e soluções analíticas mostraram que o erro da aproximação da convolução é extremamente baixo (próximo à precisão da máquina para $\epsilon$ suficientemente pequeno) e independente da malha.
• Dinâmica de Separação de Fases:
  • Simulações com diferentes escalas do núcleo newtoniano ($\eta$) demonstraram regimes variados: difusão pura, separação de fases lenta e rápida, e formação de estados estacionários não triviais.
  • Observou-se a formação de estruturas ricas e padrões complexos, especialmente quando o núcleo é uma mistura de potencial newtoniano e um mollifier (suavizador), permitindo valores positivos e negativos no núcleo.
• Potencial Regularizado: Introduziram um potencial regularizado $F_\omega$ para simulações de longo prazo. Os resultados mostraram que o uso do potencial regularizado acelera significativamente a convergência para o estado estacionário (reduzindo o custo computacional) sem perder precisão, pois remove a rigidez associada à singularidade logarítmica perto das fases puras.
• Eficiência: O método conseguiu obter soluções de alta resolução com recursos computacionais limitados (um laptop comum), demonstrando a eficácia da pré-cálculo da matriz de convolução.

5. Significado e Impacto

Este trabalho é significativo por fornecer a primeira abordagem numérica robusta e de alta resolução para o sistema de Cahn–Hilliard não local com potencial logarítmico singular e núcleo newtoniano singular em domínios gerais.

• Superação de Limitações: Resolve o problema da singularidade dupla (no potencial e no núcleo) que dificultava simulações anteriores.
• Aplicabilidade: O método é flexível para diferentes núcleos e geometrias de domínio, sendo aplicável a problemas de dinâmica de fluidos complexos, crescimento tumoral (modelo Hawkins-Daruud) e sistemas de reação-difusão.
• Ferramenta para Controle: A eficiência do método abre caminho para problemas de controle ótimo em sistemas não locais, que eram anteriormente computacionalmente inviáveis.

Em resumo, o artigo estabelece um novo padrão para a simulação numérica de modelos de separação de fases não locais, combinando rigor analítico com eficiência computacional através de técnicas espectrais avançadas.