Application of a multicomponent model of convectional reaction-diffusion to description of glucose gradients in a neurovascular unit

A supply of glucose to a nervous tissue is fulfilled by a cerebrovascular network, and further diffusion is known to occur at both an arteriolar and a microvascular level. Despite a direct relation, a blood flow dynamic and reaction-diffusion of metabolites are usually considered separately in the mathematical models. In the present study they are coupled in a multiphysical approach which allows to evaluate the effects of capillary blood flow changes on near-vessels nutrient concentration gradients evidently. Cerebral blood flow (CBF) was described by the non-steady-state Navier-Stokes equations for a non-Newtonian fluid whose constitutive law is given by the Carreau model. A three-level organization of blood–brain barrier (BBB) is modelled by the flux dysconnectivity functions including densities and kinetic properties of glucose transporters. The velocity of a fluid flow in brain extracellular space (ECS) was estimated using Darcy’s law. The equations of reaction-diffusion with convection based on a generated flow field for continues and porous media were used to describe spatial-time gradients of glucose in the capillary lumen and brain parenchyma of a neurovascular unit (NVU), respectively. Changes in CBF were directly simulated using smoothing step-like functions altering the difference of intracapillary pressure in time. The changes of CBF cover both the decrease (on 70%) and the increase (on 50%) in a capillary flow velocity. Analyzing the dynamics of glucose gradients, it was shown that a rapid decrease of a capillary blood flow yields an enhanced level of glucose in a near-capillary nervous tissue if the contacts between astrocytes end-feet are not tight. Under the increased CBF velocities the amplitude of glucose concentration gradients is always enhanced. The introduced approach can be used for estimation of blood flow changes influence not only on glucose but also on other nutrients concentration gradients and for the modelling of distributions of their concentrations near blood vessels in other tissues as well.