Journal
COMPUTATIONAL GEOSCIENCES
Volume 14, Issue 3, Pages 421-433Publisher
SPRINGER
DOI: 10.1007/s10596-009-9173-7
Keywords
Reactive transport; Porous media; Numerical simulation; Size reduction; Complementarity problems
Ask authors/readers for more resources
In this article, an approach for the efficient numerical solution of multi-species reactive transport problems in porous media is described. The objective of this approach is to reformulate the given system of partial and ordinary differential equations (PDEs, ODEs) and algebraic equations (AEs), describing local equilibrium, in such a way that the couplings and nonlinearities are concentrated in a rather small number of equations, leading to the decoupling of some linear partial differential equations from the nonlinear system. Thus, the system is handled in the spirit of a global implicit approach (one step method) avoiding operator splitting techniques, solved by Newton's method as the basic algorithmic ingredient. The reduction of the problem size helps to limit the large computational costs of numerical simulations of such problems. If the model contains equilibrium precipitation-dissolution reactions of minerals, then these are considered as complementarity conditions and rewritten as semismooth equations, and the whole nonlinear system is solved by the semismooth Newton method.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available