Stable $C^1$-conforming finite element methods for a class of nonlinear fourth-order evolution equations

Resumo Técnico: Métodos de Elementos Finitos Estáveis e Conformes em $C^1$ para uma Classe de Equações de Evolução de Quarta Ordem Não Lineares

Formulação do Problema
O artigo aborda a análise numérica de uma classe geral de equações diferenciais parciais (EDPs) de evolução de quarta ordem não lineares definidas em um domínio limitado $\Omega \subset \mathbb{R}^d$ ($d \in \{1, 2, 3\}$). A equação governante é dada por:
$$\frac{\partial u}{\partial t} - \beta_1 \Delta u + \beta_2 \Delta^2 u = \beta_3(1 - |u|^2)u - \beta_4 u \times \Delta u + \beta_5 \Delta(|u|^2 u) + \beta_6 (j \cdot \nabla)u$$
sujeita a condições de contorno de Neumann ($\frac{\partial u}{\partial n} = 0, \frac{\partial (\Delta u)}{\partial n} = 0$) e uma condição inicial $u(0) = u_0$. Aqui, $u$ pode ser escalar-valued ($m=1$) ou vetor-valued ($m=3$). A equação abrange vários modelos físicos dependendo dos coeficientes $\beta_i$ e da dimensão $d$, incluindo:

• Micromagnetismo: A equação de Landau–Lifshitz–Baryakhtar (LLB) (com ou sem torque de spin) e regularizações da equação de Landau–Lifshitz–Bloch.
• Formação de Padrões: A equação de Swift–Hohenberg.
• Separação de Fases e Difusão: Equações de Cahn–Hilliard convectivas com termos de fonte, modelos de difusão de bi-fluxo anômalo e modelos de dispersão de população generalizados.

Os autores observam que, embora existam métodos numéricos para instâncias específicas dessas equações (por exemplo, Swift–Hohenberg escalar ou Cahn–Hilliard padrão), um esquema numérico unificado e geral para o sistema completo (1.1) com todos os coeficientes potencialmente não nulos era ausente na literatura.

Metodologia
Os autores propõem e analisam métodos de elementos finitos totalmente discretos utilizando elementos conformes em $C^1$. Esta escolha é necessária devido à derivada espacial de quarta ordem ($\Delta^2 u$), que exige que o espaço de elementos finitos seja um subespaço de $H^2(\Omega)$.

• Discretização Espacial: O método emprega polinômios por partes de grau $r$ (onde $r=3$ para $d=1,2$ e $r=9$ para $d=3$ para manter a continuidade $C^1$). Elementos específicos mencionados incluem elementos Hermite $C^1$-P3 (1D), macroelementos Hsieh–Clough–Tocher (HCT) (2D) e elementos de Zenísek (3D).
• Discretização Temporal: Dois esquemas de passo de tempo linearizados são propostos:
  1. Um método de Euler semi-implícito (primeira ordem no tempo).
  2. Um método BDF2 semi-implícito (segunda ordem no tempo).
     Em ambos os esquemas, os termos não lineares são linearizados avaliando-os no passo de tempo anterior, resultando em um sistema de equações lineares em cada passo. Isso evita o custo computacional de resolver sistemas não lineares via método de Newton.

Principais Contribuições e Resultados Teóricos
O artigo fornece uma análise de erro rigorosa e uma prova de estabilidade para os esquemas propostos:

1. Estabilidade:

   • Os esquemas são provados como sendo incondicionalmente estáveis na norma $H^2$ para dimensões $d=1$ e $d=2$. Isso significa que a estabilidade ocorre sem restrições no tamanho do passo de tempo $k$ em relação ao tamanho da malha $h$, e sem assumir triangulações quase uniformes.
   • Para $d=3$, os esquemas são condicionalmente estáveis. A estabilidade é garantida se os dados iniciais forem suficientemente pequenos ou se o passo de tempo satisfizer uma restrição específica em relação à energia inicial. Os autores reconhecem que elementos conformes em $C^1$ em 3D (exigindo polinômios de grau 9) são computacionalmente proibitivos, sugerindo que métodos não conformes são mais práticos para $d=3$, embora a análise teórica seja incluída para completude.

2. Estimativas de Erro e Convergência:

   • Os autores estabelecem taxas de convergência ótimas nas normas $L^\infty(0, T; H^\beta)$ para $\beta = 0, 1, 2$.
   • Para o esquema semi-discreto, o erro é limitado por $O(h^{r+1-\beta})$.
   • Para os esquemas totalmente discretos:
     • O método de Euler semi-implícito alcança um erro de $O(h^{r+1-\beta} + k)$.
     • O método BDF2 semi-implícito alcança um erro de $O(h^{r+1-\beta} + k^2)$.
   • A análise baseia-se no desenvolvimento de uma projeção elíptica especializada adaptada ao problema não linear e na estimativa cuidadosa dos termos não lineares usando imersões de Sobolev e desigualdades de inversão.

3. Experimentos Numéricos:

   • As taxas de convergência teóricas foram verificadas usando o pacote de código aberto FreeFEM com elementos HCT em 2D.
   • Experimentos foram conduzidos na equação de Landau–Lifshitz–Baryakhtar (vetor-valued) e na equação de Cahn–Hilliard com um termo de fonte (escalar-valued).
   • Os resultados confirmaram as ordens de convergência previstas: 4ª ordem em $L^2$, 3ª ordem em $H^1$ e 2ª ordem em $H^2$ para a discretização espacial, consistentes com o uso de elementos cúbicos ($r=3$).

Significância e Alegações
O artigo afirma fornecer o primeiro esquema numérico capaz de aproximar a solução para o sistema geral de quarta ordem não linear (1.1) em sua total generalidade. As principais vantagens destacadas são:

• Linearidade: Ao linearizar os termos não lineares, o método evita a complexidade e a potencial instabilidade de resolver sistemas não lineares em cada passo de tempo, tornando-o vantajoso para aplicações como simulações micromagnéticas.
• Estabilidade: A estabilidade incondicional em $H^2$ para $d=1,2$ garante um desempenho robusto sem restrições severas de passo de tempo.
• Precisão: O uso de elementos conformes em $C^1$ permite uma precisão de ordem superior para capturar os gradientes da solução em comparação com métodos não conformes, apesar do aumento no custo da continuidade $C^1$.

Os autores mantêm um tom modesto em relação ao caso 3D, declarando explicitamente que, embora a análise teórica esteja completa, a aplicação prática de elementos conformes em $C^1$ em 3D é atualmente computacionalmente inviável devido ao alto grau polinomial exigido. O trabalho serve como uma estrutura teórica fundamental para resolver equações de evolução de quarta ordem não lineares complexas que surgem na física, ciência dos materiais e biologia.