High Order semi-implicit Rosenbrock type and Multistep methods for evolutionary partial differential equations with higher order derivatives

Título: Métodos Semi-Implicitos de Alta Ordem do Tipo Rosenbrock e Multietapas para Equações Diferenciais Parciais Evolutivas com Derivadas de Ordem Superior

Autores: Sebastiano Boscarino e Giuseppe Izzo (Universidade de Catania e Universidade de Nápoles, Itália).

---

1. Problema Abordado

O artigo foca na resolução numérica de Equações Diferenciais Parciais (EDPs) dependentes do tempo que envolvem derivadas espaciais de ordem elevada ($k \ge 2$). Exemplos citados incluem:

• Problemas de difusão de segunda ordem (ex: deslocamento miscível em meios porosos).
• Equações de convecção-difusão não lineares.
• Equações dispersivas de terceira ordem (ex: Equação de Korteweg-de Vries - KdV e suas variantes generalizadas).
• Equações de difusão de quarta ordem (ex: equações do tipo bi-harmônica para modelagem de vigas finas e separação de fases).

O Desafio:

• Métodos Explícitos: Sofrem de restrições severas de estabilidade no passo de tempo ($\Delta t$), exigindo que $\Delta t = O(\Delta x^k)$. Para ordens altas ($k \ge 2$), isso torna o passo de tempo extremamente pequeno, resultando em custos computacionais proibitivos.
• Métodos Implícitos Totais: Superam as restrições de estabilidade, mas exigem a resolução de um sistema algébrico não linear em cada passo de tempo (geralmente via iteração de Newton), o que é computacionalmente caro e complexo.

O objetivo é desenvolver esquemas que combinem a estabilidade dos métodos implícitos com a simplicidade computacional dos métodos explícitos, evitando iterações não lineares.

---

2. Metodologia Proposta

Os autores propõem uma estratégia Semi-Implicita (SI) integrada a dois frameworks: métodos do tipo Rosenbrock e métodos Multietapas (LM).

Conceito Central (Estratégia SI)

A EDP semi-discretizada é escrita na forma:
$$\frac{dU}{dt} = F(U) + B(U)U$$
Onde $F(U)$ representa termos não rígidos (explícitos) e $B(U)U$ representa os termos rígidos de alta ordem.
Na estratégia SI:

• O termo $F(U)$ e a matriz $B(U)$ são tratados explicitamente.
• Apenas o vetor $U$ dentro do termo rígido é tratado implicitamente.
• Isso transforma o sistema não linear em um sistema linear que pode ser resolvido sem iterações de Newton.

Abordagem 1: Métodos SI do Tipo Rosenbrock

• Baseados em esquemas Runge-Kutta aditivos semi-implícitos.
• Linearizam o termo rígido em torno da solução no tempo atual $t_n$.
• Utilizam a mesma matriz Jacobiana $J = B(U_n)$ para todos os estágios do método, permitindo uma única decomposição LU por passo de tempo.
• Contribuição: Derivação de condições de ordem até a terceira ordem ($p=3$) e construção de esquemas com propriedades de estabilidade $L$-estável (estabilidade estrita para problemas rígidos).

Abordagem 2: Métodos Multietapas Semi-Implicitos (SI-LM) com Preditor-Corretor

• O artigo identifica uma falha em esquemas SI-LM anteriores: o uso de um preditor puramente explícito pode reintroduzir instabilidades para derivadas de ordem superior ($k > 2$).
• Inovação: Propõe um esquema Preditor-Corretor Modificado onde ambos o preditor e o corretor são tratados de forma semi-implícita.
  • Preditor: Usa um método SI de ordem $p^*$ (ex: Euler SI).
  • Corretor: Usa um método SI de ordem $p$ (ex: BDF - Backward Differentiation Formula).
• Iteração: Para atingir a ordem de convergência esperada $p$, realiza-se um número fixo de iterações de correção ($\mu$). O artigo demonstra que $\mu \approx p - p^*$ iterações são suficientes para recuperar a ordem e a constante de erro do corretor.
• Vantagem: Evita a instabilidade do preditor explícito e mantém o custo computacional baixo (resolvendo apenas sistemas lineares).

---

3. Principais Contribuições

1. Desenvolvimento de Esquemas de Alta Ordem: Construção de métodos SI-Rosenbrock de ordem 3 e métodos SI-LM (baseados em BDF) de ordens 2, 3 e 4.
2. Estabilidade para Derivadas de Alta Ordem: Demonstração teórica e prática de que a estratégia SI, quando aplicada corretamente (especialmente no preditor dos métodos multietapas), remove as restrições severas de $\Delta t$ associadas a derivadas de ordem $k \ge 2$.
3. Análise de Condições de Ordem: Derivação detalhada das condições de ordem para os métodos Rosenbrock SI até a terceira ordem, incluindo a análise de estabilidade e regiões de estabilidade $L$-estáveis.
4. Validação Numérica: Extensa bateria de testes numéricos em equações de difusão, convecção-difusão, KdV e bi-harmônicas, demonstrando que os esquemas atingem as ordens de precisão esperadas e são estáveis para passos de tempo grandes.

---

4. Resultados Numéricos

Os autores realizaram experimentos numéricos em 1D com condições de contorno periódicas, utilizando esquemas de diferenças finitas de alta ordem (WENO) para a discretização espacial.

• Testes de Ordem de Precisão:
  • Os métodos SI-Rosenbrock ($p=3$) e SI-LM ($p=2,3,4$) mostraram taxas de convergência consistentes com a ordem teórica em normas $L_1$, $L_2$ e $L_\infty$.
  • Exemplo: Para a equação KdV generalizada e a equação bi-harmônica, os erros diminuíram conforme o refinamento da malha, confirmando a precisão de ordem 3 e 4.
• Estabilidade:
  • Os esquemas permaneceram estáveis para passos de tempo $\Delta t$ muito maiores do que os permitidos por métodos explícitos (que exigiriam $\Delta t \propto \Delta x^k$).
  • A estratégia de preditor-corretor semi-implícito (SI-PC) provou ser robusta, evitando a instabilidade observada em preditores puramente explícitos para problemas de alta ordem.
• Eficiência:
  • Os métodos multietapas SI-PC exigem a resolução de $p$ ou $p+1$ sistemas lineares por passo de tempo, o que é computacionalmente mais eficiente do que métodos Runge-Kutta de alta ordem que exigem mais estágios (e mais sistemas lineares) para atingir a mesma ordem.

---

5. Significância e Conclusão

O trabalho oferece uma solução prática e eficiente para um problema clássico na simulação numérica de EDPs: o custo computacional de métodos implícitos versus a instabilidade de métodos explícitos em problemas de alta ordem.

• Flexibilidade: A estratégia SI pode ser acoplada a diversos esquemas de discretização espacial (diferenças finitas, Galerkin descontínuo, etc.).
• Simplicidade: Elimina a necessidade de iterações de Newton, tornando a implementação mais simples e robusta.
• Aplicabilidade: Os métodos são particularmente úteis para modelagem de fenômenos físicos complexos como separação de fases, dinâmica de fluidos dispersivos e elasticidade de gradiente de tensão, onde derivadas de alta ordem são comuns.

Em resumo, os autores demonstram que é possível construir esquemas de alta ordem, estáveis e livres de iterações não lineares para EDPs evolutivas com derivadas espaciais de alta ordem, superando as limitações dos métodos tradicionais.