LARGE‐SCALE SIMULATION OF THE HUMAN ARTERIAL TREE

1 Full‐scale simulations of the virtual physiological human (VPH) will require significant advances in modelling, multiscale mathematics, scientific computing and further advances in medical imaging. Herein, we review some of the main issues that need to be resolved in order to make three‐dimensional (3D) simulations of blood flow in the human arterial tree feasible in the near future. 2 A straightforward approach is computationally prohibitive even on the emerging petaflop supercomputers, so a three‐level hierarchical approach based on vessel size is required, consisting of: (i) a macrovascular network (MaN); (ii) a mesovascular network (MeN); and (iii) a microvascular network (MiN). We present recent simulations of MaN obtained by solving the 3D Navier–Stokes equations on arterial networks with tens of arteries and bifurcations and accounting for the neglected dynamics through proper boundary conditions. 3 A multiscale simulation coupling MaN–MeN–MiN and running on hundreds of thousands of processors on petaflop computers will require no more than a few CPU hours per cardiac cycle within the next 5 years. The rapidly growing capacity of supercomputing centres opens up the possibility of simulation studies of cardiovascular diseases, drug delivery, perfusion in the brain and other pathologies.

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

[2]  Suchuan Dong,et al.  Cross-site computations on the TeraGrid , 2005, Computing in Science & Engineering.

[3]  G. Karniadakis,et al.  Spectral/hp Element Methods for CFD , 1999 .

[4]  M. Zamir On fractal properties of arterial trees. , 1999, Journal of theoretical biology.

[5]  Aaron L. Fogelson,et al.  Computational Modeling of Blood Clotting: Coagulation and Three-dimensional Platelet Aggregation , 2003 .

[6]  Richard N Upton,et al.  CEREBRAL UPTAKE OF DRUGS IN HUMANS , 2007, Clinical and experimental pharmacology & physiology.

[7]  Chi-Wang Shu,et al.  Discontinuous Galerkin Methods: Theory, Computation and Applications , 2011 .

[8]  Robert W. Dutton,et al.  Level set methods and MR image segmentation for geometric modeling in computational hemodynamics , 1998, Proceedings of the 20th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. Vol.20 Biomedical Engineering Towards the Year 2000 and Beyond (Cat. No.98CH36286).

[9]  D A Steinman,et al.  Computational blood flow modelling: errors associated with reconstructing finite element models from magnetic resonance images. , 1997, Journal of biomechanics.

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

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

[12]  I. Marshall,et al.  MRI and CFD studies of pulsatile flow in healthy and stenosed carotid bifurcation models. , 2004, Journal of biomechanics.

[13]  A. Quarteroni,et al.  On the coupling of 3D and 1D Navier-Stokes equations for flow problems in compliant vessels , 2001 .

[14]  S. Sherwin,et al.  One-dimensional modelling of a vascular network in space-time variables , 2003 .

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

[16]  Charles A. Taylor,et al.  Morphometry-Based Impedance Boundary Conditions for Patient-Specific Modeling of Blood Flow in Pulmonary Arteries , 2007, Annals of Biomedical Engineering.

[17]  Joaquim Peiró,et al.  Three-dimensional reconstruction of autologous vein bypass graft distal anastomoses imaged with magnetic resonance: clinical and research applications. , 2003, Journal of vascular surgery.

[18]  Ioannis Chatzis,et al.  Permeability and electrical conductivity of porous media from 3D stochastic replicas of the microstructure , 2000 .

[19]  S. Sherwin,et al.  Modelling the circle of Willis to assess the effects of anatomical variations and occlusions on cerebral flows. , 2007, Journal of biomechanics.

[20]  George Em Karniadakis,et al.  Distribution of WSS on the Internal Carotid Artery With an Aneurysm: A CFD Sensitivity Study , 2007 .

[21]  Juan R Cebral,et al.  Patient-specific computational modeling of cerebral aneurysms with multiple avenues of flow from 3D rotational angiography images. , 2006, Academic radiology.

[22]  F. Larsen,et al.  Persistent Arterial Hyperammonemia Increases the Concentration of Glutamine and Alanine in the Brain and Correlates with Intracranial Pressure in Patients with Fulminant Hepatic Failure , 2006, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[23]  Spencer J. Sherwin,et al.  Parallel performance of the coarse space linear vertex solver and low energy basis preconditioner for spectral/hp elements , 2009, Parallel Comput..

[24]  John G. Wood,et al.  MCDONALDʼS BLOOD FLOW IN ARTERIES: THEORETICAL, EXPERIMENTAL AND CLINICAL PRINCIPLES, 4TH EDITION , 1998 .

[25]  Ian Marshall,et al.  MRI measurement of time‐resolved wall shear stress vectors in a carotid bifurcation model, and comparison with CFD predictions , 2003, Journal of magnetic resonance imaging : JMRI.

[26]  G. Born,et al.  Aggregation of Blood Platelets by Adenosine Diphosphate and its Reversal , 1962, Nature.

[27]  S. Sherwin,et al.  Pulse wave propagation in a model human arterial network: Assessment of 1-D visco-elastic simulations against in vitro measurements , 2011, Journal of biomechanics.

[28]  D. Krause,et al.  CEREBROVASCULAR EFFECTS OF OESTROGEN: MULTIPLICITY OF ACTION , 2007, Clinical and experimental pharmacology & physiology.

[29]  R. Shadwick,et al.  Circulatory mechanics in the toad Bufo marinus. II, Haemodynamics of the arterial windkessel , 1991 .

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

[31]  Guy Marchal,et al.  Nonrigid registration for subtraction CT angiography applied to the carotids and cranial arteries. , 2007, Academic radiology.

[32]  W. Bautz,et al.  Clinical evaluation of bone-subtraction CT angiography (BSCTA) in head and neck imaging , 2006, European Radiology.

[33]  David J. Begley,et al.  The Blood-brain Barrier and Drug Delivery to the CNS , 2000 .

[34]  Céline Fouard,et al.  A Novel Three‐Dimensional Computer‐Assisted Method for a Quantitative Study of Microvascular Networks of the Human Cerebral Cortex , 2006, Microcirculation.

[35]  T F Sherman,et al.  On connecting large vessels to small. The meaning of Murray's law , 1981, The Journal of general physiology.

[36]  Alfio Quarteroni,et al.  Multiscale modelling of the circulatory system: a preliminary analysis , 1999 .

[37]  Nadia Magnenat-Thalmann,et al.  The SPHERIGON: a simple polygon patch for smoothing quickly your polygonal meshes , 1998, Proceedings Computer Animation '98 (Cat. No.98EX169).

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

[39]  L. Formaggia,et al.  Numerical modeling of 1D arterial networks coupled with a lumped parameters description of the heart , 2006, Computer methods in biomechanics and biomedical engineering.

[40]  Jing Wan,et al.  A One-dimensional Finite Element Method for Simulation-based Medical Planning for Cardiovascular Disease , 2002, Computer methods in biomechanics and biomedical engineering.

[41]  G. Karniadakis,et al.  Blood flow velocity effects and role of activation delay time on growth and form of platelet thrombi , 2006, Proceedings of the National Academy of Sciences.

[42]  H. Ehmke,et al.  Role of angiotensin II in dynamic renal blood flow autoregulation of the conscious dog , 2002, The Journal of physiology.

[43]  D. Saloner,et al.  Numerical analysis of flow through a severely stenotic carotid artery bifurcation. , 2002, Journal of biomechanical engineering.

[44]  P. Fischer,et al.  Simulation of high-Reynolds number vascular flows , 2007 .

[45]  Michael E. Papka,et al.  Simulating and visualizing the human arterial system on the TeraGrid , 2006, Future Gener. Comput. Syst..

[46]  S. Whitaker The method of volume averaging , 1998 .

[47]  E. Ritenour Doppler Ultrasound: Physics, Instrumentation and Clinical Applications , 1990 .

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

[49]  Sylvie Lorthois,et al.  Numerical Simulation of Magnetic Resonance Angiographies of an Anatomically Realistic Stenotic Carotid Bifurcation , 2005, Annals of Biomedical Engineering.

[50]  S. Sherwin,et al.  Pulse wave propagation in a model human arterial network: assessment of 1-D numerical simulations against in vitro measurements. , 2007, Journal of biomechanics.

[51]  R. M. C. So,et al.  The effect of inlet velocity profile on the bifurcation COPD airway flow , 2006, Comput. Biol. Medicine.

[52]  M. Zamir,et al.  The Physics of Pulsatile Flow , 2000, Biological Physics Series.

[53]  George E. Karniadakis,et al.  A Reconstruction Method for Gappy and Noisy Arterial Flow Data , 2007, IEEE Transactions on Medical Imaging.

[54]  Spencer J. Sherwin,et al.  High-order algorithms for vascular flow modelling , 2002 .

[55]  R. Aaslid,et al.  A transcranial Doppler method in the evaluation of cerebrovascular spasm , 2004, Neuroradiology.

[56]  Sean D. Peterson,et al.  Experimental Study of the Effect of a Skewed Inlet Flow Profile on Stenotic Flow Development , 2006 .

[57]  K. Parker,et al.  Wave propagation in a model of the arterial circulation. , 2004, Journal of biomechanics.

[58]  M. Olufsen,et al.  Numerical Simulation and Experimental Validation of Blood Flow in Arteries with Structured-Tree Outflow Conditions , 2000, Annals of Biomedical Engineering.

[59]  Ulrich Rüde,et al.  Is 1.7 x 10^10 Unknowns the Largest Finite Element System that Can Be Solved Today? , 2005, ACM/IEEE SC 2005 Conference (SC'05).

[60]  David A. Vorp,et al.  Computational modeling of arterial biomechanics: insights into pathogenesis and treatment of vascular disease. , 2003, Journal of vascular surgery.

[61]  B. Rutt,et al.  Reconstruction of carotid bifurcation hemodynamics and wall thickness using computational fluid dynamics and MRI , 2002, Magnetic resonance in medicine.

[62]  Alex M. Andrew,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science (2nd edition) , 2000 .

[63]  Paul F. Fischer,et al.  Terascale Spectral Element Algorithms and Implementations , 1999, ACM/IEEE SC 1999 Conference (SC'99).

[64]  Dongbin Xiu,et al.  Parametric uncertainty analysis of pulse wave propagation in a model of a human arterial network , 2007, J. Comput. Phys..

[65]  C. Putman,et al.  Patient-specific computational fluid dynamics modeling of anterior communicating artery aneurysms: a study of the sensitivity of intra-aneurysmal flow patterns to flow conditions in the carotid arteries. , 2006, AJNR. American journal of neuroradiology.

[66]  R. Aaslid,et al.  Noninvasive transcranial Doppler ultrasound recording of flow velocity in basal cerebral arteries. , 1982, Journal of neurosurgery.

[67]  A. Popel,et al.  A computational study of the effect of capillary network anastomoses and tortuosity on oxygen transport. , 2000, Journal of theoretical biology.

[68]  S. Sherwin,et al.  Lumped parameter outflow models for 1-D blood flow simulations: Effect on pulse waves and parameter estimation , 2008 .

[69]  R. Herfkens,et al.  Phase contrast cine magnetic resonance imaging. , 1991, Magnetic resonance quarterly.

[70]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid , 2012 .

[71]  H.M. Tufo,et al.  Terascale Spectral Element Algorithms and Implementations , 1999, ACM/IEEE SC 1999 Conference (SC'99).

[72]  S. Sherwin,et al.  Mesh generation in curvilinear domains using high‐order elements , 2002 .

[73]  Pablo J. Blanco,et al.  Multidimensional modelling for the carotid artery blood flow , 2006 .

[74]  W. Arendshorst,et al.  A novel mechanism of renal blood flow autoregulation and the autoregulatory role of A1 adenosine receptors in mice. , 2007, American journal of physiology. Renal physiology.

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