Introdução à dinâmica de fluidos criogênicos em sistemas espaciais

Fluidos criogênicos – gases liquefeitos armazenados em temperaturas abaixo de -150 °C – são o sangue vital dos modernos sistemas de propulsão espacial e suporte de vida. Hidrogênio líquido (LH2) e oxigênio líquido (LOX) alimentam os principais motores de foguetes, como Raptor da SpaceX, BE-4 da Blue Origin e Sistema de Lançamento Espacial da NASA. Gerenciar esses fluidos em ambientes térmicos extremos e microgravidade apresenta desafios de engenharia únicos. Mudança de fase, estratificação, instabilidade geiser e fervura devem ser preditos com precisão para garantir a segurança da missão e eficiência de combustível.

Ferramentas de dinâmica de fluidos computacional (CFD) como COMSOL Multiphysics com o Módulo CFD oferecem uma estrutura robusta para simular o comportamento de fluidos criogênicos. Ao conectar fluxo de fluidos, transferência de calor, termodinâmica e mudança de fase em um único ambiente, os engenheiros podem testar projetos de tanque, linhas de alimentação e estratégias de isolamento antes da criação do hardware. Este artigo fornece um guia confiável e pronto para a produção para modelar a dinâmica de fluidos criogênicos em CFD COMSOL para aplicações espaciais.

Por que os fluidos criogênicos exigem simulação especializada

Ao contrário dos líquidos comuns, os criógenos apresentam fortes variações de propriedades com temperatura e pressão. O hidrogênio líquido, por exemplo, tem uma densidade de apenas cerca de 70 kg/m3 em seu ponto de ebulição e uma condutividade térmica que muda acentuadamente perto do ponto crítico. regimes de ebulição - nucleato, transição e ebulição de filme - pode coexistir em um único tanque devido à entrada de calor localizada do motor ou radiação solar. Na microgravidade, a tensão superficial domina, criando fluxos capilares que influenciam o posicionamento líquido e a transferência de calor.

Essas complexidades tornam as correlações empíricas pouco confiáveis para o projeto.A simulação multifísica deve resolver simultaneamente:

  • Transferência de calor conjugado: condução através de paredes de tanque, convecção no criógeno, e radiação de fontes externas.
  • Dinâmica de mudança de fase:] evaporação, condensação e efeitos de gás real perto da linha de saturação.
  • Regimes de fluxo bifásico: borbulhante, estratificado, anular e nebulosidade em tubos e tanques.
  • Interação da estrutura do fluido: deformação do tanque sob tensões térmicas criogênicas e acelerações de lançamento.

A abordagem multifísica da COMSOL permite que os engenheiros combinem estes fenómenos num único modelo em vez de se iterar entre ferramentas separadas.

Capacidades Principais do CFD COMSOL para Modelação Criogênica

Equações e Modelos de Turbulência Governantes

O Módulo CFD resolve as equações de Navier-Stokes para massa, momento e energia, juntamente com um modelo de turbulência adequado. Para tanques criogênicos, o modelo k-ω SST funciona bem porque resolve fluxos de parede próximos com precisão – crítico para prever transferência de calor através da parede do tanque. Em cenários de baixa gravidade, modelos laminados ou transitórios são, por vezes, preferidos porque os níveis de turbulência são baixos e a flutuabilidade é a mistura de unidades.

Mudança de Fase e Fluxo Multifásico

O COMSOL oferece duas abordagens principais para a mudança de fase:[
]Set de Nível e Campo de Fase Os métodos de Rastreiam explicitamente a interface líquido-vapor. Estes são ideais para problemas de superfície livre, tais como perda de tanques ou impacto de jato. O Modelo de Mixtura[] é mais adequado para fluxos bifásicos dispersos (bolhas ou gotas). Para ferver, a interface Alteração de Fase Termal pode alavancar as funções definidas pelo usuário para taxas de transferência de massa com base em gradientes de temperatura.

Termodinâmica não-ideal

O recurso de Gás Real da COMSOL permite importar equações fundamentais de estado, como Peng-Robinson ou Span-Wagner, para previsões de propriedades altamente precisas. Para o hidrogénio e o oxigénio, estas equações são essenciais perto do ponto crítico, onde as suposições de gás ideal dão erros superiores a 30 %.

Bibliotecas de Materiais Personalizados

Você pode definir tabelas de propriedades de NIST REFPROP dados diretamente em COMSOL, garantindo densidade, calor específico, viscosidade e condutividade térmica são funções de temperatura e pressão. Esta é uma vantagem chave sobre correlações de propriedade genéricas.

Processo de modelagem passo a passo em COMSOL

1. Definição de Geometria

Comece com uma geometria axissimétrica 2D ou completa 3D do tanque criogénico, linha de alimentação ou cabeça de propulsor. Para um tanque esférico típico, inclua a espessura da parede, as portas de entrada/saída, os desfibriladores anti-vórtices e as camadas de isolamento. Simplifique filetes não críticos e furos de parafusos para manter a malha manejável.

2. Atribuindo Propriedades do Material

Crie um material personalizado para o criogénio utilizando tabelas NIST. Por exemplo, para LH2 a 20 K a 30 K, importe densidade em função da temperatura e pressão e defina a viscosidade dinâmica utilizando o modelo Lennard-Jones. Para a parede do tanque (por exemplo, alumínio 2219 ou aço inoxidável 316L), defina condutividade térmica, densidade e coeficiente de expansão térmica. Adicione um material de isolamento como isolamento multicamadas (MLI) com uma condutividade térmica eficaz que depende da pressão de vácuo.

3. Interfaces Físicas e Acoplamento

Adicione as seguintes interfaces do Assistente de Modelos:

  • Fluxo de Lamas ou Fluxo de Turbulentas (k-ω SST)[] para movimento de fluidos.
  • Transferência de calor em fluidos e Transferência de calor em sólidos para transferência de calor conjugado.
  • Campo de Fase (ou Conjunto de Nível) para o rastreamento de interface, se investigar dinâmica bifásica.
  • Fluxo não-isotérmico multifísico de acoplamento para ter em conta a densidade e flutuabilidade dependentes da temperatura.

Para modelos de vaporização, adicione uma Fonte de calor na interface igual ao calor latente de vaporização vezes a taxa de evaporação. Isto pode ser implementado através de uma expressão definida pelo usuário que ativa quando a temperatura líquida atinge o ponto de saturação na pressão local.

4. Condições de limite

Condições típicas para um tanque de propelente durante uma fase costeira:

  • Superfície exterior da parede do tanque: fluxo de calor da radiação solar (por exemplo, 1,3 kW/m2 a 1 UA mais infravermelho da Terra).
  • Superfície interna da parede do tanque: sem deslizamento, fluxo de calor conjugado.
  • Interface de vapor líquido: continuidade da temperatura, equilíbrio de massa com taxa de mudança de fase.
  • Saída (linha de alimentação): limite de pressão igual à pressão de ullagem do tanque menos cabeça.
  • Entrada (para enchimento): condição de limite do caudal mássico.

Para simulações de microgravidade, use uma Força de Volume para tensão superficial (baseada no modelo de força de superfície contínuo) e ajuste aceleração para zero ou um pequeno valor de micro-g (10 a 4 g).

5. Estratégia de Meshing

A malha de alta qualidade é fundamental para simulações criogênicas. Use uma malha de camada limite (pelo menos 5-10 camadas de prisma) na parede do tanque para capturar gradientes de temperatura íngremes. Na região de interface bifásico, refine elementos com um tamanho de cerca de 1 mm (para um tanque de 1 m de diâmetro) para resolver a curvatura da interface. Use uma malha triangular livre [ para modelos axissimétricos 2D ou Tetraédrio livre] com extrusão de camada limite para 3D.

Realize um estudo de convergência de malhas: execute o modelo com refinamento duplo e compare os perfis de temperatura e as posições de interface. Selecione a malha mais grossa que dá resultados dentro de 2 % da solução refinada.

6. Solucione configurações e passo de tempo

Os problemas criogênicos são frequentemente transitórios. Use uma fórmula de diferenciação para trás ] resolvedor de ordem 1 ou 2 com um passo de tempo inicial de 0,01 s. Para física acoplada, habilite a abordagem Solvedor de combinação completa com um fator de amortecimento de 0,8 para melhorar a convergência. Se o modelo correr lentamente, mude para uma abordagem Segregada[[, resolvendo primeiro o fluxo, então a transferência de calor, e então o campo de fase em cada passo.

Monitorar a convergência verificando os resíduos (inferiores a 10 a 4) e as quantidades conservadas (balanço de massa dentro de 0,1 %).

7. Pós-Processo e Análise

Saídas de chave a extrair:

  • Contornos de temperatura nas paredes do tanque e no domínio dos fluidos.
  • Taxa de desactivação (kg/s) integrada na interface.
  • Evolução da pressão no volume da ullagem (importante para o projeto estrutural do tanque).
  • Linhas de estria ou linhas de caminho mostrando zonas de recirculação que podem causar estratificação térmica.
  • Posição de interface em função do tempo — crítico para avaliar dispositivos de aquisição de líquidos.

As ferramentas do COMSOL Valores Derivados podem calcular automaticamente vazamento total de calor, massa de vapor gerada e fração vazia.

Estudo de caso: Simulando a Fervagem em um tanque de propelente criogênico

Considere um tanque LH2 esférico (rádio 1,5 m) com paredes de alumínio de 1 mm de espessura e 5 cm de MLI. O tanque está 50 % cheio, a uma pressão de 3 bar, com o fluido à temperatura de saturação (cerca de 23 K). O fluxo de calor externo é de 150 W/m2 do ambiente térmico da nave espacial. O objetivo é prever a taxa de fervura durante um período de 24 horas de costa.

Após a configuração do modelo descrito acima, a simulação mostra que o LMI reduz o vazamento de calor para cerca de 8 W/m2. A taxa de descamação estabiliza-se a 0,25 g/s após um transiente inicial de 200 s. As varreduras paramétricas revelam que a duplicação da espessura do LMI reduz o descamação em 45 %, mas a penalidade de peso pode não ser aceitável. O modelo também prevê um gradiente de temperatura de 2 K entre o fundo e o topo do tanque, o que ajuda a identificar locais onde a condensação de vapor poderia ocorrer, causando potencialmente decaimento da pressão do tanque e dificuldades de reignição.

Tais simulações informam diretamente o projeto da missão. Por exemplo, podem ajudar a dimensionar um sistema de ventilação termodinâmica (TVS) que ativamente esfria o tanque para evitar sobrepressurização.

Aplicações através do ciclo de vida da nave espacial

Projeto do sistema de alimentação de propulsores

As linhas de alimentação criogénicas devem transferir líquidos sem cavitação ou fluxo bifásico que possa matar o motor. Modelos COMSOL da linha de alimentação – incluindo curvas, junções e válvulas – podem prever quedas de pressão e detectar regiões onde a ebulição pode ocorrer devido a pontos quentes locais. O Modelo de Mixtura[] com velocidade de deslizamento captura separação de fase em alimentações de microgravidade, ajudando engenheiros a projetar dispositivos de aquisição de líquidos (LADs) como canais de tela ou palhetas.

Validação do sistema de proteção térmica

Materiais de isolamento, incluindo aerogéis e cobertores MLI, podem ser modelados como meios porosos ou sólidos em camadas com condutividades térmicas eficazes. Ao executar uma simulação de estado estacionário do tanque com aquecimento de pior caso, os engenheiros verificam que o isolamento mantém o criogênio abaixo do seu ponto de ebulição durante toda a missão.

Iniciar as operações de Chill-Down e Preencher

O enchimento de um tanque com criogênio das condições ambientais produz uma ebulição violenta e choque térmico. As simulações transitórias em COMSOL podem avaliar o tempo necessário para arrefecer as paredes do tanque para temperaturas seguras (frequentemente abaixo de 100 K) antes do enchimento principal. Estes modelos também ajudam a otimizar a taxa de enchimento para evitar a formação excessiva de vapor que pode danificar válvulas.

Gestão de Propelentes em Voo

A descamação em microgravidade provoca movimento líquido que pode desestabilizar uma nave espacial. As simulações de estrutura fluida em conjunto com a COMSOL permitem que os engenheiros projetem desorientações, diafragmas ou tanques de tensão superficial que mantêm o líquido estabilizado. O método Phase Field captura a interface em movimento durante as queimaduras do motor ou os propulsores de atitude, fornecendo forças nas paredes do tanque para análise estrutural.

Correlações de teste de solo

Antes de se poder confiar numa simulação para o voo, esta deve ser validada em função dos dados de ensaio em terra. O módulo de estimação de parâmetros da COMSOL pode ajustar parâmetros incertos (por exemplo, resistência ao contacto, factores de visão de radiação) para corresponder às leituras de termopar de um teste de subescala. Este modelo calibrado torna-se então o duplo digital para o tanque de voo.

Recursos externos para temas avançados

Para engenheiros que buscam aprofundar seu conhecimento criogênico de CFD, recomendam-se os seguintes recursos:

Instruções futuras: Gêmeos digitais e aprendizagem de máquina

A próxima fronteira na modelagem de fluidos criogênicos é a criação de gêmeos digitais em tempo real. Modelos de ordem reduzida (ROMs) podem ser construídos a partir de varreduras paramétricas COMSOL, em seguida, incorporados em software de voo de veículos para estimativa de estado. Por exemplo, uma ROM que prevê taxa de fervura em função de vazamento de calor, nível de enchimento e aceleração pode funcionar em milissegundos, permitindo a medição de propelente a bordo.

O aprendizado de máquina também está entrando em campo. As redes neurais treinadas em bancos de dados de simulação COMSOL podem prever colapso de pressão após manobras de reorientação mais rápidas do que um modelo multifísico completo. Essas abordagens híbridas mantêm a precisão ao cortar o custo computacional por ordens de magnitude.

O COMSOL LiveLink para MATLAB torna simples a interface de dados de simulação com tais ferramentas de IA, permitindo aos engenheiros melhorar iterativamente as suas ROMs com novos dados de alta fidelidade à medida que a missão avança.

Conclusão

Modelar a dinâmica de fluidos criogênicos no COMSOL CFD proporciona aos engenheiros espaciais um ambiente poderoso e unificado para enfrentar os desafios multifísicos dos fluxos de baixa temperatura. Do projeto de conceito inicial até operações de voo, essas simulações permitem melhores escolhas de isolamento, sistemas de alimentação mais eficientes e gerenciamento de propulsores mais seguros. Ao alavancar equações de gás real, modelos avançados de mudança de fase e validação completa contra dados de teste, os engenheiros podem reduzir o risco de desenvolvimento e reduzir os tempos de ciclo.

Os exemplos e fluxos de trabalho apresentados aqui oferecem um ponto de partida para a construção de modelos de CFD criogênicos confiáveis. À medida que a exploração espacial avança para um espaço mais profundo e para missões mais longas, o gerenciamento preciso de fluidos criogênicos torna-se ainda mais crítico. O desenvolvimento contínuo da COMSOL, especialmente em fluxo multifásico e acoplamento térmico, posiciona-a como uma ferramenta essencial na caixa de ferramentas do engenheiro aeroespacial.