Mathematical Analysis on Pulsatile Flow through a Catheterized Stenosed Artery

In this paper, the pulsatile flow of blood through an inclined catheterized stenosed artery is analyzed. Perturbation method is used to solve the implicit system of partial differential equations with suitable boundary conditions. Various analytical expressions axial velocity, flow rate, wall shear stress and effective viscosity have been derived with the help of MATLAB for understanding the fluid flow phenomena. The combined effect of catheterization, body acceleration, slip and inclination has been seen by plotting the graph and observed that axial velocity and flow rate increases with the increase in body acceleration, inclination angle and slip velocity while axial velocity diminishes on increasing the catheter radius. Wall shear stress increases with the increase in catheter radius and body acceleration but presence of slip velocity reduces the wall shear stress. Effective viscosity diminishes on increasing body acceleration and inclination angle, whereas slightly augmented in non-inclined stenosed artery.


Introduction
In recent years, the study of blood flow through obstructed arteries has received much attention, due to its ample importance in human cardiovascular system.
The theoretical analysis on blood flow is very useful as it plays a significant role to diagnose and understand many cardiovascular diseases such as coronary thrombosis, angina, pectoris, strokes etc.The reason behind the malfunction of cardio-vascular system is the presence of fats, cholesterol and lipoproteins at the sites of atherosclerotic lesion in the artery.By the development of atherosclerotic plaques that protrude into the lumen, arteries get narrowed and stenosed arteries are formed.Presence of atherosclerosis (stenosis), increased the resistance and therefore blood flow reduced inside the artery and also the remarkable change occurred in pressure distribution and wall shear stress.Bennett [1] has observed that the presence of wall slip at the flow boundaries decreases the apparent (effective) viscosity.Tu and Deville [2] observed the blood flow in diseased conditions.Mathematical models for blood flow through stenosed arterial segment, by taking a velocity slip condition at the constricted wall were developed [3] [4].
Our Body experiences the acceleration while running, jumping and travelling etc., due to which different health related issues arise such as headache, loss of vision, increase in pulse rate etc. Sud and Sekhon [5] studied on pulsatile flow of blood through a rigid circular tube subjected to periodic body acceleration considering blood as Newtonian fluid.In [6] [7] [8], the pulsatility of blood is discussed by treating blood as a Newtonian fluid.Liepsch [9] compared the details of flow behaviors with hemodynamic approach and observed that flow rate varies due to pulsatile nature of blood.Since heart pumps of the blood are periodic in nature, the blood flow in narrow arteries can be assumed as pulsatile.Young [10] observed that blood can be assumed as Newtonian when it passes through narrow arteries (0.02 -0.1 mm) at low shear rate (<10/s) particularly in diseased state while blood being suspension of corpuscles in plasma behaves like a non-Newtonian fluid when it flows through larger arteries at high shear rate [11] [12].
Catheter is a thin hollow flexible tube made from medical grade materials and is used in a broad range of functions.Srinivasacharya and Srikanth [13] observed that catheterization refers to a procedure in which a long, thin, flexible plastic tube (catheter) is inserted into an artery.The insertion of a catheter into a constricted tube (i.e.stenosed artery) results in an annular region between the walls of the catheter and artery.This will alter the flow field, modify the pressure distribution and increase the resistance.MacDonald [14] considered the pulsatile blood flow in a catheterized artery and obtained theoretical estimates for pressure gradient corrections for catheters, which are positioned eccentrically, as well as coaxially with the artery.Karahalios [15] discussed some possible effects of a catheter on the arterial wall.
Cardiac catheterization is a procedure used to diagnose and treat cardiovascular conditions in surgical procedures.During cardiac catheterization, coronary angiography is done by inserting a catheter in an artery or vein in your groin, neck or arm and threaded through your blood vessels to your heart.Back and Denton [16] studied the variation of wall shear stress and discussed its clinical importance in coronary angioplasty.An appropriate size of a catheter is required for each experiment in order to reduce the error due to the wave reflection at the tip of the catheter [17].Sarkar and Jayaraman [18] discussed the changed flow patterns of pulsatile blood flow in a catheterized stenosed artery.Effect of slip velocity on blood flow through a catheterized stenosed artery is studied in [19].
All above mentioned studies are based on the horizontal blood carrying vessels.Although it is well-known that many ducts in physiological systems are no horizontal, few of them have some inclination to the axis.Therefore, a gravita-Journal of Applied Mathematics and Physics tional force has been accounted there due to inclination.Chaturani and Upadhyay [20] investigated the gravity flow of fluid with couple stress along an inclined plane.Vajravelu, Sreenadh and Ramesh Babu [21] showed the theory on peristaltic transport of Herschel Bulkley fluid through an inclined tube.The study on pulsatile flow of blood through rigid inclined circular tubes under the influence of periodic body acceleration is proposed in [22].Maruti Prasad and Radha Krishnamcharya [23] have analyzed a model of steady blood flow through an inclined artery with a stenotic wall.Biswas and Paul [24] developed a mathematical model on the steady flow of blood through an inclined tapered constricted artery assuming blood as Newtonian fluid, with an axial slip velocity at the vessel wall.Pulsatile flow of blood in a catheterized inclined artery [25] with a slip velocity at the stenosed arterial wall under the influence of magnetic field is also discussed.To generalize the problem further, present model has been established to analyze the effect of catheterization, body acceleration, slip velocity and pulsatility on arterial flow characteristics in inclined stenosed artery.The detail study on the present work is discussed through the graph considering blood as Newtonian fluid.

Mathematical Formulation
Consider the pulsatile flow of blood through an inclined axially symmetric ca For solving the model we have assumed that blood flow is laminar, incompressible, one dimensional and fully developed in a rigid circular artery of radius 0 R in which a catheter of radius c R is introduced coaxially and the blood is modeled as a Newtonian fluid.It is assumed that the pulsatile flow in the artery is due to a prescribed periodic pressure gradient along the axis of the artery.It is also assumed that the length of the artery is large enough when compared to its diameter so that the entrance, end and special wall effects can be neglected.
The geometry of the stenosis which is assumed to be manifested in the arterial segment is described as [26] ( ) , in the normal artery region where ( ) R z and 0 R is the radius of the artery with and without stenosis respectively, d is the location of the stenosis, s L is the length of the stenosis and s δ denotes the maximum height of the stenosis in to the lumen such that The Navier-Stokes equations governing the fluid flow is given by, [27].
where v represents the axial velocity along z direction, t is the time, ρ is the density, p is the pressure, τ is the shear stress, β is the inclination angle and ( ) B t is the body acceleration.
Newtonian fluid can be represented by the equation where µ denotes the coefficient of viscosity of blood.
The boundary conditions are ( ) where s v is the axial slip velocity at the stenotic wall and c R is the radius of catheter.
Since the pressure gradient is the function of z and t , therefore can be represented as ( 0 R ).  ( ) where 0 A is the steady state pressure gradient, 1 A is the amplitude of the fluctuating component and both 1 0 , A A are function of z .It can be seen that the radial velocity is very small in magnitude so that it may be neglected for problem with mild stenosis.The frequency of oscillation of the pulsatile flow is denoted by p ω and defined as 2π , where p f is the pulse rate frequency.
The periodic body acceleration ( ) B t in the axial direction is given by ( ) ( ) where the amplitude of body acceleration is 0 a and φ is the phase angle of body acceleration with respect to the pressure gradient.f is assumed to be small so that wave effect can be neglected.
Let us introduce the following non-dimensional variables ( ) ( ) where α is the pulsatile Reynold's number or generalized Womersley fre- quency parameter.
Using non-dimensional variables, Equation (2) becomes Equation ( 4) becomes ( ) On substituting the value of τ in Equation ( 10), we have The boundary conditions ( 5) and ( 6) reduces to ( ) The geometry of an arterial stenosis in non-dimensional form is given by ( ) 1, in the normal artery region The non-dimensional volumetric flow rate is defined by where Effective viscosity e µ defined as Can be expressed in the non dimensional form as ( ) ( )

Method of Solution
In this paper perturbation method is used with a small parameteric value of pulsatile Reynolds number α of the series expansion to solve this system of non linear equations.Since non-dimentionalized Equation (10), Equation (11) has 2 α term which is dependent on time, therefore expanding Equation ( 10), Equa- tion (11) about 2 α .The axial velocity v, shear stress τ are expressed as follows in terms of 2 α (where 2 0 α Substituting Equation (18) in Equation ( 12), we get ( ) where ( ) ( ) ( ) Substituting Equation (18) in Equation (( 12) and ( 13)), we have ( ) Integrating Equations (( 20) and ( 21)) , with the help of Equations (( 22) and ( 23)), the expression for v 0 and v 1 are obtained as With the help of Equations (( 11) and ( 19)), the wall shear stress w τ can be obtained as From Equations (( 16) and ( 26)) the expression for volumetric flow rate ( ) The expression for effective viscosity e µ can be obtained from Equationd (( 17) and ( 26)) as

Numerical Results and Discussion
In the present investigation our objective is, to carry out the combined result of applied body acceleration and slip velocity in an inclined catheterized stenosed artery.The expression of axial velocity, flow rate, wall shear stress and effective    Figure 5 shows that on increasing the catheter radius, axial velocity decreases while slip improves the velocity profile.In Figure 6 and Figure 7, we dealt with flow rate versus axial distance and we found that flow rate increases with the increase in body acceleration, slip.Flow rate slightly increases when artery is inclined and perpendicular to the horizontal axis as compared to the normal (non-inclined) artery.The variation of wall shear stress with axial distance is     for different values of body acceleration has been seen in Figure 10 and found that on increasing body acceleration and inclination angle, effective viscosity diminishes whereas slightly augmented in normal (non-inclined) stenosed artery.

Conclusion
In the present mathematical model, pulsatile blood flow through an inclined stenosed catheterized artery with periodic body acceleration and axial slip velocity at the constricted wall has been considered.The flowing fluid is represented by a Newtonian fluid.Analytic expressions for flow variables and their variations with different flow parameters have been obtained and are represented graphically.The results based on the mathematical analysis and the subsequent numerical evaluation of the flow quantities show that the axial velocity and flow rate increases with the increase in body acceleration, inclination angle and slip velocity while axial velocity diminishes on increasing the catheter radius.Wall shear stress increases with the increase in catheter radius and body acceleration but presence of slip velocity reduces the wall shear stress.Effective viscosity diminishes on increasing body acceleration and inclination angle, whereas slightly augmented in normal (non-inclined) stenosed artery.This model concludes that slip velocity plays a very eminent role in blood flow models through an inclined stenosed artery with catheter.In future, this study may be helpful for the purpose of simulation and validation of different models in different conditions of stenosis and to investigate that which parameter has the most dominating role.

.
We have used cylindrical polar co-ordinates(  )    , , r z φ , whose origin is located on the vessel (stenosed artery) axis and r , z denote the radial and axial co-ordinates respectively.It can be shown that the magnitude of radial velocity is negligibly small and can be neglected for a low mean Reynold number flow problems in case of mild stenosis.
S. U. Siddiqui, C. Awasthi DOI: 10.4236/jamp.2017.591571878 Journal of Applied Mathematics and Physics ( ) frequency in Hz.The frequency of the body acceleration b viscosity are obtained and computed for the fixed values of F = 0.2, e = 1, ϕ = 0.2, α = 0.5, ω = 1 through the MATLAB.Figures2-5reveals the variation of axial velocity with radial distance.Figure2depicts that the magnitude of axial velocity increases on increasing slip velocity and almost double for increased value of body acceleration from B = 0 to B = 2. From Figure3and Figure4, we observed that the axial velocity is greater in inclined artery along with the increase in slip as compared to the normal artery (non-inclined) without slip and body acceleration also enhances the velocity graph in inclined position.