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

Resumo Técnico: Modelo de Erosão de Pedras Carbonatadas por CO2

1. O Problema
O artigo aborda o dano químico causado por poluentes ambientais em artefatos de pedra, especificamente em edifícios antigos e esculturas de pedra carbonatada (como mármore e calcário). O foco principal é a erosão resultante da penetração de ácido carbônico ($H_2CO_3$), formado pela reação entre dióxido de carbono ($CO_2$) atmosférico e água nos poros da pedra.

• Contexto: Aumento da concentração de $CO_2$ e poluentes (como $SO_x$ e $NO_x$) devido a atividades humanas e mudanças climáticas tem acelerado a degradação de patrimônio cultural nas últimas décadas.
• Limitações de Modelos Anteriores: Fórmulas estatísticas existentes (como o índice de erosão de Lipfert) fornecem estimativas médias e não consideram a especificidade dos materiais, a evolução temporal da degradação, fenômenos microclimáticos locais ou o efeito específico do $CO_2$ como poluente predominante.

2. Metodologia
Os autores propõem um novo modelo matemático e um algoritmo numérico para simular o processo de erosão.

• Modelo Matemático:

  • O sistema é formulado como um conjunto de equações de reação-transporte não lineares em meios porosos, governado pelo fluxo de Darcy.
  • Equações Principais: O modelo descreve a evolução da umidade ($\theta$), a concentração de ácido carbônico dissolvido ($c_a$) e a porosidade ($n$) da pedra.
  • Mecanismos Físicos:
    1. Dissolução de $CO_2$ gasoso na água dos poros, gerando ácido carbônico.
    2. Transporte de água e ácido através da matriz porosa (difusão e advecção).
    3. Reação química de dissolução do carbonato de cálcio ($CaCO_3$) pelo ácido, levando ao aumento da porosidade e, eventualmente, à perda de material.
  • Funções de Absorção: O modelo incorpora duas formulações para a função de absorção $B(s)$ (relacionando pressão capilar e conteúdo de umidade): uma simétrica e uma assimétrica, ambas derivadas da Lei de Darcy para meios não saturados.

• Abordagem Numérica:

  • Discretização: Utiliza o método de Diferenças Finitas para pontos internos do domínio computacional.
  • Gestão de Fronteiras Móveis: Para rastrear a interface de erosão (a fronteira que se move à medida que a pedra se dissolve), emprega-se o Método Level-Set combinado com o método de "pontos fantasmas" (ghost-point method).
  • Esquema Temporal: Utiliza o método de Euler Implícito para permitir passos de tempo maiores e evitar restrições severas de estabilidade.
  • Geometria: O domínio físico é embutido em uma grade cartesiana maior, permitindo a simulação de geometrias complexas e irregulares típicas de artefatos reais.

3. Contribuições Principais

• Novo Modelo Matemático: Desenvolvimento de um sistema acoplado que integra a dinâmica de fluidos em meios porosos com a cinética química de dissolução, especificamente focado na erosão induzida por $CO_2$.
• Algoritmo Robusto: Criação de um esquema numérico capaz de lidar com fronteiras móveis complexas em 1D e 2D, utilizando Level-Set para capturar a evolução da superfície da pedra com precisão.
• Flexibilidade de Parâmetros: O modelo permite a calibração com dados reais (como taxas de recessão de inscrições em lápides) e pode incorporar variáveis ambientais dinâmicas (temperatura, umidade, concentração de $CO_2$).
• Validação Comparativa: Demonstração de que diferentes formulações de funções de absorção (simétrica vs. assimétrica) produzem resultados consistentes em cenários padrão, mas oferecem flexibilidade para cenários mais complexos.

4. Resultados

• Simulações 1D e 2D: Os testes numéricos mostraram que o modelo reproduz realisticamente a evolução da frente de erosão ao longo do tempo (semanas a anos).
  • A porosidade aumenta progressivamente até atingir um limiar crítico ($n_{max}$), momento em que ocorre a perda de material.
  • O perfil de penetração de umidade e ácido segue uma progressão linear qualitativa, consistente com a literatura.
• Taxas de Erosão:
  • Em condições ambientais padrão (simulando um ano), a taxa de recuo da superfície foi de aproximadamente $2,36 \times 10^{-6}$ cm/h, alinhando-se com estimativas de áreas suburbanas ($\sim 1,7$ mm/100 anos).
  • Em condições "catastróficas" (pedra submersa/saturada), a taxa de erosão aumentou significativamente para $8,25 \times 10^{-4}$ cm/h, demonstrando a sensibilidade do modelo à umidade.
• Ordem de Convergência: O esquema numérico apresentou uma ordem de convergência global de 1 (primeira ordem no tempo e no espaço), conforme esperado para a discretização utilizada.
• Desempenho Computacional: Simulações de 1 ano em 1D levaram cerca de 30 minutos, enquanto em 2D levaram cerca de 60 horas (usando 40 núcleos), indicando viabilidade para estudos de caso.

5. Significado e Perspectivas Futuras

• Aplicação Prática: O modelo oferece uma ferramenta para prever a degradação de patrimônio cultural baseada em dados ambientais específicos do local (ex: Parque Arqueológico de Ostia Antica), permitindo estratégias de conservação personalizadas.
• Limitações Atuais: O trabalho é considerado preliminar devido à necessidade de calibração mais precisa dos parâmetros de permeabilidade e consumo ($K_w, K_a, K_n, n_{max}$) com dados experimentais mais ricos.
• Trabalhos Futuros:
  • Extensão para simulações 3D.
  • Inclusão de efeitos gravitacionais (capilaridade ascendente).
  • Modelagem de ciclos de congelamento e degelo.
  • Desenvolvimento de mecanismos de formação de crostas (interação entre erosão e deposição).
  • Aumento da ordem de precisão do esquema numérico.

Em suma, o artigo estabelece uma base matemática e computacional robusta para entender e prever a erosão química de pedras carbonatadas, superando as limitações de modelos estatísticos tradicionais ao incorporar a física e a química do processo de dissolução.