Fast Prediction of Aero-Thermal Performance During Hypersonic Vehicle Design

Resumo Técnico: Predição Rápida de Desempenho Aero-Térmico Durante o Design de Veículos Hipersônicos

Definição do Problema
O projeto de veículos hipersônicos requer a estimativa rápida de forças aerodinâmicas e fluxos de calor superficiais durante os estágios iniciais de desenvolvimento. Embora simulações de Dinâmica de Fluidos Computacional (CFD) de alta fidelidade ofereçam precisão detalhada, elas são computacionalmente proibitivas para as extensas iterações de design e otimizações de trajetória necessárias nas fases preliminares. Métodos tradicionais de nível de engenharia, como técnicas de inclinação superficial local (cunha tangente, Newtoniano modificado, Prandtl-Meyer), oferecem velocidade, mas frequentemente carecem de fidelidade na modelagem de efeitos viscosos. Especificamente, as abordagens padrão tipicamente assumem que a camada limite é totalmente laminar ou totalmente turbulenta. A suposição de um escoamento totalmente turbulento, embora conservadora, frequentemente leva à superpredição do arrasto de fricção superficial e do fluxo de calor superficial. Isso resulta em Sistemas de Proteção Térmica (TPS) superdimensionados, aumento da massa do veículo e maior consumo de combustível. Além disso, poucas ferramentas de engenharia rápidas incorporam atualmente os efeitos de transição de camada limite dentro de estruturas de inclinação local, apesar do impacto significativo dos escoamentos transicionais nas cargas aerotérmicas.

Metodologia
O artigo apresenta o HipeX, uma ferramenta de análise aerotermodinâmica leve projetada para preencher a lacuna entre métodos analíticos clássicos e CFD de alta fidelidade. O HipeX utiliza uma abordagem híbrida de nível de engenharia baseada em métodos de inclinação superficial local para estimar propriedades de escoamento invíscido e emprega um modelo de camada limite misto para abordar efeitos viscosos e transição.

• Modelagem de Escoamento Invíscido: O solver processa malhas superficiais não estruturadas (formato STL) e seleciona o modelo de escoamento apropriado com base na inclinação superficial local em relação ao vetor de velocidade incidente:

  • Teoria Newtoniana Modificada: Aplicada para choques descolados (corpos rombudos) e grandes ângulos de deflexão.
  • Relações de Choque Oblíquo: Usadas para choques aderidos e regiões de compressão.
  • Expansão de Prandtl-Meyer: Aplicada para regiões de expansão de escoamento, com uma salvaguarda que limita a aplicação a 50% do ângulo de deflexão máximo teórico para garantir estabilidade numérica.
  • Fórmula de Gaubeaud: Usada como um fallback empírico quando as soluções de Prandtl-Meyer falham.
  • As propriedades atmosféricas são determinadas via modelos padrão baseados na altitude de voo.

• Modelagem Viscosa e Térmica:

  • Compressibilidade: Corrigida usando o método da temperatura de referência de Eckert, aplicável tanto para regimes laminares quanto turbulentos.
  • Modelagem de Transição: Um número de Reynolds crítico ($Re_{cr}$) é calculado empiricamente como uma função do número de Mach para determinar o início da transição. O modelo emprega um fator de intermitência ($\gamma_{it}$) variando de 0 (laminar) a 1 (turbulento) para misturar os coeficientes de fricção superficial de fluxos laminares e turbulentos ao longo de uma região de transição definida.
  • Cálculo de Fluxo de Calor: O fluxo de calor superficial é derivado do coeficiente de fricção superficial usando a analogia de Reynolds e o número de Stanton. O modelo considera tanto a transferência de calor por convecção quanto a transferência de calor por radiação (assumindo uma emissividade de 0,4), calculando o fluxo de calor líquido e a carga térmica total integrando sobre a área molhada.

• Fluxo de Trabalho: A ferramenta aceita geometria 3D e condições de voo (Mach, altitude, Ângulo de Ataque, Deriva). Ela computa propriedades invíscidas, aplica correções viscosas, integra forças para determinar coeficientes aerodinâmicos e exporta os resultados em formato VTK. O código suporta execução paralela para análise de trajetória.

Principais Contribuições

1. Modelo de Transição Misto: A principal contribuição é a integração de um modelo de camada limite transicional em uma estrutura de inclinação local rápida. Isso permite que o solver estime a mudança de escoamento laminar para turbulento com base em números de Reynolds locais e números de Mach, em vez de depender de suposições binárias laminar/turbulento.
2. Eficiência Computacional: O HipeX é projetado para execução rápida, capaz de computar o desempenho de uma trajetória completa (aprox. 300 pontos) em uma estação de trabalho padrão em cerca de um minuto.
3. Validação contra Dados de Referência: A ferramenta é validada contra dados de CFD e referências experimentais para a configuração HEXAFLY-INT (um design de corpo de sustentação misto para cruzeiro em Mach 7–8).

Resultados
A validação foi realizada utilizando a configuração HEXAFLY-INT através de vários números de Mach (M=2 a M=7,5) e trajetórias de voo.

• Coeficientes Aerodinâmicos: Para números de Mach elevados (M=6 e M=7), as previsões do HipeX para sustentação ($C_L$), arrasto ($C_D$) e centro de pressão ($X_{cp}$) alinham-se estreitamente com os dados de CFD. A solução de transição mista geralmente fornece resultados mais precisos do que as suposições puramente laminares ou totalmente turbulentas. Em números de Mach mais baixos (M=4 e M=2), as discrepâncias aumentam, particularmente nos valores absolutos de $C_L$ e $C_D$, embora a razão sustentação-arrasto ($L/D$) permaneça razoavelmente precisa. A ferramenta identifica corretamente os limites dos métodos de inclinação local conforme o escoamento se afasta das condições fortemente hipersônicas.
• Cargas Térmicas: Em casos onde o escoamento é predominantemente laminar (alta altitude), a solução mista produz resultados semelhantes ao caso totalmente laminar. No entanto, em regimes transicionais, o modelo misto captura níveis intermediários de fluxo de calor que situam-se entre os extremos laminar e turbulento.
• Implicações para o TPS: O estudo demonstra que assumir escoamento totalmente turbulento leva à superpredição significativa do aquecimento, particularmente próximo ao bordo de ataque. Ao capturar a natureza laminar ou parcialmente turbulenta da camada limite ao longo de grande parte da trajetória, o HipeX fornece estimativas de fluxo de calor mais baixas e realistas. Isso sugere potencial para reduzir a massa do TPS sem comprometer as margens de segurança.
• Início da Transição: Comparações com o código BLITZ mostram tendências semelhantes na localização da transição, embora o HipeX tenda a prever um início de turbulência mais precoce. O modelo atualmente tem dificuldade em evitar a previsão de relaminarização em certas regiões de expansão.

Significância e Alegações
O artigo posiciona o HipeX como uma base robusta para avaliação aerotermodinâmica e otimização de design em estágio inicial. Sua significância reside em oferecer um "equilíbrio entre compensações" entre conservadorismo de design e eficiência de missão. Ao afastar-se da suposição excessivamente conservadora de escoamento totalmente turbulento, a ferramenta permite o dimensionamento mais preciso de sistemas de proteção térmica, potencialmente reduzindo a massa do veículo e o consumo de combustível.

Os autores reconhecem que o HipeX herda as limitações intrínsecas dos métodos de inclinação de superfície local, especificamente a incapacidade de prever interações de choque, zonas de separação e fortes efeitos tridimensionais. Adicionalmente, a determinação precisa do início da transição em condições reais de voo permanece desafiadora devido às dependências de rugosidade superficial e perturbações do escoamento livre. No entanto, a ferramenta obtém sucesso ao preencher a lacuna entre métodos analíticos clássicos e CFD de alta fidelidade, permitindo a avaliação rápida de trajetórias completas e grandes conjuntos de configurações geométricas essenciais para o design preliminar.