Fast Computation of Blood Flows in Arterial Models

The computation of the blood flow in arterial models requires the solution of the incompressible Navier-Stokes equations [1,2]. The geometry of these models is characterized by complex branching tubular domains. The main challenge with this kind of geometries is that the convergence rate of the pressure Poisson solver is dominated by the graph depth of the computational grid. The purpose of this paper is to present a novel technique (deflated conjugate gradients algorithm) for accelerating the solution of the pressure Poisson equation that is especially suited for elongated domains.