A Hybrid Finite-Volume Reconstruction Framework for Efficient High-Order Shock-Capturing on Unstructured Meshes

Resumo Técnico: Uma Estrutura de Reconstrução de Volumes Finitos Híbrida para Captura de Choque de Alta Ordem Eficiente em Malhas Não Estruturadas

Definição do Problema
A Dinâmica de Fluidos Computacional (CFD) para fluxos compressíveis enfrenta um compromisso persistente entre precisão e custo computacional, particularmente em malhas não estruturadas exigidas por geometrias industriais complexas. Métodos de baixa ordem são robustos e não oscilatórios próximos a choques, mas sofrem com dissipação numérica excessiva, falhando em resolver estruturas de fluxo de pequena escala em regiões suaves. Por outro lado, esquemas de alta ordem não oscilatórios, como o Weighted Essentially Non-Oscillatory (WENO) e suas variantes, oferecem precisão superior em regiões suaves, mas incorrem em um custo computacional significativo e podem apresentar dificuldades de robustez em fluxos dominados por choques fortes. As abordagens híbridas existentes frequentemente dependem de detecção a posteriori (por exemplo, o paradigma MOOD), que pode ser computacionalmente caro devido à recomputação iterativa de células "problemáticas", ou de detectores a priori que são ou muito sensíveis ou muito restritivos, levando a uma dissipação ou oscilações desnecessárias.

Metodologia
Os autores propõem uma estrutura de reconstrução híbrida de ordem arbitrária e multidimensional implementada no solver de volumes finitos de código aberto UCNS3D. A inovação central é uma estratégia de detecção a priori de etapa única que classifica o campo de fluxo local em três regiões distintas, atribuindo um método de reconstrução ideal para cada uma:

1. Regiões Suaves: É empregada uma reconstrução linear de alta ordem. Isso evita o custo computacional de esquemas não lineares onde eles não oferecem benefício de precisão e evita problemas de degeneração de ordem observados às vezes em esquemas WENO tradicionais em extremos suaves.
2. Regiões Fracamente Não Suaves: O esquema Compact Weighted Essentially Non-Oscillatory (CWENOZ) é aplicado. Ele lida com regiões com gradientes fortes ou turbulência em desenvolvimento, onde a precisão de alta ordem é necessária, mas as descontinuidades ainda não estão totalmente formadas.
3. Regiões Descontínuas: Um esquema MUSCL de segunda ordem (com limitadores de inclinação) é ativado próximo a choques reais e descontinuidades para garantir estabilidade e comportamento não oscilatório.

Componentes Principais:

• Detector de Admissibilidade Numérica (NAD) Redesenhado: Diferente da abordagem iterativa a posteriori do paradigma MOOD, este NAD utiliza um Princípio do Máximo Discreto (DMP) relaxado para classificar células em uma única etapa. Ele introduz parâmetros auxiliares ($\delta_m$ e $\delta_w$) para definir margens para as regiões MUSCL e CWENOZ, respectivamente. Isso permite o ajuste flexível para equilibrar robustez e resolução.
• Critério de Desativação: Para evitar o acionamento falso em regiões de campo distante suaves, um critério de desativação semelhante ao limitador de Venkatakrishnan é usado, garantindo que esquemas lineares de alta ordem permaneçam ativos onde o fluxo é genuinamente suave.
• Reconstrução: A estrutura utiliza um método de reconstrução de mínimos quadrados com um algoritmo Stencil-based Compact (SBC), suportando ordens polinomiais arbitrárias ($P3, P5, P7$, etc.) em malhas de elementos mistos (triângulos, quadriláteros, tetraedros, hexaedros).
• Integração de Fluxo e Tempo: O esquema utiliza solvers Riemann HLLC ou HLL para fluxos invariantes, mínimos quadrados restritos para gradientes viscosos e um integrador de tempo Runge-Kutta de Estabilidade Forte Preservada (SSP) de 4ª ordem.

Principais Contribuições

• Estratégia Híbrida: A integração de esquemas lineares, CWENOZ e MUSCL em uma estrutura unificada que minimiza o uso de reconstruções não lineares dispendiosas sem sacrificar a estabilidade.
• Eficiência: Ao restringir as reconstruções custosas de CWENOZ e MUSCL apenas às regiões necessárias, o método reduz significativamente o overhead computacional em comparação com esquemas não lineares uniformes de alta ordem.
• Flexibilidade: Os parâmetros do NAD permitem que a estrutura opere como um híbrido de três esquemas, um híbrido de dois esquemas (ex: Linear-CWENOZ ou MUSCL-CWENOZ), ou até mesmo uma abordagem de esquema único, adaptando-se à física do fluxo e aos requisitos de resolução específicos.
• Implementação: Implementação completa no solver de código aberto UCNS3D, suportando reconstruções de ordem arbitrária em malhas não estruturadas.

Resultados
A estrutura foi validada através de uma série abrangente de testes de referência 2D e 3D:

• Evolução de Vórtice 2D: Em fluxo suave, o método híbrido produziu resultados idênticos ao esquema linear de alta ordem, confirmando que não ativa erroneamente os limitadores não lineares. Alcançou a mesma precisão que o CWENOZ com um tempo de CPU significativamente menor (reduzindo o custo em quase metade em comparação ao CWENOZ).
• Problema de Shu-Osher: O método capturou com sucesso interações de choque-onda de entropia. O estudo destacou que a escolha do limitador de inclinação (Barth & Jespersen vs. MOGE) impacta significativamente o desempenho do detector, com o limitador de Barth & Jespersen fornecendo o melhor equilíbrio entre robustez e precisão para este setup híbrido específico.
• Reflexão de Mach Duplo & Degrau de Face Voltada para Frente: O método resolveu interações de choque e camadas de cisalhamento complexas. Embora simulações de ordem superior (P7) tenham mostrado um pouco mais de dissipação do que o CWENOZ puro na resolução de pequenos vórtices, isso foi atribuído à ativação do esquema MUSCL e pode ser mitigado relaxando os critérios híbridos (aumentando a região CWENOZ).
• Problema de Riemann 2D: O método híbrido demonstrou excelente robustez em várias configurações. A configuração padrão melhorou a eficiência computacional em 41,67% sobre o CWENOZ, mantendo capacidades de captura de choque comparáveis.
• Vórtice de Taylor-Green Viscoso Supersônico (3D): Em uma simulação de fluxo turbulento ($M=1.25, Re=1600$), o esquema híbrido de 6ª ordem resolveu uma gama mais ampla de números de onda e igualou soluções de referência estabelecidas (TENO e DG) melhor do que o esquema CWENOZ6 uniforme. Crucialmente, alcançou um ganho de velocidade de 2,5× em comparação com o esquema CWENOZ6 uniforme no mesmo hardware.
• Fluxo Mach 3 sobre um Cilindro: O método capturou o choque de proa e instabilidades de esteira de forma robusta, com o esquema linear de alta ordem dominando o domínio e os esquemas não lineares ativando-se apenas próximos aos choques. Isso resultou em uma redução de três vezes no custo computacional por iteração em comparação com uma discretização uniforme CWENO6.

Significância e Alegações
O artigo afirma que esta estrutura híbrida aproxima a precisão de alta ordem das simulações de CFD práticas e de escala industrial. Ao combinar o custo computacional reduzido de esquemas lineares em regiões suaves com a robustez de esquemas não lineares próximos a choques, o método oferece uma alternativa confiável aos esquemas não lineares uniformes de alta ordem. Os autores enfatizam que a estrutura mantém as ordens de precisão projetadas em regiões suaves enquanto garante soluções não oscilatórias perto de descontinuidades. A capacidade de alcançar ganhos de velocidade significativos (até 2,5× em turbulência 3D) sem comprometer a estabilidade ou a precisão sugere que esta abordagem é um caminho viável para simulações de alta fidelidade de fluxos compressíveis complexos em malhas não estruturadas.