Reduced-order modeling of blood flow for noninvasive functional evaluation of coronary artery disease

AbstractWe present a novel computational approach, based on a parametrized reduced-order model, for accelerating the calculation of pressure drop along blood vessels. Vessel lumina are defined by a geometric parametrization using the discrete empirical interpolation method on control points located on the surface of the vessel. Hemodynamics are then computed using a reduced-order representation of the parametrized three-dimensional unsteady Navier–Stokes and continuity equations. The reduced-order model is based on an offline–online splitting of the solution process, and on the projection of a finite volume full-order model on a low-dimensionality subspace generated by proper orthogonal decomposition of pressure and velocity fields. The algebraic operators of the hemodynamic equations are assembled efficiently during the online phase using the discrete empirical interpolation method. Our results show that with this approach calculations can be sped up by a factor of about 25 compared to the conventional full-order model, while maintaining prediction errors within the uncertainty limits of invasive clinical measurement of pressure drop. This is of importance for a clinically viable implementation of noninvasive, medical imaging-based computation of fractional flow reserve.

[1]  Michail I. Papafaklis,et al.  Endothelial Shear Stress and Coronary Plaque Characteristics in Humans: Combined Frequency-Domain Optical Coherence Tomography and Computational Fluid Dynamics Study , 2014, Circulation. Cardiovascular imaging.

[2]  Matthias Gutberlet,et al.  1-Year Outcomes of FFRCT-Guided Care in Patients With Suspected Coronary Disease: The PLATFORM Study. , 2016, Journal of the American College of Cardiology.

[3]  A. Quarteroni,et al.  Numerical modeling of hemodynamics scenarios of patient-specific coronary artery bypass grafts , 2017, Biomechanics and Modeling in Mechanobiology.

[4]  Clarence W. Rowley,et al.  Model Reduction for fluids, Using Balanced Proper Orthogonal Decomposition , 2005, Int. J. Bifurc. Chaos.

[5]  K. Aida-zade,et al.  Solving systems of differential equations of block structure with nonseparated boundary conditions , 2015 .

[6]  N Westerhof,et al.  Evaluation of methods for estimation of total arterial compliance. , 1995, The American journal of physiology.

[7]  Manesh R. Patel,et al.  Quality-of-Life and Economic Outcomes of Assessing Fractional Flow Reserve With Computed Tomography Angiography: PLATFORM. , 2015, Journal of the American College of Cardiology.

[8]  R. Banerjee,et al.  In Vitro Quantification of Guidewire Flow-Obstruction Effect in Model Coronary Stenoses for Interventional Diagnostic Procedure , 2007 .

[9]  Liang Zhong,et al.  Simplified Models of Non-Invasive Fractional Flow Reserve Based on CT Images , 2016, PloS one.

[10]  A. Quarteroni,et al.  Numerical solution of parametrized Navier–Stokes equations by reduced basis methods , 2007 .

[11]  Vartan Kurtcuoglu,et al.  Patient-specific three-dimensional simulation of LDL accumulation in a human left coronary artery in its healthy and atherosclerotic states. , 2009, American journal of physiology. Heart and circulatory physiology.

[12]  A. Quarteroni,et al.  Reduced Basis Methods for Partial Differential Equations: An Introduction , 2015 .

[13]  Alfio Quarteroni,et al.  Comparisons between reduced order models and full 3D models for fluid-structure interaction problems in haemodynamics , 2014, J. Comput. Appl. Math..

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

[15]  D. Comaniciu,et al.  A machine-learning approach for computation of fractional flow reserve from coronary computed tomography. , 2016, Journal of applied physiology.

[16]  William Wijns,et al.  Percutaneous coronary intervention of functionally nonsignificant stenosis: 5-year follow-up of the DEFER Study. , 2007, Journal of the American College of Cardiology.

[17]  Gianluigi Rozza,et al.  Advances in Reduced order modelling for CFD: vortex shedding around a circular cylinder using a POD-Galerkin method , 2017 .

[18]  M. Schaap,et al.  3D reconstruction techniques of human coronary bifurcations for shear stress computations. , 2014, Journal of biomechanics.

[19]  Charles-Henri Bruneau,et al.  Enablers for robust POD models , 2009, J. Comput. Phys..

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

[21]  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.

[22]  G W Hamilton,et al.  Physiologic basis for assessing critical coronary stenosis. Instantaneous flow response and regional distribution during coronary hyperemia as measures of coronary flow reserve. , 1974, The American journal of cardiology.

[23]  A. Khera,et al.  Forecasting the Future of Cardiovascular Disease in the United States: A Policy Statement From the American Heart Association , 2011, Circulation.

[24]  Hrvoje Jasak,et al.  A tensorial approach to computational continuum mechanics using object-oriented techniques , 1998 .

[25]  Pim A. L. Tonino,et al.  Study : Fractional Flow Reserve Versus Angiography in Multivessel Angiographic Versus Functional Severity of Coronary Artery Stenoses in the , 2010 .

[26]  A. Kono,et al.  Fractional flow reserve computed from noninvasive CT angiography data: diagnostic performance of an on-site clinician-operated computational fluid dynamics algorithm. , 2015, Radiology.

[27]  Charbel Farhat,et al.  Nonlinear model order reduction based on local reduced‐order bases , 2012 .

[28]  M. Darwish,et al.  The Finite Volume Method in Computational Fluid Dynamics: An Advanced Introduction with OpenFOAM® and Matlab , 2015 .

[29]  D. Ricci,et al.  Redefining the normal angiogram using population-derived ranges for coronary size and shape: validation using intravascular ultrasound and applications in diverse patient cohorts , 2007, The International Journal of Cardiovascular Imaging.

[30]  A. Quarteroni,et al.  A reduced computational and geometrical framework for inverse problems in hemodynamics , 2013, International journal for numerical methods in biomedical engineering.

[31]  Martin J. Leahy,et al.  Measurement of the blood flow rate and velocity in coronary artery stenosis using intracoronary frequency domain optical coherence tomography: Validation against fractional flow reserve , 2014, International journal of cardiology. Heart & vasculature.

[32]  Charles A. Taylor,et al.  Computational fluid dynamics applied to cardiac computed tomography for noninvasive quantification of fractional flow reserve: scientific basis. , 2013, Journal of the American College of Cardiology.

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

[34]  P. Stella,et al.  Deferral vs. performance of percutaneous coronary intervention of functionally non-significant coronary stenosis: 15-year follow-up of the DEFER trial. , 2015, European heart journal.

[35]  G. Rozza,et al.  Finite volume POD-Galerkin stabilised reduced order methods for the parametrised incompressible Navier–Stokes equations , 2017, Computers & Fluids.

[36]  Hrvoje Jasak,et al.  Error analysis and estimation for the finite volume method with applications to fluid flows , 1996 .

[37]  Vartan Kurtcuoglu,et al.  Choosing the optimal wall shear parameter for the prediction of plaque location-A patient-specific computational study in human right coronary arteries. , 2010, Atherosclerosis.

[38]  Pablo Lamata,et al.  Catheter-induced Errors in Pressure Measurements in Vessels: an In-vitro and Numerical Study , 2014, IEEE Transactions on Biomedical Engineering.

[39]  David Amsallem,et al.  Efficient model reduction of parametrized systems by matrix discrete empirical interpolation , 2015, J. Comput. Phys..

[40]  Vartan Kurtcuoglu,et al.  Patient-Specific Surgical Planning, Where Do We Stand? The Example of the Fontan Procedure , 2015, Annals of Biomedical Engineering.

[41]  Rafael Palacios,et al.  On-Demand Aerodynamics in Integrally Actuated Membranes with Feedback Control , 2017 .

[42]  Didier Sornette,et al.  Early Warning Signals of Financial Crises with Multi-Scale Quantile Regressions of Log-Periodic Power Law Singularities , 2015, PloS one.

[43]  Alison L. Marsden,et al.  Patient-Specific Multiscale Modeling of Blood Flow for Coronary Artery Bypass Graft Surgery , 2012, Annals of Biomedical Engineering.

[44]  Alfio Quarteroni,et al.  Numerical approximation of parametrized problems in cardiac electrophysiology by a local reduced basis method , 2018, Computer Methods in Applied Mechanics and Engineering.

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

[46]  Gianluigi Rozza,et al.  Fast simulations of patient-specific haemodynamics of coronary artery bypass grafts based on a POD-Galerkin method and a vascular shape parametrization , 2016, J. Comput. Phys..

[47]  N. Nguyen,et al.  A general multipurpose interpolation procedure: the magic points , 2008 .