Numerically accelerated pore-scale equilibrium dissolution

A Parallel Coupled Lattice Boltzmann-Volume of Fluid Framework for Modeling Porous Media Evolution

In this paper, we present a framework for the modeling and simulation of a subset of physical/chemical processes occurring on different spatial and temporal scales in porous materials. In order to improve our understanding of such processes on multiple spatio-temporal scales, small-scale simulations of transport and reaction are of vital importance. Due to the geometric complexity of the pore space and the need to consider a representative elementary volume, such simulations require substantial numerical resolutions, leading to potentially huge computation times. An efficient parallelization of such numerical methods is thus vital to obtain results in acceptable wall-clock time. The goal of this paper was to improve available approaches based on lattice Boltzmann methods (LBMs) to reliably and accurately predict the combined effects of mass transport and reaction in porous media. To this end, we relied on the factorized central moment LBM as a second-order accurate approach for modeling transport. In order to include morphological changes due to the dissolution of the solid phase, the volume of fluid method with the piece-wise linear interface construction algorithm was employed. These developments are being integrated into the LBM research code VirtualFluids. After the validation of the analytic test cases, we present an application of diffusion-controlled dissolution for a pore space obtained from computer tomography (CT) scans.

Influence of Micro-Pore Connectivity and Micro-Fractures on Calcium Leaching of Cement Pastes-A Coupled Simulation Approach

A coupled numerical approach is used to evaluate the influence of pore connectivity and microcracks on leaching kinetics in fully saturated cement paste. The unique advantage of the numerical model is the ability to construct and evaluate a material with controlled properties, which is very difficult under experimental conditions. Our analysis is based on two virtual microstructures, which are different in terms of pore connectivity but the same in terms of porosity and the amount of solid phases. Numerical fracturing was performed on these microstructures. The non-fractured and fractured microstructures were both subjected to chemical leaching. Results show that despite very different material physical properties, for example, pore connectivity and effective diffusivity, the leaching kinetics remain the same as long as the amount of soluble phases, i.e., buffering capacity, is the same. The leaching kinetics also remains the same in the presence of microcracks.