Study of the Effects of Artificial Dissipation and Other Numerical Parameters on Shock Wave Resolution

Este estudo investiga os efeitos de diferentes esquemas numéricos e parâmetros de malha na resolução de ondas de choque em escoamentos supersônicos invíscidos bidimensionais, demonstrando que a formulação AUSM+ proposta oferece maior robustez na eliminação de perturbações não físicas e apresenta boa concordância com resultados experimentais.

O essencial

Imagine que você é um engenheiro tentando projetar a cápsula de uma nave espacial que vai entrar na atmosfera da Terra em velocidade supersônica (mais rápida que o som). O maior desafio aqui é entender como o ar se comporta quando bate na frente da nave, criando uma onda de choque gigante, como a esteira de um barco, mas feita de ar superaquecido.

Este artigo é como um "laboratório de testes" onde os autores, Frederico e João, pegaram cinco ferramentas matemáticas diferentes (chamadas de esquemas numéricos) para simular esse cenário no computador e ver qual delas desenha a onda de choque com mais precisão.

Aqui está a explicação do que eles descobriram, usando analogias do dia a dia:

1. O Problema: O "Mapa" Distorcido

Para fazer a conta no computador, eles transformaram o espaço ao redor da nave em uma grade (como um papel milimetrado). O problema é que, perto da nave, essa grade precisa curvar muito para acompanhar o formato redondo.

• A Analogia: Imagine tentar desenhar um círculo perfeito em um papel de grade quadrado. Se você forçar o papel a curvar, as linhas da grade ficam tortas.
• O Erro: Quando as linhas da grade estão tortas, o computador começa a "alucinar". Ele vê ventos e mudanças de pressão onde não deveria haver nada (no ar livre, antes da onda de choque). É como se o GPS do carro mostrasse que você está atravessando um prédio quando você está na estrada. Isso é chamado de "perturbação não física".

2. As Ferramentas (Os Esquemas)

Os autores testaram cinco métodos diferentes para calcular o fluxo de ar:

• Beam e Warming: Um método clássico, mas que precisa de um "amortecedor" extra para não vibrar demais.
• Steger e Warming & van Leer: Métodos que tentam prever para onde o ar vai (para cima ou para baixo). O van Leer é um pouco mais refinado que o Steger.
• AUSM+ (duas versões): Métodos mais modernos. Os autores criaram duas versões diferentes de como aplicar essa fórmula no computador.

3. A Descoberta: O "Solução Mágica" (Subtração do Fluxo Livre)

Eles perceberam que, sem ajuda, quase todos os métodos criavam aquele "ruído" ou "alucinação" no ar antes da onda de choque.

• A Solução: Eles aplicaram uma técnica chamada "Subtração do Fluxo Livre".
• A Analogia: Pense em tentar ouvir uma música fraca em uma sala barulhenta. Em vez de tentar ouvir a música e ignorar o barulho, você grava o barulho da sala, toca-o de volta e inverte a fase (cancela o som). O que sobra é apenas a música.
• O Resultado: Ao subtrair matematicamente o estado do ar "calmo" (o fluxo livre) antes de calcular, o computador para de alucinar. O ruído desapareceu quase completamente (em mais de 6 ordens de magnitude!). Isso é como se o GPS parasse de mostrar prédios e mostrasse apenas a estrada.

4. O Dilema da Precisão vs. Estabilidade

Eles queriam usar métodos de "segunda ordem" (mais precisos, como uma foto em alta resolução) em vez de "primeira ordem" (menos precisos, como um desenho simples).

• O Problema: Quando tentaram usar a "alta resolução" sem ajuda, a simulação explodiu (divergiu) porque o computador ficou instável.
• A Solução 1 (Adição de Dissipação Artificial): Adicionar um "amortecedor" artificial.
  • O Efeito Colateral: Funcionou para estabilizar, mas empurrou a onda de choque para a frente, como se o amortecedor estivesse empurrando a onda. A posição da onda ficou errada.
• A Solução 2 (Limitador de Fluxo): Um "freio inteligente" que só age onde o ar está mudando bruscamente (na onda de choque) e deixa o resto fluir suavemente.
  • O Resultado: Este foi o vencedor. Ele estabilizou a simulação sem empurrar a onda de choque para o lugar errado. A posição da onda bateu perfeitamente com os dados experimentais reais.

5. O Grande Vencedor: AUSM+ (Versão 2)

Dentre todas as tentativas, uma versão específica do método AUSM+ (chamada de "Ap. 2" no texto) se destacou.

• Por que? Ela foi construída de uma forma que, naturalmente, não cria aqueles erros de "alucinação" no ar calmo, mesmo sem precisar de correções extras. É como se fosse um carro que já vem de fábrica com um sistema de suspensão perfeito para estradas tortas.

Conclusão Simples

O estudo nos ensina que, para simular o voo supersônico com precisão:

1. O computador pode "inventar" problemas se a grade de cálculo for muito curva.
2. A melhor maneira de consertar isso é "subtrair" o estado do ar calmo da equação (como cancelar o ruído de fundo).
3. Para ter precisão sem distorcer a posição da onda de choque, o uso de um "freio inteligente" (limitador) é melhor do que adicionar "amortecedores" artificiais brutos.

No fim, os resultados das melhores simulações bateram certinho com os dados de testes reais feitos em túneis de vento, provando que a matemática, quando bem ajustada, consegue prever a realidade com exatidão.

Resumo técnico

Título: Estudo dos Efeitos da Dissipação Artificial e Outros Parâmetros Numéricos na Resolução de Ondas de Choque

Autores: Frederico Bolsoni Oliveira e João Luiz F. Azevedo (Instituto Tecnológico de Aeronáutica e Instituto de Aeronáutica e Espaço, Brasil).

---

1. Problema Investigado

O artigo aborda os desafios numéricos na simulação de escoamentos supersônicos invíscidos ao redor de corpos rombudos (blunt bodies), especificamente focando na formação e resolução de ondas de choque descoladas (choques de proa ou bow shocks). O problema central reside na capacidade dos esquemas numéricos de:

• Capturar corretamente a estrutura da onda de choque sem introduzir distorções não físicas.
• Manter o estado de corrente livre (freestream) não perturbado a montante do choque.
• Lidar com as inconsistências geradas pela transformação de coordenadas em malhas curvilíneas, que podem induzir oscilações espúrias e erros na conservação da massa e momento.

O estudo é motivado pela necessidade de validar ferramentas de Dinâmica dos Fluidos Computacional (CFD) internas para o design de escudos térmicos e narizes de veículos de reentrada, onde a precisão na previsão do arrasto e da temperatura de parede é crítica.

2. Metodologia

Os autores utilizaram um código CFD interno para resolver as equações de Euler compressíveis bidimensionais em coordenadas curvilíneas gerais. O caso de teste consiste em um corpo circular 2D (diâmetro de 76,2 mm) imerso em um escoamento supersônico de ar frio (Mach 8,0).

Esquemas Numéricos Avaliados:
Foram comparados cinco esquemas de diferenças finitas:

1. Beam e Warming: Esquema implícito de fatoração aproximada com diferenças centrais de segunda ordem (requer dissipação artificial).
2. Steger e Warming: Esquema de divisão de vetor de fluxo (Flux Vector Splitting - FVS) original.
3. van Leer: Variante do FVS que garante diferenciabilidade contínua dos autovalores, melhorando a robustez.
4. AUSM+ (Advection Upstream Splitting Method Plus) - Ap.1: Reinterpretação do esquema de Liou para diferenças finitas, reconstruindo termos métricos usando a mesma função de divisão de pressão.
5. AUSM+ (Ap.2): Reinterpretação alternativa onde os termos métricos e o Jacobiano são reconstruídos via média simples dos nós adjacentes.

Técnicas de Controle Numérico:

• Dissipação Artificial: Uso do modelo não linear de Pulliam (diferenças de 2ª e 4ª ordem) para estabilizar esquemas centrados e de segunda ordem.
• Limitadores de Fluxo: Aplicação do limitador minmod para evitar oscilações em esquemas de segunda ordem e permitir transição suave para primeira ordem perto de descontinuidades.
• Subtração de Corrente Livre (Freestream Subtraction): Modificação das equações discretas para subtrair o estado de corrente livre dos fluxos, visando eliminar erros de conservação em malhas não uniformes.

3. Contribuições Principais

• Identificação da Fonte de Oscilações: O estudo demonstra que oscilações não físicas na corrente livre a montante do choque e deformações no perfil do choque são causadas principalmente por inconsistências no tratamento discreto dos termos métricos da transformação de coordenadas, e não apenas pela dissipação do esquema.
• Comparação de Reinterpretações do AUSM+: A proposta de duas formulações diferentes do esquema AUSM+ para o contexto de diferenças finitas. A formulação Ap.2 (com reconstrução por média simples dos termos métricos) mostrou-se superior, eliminando as oscilações de corrente livre mesmo sem a subtração de corrente livre.
• Análise da Dissipação Artificial: Demonstração de que a adição explícita de dissipação artificial (como no modelo de Pulliam) pode deslocar a posição da onda de choque a montante (upstream), introduzindo um viés numérico significativo em comparação com esquemas que utilizam limitadores de fluxo.
• Validação de Técnicas de Correção: Validação prática de que a subtração de corrente livre é uma solução eficaz e de baixo custo computacional para mitigar erros de conservação em malhas curvilíneas.

4. Resultados

• Comportamento de Corrente Livre: Esquemas de primeira ordem (Steger & Warming, van Leer, AUSM+ Ap.1) apresentaram oscilações significativas (até 4,4% de desvio no número de Mach) na região de corrente livre não perturbada. O esquema AUSM+ Ap.2 não apresentou essas oscilações.
• Efeito da Subtração de Corrente Livre: A aplicação da subtração de corrente livre reduziu as oscilações em pelo menos 6 ordens de magnitude para os esquemas afetados, restaurando o estado de corrente livre e eliminando deformações no choque.
• Estabilidade de Esquemas de Segunda Ordem: Esquemas de segunda ordem puros divergiram devido à falta de dissipação. A estabilização foi alcançada tanto pela adição de dissipação artificial quanto pelo uso de limitadores de fluxo.
• Posição do Choque:
  • Esquemas com dissipação artificial explícita tenderam a deslocar o choque para montante.
  • Esquemas com limitadores de fluxo (sem dissipação explícita adicional) apresentaram melhor concordância com dados numéricos e experimentais da literatura (Peery & Imlay, Lin, Holden et al.), posicionando o choque ligeiramente a jusante em relação aos esquemas com dissipação.
• Coeficiente de Pressão ($C_p$): Todos os esquemas estáveis mostraram boa concordância com dados experimentais para o coeficiente de pressão na superfície do cilindro, exceto nas regiões próximas à linha de centro (onde o escoamento é subsônico) e na linha sônica, onde o esquema de Beam e Warming e o Steger & Warming apresentaram desvios ou oscilações.

5. Significância e Conclusões

O trabalho fornece diretrizes práticas essenciais para a simulação de escoamentos supersônicos com choques fortes em malhas curvilíneas:

1. Recomendação de Correção: A subtração de corrente livre é altamente recomendada como uma correção simples e eficaz para eliminar oscilações espúrias em malhas não uniformes.
2. Escolha do Esquema AUSM+: A formulação AUSM+ Ap.2 é preferível em contextos de diferenças finitas, pois sua abordagem de reconstrução de termos métricos evita oscilações intrínsecas, mesmo na ausência de outras correções.
3. Estabilização de Alta Ordem: Para esquemas de segunda ordem, o uso de limitadores de fluxo é superior à adição explícita de dissipação artificial, pois evita o deslocamento indesejado da onda de choque e mantém maior precisão na sua localização, embora possa limitar a redução do resíduo da solução.

Em suma, o estudo destaca que a precisão na resolução de ondas de choque não depende apenas da escolha do esquema de discretização, mas criticamente do tratamento dos termos métricos da malha e das técnicas de correção de corrente livre aplicadas.