Numerical Modeling of Fluid-Structure Interaction in Arteries with Anisotropic Polyconvex Hyperelastic and Anisotropic Viscoelastic Material Models at Finite Strains

Research output: Contribution to journalResearch articleContributed

Contributors

  • Daniel Balzani - , Chair of Mechanics (OTT Professorship), Dresden Center for Computational Materials Science (DCMS) (Author)
  • Simone Deparis - , Swiss Federal Institute of Technology Lausanne (EPFL) (Author)
  • Simon Fausten - , University of Duisburg-Essen (Author)
  • Davide Forti - , Swiss Federal Institute of Technology Lausanne (EPFL) (Author)
  • Alexander Heinlein - , University of Cologne (Author)
  • Axel Klawonn - , University of Cologne (Author)
  • Alfio Quarteroni - , Swiss Federal Institute of Technology Lausanne (EPFL) (Author)
  • Oliver Rheinbach - , Freiberg University of Mining and Technology (Author)
  • Jörg Schröder - , University of Duisburg-Essen (Author)

Abstract

The accurate prediction of transmural stresses in arterial walls requires on the one hand robust and efficient numerical schemes for the solution of boundary value problems including fluid–structure interactions and on the other hand the use of a material model for the vessel wall that is able to capture the relevant features of the material behavior. One of the main contributions of this paper is the application of a highly nonlinear, polyconvex anisotropic structural model for the solid in the context of fluid–structure interaction, together with a suitable discretization. Additionally, the influence of viscoelasticity is investigated. The fluid–structure interaction problem is solved using a monolithic approach; that is, the nonlinear system is solved (after time and space discretizations) as a whole without splitting among its components. The linearized block systems are solved iteratively using parallel domain decomposition preconditioners. A simple – but nonsymmetric – curved geometry is proposed that is demonstrated to be suitable as a benchmark testbed for fluid–structure interaction simulations in biomechanics where nonlinear structural models are used. Based on the curved benchmark geometry, the influence of different material models, spatial discretizations, and meshes of varying refinement is investigated. It turns out that often-used standard displacement elements with linear shape functions are not sufficient to provide good approximations of the arterial wall stresses, whereas for standard displacement elements or F-bar formulations with quadratic shape functions, suitable results are obtained. For the time discretization, a second-order backward differentiation formula scheme is used. It is shown that the curved geometry enables the analysis of non-rotationally symmetric distributions of the mechanical fields. For instance, the maximal shear stresses in the fluid–structure interface are found to be higher in the inner curve that corresponds to clinical observations indicating a high plaque nucleation probability at such locations.

Details

Original languageEnglish
Article numbere02756
JournalInternational Journal for Numerical Methods in Biomedical Engineering
Volume32
Issue number10
Publication statusPublished - 2015
Peer-reviewedNo

External IDs

Scopus 84989223260

Keywords

Keywords

  • monolithic fluid–structure interaction, polyconvex hyperelasticity, domain decomposition, anisotropic, almost incompressible, parallel

Library keywords