The error estimate of entropy-stable discontinuous Galerkin methods for hyperbolic conservation laws

Resumo Técnico: Estimativas de Erro de Métodos Discontínuos de Galerkin Entropicamente Estáveis para Leis de Conservação Hiperbólicas

Definição do Problema
Leis de conservação hiperbólicas governam fenômenos físicos como a dinâmica de gases, onde as soluções frequentemente desenvolvem choques e descontinuidades de contato em tempo finito. Para garantir a bem-postura desses problemas, são necessárias desigualdades de entropia para selecionar a solução fraca fisicamente admissível. Embora Chen e Shu [J. Comput. Phys. 345 (2017)] tenham estabelecido um arcabouço unificado para a construção de métodos de Discontinuous Galerkin (DG) de alta ordem com estabilidade de entropia (ESDG) usando integração numérica específica e operadores de soma por partes (SBP), uma teoria de convergência rigorosa estava ausente. As análises de erro existentes para métodos DG padrão baseiam-se na suposição de que a solução numérica reside em um espaço prescrito de polinômios por partes. No entanto, para o método ESDG nodal, os valores nodais produzidos pelo esquema não correspondem, geralmente, a um polinômio por partes de grau $k$ em múltiplas dimensões espaciais. Consequentemente, as estimativas de erro clássicas não são diretamente aplicáveis, deixando uma lacuna na justificativa teórica para a convergência desses esquemas de estabilidade de entropia.

Metodologia
Este artigo estabelece estimativas de erro a priori rigorosas para métodos ESDG semi-discretos em malhas não estruturadas gerais. A análise emprega um argumento de consistência-estabilidade do tipo diferença finita formulado diretamente no nível nodal, em vez de depender de argumentos variacionais padrão.

A metodologia central envolve:

1. Hipótese de Reconstrução Polinomial: Os autores assumem que, para cada elemento, o vetor nodal admite uma reconstrução polinomial única em um espaço de dimensão finita $V_h(K)$ contendo polinômios de grau $k$. Isso permite o uso de estimativas inversas padrão e teoria de aproximação.
2. Limite a priori em $L^\infty$: A análise assume um limite no erro na norma $L^\infty$, o qual é justificado via um argumento de continuidade (método de continuação).
3. Arcabouço de Consistência-Estabilidade:
   • Consistência: Um erro de truncamento local é derivado inserindo os valores nodais da solução exata no esquema semi-discreto.
   • Estabilidade: Uma equação de erro é formulada subtraindo o esquema numérico da evolução da solução exata. O erro é medido em uma norma baseada em quadratura equivalente à norma $L^2$ quebrada no espaço de reconstrução.
   • Estimativas de Energia: A equação de erro é multiplicada pelo erro da variável de entropia e integrada. Os termos de volume são controlados usando a estrutura simétrica de SBP combinada com a simetrização de entropia dos Jacobianos de fluxo. As contribuições de interface são gerenciadas via a estabilidade de entropia do fluxo numérico. Termos não lineares são estimados usando expansões de Taylor e limites de norma de matriz escalonados pela malha.
   • Fechamento: A desigualdade diferencial resultante é fechada usando um argumento de Gronwall.

O arcabouço é estendido ao método DG Livre de Oscilações com Estabilidade de Entropia (ESOFDG), que inclui termos de amortecimento adicionais para suprimir oscilações espúrias. Os autores provam que esses termos de amortecimento não degradam a ordem de convergência.

Principais Contribuições

• Teoria de Convergência Rigorosa: O artigo fornece a primeira estimativa de erro a priori rigorosa para o método ESDG semi-discreto introduzido por Chen e Shu [6], provando uma estimativa de erro de $O(h^k)$ em uma norma baseada em quadratura para soluções suaves de tanto leis de conservação escalares quanto sistemas.
• Resolução da Lacuna da Suposição Polinomial: Ao introduzir a hipótese de reconstrução polinomial, os autores preenchem a lacuna entre a formulação de diferenciação de fluxo nodal (que viola a suposição padrão de polinômio por partes) e as técnicas clássicas de análise de erro.
• Extensão ao ESOFDG: A análise é estendida com sucesso ao método ESOFDG, demonstrando que os termos de amortecimento adicionais, projetados para controlar oscilações, preservam a taxa de convergência teórica.
• Lemas Técnicos: O artigo estabelece resultados técnicos necessários em relação às propriedades de escala dos operadores SBP, erros de truncamento local e limites nos coeficientes de amortecimento para o esquema ESOFDG.

Resultados

• Limites Teóricos: Sob as suposições de reconstrução polinomial e um limite a priori em $L^\infty$, o artigo prova que o erro na norma baseada em quadratura escala como $O(h^k)$ para o grau polinomial $k$. Este resultado é válido tanto para equações escalares quanto para sistemas de leis de conservação.
• Experimentos Numéricos: Testes numéricos na equação de Burgers (1D e 2D) e nas equações de Euler 2D (evolução de vórtice) confirmam os achados teóricos.
  • As taxas de convergência observadas geralmente alinham-se ou excedem ligeiramente a previsão teórica de $O(h^k)$.
  • Em alguns casos (ex: Burgers 1D com $k=3$), as taxas observadas aproximam-se de $O(h^{k+1})$, sugerindo que o limite teórico pode não ser agudo (sharp).
  • O esquema ESOFDG exibe taxas de convergência indistinguíveis do esquema ESDG padrão, confirmando que os termos de amortecimento não sacrificam a precisão.

Significância e Alegações
Os autores alegam que este trabalho fecha uma lacuna crítica na teoria de convergência para esquemas DG com estabilidade de entropia. Ao estabelecer uma estimativa de erro $O(h^k)$ rigorosa, o artigo fornece uma base teórica para a confiabilidade desses métodos de alta ordem na captura de soluções fisicamente admissíveis.

O artigo reconhece modestamente que, embora os experimentos numéricos sugiram que as taxas observadas possam exceder o limite teórico em até meio de ordem (potencialmente indicando superconvergência ou mecanismos delicados de cancelamento), a análise atual fornece um limite inferior sólido, embora potencialmente não agudo. Os autores observam que se uma estimativa de ordem $(k + 1/2)$ mais aguda pode ser derivada dentro do arcabouço variacional de elementos finitos permanece uma questão aberta. A principal contribuição é o estabelecimento de uma teoria de convergência rigorosa para uma classe de esquemas que anteriormente careciam de tal justificativa.