The studies in numerical approximation of partial differential equations are characterized by the necessity of managing complex geometries and their discretization. We focus our attention on two different fields where complex geometries are very common: the mathematical modeling of fluidstructure interaction problems and the family of virtual element methods. We first consider the finite element approximation of fluidstructure interaction problems described by a distributed Lagrange multiplier formulation. A crucial aspect of the method is the assembly of the finite element matrix describing the coupling between fluid and structure. This procedure consists in integrating over the solid domain both solid and fluid basis functions, which are defined on two different nonmatching grids. This can be done in two ways: the exact approach is based on a composite quadrature rule defined on the intersection between the two involved meshes, whereas the approximate approach is carried out roughly integrating on each solid element. We discuss and compare these two approaches both from the theoretical and computational point of view. Moreover, we present a preliminary study on a parallel solver for our fictitious domain formulation. As a first step towards the design of an effective solver, we consider block preconditioners, which are tested on two simplified academic problems and compared in terms of optimality, weak and strong scalability. Virtual element methods are known to tackle complex geometries without limitations on the degree of the polynomial that partially contributes to the approximated solution. An important aspect of this method is that we are not required to explicitly compute the basis functions of the VEM space since these are solutions of PDEs. Due to this fact, several quantities are not computed exactly, such as the bilinear form, for which the action on the nonpolynomial part is handled by a stabilization term, and the error, where the contribution of the nonpolynomial part is neglected. In certain cases, such as when anisotropic problems are considered, the method may show poor performance if endowed with standard stabilization terms, which have an isotropic structure. We propose a model order reduction technique constructed by means of the reduced basis method for efficiently solving the equation associated with each virtual basis function. The idea is to replace the stabilization term with an actual approximation of the nonpolynomial contribution. We show that this operation produces good results even if done in a very rough way.
Numerical approximation of PDEs on complex geometries: fluidstructure interaction problems and virtual element methods
CREDALI, FABIO
20231128
Abstract
The studies in numerical approximation of partial differential equations are characterized by the necessity of managing complex geometries and their discretization. We focus our attention on two different fields where complex geometries are very common: the mathematical modeling of fluidstructure interaction problems and the family of virtual element methods. We first consider the finite element approximation of fluidstructure interaction problems described by a distributed Lagrange multiplier formulation. A crucial aspect of the method is the assembly of the finite element matrix describing the coupling between fluid and structure. This procedure consists in integrating over the solid domain both solid and fluid basis functions, which are defined on two different nonmatching grids. This can be done in two ways: the exact approach is based on a composite quadrature rule defined on the intersection between the two involved meshes, whereas the approximate approach is carried out roughly integrating on each solid element. We discuss and compare these two approaches both from the theoretical and computational point of view. Moreover, we present a preliminary study on a parallel solver for our fictitious domain formulation. As a first step towards the design of an effective solver, we consider block preconditioners, which are tested on two simplified academic problems and compared in terms of optimality, weak and strong scalability. Virtual element methods are known to tackle complex geometries without limitations on the degree of the polynomial that partially contributes to the approximated solution. An important aspect of this method is that we are not required to explicitly compute the basis functions of the VEM space since these are solutions of PDEs. Due to this fact, several quantities are not computed exactly, such as the bilinear form, for which the action on the nonpolynomial part is handled by a stabilization term, and the error, where the contribution of the nonpolynomial part is neglected. In certain cases, such as when anisotropic problems are considered, the method may show poor performance if endowed with standard stabilization terms, which have an isotropic structure. We propose a model order reduction technique constructed by means of the reduced basis method for efficiently solving the equation associated with each virtual basis function. The idea is to replace the stabilization term with an actual approximation of the nonpolynomial contribution. We show that this operation produces good results even if done in a very rough way.File  Dimensione  Formato  

tesi_credali_definitiva.pdf
accesso aperto
Descrizione: Tesi_definitiva
Tipologia:
Tesi di dottorato
Dimensione
13.27 MB
Formato
Adobe PDF

13.27 MB  Adobe PDF  Visualizza/Apri 
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.