
2.1 Flow in (large) compliant vessels
The exchange of forces between blood and vascular wall and its resulting displacement can be eval-
uated by employing detailed three-dimensional (3D) fluid-structure interaction models [19] or sim-
plified one-dimensional (1D) blood flow modeling approaches which consider field variations only
along the main flow (axial) direction [20]. The latter methodology is less accurate (especially around
localized anatomical details), but has many competitive computational advantages and can be easily
adopted for vessel networks [21]. This makes 1D blood flow modeling the best option for describing
the fluid mechanics in reactive hyperemia, which involves propagating phenomena along the vascu-
lature with time scale of many seconds. Furthermore, through 1D blood flow modeling it is also
possible to perform wave intensity analysis which allows the quantification of pressure waveforms
travelling forward to the microcirculation and backward to the heart [20]. In arteries (and arterioles)
the blood behaviour can be generally approximated as homogeneous and Newtonian since the size of
the red blood cells carried by the plasma is considerably smaller (10 times) than the vessel diameter.
Furthermore, flow is also generally considered incompressible and laminar, with a Poiseuille velocity
profile. In 1D blood flow modeling, each compliant vessel can be treated as axisymmetric, with
blood velocity (u), pressure (P) and flow described as continuous variables along its axial direction
z. These quantities are averaged over the cross-sectional area (A) and their variation along the radial
direction is considered negligible. The Navier-Stokes equations for 1D blood flow in compliant vessels
can be expressed in terms of cross-sectional area and velocity averaged over the cross-section:
∂A
∂t +∂(Au)
∂z = 0,
∂u
∂t +u∂u
∂z +1
ρ
∂P
∂z +µu
ρ
8π
A= 0,
(1)
where tis time, µis the fluid dynamic viscosity, ρis the fluid density, while Q=Au is the (volumetric)
flow rate. It is worth noting that (1) can also be written in terms of flow rate and pressure [22] or
cross-sectional area and flow rate [23]. The mechanics underlying the vascular wall deformation
appears to be complex, mainly due to vessel visco-elastic properties and the capacity to produce
active tone for diameter regulation. To describe the interaction between blood and vessel wall,
different approaches can be used [24, 25, 26]. In the simplest case, the fluid pressure is related to the
cross-section of the vessel via a linear function with respect to the luminal diameter (D=p4A/π)
P=Pext +β(√A−√A0),(2)
where Pext is the external pressure from the surrounding tissue, βis a parameter representing the
wall elasticity and A0is the unstressed cross-section area. However, modeling reactive hyperemia
requires a vessel wall model able to describe the hyperpolarization-induced dilation and then the
resulting (compliant) structural response to hyperemic flow. Given such complexity, wall mechanics
models derived from conservation laws retaining mechanobiological features [27, 28] are preferable
over tube laws which are purely phenomenological; indeed the latter, for this specific application,
would require the introduction of several non-physical parameters. It is also noted that, if desired,
wall viscoelasticity can be integrated into the blood pressure-wall deformation law by using more
complex constitutive models [29, 30, 31] with a consequent decrease in computational efficiency. The
haemodynamic features (viscosity and density) and the geometric (diameter and length) and struc-
tural (stiffness) wall properties can be made vessel specific and can reflect different physiological and
pathological conditions such as ageing or hypertension [20, 32]. The blood flow variables described
by system (1) and (2) can be computed by employing a broad variety of numerical schemes including
finite differences, finite elements and finite volumes and can be extended to large vessel networks by
imposing mass and momentum conservation at the interface between vessels [33, 20, 34, 22].
2.2 Biphasic flow in microvessels
During reactive hyperemia, flow in the microvessels is altered from its physiological range due first
to the sudden pressure reduction induced by the upstream occlusion and then by arteriolar luminal
expansion consequent to the wall relaxation driven by ischemic tissues. Due to its morphology and
function, the downstream vasculature constitutes the site of major blood pressure drop along the
cardiovascular system. The work by Secomb [35] provides a detailed characterization of the flow
through microcirculatory networks. In these microvessels Reynolds number is <1, and therefore the
blood flow can be described with a good level of accuracy as incompressible Stokes fluid, for which the
convective component is neglected. Blood flow can still be described by (1) but most assume a rigid
microvessel wall, which implies considering only the momentum conservation equation for relating
blood flow and pressure. Furthermore, the biphasic nature of flow requires a modeling framework
that accounts for the main rheological properties of its components. Since they are concentrated in
3