An LES model with finite-rate phase change and subgrid spray based on a thermodynamically consistent four-equation multiphase model

Ce travail propose un modèle de simulation LES intégrant un changement de phase à taux fini et un spray de sous-maille, basé sur un modèle multiphasique à quatre équations thermodynamiquement cohérent, qui permet de prédire avec précision les surfaces interfaciales et les régimes d'évaporation tout en optimisant le coût de calcul.

L'essentiel

Le défi : Simuler le "chaos" d'un moteur

Imaginez que vous essayez de prédire exactement comment une goutte de carburant va se transformer en nuage de gaz à l'intérieur d'un moteur de voiture ou d'une fusée. C'est un cauchemar de physicien ! C'est un mélange de liquide qui explose en gouttelettes, de gaz qui s'échappe, de changements de température brutaux et de pressions énormes.

Si on essaie de tout calculer avec une précision absolue, l'ordinateur mettrait des années à finir la simulation. Si on simplifie trop, le résultat est faux et le moteur pourrait exploser en vrai.

L'idée des chercheurs : "Le compromis intelligent"

Les chercheurs de Stanford et de l'université de KTH ont créé un nouveau modèle mathématique (appelé LES) qui agit comme un "super-filtre".

1. La métaphore de la foule (Le modèle à 4 équations)

Imaginez que vous regardez une foule immense dans un stade.

• L'approche classique (trop lourde) : Vous essayez de suivre chaque personne, chaque mouvement de bras, chaque battement de cœur. C'est impossible à gérer.
• L'approche des chercheurs (le modèle à 4 équations) : Vous regardez la foule comme un tout. Vous supposez que dans un petit groupe de personnes, tout le monde a à peu près la même température et la même pression. Cela permet de gagner un temps fou tout en restant très précis sur le mouvement global de la "masse".

2. La métaphore du nuage de parfum (Le spray subgrid)

Dans un moteur, il y a des gouttelettes tellement minuscules qu'elles sont "invisibles" pour l'ordinateur (elles sont plus petites que les pixels de notre simulation).
Pour ne pas les perdre, les chercheurs ont inventé une sorte de "nuage de probabilité". Au lieu de dessiner chaque micro-goutte, ils calculent une "densité de surface". C'est comme si, au lieu de dessiner chaque grain de poussière dans un rayon de soleil, vous calculiez simplement la luminosité du nuage de poussière. Cela permet de savoir combien de surface est disponible pour que le carburant s'évapore.

3. La métaphore du thermostat (Le changement de phase)

Le carburant doit passer de l'état liquide à l'état gazeux (l'évaporation). Mais cela ne se fait pas instantanément.
Les chercheurs ont ajouté une règle de "limite thermodynamique". Imaginez un thermostat : si vous demandez à une pièce de passer de 0°C à 100°C en une seconde, le thermostat dira "Stop, c'est physiquement impossible !". Le modèle des chercheurs empêche l'ordinateur de faire des erreurs absurdes (comme faire apparaître du gaz là où il n'y a pas assez d'énergie) en s'assurant que le changement respecte toujours les lois de la nature.

Le résultat : Ça marche !

Pour vérifier si leur "calculateur intelligent" fonctionnait, ils l'ont testé sur un cas réel très célèbre dans le monde de l'ingénierie : le Spray A de l'ECN (un test de pulvérisation de carburant ultra-précis).

Le verdict : Les courbes dessinées par l'ordinateur collent presque parfaitement aux mesures réelles faites en laboratoire.

En résumé

Ce papier, c'est l'art de "tricher intelligemment". Les chercheurs ont trouvé comment simplifier les calculs les plus complexes de la physique (mélange de liquides, évaporation, explosions de gouttelettes) sans perdre la réalité. C'est une étape cruciale pour concevoir des moteurs de demain : plus propres, plus puissants et surtout, plus sûrs.

Résumé technique

Résumé Technique : Modèle LES avec changement de phase à taux fini et spray subgrid basé sur un modèle multiphasique à quatre équations thermodynamiquement cohérent

Problématique

La simulation numérique des systèmes multiphasiques multicomposants à grande échelle (moteurs à combustion, propulsion de fusée) est extrêmement complexe. Les défis majeurs résident dans la nécessité de capturer simultanément le mélange intra-phase, les interfaces à fort rapport de densité et les transferts de masse (évaporation, ébullition, cavitation).

Les modèles existants font souvent un compromis : soit ils utilisent le modèle d'équilibre homogène (HEM), qui suppose un changement de phase instantané (ce qui est physiquement incorrect pour les structures d'interface non résolues), soit ils utilisent des modèles de non-équilibre complet, qui sont extrêmement coûteux en ressources de calcul. Le problème central est de trouver un équilibre entre la précision physique des taux de changement de phase et l'efficacité computationnelle nécessaire pour les simulations de grande échelle (LES - Large Eddy Simulation).

Méthodologie

Les auteurs proposent un cadre de modélisation basé sur un modèle multiphasique à quatre équations, qui suppose un équilibre de pression, de température et de vitesse au niveau de la sous-maille (subgrid). Pour compenser les restrictions de ce modèle, ils introduisent deux innovations majeures :

1. Changement de phase à taux fini thermodynamiquement borné : Au lieu d'un équilibre instantané, ils utilisent un modèle basé sur la théorie cinétique de Hertz-Knudsen. Pour garantir la cohérence thermodynamique et éviter les transferts de masse non physiques (surpassement de l'équilibre), ils ont formulé une nouvelle expression où le taux de transfert est limité par l'évolution du potentiel chimique (énergie libre de Gibbs).
2. Modèle de spray subgrid ($\Sigma$) "phase-confined" : Pour prédire l'aire interfaciale de la spray non résolue par la maille, ils utilisent une équation de transport pour la densité de surface $\Sigma'$. L'innovation réside dans une formulation qui "confine" la diffusion de cette surface uniquement dans la phase gazeuse, évitant ainsi les fuites non physiques de la phase liquide à travers les interfaces.
3. Indicateurs de régime : Le modèle utilise des indicateurs basés sur la fraction volumique liquide et le nombre de Jakob pour distinguer dynamiquement les régimes de spray (dense vs dilué) et les modes de changement de phase (flashing/ébullition vs évaporation).

Contributions Clés

• Développement d'un modèle de changement de phase à taux fini qui est intrinsèquement lié à l'équilibre du potentiel chimique, garantissant qu'il ne produit jamais d'états thermodynamiques impossibles.
• Introduction d'une nouvelle forme de l'équation $\Sigma$ pour la LES, qui limite la diffusion de la surface interfaciale à la phase gazeuse, assurant une robustesse accrue lors des simulations de spray.
• Couplage efficace entre un modèle de transport de surface subgrid et un modèle de changement de phase cinétique, permettant de simuler des phénomènes complexes comme le flashing (vaporisation rapide due à une chute de pression) et l'évaporation de gouttelettes.
• Réduction du coût computationnel en utilisant le cadre à quatre équations tout en maintenant une précision de haute fidélité.

Résultats

Le modèle a été validé par deux cas d'étude basés sur les expériences du réseau de combustion de moteurs (ECN Spray A) :

• Cas non-évaporant : La simulation reproduit avec précision la densité de masse projetée (PMD), l'aire de surface projetée (PSA) et le diamètre moyen de Sauter (SMD), validant ainsi la dynamique du spray et le modèle de surface $\Sigma$.
• Cas évaporant : Le modèle capture avec succès la longueur de pénétration du liquide et de la phase gazeuse. Les résultats montrent un excellent accord avec les mesures expérimentales de pénétration, confirmant la capacité du modèle à gérer le couplage entre la dynamique du spray et le transfert de masse thermique.

Signification et Impact

Ce travail représente une avancée significative pour la simulation numérique des fluides (CFD) appliquée à l'ingénierie de combustion. En fournissant un cadre qui est à la fois thermodynamiquement cohérent et computationnellement efficace, il permet aux chercheurs de simuler des sprays complexes dans des conditions réelles (haute pression, haute température) avec une fidélité proche de la réalité, sans nécessiter la puissance de calcul prohibitive des modèles de non-équilibre complets. Cela ouvre la voie à une meilleure conception des injecteurs de carburant et des systèmes de propulsion.