⚡ electrical engineering

CFD Model Predictions of Zero Boil‑Off Unvented and Vented Tank Filling in Microgravity

Cet article présente des modèles de CFD 2D et 3D validés par des expériences de remplissage de réservoirs en microgravité afin de caractériser la relation entre les taux de remplissage, les régimes d'écoulement de type geyser et l'évolution de la pression, démontrant finalement que les geysers stables produisent des pressions de pointe plus élevées et que la modélisation 3D est essentielle pour capturer les comportements de geysers instables.

Auteurs originaux : Rebecca L. Winter, Mohammad Kassemi

Publié 2026-07-03

📖 6 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Rebecca L. Winter, Mohammad Kassemi

Article original sous licence CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Imaginez que vous essayez de remplir un ballon d'eau géant et vide flottant dans l'espace. Sur Terre, la gravité tire l'eau vers le bas, elle s'accumule donc au fond. Mais dans l'espace, il n'y a pas de « bas ». Si vous projetez un jet de liquide froid dans ce ballon flottant, il ne se contente pas de se déposer ; il se comporte comme une fontaine sauvage et imprévisible.

Ce document est une étude de simulation informatique sur ce problème précis : comment remplir un réservoir de carburant super-froid dans l'espace sans qu'il n'explose ou ne s'évapore par ébullition.

Voici une décomposition de ce que les chercheurs ont découvert, en utilisant des analogies simples :

Le gros problème : L'effet « Geyser »

Lorsque vous projetez un liquide froid dans un réservoir chaud et vide dans l'espace, le liquide ne reste pas simplement là. Il jaillit comme un geyser.

Le problème : À mesure que le liquide froid frappe l'air chaud à l'intérieur du réservoir, l'air essaie de se condenser (redevenir liquide), mais le liquide essaie aussi de bouillir. Cela crée un bras de fer qui modifie la pression à l'intérieur du réservoir.
Le danger : Si la pression devient trop élevée trop vite, le réservoir pourrait éclater. Si vous essayez de laisser sortir le gaz (ventiler) pour sauver le réservoir, vous perdez du carburant précieux. L'objectif est de remplir le réservoir sans laisser de gaz s'échapper (Zéro Évaporation/Zero Boil-Off) et sans que le réservoir n'éclate.

Les quatre « pas de danse » du liquide

Les chercheurs ont utilisé un super-ordinateur pour observer comment le liquide se comporte. Ils ont découvert que le liquide exécute quatre « pas de danse » distincts selon la vitesse à laquelle on le projette :

Le Geyser Stable Faible : Le liquide jaillit un peu et reste en place, comme une fontaine calme. Il est stable et prévisible.
Le Geyser Oscillant : Le liquide jaillit mais commence à pencher et à osciller d'avant en arrière comme une méduse ou un danseur ivre. Il est toujours stable, mais il se déplace autour du centre.
Le Détachement de Gouttelettes : Le liquide jaillit, mais à mesure qu'il monte en hauteur, il se fragmente en gouttelettes individuelles, comme un tuyau d'arrosage qui projette trop fort. Ces gouttes volent à travers le réservoir et frappent l'autre paroi.
La Colonne de Jet Pleine : Le liquide traverse tout le réservoir en un seul flux solide et ininterrompu, comme un rayon laser d'eau frappant le mur opposé.

La découverte surprenante : Lent ne signifie pas toujours sûr

Vous pourriez penser : « Si je remplis le réservoir lentement, la pression restera basse. » Les modèles informatiques ont montré que ce n'est pas toujours vrai.

L'analogie : Imaginez que vous remplissez un seau avec un tuyau d'arrosage. Si vous ouvrez l'eau en un tout petit filet (vitesse lente), l'eau peut jaillir en une ligne droite et stable qui frappe le fond et crée un énorme éclat de chaleur. Cela peut en fait provoquer un pic de pression énorme dès le début.
La réalité : Si vous ouvrez l'eau un peu plus vite, le jet peut se briser en gouttelettes. Ces gouttelettes se mélangent mieux avec l'air, refroidissant l'air plus efficacement. Cela peut en fait entraîner des pics de pression plus faibles que le simple filet lent.
La leçon : Ralentir le débit de remplissage ne garantit pas la sécurité ; cela peut parfois aggraver le pic de pression parce que le comportement du liquide change.

Le « Nombre Magique » (Nombre de Weber)

Les chercheurs ont trouvé un moyen de prédire quel « pas de danse » le liquide fera en utilisant un concept mathématique appelé le Nombre de Weber.

Considérez cela comme un panneau de limitation de vitesse pour le liquide.
Si le nombre est inférieur à 1,5, le liquide a tendance à rester calme (Geyser Stable).
Si le nombre est supérieur à 1,7, le liquide a tendance à devenir sauvage (Gouttelettes ou Jet Plein).
La zone grise : Entre 1,5 et 1,7, c'est incertain. Le liquide peut être calme ou sauvage, et l'ordinateur ne pouvait pas prédire parfaitement lequel simplement en regardant la vitesse. Cela suggère que d'autres facteurs cachés sont en jeu.

2D vs 3D : Le plat contre le réel

Les chercheurs ont modélisé cela de deux manières :

2D (Plat) : Comme regarder l'ombre du réservoir. C'est plus rapide à calculer mais cela manque les mouvements de danse « oscillants ». Le modèle pense que le liquide reste parfaitement droit, ce qui n'est pas vrai dans l'espace.
3D (Réel) : Comme regarder le réservoir réel. C'est nécessaire pour voir le geyser oscillant et les gouttelettes qui volent partout. Le modèle 2D donne parfois de mauvaises réponses sur la pression car il ne peut pas voir le liquide osciller.

Qu'en est-il des évents ?

Ils ont également examiné ce qui se passe si l'on ouvre une valve pour laisser sortir le gaz (Ventilé).

Le résultat : Si vous laissez sortir le gaz, la pression reste basse et la température reste stable. Le liquide se comporte presque de la même manière que dans le cas du réservoir scellé, mais sans les pics de pression dangereux.
Le piège : Si vous ne ventilez pas, la pression augmente, la température change et le liquide se comporte différemment. Vous ne pouvez pas utiliser une simulation « ventilée » pour prédire ce qui se passe dans un réservateur « scellé ».

Le cas du vide

Ils ont également simulé le remplissage d'un réservoir qui était presque un vide (vide d'air).

Le début : Lorsque le liquide frappe initialement le vide, il se transforme instantanément en gaz (bout de façon violente), créant un immense nuage.
La stabilisation : Après environ 30 secondes, la pression augmente suffisamment pour que le liquide cesse de bouillir et commence à se comporter comme les cas normaux de « geyser » décrits précédemment.

L'essentiel à retenir

Pour remplir un réservoir de carburant dans l'espace en toute sécurité, sans perdre de carburant ni faire exploser le réservoir, vous ne pouvez pas simplement deviner la vitesse. Vous devez comprendre exactement comment le liquide va « danser » (geyser, oscillation ou fragmentation en gouttelettes). Parfois, aller plus lentement provoque un pic de pression plus élevé, et vous avez besoin d'un modèle informatique 3D pour voir l'image complète. Le « nombre magique » (nombre de Weber) aide à prédire le comportement, mais il existe encore une zone complexe où le liquide est difficile à prédire.

Résumé technique : Prédictions de modèles CFD du remplissage de réservoirs avec et sans évents en régime de zéro évaporation (Zero Boil-Off) en microgravité

Énoncé du problème
Le stockage et le transfert sûrs et efficaces de carburants cryogéniques liquides en microgravité constituent un défi critique pour les futures missions spatiales de longue durée de la NASA. Le ravitaillement en orbite est essentiel pour les missions vers la Lune et Mars, pourtant l'absence de gravité complique la gestion des fluides. Une préoccupation majeure est l'augmentation de la pression à l'intérieur d'un réservoir pendant le remplissage, ce qui peut résulter des fuites thermiques ou du processus de remplissage lui-même, pouvant entraîner la rupture de conduites ou l'auto-inflammation du carburant. Bien que l'évent du réservoir empêche l'augmentation de la pression, cela réduit l'efficacité du transfert. Par conséquent, le remplissage « sans évent » (non ventilé) est un mode opérationnel privilégié mais complexe qui nécessite une compréhension approfondie des écoulements diphasiques, du transfert de chaleur et de la dynamique de changement de phase en microgravité. Cette étude répond au besoin de modèles prédictifs pour comprendre l'évolution de la pression transitoire et les comportements de l'interface liquide-gaz (spécifiquement la formation de « geysers ») lors du remplissage de réservoirs avec ou sans évents en microgravité.

Méthodologie
Les auteurs ont employé une approche de modélisation multidimensionnelle combinant l'analyse thermodynamique et la mécanique des fluides numérique (CFD) :

Modèle thermodynamique : Un modèle thermodynamique homogène de dimension zéro (0D) a été développé pour prédire la pression finale à l'état stationnaire. Ce modèle suppose l'homogénéité du fluide à la température et à la pression de saturation, en résolvant les équations de conservation de la masse et de l'énergie pour établir une limite théorique du système.
Modèles CFD : Des simulations CFD biphasiques ont été réalisées avec ANSYS Fluent en utilisant le schéma Volume de Fluide (VOF), une formulation explicite et un suivi d'interface abrupte.
- Dimensions : Des modèles axisymétriques 2D et paramétriques 3D ont été utilisés.
- Fluide : Le perfluoro-n-pentane (PNP) a été utilisé comme fluide de simulation, avec des propriétés dérivées de la base de données DIPPR à 253,18 K.
- Conditions : Les simulations ont couvert des taux de remplissage de 1,0 à 1,65 mL/s dans un réservoir cylindrique (hauteur de 0,094 m, rayon de 0,03 m) sous des conditions de microgravité ( $9,81 \times 10^{-6}$ m/s²).
- Scénarios : L'étude a examiné les configurations non ventilées (sans évent) et ventilées (évents latéraux et supérieurs), divers profils de vitesse d'entrée (écoulement de type plug plat et écoulement de Poiseuille parabolique), ainsi que différentes conditions initiales du réservage (vide, partiellement pré-rempli et quasi-vide).
- Validation : Des études d'indépendance de maillage et de pas de temps ont été effectuées pour garantir la précision des solutions. Les résultats de la CFD ont été validés par rapport à des données expérimentales antérieures issues d'études en tour de chute (ex. Symons et al., Govindan et Dreyer).

Contributions clés et résultats

Classification du comportement des geysers : Les modèles CFD 3D ont identifié quatre comportements de geyser distincts basés sur la vitesse d'entrée et le profil d'écoulement :
1. Geyser stable de faible amplitude (Low Stable Geyser) : Un petit monticule de liquide avec une croissance constante.
2. Geyser oscillant (Wobbly Geyser) : Une colonne de liquide stable mais s'inclinant ou se courbant de manière erratique.
3. Détachement de gouttelettes (Droplet Shedding) : Un régime instable où la colonne de liquide se fragmente en gouttelettes.
4. Colonne de jet complet (Full Jet Column) : Un régime instable où le liquide atteint la paroi opposée du réservoir sans se briser.
  L'étude a révélé que les modèles axisymétriques 2D ne pouvaient pas capturer le comportement du « geyser oscillant » en raison de sa nature asymétrique, prédisant souvent un étranglement non physique à la place.
Évolution de la pression et pics transitoires : Une découverte significative est la relation entre la stabilité du geyser et la pression transitoire. Contrairement à l'hypothèse selon laquelle des taux de remplissage plus faibles produisent toujours des pressions plus basses, les modèles ont montré que les cas de geysers stables (faible et oscillant) produisaient des pics de pression transitoire plus élevés que les cas instables. Ceci est attribué à la surface d'interface plus petite dans les geysers stables, ce qui réduit l'efficacité de la condensation de la vapeur. Les geysers instables (détachement de gouttelettes/jets) fournissent des surfaces plus grandes, facilitant une condensation rapide et atténuant les pics de pression.
Seuil du nombre de Weber : L'étude a confirmé que le nombre de Weber (ajusté pour le profil de vitesse d'entrée) prédit efficacement la transition entre les régimes stables et instables. Un seuil critique compris entre 1,5 et 1,7 a été identifié, ce qui est cohérent avec les découvertes expérimentales antérieures. Cependant, les auteurs ont noté des incohérences au sein de cette plage, suggérant que le nombre de Weber seul ne suffit peut-être pas à distinguer les geysers stables de faible amplitude des geysers oscillants, impliquant que d'autres paramètres influencent la stabilité.
Accord Thermodynamique vs CFD : Bien que le modèle thermodynamique 0D ait réussi à prédire la pression finale à l'état stationnaire (convergeant vers la limite thermodynamique), il a échoué à capturer les dépassements de pression transitoires observés dans les modèles CFD. Le modèle thermodynamique a également suggéré de manière incorrecte que des taux de remplissage plus faibles résulteraient toujours en des pressions initiales plus basses, manquant la nuance des pics de pression induits par le geyser.
Remplissage ventilé vs non ventilé et conditions initiales :
- Remplissage ventilé : Les formes de geyser dans les cas ventilés étaient presque identiques aux cas non ventilés, bien que l'évent empêche l'augmentation de la pression et de la température. Les évents latéraux entraînaient des conditions isothermes, tandis que les évents supérieurs permettaient une certaine augmentation de la température en raison de l'obstruction partielle par le liquide.
- Conditions initiales : Les réservoirs partiellement pré-remplis présentaient des augmentations de température réduites par rapport aux réservoirs vides. Le remplissage en quasi-vide présentait une phase initiale d'évaporation par flash, mais après environ 30 secondes, le comportement du geyser convergeait pour ressembler aux cas de remplissage standards.
- Effets de la gravité : Les comparaisons entre la microgravité et la gravité zéro (0g) n'ont montré aucune différence discernable dans le comportement des geysers, validant l'utilisation des données de tour de chute (0g) pour prédire les performances des missions en orbite terrestre basse et dans l'espace profond.

Signification
L'article affirme que la compréhension des comportements spécifiques des geysers est cruciale pour prédire avec précision l'évolution de la pression transitoire lors du remplissage de réservoirs non ventilés, un facteur souvent négligé par les modèles thermodynamiques simplifiés. La recherche souligne que la réduction des taux de remplissage n'est pas une stratégie universellement efficace pour prévenir la surpression, car elle peut involontairement faire passer le système dans un régime de geyser stable avec des pics de pression plus élevés. En validant les modèles CFD par rapport aux données expérimentales et en établissant les limites des approximations axisymétriques 2D pour les écoulements asymétriques, l'étude fournit un cadre plus robuste pour la conception de systèmes de gestion de fluides cryogéniques pour les futures missions spatiales. Ce travail souligne la nécessité de capturer les détails de la physique des écoulements diphasiques pour garantir la sécurité et l'efficacité des opérations de ravitaillement en orbite.