Three-dimensional Modelling of Reactive Solutes Transport in Porous Media

The simulation of solutes transport in porous media is of importance in many disciplines including CO2 sequestration simulations, the evaluation of underground water quality, nuclear waste simulations, and petroleum reservoir engineering. We present a numerical method for coupling 3D transport with equilibrium reactions. Our method is based on a fixed-point algorithm that enables us to couple different transport and chemistry modules. A splitting technique is used for transport that enables us to use different spatial and temporal approximations for different transport processes. The Newton's method is used to solve the chemical system. The reliability of the method is assessed by comparing the obtained results in a column with the one-dimensional problem.