Predicting cement microstructure and mechanical properties in hydrating cement paste with a Phase-Field model

In dit artikel wordt een aangepast faseveldmodel gepresenteerd dat de microstructuurontwikkeling en mechanische eigenschappen van hydraterend cement nauwkeurig voorspelt door thermodynamische inconsistenties op te lossen en zo een consistente koppeling tussen chemie, microstructuur en mechanisch gedrag mogelijk maakt.

De kern

De Digitale Bakkerij: Hoe Cement Hardt en Sterk Wordt

Stel je voor dat je een gigantische, onzichtbare bakkerij binnenstapt. In deze bakkerij wordt geen brood gebakken, maar cement. Maar dit is geen gewone bakkerij; het is een wereld op microscopisch niveau, waar deeltjes met elkaar dansen, oplossen en weer samenkomen om een stevige muur te bouwen.

Dit wetenschappelijke artikel vertelt het verhaal van een nieuwe manier om te voorspellen hoe die dans verloopt en hoe sterk het eindproduct wordt. De onderzoekers gebruiken een slim computermodel, een soort "digitale kristalbal", om te kijken wat er gebeurt als cement en water mengen.

Hier is hoe het werkt, vertaald in alledaags taal:

1. Het Probleem: De Oude Kaart is Verkeerd

Vroeger gebruikten wetenschappers simpele modellen om te voorspellen hoe cement hardt. Het was alsof ze probeerden een complex balletje te beschrijven door alleen te kijken naar de gemiddelde snelheid van de dansers. Ze wisten dat er gaten (poreusheid) in het cement zaten, maar hun modellen maakten vaak te veel gaten. Het resultaat? Ze dachten dat het cement zwakker was dan het in werkelijkheid is.

2. De Oplossing: De "Phase-Field" Methode (De Vloeibare Verwarming)

De onderzoekers hebben een nieuw model ontwikkeld dat ze een Phase-Field model noemen.

• De Analogie: Stel je voor dat je een bakje met warme chocoladestukjes (het cement) in een bad van koude melk (het water) doet.
• Het Oude Model: Zou zeggen: "De stukjes smelten en verdwijnen direct."
• Het Nieuwe Model: Kijkt naar de thermodynamica. Het ziet precies hoe de chocolade langzaam smelt, hoe de vloeistof beweegt, en hoe er nieuwe, stevige kristallen (het harde cement) ontstaan uit de vloeistof.

Het grote geheim van dit nieuwe model is dat het twee dingen beter doet dan de oude methoden:

1. Het herkent het verschil: Het weet precies wanneer iets oplost en wanneer iets neerslaat (vast wordt). In oude modellen werden deze processen vaak verward.
2. Het is natuurgetrouw: Het model zorgt ervoor dat de overgangen tussen de verschillende materialen glad zijn, net als in de echte wereld, in plaats van dat het eruitziet als een pixelated videogame (zoals bij oudere methoden).

3. De Dans van de Deeltjes

In hun digitale wereld laten ze zien hoe het proces verloopt:

• Stap 1: De harde cementkorrels beginnen te smelten in het water.
• Stap 2: De vloeistof wordt verzadigd met de "opgeloste" cementdeeltjes.
• Stap 3: Zodra er genoeg is, beginnen deze deeltjes weer vast te worden en vormen ze een nieuw, stevig netwerk (de "gel").
• Het Resultaat: Na verloop van tijd vullen deze nieuwe kristallen de gaten op. Het model laat zien hoe dit netwerk groeit, net als een schuim dat steeds dichter wordt, totdat het een stevige, doorlopende structuur vormt.

4. De Sterkte-Test: Van Microscopisch naar Macrowereld

Zodra het model de nieuwe structuur heeft getekend, doen de onderzoekers alsof ze erop gaan duwen en trekken.

• Ze berekenen hoe de spanning door het materiaal loopt.
• De Verrassing: Omdat hun model minder "gaten" (poreusheid) voorspelt dan de oude methoden, blijkt het cement veel sterker te zijn.
• Het is alsof je een muur bouwt: als je oude methoden dachten dat er veel lege ruimtes in de stenen zaten, dachten ze dat de muur zou instorten. Maar omdat het nieuwe model ziet dat die ruimtes eigenlijk gevuld zijn met stevige kristallen, voorspelt het dat de muur stevig blijft staan.

5. Waarom is dit belangrijk voor ons?

Dit klinkt misschien als pure theorie, maar het heeft grote gevolgen:

• Duurzamer bouwen: Als we precies weten hoe cement werkt, kunnen we minder cement gebruiken voor even sterke gebouwen. Dat scheelt enorm veel CO2-uitstoot.
• Betere voorspellingen: We kunnen nu beter voorspellen hoe lang een brug of een dam meegaat, zonder dat we decennia moeten wachten om het te zien.
• De "Digitale Tweeling": In de toekomst kunnen ingenieurs een computermodel draaien van hun specifieke betonmix, zien hoe het hardt, en de sterkte testen voordat er ook maar één baksteen wordt gelegd.

Kortom:
De onderzoekers hebben een nieuwe, slimme "bril" opgezet om naar het binnenste van cement te kijken. Met deze bril zien ze dat het cement veel steviger en minder poreus is dan we dachten. Ze hebben de digitale kaart van de bouw van onze wereld opnieuw getekend, en die kaart is veel nauwkeuriger dan ooit tevoren.

Technische samenvatting

Titel: Voorspelling van cementmicrostructuur en mechanische eigenschappen in hydraterende cementpasta met een Fase-Veldmodel

Auteurs: Alexandre Sac-Morane, Katerina Ioannidou, Manolis Veveakis, Hadrien Rattez.

---

1. Probleemstelling

Het voorspellen van de evoluerende microstructuur van hydraterend cement is essentieel voor het begrijpen en modelleren van de ontwikkeling van de mechanische eigenschappen van cementgebonden materialen. Bestaande modellen hebben echter beperkingen:

• Empirische en homogenisatiemodellen: Deze zijn vaak gebaseerd op idealisaties van de microstructuur en missen de fysieke koppeling tussen chemie en morfologie.
• Cellular Automata (CA) modellen (bijv. CEMHYD3D): Hoewel deze veel worden gebruikt, zijn ze beperkt tot een vaste voxel-resolutie en gebruiken ze semi-empirische reactiecoëfficiënten. Ze hebben de neiging om de poreuze ruimte te overschatten en kunnen geen lokale fysische beperkingen op microscopische schaal strikt naleven.
• Bestaande Fase-Veld (PF) modellen: De eerdere PF-formuleringen voor cementhydratatie bevatten fysieke inconsistenties, zoals spontane gel-precipitatie bij evenwicht en het aannemen van gelijke evenwichtsconcentraties voor oplossing en precipitatie, wat in strijd is met de gevestigde hydratatietheorie.

Het doel van dit werk is om een fysisch gebaseerd Fase-Veldmodel te ontwikkelen dat deze inconsistenties oplost en een nauwkeurige koppeling biedt tussen hydratatiechemie, microstructuurvorming en de resulterende mechanische respons.

---

2. Methodologie

A. Het Aangepaste Fase-Veld (PF) Model

De auteurs hebben een PF-model ontwikkeld dat de microstructuur-evolutie tijdens hydratatie beschrijft via drie gekoppelde vergelijkingen (opgelost met de finite-element solver MOOSE):

1. Allen-Cahn vergelijkingen: Voor de fasevariabelen van de bron (C3S, het klinker) en het product (C-S-H gel). Deze beschrijven de heterogene oplossing/precipitatie.
2. Diffusievergelijking: Voor de opgeloste stof (solute, $c$) in het water.

Kernverbeteringen ten opzichte van eerdere modellen:

• Herschreven vrije-energie potentiaal: Een nieuwe vrije-energiefunctie wordt geïntroduceerd om spontane precipitatie bij evenwicht te voorkomen.
• Distincte evenwichtsconstanten: Er worden verschillende evenwichtsconcentraties gebruikt voor de oplossing van de bron ($c_{eq, C3S}$) en de precipitatie van het gel ($c_{eq, CSH}$). Dit weerspiegelt de thermodynamische realiteit dat deze waarden fundamenteel verschillend zijn.
• Tilting-functies: Externe brontermen ($E_d/E_p$) worden toegevoegd om de vrije-energie te "hellend" te maken, wat de kinetiek van oplossing en precipitatie stuurt op basis van de lokale concentratieafwijkingen van het evenwicht.

B. Simulatieopzet

• Systeem: Twee water-cementverhoudingen (w/c) werden onderzocht: 0,3 (oververzadigd) en 0,5 (onderverzadigd).
• Domein: 2D-slices van 100 µm (om rekentijd te beperken, hoewel het model 3D-uitbreidbaar is).
• Netwerk: Een meshgrootte van 0,2 µm.
• Initieel: Sferische deeltjes met een aangepaste deeltjesgrootteverdeling (minimale straal 2,4 µm om numerieke stabiliteit te garanderen).
• Validatie: De resultaten werden vergeleken met CEMHYD3D-simulaties en experimentele data uit de literatuur.

C. Mechanische Karakterisering

Na het berekenen van de microstructuur op een bepaald hydratatietijdstip, wordt een computational homogenisatie uitgevoerd:

• De microstructuur wordt onderworpen aan trek- en schuifbelastingen.
• De fases (pore/water, bron, pasta) krijgen isotrope elastische eigenschappen toegekend (zie Tabel 2 in het artikel).
• De homogeniseerde Young's modulus ($E_h$) en Schuifmodulus ($G_h$) worden berekend uit de spanning-rekrelaties.

---

3. Belangrijkste Resultaten

Microstructuur Evolutie

• Realistische Morphologie: Het model genereert realistische porositeitsniveaus en continue fasegrenzen. In tegenstelling tot CA-modellen, die vaak een "fingering"-patroon (vingervormige groei) missen of onnauwkeurig vastleggen, toont het PF-model soepelere interfaces die beter overeenkomen met experimentele waarnemingen.
• Hydratatiekinetiek: De simulaties tonen een snelle initiële oplossingsfase gevolgd door een langzamere fase, veroorzaakt door het "opsluiten" van opgeloste stof in de pasta (diffusie-beperking). Dit komt goed overeen met experimentele data.
• Volumefracties: Het PF-model voorspelt een lagere porositeitsfractie dan CEMHYD3D. CEMHYD3D neigt de poreuze ruimte te overschatten, wat leidt tot een onnauwkeurige voorspelling van de sterkte.
• Topologie: De analyse van de "twee-punts correlatiefunctie" toont aan dat het model de overgang van geïsoleerde aggregaten naar een doorlopend netwerk (percolatie) correct vastlegt.

Mechanische Eigenschappen

• Stijfheid: De voorspelde Young's en Schuifmoduli nemen toe met de hydratatietoestand.
• Vergelijking met Literatuur:
  • De resultaten van het PF-model liggen dichter bij waarden die zijn afgeleid van echte microstructuren (SEM-observaties) dan de resultaten van CEMHYD3D.
  • De CEMHYD3D-resultaten vertonen een kunstmatige verzwakking van het materiaal door de overschatting van de poriën.
  • Het PF-model levert waarden op die vaak in de bovenste bandbreedte van de literatuurdata liggen, wat wijst op een robuustere microstructuurvoorspelling.
• Invloed van w/c: De microstructuur met w/c = 0,3 is stijver dan die met w/c = 0,5, voornamelijk door de aanwezigheid van resterende, stijvere klinker-inclusies in het eerste geval.

---

4. Bijdragen en Innovatie

1. Fysische consistentie: Het introduceren van een aangepaste vrije-energie potentiaal en distincte evenwichtsconstanten lost fundamentele thermodynamische inconsistenties op in eerdere PF-modellen voor cement.
2. Koppeling Chemie-Mechanica: Het werk demonstreert een robuuste workflow die direct van chemische reacties naar microstructuur en vervolgens naar macroscopische mechanische eigenschappen springt, zonder empirische tussenstappen.
3. Validatie: Uitgebreide validatie tegen zowel numerieke benchmarks (CEMHYD3D) als experimentele data, waarbij het PF-model superieure prestaties toont in het voorspellen van porositeitsverdeling en stijfheid.
4. Openbaarheid: De code en scripts zijn beschikbaar gesteld via GitHub, wat reproduceerbaarheid en verdere ontwikkeling mogelijk maakt.

---

5. Betekenis en Toekomstperspectief

Deze studie bevestigt dat Fase-Veldmodellen een krachtig alternatief zijn voor traditionele Cellular Automata en homogenisatiemodellen voor cement. Door de microstructuur expliciet en fysisch correct te modelleren, kunnen onderzoekers de invloed van parameters zoals deeltjesgrootteverdeling en chemische samenstelling nauwkeuriger voorspellen.

Beperkingen en Toekomst:

• Huidige simulaties zijn 2D en computatieel intensief (ca. 10 uur voor een 2D-simulatie op een standaard PC).
• Chemische krimp en visco-elastisch gedrag van de pasta zijn in deze eerste versie vereenvoudigd of genegeerd.
• Toekomstig werk: Uitbreiding naar 3D om complexe porienetwerken en tortuositeit beter te vangen, en koppeling met andere fysica zoals temperatuur, spanning en vloeistoftransport voor een volledig multiphysica-benadering.

Concluderend biedt dit model een veelbelovende basis voor het verbeteren van het begrip van cementhydratatie op microscale schaal en kan het dienen als een complementair hulpmiddel naast traditionele 28-dagen testprotocollen.