Modelling dynamic changes in blood flow and volume in the cerebral vasculature

ABSTRACT The cerebral microvasculature plays a key role in the transport of blood and the delivery of nutrients to the cells that perform brain function. Although recent advances in experimental imaging techniques mean that its structure and function can be interrogated to very small length scales, allowing individual vessels to be mapped to a fraction of 1&mgr;m, these techniques currently remain confined to animal models. In‐vivo human data can only be obtained at a much coarser length scale, of order 1mm, meaning that mathematical models of the microvasculature play a key role in interpreting flow and metabolism data. However, there are close to 10,000 vessels even within a single voxel of size 1mm3. Given the number of vessels present within a typical voxel and the complexity of the governing equations for flow and volume changes, it is computationally challenging to solve these in full, particularly when considering dynamic changes, such as those found in response to neural activation. We thus consider here the governing equations and some of the simplifications that have been proposed in order more rigorously to justify in what generations of blood vessels these approximations are valid. We show that two approximations (neglecting the advection term and assuming a quasi‐steady state solution for blood volume) can be applied throughout the cerebral vasculature and that two further approximations (a simple first order differential relationship between inlet and outlet flows and inlet and outlet pressures, and matching of static pressure at nodes) can be applied in vessels smaller than approximately 1mm in diameter. We then show how these results can be applied in solving flow fields within cerebral vascular networks providing a simplified yet rigorous approach to solving dynamic flow fields and compare the results to those obtained with alternative approaches. We thus provide a framework to model cerebral blood flow and volume within the cerebral vasculature that can be used, particularly at sub human imaging length scales, to provide greater insight into the behaviour of blood flow and volume in the cerebral vasculature. HIGHLIGHTS Four modelling assumptions for the cerebral vasculature are derived and a new model proposed.Advection can be neglected.Blood volume can be assumed to be quasi‐steady state.Static pressure can be matched at bifurcations for vessels < 1mm in diameter.A first order flow‐pressure relationship can be used for vessels <1mm in diameter.

[1]  Stephen J. Payne,et al.  Investigating the effects of a penetrating vessel occlusion with a multi-scale microvasculature model of the human cerebral cortex , 2018, NeuroImage.

[2]  Franck Plouraboué,et al.  Cerebral Blood Flow Modeling in Primate Cortex , 2010, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[3]  Anna Devor,et al.  Oxygen advection and diffusion in a three- dimensional vascular anatomical network. , 2008, Optics express.

[4]  Anders M. Dale,et al.  A vascular anatomical network model of the spatio-temporal response to brain activation , 2008, NeuroImage.

[5]  Andreas A Linninger,et al.  Hematocrit Distribution and Tissue Oxygenation in Large Microcirculatory Networks , 2015, Microcirculation.

[6]  Patrick Jenny,et al.  Vascular Graph Model to Simulate the Cerebral Blood Flow in Realistic Vascular Networks , 2009, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[7]  A. Pries,et al.  Biophysical aspects of blood flow in the microvasculature. , 1996, Cardiovascular research.

[8]  S. Čanić,et al.  Mathematical analysis of the quasilinear effects in a hyperbolic model blood flow through compliant axi‐symmetric vessels , 2003 .

[9]  A. Linninger,et al.  Cerebral Microcirculation and Oxygen Tension in the Human Secondary Cortex , 2013, Annals of Biomedical Engineering.

[10]  Stephen J Payne,et al.  A statistical model of the penetrating arterioles and venules in the human cerebral cortex , 2016, Microcirculation.

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

[12]  D. Kleinfeld,et al.  Fluctuating and sensory-induced vasodynamics in rodent cortex extend arteriole capacity , 2011, Proceedings of the National Academy of Sciences.

[13]  Michael Breakspear,et al.  Deconvolution of neural dynamics from fMRI data using a spatiotemporal hemodynamic response function , 2014, NeuroImage.

[14]  J. Alastruey,et al.  A Novel Analytical Approach to Pulsatile Blood Flow in the Arterial Network , 2016, Annals of Biomedical Engineering.

[15]  G. Langewouters,et al.  The static elastic properties of 45 human thoracic and 20 abdominal aortas in vitro and the parameters of a new model. , 1984, Journal of biomechanics.

[16]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[17]  Maiken Nedergaard,et al.  Two-photon NADH imaging exposes boundaries of oxygen diffusion in cortical vascular supply regions , 2010, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

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

[19]  Stephen J. Payne,et al.  Physiology for Engineers , 2020, Biosystems & Biorobotics.

[20]  Y. Fung,et al.  Mechanics of the Circulation , 2011, Developments in Cardiovascular Medicine.

[21]  A. Sbarbati,et al.  The microvascular system in ischemic cortical lesions , 1996, Acta Neuropathologica.

[22]  David Kleinfeld,et al.  The capillary bed offers the largest hemodynamic resistance to the cortical blood supply , 2017, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[23]  D. Kleinfeld,et al.  Topological basis for the robust distribution of blood to rodent neocortex , 2010, Proceedings of the National Academy of Sciences.

[24]  A. Huk,et al.  A Model for Transient Oxygen Delivery in Cerebral Cortex , 2009, Front. Neuroenerg..

[25]  D. Kleinfeld,et al.  Correlations of Neuronal and Microvascular Densities in Murine Cortex Revealed by Direct Counting and Colocalization of Nuclei and Vessels , 2009, The Journal of Neuroscience.

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

[27]  Gianni Pedrizzetti,et al.  Cardiovascular fluid mechanics , 2003 .

[28]  Giorgio A. Ascoli,et al.  Morphometric, geographic, and territorial characterization of brain arterial trees , 2014, International journal for numerical methods in biomedical engineering.

[29]  Patrick Jenny,et al.  Depth-dependent flow and pressure characteristics in cortical microvascular networks , 2017, PLoS Comput. Biol..

[30]  Anders M. Dale,et al.  Depth-resolved optical imaging and microscopy of vascular compartment dynamics during somatosensory stimulation , 2007, NeuroImage.

[31]  J. Mayhew,et al.  A Model of the Dynamic Relationship between Blood Flow and Volume Changes during Brain Activation , 2004, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[32]  Johannes Reichold,et al.  The microvascular system of the striate and extrastriate visual cortex of the macaque. , 2008, Cerebral cortex.

[33]  Mathieu Sellier,et al.  A computational model of hemodynamic parameters in cortical capillary networks. , 2011, Journal of theoretical biology.

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

[35]  D. Acheson Elementary Fluid Dynamics , 1990 .

[36]  Anna Devor,et al.  Quantifying the Microvascular Origin of BOLD-fMRI from First Principles with Two-Photon Microscopy and an Oxygen-Sensitive Nanoprobe , 2015, The Journal of Neuroscience.

[37]  S. Payne,et al.  Multi-scale homogenization of blood flow in 3-dimensional human cerebral microvascular networks. , 2015, Journal of theoretical biology.

[38]  S. Payne Cerebral Blood Flow and Metabolism: A Quantitative Approach , 2017 .

[39]  M. Holmes Introduction to Perturbation Methods , 1995 .

[40]  A. Hudetz,et al.  A semi-empirical model of apparent blood viscosity as a function of vessel diameter and discharge hematocrit. , 1991, Biorheology.

[41]  C. Lucas An anatomical model of the cerebral vasculature and blood flow , 2013 .

[42]  John E. W. Mayhew,et al.  A time-invariant visco-elastic windkessel model relating blood flow and blood volume , 2009, NeuroImage.

[43]  Stephen Payne,et al.  The Influence of Network Structure on the Transport of Blood in the Human Cerebral Microvasculature , 2012, Microcirculation.

[44]  S. Lorthois,et al.  Simulation study of brain blood flow regulation by intra-cortical arterioles in an anatomically accurate large human vascular network. Part II: Flow variations induced by global or localized modifications of arteriolar diameters , 2011, NeuroImage.

[45]  A. Pries,et al.  Blood viscosity in tube flow: dependence on diameter and hematocrit. , 1992, The American journal of physiology.

[46]  Rebecca J. Shipley,et al.  Multiscale Modelling of Fluid and Drug Transport in Vascular Tumours , 2010, Bulletin of mathematical biology.

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

[48]  R. Buxton,et al.  Dynamics of blood flow and oxygenation changes during brain activation: The balloon model , 1998, Magnetic resonance in medicine.

[49]  Stephen J. Payne,et al.  A generalized mathematical framework for estimating the residue function for arbitrary vascular networks , 2013, Interface Focus.

[50]  Stephen John Payne,et al.  Oxygen delivery from the cerebral microvasculature to tissue is governed by a single time constant of approximately 6 seconds , 2018, Microcirculation.

[51]  A. Pries,et al.  Blood flow in microvascular networks. Experiments and simulation. , 1990, Circulation research.

[52]  P A Robinson,et al.  Spatiotemporal hemodynamic response functions derived from physiology. , 2014, Journal of theoretical biology.

[53]  Seong-Gi Kim,et al.  Cerebral Oxygen Delivery and Consumption During Evoked Neural Activity , 2010, Front. Neuroenerg..

[54]  David Ress,et al.  Model of the Transient Neurovascular Response Based on Prompt Arterial Dilation , 2013, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[55]  Yutaka Tomita,et al.  RBC velocities in single capillaries of mouse and rat brains are the same, despite 10-fold difference in body size , 2010, Brain Research.

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

[57]  H. Duvernoy,et al.  Cortical blood vessels of the human brain , 1981, Brain Research Bulletin.

[58]  Francis Cassot,et al.  Morphometry of the human cerebral cortex microcirculation: General characteristics and space-related profiles , 2008, NeuroImage.

[59]  D. F. Young,et al.  Computer simulation of arterial flow with applications to arterial and aortic stenoses. , 1992, Journal of biomechanics.