FEM multigrid techniques for fluid-structure interaction with application to hemodynamics

We present special finite element and multigrid techniques for solving prototypical cerebral aneurysm hemodynamics problems numerically. An arbitrary Lagrangian-Eulerian (ALE) formulation is employed for this fluid-structure interaction (FSI) application. We utilize the well-known high order finite element pair Q"2P"1 for discretization in space to gain high accuracy and robustness and perform as time-stepping a fully implicit second order accurate time integrator. The resulting nonlinear discretized algebraic system is solved by an iterative Newton solver which approximates the Jacobian matrix by the divided difference approach, and the resulting linear system is solved by means of Krylov type and geometrical multigrid solvers with a Vanka-like smoother. The aim of this paper is to study the interaction of the elastic walls of an aneurysm with the geometrical shape of an implanted stent structure for prototypical 2D configurations. Preliminary results for the stent-assisted occlusion of a cerebral aneurysm and a qualitative analysis of the behavior of the elasticity of the walls vs. the geometrical details of the stent for prototypical flow situation are presented. Additionally, our approach is designed in such a way that complicated realistic constitutive relations for biomechanics applications for blood vessel simulations can be easily integrated.

[1]  J. Baranger,et al.  A formulation of Stokes's problem and the linear elasticity equations suggested by the Oldroyd model for viscoelastic flow , 1992 .

[2]  Jaroslav Hron,et al.  Numerical Study of a Modified Time-Stepping θ-Scheme for Incompressible Flow Simulations , 2006, J. Sci. Comput..

[3]  Timothy A. Davis,et al.  A combined unifrontal/multifrontal method for unsymmetric sparse matrices , 1999, TOMS.

[4]  Daniele Boffi,et al.  On the quadrilateral Q2-P1 element for the Stokes problem , 2000 .

[5]  Douglas N. Arnold,et al.  Approximation by quadrilateral finite elements , 2000, Math. Comput..

[6]  A D McCulloch,et al.  A three-dimensional finite element method for large elastic deformations of ventricular myocardium: I--Cylindrical and spherical polar coordinates. , 1996, Journal of biomechanical engineering.

[7]  Alfio Quarteroni,et al.  Modeling the cardiovascular system: a mathematical challenge , 2001 .

[8]  Rainald Löhner,et al.  Computational modelling of blood flow in side arterial branches after stenting of cerebral aneurysms , 2008 .

[9]  Alfio Quarteroni,et al.  Computational vascular fluid dynamics: problems, models and methods , 2000 .

[10]  Matthias Heil,et al.  Stokes flow in collapsible tubes: computation and experiment , 1997, Journal of Fluid Mechanics.

[11]  Miguel A. Fernández,et al.  Numerical simulation of blood flows through a porous interface , 2008 .

[12]  S. Turek,et al.  A Computational Comparison of Two FEM Solvers for Nonlinear Incompressible Flow , 2003 .

[13]  Alvaro Valencia,et al.  Blood flow dynamics and fluid–structure interaction in patient‐specific bifurcating cerebral aneurysms , 2008 .

[14]  S. Vanka.  Block-implicit multigrid solution of Navier-Stokes equations in primitive variables , 1986 .

[15]  Philip M. Gresho,et al.  On the theory of semi‐implicit projection methods for viscous incompressible flow and its implementation via a finite element method that also introduces a nearly consistent mass matrix. Part 1: Theory , 1990 .

[16]  Stefan Turek,et al.  Efficient Solvers for Incompressible Flow Problems - An Algorithmic and Computational Approach , 1999, Lecture Notes in Computational Science and Engineering.

[17]  Stefan Turek,et al.  Numerical Simulation of Laminar Incompressible Fluid-Structure Interaction for Elastic Material with Point Constraints , 2010 .

[18]  H. Damanik,et al.  A monolithic FEM approach for the log-conformation reformulation (LCR) of viscoelastic flow problems , 2010 .

[19]  A. Ouazzi.  Finite element simulation of nonlinear fluids with application to granular material and powder , 2006 .

[20]  P J Hunter,et al.  A three-dimensional finite element method for large elastic deformations of ventricular myocardium: II--Prolate spheroidal coordinates. , 1996, Journal of biomechanical engineering.

[21]  R. Rannacher,et al.  Simple nonconforming quadrilateral Stokes element , 1992 .

[22]  Rekha Ranjana Rao,et al.  A Newton-Raphson Pseudo-Solid Domain Mapping Technique for Free and Moving Boundary Problems , 1996 .

[23]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[24]  Patrick Le Tallec,et al.  Numerical analysis of a linearised fluid-structure interaction problem , 2000, Numerische Mathematik.

[25]  M. Heil.  Stokes flow in an elastic tube—a large-displacement fluid-structure interaction problem , 1998 .

[26]  S. Turek.  a Comparative Study of Time-Stepping Techniques for the Incompressible Navier-Stokes Equations: from Fully Implicit Non-Linear Schemes to Semi-Implicit Projection Methods , 1996 .

[27]  Stefan Turek,et al.  A Monolithic FEM/Multigrid Solver for an ALE Formulation of Fluid-Structure Interaction with Applications in Biomechanics , 2006 .