Invariant domain preserving limiting of time explicit and time implicit discretizations for systems of conservation laws

Resumo Técnico: Limitação de Domínios Invariantes Preservadores para Discretizações Temporais Explícitas e Implícitas

Enunciação do Problema
O artigo aborda o desafio de projetar aproximações numéricas de alta ordem (AO) para sistemas não lineares hiperbólicos de leis de conservação que preservem domínios invariantes (DI). Embora os métodos de AO capturem com precisão soluções suaves, eles frequentemente falham em satisfazer restrições físicas no nível discreto, como a positividade da densidade e da energia interna ou o cumprimento de desigualdades de entropia. Por outro lado, esquemas de baixa ordem (BO) (por exemplo, Godunov, Lax-Friedrichs) são conhecidos por preservar esses domínios invariantes, mas sofrem de difusão numérica excessiva. O objetivo é construir uma técnica de limitação que restrinja a solução de AO para permanecer dentro do conjunto convexo de estados admissíveis $D$ (por exemplo, densidade positiva, energia interna positiva, limites de entropia) sem sacrificar a conservação ou a precisão de alta ordem em regiões suaves.

Metodologia
Os autores propõem um quadro geral de limitação aplicável a qualquer discretização espacial conservativa (Volume Finito, Galerkin Descontínuo) e a uma ampla gama de esquemas de integração temporal, incluindo Runge-Kutta (RK) explícitos e implícitos, métodos multietapa e integração temporal Galerkin Descontínua (DG).

1. Conceito Central: O método limita a solução de AO em torno de uma solução de BO que é garantida como preservadora de DI. Ele generaliza o quadro de Transporte Corrigido por Fluxo (FCT) ao limitar fluxos antidifusivos (a diferença entre os resíduos de AO e BO), mas define os coeficientes de limitação para expressar a solução final como uma combinação convexa de quantidades preservadoras de DI, semelhante ao quadro de limitação convexa.
2. Coeficientes de Limitação: Diferentemente do FCT padrão, que impõe limites aos componentes da solução, esta abordagem define coeficientes $l_{\kappa,\kappa'}^{i,j}$ de modo que a solução limitada $u_h^{n+1}$ seja uma combinação convexa da solução de BO e da atualização de AO. Os coeficientes são determinados encontrando os maiores fatores de escala $k$ tais que os estados intermediários permaneçam dentro do domínio invariante $D$.
3. Procedimento Iterativo: O limitador é incorporado em um algoritmo iterativo. Começando pela solução de BO, o método refina iterativamente a solução limitada ao reaplicar o limitador à diferença entre a solução de AO e a solução limitada atual. Este processo continua até que uma tolerância seja atendida ou um número máximo de iterações seja alcançado, visando recuperar a solução de AO o mais próximo possível, mantendo as propriedades de preservação de DI.
4. Aceleração: Um fator de aceleração heurístico $\beta \geq 1$ é introduzido para escalar os coeficientes de limitação, acelerando a convergência do processo iterativo em direção à solução de AO.
5. Adaptação à Integração Temporal: O quadro adapta-se a vários esquemas temporais definindo um esquema de BO de passo único (Euler-Forward para explícito, Euler-Backward para implícito) para calcular os resíduos necessários e os fluxos antidifusivos, evitando a necessidade de recalcular soluções de BO em cada estágio de métodos multietapa ou multietapa.
6. Discretização Espacial: Derivações específicas para fluxos antidifusivos são fornecidas para esquemas de Volume Finito (VF), esquemas DG gerais e o Método de Elementos Espectrais Galerkin Descontínuos (DGSEM). Para DGSEM, o método combina limitação baseada em médias de célula com um limitador de escala linear para garantir que os valores nodais também satisfaçam as restrições do domínio invariante.

Principais Contribuições

• Generalização para Sistemas e Esquemas Implícitos: O método estende conceitos de limitação convexa para sistemas de leis de conservação e, crucialmente, aplica-se a discretizações temporais implícitas, um domínio onde abordagens anteriores de limitação convexa estavam restritas a esquemas explícitos.
• Limitação Baseada em Fluxo sem Estados "Bar": A abordagem define coeficientes de limitação para formar combinações convexas de soluções de BO e AO sem depender de "estados bar" (estados intermediários usados em algumas formulações de limitação convexa). Isso permite o uso de qualquer fluxo numérico consistente e conservativo no esquema de AO.
• Refinamento Iterativo: A introdução de um procedimento iterativo permite que a solução limitada convirja para a solução de AO, recuperando a precisão em regiões suaves enquanto preserva estritamente os domínios invariantes.
• Heurística de Aceleração: Um simples fator de escala $\beta$ é proposto para reduzir significativamente o número de iterações necessárias para a convergência.
• Ampla Aplicabilidade: A derivação é apresentada para discretizações conservativas gerais, com implementações específicas detalhadas para VF, DG e DGSEM, acopladas a integradores temporais RK explícitos/implícitos, multietapa e DG.

Resultados Numéricos
O artigo valida o método através de experimentos numéricos em problemas 1D e 2D envolvendo equações hiperbólicas escalares e as equações de Euler compressíveis:

• Problemas de Riemann (Ondas de Choque Fortes): Para as equações de Euler compressíveis, o esquema de AO sem limitação falha devido a pressão/densidade negativas. O limitador proposto resolve com sucesso ondas de choque fortes de forma nítida, mantendo a positividade e os limites de entropia tanto para passos temporais explícitos (ERK3) quanto implícitos (DIRK33), mesmo em números de CFL grandes (até 10).
• Instabilidade de Kelvin-Helmholtz: Em simulações 2D com grandes passos temporais, o esquema de BO amortecia excessivamente as pequenas escalas. O esquema limitado recupera a alta resolução do esquema de AO, preservando estruturas de pequena escala.
• Difração de Choque e Escoamento em Degrau: O limitador demonstra robustez em geometrias 2D complexas, mantendo a admissibilidade física onde o esquema de AO ilimitado falha. Ele alcança resolução comparável a esquemas explícitos com pequenos passos temporais, mas com a estabilidade de métodos implícitos.
• Equações Escalares: Testes em transporte linear, equação de Burgers e o problema KPP (fluxo não convexo) mostram que o limitador remove efetivamente oscilações espúrias e impõe o princípio do máximo e desigualdades de entropia.
• Convergência: O limitador iterativo com o fator de aceleração $\beta=2$ é mostrado a convergir para a solução de AO significativamente mais rápido do que a abordagem padrão ($\beta=1$).

Significância e Alegações
Os autores afirmam que este trabalho fornece um limitador bastante geral capaz de lidar com qualquer lei de conservação hiperbólica em malhas gerais com uma ampla gama de esquemas de integração temporal. A principal significância reside em preencher a lacuna entre a robustez de esquemas implícitos de baixa ordem e a precisão de métodos de alta ordem. Ao garantir que o esquema numérico preserve domínios invariantes no nível discreto, o método permite o uso de grandes passos temporais em simulações implícitas sem o risco de produzir soluções não físicas (por exemplo, densidade negativa). O artigo enfatiza que o limitador mantém a conservação e pode ser aplicado iterativamente para refinar a solução, oferecendo uma ferramenta prática para simulações robustas de alta ordem de problemas complexos de dinâmica de fluidos. O trabalho conclui observando que extensões futuras considerarão passos temporais Euler-Backward linearizados para problemas de estado estacionário.