🔢 mathematics

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

Questo articolo presenta un nuovo modello matematico e un algoritmo numerico basato sul metodo level-set per descrivere l'erosione delle pietre carbonatiche causata dalla dissoluzione della matrice porosa dovuta all'acido carbonico.

Autori originali: Elishan Christian Braun, Gabriella Bretti, Samuele Ferri, Maria Laura Santarelli, Matteo Semplice

Pubblicato 2026-04-01

📖 4 min di lettura🧠 Approfondimento

CC0 1.0

Autori originali: Elishan Christian Braun, Gabriella Bretti, Samuele Ferri, Maria Laura Santarelli, Matteo Semplice

Articolo originale dedicato al pubblico dominio sotto CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ Questa è una spiegazione generata dall'IA dell'articolo qui sotto. Non è stata scritta né approvata dagli autori. Per precisione tecnica, consulta l'articolo originale. Leggi il disclaimer completo

🏛️ Il Nemico Invisibile: Come l'aria "morde" la pietra

Immaginate che le statue antiche, i templi e i palazzi storici siano come dei giganti di pietra che dormono da migliaia di anni. Per secoli sono stati al sicuro, ma negli ultimi decenni qualcosa è cambiato. L'aria che respiriamo oggi non è più quella di una volta: è piena di "sporcizia" invisibile, in particolare di anidride carbonica (CO2).

Questo studio è come un filmato al rallentatore che ci mostra cosa succede quando questa aria inquinata incontra la pietra.

1. La Pioggia Acida è come un "Aceto"

Sapete cosa succede se versate dell'aceto su un guscio d'uovo? Il guscio inizia a frizzare e sciogliersi.
Nella realtà, l'anidride carbonica (CO2) che esce dalle auto e dalle fabbriche si mescola con l'umidità dell'aria (la rugiada o la pioggia leggera) e crea una sostanza chiamata acido carbonico.
Per la pietra calcarea (quella usata per il marmo e il travertino), questo acido è come un aceto molto debole ma costante. Non la scioglie in un secondo, ma la "morde" lentamente, giorno dopo giorno, anno dopo anno.

2. La Pietra è come una Spugna

La pietra non è un blocco solido e perfetto come sembra. Se la guardaste al microscopio, sembrerebbe una spugna piena di buchini microscopici (i pori).
Quando l'acido carbonico entra in questi buchini, inizia a sciogliere la "colla" che tiene uniti i granelli di pietra.

Il processo: L'acido entra nei buchini ➡️ scioglie la colla ➡️ la pietra diventa più porosa ➡️ alla fine, la pietra si sgretola e sparisce.

3. I Matematici hanno creato un "Simulatore di Degrado"

Gli autori di questo studio sono dei matematici e ingegneri che hanno deciso di non aspettare 100 anni per vedere quanto si rovina una statua. Hanno creato un modello matematico, ovvero una formula complessa che funziona come un videogioco di simulazione.

Hanno scritto un codice che dice al computer:

"Ehi, immagina una statua di marmo."
"Immagina che ci sia dell'aria con un certo livello di CO2."
"Simula come l'acqua e l'acido entrano nei buchini della pietra."
"Calcola quanto velocemente la superficie della statua si assottiglia."

4. La "Linea Magica" che si muove

Il cuore del loro lavoro è tracciare quella che chiamano la "fronte di erosione".
Immaginate di avere una torta di cioccolato. Se mettete la torta in un forno, il calore entra dall'esterno verso il centro. La linea che separa la parte già cotta da quella ancora cruda è la "fronte".
In questo caso, la "fronte" è il confine tra la parte di pietra che è ancora sana e la parte che è stata "mangiata" dall'acido. Il loro modello riesce a prevedere esattamente quanto velocemente questa linea si spinge verso l'interno della pietra.

5. Cosa hanno scoperto?

Hanno fatto due tipi di esperimenti nel loro simulatore:

Scenario Normale: Come se la pietra fosse all'aperto in una città. La pietra si consuma molto lentamente (pochi millimetri ogni 100 anni), ma è un danno costante.
Scenario "Catastrofico": Come se la pietra fosse completamente sommersa in acqua satura di acido. Qui il danno è velocissimo, come se la pietra venisse attaccata da un mostro affamato.

Hanno anche scoperto che il loro modello funziona bene sia per calcolare il danno su una singola linea (1D) che su una superficie complessa (2D), come la facciata di un edificio.

🎯 Perché è importante?

Prima di questo studio, per sapere quanto una statua era rovinata, bisognava aspettare che passasse del tempo o fare stime medie. Ora, con questo modello, possiamo:

Prevedere il futuro: Se sappiamo quanto è inquinata l'aria in una città, possiamo dire quanto durerà una statua tra 50 o 100 anni.
Proteggere meglio: Possiamo capire se è meglio coprire una statua, spostarla o trattarla con prodotti specifici per bloccare l'ingresso dell'acido.
Capire i dettagli: Non è solo una stima media; il modello vede come l'umidità e l'inquinamento agiscono nei micro-buchi della pietra.

In sintesi

Questo studio è come avere una sfera di cristallo matematica che ci permette di vedere il lento "invecchiamento" delle nostre opere d'arte causato dall'inquinamento. Ci aiuta a capire che l'aria che respiriamo oggi sta scrivendo una storia di degrado sulle pietre del passato, e ci dà gli strumenti per provare a fermarla prima che sia troppo tardi.

Titolo: Dissoluzione di pietre carbonatiche causata da inquinanti CO2: un modello di erosione

1. Il Problema

Il lavoro affronta il degrado chimico delle pietre carbonatiche (come marmo e calcare), materiali fondamentali per il patrimonio architettonico e scultoreo antico. Il fenomeno è accelerato dalla presenza di inquinanti atmosferici, in particolare l'anidride carbonica ( $CO_2$ ), che in presenza di umidità forma acido carbonico ( $H_2CO_3$ ).

Contesto: Le piogge acide e l'acidificazione ambientale (pH < 5) causata da $CO_2$ , $SO_x$ e $NO_x$ stanno erodendo monumenti a un ritmo senza precedenti negli ultimi 20 anni.
Limiti degli approcci esistenti: Gli indici di erosione statistici attuali (es. formula di Lipfert) forniscono solo stime medie basate su precipitazioni e deposizioni di gas, ma non tengono conto della specificità dei materiali, dell'evoluzione temporale locale del degrado o della dinamica microclimatica specifica legata alla $CO_2$ .
Obiettivo: Sviluppare un modello matematico che descriva la cinetica di erosione della matrice porosa dovuta alla penetrazione dell'acido carbonico, permettendo previsioni più accurate basate su dati reali (es. sensori in siti archeologici come Ostia Antica).

2. Metodologia

Gli autori propongono un nuovo modello matematico formulato come un sistema non lineare di reazione-trasporto in mezzi porosi, governato dal flusso di Darcy.

A. Il Modello Matematico
Il sistema è descritto da equazioni differenziali alle derivate parziali (PDE) che modellano:

Trasporto di umidità ( $\theta$ ): Governato da un'equazione di Richards non lineare, che descrive la diffusione dell'acqua nei pori.
Trasporto di acido carbonico ( $c_a$ ): Include diffusione (legge di Fick), trasporto convettivo e reazione chimica.
Evoluzione della porosità ( $n$ ): La reazione chimica $CaCO_3 + H_2CO_3 \rightarrow Ca^{2+} + 2HCO_3^-$ consuma la matrice solida, aumentando la porosità fino a un valore critico ( $n_{max}$ ), oltre il quale si verifica la perdita di materiale (erosione).

Il sistema è accoppiato con condizioni al contorno di tipo Neumann/Robin che descrivono l'interazione tra l'ambiente esterno (umidità $E$ , concentrazione $C$ ) e la superficie della pietra.

B. Funzioni di Assorbimento
Per descrivere la relazione tra pressione capillare e contenuto di umidità (funzione $B$ ), vengono proposte due formulazioni:

Simmetrica ( $B'$ ): Più semplice, basata su 3 parametri.
Asimmetrica ( $B'_{kP}$ ): Più complessa, che separa permeabilità e pressione capillare, permettendo una descrizione più ricca delle proprietà del materiale (adatta a dati granulometrici e test di permeabilità).

C. Metodo Numerico
Per gestire il movimento del fronte di erosione (un confine mobile che cambia nel tempo), gli autori utilizzano un approccio ibrido:

Differenze Finite: Per la discretizzazione interna del dominio.
Metodo Level-Set: Per tracciare il confine mobile. La funzione level-set $\phi(x) = n(x) - n_{max}$ definisce la superficie della pietra.
Ghost Points: Per applicare le condizioni al contorno su griglie cartesiane dove il confine non è allineato, utilizzando interpolazione di Lagrange di secondo ordine.
Schemi Temporali: Metodo di Eulero implicito per garantire stabilità con passi temporali ampi.

3. Risultati Principali

I test numerici sono stati condotti in 1D e 2D con parametri calibrati su dati letterari e sperimentali (es. tassi di recessione di 2 mm/100 anni per siti suburbani).

Confronto Modelli: Le simulazioni con le funzioni di assorbimento simmetriche e asimmetriche producono risultati molto simili per i profili di umidità e acido, confermando la validità di entrambe, sebbene quella asimmetrica offra maggiore flessibilità parametrica.
Dinamica del Fronte di Erosione:
- Il modello riproduce realisticamente l'avanzamento del fronte di erosione.
- In condizioni standard (1 anno di simulazione), il tasso di erosione calcolato è di circa $2.36 \times 10^{-6}$ cm/h, coerente con i dati letterari.
- In condizioni "catastrofiche" (pietra sommersa, saturazione completa), il tasso di erosione aumenta drasticamente a $8.25 \times 10^{-4}$ cm/h, dimostrando la sensibilità del modello all'umidità.
Convergenza: Lo schema numerico mostra un ordine di convergenza globale di 1 (coerente con la discretizzazione di primo ordine in spazio e tempo).
Simulazioni 2D: Il modello gestisce efficacemente geometrie complesse e confini mobili in due dimensioni, mostrando l'evoluzione della porosità e la perdita di materiale nel tempo.

4. Contributi Chiave

Nuovo Framework Matematico: Integrazione di chimica (dissoluzione carbonatica), fisica dei pori (flusso di Darcy/Reichards) e cinetica di erosione in un unico sistema accoppiato.
Gestione del Confine Mobile: Applicazione efficace del metodo Level-Set combinato con differenze finite e ghost points per simulare l'erosione su geometrie arbitrarie senza bisogno di rimodellare la griglia.
Flessibilità Operativa: Il modello è progettato per incorporare dati ambientali in tempo reale (concentrazione $CO_2$ , pioggia, umidità, temperatura) e caratteristiche del materiale, permettendo scenari predittivi personalizzati.
Validazione: I risultati numerici sono in accordo con i tassi di erosione riportati in letteratura per aree urbane e suburbane.

5. Significato e Prospettive Future

Questo lavoro rappresenta un passo fondamentale verso la conservazione predittiva del patrimonio culturale.

Impatto Pratico: Il modello può essere utilizzato per sviluppare algoritmi di previsione che correlano i livelli di inquinamento locale con il danno sui materiali, supportando strategie di conservazione specifiche per sito.
Limitazioni Attuali: I risultati sono preliminari a causa della necessità di calibrare meglio alcuni parametri ( $K_w, K_a, K_n, n_{max}$ ) con dati sperimentali specifici su diversi tipi di pietra.
Sviluppi Futuri:
- Estensione a 3D.
- Inclusione di effetti gravitazionali (risalita capillare).
- Modellazione della formazione di croste (interazione tra erosione e deposizione).
- Simulazione di cicli gelo-disgelo.
- Aumento dell'ordine di accuratezza dello schema numerico.

In sintesi, il paper offre uno strumento matematico robusto e flessibile per comprendere e prevedere il degrado chimico delle pietre storiche, spostando l'attenzione da stime statistiche medie a modelli dinamici basati sui processi fisici e chimici reali.