Introducción a la dinámica de fluidos criogénicos en los sistemas espaciales

Los fluidos criogénicos — gases licuados almacenados a temperaturas inferiores a -150 °C— son el sistema de propulsión espacial y soporte vital modernos. El hidrógeno líquido (LH2) y el oxígeno líquido (LOX) potencian los motores principales de cohetes como el Raptor de SpaceX, el BE-4 de Blue Origin y el sistema de lanzamiento espacial de la NASA.

Las herramientas de dinámica de fluidos computacionales (CFD) como COMSOL Multiphysics con el módulo CFD ofrecen un marco robusto para simular el comportamiento de fluido criogénico. Al acoplamiento de flujo de fluidos, transferencia de calor, termodinámica y cambio de fase en un solo entorno, los ingenieros pueden probar diseños de tanques, líneas de alimentación y estrategias de aislamiento antes de un modelo dinámico de artículo.

Por qué los fluidos criogénicos exigen simulación especializada

A diferencia de líquidos ordinarios, los criogénicos presentan fuertes variaciones de propiedades con temperatura y presión. El hidrógeno líquido, por ejemplo, tiene una densidad de sólo 70 kg/m3 en su punto de ebullición y una conductividad térmica que cambia bruscamente cerca del punto crítico. Los regímenes de la bobinación —núcleo, transición y ebullición de película— pueden coexistir en un solo tanque debido a la entrada de calor localizada del motor o radiación solar.

Estas complejidades hacen que las correlaciones empíricas no sean confiables para el diseño. La simulación multifísica debe resolver simultáneamente:

  • Conjugar la transferencia de calor: conducción a través de las paredes del tanque, la convección en el criogénico y la radiación de fuentes externas.
  • ]Modificación del cambio de fase: evaporación, condensación y efectos de gases reales cerca de la línea de saturación.
  • Dos regímenes de flujo: flujos de bublica, estratificada, anular y de neblina en tuberías y tanques.
  • Interacciones de estructura fluida: deformación de tanques bajo tensiones térmicas criogénicas y aceleraciónes de lanzamiento.

El enfoque multifísico deCOMSOL permite a los ingenieros combinar estos fenómenos en un solo modelo en lugar de iterar entre herramientas separadas.

Capacidades básicas de la CFD COMSOL para la modelación criogénica

Modelos de Ecuaciones y Turbulencias que rigen

El módulo CFD resuelve las ecuaciones de Navier‐Stokes para masa, impulso y energía, junto con un modelo de turbulencia adecuado. Para los tanques criogénicos, el modelo k‐ω SST funciona bien porque resuelve flujos de casi paredes con precisión, críticos para predecir la transferencia de calor a través de la pared del tanque.

Cambio de fase y flujo multifase

COMSOL ofrece dos enfoques principales para el cambio de fase:
Conjunto de distancia y Phase Field métodos de seguimiento de la interfaz de vacío líquido explícitamente. Estos son ideales para problemas de superficie libre como el apalancamiento de tanques o la dispersión de chorro de fase [LT:5]

Termodinámica no intradeal

La función de gas real le permite importar ecuaciones fundamentales de estado como Peng‐Robinson o Span‐Wagner para predicciones de propiedades altamente precisas. Para hidrógeno y oxígeno, estas ecuaciones son esenciales cerca del punto crítico, donde las suposiciones ideales-gas dan errores superiores al 30 %.

Bibliotecas de material personalizado

Puede definir tablas de propiedades de NIST REFPROP data directamente en COMSOL, asegurando densidad, calor específico, viscosidad y conductividad térmica son funciones de temperatura y presión. Esta es una ventaja clave sobre las correlaciones de propiedades genéricas.

Proceso de modelado de paso a paso en COMSOL

1. Definición de geometría

Comience con una geometría 3D axisymmétrica o completa del tanque criogénico, la línea de alimentación o la cabeza de empuje. Para un tanque de propulsión esférica típico, incluya el espesor de la pared, los puertos de entrada/outlet, las bafles anti-vortex y las capas de aislamiento. Simplifique los filetes no críticos y los agujeros de perno para mantener la meshing manejable.

2. Asignación de propiedades materiales

Crear un material personalizado para el criogénico utilizando tablas NIST. Por ejemplo, para LH2 a 20 K a 30 K, densidad de importación como función de temperatura y presión, y establecer la viscosidad dinámica utilizando el modelo Lennard-Jones. Para la pared del tanque (por ejemplo, aluminio 2219 o acero inoxidable 316L), definir conductividad térmica, densidad y coeficiente de expansión térmica.

3. Interfaces Físicas y Coupling

Agregue las siguientes interfaces del Asistente Modelo:

  • Flow Laminar] o Flujo turbulento (k-ω SST) para el movimiento fluido.
  • ] Transferencia de calor en fluidos y ]] Transferencia de calor en sólidos para la transferencia de calor conjugada.
  • Phase Field] (o Level Set) para el seguimiento de la interfaz, si se investiga la dinámica de dos fases.
  • N-Isothermal Flow acoplamiento multifísico para dar cuenta de la densidad y la flotabilidad dependientes de la temperatura.

Para los modelos de hervidor, agregue Fuente de calor en la interfaz igual al calor latente de los tiempos de evaporación, que se puede aplicar mediante una expresión definida por el usuario que se activa cuando la temperatura líquida alcanza el punto de saturación a la presión local.

4. Condiciones de los límites de los límites

Condiciones típicas para un tanque de propulsión durante una fase costera:

  • Superficie exterior de la pared de tanque: flujo de calor de la radiación solar (por ejemplo, 1,3 kW/m2 a 1 UA más infrarrojos terrestres).
  • Superficie interior de la pared del tanque: sin-slip, flujo de calor conjugado.
  • Interfaz líquido-vapor: continuidad de la temperatura, equilibrio de masa con velocidad de cambio de fase.
  • Outlet (línea alimentada): límite de presión igual a la presión de depósito de menos cabeza.
  • Entrada (para relleno): estado de límite de flujo de masa.

Para simulaciones de microgravedad, utilice una Fuerza de volumen] para la tensión superficial (basada en el modelo de fuerza superficial continua) y establece aceleración a cero o a un pequeño valor de microg (10−4 g).

5. Estrategia de lucha contra la pobreza

El corte de alta calidad es crítico para simulaciones criogénicas. Usa una malla de capa de límite (al menos 5-10 capas de prisma) en la pared del tanque para capturar gradientes de temperatura empinada. En la región de interfaz de dos fases, refina elementos a un tamaño de aproximadamente 1 mm (para un tanque de 1 m de diámetro) para resolver la curvatura de la interfaz.

Realizar un estudio de convergencia de malla: ejecutar el modelo con refinamiento doble y comparar perfiles de temperatura y posiciones de interfaz. Seleccione la malla más gruesa que da resultados dentro del 2 % de la solución refinada.

6. Ajustes de Solver y Pasos de Tiempo

Los problemas de calor son a menudo transitorios. Usar un Fórmula de diferenciación de fondo (BDF)] solucionador del orden 1 o 2 con un paso inicial de 0.01 s. Para la física acoplada, active el Totalmente acoplado solucionador con un factor de amortiguación de 0.8 para mejorar la convergencia[LT]

Monitor convergencia comprobando los residuos (abajo 10-4) y las cantidades conservadas (con balance de masa dentro del 0,1 %).

7. Proceso y análisis post-

Principales productos para extraer:

  • Contornos de la temperatura] en las paredes del tanque y el dominio del fluido.
  • Tasa de boil-off (kg/s) integrada sobre la interfaz.
  • Evolución de la presión] en el volumen de ullage (importante para el diseño estructural del tanque).
  • Estadísticas o líneas de ruta mostrando zonas de recirculación que pueden causar estratificación térmica.
  • Posición de la relación como función del tiempo, crítica para la evaluación de los dispositivos de adquisición de líquidos.

Los valores escritos Las herramientas de COMSOL pueden calcular automáticamente la fuga total de calor, la masa de vapor generada y la fracción de vacío.

Estudio de caso: Simulación de boil‐Off en un tanque de propellante criogénico

Considere un tanque LH2 esférico (radius 1,5 m) con paredes de aluminio de 1 mm de espesor y 5 cm de MLI. El tanque está lleno 50 %, a una presión de 3 bar, con el fluido a temperatura de saturación (aproximadamente 23 K). El flujo de calor externo es de 150 W/m2 del entorno térmico de la nave espacial. El objetivo es predecir la tasa de caldera durante un período de 24 horas.

Después de configurar el modelo como se describe anteriormente, la simulación muestra que el MLI reduce la fuga de calor a alrededor de 8 W/m2. La tasa de caldera se estabiliza a 0.25 g/s después de un transitorio inicial de 200 s. Los barridos paramétricos revelan que duplicar el espesor del MLI reduce el calentado en un 45 %, pero la pena de peso puede no ser aceptable.

Estas simulaciones informan directamente el diseño de la misión. Por ejemplo, pueden ayudar a tamaño un sistema de ventilación termodinámica (TVS) que enfría activamente el tanque para prevenir la sobrepresión.

Aplicaciones en todo el ciclo de vida de la nave espacial

Diseño de sistema de alimentación

Las líneas de alimentación criogénicas deben transferir líquido sin cavitación o flujo de dos fases que podrían morir de hambre. Los modelos COMSOL de la línea de alimentación, incluyendo curvas, uniones y válvulas, pueden predecir caídas de presión y detectar regiones donde se puede producir ebullición debido a puntos calientes locales. Modelo de la micrograbación]

Validación del sistema de protección térmica

Los materiales de aislamiento, incluyendo aerogeles y mantas MLI, pueden ser modelados como sólidos porosos o estucados con conductividades térmicas efectivas. Al ejecutar una simulación de estado estable del tanque con la calefacción de peor maletín, los ingenieros verifican que el aislamiento mantiene el criogen debajo de su punto de ebullición en toda la misión.

Lanzamiento de las operaciones de doblado y de relleno

El llenado de un tanque con criogénico de las condiciones ambientales produce cocción violenta y shock térmico. simulaciones transitorias en COMSOL pueden evaluar el tiempo necesario para enfriar las paredes del tanque a temperaturas seguras (a menudo por debajo de 100 K) antes del relleno principal. Estos modelos también ayudan a optimizar la tasa de llenado para evitar la formación excesiva de vapor que puede dañar las válvulas.

Gestión de Propellant In‐Flight

El cierre en microgravedad causa movimiento líquido que puede desestabilizar una nave espacial. Simulación de estructura de fluidos en COMSOL permite a los ingenieros diseñar baffles, diafragmas o tanques de tensión superficial que mantienen el líquido asentado. El método Phase Field captura la interfaz de movimiento durante las quemaduras del motor o los propulsores de actitud, proporcionando fuerzas estructurales.

Correlaciones de prueba de tierra

Antes de que se pueda confiar en una simulación para el vuelo, debe validarse contra los datos de prueba de tierra. El módulo de estimación de parámetros de COMSOL puede ajustar parámetros inciertos (por ejemplo, resistencia al contacto, factores de vista de radiación) para que coincidan con las lecturas de termopares de una prueba de subescala.

Recursos externos para temas avanzados

Para los ingenieros que buscan profundizar su conocimiento criogénico de la CFD, se recomiendan los siguientes recursos:

Instrucciones futuras: Gemelos digitales y aprendizaje automático

La siguiente frontera en el modelado de fluidos criogénicos es la creación de gemelos digitales en tiempo real. Los modelos de orden reducido (ROM) se pueden construir desde barridos paramétricos COMSOL, luego incrustados en software de vuelo de vehículos para la estimación del estado. Por ejemplo, un ROM que predice la tasa de desprendimiento como función de fuga de calor, nivel de llenado y aceleración puede funcionar en milisegundos, permitiendo el gauging propuls de propulsión.

Las redes neuronales entrenadas en bases de datos de simulación COMSOL pueden predecir el colapso de presión después de maniobras de reorientación más rápido que un modelo multifísico completo. Estos enfoques híbridos mantienen la precisión al reducir el costo computacional por órdenes de magnitud.

COMSOL LiveLink for MATLAB hace que sea más sencillo interactuar datos de simulación con tales herramientas de inteligencia artificial, permitiendo a los ingenieros mejorar iterativamente sus ROMs con datos de alta fidelidad frescos a medida que avanza la misión.

Conclusión

La modelación de dinámicas de fluido criogénico en COMSOL CFD proporciona a los ingenieros espaciales un entorno potente y unificado para abordar los desafíos multifísicos de flujos de baja temperatura. Desde el diseño inicial del concepto a través de operaciones de vuelo, estas simulaciones permiten mejores opciones de aislamiento, sistemas de alimentación más eficientes y una gestión de propulsión más segura.

Los ejemplos y flujos de trabajo aquí presentados ofrecen un punto de partida para la construcción de modelos criogénicos fiables de CFD. A medida que la exploración espacial empuja hacia un espacio más profundo y misiones más largas, la gestión precisa de fluido criogénico se vuelve aún más crítica. El desarrollo continuo de COMSOL, especialmente en flujo multifase y acoplamiento térmico, lo posiciona como una herramienta esencial en el cajón de herramientas del ingeniero aeroespacial.