Having in mind the modelling of marble degradation under chemical pollutants, e.g. the sulfation process, we consider governing nonlinear diffusion equations and their numer-ical 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 dis-cretizations, 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. (C) 2020 Elsevier Inc. All rights reserved.

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

Coco, A
Primo
;
2020

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 numer-ical 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 dis-cretizations, 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. (C) 2020 Elsevier Inc. All rights reserved.
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/20.500.11769/541599
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus ND
  • ???jsp.display-item.citation.isi??? 0
social impact