Novel Adaptive Methods for Hyperbolic Conservation Laws Based on New Quasi-Linear Seventh- and Ninth-Order Schemes

Resumo Técnico: Novos Métodos Adaptativos para Leis de Conservação Hiperbólicas Baseados em Novos Esquemas Quasi-Lineares de Sétima e Nona Ordem

Enunciado do Problema
O artigo aborda o desafio de resolver numericamente sistemas de leis de conservação hiperbólicas unidimensionais e bidimensionais, especificamente as equações de Euler da dinâmica de gases. Esses sistemas frequentemente desenvolvem estruturas complexas envolvendo choques, rarefações e descontinuidades de contato. Uma dificuldade primária na simulação numérica é alcançar simultaneamente (i) uma resolução nítida e não oscilatória de descontinuidades e (ii) alta precisão em regiões suaves. Métodos de ordem superior tradicionais frequentemente dependem de limitadores não lineares computacionalmente caros que são desnecessários em regiões suaves, enquanto métodos de ordem inferior podem carecer da resolução necessária para estruturas suaves.

Metodologia
Os autores propõem um novo arcabouço numérico adaptativo que particiona automaticamente o domínio computacional em regiões suaves e não suaves ("rugosas") utilizando um indicador de suavidade (SI). A estratégia, baseando-se no trabalho de [3], emprega diferentes esquemas numéricos baseados nas características locais da solução:

1. Regiões Suaves: A inovação central deste artigo é a substituição do esquema de diferenças finitas quasi-linear (QL) de quinta ordem (usado em [3]) por novos esquemas QL de sétima e nona ordem [10]. Esses esquemas são formulados como métodos de passo de tempo de três estágios (ordem espacial $r$, ordem temporal 3). Eles utilizam aproximações de diferenças centrais de ordem superior para as derivadas de fluxo e incorporam um termo de viscosidade artificial de ordem superior para estabilizar a solução. A força dessa viscosidade é controlada por um parâmetro $B$, determinado via análise de estabilidade de von Neumann.

   • Para o esquema de sétima ordem ($r=7$), o parâmetro de viscosidade ideal foi encontrado como $B^* = 1/256$.
   • Para o esquema de nona ordem ($r=9$), o parâmetro de viscosidade ideal é $B^* = 1/1024$.
   • A análise de estabilidade indica que, embora os esquemas sejam estáveis para números CFL moderados, o esquema de nona ordem admite uma região de estabilidade maior (até $\lambda = 0,75$) em comparação ao esquema de sétima ordem.

2. Regiões Não Suaves ("Rugosas"): O domínio é subdividido em áreas de contato e o restante da região rugosa.

   • Superfícies de Contato: Identificadas onde há saltos de densidade mas a pressão permanece contínua, estas áreas são tratadas com esquemas de upwind central de baixa dissipação (LDCU) de segunda ordem [6] equipados com um limitador sobrecompressivo [14].
   • Outras Áreas Rugosas (Choques): Estas regiões utilizam o limitador dissipativo Minmod2 [14].
   • Ambos os limitadores são aplicados às variáveis características locais via decomposição característica local para suprimir oscilações espúrias.

3. Dimensionalidade: Os esquemas 1-D são estendidos para 2-D usando uma construção dimensão-por-dimensão, mantendo os mesmos parâmetros de viscosidade motivados pela estabilidade.

Principais Contribuições

• Desenvolvimento de Esquemas de Ordem Superior: O artigo introduz e integra novos esquemas de diferenças finitas QL de sétima e nona ordem em um arcabouço adaptativo.
• Análise de Estabilidade: Uma análise de estabilidade de von Neumann rigorosa é conduzida para determinar os parâmetros ideais de viscosidade artificial ($B^*$) para ambos os esquemas de sétima e nona ordem, garantindo o amortecimento de oscilações de escala de grade sem comprometer a estabilidade para números CFL razoáveis.
• Refinamento da Estratégia Adaptativa: Os autores demonstram que substituir o esquema de quinta ordem por esquemas de sétima e nona ordem em regiões suaves de um algoritmo adaptativo existente aumenta significativamente a resolução sem alterar o tratamento robusto de descontinuidades.

Resultados
O desempenho dos novos esquemas adaptativos (referidos como 5ª Ordem, 7ª Ordem e 9ª Ordem) foi avaliado através de uma série de testes de referência 1-D e 2-D para as equações de Euler:

• Verificação de Precisão: Testes de convergência em soluções suaves confirmaram que os esquemas atingem suas ordens teóricas de precisão (7ª e 9ª) tanto em casos 1-D quanto 2-D.
• Interação Choque-Entropia (1-D): No problema de Titarev-Toro, os esquemas de 7ª e 9ª ordem demonstraram maior precisão na resolução de perturbações de densidade de alta frequência em comparação ao esquema de 5ª ordem.
• Onda de Detonação (Blast Wave) (1-D): Todos os esquemas resolveram com sucesso interações de choque fortes sem oscilações espúrias.
• Problemas de Riemann 2-D: Em configurações envolvendo instabilidades de jato e estruturas de vórtice complexas, os esquemas de 7ª e 9ª ordem forneceram resolução superior de características de pequena escala e instabilidades de banda lateral em comparação ao esquema de 5ª ordem.
• Problema de Implosão: Os esquemas de ordem superior reduziram a dissipação numérica, permitindo que o jato se propagasse mais longe ao longo da diagonal comparado ao esquema de 5ª ordem.
• Instabilidade de Rayleigh-Taylor (2-D): Os esquemas de 7ª e 9ª ordem resolveram mais estruturas de solução de pequena escala, indicando menor dissipação numérica.
• Custo Computacional: O tempo de CPU para os esquemas de 7ª e 9ª ordem aumentou em apenas aproximadamente 7,0% e 11,7%, respectivamente, em relação ao esquema de 5ª ordem.

Significância
O artigo afirma que os esquemas adaptativos propostos, ao incorporar os novos esquemas QL de sétima e nona ordem em regiões suaves, alcançam maior resolução e contêm uma quantidade menor de dissipação numérica em comparação com seus equivalentes utilizando o esquema QL de quinta ordem. Crucialmente, essas melhorias são alcançadas mantendo a robustez e a resolução nítida de descontinuidades proporcionadas pela estratégia adaptativa e exigindo apenas um aumento marginal no custo computacional. Os resultados sugerem que os novos esquemas oferecem um equilíbrio mais eficiente entre precisão e custo computacional para simular problemas complexos de dinâmica de gases.