Problems with localized nonhomogeneous material properties present well-known challenges for numerical simulations. In particular, such problems may feature large differences in length scales, causing difficulties with meshing and preconditioning. These difficulties are increased if the region of localized dynamics changes in time. Overlapping domain decomposition methods, which split the problem at the continuous level, show promise due to their ease of implementation and computational efficiency. Accordingly, the present work aims to further develop the mathematical theory of such methods at both the continuous and discrete levels. For the continuous formulation of the problem, we provide a full convergence analysis. For the discrete problem, we show how the described method may be interpreted as a Gauss-Seidel scheme or as a Neumann series approximation, establishing a convergence criterion in terms of the spectral radius of the system. We then provide a spectral scaling argument and provide numerical evidence for its justification.

A theoretical and numerical analysis of a Dirichlet-Neumann domain decomposition method for diffusion problems in heterogeneous media

S Bertoluzza;F Auricchio
2022

Abstract

Problems with localized nonhomogeneous material properties present well-known challenges for numerical simulations. In particular, such problems may feature large differences in length scales, causing difficulties with meshing and preconditioning. These difficulties are increased if the region of localized dynamics changes in time. Overlapping domain decomposition methods, which split the problem at the continuous level, show promise due to their ease of implementation and computational efficiency. Accordingly, the present work aims to further develop the mathematical theory of such methods at both the continuous and discrete levels. For the continuous formulation of the problem, we provide a full convergence analysis. For the discrete problem, we show how the described method may be interpreted as a Gauss-Seidel scheme or as a Neumann series approximation, establishing a convergence criterion in terms of the spectral radius of the system. We then provide a spectral scaling argument and provide numerical evidence for its justification.
2022
Istituto di Matematica Applicata e Tecnologie Informatiche - IMATI -
Finite element methods
Numerical PDEs
Multimesh methods
Thermal problems
Multiscale methods
Domain decomposition
File in questo prodotto:
File Dimensione Formato  
prod_463835-doc_181774.pdf

accesso aperto

Descrizione: A theoretical and numerical analysis of a Dirichlet-Neumann domain decomposition method for diffusion problems in heterogeneous media
Tipologia: Documento in Pre-print
Licenza: Altro tipo di licenza
Dimensione 658.54 kB
Formato Adobe PDF
658.54 kB Adobe PDF Visualizza/Apri
prod_463835-doc_181775.pdf

solo utenti autorizzati

Descrizione: A theoretical and numerical analysis of a Dirichlet-Neumann domain decomposition method for diffusion problems in heterogeneous media
Tipologia: Versione Editoriale (PDF)
Licenza: NON PUBBLICO - Accesso privato/ristretto
Dimensione 655.9 kB
Formato Adobe PDF
655.9 kB Adobe PDF   Visualizza/Apri   Richiedi una copia
TheoreticalTwoLeve_PostPrintl.pdf

Open Access dal 27/11/2023

Licenza: Creative commons
Dimensione 731.16 kB
Formato Adobe PDF
731.16 kB Adobe PDF Visualizza/Apri

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.14243/446050
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 1
  • ???jsp.display-item.citation.isi??? 1
social impact