🔬 matéria condensada

A Conservative Darcy-Brinkman-Stokes Formulation for PHREEQC-Coupled Reactive Transport with Cell-Based Surface Complexation

Este artigo apresenta uma formulação de Darcy-Brinkman-Stokes conservativa de massa implementada no GeoChemFoam e acoplada ao PHREEQC, permitindo simulações de transporte reativo precisas com complexação de superfície baseada em células tanto em estruturas de meios porosos resolvidas quanto não resolvidas.

Autores originais: Hannah P. Menke, Florian Doster, Julien Maes

Publicado 2026-10-15

📖 5 min de leitura🧠 Leitura aprofundada

CC BY 4.0

Autores originais: Hannah P. Menke, Florian Doster, Julien Maes

Artigo original sob licença CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta é uma explicação gerada por IA do artigo abaixo. Não foi escrita nem endossada pelos autores. Para precisão técnica, consulte o artigo original. Ler aviso legal completo

Profundamente no subsolo, a água não simplesmente flui através de túneis vazios como carros em uma rodovia. Em vez disso, ela se move através de um labirinto de rocha que é porosa, o que significa que é preenchida com minúsculos buracos e rachaduras interconectados. Em muitos desses ambientes subterrâneos, a rocha não é apenas um recipiente passivo; suas superfícies são quimicamente ativas. À medida que a água se move, minerais dissolvidos e contaminantes podem aderir a essas superfícies, um processo que retarda seu deslocamento e altera sua concentração. Essa interação é crítica para entender como poluentes se movem através das águas subterrâneas, como podemos capturar minerais valiosos como o lítio de salmouras e como resíduos radioativos podem ser armazenados com segurança por milhares de anos. O desafio para os cientistas é que a estrutura interna da rocha é incrivelmente complexa. Algumas partes possuem canais amplos e abertos por onde a água corre rapidamente, enquanto outras partes são feitas de um material fino, semelhante a uma esponja, com poros microscópicos que retêm uma vasta área de superfície interna. Os modelos computacionais tradicionais frequentemente têm dificuldade em capturar simultaneamente tanto o fluxo rápido nos grandes canais quanto a química lenta e pegajosa que ocorre dentro dos poros minúsculos e invisíveis.

Uma equipe de pesquisadores desenvolveu uma nova maneira de simular essas interações complexas, permitindo modelar tanto o fluxo de água quanto a adesão química de partículas em uma grade única e unificada. A abordagem deles trata a rocha não como uma coleção de buracos grandes e esponjas minúsculas separadas, mas como um material contínuo onde a quantidade de água e a quantidade de superfície reativa podem variar de um pequeno bloco do modelo para o próximo. Neste novo método, o computador divide a imagem da rocha em pequenas células. Para cada célula, o modelo calcula quanta água está realmente se movendo através dela e quanta área de superfície está disponível para os produtos químicos aderirem, mesmo que essa área de superfície seja pequena demais para ser vista na imagem original. Ao monitorar a massa total dos produtos químicos — tanto a parte que flutua na água quanto a parte que está grudada na rocha — os pesquisadores garantem que nada seja perdido ou criado por erro durante o cálculo. Essa conservação de massa é vital para fazer previsões precisas sobre quanto tempo leva para um contaminante viajar ou quanto de um mineral valioso pode ser recuperado.

A equipe testou seu novo método contra uma série de testes rigorosos para garantir que funcionasse corretamente. Eles começaram com cenários unidimensionais simples, onde a resposta poderia ser calculada à mão, confirmando que seu modelo computacional correspondia aos resultados esperados para a velocidade com que os produtos químicos se movem e o quanto eles desaceleram devido à adesão. Eles então passaram para imagens tridimensionais mais complexas de rochas reais e materiais sintéticos. Em um teste, simularam uma rocha com formato esférico suave para ver se seu método poderia medir a área de superfície com precisão. Eles descobriram que sua abordagem poderia calcular a área de superfície com um erro de menos de um por cento, mesmo quando a resolução da imagem era relativamente baixa. Em outro teste, compararam seu novo método com uma simulação muito mais detalhada e computacionalmente cara, que mapeava explicitamente cada minúsculo poro. Seu novo método, mais rápido, reproduziu os resultados da simulação detalhada com uma diferença de menos de um por cento no total de lítio capturado, provando que poderia lidar com a complexidade da rocha do mundo real sem precisar resolver cada detalhe microscópico.

A equipe então aplicou seu método a três cenários distintos do mundo real para ver como ele se comportava na prática. Primeiro, modelaram um sistema de rocha fraturada onde íons de chumbo se moviam através de rachaduras revestidas com óxido de ferro. A simulação mostrou como o chumbo desacelerava ao aderir ao revestimento, criando uma imagem clara de como o contaminante se espalharia pela rede de fraturas. Segundo, simularam a captura de manganês de uma solução de salmoura contendo uma grande quantidade de magnésio. Como o manganês e o magnésio competem pelos mesmos pontos de adesão na superfície da rocha, o modelo teve que calcular qual deles venceria. Os resultados mostraram como o manganês gradualmente deslocava o magnésio, acumulando-se na superfície da rocha em um padrão que dependia do fluxo da água. Finalmente, modelaram o movimento de estrôncio através de uma formação de calcário, um cenário relevante para a segurança de depósitos geológicos profundos de resíduos nucleares. Neste caso, o estrôncio teve que competir com o cálcio, que já estava presente na água. A simulação rastreou com sucesso como o estrônio empurrava o cálcio para fora das superfícies da rocha conforme avançava, criando uma onda distinta de troca química.

O principal achado deste trabalho é que é possível modelar o transporte reativo em formações rochosas complexas e de múltiplas escalas sem a necessidade de resolver cada poro microscópico individualmente. Ao tratar a área de superfície reativa como uma propriedade da própria célula, em vez de apenas o limite visível da rocha, os pesquisadores conseguem incluir o efeito de superfícies internas invisíveis que são cruciais para as reações químicas. Isso permite simulações precisas de como contaminantes e minerais valiosos se movem pelo subsolo, mesmo quando a rocha contém uma mistura de grandes espaços abertos e regiões minúsculas semelhantes a esponjas. O método é conservativo, o que significa que segue estritamente as leis do balanço de massa, garantindo que a quantidade total de material permaneça constante durante toda a simulação. Essa confiabilidade torna o método uma ferramenta poderosa para prever o comportamento de longo prazo de sistemas subterrâneos, desde a limpeza de aquíferos poluídos até o design de armazenamento mais seguro para resíduos radioativos. Os pesquisadores demonstraram que sua abordagem funciona em uma ampla gama de velocidades de fluxo e tipos de rocha, fornecendo uma estrutura robusta para compreender a química oculta das profundezas da Terra.

Resumo Técnico: Uma Formulação Conservativa de Darcy-Brinkman-Stokes para Transporte Reativo Acoplado ao PHREEQC

Definição do Problema
O transporte reativo em meios porosos de múltiplas escalas, como os carbonatos, apresenta um desafio significativo de modelagem. Esses sistemas frequentemente contêm macroporos abertos ao lado de regiões de matriz microporosa. Simulações baseadas em imagens podem resolver os macroporos, mas a área de superfície reativa interna das fases microporosas está frequentemente abaixo da resolução da imagem (escala sub-voxel). Os modelos clássicos de escala de Darcy dependem de suposições de upscaling constitutivas que não conseguem capturar a heterogeneidade espacial, enquanto modelos de escala de poro totalmente resolvidos (ex: Darcy-Brinkman-Stokes ou DBS explícito de fronteira) falham quando a área reativa dominante reside dentro de células computacionais grosseiras, em vez de na fronteira da malha. Os fluxos de trabalho existentes acoplados ao PHREEQC tipicamente impõem a química como condições de fluxo em interfaces resolvidas, deixando as áreas de superfície internas não resolvidas sem representação ou representadas por sumidouros volumétricos ad hoc que comprometem a conservação de massa e a definição física dos inventários em escala de célula.

Metodologia
Os autores apresentam uma formulação DBS conservadora de massa implementada no GeoChemFoam (um toolbox baseado em OpenFOAM) e acoplada ao PHREEQC via interface PhreeqcRM. A inovação central é uma abordagem de complexação de superfície baseada em célula, onde a superfície reativa é uma propriedade do volume da célula, e não da fronteira da malha.

Equações Governantes: O fluxo é governado pelas equações de Darcy-Brinkman-Stokes, onde um campo de porosidade centrado na célula ( $\epsilon$ ) varia suavemente do sólido para o poro aberto. O transporte é resolvido para a fração de água móvel física ( $\theta$ ), que governa o armazenamento e a difusão aquosa.
Definição de Área de Superfície: A área de superfície reativa ( $a_{surf}$ $a_{s u r f}$ ) é definida como uma densidade de volume bruto centrada na célula ( $m^2 m^{-3}_{bulk}$ $m^{2} m_{b u l k}^{- 3}$ ) composta por dois termos aditivos:
1. Termo de Interface ( $a_{interface}^{surf}$ ): Reconstruído da interface de fase poro-sólido resolvida via algoritmos de interface suave (ex: occupancyGaussian).
2. Termo Intrínseco ( $a_{intrinsic}^{surf}$ ): Uma contribuição de densidade de superfície por fase ( $s_k$ ) prescrita e sub-voxel, ponderada pela fração de fase geométrica ( $f_k$ ) dentro da célula. Isso permite a representação da área de superfície microporosa interna que não é resolvida pela segmentação da imagem.
Acoplamento Conservativo: O inventário de massa local é definido como $M = \theta C + a_{surf} \Gamma$ $M = θ C + a_{s u r f} Γ$ , onde $C$ $C$ é a concentração aquosa e $\Gamma$ $Γ$ é a concentração de superfície. O solver utiliza um esquema de divisão de Lie de primeira ordem:
1. Passo de Transporte: Avança as espécies aquosas de forma conservativa usando o transporte DBS de volume finito.
2. Passo de Química: Mapeia o estado aquoso atualizado, a fração de água móvel e o inventário de superfície baseado em célula para o PHREEQC. O PHREEQC resolve o equilíbrio termodinâmico local (complexação de superfície) e retorna as espécies aquosas e adsorvidas atualizadas, garantindo a conservação estrita de massa entre as fases aquosa e de superfície dentro de cada célula.
Implementação: A área de superfície reativa e as densidades de sítios são tratadas como campos fixos durante a simulação. Células com frações de água móvel abaixo de um limiar ( $\theta_{mask}$ ) são excluídas da química.

Contribuições Principais

Desacoplamento de Volumes de Fluxo e Química: A formulação distingue a porosidade numérica DBS (que governa a resistência ao fluxo) da fração de água móvel física (que governa o volume de solução química), embora ambas sejam definidas como iguais nos casos saturados apresentados.
Superfície Não Resolvida Baseada em Célula: O método integra com sucesso a área de superfície interna não resolvida como uma densidade volumétrica intrínseca, permitindo a modelagem de transporte reativo em regiões de matriz microporosa sem a necessidade de resolução explícita da rede de poros interna.
Balanço de Massa Conservativo: A abordagem garante que as espécies dissolvidas e ligadas à superfície entrem no mesmo balanço de massa local, evitando a redistribuição heurística de termos de fonte.

Resultados e Validação
A formulação foi verificada através de um conjunto abrangente de testes:

Benchmark de Retardo 1D: A resposta de advecção-difusão-sorção acoplada correspondeu às soluções analíticas e referências de domínio finito, demonstrando fatores de retardo corretos e convergência de primeira ordem em espaço e tempo.
Balanço de Massa em Sistema Fechado: Testes de balanço de massa confirmaram que a transferência de massa entre as fases aquosa e de superfície dentro do PHREEQC preserva o inventário total com resíduos abaixo de $10^{-9}$ .
Reconstrução de Área de Superfície: Algoritmos de interface suave foram testados contra geometrias analíticas (esferas e postes circulares). Os erros foram limitados a níveis sub-percentuais (ex: $<0,2\%$ para postes 2D, $<1\%$ para esferas 3D com resolução suficiente) em comparação com referências snappyHexMesh de ajuste de corpo.
Benchmark de Duas Escalas (Resolvido-para-DBS): Uma comparação de micromodelo de lítio entre uma simulação de escala de poro totalmente resolvida e a formulação DBS escalonada mostrou que a representação de área intrínseca baseada em célula reproduziu tanto os inventários aquosos quanto os adsorvidos com uma precisão de 0,6% em números de Péclet ( $Pe$ ) de 0,1 a 10. Embora a captação total tenha sido altamente precisa, o perfil de sorção longitudinal apresentou discrepâncias de forma de até 4,75%, atribuídas ao suavizamento das variações de concentração em escala de poro.
Casos de Aplicação: O modelo foi aplicado a três cenários complexos baseados em imagens:
- Rede de Fraturas 2D: Retardo de Pb $^{2+}$ por revestimentos de óxido de ferro hidratado (HFO), capturando a heterogeneidade espacial na ocupação de sítios fortes.
- Célula de Kelvin 3D: Troca competitiva Mn $^{2+}$ /Mg $^{2+}$ , demonstrando a capacidade do modelo de lidar com a competição de múltiplas espécies em uma rede microporosa.
- Calcário Estaillades 3D: Troca Sr $^{2+}$ /Ca $^{2+}$ simulando a migração de radionuclídeos, mostrando padrões de troca cromatográfica onde o Sr $^{2+}$ desloca o Ca $^{2+}$ tanto na fase aquosa quanto na fase adsorvida.

Significância e Alegações
O artigo afirma que este framework suporta a modelagem de complexação de superfície conservativa tanto em estruturas de poros resolvidas quanto não resolvidas. Ao tratar a superfície reativa como uma propriedade da célula em vez de uma fronteira de malha, o método supera as limitações de resolução dos modelos tradicionais de escala de poro em meios de múltiplas escalas. Isso é particularmente significativo para aplicações onde a seletividade e a acessibilidade da superfície controlam o desempenho, como a retenção de contaminantes, o transporte de radionuclídeos em repositórios geológicos e a recuperação de minerais críticos (ex: extração direta de lítio). Os autores observam que a precisão da representação da área intrínseca depende de caracterização externa (ex: adsorção BET) para restringir a densidade de superfície sub-voxel, e que a implementação atual assume equilíbrio local e áreas de superfície fixas, excluindo reações cinéticas e dissolução/precipitação mineral.

Mais como este