A Study for Analyzing the Effect of Overlapping Stenosis and Dilatation on Non-Newtonian Blood Flow in an Inclined Artery

Show more

1. Introduction

In general, the non-Newtonian effects in blood are depending on the magnitude of deformation rates; consequently non-Newtonian effects may exist or be enhanced at low shear rates flow regimes. The types of deformation, shear or elongation also influence the non-Newtonian effects of blood [1] . For a given rheological model [2] , the flow resistance and wall shear stress increase as the size of stenosis increases while the non- Newtonian nature of blood reduces the resistance and stress; consequently it contributes to the body protection. In this context, shear thinning seems to have the most significant role in facilitating blood flow through stenotic vessels. Under steady flow conditions blood is a predominantly shear thinning fluid that exhibits important non- Newtonian impact that was discussed by [3] [4] [5] . Blood also demonstrates yield stress which has been an issue of debates [6] [7] . A diseased artery is result of stenosis or dilatation. The word stenosis signifies the contraction of blood vessel due to deposition of fatty acids or plaque on the wall of artery. It is an abnormal and unnatural growth that develops at various locations of the cardiovascular system under diseased conditions. Through various means it is proved and well established that the constriction contributes to significant changes in the blood flow, wall shear stress, flow rate and the impedance (flow resistance) from its normal behavior. The particular cause of the occurrence of permanent dilatation is not clear. It is considered that besides complex interaction among several biological factors, some precise alterations in the hemodynamic stimuli on the arterial wall may be responsible for dilatation. After the formation of aneurysm a gradual expansion might persist on the weakened arterial wall due to the hemodynamic forces exerted by the pulsatile flow of blood. The specific alterations in the hemodynamic forces exerted on the arterial wall may play a key role to the origin and progression of the disease is widely considered hypothesis and many studies, specifically, [8] - [14] were carried out. Realizing the importance of the role of geometrical interaction with flow, several studies based on fluid dynamics through stenosed artery have been carried out to evaluate the flow pattern, wall shear stress, flow rate and impedance by many physicists and mathematicians. Newtonian blood flow in a tapered artery was studied in [15] and observed that the problem becomes more serious in the presence of an overlapping stenosis in the artery and effect of vessel tapering is another important factor. Particulate suspension blood flow through stenotic arteries, to analyze the effects of hematocrit and stenosis shape was considered in [16] and observed that the flow resistance decreases with increasing shape parameter but increases with hematocrit. Herschel-Bulkley fluid flow through an inclined tube with two mild stenoses was studied by [17] and observed that the flow resistance increases with the thickness of stenoses, yield stress and power law index. Blood flow through an overlapping constriction in arteries has been investigated by [18] and the expression for the impedance, the wall shear stress, the shear stress at the stenosis throats and at critical height of the stenosis has been derived. Shear stress at the stenosis throats and at critical height of the stenosis has been derived. Flow of blood in an inclined artery which is radially symmetric and axially non-symmetrical was discussed in [19] . The pulsatile flow through a stenosed artery filled with porous medium in the presence of transverse static magnetic field under the consideration of hematocrit dependent viscosity of blood was discussed by [20] . The pulsatile flow of blood through partially occluded elastic arteries was discussed in [21] and observed that stenosis plays a critical role in plaque growth. The blood flow through tapered elastic artery with overlapping stenosis was studied by [22] and computed the expressions of impedance resistance to flow and wall shear stress. Therefore, the knowledge of the flow field in a stenosed tube may help in proper understanding and curative measures of arterial diseases. This can be done through experimental as well as analytical studies. The simulations from the analytical studies may provide many outcomes in-vitro for understanding of the interconnection of hemodynamic parameters. The aim of the present study is to investigate the effects of an overlapping stenoses and a dilatation on the flow properties and hemodynamic indicators for flow of non-Newtonian blood that described by Herschel-Bulkley model.

2. Formulation of the Problem

It is a common clinical finding that the dilatation in the arterial wall arises distal to the constriction [23] . Also, [24] mentioned that the multiple stenoses are common in the arteries as are post-stenotic dilatations of the artery.

In the present study blood is considered as non-Newtonian in nature described by the Herschel-Bulkley model that flowing in a rigid walled unhealthy artery suffering with two stenoses which are overlapped alongside and a dilatation distal to the stenotic region. The axis of the circular tube is taken along the z-axis and the stenoses and dilatation are assumed to be symmetrical about the z-axis. It is assumed that the tube is inclined at an angle α to the horizontal as shown in Figure 1. The blood vessel geometry is determined by the radius R of the inlet and outlet non-constricted segment, whereas the radius of the smooth axisymmetric constricted and dilatation segment is given by

Figure 1. Physical model of the problem.

(1)

where, 6L is the length of stenosed segment and 2L is the length of dilatation segment. The δ_{1} is maximum thickness of the stenosis and δ_{2} is maximum dilatation in the arterial wall. In the axisymmetric coordinates system, the axis of the vessel is coincides with the z*-axis.

The function h(z) is plotted on MatLab to draw the Figure 1.

The equation of continuity for incompressible fluid is given by

. (2)

Since, flow is axisymmetric, i.e., , and there is no radial flow, i.e.,. Also, flow is developed, i.e.,.

The governing equation of motion for steady developed flow of blood in the inclined vessel is given by

(3)

where, g is acceleration due to gravity and is shear stress.

The Herschel-Bulkley fluid is defined by the constitutive equation

(4)

(5)

where, is yield stress, n and k are parameters which represents non-Newtonian effects.

The boundary conditions are given by

(6)

. (7)

3. Method of Solution

Using the following non-dimensional parameters:

, , , , , , , , , into the Equations (1) to (7), the non-di-

mensional form of these equations are given by (8) to (14), respectively

(8)

. (9)

The equation of motion is

(10)

where, is dimensionless shear stress for Herschel-Bulkley fluid and given by:

(11)

. (12)

Corresponding boundary conditions are reduced to

(13)

. (14)

Equation (10) can be written as

(15)

where,.

The solution of (15) under the boundary condition (13) gives

(16)

Using (16) in the constitutive Equation (11) for provides

. (17)

Integrating Equation (17) w.r.t. r and applying boundary condition (14), the is given by

(18)

3.1. Skin Friction ()

The skin friction from Equation (16) is given by

. (19)

3.2. Volumetric Flow Rate (Q)

(20)

Substituting the value of from Equations (11) into (20), we obtain

(21)

when, then Equation (21) reduces to

. (22)

From the Equation (22), is given by

. (23)

Using Equation (19) in Equation (23)

. (24)

3.3. Pressure Drop ()

The pressure drop across the stenosis can be obtained by integrating the Equation (24) along the length of the artery from −5 L to 10 L,

(25)

In the absence of any constriction, the pressure drop is given by

. (26)

3.4. Impedance (Resistance to Flow Ratio λ)

The resistance to flow in the diseased artery is given by

(27)

If there is no diseased portion then resistance to flow is given by

(28)

Therefore, the impedance-resistance to flow ratio may be expressed as

(29)

From Equation (27) and (28), is expressed as

(30)

The values of are computed numerically using Simpson’s one-third rule and presented through graphs.

4. Results and Discussion

The effects of various parameters pertinent to the model on velocity profile, skin-fric- tion and impedance are computed and plotted to analyze them.

Figure 2 demonstrates the effect of varying thickness of the overlapped stenosed segment keeping other parameters fixed. It is observed that the axial flow velocity enhances with the increase in the thickness of the stenosis. As per mathematical formulation of the geometry, it can be seen that in the overlapped portion of the stenoses

Figure 2. Variation in velocity profile in radial direction at k = 4, n = 1.5, α = π/6, τ_{0} = 0.5, δ_{2} = 10% of R.

having peak at the location corresponding to which and at, , therefore the projection of stenosis in the lumen at is more

as compared to that at. It is observed that flow velocity at is more as compared to that at. Also, in the dilatation region the effect of is not observed. The plug flow radius decreases significantly with the increase of maximum thickness of stenosis.

Figure 3 demonstrates the effect of varying thickness of dilatation at different locations starting from the pre constricted region i.e. at location. It is observed that the axial velocity at the locations in the stenosed region located at and

, the respective passage is and, the flow velocity sig-

nificantly increases with the reduction of the flow passage. It is plausible that the flow velocity reduces from the velocity in the normal region of the vessel. In the dilatation region, the flow velocity reduces significantly in magnitude as compared to the other locations in the region of overlapped stenosis or normal region of vessel.

The flow velocity retarded with the increase in the maximum dilatation in the vessel wall. The results in the region of stenosis and effect of maximum thickness of stenosis on the flow profile pattern are in good agreement with results already published. The effect of non-Newtonian behaviour of blood as compare to blood of Newtonian nature is plausible as the profiles are not of parabolic shape.

Figure 3. Variation in velocity profile in radial direction at k = 4, n = 1.5, α = π/6, τ_{0} = 0.5, δ_{1} = 20% of R.

In the healthy human blood, the value of yield stress ranges in between 0.01 dyne/s^{2} and 0.06 dyne/s^{2}, while in unhealthy human, suffering with myocardial infarction the yield stress may increases up to five times [25] .

The measurement of yield stresses may used in the diagnosis of pathological conditions such as diabetes [26] , systemic sclerosis [27] etc. With pressure-scanning capillary hemorheometer [28] the yield stress for the Herschel-Bulkley model was found 0.0325 dyne/s^{2}. In the present study the values of yield stress are taken 0.05, 0.025 and 0.50. The last value correlates the diseased situation.

Figure 4 shows comparative study for velocity profile in respect of yield stress of healthy and diseased arterial blood. In unhealthy situation at the yield stress τ_{0} = 0.5, the plug flow radius increases irrespective of the locations whether it is falls within region of the stenosis or in the region of dilatation.

Figure 5 depict the effect of power law index on the flow profiles. It is observed that velocity profiles enhances in all region of the vessel. This enhancement in the profiles is more in the region having lesser flow passage. The plug flow radius become widen with the increase of power law index. The flow profiles at greater value of power law index are tends to be curved in the vicinity of the axial line while for small value of power law index profiles are flatten. This shows that with increasing power law index the fluid get impetus even in the plug flow region.

Figure 4. Variation in velocity profile in radial direction at k = 4, n = 1.5, α = π/6, δ_{1} = 10%, δ_{2} = 10% of R.

Figure 5. Variation in velocity profile in radial direction at k = 4, n = 1.5, α = π/6, τ_{0} = 0.5, δ_{1} = 20% of R, δ_{2} = 10% of R.

Figure 6 shows the effect of degree of stenosis on the skin-friction along the axial direction. The skin-friction is more for lower value of maximum thickness of overlapped stenoses. Also, the skin-friction sharply increases in the region of growing thickness of stenosis up to apex of maximum thickness and reduces in the overlapped stenotic region and then again rises up to second peak of the overlapped region.

Figure 7 depict that in the region of dilatation, the skin-friction reduces with the increase of maximum expansion of dilatation in the artery. This reduction in the skin- friction with the increase in expansion of dilatation is lesser than the rise in the stenotic region.

Figure 8 demonstrates the effect of yield stress on the skin-friction. With the increase in the value of yield stress, there is significant increase in the skin-friction in the entire region. Near the extremities of the overlapped stenosis, the steepest enhancement in the skin-friction is observed. This outcome suggests that in case of unhealthy situation at the yield stress τ_{0} = 0.5 arterial wall bears more stresses which may responsible for the wall rupturing.

The typical values of the power index n for blood flow is taken to lie between 0.9 and 1.1 are used in [29] .

Figure 9 depicts the effect of power law index on the skin-friction. With the increase in power law index profound reduction in skin-friction is obtained. This proves the importance of power law index of the blood during the theoretical studies. The increase

Figure 6. Variation in skin-friction in axial direction at k = 4, n = 1.5, α = π/6, τ_{0} = 0.5, δ_{2} = 10% of R.

Figure 7. Variation in skin-friction in axial direction at k = 4, n = 1.5, α = π/6, τ_{0} = 0.5, δ_{1} = 20% of R.

Figure 8. Variation in skin-friction in axial direction at k = 4, n = 1.5, α = π/6, δ_{1} = 20% of R, δ_{2} = 10% of R.

Figure 9. Variation in skin-friction in axial direction at k = 4, τ_{0} = 0.5, α = π/6, δ_{1} = 20% of R, δ_{2} = 10% of R.

in power law index overcomes the stress on wall of the vessel consequently on the endothelial cells.

Corresponding to a fixed thickness of stenosis, the increase in maximum expansion in dilatation of the vessel, the impedance to the flow is reduces as observed in Figure 10. Also, it is plausible from Figure 10 that for a fixed value of maximum expansion in the dilatation, the impedance to the flow rises on increasing maximum thickness of the pre stenotic region. Findings of this study are in good agreement with outcomes of [30] .

Figure 11 & Figure 12 show that the increase in yield stresses increases impedance to the flow. This suggests that in case of unhealthy situation, at the yield stress τ_{0} = 0.5, the resistance to flow ratio increases hence the reduction in the transportation of blood, results in the less oxygen and nutrients supply to the tissues. It is also evident from the Figure 11 that when projection of stenosis increases in the lumen the impedance increases while from Figure 12 it is seen that when the expansion of the dilatation increases the impedance decreases.

Figure 13 shows that corresponding to a fixed expansion in dilatation; the increase in maximum thickness of stenotic region has enhanced impedance to the flow. Also, it is plausible from Figure 13 that for a fixed value of maximum thickness of the stenotic region, the impedance to the flow lowered down on increasing maximum expansion of the dilatation.

The impedance to the flow reduces with the increase of power law index as shown in Figure 14 & Figure 15. It is noted that in both the cases of n < 1 or n > 1 impedance

Figure 10. Variation in impedance with thickness of stenosis at k = 4, n = 1.5, α = π/6, τ_{0} = 0.5.

Figure 11. Variation in impedance with thickness of stenosis at k = 4, n = 1.5, α = π/6, δ_{2} = 10% of R.

Figure 12. Variation in impedance with thickness of dilatation at k = 4, n = 1.5, α = π/6, δ_{1} = 20% of R.

Figure 13. Variation in impedance with thickness of dilatation at k = 4, n = 1.5, α = π/6, τ_{0} = 0.5.

Figure 14. Variation in impedance with thickness of stenosis at k = 4, α = π/6, τ_{0} = 0.5, δ_{2} = 10% of R.

Figure 15. Variation in impedance with thickness of dilatation at k = 4, α = π/6, δ_{1} = 20% of R.

reduces with the increase of n since the pressure drop increases with the increasing value of n. Also, the impedance enhances on the advancement of thickness of stenotic region while it reduces with the increase in maximum expansion of the dilatation as observed in Figure 14 & Figure 15 respectively.

The inclination α of the vessel affects significantly on the impedance to the flow, Figure 16 & Figure 17. The impedance to flow is increases with the increase in the

Figure 16. Variation in impedance with thickness of stenosis at k = 4, n = 1.5, τ_{0} = 0.5, δ_{2} = 10% of R.

Figure 17. Variation in impedance with thickness of dilatation at k = 4, n = 1.5, τ_{0} = 0.5, δ_{1} = 20% of R.

value of α. At α = 0 the vessel is horizontal then the effect of gravitational field is uni-

form throughout. The case when α = π/6 and α = π/3 the body force term

play role to retard the flow.

5. Conclusions

§ In the dilatation region, the flow velocity reduces significantly in magnitude as com- pared to overlapped stenotic region or normal region of the vessel.

§ The effect of non-Newtonian behaviour on the flow profiles is plausible. The profiles are not of parabolic shape.

§ The increase in yield stress is due to unhealthy condition. The plug flow radius increases irrespective of the locations whether it falls within stenotic region or in the region of dilatation. On increasing power law index, the fluid gets impetus even in the plug flow region.

§ Near the extremities of the overlapped stenosis, the steepest enhancement in the skin-friction is observed.

§ With the increase in power law index, profound reduction in skin-friction is obtained. This suggests the importance of power law index of the blood during the theoretical studies.

§ In case of unhealthy situation, the yield stress increases from its normal value. Hence with the increase in yield stress, the resistance to flow ratio increases that results in reduction in the transportation of blood, oxygen and nutrients supply to the tissues.

§ The inclination α of the vessel affects significantly on the impedance to the flow.

§ The impedance to the flow reduces with the increase of power law index.

List of Symbols

g Acceleration due to gravity.

R Radius of unstenosed artery.

h(z) Radius of stenosed artery defined in the Equation (1).

6L Length of stenosed portion shown in the Figure 1.

2L Length of dilatation portion shown in the Figure 1.

n, k Non-Newtonian parameters taken in the Equation (11).

Pressure drop.

Q Volumetric flow rate.

r Radial coordinate.

w Dimensionless axial velocity.

z Axial coordinate.

Greek Letters

Inclination of artery with horizontal.

Maximum thickness of the stenosis.

Maximum thickness of the dilatation.

Resistance of flow in stenosed artery.

Resistance of flow in normal artery.

Impedance.

Viscosity.

Density of fluid.

Shear stress.

Yield stress.

Submit or recommend next manuscript to SCIRP and we will provide best service for you:

Accepting pre-submission inquiries through Email, Facebook, LinkedIn, Twitter, etc.

A wide selection of journals (inclusive of 9 subjects, more than 200 journals)

Providing 24-hour high-quality service

User-friendly online submission system

Fair and swift peer-review system

Efficient typesetting and proofreading procedure

Display of the result of downloads and visits, as well as the number of cited articles

Maximum dissemination of your research work

Submit your manuscript at: http://papersubmission.scirp.org/

Or contact jbise@scirp.org

References

[1] Sochi, T. (2010) Non-Newtonian Flow in Porous Media. Polymer, 51, 5007-5023.

https://doi.org/10.1016/j.polymer.2010.07.047

[2] Shukla, J.B., Parihar, R.S. and Rao, B. (1980) Effects of Stenosis on Non-Newtonian Flow of the Blood in an Artery. Bulletin of Mathematical Biology, 42, 283-294.

https://doi.org/10.1007/BF02460787

[3] Fisher, C. and Rossmann, R. (2009) Effect of Non-Newtonian Behavior on Hemodynamics of Cerebral Aneurysms. Journal of Biomechanical Engineering, 131, Article ID: 091004.

https://doi.org/10.1115/1.3148470

[4] Revellin, R., Rousset, F., Baud, D. and Bonjour, J. (2009) Extension of Murray’s Law Using a Non-Newtonian Model of Blood Flow. Theoretical Biology and Medical Modelling, 6, 7.

https://doi.org/10.1186/1742-4682-6-7

[5] Barnes, H.A. (1997) Thixotropy—A Review. Journal of Non-Newtonian Fluid Mechanics, 70, 1-33.

https://doi.org/10.1016/S0377-0257(97)00004-9

[6] Bodnar, T., Sequeira, A. and Prosi, M. (2011) On the Shear-Thinning and Viscoelastic Effects of Blood Flow under Various Flow Rates. Applied Mathematics and Computation, 217, 5055-5067.

https://doi.org/10.1016/j.amc.2010.07.054

[7] Lou, Z. and Yang, W.J. (1993) A Computer Simulation of the Non-Newtonian Blood Flow at the Aortic Bifurcation. Journal of Biomechanics, 26, 37-49.

https://doi.org/10.1016/0021-9290(93)90611-H

[8] Ku, D.N., Giddens, D.P., Zarins, C.K. and Glagov, S. (1985) Pulsatile flow and Atherosclerosis in the Human Carotid Bifurcation. Positive Correlation between Plaque Location and Low and Oscillating Shear Stress. Arteriosclerosis, 5, 293-302.

https://doi.org/10.1161/01.ATV.5.3.293

[9] Glagov, S., Zarins, C.K., Giddens, D.G. and Ku, D. (1988) Hemodynamics and Atherosclerosis, Insights and Perspectives Gained from Studies of Human Arteries. Archives of Pathology & Laboratory Medicine, 112, 1018-1031.

[10] Pedersen, E., Sung, H., Burlson, A. and Yoganathan, A. (1993) Two-Dimensional Velocity Measurements in a Pulsatile Flow Model of the Normal Abdominal Aortic Simulating Different Hemodynamic Conditions. Journal of Biomechanics, 26, 1237-1247.

https://doi.org/10.1016/0021-9290(93)90071-L

[11] Hademenos, G.J. (1995) The Physics of Cerebral Aneurysms. Physics Today, 48, 24-30.

https://doi.org/10.1063/1.881442

[12] Fox, J.L. and Hugh, A. (1996) Location of Atheroma: A Theory Based on Boundary Layer Separation. British Heart Journal, 28, 388-399.

https://doi.org/10.1136/hrt.28.3.388

[13] Malek, A.M., Alper, S.L. and Izumo, S. (1999) Hemodynamic Shear Stress and Its Role in Atherosclerosis. JAMA, 282, 2035-2042.

https://doi.org/10.1001/jama.282.21.2035

[14] Soliman, E. (2001) Cerebral Aneurysm. eMed. J., 2.

[15] Chakravarty, S. and Mandal, P.K. (1996) A Nonlinear Two-Dimensional Model of Blood Flow in an Overlapping Arterial Stenosis Subjected to Body Acceleration. Mathematical and Computer Modelling, 24, 43-58.

https://doi.org/10.1016/0895-7177(96)00079-9

[16] Srivastava, V.P. (2002) Particulate Suspension Blood Flow through Stenotic Arteries: Effects of Hematocrit and Stenosis Shape. Indian Journal of Pure and Applied Mathematics, 33, 1353-1360.

[17] Prasad, K.M. and Radhakrishnamacharya, G. (2008) Flow of Herschel—Bulkley Fluid through an Inclined Tube of Non-Uniform Cross-Section with Multiple Stenoses. Archives of Mechanics, 60, 161-172.

[18] Srivastava, V.P., Rastogi, R. and Vishnoi, R. (2010) A Two-Layered Suspension Blood Flow through an Overlapping Stenosis. Computers and Mathematics with Applications, 60, 432- 441.

https://doi.org/10.1016/j.camwa.2010.04.038

[19] Chakraborty, U.S., Biswas, D. and Paul, M. (2011) Suspension Model Blood Flow through an Inclined Tube with an Axially Non-Symmetrical Stenosis. Korea-Australia Rheology Journal, 23, 25-32.

https://doi.org/10.1007/s13367-011-0004-8

[20] Sharma, M.K., Bansal, K. and Bansal, S. (2012) Pulsatile Unsteady Flow of Blood through Porous Medium in a Stenotic Artery under the Influence of Transverse Magnetic Field. Korea-Australia Rheology Journal, 24, 181-189.

https://doi.org/10.1007/s13367-012-0022-1

[21] Wu, J., Liu, G., Huang, W., Ghista, D.N. and Wong, K.L. (2015) Transient Blood Flow in Elastic Coronary Arteries with Varying Degrees of Stenosis and Dilatations: CFD Modelling and Parametric Study. Computer Methods in Biomechanics and Biomedical Engineering, 18, 1835-1845.

https://doi.org/10.1080/10255842.2014.976812

[22] Nadeem, S. and Ijaz, S. (2015) Theoretical Analysis of Metallic Nanoparticles on Blood Flow through Tapered Elastic Artery with Overlapping Stenosis. IEEE Transactions on Nanobioscience, 14, 417-428.

https://doi.org/10.1109/TNB.2015.2389253

[23] Oka, S. (1981) Cardiovascular Hemorheology. Cambridge University Press, Cambridge.

[24] Moore, K.L. (1990) Clinically Oriented Anatomy. Williams and Wilkins, Baltimore.

[25] Chaturani, P. and Ponnalagarsamy, R. (1985) A Study of Non-Newtonian Aspects of Blood Flow through Stenosed Arteries and Its Applications in Arterial Diseases. Biorheology, 22, 521-531.

[26] Radtke, H., Schneider, R., Witt, R., Kiesewetter, H. and SchmidSchönbein, H. (1984) A Measuring Device to Determine a Universal Parameter for the Flow Characteristics of Blood: Measurement of the Yield Shear Stress in a Branched Capillary. Advances in Experimental Medicine and Biology, 169, 851-857.

https://doi.org/10.1007/978-1-4684-1188-1_79

[27] Picart, C., Carpentier, P.H., Galliard, H. and Piau, J.M. (1999) Blood Yield Stress in Systemic Sclerosis. American Journal of Physiology—Heart and Circulatory Physiology, 276, H771-H777.

[28] Lee, B., Xue, S., Nam, J., Lim H. and Shin S. (2011) Determination of the Blood Viscosity and Yield Stress with a pressure-Scanning Capillary Hemorheometer Using Constitutive Model. Korea-Australia Rheology Journal, 23, 1-6.

https://doi.org/10.1007/s13367-011-0001-y

[29] Sankar, D.S. and Hemalatha, K. (2007) A Non-Newtonian Fluid Flow Model for Blood Flow through a Catheterized Artery-Steady Flow. Applied Mathematical Modelling, 31, 1847-1864.

https://doi.org/10.1016/j.apm.2006.06.009

[30] Pincombe, B. and Mazumdar, J. (1997) The Effects of Post-Stenotic Dilatations on the Flow of a Blood Analogue through Stenosed Coronary Arteries. Mathematical and Computer Modelling, 25, 57-70.

https://doi.org/10.1016/S0895-7177(97)00039-3