Journal
SIAM JOURNAL ON SCIENTIFIC COMPUTING
Volume 35, Issue 4, Pages A1761-A1791Publisher
SIAM PUBLICATIONS
DOI: 10.1137/120896918
Keywords
nonlinear diffusion partial differential equation; nonlinear algebraic system; a posteriori error estimate; adaptive linearization; adaptive algebraic solution; adaptive mesh refinement; stopping criterion
Categories
Funding
- Groupement MoMaS (PACEN/CNRS)
- Groupement MoMaS (AN-DRA)
- Groupement MoMaS (BRGM)
- Groupement MoMaS (CEA)
- Groupement MoMaS (EdF)
- Groupement MoMaS (IRSN)
- ERT project Enhanced Oil Recovery and Geological Sequestration of CO2: Mesh Adaptivity, A Posteriori Error Control, and Other Advanced Techniques (LJLL/IFPEN)
Ask authors/readers for more resources
We consider nonlinear algebraic systems resulting from numerical discretizations of nonlinear partial differential equations of diffusion type. To solve these systems, some iterative nonlinear solver and, on each step of this solver, some iterative linear solver are used. We derive adaptive stopping criteria for both iterative solvers. Our criteria are based on an a posteriori error estimate which distinguishes the different error components, namely, the discretization error, the linearization error, and the algebraic error. We stop the iterations whenever the corresponding error no longer affects the overall error significantly. Our estimates also yield a guaranteed upper bound on the overall error at each step of the nonlinear and linear solvers. We prove the ( local) efficiency and robustness of the estimates with respect to the size of the nonlinearity owing, in particular, to the error measure involving the dual norm of the residual. Our developments hinge on equilibrated flux reconstructions and yield a general framework. We show how to apply this framework to various discretization schemes like finite elements, nonconforming finite elements, discontinuous Galerkin, finite volumes, and mixed finite elements; to different linearizations like fixed point and Newton; and to arbitrary iterative linear solvers. Numerical experiments for the p-Laplacian illustrate the tight overall error control and important computational savings achieved in our approach.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available