12.1 Synopsis

1D modelling of the vasculature (arterial network) represents and insightful and efficient tool for tackling problems encountered in arterial biomechanics as well as other engineering problems. In particular, 3D modelling of the vasculature is relatively expensive. 1D modelling provides an alternative in which the modelling assumptions provide a good balance between physiological accuracy and computational efficiency.To describe the flow and pressure in this network we consider the conservation of mass and momentum applied to an impermeable, deformable tube filled with an incompressible fluid, the nonlinear system of partial differential equations presented in non-conservative form is given by

∂U--+ H ∂U--= S
 ∂t     ∂x
(12.1)

     [  ]         [       ]        [         ]
      U             U   A             ( 0   )
U  =   A  ,  H =   ρ∂p- U  ,  S =   1  f-− s
                    ∂A              ρ  A

in which A is the Area (related to pressure), x is the axial coordinate along the vessel, U(x,t) the axial velocity, P(x,t) is the pressure in the tube, ρ is the density and finally f the frictional force per unit length. The unknowns in equation 12.1 are u,A and p; hence we must provide an explicit algebraic relationship to close this system. Typically, closure is provided by an algebraic relationship between A and p. For a thin elastic tube this is given by

          (          )          √-----
p = p + β  √A-−  ∘A-- ,  β =  ---πhE----
    0               0         (1− ν2)A0
(12.2)

where p0 is the external pressure, A0 is the initial cross-sectional area, E is the Young’s modulus, h is the vessel wall thickness and ν is the Poisson’s ratio. Other more elaborate pressure - area relationships are currently being implemented into the framework. Application of Riemann’s method of characteristics to equations 12.1 and 12.2 indicates that velocity and area are propagated through the system by forward and backward travelling waves. These waves are reflected and within the network by appropriate treatment of interfaces and boundaries. In the following, we will explain the usage of the blood flow solver on the basis of a single-artery problem and also on an arterial network consisting of 55 arteries.