A One-Dimensional Mathematical Model for Studying the Pulsatile Flow in Microvascular Networks

Techniques that model microvascular hemodynamics have been developed for decades. While the physiological significance of pressure pulsatility is acknowledged, most of the microcirculatory models use steady flow approaches. To theoretically study the extent and transmission of pulsatility in microcirculation, dynamic models need to be developed. In this paper, we present a one-dimensional model to describe the dynamic behavior of microvascular blood flow. The model is applied to a microvascular network from a rat mesentery. Intravital microscopy was used to record the morphology and flow velocities in individual vessel segments, and boundaries are defined according to the experimental data. The system of governing equations constituting the model is solved numerically using the discontinuous Galerkin method. An implicit integration scheme is adopted to increase computing efficiency. The model allows the simulation of the dynamic properties of blood flow in microcirculatory networks, including the pressure pulsatility (quantified by a pulsatility index) and pulse wave velocity (PWV). From the main input arteriole to the main output venule, the pulsatility index decreases by 66.7%. PWV obtained along arterioles declines with decreasing diameters, with mean values of 77.16, 25.31, and 8.30 cm/s for diameters of 26.84, 17.46, and 13.33 μm, respectively. These results suggest that the 1D model developed is able to simulate the characteristics of pressure pulsatility and wave propagation in complex microvascular networks.