One-dimensional model for propagation of a pressure wave in a model of the human arterial network: comparison of theoretical and experimental results

Abstract

Pulse wave evaluation is an effective method for arteriosclerosis screening. In a previous study, we verified that pulse waveforms change markedly due to arterial stiffness. However, a pulse wave consists of two components, the incident wave and multireflected waves. Clarification of the complicated propagation of these waves is necessary to gain an understanding of the nature of pulse waves in vivo. In this study, we built a one-dimensional theoretical model of a pressure wave propagating in a flexible tube. To evaluate the applicability of the model, we compared theoretical estimations with measured data obtained from basic tube models and a simple arterial model. We constructed different viscoelastic tube set-ups: two straight tubes; one tube connected to two tubes of different elasticity; a single bifurcation tube; and a simple arterial network with four bifurcations. Soft polyurethane tubes were used and the configuration was based on a realistic human arterial network. The tensile modulus of the material was similar to the elasticity of arteries. A pulsatile flow with ejection time 0.3 s was applied using a controlled pump. Inner pressure waves and flow velocity were then measured using a pressure sensor and an ultrasonic diagnostic system. We formulated a 1D model derived from the Navier-Stokes equations and a continuity equation to characterize pressure propagation in flexible tubes. The theoretical model includes nonlinearity and attenuation terms due to the tube wall, and flow viscosity derived from a steady Hagen-Poiseuille profile. Under the same configuration as for experiments, the governing equations were computed using the MacCormack scheme. The theoretical pressure waves for each case showed a good fit to the experimental waves. The square sum of residuals (difference between theoretical and experimental wave-forms) for each case was <10.0%. A possible explanation for the increase in the square sum of residuals is the approximation error for flow viscosity. However, the comparatively small values prove the validity of the approach and indicate the usefulness of the model for understanding pressure propagation in the human arterial network.