The jump filter in the discontinuous Galerkin method for hyperbolic conservation laws

Resumo Técnico: O Filtro de Salto (Jump Filter) no Método de Galerkin Descontínuo para Leis de Conservação Hiperbólicas

Definição do Problema
Esquemas numéricos lineares de alta ordem, particularmente o método de Galerkin Descontínuo (DG), são amplamente utilizados para resolver leis de conservação hiperbólicas devido à sua alta precisão e flexibilidade com malhas não estruturadas. No entanto, ao simular soluções contendo descontinuidades (como ondas de choque), esses esquemas frequentemente geram oscilações espúrias indesejáveis (o fenômeno de Gibbs). Essas oscilações podem levar à divergência do solver ou à criação de estruturas artificiais. Embora existam várias técnicas de captura de choque — como limitadores de Variação Total Limitada (TVB), limitadores WENO e métodos de viscosidade artificial — muitas sofrem de desvantagens, incluindo alto custo computacional, dificuldade de implementação em malhas não estruturadas ou curvas, a necessidade de decomposição característica ou dissipação numérica excessiva que degrada a precisão em regiões suaves.

Metodologia
Os autores propõem um novo filtro de salto (jump filter) projetado para operar dentro da estrutura DG para mitigar essas oscilações enquanto preserva a precisão de alta ordem e a compacidade do método.

• Formulação de Viscosidade Local: O núcleo do método é a introdução de um termo de viscosidade local baseado em informações de salto nas interfaces das células. Diferente da viscosidade espectral global, esta abordagem é localizada. O coeficiente de viscosidade, $\sigma_j$, é determinado adaptativamente pela magnitude dos saltos na solução e em suas derivadas através das fronteiras das células.
  • Em regiões suaves onde os saltos são mínimos, a viscosidade adicionada é desprezível, preservando a precisão de alta ordem.
  • Em regiões com descontinuidades fortes, a viscosidade aumenta para amortecer os modos polinomiais de alta ordem responsáveis pelas oscilações.
• Filtragem Espectral via Separação de Tempo (Time-Splitting): Para evitar a rigidez (stiffness) e o custo computacional associados à resolução implícita do termo de difusão (o que restringiria o passo de tempo), os autores empregam uma estratégia de separação de tempo. Após cada estágio do esquema DG de Runge-Kutta explícito padrão, o filtro de salto é aplicado como um passo de pós-processamento. Isso envolve resolver uma Equação Diferencial Ordinária (ODE) simples para os coeficientes modais da solução polinomial:
  $$\frac{d}{dt} u^l_j + \sigma^l_j u^l_j = 0$$
  A solução é atualizada multiplicando os coeficientes por um fator de decaimento exponencial, $e^{-\sigma^l_j \Delta t}$. Esta operação é explícita, computacionalmente barata e não requer decomposição característica.
• Integração de Limitador Híbrido: O filtro de salto também é integrado em uma estrutura de limitador híbrido. Nesta abordagem, células "problemáticas" são identificadas usando um indicador de salto. Para essas células, o filtro de salto é usado para computar uma solução polinomial de ordem inferior, que é então misturada com a solução DG de alta ordem. Isso visa reduzir ainda mais a dissipação numérica em comparação ao uso do filtro isoladamente.
• Aplicabilidade: O método é formulado tanto para leis de conservação escalares quanto para sistemas (especificamente as equações de Euler) em uma e duas dimensões. Ele opera diretamente no espaço das variáveis físicas (conservativas), evitando as complexidades da decomposição característica.

Principais Contribuições

1. Novo Filtro de Salto: Introdução de um filtro de viscosidade local que utiliza informações de salto nas interfaces das células para amortecer adaptativamente os modos de alta ordem.
2. Preservação das Propriedades DG: A análise teórica confirma que o esquema proposto mantém a conservação, a estabilidade $L^2$ e as taxas de convergência ótimas (precisão de alta ordem) para soluções suaves.
3. Eficiência Computacional: O método é projetado para baixo custo computacional. Ele não requer decomposição característica, opera inteiramente no espaço físico e mantém a compacidade do esquema DG (comunicação de dados local), tornando-o altamente adequado para ambientes de computação paralela e GPUs.
4. Integração Híbrida: Demonstração da integração perfeita do filtro em uma estrutura de limitador híbrido, suprimindo efetivamente as oscilações em choques fortes enquanto minimiza a dissipação numérica em regiões suaves.

Resultados Numéricos
O artigo valida o método através de extensos experimentos numéricos nas equações de Euler 1D e 2D, cobrindo:

• Fluxo Suave: Testes de precisão (ex: equação de Burgers, vórtice isentrópico) confirmam que o filtro não degrada a ordem de convergência do esquema DG padrão em malhas finas.
• Choques Fortes: Testes envolvendo o problema de Lax, ondas de detonação (blast waves) e reflexão de Mach dupla mostram que o filtro controla efetivamente as oscilações espúrias e captura perfis de choque nítidos sem os excessos (overshoots) típicos de esquemas de alta ordem não filtrados.
• Interações Complexas: O método apresenta bom desempenho em cenários de interação choque-vórtice, interação choque-bolha e instabilidade de Kelvin-Helmholtz, onde esquemas de ordem superior com o filtro resolvem estruturas de fluxo intrincadas (como linhas de deslizamento e camadas de mistura turbulenta) melhor do que contrapartes de ordem inferior.
• Robustez: O esquema demonstra robustez em simulações de jatos astrofísicos de alto número de Mach e mantém a positividade quando combinado com um limitador de preservação de positividade.

Significância e Alegações
Os autores alegam que o filtro de salto oferece uma alternativa robusta, precisa e eficiente aos limitadores tradicionais e métodos de viscosidade artificial. Sua principal significância reside na sua simplicidade e eficiência computacional:

• Evita a custosa decomposição característica exigida por muitos outros limitadores de alta ordem.
• Preserva a compacidade do método DG, facilitando a paralelização eficiente.
• Fornece uma maneira contínua de introduzir a dissipação necessária apenas onde ela é necessária, equilibrando a capacidade de captura de choque com a preservação da precisão de alta ordem em regiões suaves.

O artigo conclui que esta abordagem aborda efetivamente o desafio de simular leis de conservação hiperbólicas com soluções descontínuas, ofereando uma ferramenta prática para captura de choque de alta resolução. O trabalho futuro nota a exploração da implementação em malhas não estruturadas e problemas de estado estacionário.