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

Authors:
F. A. Radu;I. S. Pop
Affiliations:
UFZ-Helmholtz Center for Environmental Research, Permoserstr. 15, D-04318 Leipzig, Germany and University of Jena, Wöllnitzerstr. 7, D-07749, Jena, Germany;Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
Venue:
Journal of Computational and Applied Mathematics
Year:
2010

Citing 10
Cited 1

Mixed and hybrid finite element methods

Mixed and hybrid finite element methods
A characteristics-mixed finite element method for advection-dominated transport problems

SIAM Journal on Numerical Analysis
Mixed finite element methods for nonlinear second-order elliptic problems

SIAM Journal on Numerical Analysis
Finite Element Approximation of The Transport of Reactive Solutes in Porous Media. Part II: Error Estimates for Equilibrium Adsorption Processes

SIAM Journal on Numerical Analysis
Analysis of an Upwind-Mixed Finite Element Method for Nonlinear contaminant Transport Equations

SIAM Journal on Numerical Analysis
On Some Mixed Finite Element Methods with Numerical Integration

SIAM Journal on Scientific Computing
Mixed finite elements for the Richards' equation: linearization procedure

Journal of Computational and Applied Mathematics - Special issue: Selected papers from the 2nd international conference on advanced computational methods in engineering (ACOMEN2002) Liege University, Belgium, 27-31 May 2002
Order of Convergence Estimates for an Euler Implicit, Mixed Finite Element Discretization of Richards' Equation

SIAM Journal on Numerical Analysis
A combined finite volume–nonconforming/mixed-hybrid finite element scheme for degenerate parabolic problems

Numerische Mathematik
Error estimates for a mixed finite element discretization of some degenerate parabolic equations

Numerische Mathematik

Equivalent formulations and numerical schemes for a class of pseudo-parabolic equations

Journal of Computational and Applied Mathematics

Quantified Score

Hi-index	7.29

Visualization

Abstract

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.