A global method for coupling transport with chemistry in heterogeneous porous media

Abstract : Modeling reactive transport in porous media, using a local chemical equilibrium assumption, leads to a system of advection-diffusion PDE's coupled with algebraic equations. When solving this coupled system, the algebraic equations have to be solved at each grid point for each chemical species and at each time step. This leads to a coupled non-linear system. In this paper a global solution approach that enables to keep the software codes for transport and chemistry distinct is proposed. The method applies the Newton-Krylov framework to the formulation for reactive transport used in operator splitting. The method is formulated in terms of total mobile and total fixed concentrations and uses the chemical solver as a black box, as it only requires that on be able to solve chemical equilibrium problems (and compute derivatives), without having to know the solution method. An additional advantage of the Newton-Krylov method is that the Jacobian is only needed as an operator in a Jacobian matrix times vector product. The proposed method is tested on the MoMaS reactive transport benchmark.
Type de document :
Article dans une revue
Liste complète des métadonnées

Contributeur : Michel Kern <>
Soumis le : vendredi 18 décembre 2009 - 17:50:08
Dernière modification le : mardi 22 mars 2016 - 01:27:03
Document(s) archivé(s) le : jeudi 17 juin 2010 - 22:02:43


Fichiers produits par l'(les) auteur(s)




Amir Laila, Michel Kern. A global method for coupling transport with chemistry in heterogeneous porous media. Computational Geosciences, Springer Verlag, 2010, 14 (3), pp.465-481. 〈http://www.springerlink.com/content/933p55085742m203/?p=9daa7a95481344afbc491cb216fa1d3e&pi=14〉. 〈10.1007/s10596-009-9162-x〉. 〈inria-00442243〉



Consultations de la notice


Téléchargements de fichiers