A collocation scheme that is equivalent to discontinuous Galerkin discretizations

Resumo Técnico: Um Esquema de Colocação Equivalente a Discretizações de Galerkin Descontínuas

Enunciado do Problema
O artigo aborda a construção de discretizações estáveis e de alta ordem para equações diferenciais parciais (EDPs) em malhas não estruturadas. Embora os operadores de Somatório-Por-Partes (SBP) forneçam um quadro robusto para estabilidade, os operadores SBP multidimensionais em elementos simplex (como triângulos) podem ser complexos de construir e frequentemente carecem de expressões de forma fechada. Trabalhos anteriores de Chan introduziram um operador SBP baseado em projeção (denominado aqui como Operador de Colocação Modal ou MC) como um passo intermediário em discretizações de elementos finitos. No entanto, as propriedades específicas do uso direto deste operador em um esquema de colocação, particularmente sua relação com os métodos de Galerkin Descontínuos (DG), não haviam sido totalmente exploradas. A questão central é se uma semidiscretização baseada neste operador MC específico produz soluções equivalentes a uma discretização modal DG usando a mesma quadratura subjacente.

Metodologia
O autor define o operador MC utilizando bases polinomiais ortonormais e regras de quadratura com pontos interiores positivos. A metodologia central envolve:

1. Definição do Operador: Definir o operador MC de primeira derivada $D_{xd} = V_{xd}V^T W$, onde $V$ é a matriz de Vandermonde das funções de base ortonormais, $V_{xd}$ contém suas derivadas e $W$ é a matriz diagonal dos pesos de quadratura.
2. Prova de Equivalência Teórica: O artigo constrói uma prova demonstrando que a semidiscretização de um problema de valor inicial e de contorno (IVBP) hiperbólico não linear escalar usando o operador MC produz um vetor de solução $u(t)$ que é exatamente a projeção do vetor de solução modal DG $\tilde{u}(t)$ sobre o espaço nodal ($u(t) = V\tilde{u}(t)$). Esta prova baseia-se na decomposição do espaço de soluções no subespaço polinomial e no espaço nulo do operador.
3. Análise Espectral: O autor analisa os autovalores das matrizes do sistema resultante. Mostra-se que o espectro do operador MC consiste nos autovalores do operador DG correspondente mais um conjunto de autovalores nulos associados ao espaço nulo.
4. Tratamento da Inconsistência do Espaço Nulo: O artigo reconhece que os operadores MC são geralmente inconsistentes no espaço nulo (seu espaço nulo contém mais do que apenas vetores constantes). Para problemas não estacionários, esses modos extrínsecos são inicializados como zero e permanecem zero. Para problemas estacionários, o autor propõe Estabilização por Projeção Local (LPS) para remover a singularidade do espaço nulo sem afetar a exatidão polinomial.
5. Verificação Numérica: A teoria é verificada através de experimentos numéricos em malhas triangulares usando:
   • A equação de advecção com coeficientes constantes.
   • A equação não linear de Burgers.
   • Comparações entre discretizações MC "Padrão" e formulações "Conservadoras de Entropia" (antissimétricas).
   • Testes envolvendo regras de quadratura com pesos negativos.

Principais Contribuições

• Teorema de Equivalência: A contribuição principal é a prova de que uma discretização de colocação modal usando o operador MC é matematicamente equivalente a uma discretização modal DG, desde que compartilhem a mesma regra de quadratura. Isso vale independentemente do número de graus de liberdade (pontos de quadratura) no esquema de colocação.
• Propriedades Espectrais: O artigo estabelece que o raio espectral do operador MC é efetivamente determinado pelo método DG subjacente, com a adição de autovalores nulos. Isso implica que o agrupamento de nós (frequentemente problemático para métodos de elementos espectrais em triângulos) não impacta negativamente o raio espectral dos operadores MC da mesma forma que impacta os operadores SBP nodais tradicionais.
• Gerenciamento do Espaço Nulo: O trabalho esclarece que, embora os operadores MC sejam inconsistentes no espaço nulo, isso não é prejudicial para simulações dependentes do tempo se as condições iniciais forem projetadas corretamente. Para problemas de estado estacionário, o artigo demonstra como recuperar a consistência via dissipação numérica baseada em projeção (LPS).
• Nuance da Conservação de Entropia: O estudo revela que discretizações MC conservadoras de entropia (baseadas em divisões antissimétricas) não são inerentemente equivalentes a DG porque o resíduo pode não residir dentro do espaço polinomial. No entanto, a equivalência pode ser recuperada projetando o resíduo sobre a base polinomial.

Resultados

• Advecção com Coeficientes Constantes: Os resultados numéricos mostram que a diferença $L_2$ entre as soluções MC e modal DG está na precisão da máquina ($\approx 10^{-15}$) para vários graus polinomiais e níveis de exatidão de quadratura.
• Raios Espectrais: As distribuições de autovalores confirmam que os operadores MC exibem raios espectrais favoráveis em comparação com operadores SBP de norma mínima, que sofrem crescimento catastrófico devido ao agrupamento de nós em triângulos.
• Pesos Negativos: A equivalência mantém-se mesmo quando a regra de quadratura inclui pesos negativos, um cenário onde a teoria convencional de estabilidade energética SBP pode falhar em prever estabilidade. A estabilidade é, em vez disso, garantida pela equivalência com a solução DG energeticamente estável.
• Equação de Burgers: Para a equação não linear de Burgers, a discretização MC "Padrão" é equivalente a DG. A variante "Conservadora de Entropia" inicialmente mostra uma discrepância porque o resíduo não é abrangido pela base polinomial. No entanto, uma vez que o resíduo é projetado sobre a base, a equivalência é restaurada até a precisão da máquina.

Significado e Alegações
O artigo alega que o Operador de Colocação Modal introduzido por Chan não é meramente uma ferramenta intermediária, mas um esquema de discretização válido que é efetivamente um método DG modal. O significado reside na simplificação da construção de esquemas de alta ordem e estáveis em malhas simplex não estruturadas. Ao usar o operador MC, pode-se aproveitar a expressão de forma fechada do operador, mantendo ao mesmo tempo a estabilidade e as propriedades espectrais dos métodos DG.

O autor conclui modestamente que, embora os dois métodos sejam equivalentes, a escolha entre eles depende da aplicação. Para problemas lineares onde as matrizes de massa e rigidez são pré-calculadas, o DG pode ser preferível. No entanto, para problemas não lineares resolvidos com marchagem temporal explícita ou esquemas implícitos livres de matriz, o operador MC pode oferecer vantagens computacionais. O artigo defende fortemente o uso da colocação modal como o método preferido para a construção de operadores SBP multidimensionais.