Investigation of Shock-Capturing with Bound-Preserving Limiters for the Nonlinearly Stable Flux Reconstruction Method

Resumo Técnico: Investigação de Captura de Choque com Limitadores de Preservação de Limites para o Método de Reconstrução de Fluxo Não Linearmente Estável

Definição do Problema
Esquemas numéricos de alta ordem para fluxos compressíveis, como o método de Discontinuous Galerkin (DG) e a Reconstrução de Fluxo (FR), oferecem alta precisão, mas frequentemente lutam com a robustez na presença de choques fortes. Embora o Fluxo de Reconstrução Entropicamente Estável (ESFR) e sua extensão, o Fluxo de Reconstrução Não Linearmente Estável (NSFR), forneçam estabilidade não linear e estabilidade de entropia prováveis, essas propriedades sozinhas não garantem a preservação de limites físicos, como densidade e pressão positivas. Estratégias existentes para mitigar esses problemas, como viscosidade artificial ou limitação subcelular, podem alterar a ordem de precisão, introduzir limitações no passo de tempo ou exigir implementações complexas. Além disso, as condições específicas de Courant-Friedrichs-Lewy (CFL) necessárias para manter a positividade para os fluxos de dois pontos usados em métodos FR não foram totalmente estabelecidas.

Metodologia
O artigo investiga a aplicação de um limitador de preservação de limites modificado ao esquema NSFR para resolver as equações de Euler compressíveis. A metodologia foca em três áreas primárias:

1. Limitador de Preservação de Limites Modificado: Os autores estendem o limitador de preservação de limites de Zhang-Shu. O limitador original e suas subsequentes melhorias por Wang e Shu garantiam a positividade em pontos de quadratura específicos, mas careciam de garantias de robustez nos nós da solução, potencialmente permitindo valores de densidade ou pressão negativos nesses nós. A modificação proposta incorpora os valores mínimos de densidade e pressão nos nós da solução (além dos pontos de quadratura padrão) ao determinar os fatores de escala ($\theta_1$ e $\theta_2$) usados para limitar a solução. Isso garante a positividade em todo o elemento, incluindo os nós da solução.
2. Configuração do Esquema NSFR: O estudo utiliza o esquema NSFR com uma forma de diferenciação de fluxo de dois pontos de simetria de inclinação (skew-symmetric). A implementação emprega o fluxo de dois pontos de Chandrashekar-Ranocha (CHRA), que é conservativo de entropia, preservador de energia cinética e preservador de equilíbrio de pressão. A discretização temporal utiliza o método Runge-Kutta de Estabilidade Forte Preservadora (SSPRK3) de terceira ordem.
3. Investigação Paramétrica: O artigo analisa sistematicamente o impacto de vários parâmetros do esquema nas capacidades de captura de choque:
   • Fluxos de Dois Pontos: Comparação entre os fluxos CHRA, Chandrashekar (CH), Ismail-Roe (IR) e Kennedy-Gruber (KG).
   • Nós de Quadratura: Comparação entre nós Gauss-Legendre (GL) e Gauss-Lobatto-Legendre (GLL) para avaliação de fluxo.
   • Parâmetro de Reconstrução de Fluxo (FR) ($c$): Investigação de diferentes valores do parâmetro de correção $c$, variando desde o valor padrão de DG ($c_{DG}=0$) até o de diferença espectral ($c_{SD}$), o de Huynh's $g_2$ ($c_{HU}$), e o valor máximo estável ($c_+$), bem como valores superiores a $c_+$.

Principais Contribuições

• Robustez do Limitador Aperfeiçoada: A principal contribuição é a modificação do limitador de Zhang-Shu para considerar os valores dos nós da solução, o que se mostrou necessário para a simulação bem-sucedida de problemas 2D dominados por choques fortes sem falhas.
• Verificação da Condição CFL: O artigo estabelece numericamente que os fluxos de dois pontos usados em NSFR aderem a condições CFL específicas para a preservação da positividade. Especificamente, uma condição CFL baseada na distância mínima entre os nós de quadratura GLL (Eq. 29) mostrou-se suficiente para a preservação da positividade em vários graus polinomiais tanto em 1D quanto em 2D.
• Capacidades de Captura de Choque do NSFR: O estudo demonstra que o esquema NSFR, quando pareado com o limitador modificado, pode resolver choques fortes com soluções essencialmente livres de oscilações, particularmente ao utilizar parâmetros de reconstrução de fluxo mais fortes (ex: $c_+$ ou superior).
• Análise Comparativa: O trabalho fornece uma comparação abrangente entre o NSFR e o método Strong DG padrão, destacando a robustez superior e as vantagens de passo de tempo da abordagem NSFR.

Resultados
A metodologia proposta foi validada através de uma série de casos de teste 1D e 2D:

• Precisão: O limitador modificado preserva a alta ordem de precisão do esquema, conforme verificado pelos estudos de convergência no caso de Baixa Densidade.
• Robustez e Vantagens de Passo de Tempo: Nos problemas de Tubo de Choque de Sod e Shu-Osher, o esquema NSFR com o limitador modificado demonstrou robustez significativamente maior que o método Strong DG. Por exemplo, no Tubo de Choque de Sod, o NSFR pôde operar sem o limitador em um CFL de 0,2, enquanto o Strong DG falhou em CFLs tão baixos quanto 0,01 sem um limitador. Na interação de Choque e Vórtice Forte, o NSFR permitiu um CFL de 0,58 com o limitador, comparado a 0,27 sem ele, mantendo a precisão da solução.
• Mitigação de Oscilações: Observou-se que o aumento do parâmetro de correção de FR ($c$) efetivamente amortece oscilações espúrias e sobresaltos no rastro dos choques. Por exemplo, nos casos do Tubo de Choque de Leblanc e Reflexão de Mach Dupla, o esquema $c_+$ produziu soluções mais limpas com ruído reduzido em comparação ao esquema $c_{DG}$.
• Casos Extremos: O esquema simulou com sucesso benchmarks extremos, incluindo os casos de Jato Astrofísico de Mach 80 e Mach 2000. No caso de Mach 2000, um parâmetro de reconstrução de fluxo dez vezes maior que $c_+$ permitiu que a simulação rodasse com um CFL de 0,0005, preservando a positividade onde os métodos DG padrão falharam.
• Sensibilidade de Fluxo e Nó: O fluxo de Kennedy-Gruber (KG) foi considerado insuficientemente dissipativo para choques fortes, levando ao crescimento de entropia e oscilações. Enquanto os nós GL ofereceram uma integração mais forte, os nós GLL proporcionaram uma vantagem significativa nos passos de tempo permitidos devido à redução da sensibilidade em projeções de entropia nas fronteiras.

Significância
O artigo afirma que o esquema NSFR modificado, combinado com o limitador de preservação de limites aperfeiçoado, oferece um framework robusto para simular fluxos compressíveis não estacionários com choques fortes. A significância reside na capacidade de alcançar soluções essencialmente livres de oscilações sem depender de limitadores TVD ou viscosidade artificial, que frequentemente degradam a precisão. O estudo destaca que o esquema NSFR não apenas mantém a estabilidade de entropia, mas também fornece vantagens práticas em termos de passos de tempo permitidos e fidelidade de captura de choque. Os autores sugerem que a capacidade de escalar o parâmetro FR com base em descontinuidades locais representa uma direção promissora para pesquisas futuras, visando maximizar a eficiência do passo de tempo enquanto preserva a alta ordem de precisão em regiões suaves. O trabalho estabelece uma base para a aplicação desses esquemas a uma gama mais ampla de fluxos, incluindo problemas de choque viscoso e interações choque-turbulência.