Parallel Algorithms for the Solution of Large-Scale Fluid-Structure Interaction Problems in Hemodynamics

This thesis addresses the development and implementation of efficient and parallel algorithms for the numerical simulation of Fluid-Structure Interaction (FSI) problems in hemodynamics. Indeed, hemodynamic conditions in large arteries are significantly affected by the interaction of the pulsatile blood flow with the arterial wall. The simulation of fluid-structure interaction problems requires the approximation of a coupled system of Partial Differential Equations (PDEs) and the set up of efficient numerical solution strategies. Blood is modeled as an incompressible Newtonian fluid whose dynamics is governed by the Navier-Stokes equations. Different constituive models are used to describe the mechanical response of the arterial wall; specifically, we rely on hyperelastic isotropic and anistotropic material laws. The finite element method is used for the space discretization of both the fluid and structure problems. In particular, for the Navier-Stokes equations we consider a semi-discrete formulation based on the Variational Multiscale (VMS) method. Among a wide range of possible solution strategies for the FSI problem, here we focus on strongly coupled monolithic approaches wherein the nonlinearities are treated in a fully implicit mode. To cope with the high computational complexity of the three dimensional FSI problem, a parallel solution framework is often mandatory. To this end, we develop a new block parallel preconditioner for the coupled linearized FSI system obtained after space and time discretization. The proposed preconditioner, named FaCSI, exploits the factorized form of the FSI Jacobian matrix, the use of static condensation to formally eliminate the interface degrees of freedom of the fluid equations, and the use of a SIMPLE preconditioner for unsteady Navier-Stokes equations. In FSI problems, the different resolution requirements in the fluid and structure physical domains, as well as the presence of complex interface geometries make the use of matching fluid and structure meshes problematic. In such situations, it is much simpler to deal with discretizations that are nonconforming at the interface, provided however that the matching conditions at the interface are properly fulfilled. In this thesis we develop a novel interpolation-based method, named INTERNODES, for numerically solving partial differential equations by Galerkin methods on computational domains that are split into two (or several) subdomains featuring nonconforming interfaces. By this we mean that either a priori independent grids and/or local polynomial degrees are used to discretize each subdomain. INTERNODES can be regarded as an alternative to the mortar element method: it combines the accuracy of the latter with the easiness of implementation in a numerical code. The aforementioned techniques have been applied for the numerical simulation of large-scale fluid-structure interaction problems in the context of biomechanics. The parallel algorithms developed showed scalability up to thousands of cores utilized on high performance computing machines

Hybrid Multiscale Methods with Applications to Semiconductors, Porous Media, and Materials Science

In this work we consider two multiscale applications with tremendous computational complexity at the lower scale. First, we examine a model for charge transport in semicon- ductor structures with heterojunction interfaces. Due to the complex physical phenomena at the interface, the model at the design scale is unable to adequately capture the behavior of the structure in the interface region. Simultaneously it is computationally intractable to simulate the full heterostructure on the scale required near the interface. Second, we con- sider the problem of the simulation of fluid flow in a dynamically evolving porous medium. The evolution of the medium strongly couples the porescale flow solutions and the macro scale model, requiring a novel approach to communicate the porescale evolution to the macroscale without resorting to the intractable simulation of the fluid flow problem di- rectly on the porescale geometry. We formulate novel methods for these two applications in the multiscale framework. For the semiconductor problem we present iterative sub- structuring domain decomposition methods that decouple the interface computation from the macroscale model. For the fluid flow problem we develop a reduced order three-scale fluid flow model based on a spatial decomposition of the porescale geometry and the offline approximation of a stochastic process describing macroscale permeability paramaterized by the volume fraction of the evolved geometry

Substructuring Preconditioners for Mortar Discretization of a Degenerate Evolution Problem

x