Finite elements for the space approximation of a differential model for salts crystallization

Dit artikel presenteert een ruimte-tijdmethode op basis van eindige elementen en een impliciet-expliciet tijdsstappenprocedures voor de simulatie van zoutkristallisatie in steen, waarbij de stabiliteit, convergentie en prestaties in twee en drie dimensies worden geanalyseerd en vergeleken met eerdere eendimensionale modellen.

De kern

Zoutkristallen in steen: Een wiskundig avontuur om oude monumenten te redden

Stel je voor dat een oud stenen gebouw, zoals een kerk of een standbeeld, niet alleen wordt aangetast door de wind en de regen, maar ook door iets onzichtbaars: zout.

Wanneer het regent of de grondvochtigheid omhoog trekt, neemt de steen water op. Dit water bevat vaak opgelost zout. Als het water later verdampt (bijvoorbeeld door de zon of de wind), blijft het zout achter. Maar hier wordt het lastig: het zout verandert van een vloeibare vorm in harde kristallen.

De analogie: De opgezwollen spons
Stel je de steen voor als een enorme, microscopische spons.

1. De spons (de steen): Heeft van nature veel gaatjes (poriën).
2. Het water: Vult deze gaatjes.
3. Het zout: Lost op in het water en zwemt rond.
4. De kristallisatie: Zodra het water verdampt, groeien de zoutdeeltjes uit tot grote, harde kristallen.

Het probleem is dat deze kristallen ruimte nodig hebben. Ze groeien in de kleine gaatjes van de steen. Als ze te groot worden, duwen ze de wanden van de gaatjes uit elkaar. Het is alsof je een ballon probeert op te blazen in een glazen potje; op een gegeven moment barst de pot. Zo barst de steen uiteen, wat leidt tot afbrokkeling en schade aan ons cultureel erfgoed.

Wat doen de onderzoekers in dit artikel?

De auteurs van dit paper (een groep wiskundigen en ingenieurs uit Italië) hebben een computerprogramma ontwikkeld om precies te voorspellen hoe dit proces werkt.

1. Van een rechte lijn naar een 3D-blok
Vroeger hadden wetenschappers een simpele manier om dit te berekenen: ze keken alleen naar een rechte lijn (één dimensie). Alsof ze alleen naar de hoogte van een baksteen keken.

• Het probleem: In het echt zijn gebouwen 3D-voorwerpen. Zout kan ook van opzij binnendringen, niet alleen van onderen.
• De oplossing: De onderzoekers hebben hun model uitgebreid naar 2D (vlak) en 3D (ruimte). Ze gebruiken een techniek genaamd Finite Element Method (FEM).
  • De analogie: Stel je voor dat je een complexe vorm (zoals een standbeeld) wilt meten. In plaats van één rechte lijn te trekken, snijd je het object op in duizenden kleine, perfecte blokjes (tetraëders). Voor elk blokje berekenen ze de wiskunde. Dit geeft een veel realistischer beeld van wat er in de steen gebeurt.

2. Twee fases in het spel
Het model simuleert twee fases, net als een experiment in een laboratorium:

• Fase 1: Opdrinken (Imbibition): De steen staat in zout water. Het water zuigt omhoog (zoals een theezakje in heet water) en brengt zout mee.
• Fase 2: Drogen (Drying): Het water verdampt aan de oppervlakte. Het zout kristalliseert en de steen droogt uit.

3. De "Gevoeligheidsanalyse" (Het testen van de regels)
De onderzoekers wilden weten: Wat als we de regels een beetje veranderen?
Ze hebben gekeken naar drie belangrijke knoppen:

• Hoe groot zijn de zoutkristallen?
• Hoe snel groeien ze?
• Hoe snel verdampt het water?

Ze hebben de computer laten zien wat er gebeurt als ze deze knoppen 10% harder of zachter zetten. Het goede nieuws: het model is zeer robuust. Zelfs als je de instellingen een beetje verandert, blijft het gedrag van de steen redelijk stabiel. Dit betekent dat hun voorspellingen betrouwbaar zijn, zelfs als we niet 100% precies weten hoe snel het zout groeit in de natuur.

4. Waarom is dit belangrijk?
Met dit nieuwe, krachtige 3D-model kunnen restaurateurs en ingenieurs nu beter voorspellen:

• Waar in een muur het gevaarlijkste is voor zoutschade.
• Of een beschermende behandeling (zoals een lak of coating) echt werkt.
• Hoe lang een monument nog standhoudt voordat het begint te barsten.

Samenvattend:
De onderzoekers hebben een oude, simpele wiskundige formule (voor een rechte lijn) omgebouwd tot een modern, 3D-computerspel. In dit spel spelen ze de rol van een "digitale steen", waarbij ze elke seconde van het proces simuleren: van het opdrinken van zout water tot het barsten van de steen door kristalgroei. Hierdoor kunnen we in de toekomst betere beslissingen nemen om onze historische gebouwen te redden.

Technische samenvatting

Titel: Finite Elementen voor de Ruimtelijke Benadering van een Differentieel Model voor Zoutkristallisatie

1. Probleemstelling en Context

Het artikel richt zich op de degradatie van steenartefacten uit het culturele erfgoed, veroorzaakt door zoutkristallisatie in poreuze materialen. Wanneer stenen blootstaan aan atmosferische factoren en vocht, dringen zoutoplossingen via capillaire werking het materiaal binnen. Bij verdamping of verzadiging precipiteren de zouten, wat leidt tot kristalvorming. Dit proces kan plaatsvinden op het oppervlak (efflorescentie) of, wat schadelijker is, binnen de poriën (subflorescentie).

De kristallisatie binnen de poriën vermindert het beschikbare poreuze volume, verstoort het vloeistofnetwerk en vertraagt het watertransport, wat uiteindelijk leidt tot structurele schade. Bestaande wiskundige modellen (zoals in referentie [5]) zijn beperkt tot één dimensionale (1D) ruimtelijke representaties. Hoewel deze geschikt zijn voor het modelleren van verticale stijging, zijn ze onvoldoende voor realistische, complexe geometrieën in 2D en 3D. Het doel van dit werk is het uitbreiden van dit model naar multidimensionale domeinen en het ontwikkelen van een robuust numeriek oplosmiddel.

2. Wiskundig Model

Het model beschrijft de gekoppelde processen van vochttransport, zoutmigratie en kristallisatie in een poreus medium $\Omega$. De toestand wordt bepaald door vier variabelen:

• $\theta_l$: Het volumefractie van vloeibaar water.
• $c_i$: De concentratie van opgeloste zoutionen.
• $c_s$: De concentratie van gekristalliseerd zout.
• $n$: De porositeit van het materiaal.

De porositeit $n$ is geen onafhankelijke variabele maar koppelt aan $c_s$ via een constitutieve relatie ($n = n_0 - \gamma c_s$), wat de verstopping van poriën door kristalgroei simuleert.

De besturingsevergelijkingen vormen een niet-lineair systeem van partiële differentiaalvergelijkingen (PDE's):

1. Vochttransport: Een niet-lineaire diffusie vergelijking gedreven door capillaire krachten, waarbij de permeabiliteit afhangt van de verzadigingsgraad.
2. Zouttransport: Een convectie-diffusie vergelijking waarbij adveectie wordt veroorzaakt door de vochtkracht en diffusie door concentratiegradiënten.
3. Kristallisatiekinetiek: Een reactieterm die fungeert als een sink voor opgeloste ionen en een bron voor de vaste fase, afhankelijk van de lokale concentratie en beschikbare ruimte.

Het experimentele scenario omvat twee fasen:

• Opname (Imbibition): Wateropname via de onderkant (Dirichlet-voorwaarden) en uitwisseling met de lucht aan de bovenkant (Robin-voorwaarden).
• Droging: Geen toevoer van water of zout; verdamping aan de oppervlakken.

3. Methodologie

De kern van de bijdrage is de overgang van een 1D eindige-differentie (FD) methode naar een Finite Element Method (FEM) voor ruimtelijke discretisatie, gekoppeld aan een impliciet-expliciet tijdsstappen-algoritme.

• Ruimtelijke Discretisatie: Er wordt gebruikgemaakt van lineaire eindige elementen op een conformerend mesh (driehoeken in 2D, tetraëders in 3D). Dit maakt het mogelijk om complexe geometrieën en heterogene domeinen nauwkeurig te modelleren.
• Tijdsdiscretisatie: Een hybride aanpak wordt gebruikt:
  • Voor de niet-lineaire diffusie en transporttermen wordt een impliciete benadering gebruikt voor stabiliteit.
  • Voor de reactietermen (kristallisatie) wordt een expliciete benadering toegepast.
• Zwakte Formulering: De PDE's worden omgezet naar een zwakke vorm voor de implementatie in de FEniCS-computingomgeving.
• Validatie: De nieuwe FEM-oplosser wordt vergeleken met de bestaande 1D FD-methode (die als referentie dient) en onderworpen aan een gevoeligheidsanalyse.

4. Belangrijkste Bijdragen

1. Multidimensionale Uitbreiding: Het eerste model dat de zoutkristallisatie-dynamica uitbreidt van 1D naar 2D en 3D, waardoor realistischere simulaties van steenstructuren mogelijk zijn.
2. Ontwikkeling van een FEM-oplosser: Een stabiel en convergent numeriek schema is ontwikkeld en geïmplementeerd voor het simuleren van complexe geometrieën.
3. Gevoeligheidsanalyse: Een lokale gevoeligheidsanalyse is uitgevoerd voor de parameters $\gamma$ (specifiek volume kristallen), $K_s$ (kristallisatiesnelheid) en $K_w$ (uitwisselingscoëfficiënt). De analyse toont aan dat het model robuust is tegen kleine variaties in deze parameters.
4. Convergentieanalyse: Een experimentele analyse bevestigt de convergentieorden van de methode.

5. Numerieke Resultaten

De auteurs hebben simulaties uitgevoerd in 1D, 2D en 3D:

• 2D Simulaties: Een lange opnamefase (bijna 4444 uur) toonde aan dat de FEM-resultaten goed overeenkomen met de 1D FD-benchmark. De simulaties visualiseerden de vochtdoorlatingsfronten en de langzamere kristallisatie, met een waarneembare afname van de porositeit aan de bovenkant van de steen.
• 3D Simulaties: Een realistische experimentele opstelling (een prismatische steen) werd gesimuleerd voor zowel opname als droging. De resultaten tonen duidelijk dat het vochtfront ($\theta_l$) sneller beweegt dan het kristallisatiefront ($c_s$). Tijdens de drogingsfase verdwijnt het vrije vocht snel, terwijl de hoeveelheid neergeslagen zout en de porositeit grotendeels stabiel blijven.
• Convergentie: De experimentele convergentieanalyse toont een tijdsorde van $O(\Delta t)$ en een ruimtelijke orde van $O(h^2)$ voor de meeste variabelen ($c_i, c_s, n$). Voor de vochtvariabele $\theta_l$ is de convergentie iets trager ($O(h^2 \log(1/h))$) als gevolg van de niet-lineariteit in de diffusievergelijking. De conclusie is dat de FEM-methode superieure prestaties levert ten opzichte van de FD-methode voor dit type probleem.

6. Betekenis en Toekomstperspectief

Deze studie biedt een cruciaal instrument voor het begrijpen en voorspellen van de degradatie van cultureel erfgoed door zoutkristallisatie. Door over te stappen op 2D en 3D modellen kunnen conservatoren en ingenieurs realistischere scenario's simuleren voor complexe bouwmaterialen, in plaats van te vertrouwen op vereenvoudigde 1D benaderingen.

De robuustheid van het model en de stabiliteit van de FEM-oplosser maken het geschikt voor toekomstige uitbreidingen, zoals het modelleren van mechanische schade na kristallisatiecycli en de integratie in het "Stoneverse"-platform voor bredere toepassing in het veld van het culturele erfgoed.