Lattice Boltzmann pore-scale model for multicomponent reactive transport in porous media

[ 1] In this work, a multicomponent lattice Boltzmann ( LB) model is presented for simulating reactive transport in porous media at the pore scale. In the model, a set of distribution functions is introduced to simulate fluid flow and solute transport. The model takes into account advection, diffusion, homogeneous reactions among multiple aqueous species, and heterogeneous reactions between the aqueous solution and minerals, as well as changes in solid and pore geometry. Homogeneous reactions are described through local equilibrium mass action relations. Mineral reactions are treated kinetically through boundary conditions at the mineral surface. The LB equation for flow recovers the correct pore-scale continuity and Navier-Stokes equations. The LB equations for solute transport are modified to recover advection-diffusion equations for total concentrations at the pore scale. The model is applied to a hypothetical three-component system with two aqueous complexes and two mineral reactions in a simple pore geometry. The effects of advection, diffusion, reaction rate constants, equilibrium constants of both homogeneous and heterogeneous reactions, and chemical compositions on mineral alteration of the porous medium and solute concentration are analyzed.