# Steady Flow of Couple-Stress Fluid in Constricted Tapered Artery: Effects of Transverse Magnetic Field, Moving Catheter, and Slip Velocity

Authors: Hamzah Bakhti; Lahcen Azrar
Journal: Journal of Applied Mathematics
 (2016)

---

## Abstract

Steady flow of a couple-stress fluid in constricted tapered artery has been studied under the effects of transverse magnetic field, moving catheter, and slip velocity. With the help of Bessel’s functions, analytic expressions for axial velocity, flow rate, impedance, and wall shear stress have been obtained. It is of interest to note that these solutions can be used for different types of fluid flow in tubes and not only the case of blood. The effects of various geometric parameters, the parameters arising out of the fluid considered and the magnetic field, are discussed by considering the slip velocity, the catheter velocity, and tapering angle. The study of the above model is very important as it has direct applications in the treatment of cardiovascular diseases.


---

## Body

## 1. Introduction

Catheters are semirigid, thin tubes made from medical grade materials serving a broad range of functions and can be inserted in the body to detect and identify diseases or perform a surgical procedure inside the heart, brain, arms, legs, or lungs [1]. But the problem is how we can know that the catheter is going the right way or where is it. Most hospitals use X-ray machine, mounted above the couch to check the catheter’s progress. Low-dose X-rays are used to monitor the progress of the catheter tip. But if we want to navigate the catheter’s passage in real time, some ten to twenty X-ray images are made every second. Even though the radiation dose involved is very low, no radiation at all is much healthier for both the patient and the medical staff who run the risk of being exposed to the radiation on daily basis. Researchers developed a magnetic navigation system for medical instruments such as catheters and guide wires. The catheters with magnetic tips are steered within the patient, without the need for an electrophysiologist to maneuver the catheter or guide wire placement manually. A magnetic sensor on the tip of the medical instrument measures an external magnetic field and reports exactly where the tip of the instrument is located [2].Each year, heart disease is at the top of the list of the country’s most serious health problems. Statistics show that cardiovascular disease in some countries, like America, is the first health problem and the leading cause of death (recent statistics released by the American Heart Association). Catheters are often used to treat stenosis (partial occlusion of the blood vessel), aneurysm (dilation of the blood vessel resulting in stretching of the vessel wall), and embolism (complete occlusion of a blood vessel by a blood clot or some other particle), Figure1. The catheters may contain a balloon that can be used to stretch the blood vessel or a rinsing mechanism or laser device to remove embolism. They are also used to measure blood pressure in situ. Sometimes a catheter is used to deliberately create embolism, for example, to stop the flow of blood to a tumor in order to mortify it. An aneurysm can also be treated by means of embolisation [3].Figure 1
Catheters are used to treat stenosis (partial occlusion of the blood vessel), aneurysm (dilation of the blood vessel resulting in stretching of the vessel wall), and embolism (complete occlusion of a blood vessel by a blood clot or some other particle) (a reproduction from Delft Outlook 2004).The study of blood flow through different types of arteries is of considerable importance in many cardiovascular diseases; one of them is atherosclerosis. The blood flow through an artery has drawn the attention of researchers for a long time up to now, due to its great importance in medical sciences. Under normal conditions, blood flow in the human circulatory system depends upon the pumping action of the heart and this produces a pressure gradient throughout the arterial network.Many researchers have studied blood flow in the artery by considering blood as either Newtonian or non-Newtonian fluids; since blood is a suspension of red cells in plasma, it behaves as a non-Newtonian fluid at low shear rate. In the stenosed condition, substantial reduction in the lumen of an artery results in size effects (ratio of haematocrit to vessel diameter), which influences flow characteristics significantly. To study the size effect in the fluid flow, Stokes [4], Eringen [5], and Cowin [6] have proposed continuum model. Microcontinuum structure of the fluid is also referred to as couple-stress fluid which is proposed by Stokes. The foreground of couple-stress fluid is to introduce size dependent effects as mentioned above that is not present in other classical viscous fluids. Because of its mathematical simplicity, it has been widely used by number of researchers; Chaturani and Upadhya [7] have studied the pulsatile flow of couple-stress fluid through circular tubes. The Poiseuille flow of couple-stress fluid has been examined by Chaturani and Rathod [8].The application of Magnetohydrodynamics in physiological problems is of growing interest and it is very important from both theoretical and practical points of view. The application of Magnetohydrodynamics in physiological problems is of growing interest. The flow of blood can be controlled by applying appropriate quantity of magnetic field. Kollin [9] has coined the idea of electromagnetic field in the medical research for the first time in the year 1936. Korchevskii and Marochnik [10] have discussed the possibility of regulating the blood movement in human system by applying magnetic field. Rao and Deshikachar [11] have investigated the effect of transverse magnetic field in physiological type of flow, through a uniform circular pipe. Vardanyan [12] showed that the application of magnetic field reduces the speed of blood flow. A mathematical model for two-layer pulsatile flow of blood with microorganism in a uniform tube under low Reynolds number and magnetic effect has been studied by Rathod and Gayatri [13]. Thus, all these researchers have reported that the effect of magnetic field reduces the velocity of blood.There are many treatments available for diagnosing and treating constricted vessels. Catheterization (thin, flexible tube) is one of them, in which balloon angioplasty is a specialized form of catheterization. These procedures are widely used in the medical field for treating the atherosclerosis. Insertion of the catheter in a tube creates an annular region between inner wall of the artery and outer wall of the catheter which influences the flow field such as pressure distribution and shear stress at the wall. In view of its immense importance, the effect of the catheter on physiological parameters was discussed by the researchers [14, 15].The shapes of the stenosis in the above aforesaid studies have been considered to be radially symmetric or asymmetric. But while stenosis is maturing, it may grow up in series manner, overlapping with each other, and it would appear like x-shape. Riahi et al. [16] observed the pressure gradient force, flow velocity, impedance, and wall shear stress in the overlapping stenotic zone at the critical height and at the throats of such stenosis. Here they considered steady nature of flow. Srivastava and Mishra [17] explored the arterial blood flow through an overlapping stenosis by treating the blood as a Casson fluid. They figured impedance and shear stress for different stenosis heights. Chakravarty and Mandal [18] discussed the effects of the overlapping stenosis under low shear rate flow.The presence of red cell slip at the vessel wall was recommended theoretically by Vand [19], experimentally by Bennett [20] and Nubar [21], Chaturani and Biswas [22], and so forth. They used slip velocity at the wall in their analysis. Lately, Ponalagusamy [23, 25] has developed mathematical models for blood flow through stenosed arterial segment, by taking a slip velocity condition at the constricted wall. Thus, it is very appropriate to consider slip velocity at the wall of the stenosed artery in blood flow modeling. The catheter is actually moving into the body, so it is recommended to consider a nonnull velocity at the catheter wall; in this paper, we consider that the catheter is moving in the z-axis with a fixed velocity.In the present analysis a mathematical model for the steady blood flow through tapered stenosed artery under the influence of a moving catheter, slip velocity, and a magnetic field is presented by considering blood as a couple-stress fluid in a circular tube. It is assumed that the magnetic field along the radius of the pipe is present, no external electric field is imposed, and magnetic Reynolds number is very small. The motivation for studying this problem is to understand the blood flow in an artery under the effect of magnetic field alongside with the catheter inserted into the blood vessel also when the fatty plaques of cholesterol and artery clogging blood clots are formed in the lumen of the artery.The main aim of this work is to study these phenomena, obtain analytic expressions for axial velocity and shear stress, and also study the effect of magnetic field (Hartmann numberH) and couple-stress parameter (β) on the velocity and the effects of Hartmann number on the fluid velocity. Hence, the present mathematical model gives a simple form of velocity expression for the blood flow so that it will help not only people working in the field of physiological fluid dynamics but also the medical practitioners.


## 2. Problem Formulation

Let us consider a two-dimensional steady flow of blood through a rigid tapered stenosed tube by considering blood as an electrically conducting, incompressible, couple-stress fluid. The magnetic field is acting along the radius of the tube. The magnetic Reynolds number of the flow is assumed to be sufficiently small that the induced magnetic and electric fields can be neglected [24]. The catheter is assumed to be moving in the z-axis direction.Mathematical model of blood flow through a tapered stenosed arterial segment in the presence of a moving catheter and a magnetic field is to be built to study the impact of various geometric, Hartman, and fluid parameters on physiological parameters. The geometry of the tapered stenosed artery is shown in Figures2 and 3 and is expressed mathematically with inputs from [26] as(1)Rz=R0+ζz1-ϵnn/n-1n-1L0nL0n-1z-d-z-dn,d≤z≤d+L0R0+ζz,otherwise,where R0 is the radius of the annular region in case of nontapered artery in the nonconstricted domain, Rc is the radius of the catheter, L0 is the stenosis length, d indicates the location of stenosis, ϵ is the maximum height of the stenosis into the lumen, and ζ=tan⁡ϕ is the tapering parameter which represents the slope of the tapered vessel with ϕ being the tapering angle. ϕ<0, ϕ>0, and ϕ=0 are for converging taper, diverging taper, and no-taper, respectively.Figure 2
2D view of a catheterized stenosed artery.Figure 3
2D view of a tapered stenosed artery.The current densityJ is expressed by(2)J=σE+V×B,where E is the electric field intensity, σ is the electrical conductivity, B is the magnetic flux intensity, and V is the velocity vector. In the momentum equation, the electromagnetic force Fm is included and is defined as(3)Fm=J×B=σE+V×B×B.The conservation equations which govern the couple-stress fluid flow including a Lorentz force can be written in the following form:(4)∇·V=0,ρV·∇V=-∇p+λ+μ∇∇·V+ηΔ∇∇·V-ηΔ2V+μΔV+σV×B×B+ρf+12∇×ρl,where Δ is the Laplacian operator. For couple-stress fluid, shear stress tensor is not symmetric. The force stress tensor τ and the couple-stress tensor M that arises in the theory of couple-stress fluids are given by(5)τ=-p+λ∇VI+μ∇V+∇VT+12I×∇M+ρl,M=mI+2η∇∇×V+2η′∇∇×VT,where ρ is the density, V is velocity vector, p is the pressure, λ and μ are the viscosity coefficients, and η and η′ are couple-stress coefficients. f is a body force and l is a body couple moment. Further, the materials constants μ, λ, η, and η′ satisfy the following inequalities:(6)μ≥0,3λ+2μ≥0,η≥0,η≥η′.As the flow is steady and incompressible, in the absence of body force and body couple moment (4) reduce to(7)∇·V=0,ρV·∇V=-∇p+μΔV-ηΔ2V+σV×B×B.The flow is considered to take place under the influence of externally applied magnetic field in axial directionz. Under these assumptions, the governing equations may be written in the cylindrical coordinates system as follows.Equation of continuity is(8)∂u∂r+ur+∂w∂z=0.Equation of radial momentum is(9)ρu∂u∂r+w∂u∂z=-∂p∂r+μ∂2u∂r2+1r∂u∂r-ur2+∂2u∂z2-η∂4u∂r4+2∂4u∂r2∂z2+∂4u∂z4+2r∂3u∂r3+2r∂3u∂r∂z2-3r2∂2u∂r2-2r2∂2u∂z2+3r3∂u∂r-3r4u.Equation of axial momentum is(10)ρu∂w∂r+w∂w∂z=-∂p∂z+μ∂2w∂r2+1r∂w∂r+∂2w∂z2-η∂4w∂r4+2∂4w∂r2∂z2+∂4w∂z4+2r∂3w∂r3+2r∂3w∂r∂z2-1r2∂2u∂r2+1r3∂w∂r-σB02w,where B0 is the external transverse magnetic field. The nondimensional parameters are obtained as follows:(11)z∗=zL0,r∗=rR0,w∗=ww0,ϵ∗=ϵR0,ζ∗=ζL0R0,u∗=L0uw0R0,p∗=r02pw0L0μ,and then the dimensionless form of geometry (1) after dropping the stars is obtained as follows:(12)Rz=1+ζz1-ζnn/n-1n-1z-γ-z-γn,γ≤z≤γ+11+ζz,otherwise,where γ=L0/L1 and w0 is a typical axial velocity.Also, (8)–(10), on dropping the stars, become(13)∂u∂r+ur+∂w∂z=0,Reξ3u∂u∂r+w∂u∂z=-∂p∂r+ξ2∂2u∂r2+1r∂u∂r-ur2+ξ2∂2u∂z2-ξ2β2∂4u∂r4+2r∂3u∂r3-3r2∂2u∂r2+3r3∂u∂r-3r4u+ξ22r∂3u∂r∂z2+2∂4u∂r2∂z2+ξ2∂4u∂z4-2r2∂2u∂z2,Reξu∂w∂r+w∂w∂z=-∂p∂z+∂2w∂r2+1r∂w∂r+ξ2∂2w∂z2-1β2∂4w∂r4+2r∂3w∂r3-1r2∂2w∂r2+1r3∂w∂r+ξ4∂4w∂z4+2ξ2∂4w∂r2∂z2+1r∂3w∂r∂z2-H2w,where ξ=R0/L0, Re=ρu0R0/μ is the Reynolds number, β2=R02μ/η is the couple-stress fluid parameter, and H2=B02R02σ/μ is the Hartmann number.Under the assumption of mild stenosis that isϵ/R0≤1 and further assuming that ξ=R0/L0≤1, (13) get transformed into(14)∂p∂r=0,(15)∂p∂z=∂2w∂r2+1r∂w∂r-1β2∂4w∂r4+2r∂3w∂r3-1r2∂2w∂r2+1r3∂w∂r-H2w.Equation (15) can be written as(16)∂p∂z=Δw-1β2Δ2w-H2w,where(17)Δ=∂2∂r2+1r∂∂r.It can be seen that the pressure variation depends only on the axial variable. The pressure gradient∂p/∂z is produced by the pumping action of the heart.The corresponding nondimensional boundary conditions are as shown below:(18)w=vatr=Rz,w=gatr=Rc,∂2w∂r2-ωr∂w∂r=0atr=Rz,r=Rc,where v represents the slip velocity at the artery wall, g is the velocity of the moving boundary (catheter), and ω=η′/η is the parameter associated with the couple-stress fluid. No couple-stress effects will be present if η′=η, which is equivalent to saying that the couple-stress tensor is symmetric.


## 3. Solution of the Problem

As a solution of (16) with the pressure gradient being taken as a constant, take(19)-∂p∂z=ks.Letn12+n22=β2 and n12n22=β2H2, which are given by (20)n12=β2+β4-4β2H22,(21)n22=β2-β4-4β2H22;then, (16) is simplified to the form(22)Δ-n12Δ-n22w=β2ks.The solution of the above equation is obtained as(23)w=c1zI0n1r+c2zK0n1r+c3zI0n2r+c4zK0n2r+1H2ks,where I0 and K0 are modified Bessel function of order zero, first and second kinds, respectively. ci for i=1,2,3,4 are calculated numerically by solving the algebraic system obtained from the boundary conditions.Volumetric flow rateQ, in steady flow through a tube, is obtained by integrating the velocity profile over a cross section of the tube. The nondimensional volumetric flow rate Q across the radial distance is expressed as(24)Q=∫RcRz2rwdrwhich is obtained in the form Q=-2ksF(Rc,R(z)), where(25)FRc,Rz=d1n1RI1n1R-RcI1n1Rc+d2n1RcK1n1Rc-RK1n1R+d3n2RI1n2R-RcI1n2Rc+d4n2RcK1n2Rc-RK1n2R+R2-Rc22H2,where di=-ci/ks, for i=1,2,3,4. The pressure drop Δp across the stenosis between the sections z=0 and z=L is obtained by(26)Δp=Q2∫0L1FRc,Rzdz.The resistance to the flow (impedance) is obtained from(27)Rf=ΔpQ=12∫0L1FRc,Rzdz=12∫0d1FRc,Rzdz+∫dd+L01FRc,Rzdz+∫d+L0L1FRc,Rzdz.The dimensionless form of (20) is(28)Rf=Γ2∫0γ1FRc,Rzdz+∫γγ+11FRc,Rzdz+∫γ+11/γ1FRc,Rzdz,where Γ=L0/d.The shear stressτrz is calculated using the expression(29)τrz=τrzS+τrzA,where the symmetric and skew-symmetric parts of the stress are given by(30)τrzS=μ∂w∂r,(31)τrzA=∂mrθ∂r+mrθ+mθrr,respectively. Here,(32)mrθ=14η′r∂w∂r-η∂2w∂r2,mθr=14ηr∂w∂r-η′∂2w∂r2.Hence, the dimensionless shear stress for the artery is given by(33)τrz=n1c1I1n1r-c2K1n1r1-n124n12+n22+n2c3I1n2r-c4K1n2r1-n224n12+n22.Thus, the shear stress at the wall can be computed from (33) by taking r=R(z).


## 4. Results and Discussion

The study of blood flow through catheterized stenosed tapered artery with the presence of a transverse magnetic field involves the integration of various geometric and fluid variables, which influences the physiological parameters such as the fluid velocity, rate flow, and wall shear stress. Closed form solutions are obtained in terms of modified Bessel’s functions. The physiological dimensionless quantities such as the fluid velocity in the stenosis region and the wall shear stress at the maximum height of the stenosis are computed numerically for various values of the fluid and geometric parameters using the programMATLAB. The parameters considered are n (shape parameter), ζ (tapered parameter), Rc (catheter radius), Γ (stenosis length), β (couple-stress fluid parameter), H (the Hartmann number), v (slip velocity), and g (catheter wall velocity). The obtained results are analyzed graphically. The results obtained in this study are in good agreement with those reported in the literature.Figures4–6 illustrate the variation of axial velocity profile for different values of n, ζ, and β. It is observed from Figure 4 that the velocity profile increases for increasing the shape parameter and that actually this occurred due to the change of the stenosis height. It is observed from Figure 5 that the velocity increases by the increase in the tapered parameter ζ and the same goes for the couple-stress fluid parameter β, Figure 6. It is to be noted as β→∞, the properties of couple-stress in the fluid vanish and hence behave like a Newtonian fluid. Hence it is understood that the velocity is low in couple-stress fluid when compared to that of Newtonian fluid.Figure 4
Variation of axial velocityw with respect to n when ζ=0, ϵ=0.2, β=5, H=1, v=0.01, g=0.03, and Rc=0.1.Figure 5
Variation of axial velocityw with respect to ζ when n=5, ϵ=0.2, β=5, H=1, v=0.01, g=0.03, and Rc=0.1.Figure 6
Variation of axial velocityw with respect to β when n=5, ϵ=0.2, ζ=0, H=1, v=0.01, g=0.03, and Rc=0.1.In general, from Figure7, it can be observed that the axial velocity is decreasing. The effect of the Hartmann number H on the axial velocity is shown in Figure 8. By observing Figure 8, the fluid velocity is decreasing as the Hartmann number is increasing. Further, as H→0, the Hartmann number loses its properties and behaves like a normal blood flow without any magnetic field applied to it. We have concluded that the application of magnetic field reduces the speed of blood flow and that meets the results of Ponalagusamy and Tamil Selvi [25], also indicating that the wall shear stress increases when we increase the Hartmann number. The influence of catheter radius Rc on the axial velocity is directly proportional as shown in Figure 9. From the above result, we understand that as the catheter radius increases, the annular region gets narrowed, which leads to the rise in the obstruction to the flow.Figure 7
Variation of axial velocityw with respect to ϵ when n=5, ζ=0, β=5, H=1, v=0.01, g=0.03, and Rc=0.1.Figure 8
Variation of axial velocityw with respect to H when n=5, ζ=0, β=5, ϵ=0.2, v=0.01, g=0.03, and Rc=0.1.Figure 9
Variation of axial velocityw with respect to Rc when n=5, ζ=0, β=5, ϵ=0.2, v=0.01, g=0.03, and H=1.The axial velocity of the blood is high in case of high slip velocityv or catheter velocity g; it is observed that the slip velocity at the boundary facilitates the fluid flow and the same goes for the moving catheter, Figures 10 and 11.Figure 10
Variation of axial velocityw with respect to v when n=5, ζ=0, β=5, ϵ=0.2, Rc=0.1, g=0.03, and H=1.Figure 11
Variation of axial velocityw with respect to g when n=5, ζ=0, β=5, ϵ=0.2, v=0.01, Rc=0.1, and H=1.The shear stress at the wall is a significant physiological parameter to be considered in the blood flow study. Precise predictions of the distribution of the shear stress at the wall are particularly useful in assimilating the effect of blood flow in arteries in general. The shear stress at the wall is calculated at the maximum height of the stenosis. When we increase the height of the stenosisϵ, the wall shear stress increases, Figure 15. Furthermore, it is observed that the shear stress along the wall is reaching the maximum at the throat of the stenosis.The effect of the Hartmann numberH on the shear stress is shown in Figure 16. As observed from Figure 12, the shear stress at the wall is independent of the shape parameter n as pointed out by Ponalagusamy and Tamil Selvi [25]. The variation of the shear stress at the wall as the tapered parameter increases is depicted in Figure 13. Here, the shear stress at the wall is more significant in the untapered artery compared with the tapered artery which is diverging. Further, the converging tapered artery possesses higher wall shear stress than the untapered and tapered artery which is diverging.Figure 12
Variation of wall shear stressw with respect to n when ζ=0, ϵ=0.2, β=5, H=1, v=0.01, g=0.03, and Rc=0.1.Figure 13
Variation of wall shear stressw with respect to ζ when n=5, ϵ=0.2, β=5, H=1, v=0.01, g=0.03, and Rc=0.1.The slip velocity at the wall of the stenosed artery and the moving catheter significantly influence the shear stress at the wall, which is noticed from Figures18 and 19. Here, it is observed that gain in slip velocity (resp., moving catheter) reduces the shear stress at the wall. The influence of catheter radius on shear stress at the wall is shown in Figure 17. The increasing catheter radius narrows down the lumen of the artery thus resulting in the higher values of shear stress at the wall. As couple-stress fluid parameter β increases, the shear stress at the wall decreases. This behavior is shown in Figure 14. It is observed that the theoretical distribution of shear stress along the wall reaches a maximum at throat of the stenosis and then rapidly decreases in the diverging section. This behavior is the same in terms of the change in both the radius of the catheter Rc and the couple-stress fluid parameter β whose behavior is depicted graphically in Figures 17 and 14, respectively.Figure 14
Variation of wall shear stressw with respect to β when n=5, ϵ=0.2, ζ=0, H=1, v=0.01, g=0.03, and Rc=0.1.Figure 15
Variation of wall shear stressw with respect to ϵ when n=5, ζ=0, β=5, H=1, v=0.01, g=0.03, and Rc=0.1.Figure 16
Variation of wall shear stressw with respect to H when n=5, ζ=0, β=5, ϵ=0.2, v=0.01, g=0.03, and Rc=0.1.Figure 17
Variation of wall shear stressw with respect to Rc when n=5, ζ=0, β=5, ϵ=0.2, v=0.01, g=0.03, and H=1.Figure 18
Variation of wall shear stressw with respect to v when n=5, ζ=0, β=5, ϵ=0.2, Rc=0.1, g=0.03, and H=1.Figure 19
Variation of wall shear stressw with respect to g when n=5, ζ=0, β=5, ϵ=0.2, v=0.01, Rc=0.1, and H=1.


## 5. Conclusion

A mathematical model has been built to discuss the flow of blood through a catheterized asymmetric tapered stenosed artery with slip velocity at the stenosed wall and a moving catheter. Closed form solution is obtained and the effects of various geometric, fluid parameters and magnetic field on the axial velocity of the blood and the shear stress at the wall are studied.There is special importance of couple-stress fluids compared to the Newtonian fluids because of their wide existence such as oil, blood, and polymeric solutions. In view of what is mentioned above, an analytic approach was followed to solve the mathematical model of blood flow through stenosed tapered artery under the assumption of mild stenosis. The resultant observations are summarized as follows:(i)
As the height and the stenosis length are increasing, the obstruction to the flow of blood is increasing.(ii)
Converging tapered artery has more shear stress at the wall than the nontapered and the diverging tapered artery.(iii)
Diverging tapered artery has least shear stress at the wall.(iv)
The increment in the Hartmann number enhances the blood velocity and wall shear stress.(v)
The axial velocity is decreasing while the couple-stress fluid parametersβ and the height of the stenosis ϵ are increasing.(vi)
Three different values of the slip velocity and the velocity of the catheter at the arterial boundary and the catheter wall are considered; they are showing a significant influence on the shear stress at the wall and the axial velocity.The modeling and simulation of the above phenomena are very realistic and are expected to be very useful in predicting the behavior of physiological parameters in the diagnosis of various arterial diseases.



---

*Source: 9289684-2016-02-03.xml*