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.     

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.     

Note on wave propagation in a thin elastic tube containing a viscous fluid

This paper deals with the theoretical and experimental investigation of the propagation of pressure waves through a straight thin-walled elastic tube, filled with a streaming Newtonian fluid. The tube used was of silicon rubber and measured 200 mm in length, 8 mm in internal diameter and 1 mm in wall thickness. The velocity profiles across the diameter of the tube were measured at the beginning, at the middle and the end of the elastic pipe and compared with the results of a mathematical-computational model. The test data agree fairly well with the calculated values.     

Scaling phloem transport: elasticity and pressure-concentration waves

Earlier theoretical analyses of the rate of propagation of pressure-concentration waves in the phloem were performed without adequate attention to the elastic expansion of sieve tube walls. Here, it is shown that the rate of propagation of pressure-concentration waves in phloem sieve tubes is not significantly impeded by wall elasticity, but rather, as previously implicated, by the ratio of sap osmotic pressure to the axial drop in sap hydrostatic pressure. It is also shown that pressure-concentration waves move equally well in both the upstream and downstream directions. These results permit future models to ignore elastic effects, and lend additional theoretical support to the "osmoregulatory flow" hypothesis, which argues that efficient molecular control of the phloem is permitted by maintaining sieve sap hydrostatic pressure at a value that is spatially nearly constant, which in turn permits changes in sieve tube state to be rapidly transmitted throughout the sieve tube via pressure-concentration waves.     

Experimental validation of a time-domain-based wave propagation model of blood flow in viscoelastic vessels

Time-domain-based one-dimensional wave propagation models of the arterial system are preferable over one-dimensional wave propagation models in the frequency domain since the latter neglect the non-linear convection forces present in the physiological situation, especially when the vessel is tapered. Moreover, one-dimensional wave propagation models of the arterial system can be used to provide boundary conditions for fully three-dimensional fluid-structure interaction computations that are usually defined in the time domain. In this study, a time-domain-based one-dimensional wave propagation model in a cross-sectional area, flow and pressure (A,q,p)-formulation is developed. Using this formulation, a constitutive law that includes viscoelasticity based on the mechanical behaviour of a Kelvin body, is introduced. The resulting pressure and flow waves travelling through a straight and tapered vessel are compared to experimental data obtained from measurements in an in vitro setup. The model presented shows to be well suited to predict wave propagation through these straight and tapered vessels with viscoelastic wall properties and hereto can serve as a time-domain-based method to model wave propagation in the human arterial system.     

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.     

Experimental study of finite amplitude pulsatile pressure waves in elastic tubes

An experimental study of the propagation of pulsatile pressure waves in an elastic tube was made and results were compared to a theoretical analysis by Lou. The pressure waves were sinusoidally varying acting in a horizontal, longitudinally constrained tube containing water. The independent experimental parameters varied were the pressure wave frequency, pressure wave volume per cycle, static tube pressure and steady flow rate. The wave propagation speeds, measured by non-intrusive techniques, were found to be functions of the wave frequency and the phase angles of the wave elements as theoretically predicted by Lou.     

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.     

A thick walled viscoelastic model for the mechanics of arteries

A knowledge of the mechanics of arteries is of importance in the determination of vessel rheological properties and in the studies of blood flow and certain arterial diseases. Most existing arterial models treat only wave motions; however, other types of motion, in particular those associated with flow development and other end effects, occur in the vascular system. Thus, a model is needed which can be applied to a variety of possible types of motion.

An arterial model is described which includes the effects of thick walls, linear viscoelasticity, and wall tethering. The forms of the displacements and stresses are found independently of the exact form of the applied fluid stresses; thus, the results are applicable to a range of possible dynamical conditions. Displacements and stress states can then be found from experimental or theoretical knowledge of the blood pressure and flow. The results are applied to flow development and wave propagation regions in the arteries.     

Experiments on steady and oscillatory flows at moderate Reynolds numbers in a quasi-two-dimensional channel with a throat.

The study of steady and unsteady oscillatory static fluid pressures acting on the internal wall of a collapsible tube is essential for investigation of the complicated behavior observed when a flow is conveyed inside a tube. To examine the validity of two one-dimensional nonsteady theoretical flow models, this paper presents basic experimental observations of flow separation and reattachment and measured data on the static pressure distributions of the flow in a quasi-two-dimensional channel with a throat, together with information on the corresponding shape of the wall deflection and motion. For combinations of moderate Reynolds numbers and angles of the divergent segment of the channel, a smooth flow is separated from the wall downstream of the minimum cross section and reattached to the wall farther downstream. The measured data are compared with numerical results calculated by the two flow models.     

Transmission characteristics of axial waves in blood vessels

The elastic behavior of blood vessels can be quantitatively examined by measuring the propagation characteristics of waves transmitted by them. In addition, specific information regarding the viscoelastic properties of the vessel wall can be deduced by comparing the observed wave transmission data with theoretical predictions. The relevance of these deductions is directly dependent on the validity of the mathematical model for the mechanical behavior of blood vessels used in the theoretical analysis. Previous experimental investigations of waves in blood vessels have been restricted to pressure waves even though theoretical studies predict three types of waves with distinctly different transmission characteristics. These waves can be distinguished by the dominant displacement component of the vessel wall and are accordingly referred to as radial, axial and circumferential waves. The radial waves are also referred to as pressure waves since they exhibit pronounced pressure fluctuations. For a thorough evaluation of the mathematical models used in the analysis it is necessary to measure also the dispersion and attenuation of the axial and circumferential (torsion) waves.

To this end a method has been developed to determine the phase velocities and damping of sinusoidal axial waves in the carotid artery of anesthetized dogs with the aid of an electro-optical tracking system. For frequencies between 25 and 150 Hz the speed of the axial waves was between 20 and 40 m/sec and generally increased with frequency, while the natural pressure wave travelled at a speed of about 10 m/sec. On the basis of an isotropic wall model the axial wave speed should however be approximately 5 times higher than the pressure wave speed. This discrepancy can be interpreted as an indication for an anisotropic behavior of the carotid wall. The carotid artery appears to be more elastic in the axial than in the circumferential direction.     

Model of geometrical and smooth muscle tone adaptation of carotid artery subject to step change in pressure

Recent experimental studies have shown significant alterations of the vascular smooth muscle (VSM) tone when an artery is subjected to an elevation in pressure. Therefore, the VSM participates in the adaptation process not only by means of its synthetic activity (fibronectins and collagen) or proliferative activity (hypertrophy and hyperplasia) but also by adjusting its contractile properties and its tone level. In previous theoretical models describing the time evolution of the arterial wall adaptation in response to induced hypertension, the contribution of VSM tone has been neglected. In this study, we propose a new biomechanical model for the wall adaptation to induced hypertension, including changes in VSM tone. On the basis of Hill's model, total circumferential stress is separated into its passive and active components, the active part being the stress developed by the VSM. Adaptation rate equations describe the geometrical adaptation (wall thickening) and the adaptation of active stress (VSM tone). The evolution curves that are derived from the theoretical model fit well the experimental data describing the adaptation of the rat common carotid subjected to a step increase in pressure. This leads to the identification of the model parameters and time constants by characterizing the rapidity of the adaptation processes. The agreement between the results of this simple theoretical model and the experimental data suggests that the theoretical approach used here may appropriately account for the biomechanics underlying the arterial wall adaptation.     

General tube law for collapsible thin and thick-wall tubes

Modeling the complex deformations of cylindrical tubes under external pressure is of interest in engineering and physiological applications. The highly non-linear post-buckling behavior of cross-section of the tube during collapse attracted researchers for years. Major efforts were concentrated on studying the behavior of thin-wall tubes. Unfortunately, the knowledge on post-buckling of thick-wall tubes is still incomplete, although many experimental and several theoretical studies have been performed. In this study we systematically studied the effect of the wall thickness on post-buckling behavior of the tube. For this purpose, we utilized a computational model for evaluation of the real geometry of the deformed cross-sectional area due to negative transmural (internal minus external) pressure. We also developed an experimental method to validate the computational results. Based on the computed cross-sections of tubes with different wall thicknesses, we developed a general tube law that accounts for thin or thick wall tubes and fits the numerical data of computed cross-sectional areas versus transmural pressures.     

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.     

Pulse wave propagation in a model human arterial network: assessment of 1-D numerical simulations against in vitro measurements

A numerical model based on the nonlinear, one-dimensional (1-D) equations of pressure and flow wave propagation in conduit arteries is tested against a well-defined experimental 1:1 replica of the human arterial tree. The tree consists of 37 silicone branches representing the largest central systemic arteries in the human, including the aorta, carotid arteries and arteries that perfuse the upper and lower limbs and the main abdominal organs. The set-up is mounted horizontally and connected to a pulsatile pump delivering a periodic output similar to the aortic flow. Terminal branches end in simple resistance models, consisting of stiff capillary tubes leading to an overflow reservoir that reflects a constant venous pressure. The parameters required by the numerical algorithm are directly measured in the in vitro set-up and no data fitting is involved. Comparison of experimental and numerical pressure and flow waveforms shows the ability of the 1-D time-domain formulation to capture the main features of pulse wave propagation measured throughout the system test. As a consequence of the simple resistive boundary conditions used to reduce the uncertainty of the parameters involved in the simulation, the experimental set-up generates waveforms at terminal branches with additional non-physiological oscillations. The frequencies of these oscillations are well captured by the 1-D model, even though amplitudes are overestimated. Adding energy losses in bifurcations and including fluid inertia and compliance to the purely resistive terminal models does not reduce the underdamped effect, suggesting that wall visco-elasticity might play an important role in the experimental results. Nevertheless, average relative root-mean-square errors between simulations and experimental waveforms are smaller than 4% for pressure and 19% for the flow at all 70 locations studied.     

Model for calcium dependent oscillatory growth in pollen tubes

Experiments have shown that pollen tubes grow in an oscillatory mode, the mechanism of which is poorly understood. We propose a theoretical growth model of pollen tubes exhibiting such oscillatory behaviour. The pollen tube and the surrounding medium are represented by two immiscible fluids separated by an interface. The physical variables are pressure, surface tension, density and viscosity, which depend on relevant biological quantities, namely calcium concentration and thickness of the cell wall. The essential features generally believed to control oscillating growth are included in the model, namely a turgor pressure, a viscous cell wall which yields under pressure, stretch-activated calcium channels which transport calcium ions into the cytoplasm and an exocytosis rate dependent on the cytosolic calcium concentration in the apex of the cell. We find that a calcium dependent vesicle recycling mechanism is necessary to obtain an oscillating growth rate in our model. We study the variation in the frequency of the growth rate by changing the extracellular calcium concentration and the density of ion channels in the membrane. We compare the predictions of our model with experimental data on the frequency of oscillation versus growth speed, calcium concentration and density of calcium channels.     

Steady flow and wall compression in stenotic arteries: a three-dimensional thick-wall model with fluid-wall interactions.

Severe stenosis may cause critical flow and wall mechanical conditions related to artery fatigue, artery compression, and plaque rupture, which leads directly to heart attack and stroke. The exact mechanism involved is not well understood. In this paper a nonlinear three-dimensional thick-wall model with fluid-wall interactions is introduced to simulate blood flow in carotid arteries with stenosis and to quantify physiological conditions under which wall compression or even collapse may occur. The mechanical properties of the tube wall were selected to match a thick-wall stenosis model made of PVA hydrogel. The experimentally measured nonlinear stress-strain relationship is implemented in the computational model using an incremental linear elasticity approach. The Navier-Stokes equations are used for the fluid model. An incremental boundary iteration method is used to handle the fluid-wall interactions. Our results indicate that severe stenosis causes considerable compressive stress in the tube wall and critical flow conditions such as negative pressure, high shear stress, and flow separation which may be related to artery compression, plaque cap rupture, platelet activation, and thrombus formation. The stress distribution has a very localized pattern and both maximum tensile stress (five times higher than normal average stress) and maximum compressive stress occur inside the stenotic section. Wall deformation, flow rates, and true severities of the stenosis under different pressure conditions are calculated and compared with experimental measurements and reasonable agreement is found.     

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.     

Evolution of acoustically vaporized microdroplets in gas embolotherapy

Acoustic vaporization dynamics of a superheated dodecafluoropentane (DDFP) microdroplet inside a microtube and the resulting bubble evolution is investigated in the present work. This work is motivated by a developmental gas embolotherapy technique that is intended to treat cancers by infarcting tumors using gas bubbles. A combined theoretical and computational approach is utilized and compared with the experiments to understand the evolution process and to estimate the resulting stress distribution associated with vaporization event. The transient bubble growth is first studied by ultra-high speed imaging and then theoretical and computational modeling is used to predict the entire bubble evolution process. The evolution process consists of three regimes: an initial linear rapid spherical growth followed by a linear compressed oval shaped growth and finally a slow asymptotic nonlinear spherical bubble growth. Although the droplets are small compared to the tube diameter, the bubble evolution is influenced by the tube wall. The final bubble radius is found to scale linearly with the initial droplet radius and is approximately five times the initial droplet radius. A short pressure pulse with amplitude almost twice as that of ambient conditions is observed. The width of this pressure pulse increases with increasing droplet size whereas the amplitude is weakly dependent. Although the rise in shear stress along the tube wall is found to be under peak physiological limits, the shear stress amplitude is found to be more prominently influenced by the initial droplet size. The role of viscous dissipation along the tube wall and ambient bulk fluid pressure is found to be significant in bubble evolution dynamics.