AMG preconditioning for nonlinear degenerate parabolic equations on nonuniform grids with application to monument degradation

Authors:
M. Donatelli;M. Semplice;S. Serra-Capizzano
Affiliations:
Physics and Mathematics Department, University of Insubria, Via Valleggio, 11, Como, Italy;Department of Mathematics, University of Torino, Via C. Alberto, 10, Torino, Italy;Physics and Mathematics Department, University of Insubria, Via Valleggio, 11, Como, Italy
Venue:
Applied Numerical Mathematics
Year:
2013

Citing 6
Cited 0

Multigrid

Multigrid
A Globally Convergent Newton-GMRES Subspace Method for Systems of Nonlinear Equations

SIAM Journal on Scientific Computing
The asymptotic properties of the spectrum of nonsymmetrically perturbed Jacobi matrix sequences

Journal of Approximation Theory
On the Asymptotic Spectrum of Finite Element Matrix Sequences

SIAM Journal on Numerical Analysis
Preconditioned Implicit Solvers for Nonlinear PDEs in Monument Conservation

SIAM Journal on Scientific Computing
Analysis of Multigrid Preconditioning for Implicit PDE Solvers for Degenerate Parabolic Equations

SIAM Journal on Matrix Analysis and Applications

Quantified Score

Hi-index	0.00

Visualization

Abstract

Motivated by the modelling of marble degradation by chemical pollutants, we consider the approximation by implicit finite differences schemes of nonlinear degenerate parabolic equations in which sharp boundary layers naturally occur. The latter suggests to consider various types of nonuniform griddings, when defining suitable approximation schemes. The resulting large nonlinear systems are treated by Newton methods, while the locally Toeplitz linear systems arising from the Jacobian have to be solved efficiently. To this end, we propose the use of AMG preconditioners and we study the related convergence issues, together with the associated spectral features. We present some numerical experiments supporting our theoretical results on the spectrum of the coefficient matrix of the linear systems, alongside others regarding the numerical simulations in the case of the specific model.