Comparison of linear solvers for equilibrium geochemistry computations

Equilibrium chemistry computations and reactive transport modelling require the intensive use of a linear solver under very specific conditions. The systems to be solved are small or very small (4 × 4 to 20 × 20, occasionally larger) and are very ill-conditioned (condition number up to 10¹⁰⁰). These specific conditions have never been investigated in terms of the robustness, accuracy, and efficiency of the linear solver. In this work, we present the specificity of the linear system to be solved. Several direct and iterative solvers are compared using a panel of chemical systems, including or excluding the formation of mineral species. We show that direct and iterative solvers can be used for these problems and propose computational keys to improve the chemical solvers.