Fast simulations of patient-specific haemodynamics of coronary artery bypass grafts based on a POD-Galerkin method and a vascular shape parametrization

In this work a reduced-order computational framework for the study of haemodynamics in three-dimensional patient-specific configurations of coronary artery bypass grafts dealing with a wide range of scenarios is proposed. We combine several efficient algorithms to face at the same time both the geometrical complexity involved in the description of the vascular network and the huge computational cost entailed by time dependent patient-specific flow simulations. Medical imaging procedures allow to reconstruct patient-specific configurations from clinical data. A centerlines-based parametrization is proposed to efficiently handle geometrical variations. POD-Galerkin reduced-order models are employed to cut down large computational costs. This computational framework allows to characterize blood flows for different physical and geometrical variations relevant in the clinical practice, such as stenosis factors and anastomosis variations, in a rapid and reliable way. Several numerical results are discussed, highlighting the computational performance of the proposed framework, as well as its capability to carry out sensitivity analysis studies, so far out of reach. In particular, a reduced-order simulation takes only a few minutes to run, resulting in computational savings of 99 % of CPU time with respect to the full-order discretization. Moreover, the error between full-order and reduced-order solutions is also studied, and it is numerically found to be less than 1 % for reduced-order solutions obtained with just O(100) online degrees of freedom.

[1]  N. Nguyen,et al.  An ‘empirical interpolation’ method: application to efficient reduced-basis discretization of partial differential equations , 2004 .

[2]  Ignacio Llamas,et al.  Bender: a virtual ribbon for deforming 3D shapes in biomedical and styling applications , 2005, SPM '05.

[3]  A. Quarteroni,et al.  Model reduction techniques for fast blood flow simulation in parametrized geometries , 2012, International journal for numerical methods in biomedical engineering.

[4]  Gianluigi Rozza,et al.  Supremizer stabilization of POD–Galerkin approximation of parametrized steady incompressible Navier–Stokes equations , 2015 .

[5]  Milan Sonka,et al.  3D Slicer as an image computing platform for the Quantitative Imaging Network. , 2012, Magnetic resonance imaging.

[6]  Guang-Zhong Yang,et al.  Spiral phase velocity mapping of left and right coronary artery blood flow: Correction for through‐plane motion using selective fat‐only excitation , 2004, Journal of magnetic resonance imaging : JMRI.

[7]  Alfio Quarteroni,et al.  Cardiovascular mathematics : modeling and simulation of the circulatory system , 2009 .

[8]  Hui Ma,et al.  Parallel Transport Approach to Curve Framing , 1995 .

[9]  Gianluigi Rozza,et al.  On optimization, control and shape design of an arterial bypass , 2005 .

[10]  Amal Ahmed Owida,et al.  Numerical analysis of coronary artery bypass grafts: An over view , 2012, Comput. Methods Programs Biomed..

[11]  H Nygaard,et al.  The anastomosis angle does change the flow fields at vascular end-to-side anastomoses in vivo. , 1995, Journal of vascular surgery.

[12]  P. Nair,et al.  Reduced‐order modeling of parameterized PDEs using time–space‐parameter principal component analysis , 2009 .

[13]  Gianluigi Rozza,et al.  Efficient geometrical parametrisation techniques of interfaces for reduced-order modelling: application to fluid–structure interaction coupling problems , 2014 .

[14]  Foad Kabinejadian,et al.  Coronary artery bypass grafting hemodynamics and anastomosis design: a biomedical engineering review , 2013, Biomedical engineering online.

[15]  Jing Chen,et al.  Mixed Finite Element Formulation and Error Estimates Based on Proper Orthogonal Decomposition for the Nonstationary Navier-Stokes Equations , 2008, SIAM J. Numer. Anal..

[16]  J D Thomas,et al.  The effect of angle and flow rate upon hemodynamics in distal vascular graft anastomoses: a numerical model study. , 1991, Journal of biomechanical engineering.

[17]  G. Rozza,et al.  On the stability of the reduced basis method for Stokes equations in parametrized domains , 2007 .

[18]  Gianluigi Rozza,et al.  Patient-specific haemodynamics of coronary artery bypass grafts: clinical imaging and reduced order models , 2013 .

[19]  Charles Taylor,et al.  EXPERIMENTAL AND COMPUTATIONAL METHODS IN CARDIOVASCULAR FLUID MECHANICS , 2004 .

[20]  Gianluigi Rozza.  Shape design by optimal flow control and reduced basis techniques , 2005 .

[21]  A. Quarteroni,et al.  OPTIMAL CONTROL AND SHAPE OPTIMIZATION OF AORTO-CORONARIC BYPASS ANASTOMOSES , 2003 .

[22]  Andrea Manzoni,et al.  An efficient computational framework for reduced basis approximation and a posteriori error estimation of parametrized Navier–Stokes flows , 2014 .

[23]  Ramon Codina,et al.  Explicit reduced‐order models for the stabilized finite element approximation of the incompressible Navier–Stokes equations , 2013 .

[24]  Martin D. Buhmann,et al.  Radial Basis Functions: Theory and Implementations: Preface , 2003 .

[25]  H. Schima,et al.  Numerical study of wall mechanics and fluid dynamics in end-to-side anastomoses and correlation to intimal hyperplasia. , 1996, Journal of biomechanics.

[26]  José M. Vega,et al.  Reduced order models based on local POD plus Galerkin projection , 2010, J. Comput. Phys..

[27]  P. Fischer,et al.  Blood Flow in End-to-Side Anastomoses ∗ , 2008 .

[28]  L. Antiga.  Patient-specific modeling of ge-ometry and blood flow in large arteries , 2003 .

[29]  Ioannis A. Kakadiaris,et al.  Functional Morphology Analysis of the Left Anterior Descending Coronary Artery in EBCT Images , 2010, IEEE Transactions on Biomedical Engineering.

[30]  Bernard Haasdonk,et al.  Convergence Rates of the POD–Greedy Method , 2013 .

[31]  Gianluigi Rozza,et al.  Reduced basis approximation and a posteriori error estimation for the time-dependent viscous Burgers’ equation , 2009 .

[32]  A. Quarteroni,et al.  Shape optimization for viscous flows by reduced basis methods and free‐form deformation , 2012 .

[33]  Martin D. Buhmann,et al.  Radial Basis Functions , 2021, Encyclopedia of Mathematical Geosciences.

[34]  A. Marsden.  Optimization in Cardiovascular Modeling , 2014 .

[35]  David A. Steinman,et al.  Robust and objective decomposition and mapping of bifurcating vessels , 2004, IEEE Transactions on Medical Imaging.

[36]  Alfio Quarteroni,et al.  Boundary control and shape optimization for the robust design of bypass anastomoses under uncertainty , 2013 .

[37]  Max Gunzburger,et al.  POD and CVT-based reduced-order modeling of Navier-Stokes flows , 2006 .

[38]  Gianluigi Rozza,et al.  A Mathematical Approach in the Design of Arterial Bypass Using Unsteady Stokes Equations , 2006, J. Sci. Comput..

[39]  B. Haasdonk,et al.  REDUCED BASIS METHOD FOR FINITE VOLUME APPROXIMATIONS OF PARAMETRIZED LINEAR EVOLUTION EQUATIONS , 2008 .

[40]  J. Butany,et al.  Histological and morphometric analyses of early and late aortocoronary vein grafts and distal anastomoses. , 1998, The Canadian journal of cardiology.

[41]  Gianluigi Rozza,et al.  Reduced Order Methods for Modeling and Computational Reduction , 2013 .

[42]  Gianluigi Rozza,et al.  Shape Design in Aorto-Coronaric Bypass Anastomoses Using Perturbation Theory , 2006, SIAM J. Numer. Anal..

[43]  John W. Kirklin,et al.  Cardiac surgery: Morphology, diagnostic criteria, natural history, techniques, results, and indications , 1986 .

[44]  P. Nair,et al.  Nonintrusive reduced‐order modeling of parametrized time‐dependent partial differential equations , 2013 .

[45]  Gianluigi Rozza,et al.  Reduced basis approximation and a posteriori error estimation for Stokes flows in parametrized geometries: roles of the inf-sup stability constants , 2013, Numerische Mathematik.

[46]  Sabine Coquillart,et al.  Interactive Axial Deformations , 1993, Modeling in Computer Graphics.

[47]  F. Migliavacca,et al.  Computational modeling of vascular anastomoses , 2005, Biomechanics and modeling in mechanobiology.

[48]  Benjamin S. Kirk,et al.  Library for Parallel Adaptive Mesh Refinement / Coarsening Simulations , 2006 .

[49]  Myung-Soo Kim,et al.  Sweep‐based Freeform Deformations , 2006, Comput. Graph. Forum.

[50]  John W. Peterson,et al.  A high-performance parallel implementation of the certified reduced basis method , 2011 .

[51]  Miguel A. Fernández,et al.  Atlas-Based Reduced Models of Blood Flows for Fast Patient-Specific Simulations , 2010, STACOM/CESC.

[52]  Alfio Quarteroni,et al.  Isogeometric approximation of cardiac electrophysiology models on surfaces: An accuracy study with application to the human left atrium , 2017 .

[53]  S Glagov,et al.  Anastomotic intimal hyperplasia: mechanical injury or flow induced. , 1992, Journal of vascular surgery.

[54]  David A. Steinman,et al.  A Framework for Geometric Analysis of Vascular Structures: Application to Cerebral Aneurysms , 2009, IEEE Transactions on Medical Imaging.

[55]  T Shimono,et al.  Mr flow measurement in the internal mammary artery-to-coronary artery bypass graft: comparison with graft stenosis at radiographic angiography. , 2001, Radiology.

[56]  David A. Steinman,et al.  An image-based modeling framework for patient-specific computational hemodynamics , 2008, Medical & Biological Engineering & Computing.

[57]  R. Bishop.  There is More than One Way to Frame a Curve , 1975 .

[58]  Miguel A. Fernández,et al.  Group-wise construction of reduced models for understanding and characterization of pulmonary blood flows from medical images , 2014, Medical Image Anal..