Fast numerical solutions of patient‐specific blood flows in 3D arterial systems
- 16 March 2009
- journal article
- research article
- Published by Wiley in International Journal for Numerical Methods in Biomedical Engineering
- Vol. 26 (1) , 73-85
- https://doi.org/10.1002/cnm.1235
Abstract
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.Keywords
This publication has 14 references indexed in Scilit:
- Improving the speed and accuracy of projection-type incompressible flow solversComputer Methods in Applied Mechanics and Engineering, 2006
- Flow Imaging and Computing: Large Artery HemodynamicsAnnals of Biomedical Engineering, 2005
- Efficient pipeline for image-based patient-specific analysis of cerebral aneurysm hemodynamics: technique and sensitivityIEEE Transactions on Medical Imaging, 2005
- A Deflated Version of the Conjugate Gradient AlgorithmSIAM Journal on Scientific Computing, 2000
- Restarted GMRES preconditioned by deflationJournal of Computational and Applied Mathematics, 1996
- Incompressible Computational Fluid DynamicsPublished by Cambridge University Press (CUP) ,1993
- On the use of deflation to improve the convergence of conjugate gradient iterationCommunications in Applied Numerical Methods, 1988
- Deflation of Conjugate Gradients with Applications to Boundary Value ProblemsSIAM Journal on Numerical Analysis, 1987
- The rate of convergence of Conjugate GradientsNumerische Mathematik, 1986
- Implementation of Finite Element Methods for Navier-Stokes EquationsPublished by Springer Nature ,1981