Newton method for reactive solute transport with equilibrium sorption in porous media

We present a mass conservative numerical scheme for reactive solute transport in porous media. The transport is modeled by a convection-diffusion-reaction equation, including equilibrium sorption. The scheme is based on the mixed finite element method (MFEM), more precisely the lowest-order Raviart-Thomas elements and one-step Euler implicit. The underlying fluid flow is described by the Richards equation, a possibly degenerate parabolic equation, which is also discretized by MFEM. This work is a continuation of Radu et al. (2008) and Radu et al. (2009) [1,2] where the algorithmic aspects of the scheme and the analysis of the discretization method are presented, respectively. Here we consider the Newton method for solving the fully discrete nonlinear systems arising on each time step after discretization. The convergence of the scheme is analyzed. In the case when the solute undergoes equilibrium sorption (of Freundlich type), the problem becomes degenerate and a regularization step is necessary. We derive sufficient conditions for the quadratic convergence of the Newton scheme.

[1]  Mario Ohlberger,et al.  Adaptive finite volume approximations for weakly coupled convection dominated parabolic systems , 2002 .

[2]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[3]  Clint Dawson,et al.  Analysis of an Upwind-Mixed Finite Element Method for Nonlinear contaminant Transport Equations , 1998 .

[4]  P. Knabner,et al.  Numerical Methods for Elliptic and Parabolic Partial Differential Equations , 2003, Texts in Applied Mathematics.

[5]  Marián Slodička,et al.  Approximation of a nonlinear degenerate parabolic equation via a linear relaxation scheme , 2005 .

[6]  H. Rentz-Reichert,et al.  UG – A flexible software toolbox for solving partial differential equations , 1997 .

[7]  Sabine Attinger,et al.  A Mixed Hybrid Finite Element Discretization Scheme for Reactive Transport in Porous Media , 2008 .

[8]  Peter Knabner,et al.  Numerical simulation of contaminant biodegradation by higher order methods and adaptive time stepping , 2004 .

[9]  Martin Vohralík,et al.  A combined finite volume–nonconforming/mixed-hybrid finite element scheme for degenerate parabolic problems , 2006, Numerische Mathematik.

[10]  Eun-Jae Park,et al.  Mixed finite element methods for nonlinear second-order elliptic problems , 1995 .

[11]  Stefano Micheletti,et al.  On Some Mixed Finite Element Methods with Numerical Integration , 2001, SIAM J. Sci. Comput..

[12]  W. Jäger,et al.  Solution of doubly nonlinear and degenerate parabolic problems by relaxation schemes , 1995 .

[13]  Robert Klöfkorn,et al.  Local adaptive methods for convection dominated problems , 2002 .

[14]  Rainer Helmig,et al.  Recent advances in finite element methods for multi-phase flow processes in porous media , 2006 .

[15]  Peter Knabner,et al.  Error estimates for a mixed finite element discretization of some degenerate parabolic equations , 2008, Numerische Mathematik.

[16]  Almerico Murli,et al.  Numerical Mathematics and Advanced Applications , 2003 .

[17]  F. Radu,et al.  Mixed finite elements for the Richards' equation: linearization procedure , 2004 .

[18]  Martin Vohralík,et al.  A combined finite volume–finite element scheme for the discretization of strongly nonlinear convection–diffusion–reaction problems on nonmatching grids , 2009 .

[19]  Peter Knabner,et al.  Order of Convergence Estimates for an Euler Implicit, Mixed Finite Element Discretization of Richards' Equation , 2004, SIAM J. Numer. Anal..

[20]  M. Wheeler,et al.  A characteristics-mixed finite element method for advection-dominated transport problems , 1995 .

[21]  Ricardo H. Nochetto,et al.  Error estimates for two-phase stefan problems in several space variables, I: Linear boundary conditions , 1985 .

[22]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[23]  Clint Dawson,et al.  Upwind‐mixed methods for transport equations , 1999 .

[24]  Ricardo H. Nochetto,et al.  Error estimates for two-phase stefan problems in several space variables, II: Non-linear flux conditions , 1985 .

[25]  Sabine Attinger,et al.  Analysis of an Euler implicit‐mixed finite element scheme for reactive solute transport in porous media , 2009 .

[26]  Marián Slodicka,et al.  A Robust and Efficient Linearization Scheme for Doubly Nonlinear and Degenerate Parabolic Problems Arising in Flow in Porous Media , 2001, SIAM J. Sci. Comput..

[27]  John W. Barrett,et al.  Finite Element Approximation of the Transport of Reactive Solutes in Porous Media. Part 1: Error Estimates for Nonequilibrium Adsorption Processes , 1997 .