International Scholarly Research Notices
Volume 2015 (2015), Article ID 628605, 10 pages
http://dx.doi.org/10.1155/2015/628605

Research Article

Numerical Simulation of Nonlinear Pulsatile Newtonian Blood Flow through a Multiple Stenosed Artery

Satyasaran Changdar¹ and Soumen De²

¹Institute of Engineering & Management, Saltlake, Kolkata 700101, India
²Department of Applied Mathematics, University of Calcutta, 92 Acharya Prafulla Chandra Road, Kolkata 700009, India

Received 16 July 2015; Accepted 13 October 2015

Academic Editor: Abdelouahed Tounsi

Copyright © 2015 Satyasaran Changdar and Soumen De. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 by

where is the amplitude of the pulse.

Figure 1: Cylindrical coordinate system with multiple stenosis along the axial direction, where is the length of the stenosis and is distance of stenosis from the radial axis.

According to the above assumptions, the blood flow dynamics is governed by the equation of continuity

the 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.

Figure 2: The schematic representation of the stenosis along the axial direction as given by . The geometrical parameters are as follows: , , and .

We numerically simulate (3)–(5) subject to the following initial condition

and 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 as

Initial condition (7) and no-slip boundary condition (8) due to radial transformation (9) then become

Let 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:

where

Similarly 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 rate

and the wall shearing stress

in 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)

Figure 3: Distribution of axial velocities at different for . Axial velocity (a) without body acceleration and (b) with body acceleration.

(a)

(b)

Figure 4: Distribution of (a) axial velocities for various Reynolds numbers and (b) radial velocities for different Reynolds number.

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.

Figure 5: The distribution of radial velocity at .

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)

Figure 6: Contour distribution of (a) axial velocities for Reynolds number with body acceleration, (b) axial velocities for without body acceleration, (c) axial velocities for Reynolds number with body acceleration and (d) axial velocities for without body acceleration.

(a)

(b)

(c)

(d)

Figure 7: Contour distribution of (a) radial velocities for Reynolds number with body acceleration, (b) radial velocities for without body acceleration, (c) radial velocities for Reynolds number with body acceleration, and (d) radial velocities for without body acceleration.

(a)

(b)

(c)

(d)

Figure 8: 3D plot of (a) axial velocities without body acceleration at , (b) radial velocities without body acceleration at , (c) axial velocities with body acceleration at , and (d) radial velocities with body acceleration at .

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)

Figure 9: Distribution of (a) wall shear stress with body acceleration and (b) wall shear stress without body acceleration.

(a)

(b)

Figure 10: Distribution of (a) flux with body acceleration and (b) flux without body acceleration.

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 11: Distribution of (a) streamlines with body acceleration and (b) streamlines without body acceleration at time in the upper half segment of the artery.

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.

Figure 12: The streamline patterns of the flow at when there is no constriction.

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].

Figure 13: Validation of axial velocity at .

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.

References

D. M. Wooton and D. N. Ku, “Fluid mechanics of vascular systems, diseases, and thrombosis,” Annual Review of Biomedical Engineering, vol. 1, pp. 299–329, 1999.

K. S. Mekheimer and M. A. El Kot, “Suspension model for blood flow through catheterized curved artery with time-variant overlapping stenosis,” Engineering Science and Technology, vol. 18, no. 3, pp. 452–462, 2015.

K. S. Mekheimer, F. Salama, and M. A. El Kot, “The unsteady flow of a Carreau fluid through inclined catheterized arteries have a balloon (angioplasty) with time-variant overlapping stenosis,” Walailak Journal of Science and Technology, vol. 12, pp. 863–883, 2014.

S. Chakravarty and P. K. Mandal, “Mathematical modelling of blood flow through an overlapping arterial stenosis,” Mathematical and Computer Modelling, vol. 19, no. 1, pp. 59–70, 1994.

K. W. Lee and X. Y. Xu, “Modelling of flow and wall behaviour in a mildly stenosed tube,” Medical Engineering and Physics, vol. 24, no. 9, pp. 575–586, 2002.

K. C. Ang and J. N. Mazumdar, “Mathematical modelling of three-dimensional flow through an asymmetric arterial stenosis,” Mathematical and Computer Modelling, vol. 25, no. 1, pp. 19–29, 1997.

M. A. Ikbal, S. Chakravarty, K. K. L. Wong, J. Mazumdar, and P. K. Mandal, “Unsteady response of non-Newtonian blood flow through a stenosed artery in magnetic field,” Journal of Computational and Applied Mathematics, vol. 230, no. 1, pp. 243–259, 2009.

U. Khler, I. Marshall, M. B. Robertson, Q. Long, X. Y. Xu, and P. R. Hoskins, “MRI measurement of wall shear stress vectors in bifurcation models and comparison with CFD predictions,” Journal of Magnetic Resonance Imaging, vol. 14, no. 5, pp. 563–573, 2001.

J. S. Stroud, S. A. Berger, and D. Saloner, “Numerical analysis of flow through a severely stenotic carotid artery bifurcation,” Journal of Biomechanical Engineering, vol. 124, no. 1, pp. 9–20, 2002.

P. F. Fischer, F. Loth, S. E. Lee, S.-W. Lee, D. S. Smith, and H. S. Bassiouny, “Simulation of high Reynolds number vascular flows,” Computer Methods in Applied Mechanics and Engineering, vol. 196, no. 31-32, pp. 3049–3060, 2007.

R. M. Fearn, T. Mullin, and K. A. Cliffe, “Nonlinear flow phenomena in a symmetric sudden expansion,” Journal of Fluid Mechanics, vol. 211, pp. 595–608, 1990.

F. Durst, J. C. F. Pereira, and C. Tropea, “The plane symmetric sudden-expansion flow at low Reynolds numbers,” Journal of Fluid Mechanics, vol. 248, pp. 567–581, 1993.

T. R. Mahapatra, G. C. Layek, and M. K. Maiti, “Unsteady laminar separated flow through constricted channel,” International Journal of Non-Linear Mechanics, vol. 37, no. 2, pp. 171–186, 2002.

S. Chakravarty and A. K. Sannigrahi, “Effect of body acceleration on blood flow in an irregular stenosed artery,” Mathematical and Computer Modelling, vol. 19, no. 5, pp. 93–103, 1994.

M. G. Taylor, “The influence of the anomalous viscosity of blood upon its oscillatory flow,” Physics in Medicine and Biology, vol. 3, no. 3, pp. 273–290, 1959.

A. C. Burton, Physiology and Biophysics of the Circulation, Introductory Text, Year Book Medical Publishers, Chicago, Ill, USA, 1966.

C. Tu, M. Deville, L. Dheur, and L. Vanderschuren, “Finite element simulation of pulsatile flow through arterial stenosis,” Journal of Biomechanics, vol. 25, no. 10, pp. 1141–1152, 1992.