Mixed finite element discretization and Newton iteration for a reactive contaminant transport model with nonequilibrium sorption: convergence analysis and error estimates

F.A. Radu, I.S. Pop

Research output: Contribution to journalArticleAcademicpeer-review

12 Citations (Scopus)
1 Downloads (Pure)

Abstract

We present a numerical scheme for reactive contaminant transport with nonequilibrium sorption in porous media. The mass conservative scheme is based on Euler implicit, mixed finite elements, and Newton method. We consider the case of a Freundlich-type sorption. In this case, the sorption isotherm is not Lipschitz but just Hölder continuous. To deal with this, we perform a regularization step. The convergence of the scheme is analyzed. An explicit order of convergence depending only on the regularization parameter, the time step, and the mesh size is derived. We give also a sufficient condition for the quadratic convergence of the Newton method. Finally, relevant numerical results are presented.
Original languageEnglish
Pages (from-to)431-450
Number of pages20
JournalComputational Geosciences
Volume15
Issue number3
DOIs
Publication statusPublished - 2011

Fingerprint Dive into the research topics of 'Mixed finite element discretization and Newton iteration for a reactive contaminant transport model with nonequilibrium sorption: convergence analysis and error estimates'. Together they form a unique fingerprint.

Cite this