A Parallel Monolithic Domain Decomposition Method for Blood Flow Simulations in 3D

Abstract

We develop a parallel scalable domain decomposition method for the simulation of blood flows in compliant arteries in 3D, by using a fully coupled system of linear elasticity equation and incompressible Navier-Stokes equations. The system is discretized with a finite element method on unstructured moving meshes and solved by a Newton-Krylov algorithm preconditioned with an overlapping additive Schwarz method. We focus on the accuracy and parallel scalability of the algorithm, and report the parallel performance and robustness of the proposed approach by some numerical experiments carried out on a supercomputer with a large number of processors and for problems with millions of unknowns. © Springer-Verlag Berlin Heidelberg 2013.