When simulating phenomena in physics, engineering, or applied sciences, often one has to deal with functional equations that do not admit an analytical solution. Describing these real situations is, however, possible, resorting to one of its numerical approximations and treating the resulting mathematical representation. This thesis is placed in this context: Indeed the purpose is that of furnishing several useful tools to deal with some computational problems, stemming from discretization techniques. In most of the cases the numerical methods we analyse are the classical Qp Lagrangian FEM and the more recent Galerkin B-spline Isogeometric Analysis (IgA) approximation and Staggered Discontinuous Galerkin (DG) methods. As our model PDE, we consider classical second-order elliptic differential equations and the Incompressible Navier-Stokes equations. In all these situations the resulting matrix sequences {An}n possess a structure, namely they belong to the class of Toeplitz matrix sequences or to the more general class of Generalized Locally Toeplitz (GLT) matrix sequences, in the most general block k-level case. Consequently, the spectral analysis of the coefficient matrices plays a crucial role for an efficient and fast resolution. Indeed the convergence properties of iterative methods proposed, like multigrid or preconditioned Krylov techniques, are strictly related to the notion of symbol of the coefficient matrix sequence. In our setting the symbol is a function which asymptotically provides a reasonable approximation of the eigenvalues [singular values] of An by its evaluations of an uniform grid on its domain. These reasons, and many others, make the research of more and more efficient eigensolvers relevant and topical. In this direction, the second goal of this thesis is to provide new tools for computing the spectrum of preconditioned banded symmetric Toeplitz matrices, Toeplitz-like matrices, n-1K[p]n , nM[p]n , n-2L[p]n , coming from the B-spline IgA approximation of –u” = u, plus its multivariate counterpart for -u = u, and block and preconditioned block banded symmetric Toeplitz matrices. For all the above cases we propose new algorithms based on the classical concept of symbol, but with an innovative view on the errors of the approximation of eigenvalues by the uniform sampling of the symbol. The algorithms devised are special interpolation-extrapolation procedures performed with a high level of accuracy and only at the cost of computing of the eigenvalues of a moderate number of small sized matrices.

Spectral analysis and fast methods for structured matrix sequences and PDE discretizations / Furci, Isabella. - (2018).

Spectral analysis and fast methods for structured matrix sequences and PDE discretizations.

Furci, Isabella
2018-01-01

Abstract

When simulating phenomena in physics, engineering, or applied sciences, often one has to deal with functional equations that do not admit an analytical solution. Describing these real situations is, however, possible, resorting to one of its numerical approximations and treating the resulting mathematical representation. This thesis is placed in this context: Indeed the purpose is that of furnishing several useful tools to deal with some computational problems, stemming from discretization techniques. In most of the cases the numerical methods we analyse are the classical Qp Lagrangian FEM and the more recent Galerkin B-spline Isogeometric Analysis (IgA) approximation and Staggered Discontinuous Galerkin (DG) methods. As our model PDE, we consider classical second-order elliptic differential equations and the Incompressible Navier-Stokes equations. In all these situations the resulting matrix sequences {An}n possess a structure, namely they belong to the class of Toeplitz matrix sequences or to the more general class of Generalized Locally Toeplitz (GLT) matrix sequences, in the most general block k-level case. Consequently, the spectral analysis of the coefficient matrices plays a crucial role for an efficient and fast resolution. Indeed the convergence properties of iterative methods proposed, like multigrid or preconditioned Krylov techniques, are strictly related to the notion of symbol of the coefficient matrix sequence. In our setting the symbol is a function which asymptotically provides a reasonable approximation of the eigenvalues [singular values] of An by its evaluations of an uniform grid on its domain. These reasons, and many others, make the research of more and more efficient eigensolvers relevant and topical. In this direction, the second goal of this thesis is to provide new tools for computing the spectrum of preconditioned banded symmetric Toeplitz matrices, Toeplitz-like matrices, n-1K[p]n , nM[p]n , n-2L[p]n , coming from the B-spline IgA approximation of –u” = u, plus its multivariate counterpart for -u = u, and block and preconditioned block banded symmetric Toeplitz matrices. For all the above cases we propose new algorithms based on the classical concept of symbol, but with an innovative view on the errors of the approximation of eigenvalues by the uniform sampling of the symbol. The algorithms devised are special interpolation-extrapolation procedures performed with a high level of accuracy and only at the cost of computing of the eigenvalues of a moderate number of small sized matrices.
2018
Multilevel block GLT algebra, symbol, spectral distribution, asymptotic expansion, interpolationextrapolation algorithms, multigrid methods, preconditioning Krylov methods, Staggered DG methods, IgA approximation
Spectral analysis and fast methods for structured matrix sequences and PDE discretizations / Furci, Isabella. - (2018).
File in questo prodotto:
File Dimensione Formato  
PhD_Thesis_FurciIsabella_completa.pdf

Open Access dal 20/12/2019

Descrizione: testo completo tesi
Tipologia: Tesi di dottorato
Licenza: Non specificato
Dimensione 19.87 MB
Formato Adobe PDF
19.87 MB 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/11383/2090658
 Attenzione

L'Ateneo sottopone a validazione solo i file PDF allegati

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