Fast numerical solutions of patient‐specific blood flows in 3D arterial systems

The study of hemodynamics in arterial models constructed from patient‐specific medical images requires the solution of the incompressible flow equations in geometries characterized by complex branching tubular structures. 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. This paper presents a deflated preconditioned conjugate gradients (DPCG) algorithm for accelerating the pressure Poisson solver. A subspace deflation technique is used to approximate the lowest eigenvalues along the tubular domains. This methodology was tested with an idealized cylindrical model and three patient‐specific models of cerebral arteries and aneurysms constructed from medical images. For these cases, the number of iterations decreased by up to a factor of 16, while the total CPU time was reduced by up to 4 times when compared with the standard PCG solver. Copyright © 2009 John Wiley & Sons, Ltd.

[1]  J. Rae,et al.  Implementation of finite element Methods for Navier-Stokes equations , 1982 .

[2]  H. V. D. Vorst,et al.  The rate of convergence of Conjugate Gradients , 1986 .

[3]  R. Nicolaides.  Deflation of conjugate gradients with applications to boundary value problems , 1987 .

[4]  Lois Mansfield,et al.  On the use of deflation to improve the convergence of conjugate gradient iteration , 1988 .

[5]  Wei Jia,et al.  Numerical simulation of incompressible flows , 1989 .

[6]  K. Burrage,et al.  Restarted GMRES preconditioned by deflation , 1996 .

[7]  Frédéric Guyomarc'h,et al.  A Deflated Version of the Conjugate Gradient Algorithm , 1999, SIAM J. Sci. Comput..

[8]  C Thompson,et al.  Applied CFD techniques: An introduction based on finite element methods , 2002 .

[9]  M M Hafez.  Numerical Simulations of Incompressible Flows , 2003 .

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

[11]  Rainald Löhner.  Projective prediction of pressure increments , 2004 .

[12]  Alejandro F. Frangi,et al.  Efficient pipeline for image-based patient-specific analysis of cerebral aneurysm hemodynamics: technique and sensitivity , 2005, IEEE Transactions on Medical Imaging.

[13]  David A. Steinman,et al.  Flow Imaging and Computing: Large Artery Hemodynamics , 2005, Annals of Biomedical Engineering.

[14]  Rainald Löhner,et al.  Improving the speed and accuracy of projection-type incompressible flow solvers , 2006 .

[15]  Max Gunzburger,et al.  Incompressible Computational Fluid Dynamics: Trends and Advances , 2009 .