A parallel monolithic approach for fluid-structure interaction in a cerebral aneurysm

Ali Eken, Mehmet Sahin*

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

9 Citations (Scopus)

Abstract

A parallel fully-coupled numerical algorithm has been developed for the fluid-structure interaction problem in a cerebral artery with aneurysm. For the fluid part of the problem, an Arbitrary Lagrangian-Eulerian formulation based on the side-centered unstructured finite volume method is employed for the governing incompressible Navier–Stokes equations. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant-Kirchhoff material and the classical Galerkin finite element method is used to discretise the governing equations in a Lagrangian frame. The time integration method for the structure domain is based on the Newmark type generalized−α method while the second-order backward difference (BDF2) is used in the fluid domain. A special attention is given to construct an algorithm obeying the local/global discrete geometric conservation laws (DGCL) in order to conserve fluid volume at machine precision when the fluid domain is entirely enclosed by solid domain boundary. Therefore, a compatible kinematic boundary condition is applied at the interface between the solid and fluid domains. The parallel implementation of the present fully coupled unstructured fluid-structure solver is based on the PETSc library and a one-level restricted additive Schwarz preconditioner with a block-incomplete factorization within each partitioned sub-domains is utilized for the resulting fully coupled system. The proposed algorithm is initially validated for a pressure pulse propagating in a flexible tube and the mass conservation accuracy is tested for a thin elastic sphere filled with an incompressible fluid in a circular tube. Then the numerical method is applied to a complicated problem involving unsteady pulsatile blood flow in a cerebral artery with aneurysm as a realistic fluid-structure interaction problem encountered in biomechanics. Various hemodynamic quantities of interest like fluid velocities, blood pressure and wall shear stresses (WSS) are computed as well as the time dependent artery wall deformations.

Original languageEnglish
Pages (from-to)61-75
Number of pages15
JournalComputers and Fluids
Volume153
DOIs
Publication statusPublished - 10 Aug 2017

Bibliographical note

Publisher Copyright:
© 2017 Elsevier Ltd

Keywords

  • Cerebral aneurysm
  • Finite element method
  • Fluid-structure interaction
  • Monolithic method
  • Unstructured finite volume method

Fingerprint

Dive into the research topics of 'A parallel monolithic approach for fluid-structure interaction in a cerebral aneurysm'. Together they form a unique fingerprint.

Cite this