Computational simulation of blood flow and drug transportation in a large vasculature

When the blood flow is considered as steady Poiseuille flow, the governing Navier–Stokes equations can be reduced to the Hagen–Poiseuille law. This simplification results in a substantially reduced computational cost which is very useful in real-time blood flow simulations, in particular for large vasculatures containing thousands of blood vessels. By incorporating a convection equation we can also simulate transient drug transportation in the vasculature. In this paper we present at first the implementation of a real-time flow solver which is coupled with a complex arterial tree generated from the Constrained Constructive Optimisation (CCO) algorithm. The computational time for \(\sim \) 8200 vessels was 0.2 s. Secondly, we simulate the transient drug transportation in the vasculature. Thirdly, we model the delivery of the drug into a tiny tissue block by adopting a 3D diffusion equation. In conclusion the presented computational techniques constitute a pipeline for circulation modelling in multiple scales, and may be used in a variety of biomedical applications.

[1]  Karol Miller,et al.  Cellular automata coupled with steady‐state nutrient solution permit simulation of large‐scale growth of tumours , 2013, International journal for numerical methods in biomedical engineering.

[2]  Bernhard Preim,et al.  Analysis of vasculature for liver surgical planning , 2002, IEEE Transactions on Medical Imaging.

[3]  Kumar Mithraratne,et al.  Anatomically-based musculoskeletal modeling: prediction and validation of muscle deformation during walking , 2009, The Visual Computer.

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

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

[6]  Pablo J. Blanco,et al.  Assessing the influence of heart rate in local hemodynamics through coupled 3D‐1D‐0D models , 2010 .

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

[8]  A. Pries,et al.  Microvascular blood viscosity in vivo and the endothelial surface layer. , 2005, American journal of physiology. Heart and circulatory physiology.

[9]  Charles L. Weber,et al.  A Mathematical Model , 1987 .

[10]  Rotem Ben-Shachar,et al.  The biochemistry of acetaminophen hepatotoxicity and rescue: a mathematical model , 2012, Theoretical Biology and Medical Modelling.

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

[12]  Harvey Ho,et al.  A Computer Simulation for 3D Vasculature-Based Oxygen Distribution and Tumour Growth , 2015 .

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