A computational framework for derivative-free optimization of cardiovascular geometries

Abstract Predictive simulation and optimization-based design tools have great potential to improve the design of surgeries and interventions used in cardiovascular medicine. The present work builds upon recent advances in blood flow simulation capabilities to develop tools for optimization. A framework for coupling optimal shape design to time-accurate three-dimensional blood flow simulations in idealized cardiovascular geometries is presented. The optimization method employed is a tailored version of the surrogate management framework, a method that was developed previously for expensive functions with little or no gradient information. In this method, we employ a derivative-free approach using surrogates for increased efficiency together with mesh adaptive direct search to guarantee convergence to a local minimum. The optimization procedure has been fully automated to include the generation and parameterization of three-dimensional solid geometries, mesh generation, finite element flow computation, post processing and optimization. This methodology is demonstrated on three model problems that are representative of typical cardiovascular geometries: a stenosis, a vessel bifurcation modeled on Murray’s law, and an end-to-side anastomosis. These problems are used to illustrate the importance of the choice of cost function and the performance of the optimization algorithm. This framework for treatment optimization will be applicable to a wide range of cardiovascular surgical or catheter-based interventions in future work.

[1]  M Zamir,et al.  Shear forces and blood vessel radii in the cardiovascular system , 1977, The Journal of general physiology.

[2]  Charles A. Taylor,et al.  Comparison of CFD and MRI Flow and Velocities in an In Vitro Large Artery Bypass Graft Model , 2005, Annals of Biomedical Engineering.

[3]  Charles Audet,et al.  Analysis of Generalized Pattern Searches , 2000, SIAM J. Optim..

[4]  K. Perktold,et al.  Computer simulation of local blood flow and vessel mechanics in a compliant carotid artery bifurcation model. , 1995, Journal of biomechanics.

[5]  K. Katada,et al.  Magnitude and Role of Wall Shear Stress on Cerebral Aneurysm: Computational Fluid Dynamic Study of 20 Middle Cerebral Artery Aneurysms , 2004, Stroke.

[6]  Alfio Quarteroni,et al.  Numerical Treatment of Defective Boundary Conditions for the Navier-Stokes Equations , 2002, SIAM J. Numer. Anal..

[7]  J F Antaki,et al.  Computational flow optimization of rotary blood pump components. , 1995, Artificial organs.

[8]  Thomas J. R. Hughes,et al.  Finite Element Modeling of Three-Dimensional Pulsatile Flow in the Abdominal Aorta: Relevance to Atherosclerosis , 2004, Annals of Biomedical Engineering.

[9]  B. Kuban,et al.  Validated computation of physiologic flow in a realistic coronary artery branch. , 1997, Journal of biomechanics.

[10]  Meng Wang,et al.  Optimal Aeroacoustic Shape Design Using the Surrogate Management Framework , 2003 .

[11]  Sven Leyffer,et al.  Nonlinear programming without a penalty function , 2002, Math. Program..

[12]  Parviz Moin,et al.  Suppression of vortex-shedding noise via derivative-free shape optimization , 2004 .

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

[14]  C. Taylor,et al.  Predictive medicine: computational techniques in therapeutic decision-making. , 1999, Computer aided surgery : official journal of the International Society for Computer Aided Surgery.

[15]  D. E. Strandness,et al.  Duplex scanning in vascular disorders , 1990 .

[16]  Virginia Torczon,et al.  On the Convergence of Pattern Search Algorithms , 1997, SIAM J. Optim..

[17]  V. Torczon,et al.  RANK ORDERING AND POSITIVE BASES IN PATTERN SEARCH ALGORITHMS , 1996 .

[18]  Robert W. Dutton,et al.  A Software Framework for Creating Patient Specific Geometric Models from Medical Imaging Data for Simulation Based Medical Planning of Vascular Surgery , 2001, MICCAI.

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

[20]  Antony Jameson,et al.  Aerodynamic design via control theory , 1988, J. Sci. Comput..

[21]  John F LaDisa,et al.  Alterations in wall shear stress predict sites of neointimal hyperplasia after stent implantation in rabbit iliac arteries. , 2005, American journal of physiology. Heart and circulatory physiology.

[22]  M. D. McKay,et al.  A comparison of three methods for selecting values of input variables in the analysis of output from a computer code , 2000 .

[23]  Robert Michael Lewis,et al.  A Globally Convergent Augmented Lagrangian Pattern Search Algorithm for Optimization with General Constraints and Simple Bounds , 2002, SIAM J. Optim..

[24]  Robin Shandas,et al.  Influence of connection geometry and SVC-IVC flow rate ratio on flow structures within the total cavopulmonary connection: a numerical study. , 2002, Journal of biomechanical engineering.

[25]  C G Caro,et al.  Abnormal arterial flow pattern in untreated essential hypertension: possible link with the development of atherosclerosis. , 1990, Clinical science.

[26]  B. Rutt,et al.  Hemodynamics of human carotid artery bifurcations: computational studies with models reconstructed from magnetic resonance imaging of normal subjects. , 1998, Journal of vascular surgery.

[27]  J F Antaki,et al.  Computational fluid dynamics as a development tool for rotary blood pumps. , 2001, Artificial organs.

[28]  T. Hughes,et al.  Effect of exercise on hemodynamic conditions in the abdominal aorta. , 1999, Journal of vascular surgery.

[29]  Gregory R. Grant,et al.  Coexisting proinflammatory and antioxidative endothelial transcription profiles in a disturbed flow region of the adult porcine aorta. , 2004 .

[30]  Mauro Grigioni,et al.  Power dissipation associated with surgical operations' hemodynamics: critical issues and application to the total cavopulmonary connection. , 2006, Journal of biomechanics.

[31]  Charles A. Taylor,et al.  A coupled momentum method for modeling blood flow in three-dimensional deformable arteries , 2006 .

[32]  Pascal Frey,et al.  Fluid-structure interaction in blood flows on geometries based on medical imaging , 2005 .

[33]  A. J. Booker,et al.  A rigorous framework for optimization of expensive functions by surrogates , 1998 .

[34]  Alastair J. Martin,et al.  Computational approach to quantifying hemodynamic forces in giant cerebral aneurysms. , 2003, AJNR. American journal of neuroradiology.

[35]  F. Migliavacca,et al.  Computational fluid dynamics simulations in realistic 3-D geometries of the total cavopulmonary anastomosis: the influence of the inferior caval anastomosis. , 2003, Journal of biomechanical engineering.

[36]  Charles A. Taylor,et al.  Effects of Exercise and Respiration on Hemodynamic Efficiency in CFD Simulations of the Total Cavopulmonary Connection , 2007, Annals of Biomedical Engineering.

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

[38]  J. Peraire,et al.  Practical Three-Dimensional Aerodynamic Design and Optimization Using Unstructured Meshes , 1997 .

[39]  John E. Dennis,et al.  A framework for managing models in nonlinear optimization of computationally expensive functions , 1999 .

[40]  C. Zarins,et al.  Relative contribution of wall shear stress and injury in experimental intimal thickening at PTFE end-to-side arterial anastomoses. , 2002, Journal of biomechanical engineering.

[41]  C. R. Ethier,et al.  Accuracy of Computational Hemodynamics in Complex Arterial Geometries Reconstructed from Magnetic Resonance Imaging , 2004, Annals of Biomedical Engineering.

[42]  M. Czerny,et al.  Fluid Dynamics, Wall Mechanics, and Oxygen Transfer in Peripheral Bypass Anastomoses , 2002, Annals of Biomedical Engineering.

[43]  Robert Michael Lewis,et al.  Pattern Search Methods for Linearly Constrained Minimization , 1999, SIAM J. Optim..

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

[45]  C. D. Murray.  THE PHYSIOLOGICAL PRINCIPLE OF MINIMUM WORK APPLIED TO THE ANGLE OF BRANCHING OF ARTERIES , 1926, The Journal of general physiology.

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

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

[48]  R Pietrabissa,et al.  Multiscale modelling as a tool to prescribe realistic boundary conditions for the study of surgical procedures. , 2002, Biorheology.

[49]  G. Hulbert,et al.  A generalized-α method for integrating the filtered Navier–Stokes equations with a stabilized finite element method , 2000 .

[50]  Thomas J. R. Hughes,et al.  Finite element modeling of blood flow in arteries , 1998 .

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

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

[53]  F. Migliavacca,et al.  Computational fluid dynamic simulations of cavopulmonary connections with an extracardiac lateral conduit. , 1999, Medical engineering & physics.

[54]  C. Zarins,et al.  Carotid Bifurcation Atherosclerosis: Quantitative Correlation of Plaque Localization with Flow Velocity Profiles and Wall Shear Stress , 1983, Circulation research.

[55]  D. Giddens,et al.  Pulsatile flow in an end-to-side vascular graft model: comparison of computations with experimental data. , 2001, Journal of biomechanical engineering.

[56]  Jaime Peraire,et al.  Practical 3D aerodynamic design and optimization using unstructured meshes , 1996 .

[57]  C M Putman,et al.  Computational fluid dynamics modeling of intracranial aneurysms: effects of parent artery segmentation on intra-aneurysmal hemodynamics. , 2006, AJNR. American journal of neuroradiology.

[58]  J. Dennis,et al.  Trailing-edge noise reduction using derivative-free optimization and large-eddy simulation , 2007, Journal of Fluid Mechanics.

[59]  Charles A. Taylor,et al.  In Vivo Validation of Numerical Prediction of Blood Flow in Arterial Bypass Grafts , 2002, Annals of Biomedical Engineering.

[60]  L A Taber,et al.  An optimization principle for vascular radius including the effects of smooth muscle tone. , 1998, Biophysical journal.

[61]  Chandler Davis.  THEORY OF POSITIVE LINEAR DEPENDENCE. , 1954 .

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

[63]  A. Leuprecht,et al.  Numerical study of hemodynamics and wall mechanics in distal end-to-side anastomoses of bypass grafts. , 2002, Journal of biomechanics.

[64]  Rainald Löhner,et al.  Blood flow modeling in carotid arteries with computational fluid dynamics and MR imaging. , 2002, Academic radiology.

[65]  Christopher P. Cheng,et al.  Abdominal aortic hemodynamics in young healthy adults at rest and during lower limb exercise: quantification using image-based computer modeling. , 2006, American journal of physiology. Heart and circulatory physiology.

[66]  Charles Audet,et al.  Mesh Adaptive Direct Search Algorithms for Constrained Optimization , 2006, SIAM J. Optim..

[67]  Søren Nymand Lophaven,et al.  DACE - A Matlab Kriging Toolbox, Version 2.0 , 2002 .

[68]  Charles Audet,et al.  A Pattern Search Filter Method for Nonlinear Programming without Derivatives , 2001, SIAM J. Optim..

[69]  M. Heinkenschloss,et al.  Shape optimization in steady blood flow: A numerical study of non-Newtonian effects , 2005, Computer methods in biomechanics and biomedical engineering.

[70]  David A. Steinman,et al.  Image-Based Computational Fluid Dynamics Modeling in Realistic Arterial Geometries , 2002, Annals of Biomedical Engineering.

[71]  John F LaDisa,et al.  Stent design properties and deployment ratio influence indexes of wall shear stress: a three-dimensional computational fluid dynamics investigation within a normal artery. , 2004, Journal of applied physiology.

[72]  Robert Michael Lewis,et al.  Pattern Search Algorithms for Bound Constrained Minimization , 1999, SIAM J. Optim..

[73]  Arthur C. Taylor,et al.  Three-dimensional aerodynamic shape sensitivity analysis and design optimization using the Euler equations on unstructured grids , 1996 .

[74]  Oktay Baysal,et al.  Three-dimensional aerodynamic shape optimization using discrete sensitivity analysis , 1996 .

[75]  Michael Schäfer,et al.  A numerical approach for shape optimization of fluid flow domains , 2005 .

[76]  R.W. Dutton,et al.  Improving geometric model construction for blood flow modeling , 1999, IEEE Engineering in Medicine and Biology Magazine.

[77]  Jaime Peraire,et al.  PROGRESS TOWARDS A 3D AERODYNAMIC SHAPE OPTIMIZATION TOOL FOR THE COMPRESSIBLE, ElGE-Re NAVIER-STOKES EQUATIONS DISCRETIZED ON UNSTRUCTURED MESHES. , 1998 .

[78]  Thomas J. R. Hughes,et al.  Computational investigations in vascular disease , 1996 .

[79]  A. Jameson,et al.  Optimum Aerodynamic Design Using the Navier–Stokes Equations , 1997 .