🔢 mathematics

Dissolution of carbonate stones caused by CO2 pollutant: an erosion model

Dit artikel introduceert een nieuw wiskundig model en een numeriek algoritme op basis van de level-set-methode om de erosie van carbonaatstenen door koolzuur opgeloste CO2 in een poreus medium te simuleren.

Oorspronkelijke auteurs: Elishan Christian Braun, Gabriella Bretti, Samuele Ferri, Maria Laura Santarelli, Matteo Semplice

Gepubliceerd 2026-04-01

📖 4 min leestijd🧠 Diepgaand

CC0 1.0

Oorspronkelijke auteurs: Elishan Christian Braun, Gabriella Bretti, Samuele Ferri, Maria Laura Santarelli, Matteo Semplice

Oorspronkelijk artikel vrijgegeven aan het publieke domein onder CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ Dit is een AI-gegenereerde uitleg van het onderstaande artikel. Het is niet geschreven of goedgekeurd door de auteurs. Raadpleeg het oorspronkelijke artikel voor technische nauwkeurigheid. Lees de volledige disclaimer

Hoe CO2 oude stenen langzaam laat verdwijnen: Een wiskundig verhaal

Stel je voor dat een oud marmeren standbeeld in een stad niet wordt aangevallen door een hamer, maar door een onzichtbare, zachte regen van zure druppels. Dit is wat er gebeurt met ons cultureel erfgoed door de luchtvervuiling, en dit artikel vertelt het verhaal van hoe wetenschappers dit proces in kaart brengen met een nieuwe wiskundige formule.

Hier is het verhaal, vertaald naar alledaags taal, met een paar handige vergelijkingen.

1. Het probleem: De onzichtbare vijand

We weten dat oude gebouwen en standbeelden (zoals die in de oude Romeinse haven Ostia Antica, zie foto 1) er slechter uitzien dan ooit tevoren. Vroeger duurde het duizenden jaren voordat ze verweerden; nu gaat het in twintig jaar al.

De boosdoener is koolstofdioxide (CO2), een gas dat we allemaal uitademen en dat ook door auto's en fabrieken wordt uitgestoten. Wanneer dit gas in contact komt met regenwater, verandert het in koolzuur (een heel zwak zuur).

De analogie:
Stel je voor dat het marmer een spons is. Normaal gesproken is deze spons stevig. Maar als je er zure water op laat druppelen, lost het "lijm" op dat de steenkristallen bij elkaar houdt. De steen wordt zacht, gaat in de war en valt uiteindelijk uit elkaar. Dit noemen we erosie.

2. De oplossing: Een wiskundige "voorspeller"

De auteurs van dit artikel hebben een nieuwe wiskundig model bedacht. Dit is geen simpele schatting, maar een complexe computerrekening die precies beschrijft hoe het zure water door de poriën (de kleine gaatjes) van de steen trekt en de steen oplost.

Ze noemen dit een "reactie-transport systeem".

Transport: Hoe het water en het zuur zich verplaatsen door de steen (zoals water dat door een theezakje stroomt).
Reactie: Hoe het zuur chemisch reageert met de steen en hem oplost.

3. Hoe werkt de computer? (De "Ghost" methode)

Een van de lastigste dingen bij het simuleren van erosie is dat de vorm van de steen verandert. De randen worden steeds dunner en de steen "smelt" weg.

Om dit op de computer te berekenen, gebruiken de auteurs een slimme truc met Level-Set (een soort digitale randbewaker).

De analogie: Stel je voor dat je een tekening van een standbeeld maakt op een rooster van vierkante vakjes (zoals een schaakbord). De computer kijkt niet alleen naar de vakjes in de steen, maar ook naar de vakjes net erbuiten (de "spookpunten" of ghost points).
Door deze spookpunten slim te gebruiken, kan de computer precies berekenen waar de rand van de steen op elk moment zit, zelfs als die rand scheef loopt of een rare vorm heeft. Het is alsof de computer een onzichtbare lijn tekent die de steen scheidt van de lucht, en die lijn beweegt mee naarmate de steen wegsmelt.

4. Wat laten de proeven zien?

De auteurs hebben hun model getest in twee situaties:

Eenvoudig (1D): Alsof je naar een rechte staaf steen kijkt.
Complexe vorm (2D): Alsof je naar een blok steen kijkt.

De resultaten:

De computer liet zien hoe de vochtigheid en het zuur langzaam de diepte van de steen in trekken.
Na verloop van tijd (in de simulatie: een jaar) zag je dat de buitenkant van de steen langzaam verdween.
Ze ontdekten dat het proces in het begin langzaam gaat, maar dat de steen uiteindelijk "vol" raakt met water en de erosie versnelt.

De "Catastrofe"-test:
Ze hebben ook een extreme situatie getest: wat gebeurt er als de steen volledig onder water staat?

Het resultaat: De steen wordt razendsnel verzadigd. De erosie gaat dan veel sneller dan in de lucht. Dit laat zien hoe belangrijk het is om te weten hoe nat of droog een omgeving is.

5. Waarom is dit belangrijk voor de toekomst?

Dit model is niet alleen leuk voor wiskundigen; het is een voorspellingsmachine voor monumenten.

Scenario's: Steden kunnen straks vragen: "Wat gebeurt er met ons marmeren standbeeld als de CO2-niveaus in de stad stijgen?" of "Hoe snel zal dit gebouw verteren als het in een regenachtig klimaat staat?"
Besluitvorming: Met deze informatie kunnen restaurateurs beter plannen. Ze kunnen precies berekenen wanneer een monument moet worden schoongemaakt of beschermd, voordat het te laat is.

Conclusie

Kortom: De auteurs hebben een digitale "tijdmachine" gebouwd. In plaats van te wachten tot een standbeeld over 100 jaar weg is, kunnen ze nu in de computer zien hoe het er over 10, 50 of 100 jaar uit zal zien, afhankelijk van hoe vervuild de lucht is.

Het is een stap in de richting van het redden van onze geschiedenis, één wiskundige berekening per keer.

Probleemstelling

Het artikel adresseert een kritiek probleem op het gebied van het behoud van cultureel erfgoed: de chemische aantasting van carbonate stenen (zoals marmer en kalksteen) door atmosferische verontreiniging.

Achtergrond: Historische gebouwen en sculpturen weatheren (verwering) sneller dan ooit tevoren, mede door de toename van koolstofdioxide ( $CO_2$ ) en andere zure gassen in de atmosfeer.
Mechanisme: $CO_2$ reageert met water in de poriën van het gesteente tot koolzuur ( $H_2CO_3$ ). Dit zuur lost de calciet ( $CaCO_3$ ) op, wat leidt tot een afname van de stevigheid en uiteindelijk tot erosie van het materiaal.
Beperkingen van bestaande modellen: Bestaande methoden, zoals de statistische erosie-index van Lipfert, zijn te globaal. Ze geven gemiddelde schattingen op basis van neerslag en depositie, maar houden geen rekening met de specifieke materiaaleigenschappen, de tijdsafhankelijke evolutie van degradatie, of lokale microklimaten. Ze missen ook vaak de specifieke impact van $CO_2$ .

Methodologie

De auteurs introduceren een nieuw wiskundig model en een numeriek algoritme om dit proces te simuleren.

1. Wiskundig Model
Het model is een niet-lineair reactie-transportstelsel in poreuze media, gebaseerd op Darcy-stroming. Het beschrijft drie gekoppelde processen:

Transport van vocht: De beweging van water door de poreuze matrix, beschreven door een vergelijking die lijkt op de Richards-vergelijking.
Transport van koolzuur: Diffusie en convectie van koolzuur ( $H_2CO_3$ ) in het water, gevolgd door een chemische reactie met het gesteente.
Verandering van porositeit: De chemische reactie ( $CaCO_3 + H_2CO_3 \rightarrow Ca^{2+} + 2HCO_3^-$ ) verhoogt de porositeit ( $n$ ) van het materiaal. Wanneer de porositeit een bepaalde drempelwaarde ( $n_{max}$ ) bereikt, treedt materiaalverlies op.

Het systeem wordt gedefinieerd door partiële differentiaalvergelijkingen (PDE's) voor de waterfractie ( $\theta$ ), de concentratie koolzuur ( $c_a$ ) en de porositeit ( $n$ ). Er worden twee benaderingen voor de absorptiefunctie $B(s)$ (die capillaire druk koppelt aan vochtgehalte) getest:

Een symmetrische functie ( $B'$ ) met 3 parameters.
Een asymmetrische functie ( $B'_{kP}$ ) die permeabiliteit en capillaire druk apart beschrijft met meer parameters voor een gedetailleerdere fysieke weergave.

2. Numeriek Algoritme
Om de bewegende grens (erosiefront) nauwkeurig te volgen, gebruiken de auteurs een combinatie van methoden:

Finite Difference: Voor discretisatie binnen het domein.
Level-Set Methode: Om de grens van het materiaal te definiëren ( $\phi(x) = n(x) - n_{max}$ ).
Ghost-Point Methode: Voor de behandeling van randvoorwaarden op onregelmatige grenzen die niet perfect op het rooster vallen. Dit maakt het mogelijk om complexe geometrieën (zoals zuilen of reliëfs) te simuleren zonder dat het rooster hoeft te worden aangepast.
Tijdsintegratie: Een impliciete Euler-methode wordt gebruikt om grote tijdstappen mogelijk te maken en stabiliteit te garanderen.

Belangrijkste Bijdragen

Nieuw Wiskundig Kader: De ontwikkeling van een gekoppeld reactie-transportmodel specifiek voor de oplossing van carbonate stenen door $CO_2$ -gebaseerde verontreiniging.
Dynamische Erosiefront: Het model simuleert niet alleen de chemische concentraties, maar ook de fysieke terugtrekking van het steenoppervlak in de tijd.
Flexibiliteit: Het model kan worden gekalibreerd met real-time omgevingsdata (temperatuur, relatieve vochtigheid, $CO_2$ -concentratie) en materiaaleigenschappen (porositeit, permeabiliteit).
Numerieke Implementatie: Een robuust algoritme dat complexe geometrieën en bewegende grenzen aankan via de level-set ghost-point techniek.

Resultaten

De auteurs hebben het model getest in 1D en 2D scenario's:

Validatie: De numerieke resultaten komen overeen met literatuurwaarden voor erosiesnelheden (ongeveer 2 mm per 100 jaar voor suburbane gebieden).
Vergelijking Modellen: De resultaten van de symmetrische en asymmetrische absorptiefuncties waren zeer vergelijkbaar voor de gebruikte dataset, maar de asymmetrische versie biedt meer flexibiliteit voor toekomstige verfijning.
Erosiefront: De simulaties tonen een realistische voortgang van het erosiefront. De porositeit neemt toe totdat het materiaal breekt, terwijl de concentratie koolzuur afneemt door de chemische reactie.
Catastrofale Scenario's: In een "catastrofaal" scenario (volledige onderdompeling in water) werd een veel snellere erosiesnelheid waargenomen ( $8.251 \cdot 10^{-4}$ cm/uur versus $2.359 \cdot 10^{-6}$ cm/uur in standaardomstandigheden), wat de sensitiviteit van het model voor omgevingscondities aantoont.
Convergentie: Het numerieke schema toont een globale orde van nauwkeurigheid van 1, zoals verwacht voor een eerste-orde schema in tijd en ruimte.

Betekenis en Toekomstperspectief

Praktische Toepassing: Dit model biedt een krachtig hulpmiddel voor conservatoren om de levensduur van monumenten te voorspellen en behoudsstrategieën te ontwikkelen die specifiek zijn voor de locatie (bijv. rekening houdend met lokale vervuilingsniveaus of klimaatverandering).
Voorspellend Vermogen: Het kan verschillende scenario's simuleren, zoals drogere versus natere klimaten of hogere $CO_2$ -niveaus, om te zien hoe deze factoren de degradatie versnellen of vertragen.
Toekomstige Ontwikkelingen: De auteurs plannen om het model uit te breiden naar 3D, zwaartekrachteffecten (capillaire stijging) toe te voegen, en mechanismen voor korstvorming (afzetting) te integreren om het wisselspel tussen erosie en depositie beter te begrijpen. Ook wordt gewerkt aan een hogere orde van nauwkeurigheid voor de numerieke oplossing.

Kortom, dit onderzoek levert een fundamentele bijdrage aan de kwantitatieve analyse van de chemische verwering van steen, verschuivend van statische statistieken naar dynamische, procesgestuurde modellen.