Numerical Solution of Blood Flow and Mass Transport in an Elastic Tube with Multiple Stenoses

Alsemiry, Reima D.; Sarifuddin, undefined; Mandal, Prashanta K.; Sayed, Hamed M.; Amin, Norsarahaida

doi:https://doi.org/10.1155/2020/7609562

BioMed Research International

On this page

Abstract Introduction Results Conclusion Data Availability Disclosure Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 7609562 | https://doi.org/10.1155/2020/7609562

Numerical Solution of Blood Flow and Mass Transport in an Elastic Tube with Multiple Stenoses

Reima D. Alsemiry,^1,2 Sarifuddin,³Prashanta K. Mandal,⁴Hamed M. Sayed,^1,5and Norsarahaida Amin²

Academic Editor: Hwa-Liang Leo

Received16 Sept 2019

Accepted18 Dec 2019

Published31 Jan 2020

Abstract

The simultaneous effect of flexible wall and multiple stenoses on the flow and mass transfer of blood is investigated through numerical computation and simulations. The solution is obtained using the Marker and Cell technique on an axisymmetric model of Newtonian blood flow. The results compare favorably with physical observations where the pulsatile boundary condition and double stenoses result in a higher pressure drop across the stenoses. The streamlines, the iso-concentration lines, the Sherwood number, and the mass concentration variations along the entire wall segment provide a comprehensive analysis of the mass transport characteristics. The double stenoses and pulsatile inlet conditions increase the number of recirculation regions and effect a higher mass transfer rate at the throat, whereby more mass is expected to accumulate and cause further stenosis.

1. Introduction

Caro et al. [1] postulated that atherosclerosis, which is a narrowing of the artery as a result of plaque build-up may occur due to shear-dependent mass transfer mechanism between blood cholesterol and the arterial wall. Cholesterol exists in blood in the form of low density lipoproteins (LDLs) whose deposition along the walls of the artery is a key step in atherogenesis, which would lead to stenosis. Stenosis can affect the velocity of blood flowing through the artery, affecting blood pressure, collapsing the heart, which could in turn lead to disastrous consequences. Thus, an understanding of the behavior of local mass transport in arterial stenosis is important in the study of the formation and development of atherosclerotic lesions for appropriate assessment on the possible correlation between the site of atherosclerotic lesions and the pattern of mass transport.

Ethier [2] carried out computational modelling of mass transfer and studied its links to atherosclerosis. Other studies on mass transport and fluid flow in stenotic arteries of axisymmetric and asymmetric models have been carried out by [3–6]. In these studies, the arterial wall was considered as rigid and the artery is assumed to have single mild stenosis, in which the geometry of the stenosis is represented by the usual cosine curve along with a restriction that the ratio of the severity of stenosis and the radius of the artery is very small. In reality, this is not the case where in many medical situations, the patient is found to have multiple stenoses in the same arterial segment.

Investigations on the effect of multiple stenoses on blood flow have been carried out amongst others by [7–10]. These studies showed that from both experimental results and theoretical calculations, the total effect of a series of noncritical stenoses is approximately equal to the sum of their individual effects where they can be critical and produce symptoms of arterial insufficiency. The flow energy loss due to the presence of the stenoses, which is directly related to the pressure drop across them, increases with the number of stenoses but is not strongly dependent on the spacing between them. The authors of [11–19] have also investigated blood flow through multiple stenoses; however, these studies have not considered the mass transfer.

Another aspect to be considered in arterial blood flow is the cyclic nature of the heart pump which creates pulsatile conditions in the arteries, giving rise to unsteady flow. It is observed that most CFD models of arterial hemodynamics make the simplifying assumptions of rigid walls and fully developed inlet velocities (cf. [13–19]). But the arteries are not rigid tubes. They adapt to varying flow conditions by enlarging or shrinking. All of these physiological conditions make the modelling and consequently the solution to be almost impossible to be obtained analytically and challenging computationally. Nandakumar and Anand [20] studied steady and pulsatile flow of blood through a channel with single as well as double stenoses on the assumption that the pulsations of flow are damped in the small vessels; thus the flow is effectively steady in the capillaries and the veins while Liu and Tang [21] investigated the influence of distal stenosis on blood flow through curved arteries with two stenoses. But again, these studies on pulsatile flow have also not considered the mass transfer. In another study, Layek et al. [22] investigated the effect of multiple stenoses on the flow of Newtonian fluid in a rigid tube and opined that the disturbance created by the constrictions is mainly concentrated at the downstream of the last constriction. Considering the flow of Newtonian fluid in a two-dimensional channel having a single constriction, Layek and Midya [23] concluded that the maximum stress and the length of the recirculation region associated with two shear layers of the constriction do increase with the increasing area reduction of the constriction. They further concluded that the flow-field separates after the symmetry breaking bifurcation, and the symmetry of the flow depends on Reynolds’ number and the height of the constriction. The flow of a fluid having hematocrit-dependent viscosity past a tube with partially overlapped constriction has been investigated by Layek et al. [24]. They observed that the peak value of the wall shear stress decreases with increasing haematocrit parameter while a reverse trend is observed for the flow separation region. They also opined that the deformability of the wall does reduce the wall shear stress as compared to the rigid wall case. All these studies [22–24] ignored the flow pulsatility and/or consideration of multiple constrictions, and the mass transfer as well which plays a pivotal role in the genesis and evolution of atherosclerosis.

Based on the gap established above, with regard to studies involving mass transfer, the following work seeks to analyze the flow and mass transfer characteristics of pulsatile blood flow through an artery with double stenoses. The fluid considered is Newtonian in an axisymmetric setting, while the pair of stenoses vary in severities, lengths, and distances between them. The equation for stenoses is given in an algebraic form which could represent both moderate and severe stenoses instead of the usual cosine function which could only describe mild stenosis. The objective of the present study lies in the consideration of the transport of mass as well as momentum together through a tube with a flexible wall, resembling the flexibility of the artery in the presence of double stenoses. The flow pulsatility cannot be ruled out from the present investigation.

2. Formulation of the Problem

We consider a fully developed two-dimensional axisymmetric flow of an incompressible Newtonian fluid of density in a tube. The relevant equations of motions in vector forms are the continuity, momentum, and mass as follows:with is the material derivative, where and are the radial and axial velocity components, respectively, is the pressure, is the constant viscosity, is the mass concentration, and is the coefficient of mass diffusion.

In the cylindrical coordinate system, the corresponding equations (1)–(3) are written in a conservative form as follows:

The schematic diagram for the double stenoses is given in Figure 1, where is the radius of the artery in the stenotic region and is the radius of the artery in the nonstenotic regions. are the critical heights of the first and second stenosis respectively; is the inlet segment, is the distance between stenoses, , are the lengths of stenoses, and is the length of the arterial segment under consideration.

The equations describing the stenoses are given by the following:

The time-variant parameter is given by with representing the amplitude parameter and the angular frequency is given by , being the pulse frequency and . To the best of our knowledge equation (8) is the first equation to address double stenoses without any control on the severity of stenoses which has not been considered before.

2.1. Boundary Conditions

where is the cross-sectional average velocity of the fluid and is a constant.

3. Solution Procedure

The solution procedure involves the nondimensionalization, radial coordinate transformation, and the finite-difference Marker and Cell method (MAC) initially proposed by Harlow and Welch [25]. Sarifuddin et al. [26, 27] and Mustapha et al. [15, 16] have used the method to solve blood flow problems.

3.1. Nondimensionalization of the Equations

The nondimensional variables and parameters introduced are as follows:

Using (13), equations (4)–(7) have their respective nondimensional forms as follows (omitting bar):

The boundary conditions (9)–(12) reduce to their respective dimensionless forms:where Re is the Reynolds number, is the Schmidt number and is the Womersley number defined as follows:

3.2. Radial Coordinate Transformation

With the introduction of a radial coordinate transformation , equations (14)–(17) now become as follows:and the boundary conditions (18)–(21) become

3.3. Finite-Difference Method

The solution procedure consists of discretization of the governing equations, combining the discretized forms of the momentum and continuity equations to obtain the Poisson equation for pressure, the successive overrelaxation (SOR) method, and the pressure and velocity corrections. The schematic computational domain is given in Figure 2. The velocities and pressure are calculated at different locations of the control volume, as indicated in Figure 3. The difference equations are derived at three distinct cells, each corresponding to the continuity, axial and radial momentum equations.

The discretization of the time derivative terms is based on the first-order accurate two-level forward time differencing formula, while the convective terms in the momentum equations are discretized with a hybrid formula consisting of central differencing and second-order upwinding scheme (cf. Courant et al. [28]). The diffusive terms are discretized using second-order accurate three-point central difference formula. Thus in a finite-difference formula with and where refers to time and is the time increment. The length and width of the cell of the control volume are represented by and , respectively.

The discretized version of the continuity equation (24) at the (i, j) cell iswhere

Here and represent the respective coordinates of the center of the cell and the cell faces as shown in Figure 3, while and stand for -velocities at the top and bottom middle positions of the control volume of the continuity equation. The momentum equations (24) and (25) are written in the following forms:where are the finite-difference representation of convective and diffusive terms of the axial and radial momentum at the nth time level.

The Poisson equation for pressure is derived from equations (31)–(33) which takes the final forms:where represents the discretized form of the divergence of the velocity field at the (i, j) cell and the expressions for are the same as Mustapha et al. [15, 16]. The Poisson equation (34) for pressure is then solved using the successive overrelaxation (SOR) method to obtain the intermediate pressure field.

The increment is chosen to be along , while along is chosen to be equal to or less than a prescribed stability criterion as depicted in Figure 4, where here c is taken to be 0.05. (cf. [25, 29]).

The number of iterations is limited up to 10. The pressure and velocities then go through a correction stage to achieve better accuracy. The process is described in [15, 16, 26, 27]. When the velocity field has been obtained, the mass concentration is calculated from the respective discretized versions of equation (26) with the relevant boundary conditions (equations (27)–(30)). The values chosen for are 0.05, 2, and 3, respectively.

4. Results and Discussions

The influence of the pulsatile inlet is reflected in Figure 5 where the pressure drop in this case is higher than the one generated with the parabolic inlet. It decreases with increasing Re with a strong linear correlation between them in both cases of parabolic and pulsatile conditions (cf. Sarifuddin et al. [26]). Further, the pressure drop is seen to increase with the number of stenoses, which agrees with the experimental study of Talukder et al. [8].

The behavior of the axial and the radial velocity at the narrowest points ( and ) for different Re are shown in Figures 6(a), 6(b), 7(a), and 7(b). The axial velocity has positive values, and it is noted that the parabolic case results in higher velocity. It is also observed that the axial velocity near the wall increases with increasing Re; however, there is a cross over at and for and , respectively. Figures 7(a) and 7(b) show that the radial velocity corresponding to the pulsatile inlet assume positive values at except the value on the wall at () while negative values are observed at the narrowest point () and the flow velocity increases with increasing Re near the centerline while it is reduced near the wall with increasing Re. Both figures reveal that the velocities in the case of the parabolic inlet are negative and substantially less than that with the pulsatile inlet.

(a)

(b)

(a)

(b)

Figure 8(a) exhibits the axial velocity profiles at different locations of the stenosed arterial segment at Re = 300 for both parabolic and pulsatile inlet conditions. At (), the velocity in the parabolic case is higher than in the pulsatile case. A backflow occurs in the pulsatile case at the downstream of the narrowest point () near the wall. The curves decrease from their individual maximum at the axis as one moves away from it and finally they approach a minimum value (zero) on the wall surface. Note that the curves of the axial velocity at and at are coincident. The axial velocity at the critical height of the second stenosis is considerably higher than that of the first stenosis . Figure 8(b) shows that the radial velocity has positive values everywhere at , and , excluding the position on the wall. At and the narrowest point (), the radial velocity is observed to have all negative values. The nonzero values of the radial velocity near the wall clearly reflect the influence of the radial motion of the arterial wall in the pulsatile case.

(a)

(b)

Figure 9(a) exhibits the distribution of the wall shear stress (wss) along the arterial segment for different Re considering pulsatile as well as parabolic inlet conditions. The results show that wss for both the parabolic and pulsatile inlet conditions attain their peaks at the critical heights of the stenoses (). It is observed that separation occurs (negative values of wss) only at the downstream of the second stenosis for parabolic inlet condition. In the pulsatile case, the separation zone occurs between the two stenoses with multiple separation regions at the downstream of the second stenosis. Then, wss starts to increase slowly towards the wall surface (reattachment point).

(a)

(b)

(c)

The effect of different severities on wss is depicted in Figure 9(b). In the pulsatile case, when the two stenoses have the same severities and, , flow separation occurs at two specific places at the downstream of the first and the second stenoses with different peak values. In the case of stenoses with different severities and a larger separated region is formed at the downstream of the more severe stenoses (cf. Johnston and Kilpatrick [12]). A smaller separation region is produced in the case of pulsatile inlet condition and the peak wss is much higher than the parabolic inlet condition.

Figure 9(c) determines the effects of the length of stenoses on wss. Peak wall shear stress decreases with increasing the length of stenosis but it increases with the gap between stenosis and at this position, there is a potential that plaque would rupture whereas, at the low shear stress position, atherosclerotic development may be induced. These phenomena of separation and reattachment are due to the adverse pressure in these regions and are believed to be responsible for the malfunctioning of the cardiovascular system having atherosclerotic plaque.

Figures 10(a)–10(d) show the instantaneous patterns of streamlines governing the flow of blood through the stenoses in case of for both parabolic and pulsatile inlet conditions at Re = 300 and Re = 500. In the parabolic case, only one recirculation zone developed between the two stenoses where separation occurs at (c.f Figure 9(b)). In the case of the pulsatile inlet, a multiple recirculation region is noted in between the two stenoses (separation point ). Thus, an increase in Re and a consideration of the pulsatile flow increase the number of the recirculation region.

(a)

(b)

(c)

(d)

Figure 11(a) exhibits the profiles of mass concentration at different positions for both pulsatile and parabolic inlet conditions at Re = 300. The mass concentration at each axial position converges to zero according to the wall condition. The mass concentration gets distorted at the downstream of the stenoses which could be due to the separation of flow at the downstream of the stenoses. Along the upstream of stenoses, the flow velocity and wall shear stress (Figure 8(a)) increase as the flow gets accelerated towards the throat leading to the increase of solute concentration. It is also observed that the concentration at the throat is much higher for pulsatile flow than the parabolic one. Figure 11(b) exhibits the evolution of mass concentration at the throat of stenosis corresponding to single and double stenosis for pulsatile inlet condition. The mass concentration is higher at the throat in the case of double stenosis compared to the single one.

(a)

(b)

The mass flux to the arterial wall which is quantified through the Sherwood number over the entire stenosed arterial segment is examined and exhibited in Figures 12(a)–12(c). The Sherwood number defined by where is the local mass flux to the arterial wall and is the inlet diameter of the artery. It is observed that increases with increasing Re while distribution appreciably changes specifically at the throat, between the two stenoses and at the downstream position. The highest mass transfer is experienced at the upstream, while the minimum value occurs at the downstream of the stenoses. Note that pulsatile flow increases the Sherwood number and that it is much higher in the case of double stenosis.

(a)

(b)

(c)

The iso-concentration lines considering pulsatile as well as parabolic inlet conditions based on are displayed in Figures 13(a) and 13(b). The iso-concentration lines for parabolic and pulsatile inlets have different distributions with multiple recirculation regions nearby the downstream of the more severe stenosis in the pulsatile case. The general trend of the iso-concentration lines is that they move away from the inlet region towards the upstream of the stenosis and correspondingly impair the mass transport in this region, while they adhere to the outline of the stenosis at both the upstream and downstream ends. At this region of low wss (compare with Figure 9), cholesterol may tend to accumulate and cause more severe stenosis. This observation conforms with Schneiderman et al. [30].

(a)

(b)

5. Conclusion

The hemodynamics of the pulsatile flow and the transport of mass in an arterial segment having a couple of stenoses have been studied in relation to the distensibility of the vessel wall. Predicted results show that the pulsatile inlet and double stenoses with varying severity affect the flow characteristics significantly, especially the development of the recirculation zone and the peak value of the wall shear stress. It is also predicted that the concentration at the throat (z = 19) is much higher for pulsatile flow than the parabolic inlet condition. Moreover, the pair of stenoses contributes much to the mass concentration than the case of single stenosis. The mass flux to the arterial wall (Sherwood number) does increase with the increasing values of Re and here too, mass flux increases with the flow pulsatility and the presence of double stenoses. At the downstream, cholesterol may tend to accumulate and causes more severe stenosis. For severe stenoses, the peak value of the wall shear stress is higher in the pulsatile flow case and the iso-concentration lines show more recirculation regions nearby the downstream end and their lengths are longer. In conclusion, the results presented agree well with physical observations and provide an insight into the link between atherosclerosis, stenosis, and the pattern of mass transport.

Though the detailed knowledge of the dynamical variables is possible and provides useful elements, the mechanism of influence of the haemodynamical factors in the arterial disease is not clear. The characteristics of the red cells must be taken into consideration by including a shear-dependent viscosity in the diffusion terms in time-dependent flows highlighting the scope of further work. All these mechanical and biochemical aspects related to the biofluid dynamics are of some importance and demand further investigation. A great deal of work is needed to establish the rheological parameters for the physiological values and to understand the connection of the issues with biological facts.

Data Availability

The data on blood flow parameters used for analysis and validation purposes are from previously reported studies and datasets, which have been cited.

Disclosure

To the best of the authors' knowledge, the content of this work is correct . Each author has contributed in part to the problem formulated, numerical computations and simulations, analysis of results, and manuscript writing.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The first author wishes to acknowledge the financial support provided by the University of Taibah, Al Madinah, and the Saudi Ministry of Education. Part of this work was carried out at the University of Visva-Baharati, India, and Universiti Teknologi Malaysia. The hospitality and assistance rendered by both institutions are gratefully acknowledged.

References

C. G. Caro, J. M. Fitz-Gerald, and R. C. Schroter, “Atheroma and arterial wall shear- observation, correlation and proposal of a shear dependent mass transfer mechanism for atherogenesis,” Proceedings of the Royal Society of London. Series B. Biological Sciences, vol. 177, no. 1046, pp. 109–133, 1971.
View at: Publisher Site | Google Scholar
C. R. Ethier, “Computational modeling of mass transfer and links to atherosclerosis,” Annals of Biomedical Engineering, vol. 30, no. 4, pp. 461–471, 2002.
View at: Publisher Site | Google Scholar
L. H. Back, Y. I. Cho, D. W. Crawford, and R. F. Cuffel, “Effect of mild atherosclerosis on flow resistance in a coronary artery casting of man,” Journal of Biomechanical Engineering, vol. 106, no. 1, pp. 48–53, 1984.
View at: Publisher Site | Google Scholar
M. R. Kaazempur-Mofrad, S. Wada, J. G. Myers, and C. R. Ethier, “Mass transport and fluid flow in stenotic arteries: axisymmetric and asymmetric models,” International Journal of Heat and Mass Transfer, vol. 48, no. 21-22, pp. 4510–4517, 2005.
View at: Publisher Site | Google Scholar
A. Valencia and M. Villanueva, “Unsteady flow and mass transfer in models of stenotic arteries considering fluid-structure interaction,” International Communications in Heat and Mass Transfer, vol. 33, no. 8, pp. 966–975, 2006.
View at: Publisher Site | Google Scholar
M. A. Ikbal, S. Chakravarty, Sarifuddin, and P. K. Mandal, “Numerical simulation of mass transfer to micropolar fluid flow past a stenosed artery,” International Journal for Numerical Methods in Fluids, vol. 67, no. 11, pp. 1655–1676, 2010.
View at: Publisher Site | Google Scholar
P. E. Karayannacos, N. Talukder, R. M. Nerem, S. Roshon, and J. S. Vasko, “The role of multiple noncritical arterial stenoses in the pathogenesis of ischemia,” The Journal of Thoracic and Cardiovascular Surgery, vol. 73, no. 3, pp. 458–469, 1977.
View at: Publisher Site | Google Scholar
N. Taludker, P. E. Karayannacos, R. M. Nerem, and J. S. Vasko, “An experimental study of the fluid dynamics of multiple noncritical stenoses,” Journal of Biomechanical Engineering, vol. 99, no. 2, pp. 74–82, 1977.
View at: Publisher Site | Google Scholar
R. L. Feldman, W. W. Nichols, C. J. Pepine, and C. R. Conti, “Hemodynamic effects of long and multiple coronary arterial narrowings,” Chest, vol. 74, no. 3, pp. 280–285, 1978.
View at: Publisher Site | Google Scholar
R. L. Feldman, W. W. Nichols, C. J. Pepine, D. A. Conetta, and C. Richard Conti, “The coronary hemodynamics of left main and branch coronary stenoses,” The Journal of Thoracic and Cardiovascular Surgery, vol. 77, no. 3, pp. 377–388, 1979.
View at: Publisher Site | Google Scholar
D. Kilpatrick, S. D. Webber, and J.-P. Colle, “The vascular resistance of arterial stenoses in series,” Angiology, vol. 41, no. 4, pp. 278–285, 1990.
View at: Publisher Site | Google Scholar
P. R. Johnston and D. Kilpatrick, “Mathematical modelling of paired arterial stenoses,” in Proceedings Computers in Cardiology, pp. 229–232, Chicago, IL, USA, 1990.
View at: Google Scholar
T. S. Lee, “Numerical studies of fluid flow through tubes with double constrictions,” International Journal for Numerical Methods in Fluids, vol. 11, no. 8, pp. 1113–1126, 1990.
View at: Publisher Site | Google Scholar
N. Mustapha, S. Chakravarty, P. K. Mandal, and N. Amin, “Unsteady response of blood flow through a couple of irregular arterial constrictions to body acceleration,” Journal of Mechanics in Medicine and Biology, vol. 8, no. 3, pp. 395–420, 2008.
View at: Publisher Site | Google Scholar
N. Mustapha, P. K. Mandal, I. Abdullah, N. Amin, and T. Hayat, “Numerical simulation of generalized Newtonian blood flow past a couple of irregular arterial,” Numerical Methods for Partial Differential Equations, vol. 27, no. 4, pp. 960–981, 2011.
View at: Publisher Site | Google Scholar
N. Mustapha, P. K. Mandal, P. R. Johnston, and N. Amin, “A numerical simulation of unsteady blood flow through multi-irregular arterial stenoses,” Applied Mathematical Modelling, vol. 34, no. 6, pp. 1559–1573, 2010.
View at: Publisher Site | Google Scholar
J. C. Misra, A. Sinha, and G. C. Shit, “Mathematical modeling of blood flow in a porous vessel having double stenoses in the presence of an external magnetic field,” International Journal of Biomathematics, vol. 4, no. 2, pp. 207–225, 2011.
View at: Publisher Site | Google Scholar
V. K. Verma and P. Saraswat, “Effect of multiple stenosis on blood flow through a tube,” International Journal of Biology and Biomedical Engineering, vol. 7, pp. 753–756, 2013.
View at: Google Scholar
M. G. Rabby, S. P. Shupti, and M. M. Molla, “Pulsatile non-Newtonian laminar blood flows through arterial double stenoses,” Journal of Fluids, vol. 2014, Article ID 757902, 13 pages, 2014.
View at: Publisher Site | Google Scholar
N. Nandakumar and M. Anand, “Pulsatile flow of blood through a 2D double-stenosed channel: effect of stenosis and pulsatility on wall shear stress,” International Journal of Advances in Engineering Sciences and Applied Mathematics, vol. 8, no. 1, pp. 61–69, 2016.
View at: Publisher Site | Google Scholar
B. Liu and D. Tang, “Influence of distal stenosis on blood flow through coronary serial stenoses: a numerical study,” International Journal of Computational Methods, vol. 16, no. 3, Article ID 1842003, 2019.
View at: Publisher Site | Google Scholar
G. C. Layek, S. Mukhopadhyay, and S. A. Samad, “Oscillatory flow in a tube with multiple constrictions,” International Journal of Fluid Mechanics Research, vol. 32, no. 4, pp. 402–419, 2005.
View at: Publisher Site | Google Scholar
G. C. Layek and C. Midya, “Effect of constriction height on flow separation in a two-dimensional channel,” Communications in Nonlinear Science and Numerical Simulation, vol. 12, no. 5, pp. 745–759, 2007.
View at: Publisher Site | Google Scholar
G. C. Layek, S. Mukhopadhyay, and R. S. R. Gorla, “Unsteady viscous flow with variable viscosity in a vascular tube with an overlapping constriction,” International Journal of Engineering Science, vol. 47, no. 5-6, pp. 649–659, 2009.
View at: Publisher Site | Google Scholar
F. H. Harlow and J. E. Welch, “Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface,” Physics of Fluids, vol. 8, no. 12, pp. 2182–2189, 1965.
View at: Publisher Site | Google Scholar
C. S. Sarifuddin, S. Chakravarty, P. K. Mandal, and G. C. Layek, “Numerical simulation of unsteady generalized Newtonian blood flow through differently shaped distensible arterial stenoses,” Journal of Medical Engineering & Technology, vol. 32, no. 5, pp. 385–399, 2008.
View at: Publisher Site | Google Scholar
C. S. Layek, S. Chakravarty, and P. K. Mandal, “Effect of heat and mass transfer on non-Newtonian flow - links to atherosclerosis,” International Journal of Heat and Mass Transfer, vol. 52, no. 25-26, pp. 5719–5730, 2009.
View at: Publisher Site | Google Scholar
R. Courant, E. Isaacson, and M. Rees, “On the solution of nonlinear hyperbolic differential equations by finite differences,” Communications on Pure and Applied Mathematics, vol. 5, no. 3, pp. 243–255, 1952.
View at: Publisher Site | Google Scholar
A. Amsden and F. Harlow, The SMAC Method: A Numerical Technique for Calculating Incompressible Fluid Flows, University of California, Oakland, CA, USA, 1970.
G. Schneiderman, C. G. Ellis, and T. K. Goldstick, “Mass transport to walls of stenosed arteries: variation with Reynolds number and blood flow separation,” Journal of Biomechanics, vol. 12, no. 11, pp. 869–877, 1979.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Reima D. Alsemiry et al. 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.

PDF Download Citation

Download other formats

Order printed copies

Views

1662

Downloads

1445

Citations