Influence of the distensibility of large arteries on the longitudinal impedance: application for the development of non-invasive techniques to the diagnosis of arterial diseases
© Sahtout and Ben Salah; licensee BioMed Central Ltd. 2012
Received: 27 September 2011
Accepted: 15 March 2012
Published: 16 April 2012
This study shows that the arterial longitudinal impedance constitutes a hemodynamic parameter of interest for performance characterization of large arteries in normal condition as well as in pathological situations. For this purpose, we solved the Navier–Stokes equations for an incompressible flow using the finite element analysis method and the Arbitrary Lagrangian Eulerian (ALE) formulation. The mathematical model assumes a two-dimensional flow and takes into account the nonlinear terms in the equations of fluid motion that express the convective acceleration, as well as the nonlinear deformation of the arterial wall. Several numerical simulations of the blood flow in large vessels have been performed to study the propagation along an arterial vessel of a pressure gradient pulse and a rate flow pulse. These simulations include various deformations of the wall artery leading to parietal displacements ranging from 0 (rigid wall) to 15% (very elastic wall) in order to consider physiological and pathological cases.
The results show significant changes of the rate flow and the pressure gradient wave as a function of aosc, the relative variation in the radius of the artery over a cardiac cycle. These changes are notable beyond a critical value of aosc equal to 0.05. This critical value is also found in the evolution of the longitudinal impedance. So, above a variation of radius of 5%, the convective acceleration, created by the fluid-wall interactions, have an influence on the flow detectable on the longitudinal impedance.
The interpretation of the evolution of the longitudinal impedance shows that it could be a mean to test the performance of large arteries and can contribute to the diagnosis of parietal lesions of large arteries. For a blood vessel with a wall displacement higher than 5% similar to those of large arteries like the aorta, the longitudinal impedance is substantially greater than that obtained in the absence of wall displacement. This study also explains the effects of convective acceleration, on the shape of the decline of the pressure gradient wave and shows that they should not be neglected when the variation in radius is greater than 5%.
The elasticity of large arteries which represent the principal arterial function made that they constitute a very distensible reservoir. Their role is to transform pulsatile flow at the outlet of the heart in a continuous flow in the capillary bed. Consequently, the loss of distensibility is considered a cardiovascular risk factor . This increased rigidity is due to age [2–4] or to various diseases such as hypertension or atherosclerosis [5–8]. Several methods have been proposed for measuring the elasticity of arteries, these methods are based on several techniques such that tonometry, mechanical transducers, echotracking, ultrasonic Doppler, functional Magnetic Resonance Imaging (MRI) and photoplethysmography . However, these methodologies have limitations because some of them may apply only to superficial arteries. Others, such that ultrasonic Doppler and functional MRI, where the site of interest is the deep arteries, are mainly based on the linear theory developed by Womersley [10, 11], which is not applicable in the case of large arteries like the aorta. Indeed, the wall displacements of the large arteries are of finite amplitude and thus the convective acceleration (advection terms) introduce important nonlinear effects [12–14]. In normal physiological situations, the change in radius of these arteries is greater than 10% of the diastolic radius; this is due to the large number of elastin fibers, highly elastic, despite the presence of collagen fibers in the arterial wall. The collagen fibers having an elastic modulus of 3 to 100 times greater than that of elastin fibers are solicited from a pressure of 120mmHg [15, 16]. Beyond that pressure, and especially from 140mmHg (hypertension), we observe a hardening of the arteries and a nonlinear mechanical behavior of the wall.
In this study we are studying the effect of large deformations induced by the mechanical behavior of large arteries (radius of about 1 cm) on the longitudinal impedance. The longitudinal impedance would be an index of clinical interest to quantify the performance of local arterial function especially if it is measured using the techniques of Doppler and MRI. For this purpose, we have solved numerically the equation of Navier–Stokes using Finite Element Analysis (FEA), taking into account the fluid - structure interactions expressed by the convective acceleration terms. To be in a real situation, i.e. in physiologically conditions normal and pathological, we determined the hemodynamic quantities in the arteries where the radius varies from 0% (rigid arterial wall) up to 15% (arterial wall very elastic) of the diastolic radius.
Where W (r, z, t) and U (r, z, t) are the components of blood velocity in the longitudinal and radial directions and P (r, z, t) is the pressure. As we can see, equations (3) and (4) have nonlinear terms, corresponding to the convective accelerations. We have solved numerically these equations taking into account the nonlinear terms in order to study the influence of the distensibility on the flow. The nonlinear behavior of the arterial wall will be taken into consideration in the boundary conditions as described thereafter.
R (z, t) is the inner radius of the vessel, which is a function of z and t because of the radial deformability. The first condition expresses that longitudinal movements of the arterial wall are neglected, which is largely justified in the work of Carew  and Patel . The other boundary conditions at the wall (r = R) and at the center (r = 0) assume an axisymmetric flow and no slip on the wall.
The simulation was performed using the finite element method and a formulation specifically Arbitrary Lagrangian Eulerian (ALE) , to take into account the fluid - structure interactions. In the fluid domain, a mixed formulation velocity - pressure has been implemented. ALE-methods are frequently used to model systems where the physical domain changes with respect to time. Common examples can be found in the field of fluid–structure interactions where the domain movement is due to the force that the fluid exerts on a solid object. These systems have gained a lot of interest to describe blood flows in hemodynamic.
The longitudinal impedance is calculated from the ratios between the amplitudes of pressure gradient and the rate flow Q obtained after analyze of Fourier. The expressions of the mean component, of the fundamental and second harmonic components of the longitudinal impedance are respectively:, and .. We like to point out here that there is no linear relationship between the components of the pressure gradient and those of the flow since these quantities are the solutions of the nonlinear equations governing the flow. However, the Fourier analysis is quite justified by the fact that pressure gradient and rate flow are periodic signals.
The Fourier analysis of various pressure pulses obtained at the inlet of the blood vessel, shows that its amplitude decreases with aosc (Figure 2(a)). This indicates the role of damping played by the nonlinear elasticity of the wall. This last result is similar to the experimental results of Armentano . With regard to the wave of rate flow and of the pressure gradient at the inlet and at distance (10cm, 20cm) we found a distortion of their shape that becomes more pronounced when the distensibility of artery increases (see Figures 2(b), 4(a) and (b)). Fourier analysis shows the presence of a second harmonic whose amplitude increases with aosc (Figures 3(a) and 5), despite that the wave of the pressure at the inlet was imposed sinusoidal. The appearance of this second harmonic reflects the fact that the strong coupling fluid–structure causes significant convective acceleration. These quantities give us an answer on the shape of the pressure and of the flow wave-contour in very elastic arteries like the aorta.
The Evolution of the fundamental amplitude of longitudinal pressure gradient, (Figure 5(b)) for aosc > 0.05 is inconsistent with the conventional results of the fluid mechanics, while that for aosc < 0.05 they are in accordance. The high value of the fundamental amplitude of the pressure gradient obtained for aosc ≥ 0.10 shows the influence of the convective accelerations on the pulsatile flow. We observe two regimes of flow, one, depending on the radius and the other dominated by the elastic behavior, these two behaviors being defined by a critical value of aoscc = 0.05. This critical value would correspond to the limit of the radius variation where the assumption of linearization of Navier - Stokes equations can be applied.
The decreasing part of the mean pressure gradient (Figure 4(b)) can be explained only by the fact that the convective acceleration have a significant effect on the flow of large arteries and therefore cannot be neglected for the large strains (aosc ≥ 0.10). For aosc ≥ 0.10 we observe a reverse pressure gradient. This suggests that the convective effects caused by the elasticity of large arteries play an important role on the hemodynamic of this type of artery by reducing the wall arterial stress. This finding would result in a reduction of energy dissipation and preservation of endothelial cells [23, 24]. In summary, we can say that the influence of the convective acceleration developed on the flow is felt locally (mean pressure gradient) and at distance (pulsatile gradient pressure).
In the light of the above remarks, it seems that the nonlinear model preserves the pulse shape much better than the linear propagation mathematical model. The shape of the pressure gradient pulses, obtained from the aosc values equal to 0.10 and 0.15, reproduces with a certain similarity the pressure gradient pulses encountered in large elastic arteries like the aorta. On these physiological curves, one finds the same type of diastolic decrease (except the incisura caused by the closure of the aortic valve) which is traditionally explained by the superposition of the incident wave and reflected wave. However, several experimental studies have been designed to test the accuracy of the theoretical formulations of propagation in describing the physiological phenomenon [12–14]. These experiments have shown that many hemodynamic quantities computed by the linear propagation mathematical models are very different from the corresponding measured values and could not explain the deviation between observations and theory. To get a better response with the linear model, Reuderink and al  introduce the wall viscoelasticity and conclude that wall viscoelasticity is able to compensate the lack of convective accelerations in the linear model quite well. In other studies  this is the tapering of the tube that is introduced; however, this parameter is not involved in the large arteries. It seems that the nonlinear terms of the Navier–Stokes equations increase with the deformation, meaning that the elasticity of the wall also plays a role in the shape of the pulse decay by making it less steep.
The physical interpretation of the changes in the pressure gradient according to the variation in radius can be performed only when is reported to the rate flow, hence the calculation of the longitudinal impedance. We have therefore calculated the amplitudes of longitudinal impedance from the values of longitudinal pressure gradient and rate flow for each simulation. The fundamental amplitude of the longitudinal impedance Figure 6(a) has also a minimum at aosc = 0.05. The decreases of the impedance for aosc ≤ 0.05 reflects the fact that the flow becomes less resistant due to the increase in radius during the cardiac cycle, whereas, for aosc > 0.05, it tends to increase and to become more important than the rigid case (aosc = 0). The results found for aosc ≤ 0.05 are consistent with those observed in medium-sized arteries, which main role is to perfuse the blood downstream (see additional file 1). This last result is also encountered in large arteries in some pathology leading to a hardening of the arteries such as hypertension or atherosclerosis. We can explain the increase of the fundamental longitudinal impedance, for values of aosc ≥ 0.05 by the fact that a part of the volume of fluid remained stored inside the blood vessel (see additional file 2 and 3). We have seen, above, the importance of convective effects for these values of aosc, which also contributes to mass transport. Normally, the increases of the fundamental pressure gradient amplitude would have led to a significant increase in the fundamental rate flow amplitude. However, we observe the opposite. In fact, we see that the pulsatile rate flow diminishes with aosc at a distance of 20 cm, while it was practically constant at the entrance. These results are very remarkable for values of aosc equal to 0.10 and 0.15, encountered in very elastic arteries. Thus, knowledge of the longitudinal impedance allows us to assess the ability of the artery to store a certain volume of blood.
As shown in Figure 6(b) the mean amplitude of the longitudinal impedance is much smaller than the pulsatile component (fundamental and second harmonic). In fact, different studies [27, 28] on the "sleeve effect" and on the hemodynamic effects of vein grafts show that it is the pulsatile components that are most useful to quantify the performance of the local arterial function. However, their studies are based on Womersley's solutions applicable only in small arteries where the parietal displacements (aosc) are well below 0.05. Our focus being to study the role of wall elasticity on the local hemodynamic and chiefly at the level of large arteries, we have not, in our simulation, used restrictive assumptions on the convective acceleration terms which are the nonlinear terms of the Navier - Stokes. In [27, 28] the authors have also shown that the longitudinal impedance, increases with the frequency, we observe the same trend for the second harmonic at z = 0.10m and aosc < 0.05, see Figure 5(a). Nevertheless, for the other cases (z = 0.10m–aosc > 0.07 and z = 0.20m–0 ≤ aosc ≤ 0.15) our simulation shows that the second harmonic of the longitudinal impedance is lower than the fundamental and decreases when aosc increases. This is due to the fact that the rate flow and the pressure gradient vary in the same proportion and thus, the convective effects are less important at frequencies above the fundamental. Although the artery longitudinal impedance is usually used as intermediate parameter in calculations, our study has shown a correlation between the longitudinal impedance and the mechanical properties of the wall artery, meaning that it can constitute a clinical important clue for the screening of diseases caused by a change in the mechanical behavior of arteries. Indeed, several studies have shown the incidence of arterial stiffness on cardiovascular disease. Our previous work [29, 30] have shown that the pressure gradient and rate flow can be determined from the rate flow and therefore from velocimetry data at the center of blood vessel, thus, it would be interesting to access this index from non-invasive techniques like ultrasonic Doppler or MRI.
In our numerical model we have been subject to certain limitations as the frequency range and reflections. We do not consider multiple reflections caused by the structure of the arterial system because we want to show the response of the elasticity of the arteries on the flow [13, 31, 32]. Indeed, previous work [30, 33–35] shows that the extraction of the incident wave would assess the performance of the heart and arterial system. Regarding frequency, the frequency of 1 Hz is often used in the literature because it represents the fundamental frequency of blood pressure that determines the shape and amplitude of the waveform flow.
The main results concern the effects of the elasticity and the loss of elasticity of the arterial wall on the flow waveform and the hemodynamic.
The study of the waveform is useful for the hemodynamic signal processing. Our results show that the variation in the radius of the artery during a cardiac cycle effects the hemodynamic. We observe two types of behavior defined by a critical point corresponding to a variation of the radius equal to 5% (aosc). This critical point represents the limit of deformation of the arterial wall beyond which the convective acceleration terms play a significant role on the flow.
The longitudinal impedance was calculated to interpret the effects of convective acceleration on the flow. The results show that for aosc > 0.05 the convective acceleration decelerate the flow so that a portion of blood volume is retained in the upstream of the artery. The longitudinal impedance would, therefore, be a clinically useful parameter to assess the elasticity of the arteries, in order to detect and localize vascular diseases affecting the wall of the elastic arteries like the aorta. Thus, determination of the longitudinal impedance, for a functional exploration of the arterial system, could be performed by non-invasive techniques such as ultrasonic Doppler velocimetry.
- Safar M: Paroi artérielle et vieillissement vasculaire. Éditions scientifiques et médicales Elsevier SAS. 2002, PpGoogle Scholar
- Learoyd BM, Taylor MG: Alterations with age in the viscoelastic properties of human arterial walls. Circ Res. 1966, 18: 278-292. 10.1161/01.RES.18.3.278.View ArticleGoogle Scholar
- Isnard RN, Pannier BM, Laurent S, London GM, Diebold B, Safar ME: Pulsatile diameter and elastic modulus of the aortic arch in essential hypertension: a non-invasive study. Am Coll Cardiol. 1989, 13: 399-405. 10.1016/0735-1097(89)90518-4.View ArticleGoogle Scholar
- Lang RM, Cholley BP, Korcarz C, Marcus RH, Shroff SG: Measurements of regional elastic properties of the human aorta: a new application of transesophageal echocardiography with automated border detection and calibrated subclavian pulse tracings. Circulation. 1994, 90: 1875-1882. 10.1161/01.CIR.90.4.1875.View ArticleGoogle Scholar
- Farrar DJ, Bond MG, Riley WA, Sawyer JK: Anatomic correlates of aortic pulse wave velocity and carotid artery elasticity during atherosclerosis progression and regression in monkeys. Circulation. 1991, 83: 1754-1763. 10.1161/01.CIR.83.5.1754.View ArticleGoogle Scholar
- Ting CT, Brin KP, Lin SJ, Wang SP, Chang MS, Chiang BN, Yin FCP: Arterial hemodynamic in human hypertension. J Clin Invest. 1986, 78: 1462-1471. 10.1172/JCI112737.View ArticleGoogle Scholar
- Liu Z, Ting CT, Zhu S, Yin FCP: Aortic compliance in human hypertension. Hypertension. 1989, 14: 129-136. 10.1161/01.HYP.14.2.129.View ArticleGoogle Scholar
- O'Rourke MF: Arterial stiffness, systolic blood pressure, and logical treatment of arterial hypertension. Hypertension. 1990, 15: 339-347. 10.1161/01.HYP.15.4.339.View ArticleGoogle Scholar
- Laurent S, Cockroft J, Van Bortel L, Boutouryie P, Giannattasio C, Hayoz D, Pannier B, Vlachopoulos C, Wilkinson I, Struijker-Boudier H: Expert consensus document on arterial stiffness: methodological issues and clinical applications. Eur Heart J. 2006, 27: 2588-2605. 10.1093/eurheartj/ehl254.View ArticleGoogle Scholar
- Womersley JR: Oscillatory Motion of a Viscous Liquid in a Thin-Walled Elastic Tube–I. The Linear Approximation for Long wave. Phil Mag. 1955, 46: 199-221.View ArticleMathSciNetMATHGoogle Scholar
- Womersley JR: Oscillatory Flow in Arteries: The Constrained elastic tube as a Model of Arterial Wall and Pulse Transmission. Phys Med. 1957, 46: 178-187.View ArticleGoogle Scholar
- Milnor WR, Nichols WW: A new method of measuring propagation coefficients and characteristic impedance in blood vessels. Circ Res. 1975, 36: 631-639. 10.1161/01.RES.36.5.631.View ArticleGoogle Scholar
- Milnor WR, Bertram CD: The relation between arterial viscoelasticity and wave propagation in the canine femoral artery in vivo. Circ Res. 1978, 43: 870-879. 10.1161/01.RES.43.6.870.View ArticleGoogle Scholar
- Li JK, Melbin J, Rifflin RA, Noordergraaf A: Pulse wave propagation. Circ Res. 1981, 49: 442-452. 10.1161/01.RES.49.2.442.View ArticleGoogle Scholar
- Cox RH: Passive mechanics and connective tissue composition of canine arteries. Am J Physiol. 1978, 234: 533-541.Google Scholar
- Amentano RL: Determination in vivo arterial hemodynamic characteristics, application to hypertension. PhD. Thesis. 1999, Paris 7,Google Scholar
- Coussot P, Grossiord JL: Understand the rheology of the circulation of blood to the taking of concrete. 2002, EDP Sciences, groupe Français de RhéologieGoogle Scholar
- Flaud P, Geiger D, Oddou C, Quemada D: Pulsed flow in viscoelastic pipe. Application to the study of blood flow. J Phys. 1974, 35: 869-882. 10.1051/jphys:019740035011086900.View ArticleGoogle Scholar
- Carew TE, Vaishnav RN, Patel DJ: Compressibility of the arterial wall. Circ Res. 1968, 23: 61-68. 10.1161/01.RES.23.1.61.View ArticleGoogle Scholar
- Patel DDJ, Greenfield JC, Fry DL: In vivo pressure –Length- radius relationship of certain blood vessels in man and dog. Pulsatile Blood Flow. Edited by: Attinger EO. 1964, McGraw-Hill, New York, 293-302.Google Scholar
- Cox RH: Passive mechanics and connective tissue composition of canine arteries. Am J Physiol. 1978, 2: 533-541.Google Scholar
- Donea J, Huerta A, Ponthot JP, Rodriguez-Ferran A: Arbitrary Lagrangian - Eulerian Methods. Chapter 14. 2004, John Wiley & Sons, New YorkGoogle Scholar
- Caro CG, FitzGerald JM, Schroter RC: Atheroma and arterial wall shear observations, correlation and proposal of a shear dependant mass transfer mechanism for artherogenesis. Proc. Royal Soc. London Ser. B. 1971, 177: 109-111. 10.1098/rspb.1971.0019.View ArticleADSGoogle Scholar
- Ku DN, Giddens DP, Zarins CK, Galgov S: Pulsatile flow and atherosclerosis in the human carotid bifurcation. Arteriosclerosis. 1985, 5: 293-302. 10.1161/01.ATV.5.3.293.View ArticleGoogle Scholar
- Reuderink P, Hoogstraten HW, Sipkma P, Hillen B, Westerhof N: Linear and nonlinear one-dimensional models of pulse wavetransmission at high Womersley numbers. J Biomech. 1989, 22: 819-827. 10.1016/0021-9290(89)90065-1.View ArticleGoogle Scholar
- Belardinelli E, Cavalcanti S: Anew nonlinear two-dimensional model of blood motion in tapered and elastic vessels. Comput Biol Med. 1991, 21: 1-3. 10.1016/0010-4825(91)90030-D.View ArticleGoogle Scholar
- Jager GN, Westerhof N, Noordergaaf A: Oscillatory flow Impedance in electrical analogy of arterial system: Representation of sleeve effect and non-Newtonian properties of flow. Circ Res. 1965, 16: 121-133. 10.1161/01.RES.16.2.121.View ArticleGoogle Scholar
- Skelly CL, Meyerson S: Hemodynamics of venous grafts: Evaluation and meaning. Ann Chir Vasc. 2001, 15: 110-122.View ArticleGoogle Scholar
- Sahtout W, Ben Salah R: Influence of convective effects on hemodynamics of large arteries in normal and pathological conditions: Influence of convective effects on hemodynamics of large arteries in normal and pathological conditions. Rev Rhéol. 2010, 18: 1-10.Google Scholar
- Abdessalem KB, Sahtout W: Numerical simulation of non-invasive determination of the propagation coefficient in arterial system using two measurements sites. Eur Phys J Appl. 2007, 40: 211-219. 10.1051/epjap:2007151.View ArticleADSGoogle Scholar
- Newman DL, Greendale SE: Analysis of forward and backward pressure waves by a total occlusion method. Med Bioeng Comp. 1980, 80: 240-245.Google Scholar
- Campbell KB, Lewis CL, et al: Pulse reflection sites and effective length of the arterial system. Am J Physiol. 1989, 256: 1684-1689.Google Scholar
- Rogova I: Propagation of waves in hemodynamics blood: Application to indirect assessment of pathophysiological parameters. 1998, PhD thesis, Paris 7Google Scholar
- Khir AW, O’Brien A, et al: Determination of wave speed and wave separation in the arteries. J Biomech. 2001, 34: 1145-1155. 10.1016/S0021-9290(01)00076-8.View ArticleGoogle Scholar
- David S, Berger JK, Noordergraaf A: Differential effects of wave reflections and peripheral resistance on aortic blood pressure: a model-based study. Am J Physiol Soc. 1994, 94: 1626-1642.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.