The semismooth Newton method for multicomponent reactive transport with minerals

Abstract In this article a method is proposed for the efficient simulation of reactive systems or reactive transport problems including (equilibrium) mineral precipitation–dissolution. The difficulty lies in the fact that for larger systems it is usually a priorily not known for which mineral the fluid is saturated and for which it is undersaturated. Currently wide-spread algorithms use some trial-and-error strategy requiring repeated computations, or an approximation of equilibrium precipitation–dissolution by a kinetic description. In this article we propose to formulate the problem as a so-called complementarity problem (CP) and to solve it with the semismooth Newton method, a solution strategy well known in the field of optimization theory. The CP formulation of the mineral reactions is then combined with a reformulation for the full multicomponent reactive transport problem which leads to a reduction of the number of unknowns. The reactive transport problem is tackled in the sense of a one step (global implicit) method.

[1]  Craig M. Bethke,et al.  Geochemical Reaction Modeling: Concepts and Applications , 1996 .

[2]  Peter Knabner,et al.  A reduction scheme for coupled multicomponent transport‐reaction problems in porous media: Generalization to problems with heterogeneous equilibrium reactions , 2007 .

[3]  V. S. Tripathi,et al.  A critical evaluation of recent developments in hydrogeochemical transport models of reactive multichemical components , 1989 .

[4]  C. Steefel,et al.  Approaches to modeling of reactive transport in porous media , 1996 .

[5]  T. J. Wolery,et al.  Calculation of chemical equilibrium between aqueous solution and minerals: the EQ3/6 software package. [In FORTRAN extended 4. 6 for CDC6600 and 7600] , 1979 .

[6]  Robert Mosé,et al.  New efficient algorithm for solving thermodynamic chemistry , 2002 .

[7]  Christian Kanzow,et al.  SOLUTION OF REACTIVE TRANSPORT PROBLEMS INCLUDING MINERAL PRECIPITATION-DISSOLUTION REACTIONS BY A SEMISMOOTH NEWTON METHOD1 , 2009 .

[8]  Peter Knabner,et al.  A parallel global-implicit 2-D solver for reactive transport problems in porous media based on a reduction scheme and its application to the MoMaS benchmark problem , 2010 .

[9]  L. Gardner,et al.  Geochemical Reaction Modeling , 2000 .

[10]  Fredrik Edvard Saaf.  A study of reactive transport phenomena in porous media , 1997 .

[11]  Jesús Carrera,et al.  A formulation for decoupling components in reactive transport problems , 2004 .

[12]  Michael Ulbrich,et al.  A mesh-independence result for semismooth Newton methods , 2004, Math. Program..

[13]  Liqun Qi,et al.  Convergence Analysis of Some Algorithms for Solving Nonsmooth Equations , 1993, Math. Oper. Res..

[14]  F. Facchinei,et al.  Finite-Dimensional Variational Inequalities and Complementarity Problems , 2003 .

[15]  Francisco Facchinei,et al.  A Theoretical and Numerical Comparison of Some Semismooth Algorithms for Complementarity Problems , 2000, Comput. Optim. Appl..

[16]  J. J. Morgan,et al.  Aquatic Chemistry: Chemical Equilibria and Rates in Natural Waters , 1970 .

[17]  Christian Kanzow,et al.  Inexact semismooth Newton methods for large-scale complementarity problems , 2004, Optim. Methods Softw..

[18]  Joachim Hoffmann,et al.  Comparison of numerical methods for simulating strongly nonlinear and heterogeneous reactive transport problems—the MoMaS benchmark case , 2010 .

[19]  Kazufumi Ito,et al.  On a semi-smooth Newton method and its globalization , 2009, Math. Program..

[20]  Larry W. Lake,et al.  Effect of partial local equilibrium on the propagation of precipitation/dissolution waves , 1993 .

[21]  Peter C. Lichtner,et al.  Continuum formulation of multicomponent-multiphase reactive transport , 1996 .

[22]  Peter Knabner,et al.  A new numerical reduction scheme for fully coupled multicomponent transport‐reaction problems in porous media , 2005 .

[23]  Peter C. Lichtner,et al.  Continuum model for simultaneous chemical reactions and mass transport in hydrothermal systems , 1985 .

[24]  William D. Burgos,et al.  A general paradigm to model reaction‐based biogeochemical processes in batch systems , 2003 .

[25]  Michel Kern,et al.  Reactive transport benchmark of MoMaS , 2010 .

[26]  Christian Kanzow,et al.  The semismooth Newton method for the solution of reactive transport problems including mineral precipitation-dissolution reactions , 2011, Comput. Optim. Appl..

[27]  Kazufumi Ito,et al.  The Primal-Dual Active Set Strategy as a Semismooth Newton Method , 2002, SIAM J. Optim..

[28]  Edward W. Bolton,et al.  Long-term flow/chemistry feedback in a porous medium with heterogenous permeability: Kinetic control of dissolution and precipitation , 1999 .