Preconditioned implicit solvers for nonlinear PDEs in monument conservation

Article, Preprint English OPEN
Semplice, M. (2010)

Mathematical models for the description, in a quantitative way, of the damages induced on the monuments by the action of specific pollutants are often systems of nonlinear, possibly degenerate, parabolic equations. Although some the asymptotic properties of the solutions are known, for a short window of time, one needs a numerical approximation scheme in order to have a quantitative forecast at any time of interest. In this paper a fully implicit numerical method is proposed, analyzed and numerically tested for parabolic equations of porous media type and on a systems of two PDEs that models the sulfation of marble in monuments. Due to the nonlinear nature of the underlying mathematical model, the use of a fixed point scheme is required and every step implies the solution of large, locally structured, linear systems. A special effort is devoted to the spectral analysis of the relevant matrices and to the design of appropriate iterative or multi-iterative solvers, with special attention to preconditioned Krylov methods and to multigrid procedures. Numerical experiments for the validation of the analysis complement this contribution.
  • References (11)
    11 references, page 1 of 2

    [ADDN04] D. Aregba Driollet, F. Diele, and R. Natalini. A mathematical model for the SO2 aggression to calcium carbonate stones: numerical approximation and asymptotic analysis. SIAM J. Appl. Math., 64(5):1636{1667, 2004.

    [AFNT07] G. Al , V. Furuholt, R. Natalini, and I. Torcicollo. A mathematical model of sulphite chemical aggression of limestones with high permeability. I. Modeling and qualitative analysis. Transp. Porous Media, 69(1):109{122, 2007.

    A.E. Berger, H. Brezis, and J.C.W Rogers. A numerical method for solving the problem ut f (u) = 0. RAIRO numerical analysis, 13:297{312, 1979.

    [BLTR00] R. Bugini, M. Laurenzi Tabasso, and M. Realini. Rate of formation of black crusts on marble. A case study. J. Cultural Heritage, pages 111{116, 2000.

    H. Brezis and A. Pazy. Convergence and approximation of semigroups of nonlinear operators in Banach spaces. J. Functional Analysis, 9:63{74, 1972.

    [CGN+09] F. Clarelli, C. Giavarini, R. Natalini, C. Nitsch, and M.L Santarelli. Mathematical models for the consolidation processes in stones. In Proc. of \International Symposium: Stone Consolidation in Cultural Heritage - research and practice". Lisbona, May 2008., 2009. to appear.

    M.G. Crandall and T.M. Liggett. Generation of Semi-Groups of non linear transformations on general Banach spaces. Amer. J. Math., 93:265{298, 1971.

    [CNPS07] F. Cavalli, G. Naldi, G. Puppo, and M. Semplice. High-order relaxation schemes for non linear degenerate di usion problems. SIAM Journal on Numerical Analysis, 45(5):2098{2119, 2007.

    [GKPC89] K.L. Gauri, N.P. Kulshreshtha, A.R. Punuru, and A.N. Chowdhury. Rate of decay of marble in laboratory and outdoor exposure. J. Mater. Civil Eng., pages 73{85, 1989.

    U. Trottenberg, C. W. Oosterlee, and A. Schuller. Multigrid. Academic Press Inc., San Diego, CA, 2001. With contributions by A.

  • Metrics
    No metrics available
Share - Bookmark