Topological cell-openness index for porous materials

Cet article propose un indice de perméabilité topologique des cellules ($\tau$) fondé sur les nombres de Betti, en tant que métrique complémentaire ou alternative à la pycnométrie gazeuse pour caractériser les proportions de cellules ouvertes versus fermées dans les matériaux poreux, tout en démontrant sa corrélation avec des grandeurs physiques et son utilité pour estimer les tailles de caractéristiques.

L'essentiel

Imaginez que vous ayez une éponge. Certaines éponges sont pleines de trous qui tous communiquent avec l'extérieur, laissant l'eau s'écouler directement à travers. D'autres ont des trous, mais beaucoup d'entre eux sont piégés à l'intérieur, comme de petites bulles scellées dans du verre, de sorte que l'eau ne peut ni entrer ni sortir.

Depuis longtemps, les scientifiques disposent d'une méthode standard pour mesurer à quel point une éponge est « ouverte ». Ils l'appellent la pycnométrie gazeuse. Imaginez cela comme souffler dans l'éponge avec une paille. Si l'air peut entrer, le trou est « ouvert ». Si l'air ne peut pas entrer, le trou est « fermé ». Cette méthode vous donne un seul chiffre : le pourcentage d'espace ouvert. C'est la référence absolue de l'industrie.

Cependant, les auteurs de cet article, Michał Bogdan et Paweł Dłotko, ont remarqué un problème. Imaginez une éponge où 99 % des trous sont ouverts vers l'extérieur, mais où le 1 % restant est en fait un amas de minuscules bulles isolées piégées à l'intérieur du réseau ouvert. Le test standard « soufflez dedans » dirait : « Super ! C'est 99 % ouvert ! » et s'arrêterait là. Il manque le fait que la partie ouverte est en réalité un réseau désordonné et déconnecté plutôt qu'une seule autoroute fluide.

Pour résoudre ce problème, les auteurs ont créé un nouvel outil appelé l'Indice d'Ouverture des Cellules (τ).

Le Nouvel Outil : Compter les Boucles et les Îles

Au lieu de simplement souffler de l'air, les auteurs utilisent une branche des mathématiques appelée Analyse Topologique des Données. Vous pouvez y voir une façon super-intelligente de compter les formes et les connexions dans une image 3D du matériau.

Ils utilisent un concept appelé les nombres de Betti, qui semblent compliqués mais sont en réalité de simples compteurs pour des formes spécifiques :

• Compter les Îles (0D) : Combien y a-t-il de fragments séparés de trous ?
• Compter les Boucles (1D) : Combien d'anneaux ou de formes de beignets pouvez-vous former en marchant à travers les trous ?
• Compter les Grottes (2D) : Combien y a-t-il de bulles complètement enfermées ?

Les auteurs combinent ces comptes dans leur nouvel indice, τ.

• Si τ est 0, le matériau ressemble à un sac de billes : chaque trou est une île fermée et séparée. Rien ne se connecte.
• Si τ est 1, le matériau ressemble à un nid d'abeilles parfait : chaque trou est connecté à tous les autres dans un seul réseau géant et ouvert.

Pourquoi est-ce mieux que l'ancienne méthode ?

L'article montre que, bien que l'ancienne méthode (pycnométrie gazeuse) et la nouvelle méthode (τ) s'accordent généralement, elles divergent parfois d'une manière très intéressante.

Imaginez deux éponges qui toutes deux donnent un résultat de « 99 % ouvertes » selon l'ancienne méthode.

• L'Éponge A est un réseau parfait et interconnecté.
• L'Éponge B ressemble à un réseau, mais elle est en fait composée de 50 réseaux séparés qui touchent tous le bord de l'éponge sans se toucher entre eux.

L'ancienne méthode voit les deux comme « 99 % ouvertes ». La nouvelle méthode (τ) voit l'Éponge A comme « très ouverte » (score élevé) et l'Éponge B comme « moins ouverte » (score plus bas) car elle détecte que le réseau est brisé en pièces déconnectées. C'est comme la différence entre une ville avec un seul système d'autoroutes géant et une ville avec 50 impasses séparées qui touchent toutes, par hasard, les limites de la ville.

Lire la « Empreinte Digitale » du Matériau

Les auteurs ont également découvert qu'en observant comment ces comptes de formes changent lorsqu'ils « zooment » et « dézooment » sur l'image (un processus appelé filtrage), ils peuvent deviner la taille physique des trous.

Pensez-y comme écouter une chanson. Si vous connaissez le rythme et les notes, vous pouvez deviner la taille des instruments qui les jouent.

• Ils ont constaté que les « pics » et les « creux » dans leurs graphiques de comptage de formes correspondent à la taille des trous, à la distance entre les trous et à l'épaisseur des parois solides entre eux.
• Cela a très bien fonctionné pour les matériaux avec des trous fermés et isolés (comme un bloc de fromage suisse où les trous ne se touchent pas).
• C'était un peu plus délicat pour les réseaux ouverts et désordonnés, mais cela a tout de même fourni des indices utiles.

Est-ce que cela compte pour la vie réelle ?

Les auteurs ont testé si leur nouveau chiffre (τ) pouvait prédire la capacité d'un matériau à transporter de la chaleur ou des fluides.

• Fluides (Perméabilité) : Dans des modèles 2D, ils ont trouvé une relation très forte et claire entre leur nouvel indice et la facilité avec laquelle un fluide traverse le matériau.
• Chaleur (Conductivité thermique) : Dans des modèles 3D, leur nouvel indice s'est avéré légèrement meilleur pour prédire la façon dont la chaleur se déplace à travers le matériau par rapport à l'ancienne méthode.

La Conclusion

L'article ne prétend pas que cela guérira des maladies ou construira de nouveaux fusées immédiatement. Au contraire, il propose un simple « deuxième avis » basé sur les mathématiques pour mesurer les matériaux poreux.

Si vous analysez une éponge, une roche ou une mousse, l'ancienne méthode vous dit combien d'espace est ouvert. La nouvelle méthode des auteurs vous dit à quel point cet espace ouvert est bien connecté. Ils suggèrent que chaque fois que vous disposez d'une image 3D de haute qualité d'un matériau, vous devriez rapporter les deux chiffres : l'ancien (pour la tradition) et le nouveau (pour capturer les parties cachées et déconnectées que l'ancienne méthode manque).

Résumé technique

Résumé technique : Indice de perméabilité topologique des cellules pour les matériaux poreux

Énoncé du problème
Caractériser l'ouverture des matériaux poreux est essentiel pour comprendre leurs propriétés de transport, pourtant les normes actuelles reposent fortement sur la pycnométrie gazeuse. Cette méthode expérimentale fournit un paramètre unique, la « fraction de pores ouverts » ($\phi_0$), représentant la fraction volumique des pores connectés à la frontière de l'échantillon. Cependant, $\phi_0$ présente des limites : il nécessite la connaissance de la densité réelle du matériau, éprouve des difficultés avec les échantillons où les pores ouverts ne sont pas connectés à la surface, et ne peut pas résoudre les discontinuités structurelles internes au sein d'une phase par ailleurs « ouverte » (par exemple, plusieurs composants touchant la frontière mais déconnectés les uns des autres). Avec l'accessibilité croissante de l'imagerie microscopique 3D (par exemple, la micro-tomodensitométrie aux rayons X), il existe un besoin de méthodes complémentaires basées sur l'image, capables de capturer des informations de connectivité topologique au-delà des simples fractions volumiques.

Méthodologie
Les auteurs proposent une méthode basée sur l'image utilisant l'analyse de données topologiques (TDA), spécifiquement l'homologie persistante, pour estimer la proportion de cellules ouvertes et fermées dans des microstructures binaires 2D et 3D voxelisées.

1. Génération et sources de données : L'étude utilise des microstructures binaires synthétiques générées via un processus stochastique impliquant des pores sphériques/circulaires sur une grille, connectés par des canaux avec des probabilités variables ($p$). Quatre familles de jeux de données ont été créées : des trous partiellement connectés (2D/3D), des cellules purement fermées (2D), des cellules purement ouvertes (2D) et des cellules mixtes 3D pour l'analyse du transfert de chaleur. De plus, un jeu de données 2D fourni de l'extérieur, composé de 31 images poreuses avec des valeurs de perméabilité connues, a été utilisé.
2. Pycnométrie numérique ($\phi_0$) : Un analogue numérique de la pycnométrie gazeuse a été mis en œuvre en identifiant les composantes connexes de la phase poreuse touchant la frontière du domaine et en calculant leur fraction volumique par rapport au volume total des pores.
3. Analyse topologique : Les auteurs calculent les nombres de Betti ($\beta_0, \beta_1, \beta_2$) sur les ensembles de sous-niveaux de la transformée de distance signée (SDT) de la structure binaire.
   • $\beta_0$ : Nombre de composantes connexes.
   • $\beta_1$ : Nombre de boucles indépendantes (1D) ou de tunnels.
   • $\beta_2$ : Nombre de cavités fermées (en 3D).
   • L'analyse filtre le bruit en ne conservant que les caractéristiques ayant une persistance $\ge 1.5$.
4. L'indice d'ouverture des cellules ($\tau$) : Un nouvel indice scalaire $\tau \in [0, 1]$ est défini à partir des nombres de Betti à une valeur de filtration spécifique ($t=-1$) :
   • Pour le 2D : $\tau = \frac{\beta_1}{\beta_0 + \beta_1}$
   • Pour le 3D : $\tau = \frac{\beta_1 + \beta_2}{\beta_0 + \beta_1 + \beta_2}$
   • $\tau = 0$ indique un système de pores déconnectés et fermés ; $\tau = 1$ indique un réseau entièrement connecté et ouvert.
5. Simulations physiques : La conductivité thermique effective a été calculée pour des structures 3D en résolvant l'équation de diffusion stationnaire avec une conductivité isotrope par morceaux. Les données de perméabilité ont été extraites du jeu de données de référence externe.

Contributions et résultats clés

• Définition de $\tau$ : L'article introduit $\tau$ comme un complément topologique à $\phi_0$. Bien que $\phi_0$ et $\tau$ soient fortement corrélés, $\tau$ capture des nuances de connectivité que $\phi_0$ manque. Plus précisément, dans les régimes où $\phi_0$ se sature près de 1 (indiquant une grande ouverture), $\tau$ peut discriminer entre un réseau véritablement bien connecté et un réseau composé de plusieurs composants qui touchent individuellement la frontière mais sont mutuellement déconnectés.
• Corrélation avec les propriétés physiques :
  • Perméabilité 2D : Dans un jeu de données de 31 images poreuses, $\tau$ a montré une relation parabolique forte avec $\log_{10}(\text{perméabilité})$ ($R^2 = 0.95$), malgré une variation de $\tau$ inférieure à 1 % sur l'ensemble du jeu de données.
  • Conductivité thermique 3D : Sur 250 structures synthétiques 3D, $\tau$ a montré une corrélation positive modérée avec la trace du tenseur de conductivité thermique effective ($\text{tr}(K)$). Dans tous les sous-jeux de données, $\tau$ s'est révélé être un prédicteur légèrement meilleur de $\text{tr}(K)$ que $\phi_0$.
• Estimation de l'échelle de longueur via les courbes de Betti : Les auteurs démontrent que les courbes de Betti ($\beta_i(t)$) dérivées des filtrations SDT peuvent retrouver des longueurs géométriques caractéristiques :
  • Cellules fermées : Les points de gradient maximal dans $\beta_0$ et $\beta_1$ ont estimé avec succès le rayon moyen des pores ($r$), la demi-largeur inter-poreuse ($d/2$) et la demi-épaisseur du noyau solide ($h/2$) avec une grande précision ($R^2$ allant jusqu'à 0.95).
  • Cellules ouvertes : Bien que la méthode fonctionne pour des structures périodiques idéalisées, son pouvoir prédictif pour la largeur de gorge ($w/2$) dans des jeux de données de cellules ouvertes réalistes et bruyants était plus faible, probablement en raison de l' chevauchement des échelles de longueur qui obscurcit les signatures topologiques.

Signification et affirmations
Les auteurs présentent ce travail comme une tentative « délibérément réductrice » pour déterminer si un nombre unique et interprétable ($\tau$) peut compléter la métrique classique de pycnométrie gazeuse ($\phi_0$). Ils ne prétendent pas que $\tau$ remplace la pycnométrie, mais plutôt qu'il offre une perspective topologique qui résout des ambiguïtés structurelles invisibles pour les mesures basées sur le volume.

L'article affirme que :

1. $\tau$ fournit une information structurelle véritable concernant la connectivité interne que $\phi_0$ ne peut pas résoudre.
2. Le décalage entre $\tau$ et $\phi_0$ dans les systèmes hautement ouverts n'est pas une erreur, mais une caractéristique indiquant des sous-réseaux déconnectés.
3. Les courbes de Betti offrent une voie pour estimer directement les tailles de caractéristiques (rayon des pores, largeur de gorge) à partir de signatures topologiques, bien que cela soit plus robuste dans les systèmes à cellules fermées.
4. $\tau$ est corrélé aux propriétés de transport (perméabilité et conductivité thermique), suggérant qu'il encode des signatures structurelles physiquement significatives.

Les auteurs concluent que là où une imagerie 3D de haute qualité est disponible, $\tau$ devrait être rapporté aux côtés de $\phi_0$ pour fournir une caractérisation plus complète des microstructures de matériaux poreux.