Having in mind the modelling of marble degradation under chemical pollutants, e.g. the sulfation process, we consider governing nonlinear diffusion equations and their numerical approximation. The space domain of a computation is the pristine marble object. In order to accurately discretize it while maintaining the simplicity of finite difference discretizations, the domain is described using a level-set technique. A uniform Cartesian grid is laid over a box containing the domain, but the solution is defined and updated only in the grid nodes that lie inside the domain, the level-set being employed to select them and to impose accurately the boundary conditions. We use a Crank-Nicolson scheme in time, while for the space variables the discretization is performed by a standard Finite-Difference scheme for grid points inside the domain and by a ghost-cell technique on the ghost points (by using boundary conditions). The solution of the large nonlinear system is obtained by a Newton-Raphson procedure and a tailored multigrid technique is developed for the inner linear solvers. The numerical results, which are very satisfactory in terms of reconstruction quality and of computational efficiency, are presented and discussed at the end of the paper.

A level-set multigrid technique for nonlinear diffusion in the numerical simulation of marble degradation under chemical pollutants

Semplice M.;Serra Capizzano S.
2020-01-01

Abstract

Having in mind the modelling of marble degradation under chemical pollutants, e.g. the sulfation process, we consider governing nonlinear diffusion equations and their numerical approximation. The space domain of a computation is the pristine marble object. In order to accurately discretize it while maintaining the simplicity of finite difference discretizations, the domain is described using a level-set technique. A uniform Cartesian grid is laid over a box containing the domain, but the solution is defined and updated only in the grid nodes that lie inside the domain, the level-set being employed to select them and to impose accurately the boundary conditions. We use a Crank-Nicolson scheme in time, while for the space variables the discretization is performed by a standard Finite-Difference scheme for grid points inside the domain and by a ghost-cell technique on the ghost points (by using boundary conditions). The solution of the large nonlinear system is obtained by a Newton-Raphson procedure and a tailored multigrid technique is developed for the inner linear solvers. The numerical results, which are very satisfactory in terms of reconstruction quality and of computational efficiency, are presented and discussed at the end of the paper.
2020
https://arxiv.org/abs/1902.07029
Coco, A.; Semplice, M.; Serra Capizzano, S.
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11383/2100138
 Attenzione

L'Ateneo sottopone a validazione solo i file PDF allegati

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 5
  • ???jsp.display-item.citation.isi??? 3
social impact