Consistent multiple-relaxation-time lattice Boltzmann method for the volume averaged Navier-Stokes equations

Cet article présente une méthode de Boltzmann sur réseau à temps de relaxation multiple qui, grâce à une équation d'état provisoire et un terme de pénalité, résout de manière cohérente et précise les équations de Navier-Stokes moyennées en volume pour les écoulements multiphasiques, éliminant ainsi les vitesses parasites et les incohérences des schémas traditionnels basés sur la densité.

L'essentiel

🌊 Le Secret pour Simuler les Fluides dans un Monde "Poreux"

Imaginez que vous essayez de prédire comment l'eau coule à travers une éponge, ou comment le sang circule dans un réseau de petits vaisseaux, ou encore comment le sable bouge dans un lit fluidisé. C'est un casse-tête énorme pour les mathématiciens et les ingénieurs.

Pour résoudre ce problème, ils utilisent une méthode appelée Lattice Boltzmann (LB). On peut imaginer cette méthode comme une armée de millions de petits robots qui sautent d'une case à l'autre sur une grille virtuelle. En suivant leurs mouvements, on peut reconstituer le comportement global de l'eau ou de l'air.

Mais il y a un gros problème quand on ajoute des obstacles (comme les trous d'une éponge) : les robots se trompent, ils créent des courants fantômes et des erreurs qui faussent tout le calcul. C'est comme si, dans votre simulation, l'eau commençait à couler toute seule vers le haut sans aucune raison !

C'est ici que l'équipe de l'Université Tsinghua (en Chine) intervient avec une nouvelle recette de cuisine mathématique.

1. Le Problème : Les "Fantômes" dans la Simulation

Dans les anciennes méthodes, les robots LB avaient du mal à comprendre la différence entre l'eau pure et l'eau mélangée à des solides (comme dans une éponge).

• L'analogie : Imaginez que vous essayez de compter les voitures sur une autoroute, mais que certaines parties de la route sont bouchées par des travaux. Les anciennes méthodes de calcul confondaient la densité des voitures avec la densité des travaux. Résultat ? Des erreurs de calcul qui faisaient accélérer ou freiner les voitures sans raison (ce qu'on appelle des "vitesses fantômes").

2. La Solution : Une Nouvelle Recette (MRTLB-VANSE)

Les auteurs ont créé une nouvelle méthode, qu'on pourrait appeler le "Chef Cuisinier Multi-Relaxation". Voici comment ils ont réglé les problèmes, point par point :

• Découpler les ingrédients (La "Loi d'État Provisionnelle") :
  Avant, la quantité d'eau et la quantité de "trous" (porosité) étaient liées trop étroitement. Si les trous changeaient, l'eau changeait de poids, ce qui était faux.

  • L'analogie : Imaginez que vous faites une sauce. Avant, si vous ajoutiez des pâtes (les solides), la sauce devenait automatiquement plus lourde, même si vous n'aviez pas ajouté d'eau. La nouvelle méthode dit : "Attends, la sauce et les pâtes sont deux choses différentes". Ils ont créé une règle temporaire qui permet de gérer les pâtes et la sauce séparément, évitant ainsi que la sauce ne devienne bizarrement lourde ou légère.

• Le "Correcteur de Stress" (Le terme de pénalité) :
  Même avec la nouvelle règle, les robots faisaient encore de petites erreurs quand ils calculaient la friction (la viscosité).

  • L'analogie : C'est comme si vous conduisiez une voiture avec un volant qui tremble légèrement. Les auteurs ont ajouté un "amortisseur" intelligent dans leur calcul. Ce n'est pas une force physique réelle, mais une correction mathématique qui annule les tremblements des robots, rendant la conduite (l'écoulement du fluide) parfaitement stable et lisse.

• Le Chef d'Orchestre (MRT) :
  Au lieu de laisser chaque robot décider seul de sa vitesse (ce qui crée du chaos), ils utilisent une méthode où tous les robots écoutent un chef d'orchestre (l'opérateur MRT).

  • L'analogie : C'est la différence entre une foule qui court dans tous les sens et une troupe de danseurs synchronisés. Le chef d'orchestre assure que même si la musique (le fluide) change de rythme, les danseurs restent en harmonie, même dans des situations très difficiles (comme des trous très petits ou très grands).

3. Les Résultats : Pourquoi c'est génial ?

Les chercheurs ont testé leur nouvelle méthode avec des scénarios extrêmes :

• Des éponges avec des trous géants et minuscules : L'ancienne méthode s'effondrait (les calculs divergeaient). La nouvelle méthode a tenu bon.
• Des fluides qui bougent vite et lentement : La nouvelle méthode reste précise, là où les anciennes devenaient imprécises.
• Zéro vitesse fantôme : Dans les zones où il ne devrait pas y avoir d'écoulement, l'eau reste parfaitement immobile.

En résumé :
Ce papier présente une nouvelle façon de simuler les fluides qui circulent dans des milieux complexes (comme le sol, les poumons, ou les réacteurs chimiques). En utilisant une astuce mathématique pour séparer les ingrédients et un chef d'orchestre pour corriger les erreurs, ils ont créé un simulateur beaucoup plus fiable, plus stable et plus précis que ceux utilisés jusqu'à présent.

C'est comme passer d'une carte dessinée à la main, pleine de ratures, à un GPS haute précision qui vous guide parfaitement, même dans les ruelles les plus étroites et les plus embouteillées.

Résumé technique

1. Problématique

Les équations de Navier-Stokes moyennées en volume (VANSE) sont fondamentales pour la simulation de systèmes multiphasiques fluide-solide à l'échelle macroscopique (par exemple, lits fluidisés, écoulements dans les milieux poreux). Bien que la méthode de Boltzmann sur réseau (LBM) soit un outil puissant pour modéliser ces écoulements, les schémas basés sur la densité (density-based) existants souffrent de deux limitations majeures lorsqu'ils sont appliqués aux VANSE :

• Incohérence numérique : Les schémas traditionnels (notamment ceux utilisant l'opérateur de collision BGK à temps de relaxation unique, SRT) génèrent des erreurs numériques significatives, notamment des termes d'erreur d'ordre élevé dans le tenseur des contraintes visqueuses et des gradients de pression incorrects.
• Vitesses parasites (Spurious Velocities) : Dans les champs de fraction de vide (porosité) présentant de forts gradients ou des discontinuités, les schémas précédents produisent des accélérations non physiques, générant des vitesses parasites qui dépassent souvent les limites de compressibilité faible du LBM, rendant la simulation instable.

2. Méthodologie

Les auteurs proposent une nouvelle méthode : une méthode LBM à temps de relaxation multiple (MRT) cohérente pour les VANSE (MRTLB-VANSE), intégrée dans un cadre basé sur la densité. La méthodologie repose sur trois innovations clés :

• Découplage de la fraction de vide et de la densité :
  Au lieu d'utiliser l'équation d'état standard $p = c_s^2 \phi \rho$, les auteurs introduisent une équation d'état provisoire $p = c_s^2 \kappa \rho$ (où $\kappa$ est une constante proche de la fraction de vide minimale). Cela permet de modifier la distribution d'équilibre de la densité (en utilisant un développement de Hermite d'ordre 3) pour découpler la fraction de vide $\phi$ de la densité $\rho$ dans le terme de gradient de pression.

  • Conséquence : La force de correction de pression ne nécessite plus la discrétisation du gradient de fraction de vide ($\nabla \phi$), éliminant ainsi les oscillations non physiques aux interfaces discontinues.

• Opérateur de collision MRT et terme source de pénalité :
  L'utilisation de l'opérateur MRT (Multiple Relaxation Time) permet de traiter séparément les différents modes de relaxation. Pour corriger les erreurs résiduelles dans le tenseur des contraintes visqueuses (qui violeraient l'invariance galiléenne), un terme source de pénalité ($C$) est ajouté dans l'espace des moments. Ce terme compense spécifiquement les erreurs numériques introduites par l'équation d'état provisoire et la force de correction.

• Analyse de Chapman-Enskog :
  Une analyse rigoureuse de Chapman-Enskog a été menée pour démontrer mathématiquement que le schéma proposé récupère exactement les équations VANSE complètes avec une précision d'ordre 2, garantissant la cohérence entre l'échelle mésoscopique (LBM) et l'échelle macroscopique.

3. Contributions Clés

• Développement d'un schéma MRTLB cohérent : Création d'un algorithme LBM basé sur la densité capable de résoudre les VANSE complètes sans les incohérences des schémas SRT précédents.
• Élimination des vitesses parasites : La nouvelle distribution d'équilibre et la force de correction modifiée réduisent drastiquement les vitesses parasites dans les champs de porosité discontinus, un problème critique pour les simulations réalistes.
• Robustesse pour les grands gradients : Le schéma reste stable et précis même lorsque la fraction de vide varie fortement dans l'espace et le temps, là où les méthodes antérieures divergent ou échouent.
• Validation complète : L'article fournit une validation théorique (analyse) et numérique exhaustive, incluant des solutions analytiques, des tests d'écoulement sans flux, et la méthode des solutions manufacturées (MMS).

4. Résultats

Les résultats numériques valident l'efficacité de la méthode MRTLB-VANSE à travers plusieurs benchmarks :

• Écoulements poreux uniformes : La méthode reproduit avec précision les profils de vitesse analytiques pour les écoulements de Poiseuille et de Couette dans des milieux poreux, confirmant la capacité à modéliser l'échelle REV (Volume Élémentaire Représentatif).
• Écoulements non uniformes et discontinus :
  • Dans des tests d'écoulement 1D avec des gradients de porosité, la méthode correspond parfaitement aux solutions analytiques.
  • Dans des tests d'écoulement discontinu (sans flux), la méthode MRTLB-VANSE génère des vitesses parasites nettement inférieures à celles des méthodes de Zhang, Bukreev et Blais. Elle fonctionne même pour des fractions de vide très faibles ($\phi < 0.5$), là où d'autres schémas divergent.
• Convergence d'ordre 2 : La méthode des solutions manufacturées (MMS) confirme que l'erreur de vitesse décroît avec un ordre de convergence de 2 par rapport à l'espacement du réseau, tant pour des écoulements à divergence nulle que non nulle, et pour divers viscosités.
• Comparaison avec les méthodes SRT : Contrairement à la méthode basée sur la pression de Fu et al. (qui utilise un opérateur SRT), la méthode MRTLB-VANSE maintient sa précision et sa stabilité même à haute viscosité (temps de relaxation < 0.5), évitant l'instabilité numérique inhérente aux opérateurs SRT dans ces conditions.

5. Signification et Perspectives

Ce travail représente une avancée significative pour la simulation numérique des écoulements multiphasiques complexes. En résolvant le problème de l'incohérence et des vitesses parasites dans les schémas LBM basés sur la densité, la méthode proposée élargit considérablement le domaine d'application des modèles VANSE.

• Impact industriel : La méthode est prête pour être intégrée dans des modèles industriels complexes impliquant des interactions fluide-solide, la dissolution minérale, la combustion, ou la fracturation.
• Extensibilité : Les auteurs soulignent que l'extension de cette méthode en 3D est directe. De plus, elle ouvre la voie au couplage avec d'autres modèles physiques (champs de concentration, température, modèles multiphasiques immiscibles comme le modèle de gradient de couleur ou pseudo-potentiel) pour simuler des processus liquide-gaz-solide complexes.

En résumé, cette étude fournit un cadre numérique robuste et précis pour la simulation de systèmes multiphasiques à grande échelle, surmontant les limitations historiques des méthodes LBM appliquées aux équations moyennées en volume.