Topological cell-openness index for porous materials

Dit artikel stelt een topologische index voor celopenheid ($\tau$) voor, gebaseerd op Betti-getallen, als een complementaire of alternatieve maatstaf voor gaspycnometrie voor het karakteriseren van de verhouding tussen open en gesloten cellen in poreuze materialen, terwijl tevens de correlatie met fysische grootheden en de bruikbaarheid voor het schatten van kenmerkende afmetingen wordt aangetoond.

De kern

Stel je een spons voor. Sommige sponzen zitten vol gaten die allemaal met de buitenkant verbonden zijn, waardoor water er direct doorheen kan stromen. Anderen hebben ook gaten, maar veel daarvan zitten opgesloten in het binnenste, als kleine belletjes die in glas zijn verzegeld, zodat water er niet in of uit kan.

Al lang hebben wetenschappers een standaardmethode om te meten hoe "open" een spons is. Ze noemen dit gaspyknometrie. Denk hierbij aan het blazen in de spons met een rietje. Als de lucht erin kan, is het gat "open". Als de lucht er niet in kan, is het gat "gesloten". Deze methode geeft je één getal: het percentage open ruimte. Het is de industriestandaard.

De auteurs van dit artikel, Michał Bogdan en Paweł Dłotko, merkten echter een probleem op. Stel je een spons voor waarbij 99% van de gaten open is naar de buitenkant, maar de resterende 1% eigenlijk een hoopje tiny, geïsoleerde belletjes zijn die binnenin het open netwerk zijn opgesloten. De standaardtest "blazen erin" zou zeggen: "Geweldig! Het is 99% open!" en daar bij stopten. Het mist het feit dat het open deel eigenlijk een rommelig, niet-verbonden web is in plaats van één gladde snelweg.

Om dit op te lossen, hebben de auteurs een nieuw hulpmiddel ontwikkeld genaamd de Cell-Openness Index (τ).

Het nieuwe hulpmiddel: tellen van lussen en eilanden

In plaats van alleen maar lucht in te blazen, gebruiken de auteurs een tak van de wiskunde genaamd Topological Data Analysis. Je kunt dit zien als een superslimme manier om vormen en verbindingen te tellen in een 3D-afbeelding van het materiaal.

Ze maken gebruik van een concept genaamd Betti-getallen, die ingewikkeld klinken maar eigenlijk gewoon tellers zijn voor specifieke vormen:

• Eilanden tellen (0D): Hoeveel losse brokken gaten zijn er?
• Lussen tellen (1D): Hoeveel ringen of donut-vormen kun je maken door door de gaten te lopen?
• Grotten tellen (2D): Hoeveel volledig ingesloten belletjes zijn er?

De auteurs combineren deze tellingen in hun nieuwe index, τ.

• Als τ 0 is, is het materiaal als een zak met marbles: elk gat is een apart, gesloten eiland. Niets is verbonden.
• Als τ 1 is, is het materiaal als een perfect honingraatpatroon: elk gat is verbonden met elk ander gat in één groot, open netwerk.

Waarom is dit beter dan de oude manier?

Het artikel laat zien dat terwijl de oude methode (gaspyknometrie) en de nieuwe methode (τ) meestal overeenkomen, ze soms op een zeer interessante manier van mening verschillen.

Stel je twee sponzen voor die volgens de oude methode allebei "99% open" testen.

• Spons A is een perfect, onderling verbonden web.
• Spons B ziet eruit als een web, maar is eigenlijk gemaakt van 50 aparte webben die allemaal de rand van de spons raken, maar elkaar niet raken.

De oude methode ziet beide als "99% open". De nieuwe methode (τ) ziet Spons A als "zeer open" (hoge score) en Spons B als "minder open" (lagere score), omdat het detecteert dat het netwerk is opgebroken in niet-verbonden stukken. Het is als het verschil tussen een stad met één groot snelwegennetwerk en een stad met 50 aparte doodlopende straten die toevallig allemaal de stadsgrenzen raken.

De "vingerafdruk" van het materiaal lezen

De auteurs ontdekten ook dat ze door te kijken naar hoe deze vormtellingen veranderen terwijl ze in en uitzoomen op de afbeelding (een proces genaamd filtratie), ze de fysieke grootte van de gaten kunnen raden.

Denk hierbij aan het luisteren naar een liedje. Als je het ritme en de noten kent, kun je de grootte van de instrumenten die ze spelen raden.

• Ze ontdekten dat de "pieken" en "dalen" in hun grafieken voor vormtelling corresponderen met de grootte van de gaten, de afstand tussen de gaten en de dikte van de vaste wanden ertussen.
• Dit werkte zeer goed voor materialen met gesloten, geïsoleerde gaten (zoals een blok Zwitserse kaas waarbij de gaten elkaar niet raken).
• Het was iets lastiger voor open, rommelige netwerken, maar leverde toch nuttige aanwijzingen op.

Maakt dit iets uit voor het echte leven?

De auteurs hebben getest of hun nieuwe getal (τ) kon voorspellen hoe goed een materiaal warmte of vloeistoffen verplaatst.

• Vloeistoffen (Permeabiliteit): In 2D-modellen vonden ze een zeer sterke, duidelijke relatie tussen hun nieuwe index en hoe gemakkelijk vloeistof door het materiaal stroomt.
• Warmte (Warmtegeleiding): In 3D-modellen was hun nieuwe index iets beter in het voorspellen van hoe goed warmte door het materiaal beweegt, vergeleken met de oude methode.

De conclusie

Het artikel beweert niet dat dit ziektes zal genezen of direct nieuwe raketten zal bouwen. In plaats daarvan stelt het een eenvoudige, wiskundige "tweede mening" voor voor het meten van poreuze materialen.

Als je een spons, een rots of een schuim analyseert, vertelt de oude methode je hoeveel ruimte open is. De nieuwe methode van de auteurs vertelt je hoe goed die open ruimte is verbonden. Ze suggereren dat wanneer je een hoogwaardige 3D-afbeelding van een materiaal hebt, je beide getallen moet rapporteren: het oude (voor de traditie) en het nieuwe (om de verborgen, niet-verbonden stukken te vangen die de oude methode mist).

Technische samenvatting

Technische Samenvatting: Topologische Index voor Cel-openheid bij Poreuze Materialen

Probleemstelling
Het karakteriseren van de openheid van poreuze materialen is cruciaal voor het begrijpen van transporteigenschappen, maar huidige standaarden vertrouwen sterk op gaspycnometrie. Deze experimentele methode levert een enkele parameter op, de "fractie open poriën" ($\phi_0$), die de volumefractie van poriën voorstelt die verbonden zijn met de samplegrens. Echter, $\phi_0$ heeft beperkingen: het vereist kennis van de ware dichtheid van het materiaal, heeft moeite met samples waarbij open poriën niet met het oppervlak verbonden zijn, en kan interne structurele onderbrekingen binnen een anderszins "open" fase niet oplossen (bijvoorbeeld meerdere componenten die de grens raken maar onderling niet verbonden zijn). Met de toenemende toegankelijkheid van 3D micro-imaging (bijvoorbeeld X-ray micro-CT) is er behoefte aan complementaire, op beelden gebaseerde methoden die topologische connectiviteitsinformatie vastleggen die verder gaat dan eenvoudige volumefracties.

Methodologie
De auteurs stellen een op beelden gebaseerde methode voor die Topologische Data-analyse (TDA), specifiek persistente homologie, gebruikt om de verhouding van open en gesloten cellen in binaire 2D- en 3D-gevoxeliseerde microstructuren te schatten.

1. Data-generatie en Bronnen: De studie maakt gebruik van synthetische binaire microstructuren gegenereerd via een stochastisch proces waarbij bolvormige/cirkelvormige poriën op een rooster worden geplaatst, verbonden door kanalen met variërende kansen ($p$). Vier families datasets werden gecreëerd: gedeeltelijk verbonden gaten (2D/3D), puur gesloten cellen (2D), puur open cellen (2D), en gemengde 3D-cell voor warmteoverdrachtsanalyse. Daarnaast werd een extern aangeleverde 2D-dataset van 31 poreuze afbeeldingen met bekende permeabiliteitswaarden gebruikt.
2. Numerieke Pycnometrie ($\phi_0$): Een digitaal equivalent van gaspycnometrie werd geïmplementeerd door samenhangende componenten van de poriefase die de domeingrens raken te identificeren en hun volumefractie te berekenen ten opzichte van het totale porievolume.
3. Topologische Analyse: De auteurs berekenen Betti-getallen ($\beta_0, \beta_1, \beta_2$) op sub-niveauverzamelingen van de Signed Distance Transform (SDT) van de binaire structuur.
   • $\beta_0$: Aantal samenhangende componenten.
   • $\beta_1$: Aantal onafhankelijke lussen (1D) of tunnels.
   • $\beta_2$: Aantal ingesloten holtes (in 3D).
   • De analyse filtert ruis weg door alleen kenmerken met persistentie $\ge 1.5$ te behouden.
4. De Cel-openheidsindex ($\tau$): Een nieuwe scalair index $\tau \in [0, 1]$ wordt gedefinieerd op basis van de Betti-getallen bij een specifieke filtratiewaarde ($t=-1$):
   • Voor 2D: $\tau = \frac{\beta_1}{\beta_0 + \beta_1}$
   • Voor 3D: $\tau = \frac{\beta_1 + \beta_2}{\beta_0 + \beta_1 + \beta_2}$
   • $\tau = 0$ geeft een systeem van niet-verbonden, gesloten poriën aan; $\tau = 1$ geeft een volledig verbonden, open netwerk aan.
5. Fysische Simulaties: Effectieve thermische geleidbaarheid werd berekend voor 3D-structuren door de stationaire diffusievergelijking op te lossen met stuksgewijs isotrope geleidbaarheid. Permeabiliteitsdata werd ontleend aan de externe referentiedataset.

Belangrijkste Bijdragen en Resultaten

• Definitie van $\tau$: Het artikel introduceert $\tau$ als een topologische aanvulling op $\phi_0$. Hoewel $\phi_0$ en $\tau$ sterk gecorreleerd zijn, vangt $\tau$ connectiviteitsnuances die $\phi_0$ mist. Specifiek kan $\tau$ in regimes waar $\phi_0$ verzadigt bij 1 (wat hoge openheid aangeeft), onderscheid maken tussen een echt goed verbonden netwerk en een netwerk bestaande uit meerdere componenten die individueel de grens raken maar onderling niet verbonden zijn.
• Correlatie met Fysische Eigenschappen:
  • 2D Permeabiliteit: In een dataset van 31 poreuze afbeeldingen toonde $\tau$ een sterke parabolische relatie met $\log_{10}(\text{permeabiliteit})$ ($R^2 = 0.95$), ondanks dat $\tau$ in de dataset met minder dan 1% varieerde.
  • 3D Thermische Geleidbaarheid: Over 250 synthetische 3D-structuren toonde $\tau$ een matige positieve correlatie met de trace van de effectieve thermische geleidbaarheidstensor ($\text{tr}(K)$). In alle sub-datasets was $\tau$ een iets betere voorspeller van $\text{tr}(K)$ dan $\phi_0$.
• Schaalindicatie via Betti-curven: De auteurs tonen aan dat Betti-curven ($\beta_i(t)$) afgeleid van SDT-filtraties karakteristieke geometrische lengtes kunnen herstellen:
  • Gesloten Cellen: De punten van maximale gradiënt in $\beta_0$ en $\beta_1$ schatten de gemiddelde porieradius ($r$), de halve breedte tussen poriën ($d/2$) en de halve dikte van de vaste kern ($h/2$) met hoge nauwkeurigheid ($R^2$ tot 0.95).
  • Open Cellen: Hoewel de methode werkt voor geïdealiseerde periodieke structuren, was de voorspellende kracht voor de keelbreedte ($w/2$) in ruizige, realistische open-celdatasets zwakker, waarschijnlijk door overlappende lengteschalen die topologische signatuur verduisteren.

Betekenis en Beweringen
De auteurs positioneren dit werk als een "opzettelijk reducerende" poging om te bepalen of een enkel, interpreteerbaar getal ($\tau$) de klassieke gaspycnometrie-maatstaf ($\phi_0$) kan aanvullen. Zij claimen niet dat $\tau$ pycnometrie vervangt, maar dat het een topologisch perspectief biedt dat structurele ambiguïteiten oplost die onzichtbaar zijn voor op volume gebaseerde metingen.

Het artikel beweert dat:

1. $\tau$ echte structurele informatie biedt over interne connectiviteit die $\phi_0$ niet kan oplossen.
2. Het verschil tussen $\tau$ en $\phi_0$ in sterk open systemen geen fout is, maar een kenmerk dat niet-verbonden sub-netwerken aangeeft.
3. Betti-curven een weg bieden om karakteristieke feature-groottes (porieradius, keelbreedte) direct te schatten uit topologische signatuur, hoewel dit robuuster is in gesloten-celsystemen.
4. $\tau$ correleert met transporteigenschappen (permeabiliteit en thermische geleidbaarheid), wat suggereert dat het fysiek betekenisvolle structurele signatuur codeert.

De auteurs concluderen dat waar hoogwaardige 3D-imaging beschikbaar is, $\tau$ naast $\phi_0$ gerapporteerd moet worden om een vollediger karakterisering van microstructuren van poreuze materialen te bieden.