An optimal control approach to determine resistance‐type boundary conditions from in‐vivo data for cardiovascular simulations

The choice of appropriate boundary conditions is a fundamental step in computational fluid dynamics (CFD) simulations of the cardiovascular system. Boundary conditions, in fact, highly affect the computed pressure and flow rates, and consequently haemodynamic indicators such as wall shear stress, which are of clinical interest. Devising automated procedures for the selection of boundary conditions is vital to achieve repeatable simulations. However, the most common techniques do not automatically assimilate patient-specific data, relying instead on expensive and time-consuming manual tuning procedures. In this work, we propose a technique for the automated estimation of outlet boundary conditions based on optimal control. The values of resistive boundary conditions are set as control variables and optimized to match available patient-specific data. Experimental results on four aortic arches demonstrate that the proposed framework can assimilate 4D-Flow MRI data more accurately than two other common techniques based on Murray's law and Ohm's law. This article is protected by copyright. All rights reserved.

[1]  A. Marsden,et al.  A comparison of outlet boundary treatments for prevention of backflow divergence with relevance to blood flow simulations , 2011 .

[2]  Luca Dedè,et al.  Optimal flow control for Navier–Stokes equations: drag minimization , 2007 .

[3]  Evan M Masutani,et al.  4D Flow Vorticity Visualization Predicts Regions of Quantitative Flow Inconsistency for Optimal Blood Flow Measurement. , 2020, Radiology. Cardiothoracic imaging.

[4]  Hang Si,et al.  TetGen, a Delaunay-Based Quality Tetrahedral Mesh Generator , 2015, ACM Trans. Math. Softw..

[5]  Fuyou Liang,et al.  A multi‐scale model of the coronary circulation applied to investigate transmural myocardial flow , 2018, International journal for numerical methods in biomedical engineering.

[6]  Giancarlo Pennati,et al.  Patient‐specific parameter estimation in single‐ventricle lumped circulation models under uncertainty , 2017, International journal for numerical methods in biomedical engineering.

[7]  P. Serruys,et al.  Strain distribution over plaques in human coronary arteries relates to shear stress. , 2008, American journal of physiology. Heart and circulatory physiology.

[8]  F. N. van de Vosse,et al.  The influence of boundary conditions on wall shear stress distribution in patients specific coronary trees. , 2011, Journal of biomechanics.

[9]  Frans N van de Vosse,et al.  MRI-based quantification of outflow boundary conditions for computational fluid dynamics of stenosed human carotid arteries. , 2010, Journal of biomechanics.

[10]  Charles A. Taylor,et al.  Tuning Multidomain Hemodynamic Simulations to Match Physiological Measurements , 2010, Annals of Biomedical Engineering.

[11]  A. Veneziani,et al.  Uncertainty quantification for data assimilation in a steady incompressible Navier-Stokes problem , 2013 .

[12]  T. Pedley The Fluid Mechanics of Large Blood Vessels: Contents , 1980 .

[13]  Patrick Segers,et al.  The impact of simplified boundary conditions and aortic arch inclusion on CFD simulations in the mouse aorta: a comparison with mouse-specific reference data. , 2011, Journal of biomechanical engineering.

[14]  A. Quarteroni Numerical Models for Differential Problems , 2009 .

[15]  Gianluigi Rozza,et al.  Reduced order methods for parametric optimal flow control in coronary bypass grafts, toward patient‐specific data assimilation , 2019, International journal for numerical methods in biomedical engineering.

[16]  D. Mikhailidis,et al.  Influence of Oscillating Flow on LDL Transport and Wall Shear Stress in the Normal Aortic Arch , 2009, The open cardiovascular medicine journal.

[17]  Franck Nicoud,et al.  Data assimilation for identification of cardiovascular network characteristics , 2017, International journal for numerical methods in biomedical engineering.

[18]  Stavroula Balabani,et al.  Patient-specific haemodynamic simulations of complex aortic dissections informed by commonly available clinical datasets. , 2019, Medical engineering & physics.

[19]  Pablo J. Blanco,et al.  Boundary control in computational haemodynamics , 2017, Journal of Fluid Mechanics.

[20]  A Noordergraaf,et al.  Estimation of total systemic arterial compliance in humans. , 1990, Journal of applied physiology.

[21]  C D Murray,et al.  The Physiological Principle of Minimum Work: I. The Vascular System and the Cost of Blood Volume. , 1926, Proceedings of the National Academy of Sciences of the United States of America.

[22]  L. Chua,et al.  Computational model of blood flow in the aorto-coronary bypass graft , 2005, Biomedical engineering online.

[23]  Max Gunzburger,et al.  Perspectives in flow control and optimization , 1987 .

[24]  Jeffrey C. Lagarias,et al.  Convergence Properties of the Nelder-Mead Simplex Method in Low Dimensions , 1998, SIAM J. Optim..

[25]  Alessandro Veneziani,et al.  Patient‐specific CFD modelling in the thoracic aorta with PC‐MRI–based boundary conditions: A least‐square three‐element Windkessel approach , 2018, International journal for numerical methods in biomedical engineering.

[26]  Alison L. Marsden,et al.  SimVascular: An Open Source Pipeline for Cardiovascular Simulation , 2017, Annals of Biomedical Engineering.

[27]  J. Pepper,et al.  On the choice of outlet boundary conditions for patient-specific analysis of aortic flow using computational fluid dynamics. , 2017, Journal of biomechanics.

[28]  Anders Logg,et al.  Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book , 2012 .

[29]  M. Markl,et al.  4D flow cardiovascular magnetic resonance consensus statement , 2015, Journal of Cardiovascular Magnetic Resonance.

[30]  Alberto M. Gambaruto,et al.  Patient-specific Blood Flow Simulations: Setting Dirichlet Boundary Conditions for Minimal Error with Respect to Measured Data , 2014 .

[31]  J Tiago,et al.  A velocity tracking approach for the data assimilation problem in blood flow simulations , 2016, International journal for numerical methods in biomedical engineering.

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

[33]  Magne Nordaas,et al.  Variational data assimilation for transient blood flow simulations: Cerebral aneurysms as an illustrative example , 2016, International journal for numerical methods in biomedical engineering.

[34]  L. Antiga,et al.  Outflow conditions for image-based hemodynamic models of the carotid bifurcation: implications for indicators of abnormal flow. , 2010, Journal of biomechanical engineering.

[35]  Bradley D. Allen,et al.  4D flow imaging with MRI. , 2014, Cardiovascular diagnosis and therapy.

[36]  Martin Stoll,et al.  All-at-once solution of time-dependent Stokes control , 2013, J. Comput. Phys..

[37]  Justin S Tran,et al.  Automated Tuning for Parameter Identification and Uncertainty Quantification in Multi-scale Coronary Simulations. , 2017, Computers & fluids.

[38]  M. Salvetti,et al.  Validation of Numerical Simulations of Thoracic Aorta Hemodynamics: Comparison with In Vivo Measurements and Stochastic Sensitivity Analysis , 2018, Cardiovascular Engineering and Technology.

[39]  Ivan P. Gavrilyuk,et al.  Lagrange multiplier approach to variational problems and applications , 2010, Math. Comput..

[40]  Charles A. Taylor,et al.  Patient-specific modeling of cardiovascular mechanics. , 2009, Annual review of biomedical engineering.

[41]  Gianluigi Rozza,et al.  Reduced basis approximation of parametrized optimal flow control problems for the Stokes equations , 2015, Comput. Math. Appl..

[42]  Charles A. Taylor,et al.  Outflow boundary conditions for three-dimensional finite element modeling of blood flow and pressure in arteries , 2006 .

[43]  Anders Logg,et al.  The FEniCS Project Version 1.5 , 2015 .

[44]  G. Karniadakis,et al.  Model inversion via multi-fidelity Bayesian optimization: a new paradigm for parameter estimation in haemodynamics, and beyond , 2016, Journal of The Royal Society Interface.

[45]  Timothy J. Pedley,et al.  The fluid mechanics of large blood vessels , 1980 .

[46]  Charles A. Taylor,et al.  On Coupling a Lumped Parameter Heart Model and a Three-Dimensional Finite Element Aorta Model , 2009, Annals of Biomedical Engineering.

[47]  Patrick Amestoy,et al.  A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling , 2001, SIAM J. Matrix Anal. Appl..

[48]  T. Wonnacott,et al.  Relation between diameter and flow in major branches of the arch of the aorta. , 1992, Journal of biomechanics.

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

[50]  J-F Gerbeau,et al.  A methodological paradigm for patient‐specific multi‐scale CFD simulations: from clinical measurements to parameter estimates for individual analysis , 2014, International journal for numerical methods in biomedical engineering.

[51]  Jeff D. Eldredge,et al.  An ensemble Kalman filter approach to parameter estimation for patient-specific cardiovascular flow modeling , 2020 .

[52]  G. Karniadakis,et al.  Outflow Boundary Conditions for Arterial Networks with Multiple Outlets , 2008, Annals of Biomedical Engineering.

[53]  Fergal J Boyle,et al.  Computational fluid dynamics analysis of balloon-expandable coronary stents: influence of stent and vessel deformation. , 2014, Medical engineering & physics.

[54]  Berend E. Westerhof,et al.  The arterial Windkessel , 2009, Medical & Biological Engineering & Computing.

[55]  R. Krams,et al.  Large variations in absolute wall shear stress levels within one species and between species. , 2007, Atherosclerosis.

[56]  Tobias Schaeffter,et al.  A flexible framework for sequential estimation of model parameters in computational hemodynamics , 2020, Adv. Model. Simul. Eng. Sci..

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

[58]  Arnold P. G. Hoeks,et al.  Wall shear stress as measured in vivo: consequences for the design of the arterial system , 2008, Medical & Biological Engineering & Computing.

[59]  Marina Piccinelli,et al.  Applications of variational data assimilation in computational hemodynamics , 2012 .

[60]  Kazufumi Ito,et al.  Lagrange multiplier approach to variational problems and applications , 2008, Advances in design and control.

[61]  P. Vincent,et al.  Blood flow in the rabbit aortic arch and descending thoracic aorta , 2011, Journal of The Royal Society Interface.

[62]  Justin S. Tran,et al.  Effect of Wall Elasticity on Hemodynamics and Wall Shear Stress in Patient-Specific Simulations in the Coronary Arteries. , 2019, Journal of biomechanical engineering.