Fast blood-flow simulation for large arterial trees containing thousands of vessels

Abstract Blood flow modelling has previously been successfully carried out in arterial trees to study pulse wave propagation using nonlinear or linear flow solvers. However, the number of vessels used in the simulations seldom grows over a few hundred. The aim of this work is to present a computationally efficient solver coupled with highly detailed arterial trees containing thousands of vessels. The core of the solver is based on a modified transmission line method, which exploits the analogy between electrical current in finite-length conductors and blood flow in vessels. The viscoelastic behaviour of the arterial-wall is taken into account using a complex elastic modulus. The flow is solved vessel by vessel in the frequency domain and the calculated output pressure is then used as an input boundary condition for daughter vessels. The computational results yield pulsatile blood pressure and flow rate for every segment in the tree. This solver is coupled with large arterial trees generated from a three-dimensional constrained constructive optimisation algorithm. The tree contains thousands of blood vessels with radii spanning ~1 mm in the root artery to ~30 μm in leaf vessels. The computation takes seconds to complete for a vasculature of 2048 vessels and less than 2 min for a vasculature of 4096 vessels on a desktop computer.

[1]  Nicolas P Smith,et al.  Structural morphology of renal vasculature. , 2006, American journal of physiology. Heart and circulatory physiology.

[2]  Peter J. Hunter,et al.  A Hybrid 1D and 3D Approach to Hemodynamics Modelling for a Patient-Specific Cerebral Vasculature and Aneurysm , 2009, MICCAI.

[3]  D. Bergel,et al.  The dynamic elastic properties of the arterial wall , 1961, The Journal of physiology.

[4]  A. J. Pullan,et al.  Geometric modeling of the human torso using cubic hermite elements , 2007, Annals of Biomedical Engineering.

[5]  A. Avolio,et al.  Multi-branched model of the human arterial system , 1980, Medical and Biological Engineering and Computing.

[6]  P. Cahill,et al.  Pulsatile Flow–Induced Angiogenesis: Role of Gi Subunits , 2002, Arteriosclerosis, thrombosis, and vascular biology.

[7]  K. Burrowes,et al.  Contribution of serial and parallel microperfusion to spatial variability in pulmonary inter- and intra-acinar blood flow. , 2010, Journal of applied physiology.

[8]  H. Pain,et al.  The physics of vibrations and waves , 1968 .

[9]  W. Fleig,et al.  Hepatic arterial flow volume and reserve in patients with cirrhosis: use of intra-arterial Doppler and adenosine infusion. , 1999, Gastroenterology.

[10]  A Noordergraaf,et al.  Arterial viscoelasticity: a generalized model. Effect on input impedance and wave travel in the systematic tree. , 1970, Journal of biomechanics.

[11]  D. Bergel.  The static elastic properties of the arterial wall , 1961, The Journal of physiology.

[12]  N. Stergiopulos,et al.  Assessment of distributed arterial network models , 1997, Medical and Biological Engineering and Computing.

[13]  I. Puri,et al.  Field-induced self-assembled ferrofluid aggregation in pulsatile flow , 2005 .

[14]  Harvey Ho,et al.  Modeling the hepatic arterial buffer response in the liver. , 2013, Medical engineering & physics.

[15]  P. J. Hunter,et al.  Generation of an Anatomically Based Three-Dimensional Model of the Conducting Airways , 2000, Annals of Biomedical Engineering.

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

[17]  M Zamir,et al.  Branching characteristics of human coronary arteries. , 1986, Canadian journal of physiology and pharmacology.

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

[19]  A. Henney,et al.  The virtual liver: a multidisciplinary, multilevel challenge for systems biology , 2012, Wiley interdisciplinary reviews. Systems biology and medicine.

[20]  J K Raines,et al.  A computer simulation of arterial dynamics in the human leg. , 1974, Journal of biomechanics.

[21]  H. Pain,et al.  The Physics of Vibrations and Waves: Pain/The Physics of Vibrations and Waves , 2005 .

[22]  Alfio Quarteroni,et al.  On the physical consistency between three-dimensional and one-dimensional models in haemodynamics , 2013, J. Comput. Phys..

[23]  Alicia M. Williams.  The Hydrodynamics of Ferrofluid Aggregates , 2008 .

[24]  Tobias Preusser,et al.  Analysis and Algorithmic Generation of Hepatic Vascular Systems , 2012, International journal of hepatology.

[25]  J. Womersley.  Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known , 1955, The Journal of physiology.

[26]  Subin Solomen,et al.  Biomechanics of Circulation , 2018 .

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

[28]  Ishwar K. Puri,et al.  Particle Transport in Therapeutic Magnetic Fields , 2014 .

[29]  M J Ackerman,et al.  The Visible Human Project , 1998, Proc. IEEE.

[30]  Stephane Cotin,et al.  Real-Time Modeling of Vascular Flow for Angiography Simulation , 2007, MICCAI.

[31]  W. Schreiner,et al.  Computer-optimization of vascular trees , 1993, IEEE Transactions on Biomedical Engineering.

[32]  W Schreiner,et al.  A three-dimensional model for arterial tree representation, generated by constrained constructive optimization , 1999, Comput. Biol. Medicine.