Introduction à la dynamique des fluides cryogéniques dans les systèmes spatiaux

Les fluides cryogéniques – gaz liquéfiés stockés à des températures inférieures à -150 °C – sont le noyau vital des systèmes modernes de propulsion spatiale et de survie. L'hydrogène liquide (LH2) et l'oxygène liquide (LOX) alimentent les principaux moteurs de fusées comme SpaceXs Raptor, Blue Origin , BE-4 et NASAS Space Launch System. La gestion de ces fluides dans des environnements de microgravité et de chaleur extrême présente des défis techniques uniques.

Des outils comme COMSOL Multiphysique avec le module CFD offrent un cadre solide pour simuler le comportement des fluides cryogéniques. En accouplement le flux de fluide, le transfert de chaleur, la thermodynamique et le changement de phase dans un seul environnement, les ingénieurs peuvent tester les conceptions des réservoirs, les lignes d'alimentation et les stratégies d'isolation avant la construction du matériel.

Pourquoi les fluides cryogéniques exigent une simulation spécialisée

Contrairement aux liquides ordinaires, les cryogènes présentent de fortes variations de propriété avec la température et la pression. L'hydrogène liquide, par exemple, a une densité d'environ 70 kg/m3 à son point d'ébullition et une conductivité thermique qui change fortement près du point critique. Les régimes d'ébullition – nucléés, transitions et ébullitions de film – peuvent coexister dans un seul réservoir en raison de l'apport calorifique localisé du moteur ou du rayonnement solaire.

Ces complexités rendent les corrélations empiriques peu fiables pour la conception. La simulation multiphysique doit résoudre simultanément:

  • Transfert de chaleur conjugué: conduction à travers les parois des réservoirs, convection dans le cryogène, et radiations provenant de sources externes.
  • Dynamique du changement de phase: évaporation, condensation et gaz réel près de la ligne de saturation.
  • Régimes de débit en deux phases : débits bubble, stratifié, annulaire et de brouillard dans les tuyaux et les réservoirs.
  • Interaction de la structure du fluide: déformation du réservoir sous des contraintes thermiques cryogéniques et accélérations de lancement.

COMSOL="s approche multiphysique permet aux ingénieurs de combiner ces phénomènes dans un modèle unique plutôt que d'aller entre des outils séparés.

Capacités de base de COMSOL CFD pour la modélisation cryogénique

Équations et modèles de turbulence

Le module CFD résout les équations Navier-Stokes pour la masse, l'impulsion et l'énergie, couplées à un modèle de turbulence approprié. Pour les réservoirs cryogéniques, le modèle k‐↓ SST fonctionne bien parce qu'il résout avec précision les débits de parois proches, critique pour prédire le transfert de chaleur à travers la paroi du réservoir.

Changement de phase et flux multiphase

COMSOL propose deux approches principales pour le changement de phase:[
Niveau Set[ et Les méthodesPhase Field[ suivent explicitement l'interface liquide-vapor. Elles sont idéales pour les problèmes de surface libre tels que le sloshing de réservoir ou l'impingement de jet.Mixture Model[ est mieux adapté aux flux biphasés dispersés (boucle ou gouttelette).Pour l'ébullition, l'interface Thermal Phase Change[ peut utiliser des fonctions définies par l'utilisateur pour les taux de transfert de masse basés sur des gradients de température.

Thermodynamique non idéale

COMSOL="s La fonction Gaz réel vous permet d'importer des équations fondamentales d'état telles que Peng‐Robinson ou Span‐Wagner pour des prédictions de propriétés très précises.Pour l'hydrogène et l'oxygène, ces équations sont essentielles près du point critique, où les hypothèses de gaz idéal donnent des erreurs dépassant 30 %.

Bibliothèques de matériel personnalisé

Vous pouvez définir des tables de propriétés à partir de NIST REFPROP data directement dans COMSOL, garantissant la densité, la chaleur spécifique, la viscosité et la conductivité thermique sont des fonctions de température et de pression.

Processus de modélisation étape par étape dans le COMSOL

1. Définition de la géométrie

Pour un réservoir à propergol sphérique typique, inclure l'épaisseur de paroi, les ports d'entrée/sortie, les chicanes antivortex et les couches isolantes. Simplifiez les filets non critiques et les trous de boulons pour garder le maillage gérable.

2. Attribution des propriétés du matériau

Créer un matériau personnalisé pour le cryogène à l'aide de tables NIST. Par exemple, pour le LH2 à 20 K à 30 K, importer la densité en fonction de la température et de la pression, et régler la viscosité dynamique à l'aide du modèle Lennard‐Jones. Pour la paroi du réservoir (par exemple, aluminium 2219 ou acier inoxydable 316L), définir la conductivité thermique, la densité et le coefficient de dilatation thermique.

3. Interfaces physiques et couplage

Ajouter les interfaces suivantes à partir de l'assistant modèle:

  • Flux de laminage ou Flux de turbulence (k‐= SST)[ pour le mouvement du fluide.
  • Transfert de chaleur dans les fluides et Transfert de chaleur dans les solides pour le transfert de chaleur conjugué.
  • Champ de phase (ou ensemble de niveaux) pour le suivi de l'interface, si l'on étudie la dynamique en deux phases.
  • Raccordement multiphysique avec flux non isotherme pour tenir compte de la densité et de la flottabilité dépendantes de la température.

Pour les modèles à ébullition, ajouter un Source de chaleur à l'interface égale à la chaleur latente de la vaporisation multipliée par le taux d'évaporation. Ceci peut être implémenté par une expression définie par l'utilisateur qui active lorsque la température du liquide atteint le point de saturation à la pression locale.

4. Conditions de la frontière

Conditions typiques pour un réservoir de propergol pendant une phase côtière:

  • Surface extérieure de la paroi de la citerne: flux de chaleur provenant du rayonnement solaire (par exemple, 1,3 kW/m2 à 1 UA plus infrarouge de la Terre).
  • Surface intérieure de la paroi du réservoir : sans glissement, flux de chaleur conjugué.
  • Interface liquide-vapor : continuité de la température, bilan massique avec vitesse de variation de phase.
  • Débit (ligne d'alimentation): limite de pression égale à la pression de l'éclaboussure du réservoir moins la tête.
  • Entrée (pour le remplissage): condition limite du débit massique.

Pour les simulations de microgravité, utilisez une Force volumétrique pour la tension de surface (basée sur le modèle de force de surface continuelle) et fixez l'accélération à zéro ou à une petite valeur micro-g (10−4 g).

5. Stratégie de mesure

Utiliser un maillage de couche limite (au moins 5 à 10 couches de prisme) à la paroi du réservoir pour capter les gradients de température abrupts. Dans la région d'interface en deux phases, raffiner les éléments à une taille d'environ 1 mm (pour un réservoir de 1 m de diamètre) pour résoudre la courbure de l'interface. Utiliser un maillage Triangulaire libre pour les modèles axisymétriques 2D ou ]Tétraédral libre avec extrusion de couche limite pour la 3D.

Effectuez une étude de convergence des mailles : exécutez le modèle avec un double raffinement et comparez les profils de température et les positions d'interface. Sélectionnez le maillage plus grossier qui donne des résultats dans les 2 % de la solution raffinée.

6. Paramètres du solvant et pas de temps

Les problèmes cryogéniques sont souvent transitoires.Utilisez un formule de différenciation en arrière (BDF) solveur de l'ordre 1 ou 2 avec une étape initiale de 0,01 s. Pour la physique couplée, activez le solveur complètement couplé avec un facteur d'amortissement de 0,8 pour améliorer la convergence. Si le modèle fonctionne lentement, passez à une approche ségréguée, résolvant d'abord le flux, puis le transfert de chaleur, puis le champ de phase dans chaque étape temporelle.

Surveiller la convergence en vérifiant les résidus (moins de 10 à 4) et en conservant les quantités (balance massique dans un rayon de 0,1 %).

7. Traitement et analyse postérieurs à la mise en oeuvre

Principaux résultats à extraire :

  • Des contours température[ sur les parois des réservoirs et le domaine fluide.
  • Taux de décompression (kg/s) intégré sur l'interface.
  • Évolution de la pression[ dans le volume de l'échelle (important pour la conception de la structure du réservoir).
  • Streamlines ou lignes de chemin montrant des zones de recirculation qui peuvent causer une stratification thermique.
  • Position d'interface[ en fonction du temps – critique pour l'évaluation des dispositifs d'acquisition de liquides.

Les outils COMSOL=Valeurs dérivées peuvent calculer automatiquement la fuite totale de chaleur, la masse de vapeur générée et la fraction vide.

Étude de cas: Simulation de la bouillure dans un réservoir de propergol cryogénique

Considérez un réservoir sphérique LH2 (rayon 1,5 m) avec des parois en aluminium d'une épaisseur de 1 mm et 5 cm de MLI. Le réservoir est plein à 50 %, à une pression de 3 bar, avec le fluide à température de saturation (environ 23 K).

Après avoir mis en place le modèle décrit ci-dessus, la simulation montre que le MLI réduit la fuite de chaleur à environ 8 W/m2. Le taux d'ébullition s'est stabilisé à 0,25 g/s après un passage initial de 200 s. Les balayages paramétriques révèlent que le doublement de l'épaisseur du MLI réduit l'ébullition de 45 %, mais la pénalité pour poids peut ne pas être acceptable. Le modèle prévoit également un gradient de température de 2 K entre le fond et le dessus du réservoir, ce qui aide à identifier les endroits où la condensation de vapeur pourrait se produire, ce qui pourrait entraîner une dégradation de la pression du réservoir et des difficultés de réinflammation.

Ces simulations permettent de concevoir directement la mission. Par exemple, elles peuvent aider à dimensionner un système de ventilation thermodynamique (TVS) qui refroidit activement le réservoir pour éviter la surpressurisation.

Applications dans le cycle de vie des engins spatiaux

Conception du système d'alimentation en propergol

Les modèles COMSOL de la ligne d'alimentation – y compris les virages, les jonctions et les vannes – peuvent prédire des baisses de pression et détecter des régions où l'ébullition peut survenir en raison de points chauds locaux. Le Modèle de variance avec vitesse de glissement capture la séparation de phase dans les flux de microgravité, aidant les ingénieurs à concevoir des dispositifs d'acquisition de liquides (LAD) tels que des canaux d'écran ou des vanes.

Validation du système de protection thermique

Les matériaux d'isolation, y compris les aérogels et les couvertures MLI, peuvent être modélisés en milieux poreux ou en couches solides avec conductivité thermique efficace. En effectuant une simulation à l'état stable du réservoir avec chauffage dans le pire des cas, les ingénieurs vérifient que l'isolation maintient le cryogène sous son point d'ébullition tout au long de la mission.

Lancer les opérations de refroidissement-décharge et de remplissage du paillasson

Le remplissage d'un réservoir avec du cryogène dans des conditions ambiantes produit un violent ébullition et un choc thermique. Les simulations transitoires dans COMSOL peuvent évaluer le temps nécessaire pour refroidir les parois du réservoir à des températures sûres (souvent inférieures à 100 K) avant le remplissage principal.

Gestion des propergols en vol

Les simulations de structure de fluide en combinaison dans COMSOL permettent aux ingénieurs de concevoir des bacs, des diaphragmes ou des réservoirs de tension de surface qui maintiennent le liquide en position stabilisée. La méthode ]Phase Field capture l'interface mobile pendant les brûlures du moteur ou les propulseurs d'assiette, fournissant des forces sur les parois du réservoir pour l'analyse de structure.

Correlations des essais au sol

Avant de pouvoir faire confiance à une simulation pour le vol, elle doit être validée en fonction des données d'essai au sol. Le module d'estimation des paramètres du COMSOL peut ajuster des paramètres incertains (par exemple, résistance au contact, facteurs de vision du rayonnement) pour correspondre aux valeurs de thermocouple obtenues lors d'un essai à l'échelle inférieure.

Ressources externes pour les sujets avancés

Pour les ingénieurs qui cherchent à approfondir leurs connaissances en matière de DFC cryogénique, les ressources suivantes sont recommandées :

Orientations futures : Jumelles numériques et apprentissage automatique

La prochaine frontière de la modélisation cryogénique des fluides est la création de jumeaux numériques en temps réel. Des modèles à commande réduite (ROM) peuvent être construits à partir de balayages paramétriques COMSOL, puis intégrés dans le logiciel de vol du véhicule pour l'estimation de l'état. Par exemple, un ROM qui prédit le taux d'ébullition en fonction de la fuite de chaleur, le niveau de remplissage et l'accélération peuvent fonctionner en millisecondes, permettant le jaugeage à bord du propergol.

Les réseaux neuronaux formés sur les bases de données de simulation COMSOL peuvent prédire l'effondrement de la pression après des manœuvres de réorientation plus rapides qu'un modèle multiphysique complet. Ces approches hybrides maintiennent la précision tout en réduisant le coût de calcul par ordre de grandeur.

COMSOL , LiveLink pour MATLAB permet d'interfacer facilement les données de simulation avec ces outils d'IA, permettant aux ingénieurs d'améliorer leurs ROM de façon itérative avec des données fraîches de haute fidélité au fur et à mesure que la mission progresse.

Conclusion

La modélisation de la dynamique cryogénique des fluides dans le COMSOL CFD offre aux ingénieurs de l'espace un environnement puissant et unifié pour relever les défis multiphysiques des flux à basse température. De la conception initiale du concept aux opérations de vol, ces simulations permettent de mieux choisir l'isolation, de mieux utiliser les systèmes d'alimentation et de gérer plus efficacement les propergols.

Les exemples et les flux de travail présentés ici offrent un point de départ pour construire des modèles CCD cryogéniques fiables. L'exploration spatiale poussant vers un espace plus profond et des missions plus longues, la gestion exacte des fluides cryogéniques devient encore plus critique.