Applied Bionics and Biomechanics
Volume 2015 (2015), Article ID 406195, 12 pages
http://dx.doi.org/10.1155/2015/406195

Research Article

Biorheological Model on Flow of Herschel-Bulkley Fluid through a Tapered Arterial Stenosis with Dilatation

S. Priyadharshini and R. Ponalagusamy

Department of Mathematics, National Institute of Technology, Tiruchirappalli, Tamilnadu 620015, India

Received 19 September 2014; Accepted 18 February 2015

Academic Editor: Cecilia Laschi

Copyright © 2015 S. Priyadharshini and R. Ponalagusamy. 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 analysis of blood flow through a tapered artery with stenosis and dilatation has been carried out where the blood is treated as incompressible Herschel-Bulkley fluid. A comparison between numerical values and analytical values of pressure gradient at the midpoint of stenotic region shows that the analytical expression for pressure gradient works well for the values of yield stress till 2.4. The wall shear stress and flow resistance increase significantly with axial distance and the increase is more in the case of converging tapered artery. A comparison study of velocity profiles, wall shear stress, and flow resistance for Newtonian, power law, Bingham-plastic, and Herschel-Bulkley fluids shows that the variation is greater for Herschel-Bulkley fluid than the other fluids. The obtained velocity profiles have been compared with the experimental data and it is observed that blood behaves like a Herschel-Bulkley fluid rather than power law, Bingham, and Newtonian fluids. It is observed that, in the case of a tapered stenosed tube, the streamline pattern follows a convex pattern when we move from to and it follows a concave pattern when we move from to . Further, it is of opposite behaviour in the case of a tapered dilatation tube which forms new information that is, for the first time, added to the literature.

1. Introduction

Blood flow through a stenosed artery is one of the important areas of research because a stenosed artery affects the entire cardiovascular system. Aortic stenosis causes chest pain and decreased blood flow to the brain resulting in loss of consciousness and heart failure which increases the risk of death. It is well known that fluid dynamical factors play a pivotal role in the formation and development of stenosis. Young [1] and Young and Tsai [2] studied the effects of stenosis on blood flow through arteries. Several investigators [3–10] analyzed the blood flow through a stenosed artery and have shown that the physical parameters affect the blood flow. Pulsatile flow of blood through a stenosed porous medium under the influence of periodic body acceleration considering blood as a Newtonian fluid has been studied by El-Shahed [11]. El-Shehawey et al. [12] have examined the pulsatile flow of blood through a tube considering blood as a Newtonian fluid taking into account the body acceleration and porosity of the tube. Sharma et al. [13] investigated the effects of radial variation of hematocrit and magnetic field on the flow of blood as a Newtonian fluid through a porous medium in a stenosed artery.

Viscoplastic materials are concentrated suspensions of solid particles or macromolecules and are classified as generalized Newtonian fluids. They flow like liquids when subjected to a stress above a critical value but respond as elastic or inelastic solids below this critical stress. According to the von Mises yield criterion, flow is assumed to occur when the second invariant of the stress exceeds the so-called yield stress [4]. It is understood that the important time-independent non-Newtonian fluid possessing a fluid behavior index (power law index) and yield values is the Herschel-Bulkley fluid, which has pivotal applications in polymer processing industries [12], developing blood oxygenators, and biomechanics [4]. Further, Herschel-Bulkley fluids include both shear thinning and shear thickening materials. The practical examples of such materials are greases, colloidal suspensions, starch pastes, tooth pastes, paints, and blood flow in an artery. These fluids have been useful as lubricant in roller bearing [13].

The non-Newtonian behavior of blood has been considered and studied by [14–17]. Chaturani and Samy [18] investigated the effects of non-Newtonian nature of blood treating it as a Casson’s fluid and pulsatility on flow through a stenosed tube. The two-dimensional flow of power law fluid in stenosed arteries has been studied and the effect of power law index on the flow separation and reattachment point has been thoroughly investigated [19]. Nadeem et al. [20] and Ismail et al. [21] have investigated blood flow through a tapered artery with a stenosis assuming the blood as a non-Newtonian power law fluid model. They analyzed the influences of different parameters (power law index, flow rate, stenosis shape, and stenosis height) in different types of tapered arteries (converging tapered, diverging tapered, and nontaperted artery). Pincombe et al. [22] proposed a fully developed one-dimensional casson flow through a stenosed artery with multiple abnormal segments. They have studied the effects of multiple stenoses and poststenotic dilatation on non-Newtonian blood flow in small arteries. Scott Blair and Spanner [23] have suggested that blood obeys Casson’s model only for moderate shear rate flows and that there is no difference between Casson’s and Herschel-Bulkley plots over the range where Casson’s plot is valid (for blood). Furthermore, Sacks et al. [24] have experimentally pointed out that blood shows the behavior characteristic of a combination of Bingham-plastic and pseudoplastic fluid-Herschel-Bulkley fluid with the fluid behavior index greater than unity. In view of the experimental observation [24] and suggestion made in [23], it is pertinent to consider the behavior of blood as a Herschel-Bulkley fluid.

The non-Newtonian aspects of blood flow through stenosed arteries have been studied by [25] treating blood as a Herschel-Bulkley fluid. Biswas and Laskar [26] have investigated the steady flow of blood as a Herschel-Bulkley fluid through a stenosed artery. In these studies, the combined effects of the rheology of blood as Herschel-Bulkley fluid model, stenosis height, dilatation depth, and tapering on the flow of blood have not been investigated. Hence, the aim of the present paper is to analyze the flow of Herschel-Bulkley fluid in a tapered artery with stenosis and dilatation (Figure 2). The expressions for velocity, wall shear stress, and flow resistance have been derived. The effects of parameters such as power law index, shear dependent nonlinear viscosity, stenotic height, taper angle, dilatation depth, and the yield stress on physiologically important quantities, namely, wall shear stress and flow resistance, are presented graphically.

2. Formulation of the Problem

Consider the steady and axially symmetric flow of an incompressible Herschel-Bulkley fluid lying in a tube having length (Figures 1 and 2). We take the cylindrical coordinate system in such a way that , , and are the velocity components in , , and directions, respectively. The equations governing the two-dimensional steady incompressible Herschel-Bulkley fluid are

Figure 1: Geometry of an axially nonsymmetrical artery with stenosis and dilatation.

Figure 2: Geometry of the tapered artery with stenosis and dilatation for different taper angle.

In the above equations the extra stress tensor for Herschel-Bulkley fluid is defined as

where is the consistency index, is the power law index (or fluid behaviour index), is the yield stress, and

where , , is the rate of strain tensor component.

We introduce the nondimensional variables

where is the average velocity of flow of Newtonian fluid, , , is the radius of the normal artery, Re is the Reynolds number, is the radius of the abnormal artery, and is the viscosity of Newtonian fluid. By assuming

the cases of mild stenosis and mild dilatation , (1) with the help of (2) and (4) take the form

where . The corresponding boundary conditions are

The equations describing the geometry of the wall are

where is the maximum distance the th abnormal segment projects into the lumen and is negative for aneurysms and positive for stenosis, is the radius of the artery, and , where is the taper angle. For converging tapering becomes greater than 0, indicates the diverging tapering, and for the case of nontapered artery, is the length of the th abnormal segment, denotes the distance from the origin to the commencement of the th abnormal segment and is given by

indicates the distance between the origin of the flow region and the end of the th abnormal segment and is given by

and is the distance separating the start of the th abnormal segment from the end of the th or from the start of the segment if [22].

3. Solution of the Problem

The exact solution for velocity field satisfying the boundary conditions can be written as

where .

The plug core velocity is given by

where is the radius of the plug core region and .

Multiplying (11) by and integrating with respect to , the stream function is obtained as

The volumetric flow rate is defined as

The total flow rate is defined as

Using (11), (12), and (15), we get

The shear stress at the wall of the tapered arterial stenosis with dilatation (wall shear stress ) is defined as

The flow resistance is defined as

where is any point of cross section of nonuniform tube along the axial direction.

Case 1. For any value of yield stress , (16) can be rewritten as

where . For , one can numerically compute the value of (pressure gradient) from (19) for different values of the parameters. Equation (19) has been numerically solved for using Newton-Raphson method.

Case 2. For small value of yield stress , the expression for pressure gradient can be obtained as

Using (17) and (20), the wall shear stress is obtained as

Substituting (20) into (18), the analytical expression for flow resistance is obtained as

where

Considering the number of abnormal segments within an arterial segment as shown in Figure 1, we define as the starting point and as the ending point of each portion. Taking this into account (22) can be rewritten as

4. Discussion

A comparison between numerical values and analytical values of pressure gradient at the midpoint of stenotic region shows that, up to , the maximum error is less than 1.4% and for dilatation region the maximum error is less than 6%. This is illustrated in Tables 1 and 2. This implies that the analytical expression for pressure gradient works well for the values of yield stress till 2.4.

Table 1: A comparison between numerical value of pressure gradient and approximate value of pressure gradient at the midpoint of stenotic region for different values of yield stress taking = 0.8, = 1.2, , and .

Numerical value of
17.9423
18.1099
19.1145
20.4511
21.7841
23.1127
24.4368
25.7559

Analytical value of
17.9424
18.1104
19.124
20.4895
21.8712
23.2692
24.6834
26.1137

Table 2: A comparison between numerical value of pressure gradient and approximate value of pressure gradient at the midpoint of dilatation region for different values of yield stress taking = 0.8, = 1.2, , and .

Numerical value of
5.0954
5.2106
5.9001
6.8133
7.7185
8.6152
9.5037
10.3845

Analytical value of
5.0956
5.2116
5.9163
6.8799
7.8709
8.8892
9.9349
11.008

A comparative study of velocity profiles for fluids such as Newtonian, power law, Bingham-plastic, and Herschel-Bulkley fluids is represented graphically in Figure 3. From Figure 3, it is observed that the velocity of Herschel-Bulkley fluid agrees with experimental values compared to that of the other fluids.

Figure 3: Comparison of velocity profiles for various fluids with experimental results.

The variation of wall shear stress (WSS) with respect to axial distance for the case of a converging tapered, not tapered, and diverging tapered arterial stenosis with dilatation is displayed in Figures 4–8. WSSs of fluids such as Newtonian, power law, Bingham-plastic, and Herschel-Bulkley fluids are compared in Figure 4. It is important to note that WSS increases in the upstream of the stenotic region ( to 2.5), reaches maximum at the midpoint (), and decreases in the downstream of region ( to 3), while, in the dilatation region, WSS decreases as varies from 4 to 4.5, reaches minimum at the midpoint (), and increases in the region ( to 5). In the case of stenosis, increase is more for converging tapered artery () as compared to the case of not tapered () and diverging tapered () artery. It is observed from the view of variation of WSS around the midpoint of stenotic region that the effect of the presence of stenosis is higher on the rheology of blood as Bingham fluid model in comparison with the rheology of blood as Newtonian, Hershel-Bulkley, and power law fluid models, respectively. It is important to observe from Figure 5 that the power law index () plays a significant role in stenotic region ( to 3) since the percentage of variation in WSS is higher for stenosis as compared to the case of dilatation.

Figure 4: Axial variation of wall shear stress for Newtonian, power law, Bingham-plastic, and Herschel-Bulkley fluids with different values of .

Figure 5: Axial variation of wall shear stress () for different values of power law index () taking , , and .

Figure 6: Axial variation of wall shear stress () for different values of consistency index () taking and .

Figure 7: Variation of wall shear stress () with respect to axial distance for different values of taking , , and .

Figure 8: Axial variation of wall shear stress () for different values of by taking , , and .

Axial variation of WSS with respect to yield stress in the case of converging tapered, not tapered, and diverging tapered arterial stenosis with dilatation is displayed in Figure 7. Increase in yield stress causes wall shear stress to increase and the variation is more in the stenotic region than in the dilatation region. The effect of stenotic height on WSS has been investigated in Figure 8. As stenotic height increases, WSS increases in the stenotic region while it decreases in the dilatation region. When there is no stenosis, WSS increases linearly with respect to the axial distance. It is observed that the stenotic height plays a predominant role in increasing the WSS. The variation is more in the case of converging tapered than not tapered and diverging tapered arteries.

Figures 9–13 are prepared to see the variation of resistance to flow with respect to the axial distance in the case of converging tapered, not tapered, and diverging tapered arterial stenosis with dilatation. A comparative study of flow resistance for Newtonian, power law, Bingham-plastic, and Herschel-Bulkley fluids is depicted in Figure 9. Flow resistance increases significantly in the stenotic region ( to 3): the increase is more for Herschel-Bulkley fluid and comparatively less for Newtonian fluid. Flow resistance decreases with the axial distance in the dilatation region. The variation of flow resistance for power law and Bingham-plastic is lesser when compared with Herschel-Bulkley and greater when compared with Newtonian fluid. Figures 10 and 11 depict that flow resistance increases as and increases. Increase in power law index causes flow resistance to increase significantly as compared to the consistency index (). The effect of yield stress on flow resistance having other parameters fixed has been studied from Figure 12. Flow resistance increases as yield stress increases and the variation caused by yield stress is less compared to other parameters. Figure 13 shows that the flow resistance increases with stenotic height and its increase is more in the case of converging tapered artery.

Figure 9: Variation of flow resistance () for various fluids with respect to axial distance taking .

Figure 10: Axial variation of flow resistance () for different values of taking , , and .

Figure 11: Axial variation of flow resistance () for different values of taking , , and .

Figure 12: Axial variation of flow resistance () for different values of yield stress () taking , , and .

Figure 13: Variation of flow resistance () with respect to axial distance for different values of taking , , and .

The effects of consistency index (), power law index (), and yield stress on the stream line pattern have been examined and illustrated in Figures 14–19. In the case of tapered stenosed tube (Figures 14–16), the non-Newtonian behaviour of blood plays a predominant role in the formation of trapping bolus. Increase in or does not cause a significant change in the stream line pattern. Increase in yield stress leads to a significant increase in the size of trapping bolus. It is observed that the parameters and are weak parameters in the sense that these parameters bring a small change in the stream line pattern in comparison with the yield stress. Figures 14–19 reveal that, in the case of a tapered stenosed tube, the stream line pattern follows a convex pattern when we move from to and it follows a concave pattern when we move from to . Further, it is of opposite behaviour in the case of a tapered dilatation tube. In the case of dilatation, the variation in the stream line pattern corresponding to change in parameters is less due to lower pressure gradient. This has been illustrated in Figures 17–19.

(a)

(b)

(c)

Figure 14: Stream lines for (a) Newtonian fluid, (b) Bingham-plastic fluid (), and (c) Bingham-plastic fluid () taking and in the stenotic region.

(a)

(b)

(c)

Figure 15: Stream lines for (a) power law fluid, (b) Herschel-Bulkley fluid (), and (c) Herschel-Bulkley fluid () taking and in the stenotic region.

(a)

(b)

Figure 16: Stream lines for different values of : (a) and (b) taking , , , , and in the stenotic region.

(a)

(b)

Figure 17: Stream lines for (a) Newtonian fluid and (b) Bingham-plastic fluid taking and in the dilatation region.

(a)

(b)

Figure 18: Stream lines for (a) power law fluid and (b) Herschel-Bulkley fluid taking and in the dilatation region.

(a)

(b)

Figure 19: Stream lines for different values of : (a) and (b) taking , , , and in the dilatation region.

5. Conclusion

This work presents a model of flow of an incompressible Herschel-Bulkley fluid through a tapered artery with stenosis and dilatation. In this paper, we conclude the following.

(i)Expressions for velocity profile, wall shear stress, and flow resistance are derived.(ii)A comparison between numerical values and analytical values of pressure gradient at the midpoint of stenotic region shows that, up to , the maximum error is less than 1.4% and, for dilatation region, the maximum error is less than 6%. This implies that the analytical expression for pressure gradient works well for the values of yield stress till 2.4.(iii)Effects of parameters such as power law index, consistency index, yield stress, stenotic height, dilatation depth, and taper angle on the above mentioned physiologically important quantities are studied.(iv)For given value of power law index (), Herschel-Bulkley fluid has greater wall shear stress than the power law fluid.(v)It is important to note that increase in yield stress leads to increase in wall shear stress and resistance to flow.(vi)Flow resistance increases significantly as the stenotic height increases for given .(vii)It is observed that the parameters and are weak parameters in the sense that these parameters bring a small change in the stream line pattern in comparison with the yield stress and the stream line pattern for tapered dilatation tube is of opposite behaviour as compared to tapered stenosed tube.(viii)From the present work, results for power law (taking ), Bingham-plastic (taking ), and Newtonian fluids (taking and ) can be obtained.(ix)Results illustrated through graphs show the effects of multiple diseased portions of artery in close proximity to each other (a poststenotic dilatation) on the increase of flow resistance causing the reduction of blood flow.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

D. F. Young, “Effect of a time dependent stenosis of flow through a tube,” Journal of Engineering for Industry, vol. 90, pp. 248–254, 1968.

D. F. Young and F. Y. Tsai, “Flow characteristics in models of arterial stenoses: I. Steady flow,” Journal of Biomechanics, vol. 6, no. 4, pp. 395–410, 1973.

D. F. Young, “Fluid mechanics of arterial stenoses,” Journal of Biomechanical Engineering, vol. 101, no. 3, pp. 157–175, 1979.

R. Ponalagusamy, Blood flow through stenosed tube [Ph.D. thesis], IIT, Bombay, India, 1986.

D. B. Clegg and G. Power, “Flow of a Bingham fluid in a slightly curved tube,” Applied Scientific Research: Section A, vol. 12, no. 2, pp. 199–212, 1963.

J. A. Greenwood and J. J. Kauzlarich, “EHD lubrication with Herschel-Bulkley model greases,” ASLE Transactions, vol. 4, pp. 269–278, 1972.

P. Chaturani and R. Ponnalagarsamy, “Analysis of pulsatile blood flow through stenosed arteries and its applications to cardiovascular diseases,” in Proceedings of the 3rd National Conference on Fluid Mechanics and Fluid Power, pp. 463–468, 1984.

R. Ponalagusamy, R. T. Selvi, and A. K. Banerjee, “Mathematical model of pulsatile flow of non-Newtonian fluid in tubes of varying cross-sections and its implications to blood flow,” Journal of the Franklin Institute, vol. 349, no. 5, pp. 1681–1698, 2012.

J. C. Misra, M. K. Patra, and S. C. Misra, “A non-Newtonian fluid model for blood flow through arteries under stenotic conditions,” Journal of Biomechanics, vol. 26, no. 9, pp. 1129–1141, 1993.

V. P. Srivastava, “Arterial blood flow through a non symmetric stenosis with applications,” Japanese Journal of Applied Physics, vol. 34, pp. 6539–6545, 1995.

M. El-Shahed, “Pulsatile flow of blood through a stenosed porous medium under periodic body acceleration,” Applied Mathematics and Computation, vol. 138, no. 2-3, pp. 479–488, 2003.

E. F. Elshehawey, E. M. Elbarbary, N. A. S. Afifi, and M. El-Shahed, “Pulsatile flow of blood through a porous medium under periodic body acceleration,” International Journal of Theoretical Physics, vol. 39, no. 1, pp. 183–188, 2000.

M. K. Sharma, K. Bansal, and S. Bansal, “Pulsatile unsteady flow of blood through porous medium in a stenotic artery under the influence of transverse magnetic field,” Korea Australia Rheology Journal, vol. 24, no. 3, pp. 181–189, 2012.

P. K. Mandal, “An unsteady analysis of non-Newtonian blood flow through tapered arteries with a stenosis,” International Journal of Non-Linear Mechanics, vol. 40, no. 1, pp. 151–164, 2005.

P. K. Mandal, S. Chakravarty, A. Mandal, and N. Amin, “Effect of body acceleration on unsteady pulsatile flow of non-Newtonian fluid through a stenosed artery,” Applied Mathematics and Computation, vol. 189, no. 1, pp. 766–779, 2007.

R. Ponalagusamy, “Mathematical analysis on effect of non-Newtonian behavior of blood on optimal geometry of microvascular bifurcation system,” Journal of the Franklin Institute, vol. 349, no. 9, pp. 2861–2874, 2012.

R. Ponalagusamy, “Pulsatile flow of Herschel-Bulkley fluid in tapered blood vessels,” in Proceedings of the International Conference on Scientific Computing (CSC '13), and World Congress in Computer Science, Computer Engineering, and Applied Computing (WORLDCOMP '13), pp. 67–73, Las Vegas, Nev, USA, July 2013.

P. Chaturani and R. Ponalagusamy, “Pulsatile flow of Casson's fluid through stenosed arteries with applications to blood flow,” Biorheology, vol. 23, no. 5, pp. 499–511, 1986.

P. Chaturani and R. Ponalagarsamy, “Dilatency effects of blood on flow through arterial stenosis,” in Proceedings of the 28th Congress of the Indian Society of Theoretical and Applied Mechanics, pp. 87–96, 1983.

S. Nadeem, N. S. Akbar, A. A. Hendi, and T. Hayat, “Power law fluid model for blood flow through a tapered artery with a stenosis,” Applied Mathematics and Computation, vol. 217, no. 17, pp. 7108–7116, 2011.

Z. Ismail, I. Abdullah, N. Mustapha, and N. Amin, “A power-law model of blood flow through a tapered overlapping stenosed artery,” Applied Mathematics and Computation, vol. 195, no. 2, pp. 669–680, 2008.

B. Pincombe, J. Mazumdar, and I. Hamilton-Craig, “Effects of multiple stenoses and post-stenotic dilatation on non-Newtonian blood flow in small arteries,” Medical and Biological Engineering and Computing, vol. 37, no. 5, pp. 595–599, 1999.

G. W. Scott Blair and D. C. Spanner, An Introduction to Biorheology, Elsevier Scientific Publishing, Amsterdam, The Netherlands, 1974.

A. H. Sacks, K. R. Raman, J. A. Burnell, and E. G. Tickner, “Auscultatory versus direct pressure measurements for newtonian fluids and for blood in simulated arteries,” VIDYA Report 119, 1963.

P. Chaturani and R. Ponnalagarsamy, “A study of non Newtonian aspects of blood flow through stenosed arteries and its applications in arterial diseases,” Biorheology, vol. 22, no. 6, pp. 521–531, 1985.

D. Biswas and R. B. Laskar, “Steady flow of blood through a stenosed artery—a non-Newtonian fluid model,” Assam University Journal of Science and Technology, vol. 7, pp. 144–153, 2011.