Abstract
As hemodynamics is widely believed to correlate with anastomotic stenosis in coronary bypass surgery, this paper investigates the flow characteristics and distributions of the hemodynamic parameters (HPs) in a coronary bypass model (which includes both proximal and distal anastomoses), under physiological flow conditions. Disturbed flows (flow separation/reattachment, vortical and secondary flows) as well as regions of high oscillatory shear index (OSI) with low wall shear stress (WSS), i.e., high-OSI-and-low-WSS and low-OSI-and-high-WSS were found in the proximal and distal anastomoses, especially at the toe and heel regions of distal anastomosis, which indicate highly suspected sites for the onset of the atherosclerotic lesions. The flow patterns found in the graft and distal anastomoses of our model at deceleration phases are different from those of the isolated distal anastomosis model. In addition, a huge significant difference in segmental averages of HPs was found between the distal and proximal anastomoses. These findings further suggest that intimal hyperplasia would be more prone to form in the distal anastomosis than in the proximal anastomosis, particularly along the suture line at the toe and heel of distal anastomosis.
Similar content being viewed by others
Abbreviations
- B n :
-
Fourier coefficients
- D G :
-
diameter of graft (4 mm)
- f :
-
frequency in Hz
- i :
-
\( {\sqrt { - 1} }, \) unit imaginary number
- J 0 :
-
the first kind of Bessel function of order 0
- J 1 :
-
the first kind of Bessel function of order 1
- N :
-
the total number of Fourier transform terms
- OSI:
-
oscillatory shear index, defined as \( {\text{OSI}} = \frac{1} {2}{\left( {1 - \frac{{{\left| {{\int\limits_0^T {\tau _{{\text{w}}} {\text{d}}t} }} \right|}}} {{{\int\limits_0^T {{\left| {\tau _{{\text{w}}} } \right|}{\text{d}}t} }}}} \right)} \)
- p :
-
static pressure (Pa)
- Q G :
-
mean flow rate in the graft during the pulsatile flow cycle (m3/s)
- Q(t) :
-
flow rate at time t (m3/s)
- r :
-
radial location (m)
- R :
-
radius of the aorta (12.5 mm)
- Re :
-
Reynolds number of aorta, defined as Re = 2ρuR/μ
- t :
-
time (s)
- T :
-
time of a period (s)
- u i :
-
velocity in i direction (m/s), i = 1, 2, 3 for x, y, z directions in Cartesian coordinate respectively
- u(r, t) :
-
the distribution of axial velocity (m/s) at different radial location (r) and time (t)
- WSS:
-
time-averaged wall shear stress (Pa), defined as \( {\text{WSS}} = \frac{1} {T}{\int\limits_0^T {{\left| {\tau _{{\text{w}}} } \right|}} }{\text{d}}t \)
- WSSG:
-
normalized time-averaged spatial wall shear stress gradient, defined as \( {\text{WSSG}} = \frac{1} {T}\frac{{D_{{\text{G}}} }} {{\tau _{{\text{G}}} }}{\int\limits_0^T {{\sqrt {{\left( {\frac{{\partial \tau _{x} }} {{\partial x}}} \right)}^{2} + {\left( {\frac{{\partial \tau _{y} }} {{\partial y}}} \right)}^{2} + {\left( {\frac{{\partial \tau _{z} }} {{\partial z}}} \right)}^{2} } }} }{\text{d}}t \)
- x i :
-
location in Cartesian coordinate (m), i = 1, 2, 3 for x, y, z directions, respectively
- α :
-
Womersley number, defined as \( \alpha = R{\sqrt {\omega /\upsilon } } \)
- μ :
-
dynamic viscosity of the working fluid (Pa s)
- ρ :
-
density of the working fluid (kg/m3)
- τ G :
-
Poiseuille type wall shear stress at the graft corresponding to the mean flow rate in the graft (Pa), defined as \( \tau _{{\text{G}}} = \frac{{32\mu Q_{{\text{G}}} }} {{\pi D^{{\text{3}}}_{{\text{G}}} }} \)
- τ i :
-
wall shear stress in Cartesian coordinate (Pa), i = x, y, z for x, y, z directions, respectively
- τ w :
-
wall shear stress (Pa), defined as \( \tau _{{\text{w}}} = \left. {\mu (\partial u/\partial n)} \right|_{{{\text{wall}}}} , \) where ∂u/∂n| wall is the normal velocity gradient at the wall
- υ :
-
kinematic viscosity of the working fluid (m2/s)
- ω :
-
angular frequency in radian per second of the oscillatory motion, defined as ω = 2πf
References
Asakura T, Karino T (1990) Flow patterns and spatial distribution of atherosclerotic lesions in human coronary arteries. Circ Res 66(4):1045–1066
Bassiouny HS, White S, Glagov S, Choi E, Giddens DP, Zarins CK (1992) Anastomotic intimal hyperplasia: Mechanical injury or flow induced. J Vasc Surg 15:708–717
Buchanan JR, Kleinstreuer C, Hyun S, Truskey GA (2003) Hemodynamics simulation and identification of susceptible sites of atherosclerotic lesion formation in a model abdominal aorta. J Biomech 36:1185–1196
Buchanan JR, Kleinstreuer C, Truskey GA, Lei M (1999) Relation between non-uniform hemodynamics and sites of altered permeability and lesion growth at the rabbit aorto-celiac junction. Atherosclerosis 143:27–40
Chua LP, Zhang JM, Yu SCM, Ghista DN, Tan YS (2005) Numerical study on the pulsatile flow characteristics of proximal anastomotic model. Proc Inst Mech Eng H 219(5):361–379
Chua LP, Zhang JM, Zhou TM (2005) Numerical study of a complete anastomosis model for the coronary artery bypass. Int Commun Heat Mass 32:473–482
Cole JS, Watterson JK, O’Reilly MJG (2002) Numerical investigation of the haemodynamics at a patched arterial bypass anastomosis. Med Eng Phys 24(6):393–401
Eckstein FS, Bonilla LF, Englberger L, Immer FF, Berg TA, Schmidli J, Carrel TP (2002) The St Jude Medical symmetry aortic connector system for proximal vein graft anastomoses in coronary artery bypass grafting. J Thorac Cardiovasc Surg 123(4):777–782
Filippone A (2004) 3-D flow separation. Advanced topics in aerodynamics [online]. Available: http://aerodyn.org/Wings/3dsepar.html
Galjee MA, Rossum ACV, Doesburg T, Eenige MJ, Visser CA (1996) Value of magnetic resonance imaging in assessing patency and function of coronary artery bypass graft. Circulation 93(4):660–666
Ganong WF (1973) Review of medical physiology. Los Altos, Calif, p 453
Hofer M, Rappitsch G, Perktold K, Trubel W, Schima H (1996) Numerical study of wall mechanics and fluid dynamics in end-to-side anastomoses and correlation to intimal hyperplasia. J Biomech 29:1297–1308
Kajiya F, Tsujioka K, Ogasawara Y, Wada Y, Matsuoka S, Kanazawa S, Hiramatsu O, Tadaoka SI, Goto M, Fujiwara T (1987) Analysis of flow characteristics in poststenotic regions of the human coronary artery during bypass graft surgery. Circulation 76(5):1092–1100
Kassab GS, Navia JA (2006) Biomechanical considerations in the design of graft: the homeostasis hypothesis. Annu Rev Biomed Eng 8:499–535
Kleinstreuer C, Hyun S, Buchanan JR, Longest PW, Archie JP, Truskey GA (2001) Hemodynamic parameters and early intimal thickening in branching blood vessels. Crit Rev Biomed Eng 29(1):1–64
Lee D, Su JM, Liang HY (2001) A numerical simulation of steady flow fields in a bypass tube. J Biomech 34:1407–1416
Leuprecht A, Perktold K, Prosi M, Berk T, Trubel W, Schima H (2002) Numerical study of hemodynamics and wall mechanics in distal end-to-side anastomoses of bypass grafts. J Biomech 35:225–236
Longest PW, Kleinstreuer C (2003) Numerical simulation of wall shear stress conditions and platelet localization in realistic end-to-side arterial anastomoses. J Biomech Eng T ASME 125:671–681
Loth F, Jones SA, Zarins CK, Giddens DP, Nassar RF, Glagov S, Bassiouny HS (2002) Relative contribution of wall shear stress and injury in experimental intimal thickening at PTFE end-to-side arterial anastomoses. J Biomech Eng T ASME 124:44–50
Lu B, Dai R, Jing B, Bai H, He S, Zhang N, Wu Q, Zhu X (2001) Evaluation of coronary artery bypass graft patency using three dimensional reconstruction and flow study of electron beam tomography. Chin Med J 114(5):466–472
Moore JA, Steinman DA, Prakash S, Johnston KW, Ethier CR (1999) A numerical study of blood flow patterns in anatomically realistic and simplified end-to-side anastomoses. J Biomech Eng T ASME 121:265–272
Mowat DH, Haites NE, Rawles JM (1983) Aortic blood velocity in healthy adults using simple ultrasound techniques. Cardiovasc Res 17:75–80
Peacock J, Jones T, Tock C, Lutz R (1998) The onset of turbulence in physiological pulsatile flow in a straight tube. Exp Fluids 24:1–9
Roache PJ (1998) Verification of codes and calculations. AIAA J 36:696–702
Sankaranarayanan M., Chua LP, Ghista DN, Tan YS (2005) Computational model of blood flow in the aorto-coronary bypass graft, Biomed Eng On-line J 4(14):13
Shah PJ, Gordon I, Fuller J, Seevanayagam S, Rosalion A, Tatoulis J, Raman JS, Buxton BF (2003) Factors affecting saphenous vein graft patency: clinical and angiographic study in 1402 symptomatic patients operated on between 1977 and 1999. J Thorac Cardiovasc Surg 126(6):1972–1977
Sottiurao VS, Yao JST, Batson RC, Lim SS, Jones R, Nakamura YA (1989) Distal anastomotic intimal hyperplasia: histopathologic character and biogenesis. Ann Vasc Surg 3:26–33
Tang D, Yang C, Kobayashi S, Ku DN (2001) Steady flow and wall compression in stenotic arteries: a three-dimensional thick-wall model with fluid-wall interactions. J Biomech Eng T ASME 123:548–557
Womersley JR (1955) Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known. J Physiol 127:553–563
Acknowledgments
The financial support of A*STAR Project 0221010023 is gratefully acknowledged.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Zhang, JM., Chua, L.P., Ghista, D.N. et al. Numerical investigation and identification of susceptible sites of atherosclerotic lesion formation in a complete coronary artery bypass model. Med Biol Eng Comput 46, 689–699 (2008). https://doi.org/10.1007/s11517-008-0320-4
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11517-008-0320-4