Monolithic kinetic algorithm for heterogeneous porous media systems using a continuous one-domain approach

Résumé Technique : Algorithme Cinétique Monolithique pour les Systèmes de Milieux Poreux Hétérogènes

Énoncé du Problème
La simulation de l'écoulement des fluides dans les milieux poreux hétérogènes, en particulier dans les systèmes couplant des régions de fluide libre avec des structures poreuses, présente des défis importants. Les approches traditionnelles reposent souvent sur des formulations à « deux domaines », qui résolvent des équations de conservation distinctes pour les régions de fluide libre et poreuses et nécessitent des conditions aux limites d'interface explicites (par exemple, les conditions de Beavers-Joseph) pour les coupler. Alternativement, les approches à « un seul domaine » utilisent un ensemble unique d'équations mais ont historiquement éprouvé des difficultés de cohérence dans les milieux hétérogènes.

Les formulations existantes de la méthode de Boltzmann sur réseau (LBM) pour les milieux poreux souffrent souvent d'une mise à l'échelle incohérente de la porosité dans les équations macroscopiques récupérées. Spécifiquement, de nombreux modèles précédents échouent fréquemment à mettre à l'échelle correctement le gradient de pression et les termes de contrainte visqueuse en fonction de la porosité lorsque celle-ci varie spatialement. De plus, de nombreuses formulations expriment la contrainte visqueuse de Brinkman en termes de gradients de la vitesse intrinsèque plutôt qu'en termes de la vitesse superficielle (Darcy), ce qui entraîne des inexactitudes dans les domaines hétérogènes. Par ailleurs, les modèles LBM standards manquent souvent de la flexibilité nécessaire pour ajuster indépendamment la viscosité volumique effective ou la vitesse du son, ce qui limite l'efficacité computationnelle et la stabilité.

Méthodologie
Les auteurs proposent un modèle de Boltzmann sur réseau monolithique à domaine unique basé sur une approche continue à un seul domaine. La méthodologie repose sur les piliers suivants :

1. Modèle Macroscopique : Les hydrodynamiques cibles sont dérivées par la méthode de la moyenne volumique, résultant en un modèle de Darcy-Forchheimer-Brinkman généralisé. Ce modèle est valide sur l'ensemble du domaine de calcul, traitant la transition entre les régions de fluide libre et poreuses comme une variation continue des propriétés physiques (porosité $\epsilon$). Les équations de conservation comprennent :

   • Une équation de continuité pour la densité moyenne volumique ($\epsilon\rho$).
   • Une équation de quantité de mouvement incorporant la traînée de Darcy, la traînée inertielle de Forchheimer et un tenseur de contrainte visqueuse effectif ($\sigma_{eff}$).
   • Crucialement, la contrainte visqueuse effective est exprimée en termes de gradients de la vitesse superficielle ($v$), et non de la vitesse intrinsèque ($u = v/\epsilon$), garantissant la cohérence avec la formulation de Brinkman pour une porosité variable.

2. Solveur Cinétique : Le modèle utilise une équation cinétique généralisée unique sur des réseaux standards (par exemple, D2Q9, D3Q19).

   • Distribution d'Équilibre : Une distribution d'équilibre discrète sous forme de produit est employée pour assurer les moments corrects et l'invariance de Galilée. Cela inclut des termes de correction spécifiques ($\Phi_{ext}$ et $\Phi_{\eta}$) pour restaurer l'isotropie dans le tenseur de contrainte visqueuse et permettre un ajustement indépendant de la viscosité volumique.
   • Schéma de Force : Un schéma de force généralisé cohérent (méthode Relax-to-Force) est utilisé pour incorporer la traînée de Darcy-Forchheimer, les forces de corps et les corrections nécessaires pour les gradients de porosité.
   • Mise à l'échelle de la Porosité : Les champs localement conservés sont définis comme la densité et la quantité de mouvement moyennes volumiques, incorporant naturellement la mise à l'échelle de la porosité dans les moments cinétiques.
   • Évaluation de la Vitesse : Un algorithme itératif est proposé pour résoudre la relation implicite entre la vitesse macroscopique et les termes de force, gérant la traînée non linéaire de Darcy-Forchheimer et les corrections visqueuses.

3. Conditions aux Limites : Les auteurs adaptent les conditions aux limites de nœuds mouillés standards (schéma d'équilibre, extrapolation de non-équilibre et schéma de Tamm-Mott-Smith) à la formulation de milieu poreux. Une procédure spécifique est dérivée pour calculer la densité du fluide intrinsèque aux frontières en utilisant la continuité de la masse, tenant compte de la présence du terme de force et des variations de porosité.

Contributions Clés

• Formulation à Domaine Unique Cohérente : L'article présente la première formulation LBM qui récupère un tenseur de contrainte visqueuse de Brinkman avec une mise à l'échelle de porosité pleinement cohérente, exprimée en termes de vitesse superficielle. Cela permet au modèle de gérer une porosité variant spatialement sans conditions d'interface artificielles.
• Correction des Termes de Pression et Visqueux : Les auteurs dérivent et implémentent des corrections de force spécifiques pour rectifier les incohérences dans le gradient de pression et le tenseur de contrainte visqueuse qui surviennent dans les milieux hétérogènes, garantissant que le modèle récupère la limite macroscopique correcte.
• Ajustement Indépendant des Paramètres : Le modèle permet de spécifier la vitesse du son et la viscosité volumique effective indépendamment de la viscosité de cisaillement. Cela améliore l'efficacité computationnelle (en permettant des vitesses du son plus élevées) et la stabilité numérique (en ajustant la viscosité volumique).
• Algorithme Monolithique : L'approche utilise une seule équation cinétique et un seul solveur pour l'ensemble du domaine, simplifiant la mise en œuvre pour des géométries complexes impliquant à la fois du fluide libre et des milieux poreux.

Résultats et Validation
Le modèle a été validé par une série de tests de référence et d'études numériques :

• Analyse de Cohérence : Une analyse spectrale temporelle linéaire a confirmé que le modèle récupère la vitesse du son effective correcte et les taux de dissipation pour les modes de cisaillement et acoustiques. L'analyse a démontré que l'équilibre en forme de produit avec les corrections proposées préserve l'invariance de Galilée et l'isotropie, contrairement aux formulations d'équilibre de second ordre précédentes.
• Écoulements de Référence :
  • Écoulement de Couette : Les simulations dans un canal poreux homogène sous l'effet de la gravité montrent un excellent accord avec les solutions numériques de référence pour les profils de vitesse et de densité à travers divers nombres de Reynolds et de Darcy.
  • Écoulement de Poiseuille : Le modèle capture avec précision l'écoulement dans des canaux homogènes et, de manière cruciale, dans des canaux présentant des profils de porosité hétérogènes. Il prédit correctement les effets de canalisation de l'écoulement où le fluide circule préférentiellement à travers les régions de haute porosité.
• Démonstrations Numériques :
  • Écoulement autour d'un Cylindre Perméable : Le modèle a simulé avec succès l'écoulement confiné autour d'un cylindre poreux, identifiant cinq régimes de flux distincts basés sur le nombre de Reynolds et la porosité du cœur. Il a capturé des phénomènes tels que la stabilisation du sillage, le détachement de la recirculation et les instabilités de sillage lointain, montrant un bon accord qualitatif avec la littérature.
  • Cavité Entraînée par un Couvercle avec Obstacle Poreux : Le modèle a simulé une cavité avec un obstacle poreux interne, démontrant comment la perméabilité affecte la formation de vortex, la pénétration du flux et les motifs de recirculation. Les résultats ont mis en évidence la transition entre un comportement solide et un comportement perméable à mesure que la porosité augmente.

Signification et Revendications
Les auteurs affirment que ce travail fournit un cadre robuste et unifié pour simuler le transport dans les milieux poreux hétérogènes. En éliminant le besoin de conditions aux limites d'interface explicites et en assurant une mise à l'échelle cohérente de tous les termes hydrodynamiques, le modèle est particulièrement adapté pour :

• Les études de conception préliminaire à l'échelle d'un dispositif pour des structures poreuses d'ingénierie avec une porosité spatialement variable.
• La simulation d'interactions complexes entre fluides libres et corps perméables (par exemple, piles à combustible, systèmes de filtration, combustion poreuse).
• La capture des phénomènes de transport dans des systèmes où la transition entre l'écoulement libre et le milieu poreux est continue ou là où les interfaces nettes sont difficiles à résoudre.

L'article conclut que bien que l'étude actuelle se concentre sur la validation en 2D, la formulation est intrinsèquement multidimensionnelle et prête pour une extension aux géométries 3D. Des travaux futurs sont identifiés comme nécessaires pour inclure le bilan énergétique (transfert de chaleur), les mélanges multicomposants et les écoulements multiphasiques afin d'aborder pleinement le transport couplé dans des systèmes d'ingénierie réalistes.