Preconditioning a coupled model for reactive transport in porous media
This work addresses the computational challenge of coupling transport and chemistry in porous media simulations, but the contribution is incremental as it applies known preconditioning strategies to a specific benchmark.
The paper develops numerical methods for reactive transport in porous media that separate transport and chemistry at the software level while maintaining tight numerical coupling. The approach is applied to the 1D MoMaS benchmark, demonstrating the effectiveness of block preconditioning techniques.
We study numerical methods for solving reactive transport problems in porous media that allow a separation of transport and chemistry at the software level, while keeping a tight numerical coupling between both subsystems. After recalling how each individual subsystem can be solved, we give a formulation that eliminates the local chemical concentrations, and keeps the total concentrations as unknowns. The block structure of the model is exploited both at the nonlinear level, by eliminating some unknowns, and at the linear level by using block Gauss-Seidel or block Jacobi preconditioning. The methods are applied to the easy 1D case of the MoMaS benchmark.