Gas phase appearance and disappearance as a problem with complementarity constraints

The modeling of migration of hydrogen produced by the corrosion of the nuclear waste packages in an underground storage including the dissolution of hydrogen involves a set of nonlinear partial differential equations with nonlinear complementarity constraints. This article shows how to apply a modern and efficient solution strategy, the Newton-min method, to this geoscience problem and investigates its applicability and efficiency. In particular, numerical experiments show that the Newton-min method is quadratically convergent for this problem.

[1]  Peter Knabner,et al.  Fully coupled generalized hybrid-mixed finite element approximation of two-phase two-component flow in porous media. Part I: formulation and properties of the mathematical model , 2013, Computational Geosciences.

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

[3]  Patrick T. Harker,et al.  Finite-dimensional variational inequality and nonlinear complementarity problems: A survey of theory, algorithms and applications , 1990, Math. Program..

[4]  Ibtihel Ben Gharbia,et al.  Nonconvergence of the plain Newton-min algorithm for linear complementarity problems with a P-matrix , 2012, Math. Program..

[5]  F. Clarke.  Optimization And Nonsmooth Analysis , 1983 .

[6]  Stanislav N. Antontsev,et al.  Homogenization of Immiscible Compressible Two-Phase Flow in Porous Media: Application to Gas Migration in a Nuclear Waste Repository , 2010, Multiscale Model. Simul..

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

[8]  Jérôme Jaffré,et al.  Henry’ Law and Gas Phase Disappearance , 2009, 0904.1195.

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

[10]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[11]  Farid Smaï,et al.  Two-phase, partially miscible flow and transport modeling in porous media; application to gas migration in a nuclear waste repository , 2008, 0802.4389.

[12]  Rainer Helmig,et al.  A new approach for phase transitions in miscible multi-phase flow in porous media , 2011 .

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

[14]  Ibtihel Ben Gharbia,et al.  An Algorithmic Characterization of P-Matricity , 2013, SIAM J. Matrix Anal. Appl..

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

[16]  Peter Knabner,et al.  Fully coupled generalised hybrid-mixed finite element approximation of two-phase two-component flow in porous media. Part II: numerical scheme and numerical results , 2012, Computational Geosciences.

[17]  Anahita Abadpour,et al.  Asymptotic Decomposed Model of Two-Phase Compositional Flow in Porous Media: Analytical Front Tracking Method for Riemann Problem , 2010 .

[18]  Barbara Wohlmuth,et al.  Semismooth Newton methods for variational problems with inequality constraints , 2010 .

[19]  G. Chavent.  Mathematical models and finite elements for reservoir simulation , 1986 .

[20]  S. Kräutle,et al.  The semismooth Newton method for multicomponent reactive transport with minerals , 2011 .