Compressible and immiscible fluids with arbitrary density ratio

Cet article présente une nouvelle théorie fondée sur la thermodynamique pour la modélisation de fluides compressibles et immiscibles présentant des rapports de densité arbitrairement élevés, en remédiant aux limites des équations traditionnelles de Navier-Stokes et d'Euler par la dérivation d'équations d'évolution de la densité qui rendent compte correctement de la véritable évolution de la quantité de mouvement.

L'essentiel

Imaginez que vous essayiez de décrire une danse entre deux partenaires très différents : un éléphant lourd et aux mouvements lents (l'eau) et une plume légère et rapide (l'air). Dans le monde de la physique des fluides, ces deux-là sont « immiscibles », ce qui signifie qu'ils ne se mélangent pas comme l'huile et l'eau, mais qu'ils possèdent une frontière floue là où ils se rencontrent.

Pendant longtemps, les scientifiques ont lutté pour écrire les « règles de la danse » (les équations) pour des systèmes où la différence de densité est immense, comme un cas où l'éléphant serait 1 000 fois plus lourd que la plume. Les anciennes règles présentaient une faille majeure : elles traitaient le poids de l'éléphant comme s'il s'agissait d'un nombre constant et immuable, même précisément à la frontière floue où l'éléphant se transforme en plume. C'est comme essayer de décrire une personne en train de se transformer en fantôme en disant : « Elle reste un être humain solide tout au long du processus », ce qui n'a aucun sens.

Voici une décomposition simple de ce que ce papier propose de corriger.

1. Le problème avec les anciennes règles

La méthode traditionnelle pour calculer le mouvement des fluides (en utilisant les équations de Navier-Stokes et d'Euler) repose sur un raccourci appelé l'approximation de Boussinesq.

• L'analogie : Imaginez que vous poussez un chariot de courses. Si le chariot est rempli de briques, il est lourd. S'il est vide, il est léger. Les anciennes règles supposaient que si vous poussez un chariot qui est partiellement rempli de briques et partiellement d'air, le poids du chariot ne change jamais pendant que vous le poussez. Cela suppose simplement que le poids est une moyenne fixe.
• La faille : En réalité, à mesure que le chariot se déplace et que les briques se déplacent (diffusion), le poids change. Les anciennes règles ignoraient le fait que la « quantité de mouvement » (masse $\times$ vitesse) change parce que la masse elle-même change à la frontière. Elles supposaient également que l'air et l'eau occupent un espace parfaitement prévisible, ignorant que lorsqu'ils se mélangent à la frontière, l'espace qu'ils occupent peut en fait osciller et changer en raison de la pression.

2. La nouvelle approche : La minimisation de l'énergie

Au lieu de deviner comment la densité change, l'auteur part d'un principe fondamental : la nature cherche toujours à utiliser le moins d'énergie possible.

• L'analogie : Pensez à une balle roulant le long d'une colline. Elle ne se soucie pas des étapes spécifiques qu'elle emprunte ; elle veut simplement atteindre le bas (l'énergie la plus basse). L'auteur utilise ce concept de « colline d'énergie » pour dériver de nouvelles règles sur la façon dont l'eau et l'air interagissent.
• L'innovation clé : L'auteur introduit un concept appelé « Volume Excédentaire ».
  • Imaginez que vous ayez un seau d'eau et un seau d'air. Si vous les versez ensemble, vous pourriez vous attendre à ce que le volume total soit exactement la somme des deux seaux. Mais au niveau microscopique, lorsqu'ils se rencontrent, les molécules peuvent se serrer plus étroitement ou plus lâchement, créant un peu d'espace « supplémentaire » ou « manquant ».
  • Les anciennes règles supposaient que cet espace supplémentaire était nul partout. Ce papier affirme : « Non, cet espace supplémentaire existe, il change d'un endroit à l'autre, et il affecte la densité. »

3. Les nouvelles « Règles de la Danse » (Les Résultats)

En tenant compte de ce changement de « l'espace supplémentaire » et en utilisant la minimisation de l'énergie, l'auteur dérive un nouvel ensemble d'équations qui accomplissent trois choses principales :

A. Une nouvelle façon d'entendre le son (Vitesse du son)
Le papier montre que la vitesse du son n'est pas seulement un nombre aléatoire ; elle découle directement de la façon dont l'énergie du fluide change lorsqu'il est comprimé.

• La métaphore : Pensez au son comme à une ondulation dans une foule. La vitesse de cette ondulation dépend de la manière dont les gens (les molécules) sont serrés et de l'énergie qu'ils possèdent. La nouvelle formule calcule cette vitesse naturellement, sans avoir besoin qu'on lui indique la valeur au préalable. Elle suggère même que, dans un gaz, la vitesse du son est approximativement la même que la vitesse moyenne à laquelle les molécules de gaz rebondissent.

B. Une nouvelle règle pour la pression et la vitesse (La loi de Bernoulli)
Vous avez probablement déjà entendu parler du principe de Bernoulli : « Lorsqu'un fluide se déplace plus vite, sa pression chute. »

• Le rebondissement : L'ancienne règle fonctionne très bien pour l'eau circulant dans un tuyau, mais elle échoue lorsqu'il y a un saut de densité énorme (comme de l'eau frappant de l'air). L'auteur crée une Loi de Bernoulli Généralisée.
• La métaphore : Imaginez une rivière se jetant dans une cascade. L'ancienne règle dit que l'énergie reste la même. La nouvelle règle dit : « Attendez, alors que l'eau se transforme en brume (air), une partie de l'énergie est perdue ou transformée parce que la "densité/étouffement" de l'eau est en train de changer. » La nouvelle équation tient compte de ce transfert d'énergie, ce qui la rend précise même lorsque le fluide change de nature, passant de lourd à léger.

C. Le « Bossage » de la densité
C'est peut-être le résultat le plus visuel.

• La vieille vision : Si l'on regarde la frontière entre l'eau et l'air, les anciens modèles disaient que la densité descendrait simplement de manière fluide, comme une rampe descendante.
• La nouvelle vision : Les mathématiques de l'auteur prédisent un bossage. En traversant la frontière, la densité augmente légèrement avant de chuter vers le niveau de l'air.
• La métaphore : Imaginez une foule de personnes (l'eau) essayant de passer par une porte vers un couloir vide (l'air). En passant la porte, elles pourraient se serrer un peu plus étroitement pendant une fraction de seconde avant de s'éparpiller. La nouvelle théorie prédit ce « bossage de compactage », ce qui correspond à ce que les simulations informatiques avancées (appelées Théorie de la Fonctionnelle de la Densité) ont déjà observé, mais que les anciens modèles simples avaient manqué.

Résumé

Ce papier propose une nouvelle façon d'écrire les lois de la physique pour les fluides qui ont des poids très différents (comme l'eau et l'air).

1. Il cesse de prétendre que le poids est constant à la frontière.
2. Il admet que l'espace occupé par les molécules change (volume excédentaire).
3. Il utilise le principe de « l'énergie la plus basse » pour dériver de nouvelles règles qui expliquent comment le son voyage, comment la pression change, et pourquoi la densité à la frontière eau-air présente en fait un petit « creux » (ou bossage).

L'auteur affirme que ce nouveau cadre fonctionne pour n'importe quel mélange de fluides, peu importe la différence de poids, ouvrant la voie à des simulations informatiques plus précises de phénomènes tels que la chute de la pluie, le déferlement des vagues ou la montée des bulles.

Résumé technique

Résumé technique : Fluides compressibles et immiscibles avec des rapports de densité arbitraires : Théorie de l'interface diffuse

Énoncé du problème
Les modèles actuels de mécanique des fluides numérique (CFD) peinent à simuler avec précision les systèmes de fluides immiscibles présentant des rapports de densité élevés, tels que l'eau et l'air (rapport $\approx 1000$). Les approches existantes reposent sur deux méthodes principales, lesquelles possèdent des limitations significatives :

1. L'approximation de Boussinesq : Elle suppose une densité ($\rho$) constante au sein de l'interface, utilisant l'équation $\rho(du/dt) = -\nabla p + \dots$. Cette formulation ne parvient pas à rendre compte de l'évolution « réelle » de la quantité de mouvement $d(\rho u)/dt$, négligeant les changements de densité et le transfert de quantité de mouvement qui en résulte. De plus, elle ne peut pas générer naturellement les forces de gravité et de flottabilité sans ajouts ad hoc, car celles-ci découlent des différences de densité.
2. La méthode d'interpolation : Elle définit la densité comme une fonction linéaire de la concentration volumique ($\rho = \rho_w \phi + \rho_a(1-\phi)$). Cette approche crée un conflit fondamental entre l'incompressibilité ($d\rho/dt = 0$) et la diffusion ($d\phi/dt \neq 0$). Elle suppose également un volume excédentaire uniforme (souvent nul) à travers le domaine, ce qui contredit la réalité physique où le volume excédentaire local dépend de la pression locale. Par conséquent, les équations de continuité standards et l'interpolation linéaire échouent à capturer les profils de densité non monotones observés à l'interface eau-air.

Méthodologie
Les auteurs proposent un nouveau cadre théorique fondé sur le principe thermodynamique de minimisation de l'énergie, distinct de la dérivation traditionnelle des équations de Navier-Stokes et d'Euler. La méthodologie comprend :

• Reformulation du volume et de la densité : Au lieu de supposer un volume excédentaire nul, le modèle introduit un volume excédentaire ($v_e$) spatialement variable au sein de l'interface. Ce volume est redistribué entre les composants du fluide, conduisant à des densités partielles ($\rho'_w, \rho'_a$) qui ne sont plus constantes mais dépendent de la pression locale et du mélange. Cela nécessite deux équations d'évolution couplées : une pour la concentration volumique ($\phi$) et une pour la densité ($\rho$).
• Principe de dissipation-conservation : Les auteurs établissent un théorème reliant le premier principe de la thermodynamique (conservation) au second principe (dissipation). En définissant l'énergie totale du système ($E$) comme la somme de l'énergie libre ($F$), de l'énergie de pression ($P$), de l'énergie cinétique ($K$) et de l'énergie libre des parois ($W$), ils dérivent les équations cinétiques en minimisant la dissipation d'énergie.
• Relation de Gibbs-Duhem modifiée : Pour résoudre les incompatibilités mathématiques entre les gradients spatiaux et les dérivées fonctionnelles dans les modèles existants, une relation de Gibbs-Duhem isotherme modifiée est dérivée. Elle rend explicitement compte des changements d'entropie dépendants de la vitesse et du travail inertiel, menant à une définition rigoureuse de la pression qui inclut des contributions de l'énergie libre et de l'entropie de vitesse.
• Définition de l'énergie cinétique : L'énergie cinétique est formulée comme $K = \int \rho u \cdot u \, d\Omega$ (densité de quantité de mouvement $\zeta = \rho u$) plutôt que la forme conventionnelle $\frac{1}{2}\rho u^2$. Cette distinction permet de séparer le bilan de quantité de mouvement (changeant à la fois $u$ et $\rho$) de la dissipation de vitesse (changeant uniquement $u$).

Contributions clés
La dérivation produit une équation d'évolution de la densité généralisée et trois avancées théoriques majeures :

1. Loi de Bernoulli généralisée : Le principe de Bernoulli standard ($p + \frac{1}{2}\rho u^2 = \text{constante}$) est montré comme étant un cas particulier. La forme généralisée tient compte des changements d'entropie de vitesse et de l'évolution de la densité à l'interface : $p + \rho u^2 = \text{constante}$ (sous des conditions spécifiques d'état stationnaire non visqueux), reconnaissant que des sauts d'énergie se produisent à travers les interfaces en raison des différences de densité.
2. Formulation de la vitesse du son universelle : Une expression générale de la vitesse du son est dérivée basée sur la dérivée de l'énergie libre par rapport à la densité : $c = \sqrt{d\Psi/d\rho}$. Pour les gaz parfaits, cela s'approche de la vitesse thermique moyenne des molécules.
3. Équation d'état (EOS) généralisée : La théorie démontre que l'équation d'évolution de la densité génère naturellement une EOS pour les fluides (par exemple, $p \sim \rho \ln \rho$ pour les gaz parfaits) sans nécessiter une EOS prescrite, contrastant avec les méthodes classiques d'interface diffuse.

Résultats et validation
Le papier valide la théorie à travers trois études de cas :

• Propagation du son : L'équation de la vitesse du son dérivée est cohérente avec l'analyse par perturbation et prédit correctement la relation de dispersion des ondes sonores dans l'air.
• EOS des gaz parfaits : À l'état stationnaire, le modèle reproduit la loi des gaz parfaits, confirmant que l'équation d'évolution de la densité capture intrinsèquement les relations thermodynamiques d'état.
• Dynamique de l'interface eau-air : Une solution analytique pour une interface eau-air plane en mouvement dans un tube démontre que le modèle prédit un profil de densité non monotone. Plus précisément, la densité augmente de la valeur de l'eau de masse volumique jusqu'à un pic maximal au centre de l'interface avant de diminuer vers la valeur de l'air, ce qui est conforme aux résultats de la théorie de la fonctionnelle de la densité (DFT) [11, 12] et contredit la diminution monotone prédite par les méthodes d'interpolation linéaire. Le modèle reproduit également correctement l'équation de Young-Laplace pour les interfaces courbes.

Signification
Le document affirme fournir un cadre général pour les fluides immiscibles avec des rapports de densité arbitrairement élevés, applicable tant aux régimes compressibles qu'incompressibles. En traitant l'évolution de la densité comme une variable dynamique couplée au volume excédentaire et à la dissipation visqueuse, le modèle résout les contradictions inhérentes à l'approximation de Boussinesq et aux méthodes d'interpolation linéaire. Le travail suggère que la dissipation visqueuse altère intrinsèquement la densité, remettant en question l'hypothèse stricte d'incompressibilité dans les écoulements dissipatifs. Les auteurs postulent que cette théorie ouvre une nouvelle fenêtre pour la CFD, offrant une description physiquement cohérente du transfert de quantité de mouvement et d'énergie aux interfaces fluide-fluide là où les rapports de densité sont extrêmes.