Jtam-A4.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 51, 2, pp. 409-418, Warsaw 2013 TWO-PHASE FLOW IN A CATHETERIZED ARTERY WITH ATHEROSCLEROSIS Jessica M. Ledesma, Daniel N. Riahi, Ranadhir Roy University of Texas-Pan American, Department of Mathematics, Edinburg, Texas, USA e-mail: driahi@utpa.edu We consider the problem of blood flow in a catheterized artery and in the presence of atherosclerosis,which is chosenbasedontheavailable experimentaldata.Theatherosclerosis is a conditionwhere an arterywall thickens as a result of fattymaterials such as cholesterol. The use of catheter is of immense importance as a standard tool for diagnosis and treatment in a patience whose artery is affected adversely by the presence of atherosclerosiswithin the artery. The blood flow in the arterial tube is represented by a two-phase model composing a suspension of erythrocytes (red cells) in plasma. The coupled differential equations for both fluid (plasma) and particles (red cells) are solved theoretically subjected to reasonable modeling and approximations. The important quantities such as plasma speed, velocity of red cells, blood pressure force, impedance (blood flow resistance) and the wall shear stress are computed for different values of the catheter size, axial location of atherosclerosis and the hematocrit due to the red cells-plasma combination of the blood flow system. Key words: arterial flow, blood flow, impedance, atherosclerosis, pressure, shear stress 1. Introduction Diseases in the blood vessels and in the heart, such as heart attack and stroke, are the major causes ofmortality worldwide. The underlying cause for these events is the formation of lesions, known as atherosclerosis. These lesions and plaques can grow and occlude the artery and hence prevent blood supply to the distal bed. Plaques with calcium in them can also rupture and initiate the formation of blood clots (thrombus). The clots can form as emboli and occlude the smaller vessels that can also result in interruption of blood supply to the distal bed. Plaques for- med in coronary arteries can lead to heart attacks and clots in the cerebral circulation can result in the stroke. There are a number of risk factors for the presence of atherosclerotic lesions. The common sites for the formation anddevelopment of atherosclerosis include the coronary arteries, the branching of the subclavian and common carotids in the aortic arch, the bifurcation of the common carotid to internal and external carotids especially in the carotid sinus region distal to the bifurcation, the renal arterial branching in the descending aorta and in the ileofemoral bifurcations of the descending aorta. The common feature in the location for the development of the lesion is the presence of curvature, branching, and bifurcation present in these sites. The fluid dynamics at these sites can be anticipated to be vastly different from other segments of the arteries that are relatively straight and devoid of any branching segments. Hence, several investigators have attempted to link the fluid dynamically induced stress with the formation of atherosclerotic lesion in the human circulation. By assuming the artery to be circularly cy- lindrical in shape, Mishra and Panda (2005) studied the flow of blood in stenosed artery for the Casson type fluid. Young and Tsai (1973) discussed some characteristics of flow of blood in stented arteries. The blood vessels carry blood from the heart to all the organs and tissues of the body including brain, kidneys, gut, muscles, and the heart itself. Venkateswarlu and Rao (2004) studied an assumed oscillatory form of the blood flow through an indented tube in the 410 J.M. Ledesma et al. presence of steady single stenosis with a very simple shape. They used the so-called Einstein model for the viscosity of the blood but for a variable volume flow rate and the prescribed value for the magnitude of the pressure gradient. Srivastava et al. (2010) studied arterial blood flow through an overlapping stenosis (Mishra and Panda, 2005) by using a Casson type fluid flow. They calculated the impedance and shear stress for different stenosis height. Riahi et al. (2011) investigated arterial blood flow in the presence of an overlapping stenosis using the variable viscosity model due to Einstein for the blood flow. All the investigations described above were for the cases where no catheter was inserted into the artery, although there have also been a number of studies of the blood flow systems in catheterized arteries (Kanai et al., 1970; Back, 1994; Back et al., 1996; Srivastava and Rastogi, 2010). Inall the studies thathavebeencarried out so far and investigated arterial bloodflowsystem, and some of whichwere listed in the previous paragraph, the shape of stenosis was based on the assumed analytical description due to some forms of certain functions. In the present study, we apply for the first time an experimentally based form (Back et al., 1984) for the atherosclerosis shape in the artery (Fig. 1), where the blood is represented by a two-phase macroscopic model. InFig. 1,weprovide a dimensional shape function R(z) versus thedimensional value of the axial variable of an artery based on the actual dimensional data determined from the experimentally values of the cross-sectional area of the artery of a human (Back et al., 1984). Here, as in more realistic cases, the blood flow is composed of a suspension of erythrocytes (red cells) in plasma. Fig. 1. Experimentally collected data for the atherosclerosis (Back et al., 1984) 2. Formulation and analysis We consider the problem of axisymmetric flow of blood in a catheterized artery in the form of a circular cylindrical annulus tube with the outer radius R0 (radius of the artery) and the inner radius r1 (radius of the catheter) and in the presence of an atherosclerosis whose shape (Fig. 1) is determined from the experimentally collected data (Back et al., 1984). The artery length is assumed to be sufficiently large in comparison to its radius so that the end effects can be neglected. The two-phase flow system in a catheterized artery is based on the original governing equ- ations for themass conservation andmomentum (Batchelor, 1970) for both fluid plasma and the suspended particles (red cells) as their steady axisymmetric form in the cylindrical coordinate system with the axial direction along the co-axial direction of the catheterized artery are given by (Srivastava and Rastogi, 2010) (1−C)ρf ( uf ∂uf ∂z +vf ∂uf ∂r ) =−(1−C) ∂P ∂z +(1−C)µs∇2uf +CS′(up−uf) (1−C)ρf ( uf ∂vf ∂z +vf ∂vf ∂r ) =−(1−C) ∂P ∂r +(1−C)µs∇2vf +CS′(vp−vf) 1 r ∂ ∂r [(1−C)vf]+ ∂ ∂z [(1−C)uf] = 0 (2.1) Two-phase flow in a catheterized artery with atherosclerosis 411 and Cρp ( up ∂up ∂z +vp ∂up ∂r ) =− ∂P ∂z +CS′(uf −up) Cρp ( up ∂vp ∂z +vp ∂vp ∂r ) =− ∂P ∂r +CS′(vf −vp) 1 r ∂ ∂r (Crvp)+ ∂ ∂z (Cup)= 0 (2.2) Here ∇2 ≡ [(1/r)(∂/∂r)(r∂/∂r)+∂2/∂z2] is the Laplacian operator, r (r1 ¬ r ¬R0) and z are the cylindrical coordinates with the axial variable z along the tube axis and the radial variable r along the direction perpendicular to the tube axis, subscripts f and p refer to fluid (plasma) and particle (erythrocyte) quantities, respectively, u and v are the axial and radial velocity components, respectively, ρ is density, P is pressure, C is the volume fraction density of the particles, refers here as the hematocrit % in the blood, and the expressions for the viscosity of suspension ?s and the drag coefficient of interaction S′ have been chosen to be (Srivastava, 1996; Srivastava and Srivastava, 2009) µs = µ0 1−mC m=0.07exp [ 2.49C+ 1107 T exp(−1.69C) ] S′ =4.5 µ0 a20 4+3 √ 8C−3C2+3C (2−3C)2 (2.3) where µ0 is the plasma viscosity, a0 is the radius of a red cell and T is absolute temperature measured in Kelvin. Based on the reasonable suggestion by Charm and Kurland (1974), the expression for plasma viscosity given by (2.3)1,2 is accurate up to 60% hematocrit (C = 0.6), and expression (2.3)1,2 was derived first by Tam (1969) representing classical Stokes drag valid for a small particle Reynolds number. We consider governing equations (2.1)-(2.3) for the blood flow in the axisymmetric form and use the cylindrical coordinate system with r as the radial variable, z as the axial variable and with the z-axis along the axis of cylindrical artery tube, where a catheter in the form of a tube with a small radius but along the axis of artery is placed in the artery. The inside boundary of the artery is partially structured along a distance L0 due to the presence of atherosclerosis (Fig. 2). In Fig. 2, where the flow system and the geometry is shown in the cylindrical annulus, the catheterized arterial tube is given over a distance L= 2d+L0 in the axial direction, δ is the maximum height of atherosclerosis into the lumen, which appears at a particular location in the axial direction. In particular, we refer to a location corresponding to a value very close to the maximum height δ of the atherosclerosis as the critical height such as the location at a distance z= d+L0/2 from the origin of the coordinate system. Fig. 2. Flow geometry in a catheterized artery with atherosclerosis 412 J.M. Ledesma et al. We now first non-dimensionalize governing equations (2.1)-(2.3) using U, L0, R0, δ and µ0UL0/δ 2 as scales for thevelocity, axial length, radial length, rate of radial change andpressure, respectively,where U is themaximumvelocity for theunidirectional flow ina cylindrical annulus (White, 1991).Next,we simplify thedimensionless formsof governing equations (2.1)-(2.3) under the reasonable conditions formild atherosclerosiswith δ/R0 ≪ 1, unidirectional flowassumption (White, 1991)where the axial velocity component dominates over the radial velocity component, and subjected theassumptions that the inertial terms ingoverningmomentumequations (2.1)1,2, (2.2)1,2 are small and Re(δ/L0)≪ 1. Under these conditions and assumptions (Srivastava and Rastogi, 2010), the pressure is only a function of z and (2.1)-(2.3) lead to simpler equations. These simpler equations in non-dimensional forms are given below using the same symbols for the variables as their dimensional ones for simplicity of notations (1−C) ∂P ∂z = 1−C 1−mC 1 r ∂ ∂r ( r ∂uf ∂r ) +CSβ2(up−uf) dP dz =Sβ2(uf −up) S=4.5 4+3 √ 8C−3C2+3C (2−3C)2 (2.4) where β= δ/a0. Equations (2.4) are subjected to the following no slip boundary conditions uf =0 on r= r1 and uf =0 on r=R(z) (2.5) Using (2.4)2.3 for (uf −up) in (2.4)1 and integrating twice with respect to r andmaking use of the boundary conditions given in (2.5), we find uf =− 1−mC 4(1−C) dp dz [ R2−r2+ (R2−r21)ln r R ln R r1 ] (2.6) The expression for the axial velocity for the red cells is then found from (2.4)2,3 in terms of the axial velocity for the plasma in the form up =uf − 1 Sβ2 dp dz (2.7) Since both expressions for the axial velocity of plasma and red cells given by (2.6) and (2.7) are in terms of the unknown pressure gradient (dP/dz), we obtain an expression for the pressure gradient by assuming a prescribed volume flow rate in the annulus given by Q=2π R ∫ r1 r[(1−C)uf +Cup] dr (2.8) Using (2.6) and (2.7) in (2.8) and solving for the pressure gradient, we find dP dz =− Q π 16 (R2−r21) [ 16c Sβ2 +2 R2−r21 1−C (1−mC) ( R2+r21 R2−r21 − 1 ln R r1 )] (2.9) The flow resistance, referred to as the impedance λ, is given by λ= ∆P Q (2.10) where ∆p is the pressure drop across the length 1+2b given by ∆P =P(0)−P(1+2b) = 1+2b ∫ 0 dP dz dz (2.11) and b= d/L0. Two-phase flow in a catheterized artery with atherosclerosis 413 From (2.6), we find the wall shear stress τw to be τw =− ∂uf ∂r ∣ ∣ ∣ ∣ r=R(z) =− R(1−mC) 2(1−C) dP dz − 1−mC 4R(1−C) dP dz r21 −R 2 ln R r1 (2.12) and the shear stress force F over the surface of the artery from z = 0 to z = 1+2b is then given by F =2π 1+2b ∫ 0 τw dz (2.13) The value of the wall shear stress at a different location where the atherosclerosis is located can then be found from (2.12) for a specific z value. 3. Results and discussion We carried out numerical calculations of various expressions obtained in the previous section for several different values of the hematocrit parameter C, the radius of the catheter and over a range of values of the axial and possibly radial variables. For all the calculations, we set the volume flow rate Q=1, b=0.5 and β= δ/0.004, where δ=1−minR. Figures 3-5 present results for dP/dz (axial rate change of the blood pressure in the cathe- terized artery). Figure 3 presents the pressure gradient versus the axial variable for the catheter radius=0.2 and for several values of the hematocrit parameter. It can be seen from this figure that the blood pressure gradient is negative implying that the blood pressure force is in the direction of the positive z-axis. The blood pressure force does not vary with respect to the axial variable at the axial locations outside the atherosclerosis zone. However, the magnitude of the blood pressure force increases with the atherosclerosis effect in the atherosclerosis zone. Themagnitude of the pressure force also increases with the hematocrit effect. These results are physically and bio-medically reasonable since higher percentage of the blood cells amount in the plasma as well as more severity of the atherosclerosis can intensify the blood pressure force in the artery. Fig. 3. Axial rate of change of blood pressure versus axial variable for catheter radius=0.2 and several values of hematocrit parameters In Figs. 4 and 5, similar results to those in Fig. 3 are presented but for higher values of the catheter radius. It can be seen that the magnitude of the blood pressure force in the artery increases with the catheter radius, which again makes sense physically since higher catheter radius implies smaller annulus gap leading to a higher value of such a force. 414 J.M. Ledesma et al. Fig. 4. The same as in Fig. 3 but for catheter radius=0.4 Fig. 5. The same as in Fig. 3 but for the catheter radius=0.5 Figure 6 presents the axial velocity of the blood plasma versus the axial variable for the hematocrit parameter=0.4, the radial variable=0.6 and for two different values of the catheter radius. It can be seen from this figure that the plasma velocity is positive, which makes sense since the blood pressure force is in the positive direction of the axis of the catheterized artery system. The plasma velocity is constant outside the atherosclerosis zone, while it is variable and itsmagnitude increases with the atherosclerosis effect. This result is reasonable since higher atherosclerosis effect increases the strength of the pressure force leading to higher blood plasma speed. The plasma velocity also increases with the catheter radius, which is reasonable since higher catheter radius decreases the annulus gap leading to higher plasma speed.Our additional generated data for the plasma velocity versus different values of the hematocrit percentage indicated insignificant changes of the plasma speed with respect to the hematocrit effect. Fig. 6. Axial velocity of plasma versus the axial variable for C =0.4, r=0.6 and r1 =0.4 Figure 7 present the values of the plasma velocity versus the radial variable for the catheter radius=0.2, the axial value of z=1 and the value of 0.4 for the hematocrit parameter. It can be seen from this figure that the plasma velocity satisfied its zero no-slip conditions at the two Two-phase flow in a catheterized artery with atherosclerosis 415 boundaries of the catheterized artery system, while it has amaximum value at some location in the annulus away from the boundary. Fig. 7. Axial velocity of plasma versus ,r (r1 ¬ r¬R) for r1 =0.2, z=1 and C =0.4 Figure 8 presents the axial velocity of the red cells versus the axial variable for a fixed radial location (r=0.6), for thevalueof 0.4 for thehematocrit parameter and for twovalues 0.2 and0.4 of the catheter radius. As can be seen from this figure, the red cells velocity is higher for higher value of the catheter radius. The speed of the red cells is constant outside the atherosclerosis zone, but it varies significantly in the atherosclerosis zone and increases significantly at the locations where the atherosclerosis is more severe. Fig. 8. Axial velocity of red cells versus the axial variable for r=0.6,C =0.4 and two different r1 values 0.2 and 0.4. Fig. 9. Axial velocity of the red cells versus the radial variable for z=1,C =0.4 and r1 =0.2 Figure 9 presents the speed of the red cells versus the radial variable for z=1,C =0.4 and r1 =0.2. It can be seen from this figure that the highest speed of cells is near mid-region of the annulus, the rate of increase of the speed with respect to the radial variable is higher near the catheter, while the rate of decrease of the speed is smaller near the artery wall. Our additional generated data for the variation of the velocity of the red cells with respect to the hematocrit parameter indicated that the magnitude of such a velocity increases with C, even though such increase is not very significant. 416 J.M. Ledesma et al. Figure 10 presents the impedance versus the hematocrit parameter for values of 0.2, 0.4 and 0.6 of the catheter radius. It can be seen from this figure that the impedance increases with respect to the hematocrit parameter especially noticeably for C < 0.2 and C > 0.4. Thus, the flow resistance an increases with increase in the hematocrit character of the blood especially for a relatively small percentage of the red cells or a large percentage of such cells in the plasma. The flow resistance also is higher if the catheter radius is larger. Fig. 10. Impedance versus C for r1 =0.2, 0.4 and 0.6 Figure 11 presents the wall shear stress versus the axial variable for the value of 0.4 for the hematocrit parameter and for two values 0.2 and 0.4 of the other catheter radius. It can be seen from this figure that the shear stress increases with the catheter radius. The shear stress is constant outside the atherosclerosis zone, while it varies with respect to the axial variable and increases significantly in the atherosclerosis zone, and such an increase intensifies if the atherosclerosis becomes more severe. Although the shear stress was found to increase with the hematocrit parameter, such an increase was found to be weak. Fig. 11.Wall shear stress versus z for C =0.4 and r1 =0.2, 0.4 Fig. 12.Wall shear stress versus catheter radius for z=1 and C =0.4 Figure 12 presents the shear stress versus the radius of the catheter for z=1 and C =0.4. It can be seen from this figure that the shear stress increases with the radius of the catheter, and its rate of increase with respect to the catheter radius also increases with the catheter radius. Two-phase flow in a catheterized artery with atherosclerosis 417 4. Conclusion We investigated the catheterized arterial blood flow in the presence of the atherosclerosis effect. Our formulation was based on a two-phase blood fluid flow model of combination of plasma and red cells.We calculated steady cases for important quantities, such as the pressure gradient force, plasma velocity, red cells velocity, impedance and wall shear stress in the presence of the steady atherosclerosis effect, which is taken into account from the available experimental results for a human. These quantities were evaluated at several values of the hematocrit parameter C, the axial variable z or the radial variable r and the catheter radius. We found, in particular, that the pressure gradient force, plasma velocity, red cells velocity andwall stress are significant in atherosclerosis zone, and they become stronger if the effect of the atherosclerosis is stronger. These conditions could be relevant to the corresponding conditions for certain patients, which can require higher care and attention. The extension of the present paper to the cases of unsteady arterial flow in the presence of the pulse frequency, whose calculated results turn out to depend on the results of the present paper, will be presented in a forthcoming paper. Another important extension of the present study can be formoremedically realistic finite systemswith cases conforming tomoremedically generated data in order to identify the components of arterial blood flow diseases, which can improve health conditions of the corresponding patients. Acknowledgement This researchwas supported by UTPA-FRC2011-2012 grant. References 1. Back L.H., 1994, Estimatedmean flow resistance during coronary artery catheterization, Journal of Biomechanics, 27, 169-175 2. Back L.H., Cho Y.I., Crawford D.W., Cuffel R.F., 1984, Effect of mild atherosclerosis on flow resistance in a coronary artery casting of man,Transactions of the ASME, 106, 48-53 3. Back L.H., Kwack E.Y., Back M.R., 1996, Flow rate-pressure drop relation to coronary an- gioplasty: Catheter obstruction effect, Journal of Biomedical Engineering, 118, 83-89 4. Batchelor G., 1970,An Introduction to Fluid Dynamics, Cambridge University Press, UK 5. Charm S.E., Kurland G.S., 1974,Blood Flow and Microcirculation, JohnWiley, NewYork 6. Kanai H., Lizuka M., Sakamotos K., 1970, One of the problem in the measurement of blood pressure by catheterization:Wave reflection at the tip of catheter,Med. Biol. Engr., 28, 483-496 7. Mishra B.K., Panda T.C., 2005, Non-Newtonian model of blood flow through an arterial ste- nosis,Acta Ciencia Indica, 31, M(2), 341-348 8. Riahi D.N., Roy R., Cavazos S., 2011, On arterial blood flow in presence of an overlapping stenosis,Mathematical and Computer Modelling, 54, 2999-3006 9. Srivastava V.P., 1996, Two-phase model of blood flow through stenosed tubes in the presence of a peripheral layer: Applications, Journal of Biomechanics, 29, 1377-1382 10. Srivastava V.P., Rastogi R., 2010, Blood flow through a stenosed catheterized artery: Effects of hematocrit and stenosis shape,Computer and Mathematics with Applications, 59, 1377-1385 11. SrivastavaV.P., Rastogi R.,Mishra S., 2010,Non-Newtonian arterial blood flow through an overlapping stenosis,Applications and Applied Math.: An. Int. Journal, 5, 1, 225-238 12. Srivastava V.P., Srivastava R., 2009, Particulate suspension blood flow through a narrow catheterized artery,Computer and Mathematics with Applications, 58, 227-238 418 J.M. Ledesma et al. 13. Tam C.K.W., 1969, The drag on a cloud of spherical particles in low Reynolds number flows, Journal of Fluid Mechanics, 38, 537-546 14. Venkateswarlu K., Rao J.A., 2004, Numerical solution of unsteady blood flow through an indented tube with atherosclerosis, Indian Journal of Biochemistry and Biophysics, 41, 241-245 15. White F.M., 1991,Viscous Fluid Flow, Second Edition,McGraw-Hill, Inc., NewYork 16. Young D.F., Tsai F.Y., 1973, Flow characteristics in model of arterial stenosis-steady flow, Journal of Biomechanics, 6, 395-410 Dwufazowy przepływ krwi w tętnicy ze zmianami miażdżycowymi po wykonanym zabiegu koronarografii Streszczenie W pracy omówiono problem przepływu krwi w tętnicy po zabiegu koronarografiiw obecności zmian miażdżycowych, opierając się na osiągalnych danych z badań klinicznych. Miażdżycą nazywamy stan, w którym ściana tętnicy pogrubia się dowewnątrzwskutek odkładania się tłuszczy, głównie cholesterolu. W standardowej metodzie leczenia miażdżycy stosuje się zabieg koronarografii polegający na wprowa- dzeniu cewnika do upośledzonej tętnicy.W pracy opisano przepływ krwi w przekroju tętnicy za pomocą dwufazowegomodelu odzwierciedlającego zawiesinę czerwonych ciałek krwiw osoczu. Sprzężone, różnicz- kowe równania przepływu płynu (osocza) i ruchu cząstek (czerwonych ciałek) rozwiązano analitycznie w stopniu akceptowalnie przybliżonym. Tak istotne wielkości, jak prędkość przepływu osocza, prędkość czerwonych ciałek, ciśnienie krwi, impedancja (opory przepływu) oraz naprężenia ścinające w ścianie tętnicy obliczono dla różnych rozmiarów cewnika, osiowego rozkładu złogówmiażdżycowych oraz hema- tokrytu wywołanego dwufazową kombinacją czerwone ciałka-osocze w badanym układzie krwionośnym. Manuscript received June 9, 2012; accepted for print July 25, 2012