25.04.13, 16:00 – 17:30 Uhr
- Michel Kern
INRIA Paris-Rocquencourt, Maison de la Simulation
Chemical reactions between species in solution and the rock matrix lead to a problem coupling partial differential equations describing flow and transport, and algebraic or ordinary differential equations for the chemical reactions. After discretization, one obtains a system of non-linear equations whose size is the number of grid cell times the number of chemical species.The talk will discuss numerical issues in the numerical simulation of reactive transport. Starting from the usual "operator splitting" formulation, we show that the Newton-Krylov method (where the linear system at each Newton iteration is solved by an iterative method) leads to a numerically strongly coupled method, with a software implementation respecting the different nature of transport and chemical modules. The behavior of the method is illustrated on a synthetic test case inspired from CO2 storage.
Pfaffenwaldring 61, Raum U1.003 (MML), Universität Stuttgart