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

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.