Preconditioning a coupled model for reactive transport in porous media

Laila Amir (1) , Michel Kern (2, 3)

1 FSTG - Faculté des Sciences et Techniques Marrakech
2 SERENA - Simulation for the Environment: Reliable and Efficient Numerical Algorithms
3 CERMICS - Centre d'Enseignement et de Recherche en Mathématiques, Informatique et Calcul Scientifique

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.

Mots clés: Reactive transport; Precondtioning; Porous media; Non-linear systems

Domaines: Analyse numérique [math.NA]; Hydrologie

Laila Amir, Michel Kern. Preconditioning a coupled model for reactive transport in porous media. International Journal of Numerical Analysis and Modeling, 2019, 16 (1), pp.18-48. ⟨hal-01327307v4⟩