English

Iterative solution of piecewise linear systems for the numerical solution of obstacle problems

Numerical Analysis 2012-06-21 v1

Abstract

We investigate the use of piecewise linear systems, whose coefficient matrix is a piecewise constant function of the solution itself. Such systems arise, for example, from the numerical solution of linear complementarity problems and in the numerical solution of free-surface problems. In particular, we here study their application to the numerical solution of both the (linear) parabolic obstacle problem and the obstacle problem. We propose a class of effective semi-iterative Newton-type methods to find the exact solution of such piecewise linear systems. We prove that the semiiterative Newton-type methods have a global monotonic convergence property, i.e., the iterates converge monotonically to the exact solution in a finite number of steps. Numerical examples are presented to demonstrate the effectiveness of the proposed methods.

Keywords

Cite

@article{arxiv.0912.3222,
  title  = {Iterative solution of piecewise linear systems for the numerical solution of obstacle problems},
  author = {Luigi Brugnano and Alessandra Sestini},
  journal= {arXiv preprint arXiv:0912.3222},
  year   = {2012}
}

Comments

22 pages, 5 figures

R2 v1 2026-06-21T14:24:46.557Z