Abstract
An appropriate nonlinear blood flow model under the influence of periodic body acceleration through a multiple stenosed artery is investigated with the help of finite difference method. The arterial segment is simulated by a cylindrical tube filled with a viscous incompressible Newtonian fluid described by the Navier-Stokes equation. The nonlinear equation is solved numerically with the proper boundary conditions and pressure gradient that arise from the normal functioning of the heart. Results are discussed in comparison with the existing models.
1. Introduction
At present the investigation of blood flow analysis in a stenosed artery is very important in the medical domain because of the fact that many of the diseases such as heart attacks and strokes are related to blood flow and the physical characteristic of vessel wall. Nowadays the leading causes of the death in the world are due to heart diseases such as atherosclerosis. Atherosclerosis involves an accumulation of low-density lipoprotein in the wall of large arteries, typically where the wall shear rate is low and oscillatory [1].
Investigation of blood flow modeling through arterial multistenosis is very challenging. Accuracy of the simulation depends mainly on suitable numerical approach, realistic model geometry, and boundary conditions. Many investigators have focused their attention on blood flow through stenosed arteries with single stenosis by Mekheimer [2, 3], Chakravarty and Mandal [4], Lee and Xu [5], who pointed out that the mathematical model becomes more accurate in the presence of an overlapping stenosis instead of a mild one. Ang and Mazumdar [6] studied asymmetric arterial blood flow with numerical solution in three dimensions, and Ikbal et al. [7] have worked on unsteady response of non-Newtonian blood flow in magnetic field without considering periodic body acceleration. Khler et al. [8] studied the wall shear stress with the help of magnetic resonance imaging (MRI) measurements of the velocity field and compared them with simulation outputs. Stroud et al. [9] have studied a 2D plaque model using modeling and simulation while Fischer et al. [10] worked on numerical method for the computational study of arterial blood flow with turbulence. The asymmetric flows in a symmetric sudden expansion channel have been studied using experimental and numerical techniques by Fearn et al. [11] and Durst et al. [12]. Mahapatra et al. [13] investigated unsteady laminar separated flow through constricted channel using finite difference technique in staggered grid distribution and suggested that the critical value of Reynolds number depends on the area reduction and the length of the constriction. Chakravarty and Sannigrahi [14] solved blood flow model with body acceleration but they do not consider the nonlinear terms in the model. Blood shows a non-Newtonian behaviour at low shear rates in tubes of smaller diameters, and Taylor [15] suggested that at high shear rates commonly found in larger arteries blood behaves like a Newtonian fluid.
With the above motivation in our mind we have worked on numerical simulations of nonlinear pulsatile unsteady Newtonian blood flow in a rigid cylindrical tube through cosine-shape stenosis under the influences of periodic body acceleration. It appears that a few studies address the issue of that nonlinear terms present in Navier-Stokes equation govern blood flow with periodic body acceleration associated with an atherosclerotic plaque. The numerical solutions are obtained of the nonlinear model using appropriate finite difference method. A comparison of the axial and radial velocities, wall shear stresses, flux rates, and the streamlines with other existing model [7] has been studied also.
2. Physical Assumptions and Mathematical Model
The segment of the artery is modeled as an axisymmetric cylindrical tube with radius . The blood is modeled as a homogeneous incompressible viscous unsteady Newtonian fluid of density and kinematic viscosity . Therefore, the blood flow is governed by the incompressible Navier-Stokes equation.
A cylindrical coordinate system is chosen, where is the coordinate in the radial and circumferential direction, while is taken along the axis of artery as shown in Figure 1. The velocity components in the axial and radial directions are and , respectively. The flow is driven by a prescribed dimensional oscillatory axial pressure gradient given by [16]where and are the steady component of the pressure gradient and pulsatile component, respectively. The frequency , and is the heart pulse frequency. The pressure gradient in the radial direction is negligibly small as the lumen radius of artery is small compared to pressure wave length so that . Also because of the human body acceleration, the axial flow is subject to an external force . For the present model, we consider the periodic acceleration force given bywhere is the amplitude of the pulse.
According to the above assumptions, the blood flow dynamics is governed by the equation of continuitythe momentum equation in the radial direction (flow velocity )and the axial direction (flow velocity )In (4) and (5), is the Reynolds number and is the average velocity of the blood.
Finally, to model multiple axisymmetric stenosis, we define the following mathematical function :where is a dimensionless constant and the geometry of this axisymmetric stenosis in the cross section of the artery is shown in Figure 2.
We numerically simulate (3)–(5) subject to the following initial conditionand the no-slip boundary conditions
3. Numerical Simulation: Computational Method
We use the finite difference scheme to study the dynamics of blood flow through the cylindrical shape artery. To employ this method, first we transform our cylindrical domain into the rectangular domain by using the following radial transformation:Under this transformation, the equation of continuity (3) and the equations of motion in the radial direction (4) and axial direction (5), respectively, are rewritten asInitial condition (7) and no-slip boundary condition (8) due to radial transformation (9) then becomeLet us first apply the finite difference discretization scheme to solve nonlinear model equations (10)–(12). We use the central difference approximation to discretize the spatial derivatives and the explicit forward finite difference approximation to discretize the time derivative in the following manner:whereSimilarly we approximate all the partial derivatives of .
The axial velocity is obtained from (10) and (12) by applying the above finite difference scheme at any point in the domain of interest at any time with the help of the following discretize initial and boundary conditions (discretization of (13)):subject to the input pressure gradient and external force from relations (1) and (2). The radial velocity is obtained from (10) and (11).
Finally, we determine the volumetric flow rateand the wall shearing stressin the rectangular domain with the help of transformation (9), where is the viscosity. The discretize version of and is given by the following equations:
4. Simulation Results and Discussions
In this section, we shall discuss the numerical simulation of the nonlinear equations to study the influence of stenoses and body acceleration on the blood flow for different values of the physical parameters. The simulation parameters are as follows [7, 14]: , , , Re = 400, 600, and 800, , and . The results obtained for axial velocity by solving explicit finite difference scheme with various grid sizes are taken in order to achieve the convergence and stability. We perform the experiments for grid size and with and . The results are found to be very similar in both cases.
Figures 3(a) and 3(b) represent the behavior of the axial velocity profile of the blood at time without and with body acceleration, respectively. Both figures are drawn for at different . The comparative study between the figures (Figures 3(a) and 3(b)) reveals that the body acceleration enhances the axial velocity. The curves in these two figures reveal that the velocity profile is constant for and then velocity decreases and finally goes to zero on the constricted arterial wall. The values of depend on both the body acceleration and the stenosed zone. Also in the stenosed zone of the artery , the velocity is low and in presence of body acceleration, the velocity sharply decreases in this stenosed zone. Figure 4(a) shows the results for the distribution of axial velocity over the stenosed artery for three different Reynolds numbers. We can say from this figure that as enters into the stenosed zone the axial velocity starts decreasing from its maximum value in nonstenosed zone until the first constriction attaint its maximum value near and then it gradually increases up to and again it starts decreasing till the maximum height of the second stenosis and gradually starts increasing thereafter and finally again flows with maximum velocity in nonstenosis region. The three curves here indicate that the axial velocity increases in the constricted part of the artery as Reynolds number increases under the influence of body acceleration. Figure 4(b) represents the results for the distribution of radial velocity in the multiple stenosed artery for three different Reynolds numbers. From the figure we can say that the direction of radial velocity is negative in the stenosis zone due to presence of multiple stenosis. Thus the multiple stenosis and the Newtonian characteristics of the flowing blood affect the axial velocity profile which can be estimated by the relevant curves of the present figure.
(a)
(b)
(a)
(b)
The curves in Figure 5 describe the nature of radial velocity for three different Reynolds numbers. The velocity initially starts with zero and continues till and then it decreases gradually in negative direction.
To test the effects of body acceleration on axial and radial velocities profile several simulations have been carried out using the contour plot as shown in Figures 6(a)–6(d) and Figures 7(a)–7(d) respectively. The velocity profiles and are shown in different region with different colors representing the value of velocity with the help of color bar. One can see from these plot that velocity profile is divided into different layers due to the constriction of the artery and changes in the plot also can be observed in case of no body acceleration. We have shown the distribution of axial and radial velocities in Figures 8(a) and 8(d) in entire upper half segment of the artery using 3D plot. Figure 8(a) shows the axial velocity profile of the flow without body acceleration at time for . One can see that the velocity profile is constant in the nonstenosed zone and varies over the constricted area. In case of body acceleration the profile is shown in Figure 8(c).
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
Wall shear stress plays an important role in the creation and propagation of arteriosclerosis. If the wall shear stress is high then it may damage the arterial wall and is the main cause of the intimal thickening. On the other hand, the plaque formation in an artery is created in the regions of low arterial wall shear stress. Atherosclerotic lesions are associated with low and high wall shear stress. So it is important to study the wall shear stress distribution in the multistenosed artery. Figures 9(a) and 9(b) show the distribution of wall shear stress on the arterial segment for three different Reynolds numbers. The wall shear stress increases rapidly near to the peak of the constriction. Here the effects of Reynolds number can be observed from the figure. The wall shear stress increases as Reynolds number increases. Figures 10(a) and 10(b) show the distribution of flux over the stenosed artery for different Reynolds number. One can conclude that flux decreases near the picks of the stenosis.
(a)
(b)
(a)
(b)
The streamlines of the blood flow in the artery with multistenosis are found in the transformed rectangular domain with grid in the upper half zone and same grid also taken for lower half portion. We have plotted the different types of streamlines in Figures 11(a) and 11(b) at . All the streamlines follow the straight line path near the axis which gradually get perturbed more towards the wall of the stenosed artery. It is interesting to observe that several flow lines are attracted towards the stenotic wall upstream with the formation of circulation zones while others pass through the constricted region directly following the main stream.
(a)
(b)
Figure 12 shows the streamline patterns at when the artery is free from stenosis. It is observed that the lines are parallel to the axial direction.
For the purpose of model validation, the axial velocity profile is compared with [7], and as shown in Figure 13, the results are found to be in good agreement though their studies were based on the stenotic blood flow in which the streaming blood was treated as non-Newtonian fluid in magnetic field. Also the result agrees qualitatively well for Newtonian fluid with Tu et al. [17].
5. Conclusions
A nonlinear mathematical model for blood flow in a multiple stenosed arterial segment has been developed under the influence of body acceleration. The numerical simulation of blood flow is investigated in this study. As the Reynolds number increases, the wall shear stress increases. The multiple stenosis has significant effect on the wall shear stress in such a way that it develops more at the constricted locations than all other sites of the artery. The streamline pattern shows the distinct boundary layer characteristics in the arterial segment. This is validated by the flow visualisation observed by the simulation studies of [7]. The results obtained here would help researchers greatly in gaining better insight into blood flow models through the multistenosis artery under the influence of body acceleration.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.