PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of jbiolphysspringer.comThis journalToc AlertsSubmit OnlineOpen Choice
 
J Biol Phys. Jan 2011; 37(1): 51–60.
Published online Aug 17, 2010. doi:  10.1007/s10867-010-9196-2
PMCID: PMC3006462
The design of a reflectionless arterial prosthesis
Guy Richard Kolcorresponding author and Paul Woafo
Laboratory of Modelling and Simulation in Engineering and Biological Physics, P.O. Box 812, Yaoundé, Cameroon
Guy Richard Kol, kolguy_r/at/yahoo.fr.
corresponding authorCorresponding author.
Received April 17, 2010; Accepted July 15, 2010.
Abstract
We propose a new technique to characterize a reflectionless arterial prosthesis. The corresponding transmission and reflection coefficients are determined from the geometric and the elastic properties of the arterial wall, and the interaction between the latter and the prosthesis are studied accordingly.
Keywords: Blood waves, Transmission coefficient, Reflection coefficient, Prosthesis
Elastic tubes are widely used nowadays in fluid flow, and their importance no longer needs to be demonstrated. From the biological field up to the industrial domain, different types of elastic tubes are encountered and while many studies have been done, there still exist unresolved questions.
The first studies on elastic tubes presented a linear approach of the fluid flow therein with differences arising only in the consideration of the tube wall contribution [1, 2]. The nonlinear studies that came after brought out the importance of the fluid convective term and particularly took into account the dynamical properties of the tube wall where the nonlinear stress–strain relationship is considered [36]. Most of the works that followed contributed to improve the modeling, in a more realistic way, of elastic tubes [79] and took into account viscosity, tapering aspect of the tubes [5, 10], pathological cases of vessels and prostheses [11].
Some studies on blood waves focused on the study of the behaviour of waves in the regions of discontinuities [12, 13]. In general, any discontinuity zone present in an artery gives rise to the phenomena of wave transmission and reflection. These discontinuities are in general due to the change in the elastic properties of the arterial wall, diseases, bifurcations and prostheses.
The objective of this paper is to determine the expression of the transmission and reflection coefficients in the prosthesis using a nonlinear model of blood flow. This theory is based on the separation of pressure waves into their forward and backward running components proximal and distal to the prosthesis. We then use the same approach to investigate the characteristics of a reflectionless prosthesis.
In the study, we consider that the fluid in the vessel is Newtonian, viscous, homogeneous and incompressible. Furthermore, we take into account the variation of the radius and elasticity of the wall. The rest of the paper is therefore organized as follows: in Section 2, the physical and mathematical problems are formulated. In Section 3, the nonlinear model of the flow is used to determine the transmission and reflection coefficients of the prosthesis, and the conditions for a reflectionless prosthesis are derived accordingly. The last section is devoted to concluding remarks.
Physical model and governing equations
The artery is assumed to be conical, with a nonlinear stress–strain relation. The geometry of the artery with the presence of a prosthesis is illustrated in Fig. 1. The index i (i = 1,2 or 3) represents each portion, where i = 2 stands for the prosthesis and i = 1 and i = 3 represent the leftmost and rightmost portions of the vessel, respectively. For the present discrete case, we apply the averaging procedure, used in the literature, to the Navier–Stokes equations [5] and we derive the following one-dimensional model
equation M1
1
equation M2
2
where z is the axial coordinate, t is time, Pi is the pressure inside the vessel, Wi is the axial component of the fluid velocity, Ri is the radius of the vessel, ρ is the fluid density and νi is kinematic viscosity.
Fig. 1
Fig. 1
Geometry of the artery with prosthesis, where a stented vessel is shown in a cylindrical coordinate system (r, radial coordinate; z, longitudinal coordinate)
For the wall dynamics, we use the second law of Newton on a portion of the vessel wall and we obtain the relation [14]:
equation M3
3
where equation M4 is the pressure outside the vessel, δRi denotes the vessel deformation, ρ0i is the wall density, and Hi and hi are the effective inertial thickness and the thickness of the wall, respectively. γi is the viscoelasticity coefficient due to the absorption effects of the material. βi is the shear modulus, and [var phi]i is another coefficient of viscoelasticity. σti is the approximated exponential function for the stress–strain relationship given in [15, 16]. This approximated relation is widely known as [3, 4, 6, 7]
equation M5
4
where equation M6 is the stationary radius of the vessel, Ei is the Young modulus and ai represents the nonlinear coefficient of elasticity in the portion i .
The vessel radius Ri(z,t) in the portion i is described as
equation M7
5
As vessel segments are tapered elastic tubes with increasing rigidity away from the heart, they represent general cases of thin-walled elastic tubes with variable Young modulus. Many authors [4, 7, 8] showed that the stationary radius and the variation of elasticity along the wall will obey the relations
equation M8
6
with f (z) =  mi(z  z0i), g(z) = ni(z  z0i). Here, f (z) denotes the decrease of the tube radius and g(z) the increase of the tube wall rigidity along the tube originating from the heart. R0i is the radius at the entrance of the portion i and z0i the abscissa at the entrance of the vessel segment, in the same portion. The positive coefficients mi and ni characterize the rate of decrease and increase of the radius and wall rigidity in the portion i, respectively.
As we consider weak wall displacement, we introduce the approximation:
equation M9
7
where equation M10 and equation M11 are the normalized cross sectional area of the tube and the cross sectional area at the entrance of the ith portions, respectively. R0, w0 and E0 are, respectively, the reference values of radius, velocity and Young modulus at the entrance of the vessel. We define the viscosity coefficients equation M12 and equation M13 as:
equation M14
8
and the dimensionless quantities
equation M15
9
and we also include the following coefficients
equation M16
10
We assume that the first and last portions (i = 1 and i = 3) are made up of the same material properties.
So, Eqs. 13 thus illustrate mass conservation,
equation M17
11
momentum conservation,
equation M18
12
and wall dynamics,
equation M19
13
where
equation M20
with equation M21 being the relative stiffness of the vessel. The flow pulsation has a dimensionless period which is given by the inverse of the Strouhal number equation M22.
The linear theory for separation of forward and backward waves
Following the mathematical procedure of [3, 17], we set
equation M23
14
equation M24
15
equation M25
16
where ϵ is a small parameter measuring the weakness of dispersion and nonlinearity. Also weak are the parameters modelling the taper effect of the wall, the compliance variation and the blood viscosity which become
equation M26
17
Next, the coordinate system is changed to slowly varying coordinates
equation M27
18
where ci represents the wave group velocity and is assumed to be the velocity of the linear wave as it enters the medium. ξi is an independent variable representing the forward wave whereas ηi stands for the backward wave.
The approximation of the flow motion at the first order, ϵ1, gives:
equation M28
19
equation M29
20
equation M30
21
In the present work, the study will be limited to this first order of perturbation. The second order, useful to obtain the analytic expression of the wave, is not used here. So, we consider the waves to be separated explicitly as forward and backward elements (where f and b are, respectively, used as superscripts) [18, 19].
equation M31
22
equation M32
23
equation M33
24
Introducing Eqs. 2224 into Eqs. 1921 leads to the following relation
equation M34
25
and
equation M35
26
with
equation M36
27
Transmission and reflection coefficients
In this part, we propose a new technique for the estimation of the transmission, T, and the reflection, B, coefficients in the prosthesis proportional to the local transmission, Ti, and reflection, Bi, coefficients at the interface of the prosthesis.
For the computation of the characteristic coefficients, we use the impedance of the system given by the ratio of the pressure and flow at the interfaces [16, 17]
equation M37
28
The local transmission, Ti, and reflection, Bi, coefficients are defined as
equation M38
29
and are expressed in terms of the flows. Of course, in the literature, this is not usual, but from the conservation law equation M39 the relation Ti + Bi = 1 is not verified when the pressures are mainly used. This also brings out the difference between the method used in this paper and the one used by Stergiopulos et al. [18].
If we consider at the interfaces the conservation of mass and energy, it will appear that, at the first interface,
equation M40
30
while at the second interface
equation M41
31
Inserting Eqs. 2527 into Eqs. 30 and 31 helps to obtain the local transmission coefficient T1 (resp. T2) and the local reflection coefficient B1 (resp. B2), respectively, at the first and second interfaces
equation M42
32
equation M43
33
with Lv as the length of the left portion of the vessel and Lp as the length of the prosthesis.
In this work, we suppose that the distal forward wave, equation M44, is the sum of two waves: the part of the proximal forward wave that became transmitted through the prosthesis and the part of the distal backward wave that became reflected by the prosthesis. And we obtain a similar equation for the proximal backward running wave, equation M45. These equations are given by
equation M46
34
where T and B are the global transmission and reflection coefficients due to the overall prosthesis. Using Eq. 34, we obtain the expressions of the transmission and reflection coefficients, as
equation M47
35
Inserting relation Eq. 29 into Eq. 35 helps to derive the transmission, T, and reflection, B, coefficients in the prosthesis depending on the local transmission, Ti, and reflection, Bi, coefficients at the interfaces
equation M48
36
equation M49
37
Condition for a reflectionless arterial prosthesis
From Eq. 37, we get the conditions for a reflectionless arterial prostheses. The reflectionless arterial prosthesis is obtained when the reflection coefficient is equal to 0 and the transmission coefficient equal to 1. This gives us the following condition (with B1  0):
equation M50
38
with
equation M51
39
and
equation M52
40
Using this condition (see relation 38), the optimal thickness of the prosthesis is derived as
equation M53
41
The manufacturing of reflectionless arterial prostheses being one of the primary objectives of research in this field, our result gives a new mathematical condition for the reflectionless arterial prosthesis. This condition takes into account the geometric and elastic properties of the vessel wall.
For numerical calculation, we have used the characteristic parameters for the femoral artery of a dog [6]: Lv = 4.5 cm, equation M54 dyn/cm2, n1 = 0.067 cm  1, m1 = 0.080 cm  1 and the thickness of the vessel is 0.018 cm. We consider the prosthesis parameters to be Lp = 4.0 cm, n2 = 0.069 cm  1 and m2 = 0.089 cm  1.
Figure 2 shows that the transmission coefficient increases to the value 1 on the positive side of the ordinate and the reflection coefficient increases to the value 0 on the negative side, as the ratio of the product of the Young modulus and the wall thickness is increased and for R01/R02 = 0.86. This behaviour is comparable with the effects of the severity of an extended stenosis as obtained by Stergiopulos et al. [18]. The negative values obtained in the case of the reflection coefficient are due to the fact that the reflected waves velocity is negative (see Eq. 26).
Fig. 2
Fig. 2
The dependence of the transmission and reflection coefficients on the ratio of the product of Young modulus and wall thickness of the stented vessel for R01/R02 = 0.86
In Fig. 3, the ratio of thickness (Eq. 41) is plotted as the ratio of the radius varies. When the ratio of the radii increases, the thickness ratio decreases. This figure illustrates a reflectionless arterial prosthesis. For example, in the case of a Palmaz stent with Young modulus equation M55 dyn/cm2 [20], inserted in the femoral artery, the ratio of the thickness of the stent to that of the vessel is equation M56. Using this condition, the value of the transmission coefficient is approximately equal to 1 while the reflection coefficient is equal to 0, when the ratio of the radius is R01/R02 = 0.94. Taking into account the thickness of the femoral artery (h01 = 0.018 cm), the thickness of the stent is obtained as h02 = 0.173 cm.
Fig. 3
Fig. 3
Dependence of the ratio of thickness on that of the radius in the case of the reflectionless prosthesis
In Fig. 4, we show that the variation of the length of the left portion of the vessel brings about modifications in the characteristics of the perfect prosthesis. When the length of the vessel increases, the characteristics of the reflectionless prosthesis decrease. As a whole, these results present the interactions between the vessel wall and the prosthesis. Such a behaviour would also be observed for the variation of the elasticity of the vessel. These results are further confirmed by the case of Fig. 5, where the variation of the prosthesis length is taken into account. So, when the length of the prosthesis increases, the characteristics of the reflectionless prosthesis increase. In fact, if there is a perturbation at the entrance of the prothesis, that prothesis, due to its characteristics will overcome the perturbation for a better blood flow. That is why increasing the length of the prothesis deeply modifies the condition for a reflectionless prothesis as depicted in Fig. 5.
Fig. 4
Fig. 4
Influence of the left vessel length on the reflectionless prothesis condition
Fig. 5
Fig. 5
Influence of the prothesis length on the reflectionless prothesis condition
Conclusion
In this study, the concept of wave separation has been used to derive the reflection and transmission coefficients of blood waves created by a prosthesis inserted in an artery. This has led to the derivation of conditions for a reflectionless prosthesis. These conditions are useful for the design of a reflectionless arterial prosthesis. In our methodology, we take explicitly into account the interactions between the vessel wall and the prosthesis at the entrance and the exit of the prosthesis.
1. Fung YC. New York. New York: Springer; 1984. Biodynamics: Circulation.
2. Pedley TJ. Perspectives in Fluid Mechanics. Cambridge: Cambridge University Press; 2000.
3. Yomosa S. Solitary waves in large blood vessels J. Phys. Soc. Jpn. 1987. 56506–520.5208931161987JPSJ...56..506Y. doi: 10.1143/JPSJ.56.506. [Cross Ref]
4. Noubissie S, Woafo P. Dynamics of solitary blood waves in arteries with prostheses. Phys. Rev. E. 2003;67:419111–419118. doi: 10.1103/PhysRevE.67.041911. [PubMed] [Cross Ref]
5. Demiray H. Weakly nonlinear waves in a viscous fluid contained in a viscoelastic tube with variable cross-section Eur. J. Mech. 2005. 24337–347.3471069.740162164364. doi: 10.1016/j.euromechsol.2004.12.002. [Cross Ref]
6. Ntchantcho R, Noubissie S, Woafo P. Numerical simulation of solitary blood waves in an elastic tube subjected to a localised deformation Commun. Nonlinear Sci. Numer. Simulat. 2007. 121572–1583.15831118.3503523326472007CNSNS..12.1572N. doi: 10.1016/j.cnsns.2006.03.009. [Cross Ref]
7. Paquerot JF, Remoissenet M. Dynamics of nonlinear blood pressure waves in large arteries Phys. Lett. A 1994. 19477–82.821994PhLA..194...77P. doi: 10.1016/0375-9601(94)00729-9. [Cross Ref]
8. Noubissie S, Woafo P. Dynamics of solitary waves through taper-thin elastic tube with localised deformation Phys. Scr. 2004. 69249–256.2561067.761052004PhyS...69..249N. doi: 10.1238/Physica.Regular.069a00249. [Cross Ref]
9. Paquerot JF, Lambrakos SG. Monovariable representation of blood flow in large elastic arteries Phys. Rev. E 1994. 493432–3439.34391994PhRvE..49.3432P. doi: 10.1103/PhysRevE.49.3432. [PubMed] [Cross Ref]
10. Duan WS, Wang BR, Wei RJ. The decay of soliton in small blood arteries J. Phys. Soc. Jpn. 1994. 65945–947.9471996JPSJ...65..945D. doi: 10.1143/JPSJ.65.945. [Cross Ref]
11. Tortoriello A, Pedrizzetti G. Flow-tissue interaction with compliance mismatch in a model stented artery. J. Biomech. 2004;37:1–11. doi: 10.1016/S0021-9290(03)00259-8. [PubMed] [Cross Ref]
12. Khir AW, Parker KH. Measurements of wave speed and reflected waves in elastic tubes and bifurcations. J. Biomech. 2002;35:775–783. doi: 10.1016/S0021-9290(02)00025-8. [PubMed] [Cross Ref]
13. Khir AW, Parker KH. Wave intensity in the ascending aorta: effects of arterial occlusion. J. Biomech. 2005;38:647–655. doi: 10.1016/j.jbiomech.2004.05.039. [PubMed] [Cross Ref]
14. Formaggia L, Lamponi D, Quarteroni A. One-dimensional models for blood flow in arteries J. Eng. Math. 2003. 47251–276.2761070.760592038983. doi: 10.1023/B:ENGI.0000007980.01347.29. [Cross Ref]
15. Akio S, Masamitsu H, Yoshihisa U. Pressure pulse wave for blood flow in the aorta from the viewpoint of the nonlinear Toda lattice. Phys. Lett. A. 1996;221:395–399. doi: 10.1016/0375-9601(96)00563-4. [Cross Ref]
16. Xie J, Zhou J, Fung YC. Bending of blood vessel wall: stress–strain laws of the intima-media and adventitial layers. J. Biomech. Eng. 1995;117:136–145. doi: 10.1115/1.2792261. [PubMed] [Cross Ref]
17. Duan WS, Wang BR, Wei RJ. Reflection and transmission blood waves due to arterial branching Phys. Rev. E 1997. 551773–1778.17781997PhRvE..55.1773D. doi: 10.1103/PhysRevE.55.1773. [Cross Ref]
18. Stergiopulos N, Spiridon M, Pythoud F, Meister JJ. On the wave transmission and reflection properties of stenoses. J. Biomech. 1996;29:31–38. doi: 10.1016/0021-9290(95)00023-2. [PubMed] [Cross Ref]
19. Li JK, editor. Dynamics of the vascular system. Singapore: World Scientific; 2004.
20. Berry JL, Manoach E, Mekkaoui C, Rolland PH, Moore JE, Rachev A. Hemodynamics and wall mechanics of a compliance matching stent: in vitro and in vivo analysis. J. Vasc. Interv. Radiol. 2002;13:97–105. doi: 10.1016/S1051-0443(07)60015-3. [PubMed] [Cross Ref]
Articles from Journal of Biological Physics are provided here courtesy of
Springer Science+Business Media B.V.