Numerical homogenization of a heat transfer problem with phase transitions in heterogeneous frozen soils

Dans le vaste Grand Nord gelé, le sol n'est pas simplement de la terre ; c'est une fondation structurelle qui soutient les plus grandes villes, les pipelines et les routes du monde. Pendant des siècles, ce pergélisol est resté solide, son état gelé agissant comme un socle rigide pour les constructions humaines. Cependant, à mesure que le climat se réchauffe, ce sol gelé commence à dégeler. Lorsque la glace contenue dans le sol fond, le sol perd sa résistance, entraînant un affaissement inégal et la fissuration des bâtiments. Prédire exactement jusqu'où ce dégel ira est une question de sécurité et d'ingénierie, mais le sol lui-même présente une complication majeure. Le vrai sol n'est pas un bloc de glace uniforme ; c'est un mélange chaotique de couches, de lentilles et de poches de matériaux différents, dont les propriétés changent sur des distances aussi petites que quelques centimètres. Pour simuler la fonte d'un sol aussi complexe avec une précision parfaite, un ordinateur devrait suivre chaque couche individuelle, une tâche qui nécessite plus de puissance de calcul que ce qui est actuellement disponible pour des prévisions à grande échelle et à long terme.

Pour résoudre ce problème, les chercheurs ont recours à une technique appelée homogénéisation numérique. Imaginez que vous essayiez de prédire le flux d'eau à travers une éponge. Au lieu de modéliser chaque minuscule trou et chaque fibre, vous traitez l'éponge comme un matériau unique ayant une capacité « moyenne » à laisser passer l'eau. C'est l'essence même de l'homogénéisation : remplacer une réalité désordonnée et détaillée par un modèle simplifié et moyen, assez rapide pour être exécuté par un ordinateur. Cependant, lorsque le sol fond, cette moyenne devient un piège. La limite entre le sol gelé et le sol dégelé — le « front » — ne s'arrête pas proprement au bord du bloc de calcul de l'ordinateur ; elle le traverse de part en part. Cela crée un dilemiment : si un bloc de sol est à moitié gelé et à moitié dégelé, comment l'ordinateur doit-il décider de ses propriétés ? Doit-il se comporter davantage comme de la glace, davantage comme de l'eau, ou quelque chose entre les deux ? Le choix fait ici détermine si le modèle prédit qu'un bâtiment va s'enfoncer ou rester stable.

Dans une étude récente, Timur Zakharov, de l'Université fédérale du Nord-Est, a abordé ce dilemme spécifique en testant trois façons différentes de gérer ces blocs de sol « mixtes ». La recherche s'est concentrée sur un bâtiment chauffé posé sur un sol gelé, un scénario où la chaleur provenant de la structure fait fondre lentement le sol en dessous. L'équipe a créé un jumeau numérique de cet environnement, utilisant un motif aléatoire et réaliste de conductivité du sol pour imiter le chaos naturel du terrain. Ils ont ensuite lancé la simulation sur une grille grossière, où chaque bloc de sol était beaucoup plus grand que les couches réelles de terre et de glace, forçant le front de changement de phase à trancher au milieu de ces blocs. Ils ont comparé trois méthodes distinctes pour décider de ce qui se passe à l'intérieur de ces blocs tranchés.

La première méthode, connue sous le nom d'approche par table, pré-calcule une liste de propriétés pour chaque état possible d'un bloc de sol et cherche la réponse en fonction de la température moyenne du bloc. La deuxième méthode, un schéma nodal, décide de l'état du sol aux coins des blocs et fait la moyenne de ces décisions. La troisième méthode, appelée h-lissage (h-smoothing), calcule la fraction exacte du bloc qui est fondue en se basant sur un gradient de température linéaire, mesurant essentiellement la tranche du bloc qui s'est transformée en eau. Les chercheurs ont exécuté ces simulations en deux et trois dimensions, comparant les résultats à une solution de « référence » hautement détaillée qui suivait chaque couche de sol.

Les résultats ont révélé un compromis surprenant qui remet en question la manière dont les ingénieurs pourraient choisir leurs outils. Lorsque l'objectif est de prédire la température du sol, la méthode du h-lissage s'est avérée être la plus précise. Elle a suivi le flux de chaleur avec la plus haute précision, surtout sur les grilles grossières où les autres méthodes peinaient. Cependant, lorsque les chercheurs ont examiné la position du front de fonte — la limite réelle entre le sol gelé et le sol dégelé — l'histoire change complètement. La méthode du h-lissage, bien que performante pour la température, a systématiquement surestimé la profondeur du dégel. Elle prédisait que le sol fondrait plus profondément qu'il ne le faisait réellement.

La raison de cette divergence réside dans les lois de la conservation de l'énergie. La méthode du h-lissage, dans son effort pour être précise quant à la forme de la zone de fusion, a involontairement ajouté une petite quantité de chaleur supplémentaire au système à chaque étape de temps. Au cours d'une année, ce surplus minuscule s'est accumulé, provoquant un dégel du sol plus important que celui du sol réel. En revanche, les deux autres méthodes, particulièrement l'approche par table et le schéma nodal, étaient « conservatrices ». Elles n'ajoutaient ni ne perdaient d'énergie, et en conséquence, elles ont prédit la profondeur du dégel avec beaucoup plus de précision, même si leurs lectures de température étaient légèrement moins précises. L'étude a montré que si un chercheur a besoin de connaître la température exacte du sol, le h-lissage est le choix supérieur. Mais si l'objectif est de savoir à quelle profondeur la fondation va s'enfoncer à cause du dégel, les méthodes conservatrices sont bien plus fiables.

La recherche a également exploré ce qui se passe lorsque la température à laquelle le sol fond n'est pas la même partout, mais varie selon les différentes couches. Dans ce scénario plus complexe, la méthode basée sur la table, qui utilise une recherche détaillée des propriétés, a surpassé les méthodes plus simples à deux tenseurs pour prédire la position du front de dégel. Cela suggère que pour des conditions de sol hautement variables, une carte plus détaillée des états du sol pré-calculée est nécessaire. En fin de compte, l'étude fournit un guide clair aux ingénieurs et aux scientifiques : il n'existe pas de méthode unique « idéale ». Le choix dépend entièrement de la question posée. Si la priorité est la précision thermique, une méthode l'emporte ; si la priorité est la limite physique du dégel, une autre prend la tête. En clarifiant ces différences, ce travail aide à garantir que les modèles utilisés pour protéger les bâtiments dans l'Arctique ne soient pas seulement rapides, mais aussi dignes de confiance.