Linear propagation of pulsatile waves in viscoelastic tubes

An experimental and theoretical analysis is made of pulsatile wave propagation in deformable latex tubes as a model of the propagation of pressure pulses in arteries. A quasi one-dimensional linear model is used in which, in particular, attention is paid to the viscous phenomena in fluid and tube wall. The agreement between experimental and theoretical results is satisfactory. It appeared that the viscoelastic behaviour of the tube wall dominates the damping of the pressure pulse. Several linear models are used to describe the wall behaviour. No significant differences between the results of these models were found.     

Influence of geometric taper on the derivation of the true propagation coefficient using a three point method

We studied the effect of geometric taper on the derivation of the true propagation coefficient from three pressures determined 10 cm apart ('three-point method'). For this purpose the true propagation coefficients of a uniform latex tube (length 50 cm, outer diameter 12.73 mm, Womersley phase velocity 6.23-6.42 ms-1 (1-10 Hz), Womersley damping coefficient 0.05-0.14 m-1 (1-10 Hz) and of a tapered latex tube (length 50 cm, outer diameter varying from 15.88 to 9.45 mm, in the middle section with same properties as the uniform tube) were determined. The real part of the true propagation coefficient (the damping coefficient) was compared with apparent damping, and with the damping coefficient calculated using Womersley's theory. The imaginary part of the true propagation coefficient (the phase coefficient) was expressed in terms of phase velocity. True phase velocity was compared with measurements of apparent phase velocity, foot-to-foot velocity, and calculations of phase velocity parameters Womersley's theory and the Moens-Korteweg equation. The results show that in the uniform tube the three-point propagation coefficient is in agreement with all other estimates. Taper causes an error in the three-point propagation coefficient. At some frequencies the damping is reversed to amplification (values up to -2 m-1) and the phase velocity may be both overestimated or underestimated (up to 50%). The overestimation of true damping as reported in the literature cannot be explained from vessel taper.     

Wave Propagation in a Viscous Fluid Contained in an Orthotropic Elastic Tube

The problem of pressure wave propagation through a viscous fluid contained in an orthotropic elastic tube is considered in connection with arterial blood flow. Solutions to the fluid flow and elasticity equations are obtained for the presence of a reflected wave. Numerical results are presented for both isotropic and orthotropic elastic tubes. In particular, the pressure pulse, flow rate, axial fluid velocity, and wall displacements are plotted vs. time at various stations along the ascending aorta of man. The results indicate an increase in the peak value of the pressure pulse and a decrease in the flow rate as the pulse propagates away from the heart. Finally, the velocity of wave propagation depends mainly on the tangential modulus of elasticity of the arterial wall, and anisotropy of the wall accounts in part for the reduction of longitudinal movements and an increase in the hydraulic resistance.     

Theoretical analysis of pressure pulse propagation in arterial vessels.

An original mathematical model of viscous fluid motion in a tapered and distensible tube is presented. The model equations are deduced by assuming a two-dimensional flow and taking into account the nonlinear terms in the fluid motion equations, as well as the nonlinear deformation of the tube wall. One distinctive feature of the model is the formal integration with respect to the radial coordinate of the Navier-Stokes equations by power series expansion. The consequent computational frame allows an easy, accurate evaluation of the effects produced by changing the values of all physical and geometrical tube parameters. The model is employed to study the propagation along an arterial vessel of a pressure pulse produced by a single flow pulse applied at the proximal vessel extremity. In particular, the effects of the natural taper angle of the arterial wall on pulse propagation are investigated. The simulation results show that tapering considerably influences wave attenuation but not wave velocity. The substantially different behavior of pulse propagation, depending upon whether it travels towards the distal extremity or in the opposite direction, is observed: natural tapering causes a continuous increase in the pulse amplitude as it moves towards the distal extremity; on the contrary, the reflected pulse, running in the opposite direction, is greatly damped. For a vessel with physical and geometrical properties similar to those of a canine femoral artery and 0.1 degree taper angle, the forward amplification is about 0.9 m-1 and the backward attenuation is 1.4 m-1, so that the overall tapering effect gives a remarkably damped pressure response. For a natural taper angle of 0.14 degrees the perturbation is almost extinct when the pulse wave returns to the proximal extremity.     

Wave Propagation through a Viscous Incompressible Fluid Contained in an Initially Stressed Elastic Tube

To have a better understanding of the flow of blood in arteries a theoretical analysis of the pressure wave propagation through a viscous incompressible fluid contained in an initially stressed tube is considered. The fluid is assumed to be Newtonian. The tube is taken to be elastic and isotropic. The analysis is restricted to tubes with thin walls and to waves whose wavelengths are very large compared with the radius of the tube. It is further assumed that the amplitude of the pressure disturbance is sufficiently small so that nonlinear terms of the inertia of the fluid are negligible compared with linear ones. Both circumferential and longitudinal initial stresses are considered; however, their origins are not specified. Initial stresses enter equations as independent parameters. A frequency equation, which is quadratic in the square of the propagation velocity is obtained. Two out of four roots of this equation give the velocity of propagation of two distinct outgoing waves. The remaining two roots represent incoming waves corresponding to the first two waves. One of the waves propagates more slowly than the other. As the circumferential and/or longitudinal stress of the wall increases, the velocity of propagation and transmission per wavelength of the slower wave decreases. The response of the fast wave to a change in the initial stress is on the opposite direction.     

Periodic flow of a viscous liquid in a thick-walled elastic tube

Arterial blood flow is analyzed on the basis of a realistic model consisting of a viscous liquid contained in a thick-walled viscoelastic tube. Approximate forms of the Navier-Stokes and continuity equations are derived for this model and solved in conjunction with the equations of motion of an elastic solid. Expressions are found for the displacement of the tube wall, velocity distribution, volume flow rate and phase velocity of the pressure wave. Changes in the shape of the pressure wave caused by damping and dispersion are determined, and the effect of viscoelasticity is assessed. Numerical results are presented which correspond to observed parameters of the circulatory systems of living animals.     

Heterogeneous mechanics of the mouse pulmonary arterial network

Individualized modeling and simulation of blood flow mechanics find applications in both animal research and patient care. Individual animal or patient models for blood vessel mechanics are based on combining measured vascular geometry with a fluid structure model coupling formulations describing dynamics of the fluid and mechanics of the wall. For example, one-dimensional fluid flow modeling requires a constitutive law relating vessel cross-sectional deformation to pressure in the lumen. To investigate means of identifying appropriate constitutive relationships, an automated segmentation algorithm was applied to micro-computerized tomography images from a mouse lung obtained at four different static pressures to identify the static pressure–radius relationship for four generations of vessels in the pulmonary arterial network. A shape-fitting function was parameterized for each vessel in the network to characterize the nonlinear and heterogeneous nature of vessel distensibility in the pulmonary arteries. These data on morphometric and mechanical properties were used to simulate pressure and flow velocity propagation in the network using one-dimensional representations of fluid and vessel wall mechanics. Moreover, wave intensity analysis was used to study effects of wall mechanics on generation and propagation of pressure wave reflections. Simulations were conducted to investigate the role of linear versus nonlinear formulations of wall elasticity and homogeneous versus heterogeneous treatments of vessel wall properties. Accounting for heterogeneity, by parameterizing the pressure/distention equation of state individually for each vessel segment, was found to have little effect on the predicted pressure profiles and wave propagation compared to a homogeneous parameterization based on average behavior. However, substantially different results were obtained using a linear elastic thin-shell model than were obtained using a nonlinear model that has a more physiologically realistic pressure versus radius relationship.     

Periodic flow of a viscous liquid in a thick-walled elastic tube

Arterial blood flow is analyzed on the basis of a realistic model consisting of a viscous liquid contained in a thick-walled viscoelastic tube. Approximate forms of the Navier-Stokes and continuity equations are derived for this model and solved in conjunction with the equations of motion of an elastic solid. Expressions are found for the displacement of the tube wall, velocity distribution, volume flow rate and phase velocity of the pressure wave. Changes in the shape of the pressure wave caused by damping and dispersion are determined, and the effect of viscoelasticity is assessed. Numerical results are presented which correspond to observed parameters of the circulatory systems of living animals. This research was partially supported by the National Science Foundation; it was done in part by D. K. Whirlow in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Carnegie Institute of Technology.     

Wave Propagation through a Newtonian Fluid Contained within a Thick-Walled,Viscoelastic Tube

The propagation of harmonic pressure waves through a Newtonian fluid contained within a thick-walled, viscoelastic tube is considered as a model of arterial blood flow. The fluid is assumed to be homogeneous and Newtonian, and its motion to be laminar and axisymmetric. The wall is assumed to be isotropic, incompressible, linear, and viscoelastic. It is also assumed that the motion is such that the convective acceleration is negligible. The motion of the fluid is described by the linearized form of the Navier-Stokes equations and the motion of the wall by classical elasticity theory. The frequency dependence of the wall mechanical properties are represented by a three parameter, relaxation-type model. Using boundary conditions describing the continuity of stress and velocity components in the fluid and the wall, explicit solutions for the system of equations of the model have been obtained. The longitudinal fluid impedance has been expressed in terms of frequency and the system parameters. The frequency equation has been solved and the propagation constant also expressed in terms of frequency and system parameters. The results indicate that the fluid impedance is smaller than predicted by the rigid tube model or by Womersley''s constrained elastic tube model. Also, the velocity of propagation is generally slower and the transmission per wavelength less than predicted by Womersley''s elastic tube model. The propagation constant is very sensitive to changes in the degree of wall viscoelasticity.     

Energy dissipation and pulse wave attenuation in the canine carotid artery

This paper presents methods for calculating, for a segment of artery in vivo, (1) viscous and viscoelastic energy dissipation as a function of time, and (2) the viscous and viscoelastic components of the frequency-dependent attenuation coefficient. The calculations require measurement of arterial diameter and of intra-arterial pressure and flow-rate at two sites along the vessel. Viscous energy dissipation is calculated from the radius-dependent velocity shear in the lumen given by linear theory from the pressure measurements. The attenuation coefficient for a given harmonic of heart rate is calculated as half the sum of the viscous and viscoelastic components at that frequency of the energy dissipated per unit length by forwardtravelling waves, divided by the forward-wave flow work input to the segment at that frequency. Measurements in canine carotid arteries indicate that wall viscoelasticity contributes relatively little to energy dissipation per cardiac cycle and pulse wave attenuation.     

Wave Propagation through a Viscous Fluid Contained in a Tethered, Initially Stressed, Orthotropic Elastic Tube

To give a realistic representation of the pulse propagation in arteries a theoretical analysis of the wave propagation through a viscous incompressible fluid contained in an initially stressed elastic tube is considered. The tube is assumed to be orthotropic and its longitudinal motion is constrained by a uniformly distributed additional mass, a dashpot and a spring. The fluid is assumed to be Newtonian. The analysis is restricted to propagation of small amplitude harmonic waves whose wavelength is large compared to the radius of the vessel. Elimination of arbitrary constants from the general solutions of the equations of motion of the fluid and the wall gives a frequency equation to determine the velocity of propagation. Two roots of this equation give the velocity of propagation of two distinct outgoing waves. One of the waves propagates slower than the other. The propagation properties of s lower waves are very slightly affected by the degree of anisotropy of the wall. The velocity of propagation of faster waves decreases as the ratio of the longitudinal modulus of elasticity to the circumferential modulus decreases; transmission of these waves is very little affected. The influence of the tethering on the propagation velocity of slower waves is negligibly small; transmission of these waves is seriously affected. In tethered tubes faster waves are completely attenuated.     

Comparative study of viscoelastic arterial wall models in nonlinear one-dimensional finite element simulations of blood flow

It is well known that blood vessels exhibit viscoelastic properties, which are modeled in the literature with different mathematical forms and experimental bases. The wide range of existing viscoelastic wall models may produce significantly different blood flow, pressure, and vessel deformation solutions in cardiovascular simulations. In this paper, we present a novel comparative study of two different viscoelastic wall models in nonlinear one-dimensional (1D) simulations of blood flow. The viscoelastic models are from papers by Holenstein et al. in 1980 (model V1) and Valdez-Jasso et al. in 2009 (model V2). The static elastic or zero-frequency responses of both models are chosen to be identical. The nonlinear 1D blood flow equations incorporating wall viscoelasticity are solved using a space-time finite element method and the implementation is verified with the Method of Manufactured Solutions. Simulation results using models V1, V2 and the common static elastic model are compared in three application examples: (i) wave propagation study in an idealized vessel with reflection-free outflow boundary condition; (ii) carotid artery model with nonperiodic boundary conditions; and (iii) subject-specific abdominal aorta model under rest and simulated lower limb exercise conditions. In the wave propagation study the damping and wave speed were largest for model V2 and lowest for the elastic model. In the carotid and abdominal aorta studies the most significant differences between wall models were observed in the hysteresis (pressure-area) loops, which were larger for V2 than V1, indicating that V2 is a more dissipative model. The cross-sectional area oscillations over the cardiac cycle were smaller for the viscoelastic models compared to the elastic model. In the abdominal aorta study, differences between constitutive models were more pronounced under exercise conditions than at rest. Inlet pressure pulse for model V1 was larger than the pulse for V2 and the elastic model in the exercise case. In this paper, we have successfully implemented and verified two viscoelastic wall models in a nonlinear 1D finite element blood flow solver and analyzed differences between these models in various idealized and physiological simulations, including exercise. The computational model of blood flow presented here can be utilized in further studies of the cardiovascular system incorporating viscoelastic wall properties.     

On the propagation of a disturbance through a viscous liquid flowing in a distensible tube of appreciable mass

The velocity of propagation of a disturbance wave in a liquid flowing in a distensible tube is computed. The mathematical model is more general than those used in previous analyses: the tube wall properties are realistic; the convective part of the axial inertia forces is taken into account; radial inertia forces of both the fluid and tube wall are present; viscous stresses are present. Four parameters influencing the velocity of propagation are obtained and discussed. Curves are plotted illustrating the effects of the parameters. Contrary to the results of previous analyses, viscous effects are shown to be appreciable in blood flow. It is also shown that radial inertia effects can be important in laboratory set-ups. The material presented in this paper was adapted from the Ph.D. thesis written by the author at Harvard University.     

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

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.     

Examination of elastic non-uniformity in the arterial system using a hydraulic model

Pressure wave propagation has been examined in a model artery with spatially varying compliance. Although results were affected by viscous losses, appropriate allowance for such losses produced agreement between experimental findings and predictions of linear wave transmission theory. Particularly, the ability of non-uniformity of the tube wall to generate amplification of the pressure wave was confirmed. However, extrapolation to the physiological situation suggests that reflections from discrete sites in peripheral beds have a greater effect on pressure wave propagation than does elastic non-uniformity of major vessels. A theoretical analysis has demonstrated that the effects of elastic non-uniformity can be interpreted as the integrated effects of infinitesimal reflections from each progressive increment in wall stiffness.     

Mechanics of blood flow

The historical development of the mechanics of blood flow can be traced from ancient times, to Leonardo da Vinci and Leonhard Euler and up to the present times with increasing biological knowledge and mathematical analysis. In the last two decades, quantitative and numerical methods have steadily given more complete and precise understanding. In the arterial system wave propagation computations based on nonlinear one-dimensional modeling have given the best representation of pulse wave propagation. In the veins, the theory of unsteady flow in collapsible tubes has recently been extensively developed. In the last decade, progress has been made in describing the blood flow at junctions, through stenoses, in bends and in capillary blood vessels. The rheological behavior of individual red blood cells has been explored. A working model consists of an elastic membrane filled with viscous fluid. This model forms a basis for understanding the viscous and viscoelastic behavior of blood.     

On the velocity of propagation of pressure waves in an incompressible viscous fluid enclosed in a tube with an elastomeric wall

Rashevsky's treatment of the flow of an incompressible viscous fluid in an elastic distensible tube is applied to the same problem, except that the wall of the tube is assumed to be elastomeric. As a result the velocity of propagation is obtained in terms of the elastomeric constants of the wall, the thickness and density of the wall, the viscosity of the fluid, and the radius of the tube.     

Inflation waves induced by axial acceleration of the aorta

The phenomenon of high-amplitude inflation waves resulting from a sharp axial acceleration of the aorta, as may occur in road accidents, is investigated theoretically. The aorta is modeled as an axisymmetric tapered membranic shell (tube) made of an incompressible, nonlinear viscoelastic material with cylindrical orthotropy. It is filled with an inviscid, incompressible fluid whose flow is considered as quasi-one dimensional along the tube axis. The equations of motion of the tube and of the fluid are solved numerically, by using a two-step explicit scheme, for several axial acceleration profiles. The solutions shows that an inflation wave is generated and it propagates in opposite direction to that of the acceleration. The wall stresses, deformations and their time derivatives as well as fluid velocity and pressure are determined along the tube at different time intervals. Peak axial and circumferential stresses are high, with the latter far exceeding the former. These stresses may cause rupture of the aorta.     

Nonaxisymmetric wave propagation through a viscous fluid in a viscoelastic tube.

The eigensolution of a nonaxisymmetric linear wave propagation through a viscous fluid in a viscoelastic pipe is obtained. The solution can serve as a foundation for studying nonaxisymmetric blood flow. Some leading eigenvalues for two special cases are computed. The result shows that nonaxisymmetric modes decay fast during their propagation.     

A viscoelastic model for use in predicting arterial pulse waves

In nonlinear mathematical models of the arterial circulation, the viscoelasticity of the vessel walls has generally been neglected or only taken into account in a highly approximate manner. A new method is proposed to simulate the nonlinear viscoelastic properties of the wall material with the aid of a convolution integral of the creep function and the pressure history. With this simulation it is possible to properly describe the measured characteristics of arterial viscoelasticity. Moreover, it is utilized in a mathematical model of arterial pulse propagation to study the influence of the internal wall friction on the shape, amplitude and mean value of pressure and flow pulses. The corresponding predictions are in much better agreement with in-vivo measurements, especially for the distal part of the circulation, than those obtained without viscoelasticity.