Local Thermal Non-Equilibrium Models in Porous Media: A Comparative Study of Conduction Effects

Cette étude compare les modèles de non-équilibre thermique local (LTNE) à une référence résolue au niveau des pores pour des milieux poreux purement conductifs, démontrant que les modèles à l'échelle REV utilisant des paramètres effectifs basés sur l'homogénéisation capturent avec précision les effets de résistance interfaciale, tandis que les modèles à double réseau avec une résolution spatiale fixe présentent une plus grande déviation.

L'essentiel

Imaginez un matériau poreux, comme une éponge ou une roche, comme une ville animée. Cette ville compte deux types de résidents : les grains solides (les bâtiments) et le fluide (l'eau ou l'air circulant dans les rues).

Lorsque vous chauffez un côté de cette ville, vous voulez savoir comment la température se propage. La grande question que pose cet article est : les bâtiments et l'eau à l'intérieur chauffent-ils exactement à la même vitesse, ou prennent-ils du retard l'un par rapport à l'autre ?

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

1. Les deux façons de penser la chaleur

Par le passé, les scientifiques supposaient généralement l'Équilibre Thermique Local (ETL).

• L'analogie : Imaginez une pièce remplie de personnes se tenant par la main. Si une personne a chaud, tout le monde le ressent instantanément. Dans ce modèle, les « bâtiments » et l'« eau » sont si parfaitement connectés qu'ils ont toujours exactement la même température à n'importe quel endroit. C'est comme s'ils partageaient un seul cerveau.

Cependant, les chercheurs savaient que ce n'est pas toujours vrai. Parfois, la connexion entre le bâtiment et l'eau est « collante » ou lente. C'est le Non-Équilibre Thermique Local (NETL).

• L'analogie : Imaginez que les personnes sont dans des pièces séparées avec des portes épaisses et isolées entre elles. Si vous chauffez l'eau dans le couloir, les bâtiments peuvent rester froids pendant un moment, car la chaleur doit lutter pour traverser la porte. L'eau devient chaude, mais le bâtiment reste froid un certain temps. Ils ont des températures différentes au même endroit.

2. Les trois « cartes » utilisées pour prédire la chaleur

Pour déterminer quand ce « retard » se produit et comment le prédire, l'équipe a comparé trois façons différentes de dessiner une carte de cette ville :

• Carte A : La vue « au niveau de la rue » (Modèle résolu aux pores)

  • Ce que c'est : C'est la carte la plus détaillée. Elle dessine chaque bâtiment et chaque rue. Elle voit la forme exacte de la roche et de l'eau.
  • Le hic : Elle est incroyablement lente et coûteuse en calcul, comme essayer de simuler chaque grain de sable d'une plage. Les chercheurs l'ont utilisée comme leur « référence or » pour vérifier si les autres cartes étaient correctes.

• Carte B : La vue « quartier » (Modèle à double réseau)

  • Ce que c'est : Au lieu de dessiner chaque rue, cette carte simplifie la ville en un réseau de points (représentant les bâtiments et les poches d'eau) reliés par des lignes (représentant les connexions entre eux).
  • Le hic : Elle est plus rapide, mais elle a une résolution fixe. C'est comme regarder une ville à travers une grille de fenêtres ; vous ne pouvez pas zoomer plus près que la taille de la fenêtre. L'article a constaté que, comme cette grille est fixe, elle manque parfois les changements brusques de température se produisant juste aux bords.

• Carte C : La vue « aérienne » (Modèle à l'échelle du VRE)

  • Ce que c'est : C'est une carte de haut niveau, moyennée. Elle ne voit pas les bâtiments individuels ; elle voit des « blocs » de la ville. Elle utilise les mathématiques pour deviner le comportement moyen de tout le bloc.
  • Le hic : Pour que cela fonctionne, vous devez deviner les « propriétés moyennes » du bloc. Si vous vous trompez, toute la carte est fausse.

3. La grande expérience

Les chercheurs ont effectué des simulations sur ordinateur pour voir comment la chaleur se déplaçait dans cette « ville » dans deux conditions différentes :

Scénario 1 : La porte ouverte (faible résistance)

• Le montage : La connexion entre l'eau et la roche était parfaite (comme une porte grande ouverte). La chaleur circulait librement.
• Le résultat : La « porte ouverte » signifiait que l'eau et la roche chauffaient instantanément ensemble. L'hypothèse ETL (le cerveau unique) fonctionnait parfaitement. Les trois cartes donnaient presque la même réponse. Le « retard » n'existait pas.

Scénario 2 : La porte isolée (haute résistance)

• Le montage : La connexion était bloquée ou « collante » (comme une porte épaisse et isolée). La chaleur avait du mal à passer de l'eau à la roche.
• Le résultat : Maintenant, l'eau devenait chaude, mais la roche restait froide un certain temps. L'hypothèse ETL échouait complètement.
  • La Carte au niveau de la rue montrait le retard exact.
  • La Carte aérienne (si calculée correctement en utilisant une méthode mathématique spécifique appelée homogénéisation) correspondait très bien à la Carte au niveau de la rue.
  • La Carte quartier était correcte, mais parce que ses « fenêtres » avaient une taille fixe, elle lissait un peu trop les différences nettes.

4. La conclusion clé

La découverte la plus importante concerne la façon dont vous calculez la carte « aérienne ».

• Certaines anciennes méthodes de calcul des propriétés moyennes pour la carte aérienne ignoraient la « porte collante ». Elles supposaient que le transfert de chaleur était toujours parfait. Lorsque les chercheurs ont utilisé ces anciennes formules, la carte aérienne n'a pas réussi à montrer le retard entre l'eau et la roche.
• Cependant, lorsqu'ils ont utilisé une méthode mathématique spécifique et plus avancée (l'homogénéisation) qui prenait en compte la « porte collante » (la résistance interfaciale), la carte aérienne est devenue incroyablement précise. Elle correspondait presque parfaitement à la vue détaillée au niveau de la rue, même si elle était beaucoup plus simple.

Résumé

• Si la connexion est parfaite : Vous pouvez utiliser des modèles simples ; tout chauffe ensemble.
• Si la connexion est lente/collante : Vous devez utiliser des modèles qui permettent à l'eau et à la roche d'avoir des températures différentes.
• Le meilleur raccourci : Si vous devez modéliser un système immense (comme un aquifère entier ou une pile à combustible) et que vous ne pouvez pas simuler chaque grain, utilisez le modèle « aérien », mais assurez-vous d'utiliser les mathématiques spécifiques qui tiennent compte de la résistance entre les matériaux. Si vous faites cela, votre modèle simple sera aussi précis que celui qui est super détaillé.

Note : L'article indique explicitement que cette étude n'a examiné que la chaleur se déplaçant à travers des matériaux immobiles (conduction). Ils n'ont pas examiné la chaleur se déplaçant avec de l'eau en écoulement (convection), ce qu'ils disent étudier dans un article futur.

Résumé technique

1. Énoncé du Problème

La modélisation du transfert de chaleur dans les milieux poreux repose souvent sur l'hypothèse de l'Équilibre Thermique Local (ETL), qui postule que les phases fluide et solide partagent la même température instantanément au sein d'un volume de contrôle. Bien que valable pour de nombreux scénarios, cette hypothèse échoue dans les systèmes caractérisés par :

• De forts gradients de température.
• Des différences significatives de propriétés thermiques entre les phases.
• De fortes résistances thermiques interfaciales (faibles coefficients de transfert de chaleur).
• Une dynamique rapide.

Lorsque l'ETL échoue, des modèles de Non-Équilibre Thermique Local (NETL) sont nécessaires pour résoudre les températures distinctes des phases fluide ($T_f$) et solide ($T_s$) ainsi que l'échange de chaleur entre elles. Cependant, il manque d'études comparatives complètes validant les modèles NETL à l'échelle du continuum (Dual-Réseau et Échelle du Volume Élémentaire Représentatif - VER) par rapport à des références haute fidélité résolues au niveau des pores, notamment concernant la précision des formulations de paramètres effectifs (par exemple, homogénéisation par rapport aux corrélations empiriques).

2. Méthodologie

L'étude compare trois approches de modélisation à l'échelle du continuum sur un système entièrement saturé (un fluide, un solide) en considérant la conduction thermique pure (négligeant la convection et le transport de masse).

A. Les Trois Classes de Modèles

1. Modèle Résolu aux Pores (Référence) :

   • Approche : Résout les équations de bilan d'énergie séparément pour les domaines fluide ($\Omega_f$) et solide ($\Omega_s$) sur une grille résolvant la géométrie exacte à l'échelle des pores.
   • Interface : Résout explicitement l'interface nette fluide-solide ($\Gamma_{fs}$).
   • Implémentation NETL : Utilise un coefficient de transfert de chaleur fini ($h$) pour permettre un saut de température à travers l'interface ($\lambda_f \nabla T_f \cdot n = h(T_s - T_f)$).
   • Résolveur : Méthode des Éléments Finis (MEF) utilisant Netgen/NGSolve.

2. Modèle Dual-Réseau :

   • Approche : Représente le milieu poreux comme deux réseaux interconnectés (fluide et solide) composés de corps idéalisés (sphères) et de gorges (cylindres).
   • Interface : Ne résout pas l'interface nette ; utilise plutôt des gorges interfaciales avec une transmissivité thermique basée sur des facteurs de forme et des surfaces effectives.
   • Résolveur : Méthode des Volumes Finis aux Éléments (CVFE) utilisant DuMux.

3. Modèle à l'Échelle du VER (Volume Élémentaire Représentatif) :

   • Approche : Utilise des équations de bilan moyennées sur un volume de contrôle où les températures fluide et solide se chevauchent.
   • Implémentation NETL : Résout deux équations d'énergie couplées avec un terme d'échange de chaleur interfacial ($q_{cond} = \lambda_{eff,I} a_{fs} (T_s - T_f)$).
   • Paramètres Effectifs : L'étude évalue trois formulations différentes pour les conductivités thermiques effectives et les termes interfaciaux :
     • Nuske et al. : Mise à l'échelle simple par fraction volumique.
     • Nakayama et al. : Inclut la tortuosité et les rapports de conductivité.
     • Théorie de l'Homogénéisation : Dérive mathématiquement les paramètres effectifs à partir des équations à l'échelle des pores, incorporant explicitement le coefficient de transfert de chaleur interfacial ($h$).
   • Résolveur : Méthode des Volumes Finis (MVF) utilisant PorePy.

B. Configuration de la Simulation

• Géométrie : Un domaine 3D périodique composé de cellules de référence cubiques contenant une sphère connectée à six tubes.
• Matériaux : Granit (solide) et Eau (fluide).
• Conditions aux Limites : Configuration quasi-1D avec un échelon de température sur la limite gauche (fluide chauffé de 10 K, solide isolé).
• Cas de Test :
  • Cas 1 (Faible Résistance) : Coefficient de transfert de chaleur élevé ($h \approx 8,3 \times 10^7$ W/m²K), mimant les interfaces eau-granit standard.
  • Cas 2 (Haute Résistance) : Coefficient de transfert de chaleur faible ($h = 100$ W/m²K), mimant des interfaces à haute résistance (par exemple, cuivre-hélium) pour induire de forts effets NETL.

3. Contributions Clés

• Comparaison Systématique : Fournit une référence rigoureuse des modèles Dual-Réseau et à l'échelle du VER par rapport à une référence résolue aux pores pour des conditions NETL.
• Analyse des Formulations de Paramètres : Évalue l'efficacité des différentes formulations de paramètres effectifs. Elle démontre que les approches basées sur l'homogénéisation sont supérieures pour capturer le comportement NETL car elles intègrent mathématiquement le coefficient de transfert de chaleur interfacial, tandis que les formulations empiriques (Nuske, Nakayama) échouent souvent à prédire le NETL dans des systèmes purement conductifs sans dépendance explicite à $h$.
• Impact de la Résolution : Met en évidence le rôle critique de la résolution spatiale. Le modèle Dual-Réseau a montré des écarts par rapport à la référence en raison de sa résolution spatiale fixe, qui ne peut pas être affinée pour capturer les gradients de température raides près des limites, contrairement au modèle à l'échelle du VER où l'affinement de la grille est possible.

4. Résultats Clés

• Faible Résistance Interfaciale (Haut $h$) :

  • Le système se comporte comme un ETL ; les températures fluide et solide sont presque identiques dans tous les modèles.
  • Le modèle à l'échelle du VER (utilisant des paramètres homogénéisés) correspond étroitement à la référence résolue aux pores.
  • Le modèle Dual-Réseau a montré des températures légèrement plus basses et des gradients moins raides par rapport à la référence, attribués à sa résolution spatiale fixe et à la diffusion numérique.

• Haute Résistance Interfaciale (Bas $h$) :

  • Effets NETL : Des différences de température significatives entre les phases fluide et solide ont été observées, en particulier près de la limite chauffée.
  • Performance des Modèles :
    • Le modèle à l'échelle du VER avec paramètres homogénéisés a capturé avec précision le comportement NETL et le saut de température, correspondant étroitement à la référence résolue aux pores (bien que légèrement plus rapide en raison de la diffusion numérique).
    • Le modèle Dual-Réseau a de nouveau montré un gradient moins raide et des températures plus basses en raison de son incapacité à affiner le maillage à l'interface.
    • ETL vs NETL : Les modèles ignorant la résistance interfaciale (ETL) n'ont pas réussi à capturer le saut de température, confirmant la nécessité de modèles NETL pour les scénarios à haute résistance.

• Formulations de Paramètres Effectifs :

  • Dans le cas à haute résistance, les formulations de Nuske et Nakayama (qui n'incluent pas intrinsèquement $h$ dans la limite conductive) n'ont pas réussi à produire des effets NETL, prédisant à la place un ETL.
  • Seule l'approche basée sur l'homogénéisation a réussi à capturer le comportement NETL car le terme interfacial mis à l'échelle ($H = h \cdot a_{fs}$) dépend explicitement du coefficient de transfert de chaleur.

5. Signification et Conclusion

• Validation de l'Homogénéisation : L'étude confirme que la théorie de l'homogénéisation fournit les paramètres effectifs les plus fiables pour les modèles NETL à l'échelle du VER, car elle fait le lien rigoureux entre la physique à l'échelle des pores (incluant la résistance interfaciale) et l'échelle du continuum.
• Limites des Modèles Dual-Réseau : Bien qu'utiles pour les insights à l'échelle des pores, les modèles Dual-Réseau sont limités par leur résolution spatiale fixe. Ils peinent à résoudre les gradients thermiques nets près des limites par rapport aux modèles de continuum qui permettent l'affinement de la grille.
• Implications Pratiques : Pour les applications impliquant une forte résistance interfaciale (par exemple, certains systèmes géothermiques ou cryogéniques), s'appuyer sur des paramètres effectifs empiriques sans considérer le coefficient de transfert de chaleur peut conduire à des erreurs significatives (prédire un ETL alors que le NETL existe).
• Travaux Futurs : Les auteurs notent que la convection est un facteur dominant dans la plupart des applications pratiques et prévoient d'étendre cette étude comparative pour inclure le transport de masse et le transfert de chaleur convectif dans un article ultérieur.

En résumé, le papier établit que les modèles à l'échelle du VER utilisant des paramètres effectifs homogénéisés constituent l'approche de continuum la plus robuste pour simuler le NETL dans les milieux poreux, à condition que le coefficient de transfert de chaleur interfacial soit explicitement pris en compte dans la formulation du modèle.