Analysis of heat transfer and water flow with phase change in saturated porous media by bond-based peridynamics

Cet article présente et valide un cadre de peridynamique basé sur les liaisons pour modéliser avec précision le transfert de chaleur couplé et l'écoulement d'eau sous pression avec changement de phase dans les milieux poreux saturés, offrant une approche non locale robuste pour prédire les interfaces de phase et les distributions thermodynamiques dans des scénarios complexes tels que le gel et le dégel des sols.

L'essentiel

Imaginez une éponge complètement imbibée d'eau. Maintenant, imaginez que cette éponge se trouve dans un congélateur, et que l'eau à l'intérieur se transforme lentement en glace. En gelant, l'eau se dilate, l'écoulement de l'eau liquide change, et la température évolue de manière complexe. C'est ce type de problème que les scientifiques rencontrent lorsqu'ils étudient des phénomènes comme le sol gelé (pergélisol) ou la formation de glace dans le sol.

Cet article présente une nouvelle méthode pour simuler et prédire exactement ce qui se passe à l'intérieur de cette « éponge » lorsque l'eau gèle et dégèle, tout en y circulant. Voici un résumé simple de leur travail :

Le Problème : Le « Nœud Inextricable » des Mathématiques

Traditionnellement, les scientifiques utilisent des mathématiques qui examinent un matériau point par point, comme si l'on regardait un seul pixel sur un écran. Cela fonctionne bien pour des objets lisses. Mais lorsque l'eau se transforme en glace, les choses se compliquent :

• Le Problème de la Frontière : La ligne entre l'eau liquide et la glace solide est une cible mobile. C'est comme essayer de tracer une ligne sur une feuille de papier qui continue de bouger et de changer de forme.
• Le « Craquement » : Lorsque l'eau gèle, ses propriétés changent instantanément. Les mathématiques traditionnelles peinent à gérer ces « sauts » soudains ou ces bords nets, provoquant souvent l'arrêt ou des réponses erronées de la simulation informatique.

La Solution : La « Surveillance de Quartier » (Péridynamique)

Les auteurs proposent d'utiliser une méthode appelée Péridynamique à base de liaisons. Au lieu d'examiner un point isolé, imaginez que chaque particule minuscule de l'éponge est une personne dans un quartier.

• L'Horizon : Chaque personne possède un « horizon » (un cercle autour d'elle). Elle ne peut parler et interagir qu'avec ses voisins situés dans ce cercle.
• Les Liaisons : Si deux voisins sont proches, ils sont connectés par une « liaison ».
• La Magie : Dans ce modèle, si une liaison se brise (comme lorsque la glace se forme et bloque l'écoulement de l'eau), les mathématiques ne s'effondrent pas. Le système cesse simplement d'envoyer des messages à travers cette liaison brisée. Cela le rend incroyablement efficace pour gérer les fissures, les fronts de glace mobiles et les changements soudains sans se perdre.

Ce Qu'ils Ont Fait : Les Expériences sur l'« Éponge »

L'équipe a construit un modèle informatique basé sur cette idée de « quartier » pour suivre trois phénomènes se produisant simultanément :

1. Le transfert de chaleur : Comment le froid se propage.
2. Le mouvement de l'eau : Comment le liquide s'écoule à travers l'éponge.
3. Le changement de phase : Comment l'eau se transforme en glace et inversement.

Ils ont testé leur nouveau modèle de trois manières :

1. Le Test 1D (Le Long Couloir) : Ils ont simulé une bande longue et étroite de sol gelé. Ils ont comparé leurs résultats à une « référence absolue » mathématique connue (une solution exacte). Leur modèle correspondait parfaitement, prouvant qu'il pouvait gérer correctement le gel.
2. Le Test d'Écoulement (La Rivière) : Ils ont simulé l'écoulement de l'eau à travers le matériau sans gel. Là encore, leurs résultats correspondaient parfaitement aux mathématiques connues.
3. Le Test Complexe (L'Île Gelée) : C'était le grand défi. Ils ont créé une simulation 2D d'une « île » de glace gelée à l'intérieur d'une éponge plus chaude remplie d'eau. Ils ont comparé leurs résultats à une méthode standard très populaire appelée Méthode des Éléments Finis (MEF).
   • Le Résultat : Leur modèle correspondait à la méthode standard lorsque les conditions étaient calmes.
   • Le Super-pouvoir : Lorsqu'ils ont augmenté la pression de l'eau pour faire couler l'eau très vite, la méthode standard (MEF) s'est perdue et a échoué. Leur nouveau modèle de « quartier » a continué de fonctionner parfaitement, gérant l'écoulement à grande vitesse et la fonte de la glace sans broncher.

Pourquoi Cela Compte (Selon l'Article)

Les auteurs expliquent que cette simulation réussie est une première étape cruciale. En suivant avec précision comment la chaleur et l'eau se déplacent ensemble tandis que la glace se forme et fond, ils posent les fondations d'un modèle plus complexe. Ce futur modèle pourrait nous aider à comprendre :

• Comment se comporte le pergélisol (sol gelé en permanence).
• Le phénomène du gélivage, où le sol gelé se soulève et endommage les routes, les bâtiments et les mines.

En bref, l'article présente un nouveau système robuste de « surveillance de quartier » pour les mathématiques, capable de gérer les frontières mouvantes et désordonnées de l'eau gelant dans le sol mieux que les anciennes méthodes, en particulier lorsque l'eau se déplace rapidement.

Résumé technique

Résumé technique : Analyse du transfert de chaleur et de l'écoulement de l'eau avec changement de phase dans des milieux poreux saturés par la peridynamie basée sur les liens

Énoncé du problème
L'article aborde les défis physiques et mathématiques complexes associés à la modélisation du transport de chaleur et de masse dans des milieux poreux saturés où se produit un changement de phase (par exemple, le gel et le dégel des sols). Ces processus impliquent de fortes non-linéarités physiques, telles que des changements rapides des propriétés physiques et mécaniques entre les phases, et des discontinuités géométriques aux interfaces entre phases. Les méthodes numériques traditionnelles basées sur des formulations différentielles locales, telles que la méthode des éléments finis (MEF), peinent à gérer de manière robuste ces discontinuités et ces fortes non-linéarités, en particulier lorsqu'il s'agit d'interfaces évolutives et de forts gradients. Les auteurs visent à développer une approche non locale capable de prédire avec précision l'emplacement des interfaces de phase et de calculer les distributions de température et de pression dans des conditions d'écoulement de l'eau piloté par la pression.

Méthodologie
Les auteurs développent un modèle non local basé sur la peridynamie basée sur les liens (PD). Dans ce cadre, les équations aux dérivées partielles classiques sont remplacées par des équations intégrales-différentielles, permettant le traitement des discontinuités sans singularités.

1. Formulation théorique :

   • Équations gouvernantes : L'étude formule la conservation de la masse pour l'eau (en considérant les phases liquide et glace) et la conservation de l'énergie (prenant en compte la conduction thermique, la convection et la chaleur latente). Ces équations sont dérivées des équations locales classiques (loi de Darcy, loi de Fourier et bilan énergétique) et reformulées dans un contexte PD non local.
   • Modélisation du changement de phase : Le modèle intègre une courbe caractéristique de gel du sol (relation de type Weibull) pour définir la saturation en eau liquide en fonction de la température. La chaleur latente est traitée comme un terme de capacité thermique spécifique dépendant de la température.
   • Classification des liens : Le milieu poreux est discrétisé en particules connectées par des « liens de transport » (liens-t). Les liens sont classés en fonction des phases des particules connectées :
     • Liens-t liquides : Transfèrent à la fois la chaleur et l'eau.
     • Liens-t interfaciaux : Relient différentes régions (liquide/solide) et transfèrent à la fois la chaleur et l'eau.
     • Liens-t solides (imperméables) : Relient les particules au sein de la région solide (glace). Ceux-ci sont définis comme imperméables, transférant uniquement la chaleur, simulant ainsi efficacement la perméabilité négligeable du sol gelé.
   • Couplage : Le modèle couple l'écoulement de l'eau piloté par la pression avec le transfert de chaleur transitoire. Le flux d'eau est calculé sur la base des différences de potentiel d'écoulement à travers les liens, et ce flux entraîne le transport de chaleur convectif.

2. Implémentation numérique :

   • Le domaine est discrétisé en une grille uniforme de particules.
   • L'équation de conservation de la masse est résolue pour déterminer le champ de pression et les vecteurs de flux d'eau.
   • L'équation de conservation de l'énergie est résolue pour mettre à jour le champ de température, intégrant à la fois les termes conductifs et convectifs dérivés des interactions entre liens.
   • Une attention particulière est portée au traitement des termes d'auto-interaction (où une particule interagit avec elle-même dans la formulation intégrale) en utilisant des techniques de moyennage pour les plus proches voisins.

Contributions clés

• Couplage novateur : Ce travail présente la première formulation peridynamique basée sur les liens qui couple le transfert de chaleur transitoire avec le changement de phase et l'écoulement de l'eau piloté par la pression dans des milieux poreux saturés. Les développements PD précédents avaient abordé la conduction thermique ou l'écoulement des fluides séparément, ou les avaient couplés sans changement de phase.
• Suivi d'interface : La formulation non locale gère naturellement l'interface mobile entre les phases liquide et solide sans nécessiter de reconstruction explicite du maillage ou d'algorithmes de suivi d'interface typiques des méthodes locales.
• Gestion des forts gradients : La formulation est conçue pour rester stable dans des conditions de forts gradients de pression et de fortes vitesses d'eau, des scénarios où les méthodes différentielles locales souffrent souvent d'instabilité numérique ou de divergence.

Résultats et vérification
La méthodologie proposée a été vérifiée par trois cas d'essai distincts :

1. Transfert de chaleur 1D avec changement de phase : Les résultats PD pour un front de gel dans un domaine semi-infini ont été comparés à une solution analytique existante. Le modèle PD a montré un accord très étroit, tant pour la position du front de phase que pour la distribution de température au fil du temps.
2. Transfert de chaleur 1D avec écoulement d'eau : Un problème convectif-conductif a été résolu et comparé à une solution analytique. Les résultats ont démontré une mise en œuvre précise du transport de chaleur convectif.
3. Problème couplé 2D (Inclusion de gel) : Un problème 2D impliquant le dégel d'une inclusion gelée dans un milieu poreux sous écoulement piloté par la pression a été simulé. Les résultats PD ont été comparés aux simulations par la méthode des éléments finis (MEF) (utilisant Comsol Multiphysics).
   • Température : Les résultats PD et MEF ont montré un bon accord pour les distributions de température.
   • Pression : Le modèle PD a suivi avec succès les variations de pression lors de la fonte de l'inclusion gelée. Notamment, le modèle PD a capturé des valeurs de pression négatives extrêmes au sein de l'inclusion gelée (en raison des changements de densité dans des conditions de volume presque constant) que le modèle MEF n'a pas réussi à suivre, entraînant une divergence de la solution MEF à des gradients de pression plus élevés.
   • Écoulement à haute vitesse : Un quatrième cas d'essai impliquant une « fissure » dans un corps gelé soumis à un écoulement d'eau chaude à haute vitesse a démontré que le modèle PD pouvait gérer le transfert de chaleur dominé par la convection et le changement de phase rapide là où les solutions MEF devenaient instables ou divergentes.

Signification et affirmations
L'article affirme que l'approche peridynamique basée sur les liens développée fournit une alternative précise et robuste aux méthodes locales pour la modélisation des processus thermo-hydrauliques couplés avec changement de phase. Les exercices de vérification démontrent que la méthodologie capture correctement la physique du transport de chaleur et d'eau en présence de discontinuités évolutives.

Les auteurs positionnent ce travail comme une étape cruciale vers un modèle complet thermo-hydro-mécanique (THM). Un tel modèle permettrait de décrire les comportements hydrologiques complexes dans les sols de pergélisol et des phénomènes comme le soulèvement par le gel, qui constituent des défis majeurs en génie civil et minier. L'article note modestement que, bien que le modèle actuel gère efficacement le couplage hydro-thermique, des travaux futurs seront nécessaires pour intégrer la déformation mécanique et les implications physiques spécifiques des gradients de pression extrêmes observés dans les zones gelées.