Influence of model boundary conditions on blood flow patterns in a patient specific stenotic right coronary artery
© Liu et al.; licensee BioMed Central Ltd. 2015
Published: 9 January 2015
In literature, the effect of the inflow boundary condition was investigated by examining the impact of the waveform and the shape of the spatial profile of the inlet velocity on the cardiac hemodynamics. However, not much work has been reported on comparing the effect of the different combinations of the inlet/outlet boundary conditions on the quantification of the pressure field and flow distribution patterns in stenotic right coronary arteries.
Non-Newtonian models were used to simulate blood flow in a patient-specific stenotic right coronary artery and investigate the influence of different boundary conditions on the phasic variation and the spatial distribution patterns of blood flow. The 3D geometry of a diseased artery segment was reconstructed from a series of IVUS slices. Five different combinations of the inlet and the outlet boundary conditions were tested and compared.
The temporal distribution patterns and the magnitudes of the velocity, the wall shear stress (WSS), the pressure, the pressure drop (PD), and the spatial gradient of wall pressure (WPG) were different when boundary conditions were imposed using different pressure/velocity combinations at inlet/outlet. The maximum velocity magnitude in a cardiac cycle at the center of the inlet from models with imposed inlet pressure conditions was about 29% lower than that from models using fully developed inlet velocity data. Due to the fact that models with imposed pressure conditions led to blunt velocity profile, the maximum wall shear stress at inlet in a cardiac cycle from models with imposed inlet pressure conditions was about 29% higher than that from models with imposed inlet velocity boundary conditions. When the inlet boundary was imposed by a velocity waveform, the models with different outlet boundary conditions resulted in different temporal distribution patterns and magnitudes of the phasic variation of pressure. On the other hand, the type of different boundary conditions imposed at the inlet and the outlet did not have significant effect on the spatial distribution patterns of the PD, the WPG and the WSS on the lumen surface, regarding the locations of the maximum and the minimum of each quantity.
The observations from this study indicated that the ways how pressure and velocity boundary conditions are imposed in computational models have considerable impact on flow velocity and shear stress predictions. Accuracy of in vivo measurements of blood pressure and velocity is of great importance for reliable model predictions.
The coronary arteries are a common and important site of atherosclerotic lesion development . It is believed that mechanical stresses resulting from intravascular pressure and flow may contribute to the pathogenesis of atherosclerosis. The curvature and the presence of stenosis in coronary arteries result in considerable change to the flow field. Flow disturbances associated with a stenosis could contribute to modifications in the rate of plaque development, the direction of plaque extension and the composition of the plaque [2–4]. A detailed hemodynamic evaluation of disturbed flow and the spatial and temporal flow distribution patterns may give additional insight to understanding the progression of atherosclerosis and may have useful clinical value. However, it is difficult to measure local flow patterns and mechanical forces in vivo with sufficient accuracy. Computational models become a useful tool in that regard.
In recent decades, much progress had been made in research coupling medical imaging and computational fluid dynamics (CFD) to study cardiovascular hemodynamics [1, 4–9]. The methods developed could provide investigators with powerful new tools, rivaling and even surpassing experimental fluid mechanics methods to investigate the mechanisms of disease and to design medical devices and therapeutic interventions . Torii et al. investigated the effects of wall compliance on coronary hemodynamics by coupling FSI analysis of the human RCA in conjunction with physiological velocity and pressure waveforms . Johnston et al. compared models of Newtonian and non-Newtonian flow in four healthy right coronary arteries reconstructed from biplane angiograms . Zeng et al. examined the effects of cardiac motion on right coronary artery hemodynamics . Myers et al. investigated the effects of arterial geometry, pulsatile flow condition and dynamic vessel motion on flow patterns in a model with a patient specific right coronary artery geometry . Tang et al. provided a review of the key steps and findings of their group's in image-based models for human carotid and coronary plagues .
In general the solution of a system of partial differential equations is only uniquely determined if the boundaries of the domain are imposed with necessary conditions. Therefore, to perform the 3D simulations of blood flow in human arteries, proper boundary conditions need to be specified at the inlet and the outlet boundaries when solving the governing Navier-Stokes equations. In addition to the consideration of the well-posedness of the computational problem, more natural and realistic conditions relating pressure and velocity waveforms for image-based modeling may enable researchers to better quantify and predict the local pressure fields and flow distribution in the upstream and downstream domains of the patient specific stenotic artery. On the other hand, improperly imposed boundary conditions may result in misleading hemodynamic information. Usually the velocity profile or the pressure profile is explicitly prescribed at the inlet. There is a variety of approaches for the outlet boundary condition. Three of the commonly used simple outlet conditions are prescribed velocity or pressure profile, or zero normal traction. Different combinations of the above mentioned boundary conditions at the inlet and the outlet have been applied by many researches [4, 5, 9–15].
Several authors have studied the influence of the flow waveform and the spatial profile of the inlet velocity on CFD predictions of blood flow in human arteries [1, 16–18]. Mark et al. demonstrated the importance of unsteadiness in a left anterior descending coronary artery flow phantom . Ethier et al. examined the flow waveform effects on end-to-side anastomotic flow patterns . Myers et al. investigated the effects of the inlet velocity profile and flow waveform on the flow patterns on the hemodynamics in the proximal, medial, and distal arterial regions of human right coronary artery. In their study, a zero-traction condition was imposed at the model outlet and a velocity profile was imposed at the inlet. They compared three different inlet velocity profiles: a fully developed inlet velocity profile, a "blunt" inlet velocity profile, and a fully developed Dean-type velocity profile. They also used two flow waveforms: sinusoidal and physiological waveforms in the study. They concluded that changes in the inlet velocity profiles did not produce significant changes in the arterial velocity and wall shear stress patterns and the waveform shape is not a significant factor in the right coronary time-averaged WSS features . Morbiducci et al. studied the influence of assumptions regarding velocity profiles at the inlet section of the ascending aorta. They concluded that the plausibility of the assumption of idealized velocity profiles as inlet boundary conditions can lead to misleading representations of the aortic hemodynamics both in terms of disturbed shear and bulk flow structures . However, not much work has been reported on comparing the effect of the different combinations of the inlet/outlet boundary conditions on the quantification of the pressure field and flow distribution patterns. The aim of this work is to investigate the influence of the different types of boundary conditions imposed at the inlet/outlet on the spatial and the temporal distribution patterns of the blood flow in a patient-specific stenotic right coronary artery.
Patient-specific plaque geometry and flow data
The flow model and numerical method
where η0 = 0.056 Pa∙s is the zero shear rate viscosity, η∞ = 0.00345 Pa∙s is the infinite shear rate viscosity, is the shear rate, λ = 3.313s is a parameter, and n = 0.3568 is a dimensionless parameter. In the computations, the blood density ρ was assumed to be constant at 1050 kg/m 3 .
V-NS0 model: A fully developed velocity profile with pulse waveform (Figure 2(a)) was imposed at the inlet and a zero normal stress condition was imposed at the outlet.
P-V model: A time dependent pressure with waveform Pin (Figure 2(b)) was imposed at the inlet. A fully developed velocity profile with pulse waveform (Figure 2(a)) was multiplied by a scalar based on the sizes of the inlet/outlet cross sections such that the flow rate at the outlet is equal to that at the inlet. Then this adjusted velocity waveform was imposed at the outlet.
V-P model: A fully developed velocity profile with pulse waveform (Figure 2(a)) was imposed at the inlet and a time dependent pressure with waveform Pout (as shown in Figure 2(b)) was imposed at the outlet.
P-P model: A time dependent pressure with waveform Pin and Pout (as shown in Figure 2(b)) was imposed at the inlet and the outlet boundary, respectively.
P-NS0 model: A time dependent pressure with waveform Pin (Figure 2(b)) was imposed at the inlet and a zero normal stress condition was imposed at the outlet.
In the process of the computer simulations using the models with different combinations of the inlet/outlet boundary conditions, the computation of the P-V model was first carried out. This model applied the real data measured from the patients at the boundaries and thus its simulation results were used for the comparison with other models. The blood pressure at the outlet boundary obtained from P-V model was averaged on the lumen cross section at each time during the cardiac cycle, resulting in a pressure waveform named Pout (see Pout in Figure 2(b)). This pulse pressure waveform was then used as the boundary condition at the outlet for the V-P model and the P-P model. In all models imposed with pressure boundary condition at either the inlet or the outlet the pressure is assumed as uniform spatially over the inlet or the outlet cross section.
Validation of numerical solutions
Global maximums of the mean PD and the mean WPG, and the global maximum and minimum of the mean WSS for five models.
Max PDmean (mmHg)
Max WPGmean (N/m 3 )
Max WSSmean (N/m 2 )
Min WSSmean (N/m 2 )
BCs have a notable effect on the velocity magnitude in upstream region
Figure 5(1a) indicates that the velocity magnitude of the P-V/P-P models is much lower than that of the V-P/V-NS0 models at the center of the inlet lumen cross section, while the velocity magnitude of the P-NS0 model is unrealistically high. The maximum velocity magnitude (0.708 m/s) in a cardiac cycle from P-V/P-P models was about 29% lower than that from V-P/V-NS0 models which used experimental velocity data. Both P-V and P-P models underestimate the blood velocity magnitude along the center line of the artery. With the existence of the secondary flow and the flow shifting in a stenotic curved artery, the blood velocity magnitude is highly non-uniformly distributed. The underestimate of the velocity magnitude of the P-V/P-P models on the center line varies significantly in different regions of the artery. Compared with the V-P/V-NS0 model, the velocity magnitude at the center of the 13th lumen slice resulted from the P-V model or the P-P model is about 8% or 10% lower, respectively. These differences gradually diminish downstream along the artery.
The underestimation of the P-V/P-P models on the velocity magnitude at the center line in the upstream may be explained by the contour plots in Figure 6. All four models have the same flow rate during the cardiac cycle but with two different spatial distribution patterns of the velocity magnitude, one is blunter and the other is fully developed. The blunter inlet flow distribution of the P-V/P-P models results in a lower velocity magnitude at the center.
Effect of the BCs on the WSS
Figure 5(2a&2b) show that the WSS at the inner wall resulted from different model is quite different, especially near the inlet of the artery. Compared with the V-P/V-NS0 models, the maximum wall shear stress at inlet in a cardiac cycle resulted from P-V model or P-P model is higher by 29% or 25%, respectively. The plots of the time average mean WSS along the inner curve and the outer curve in Figure 8(3a&3b) can further confirm the above differences observed. The V-P and the V-NS0 models underestimate the mean WSS by 8%, on average over the lumen wall of the upper stream segment, compared with the P-V model. There is no significant difference between the mean WSS of the P-P and the P-V models in this region (1.5% on average). The cause of the errors in this region maybe the assumption of the inlet fully developed spatial velocity profile on the V-P and the V-NS0 models. The blunter inlet flow distribution of the P-V/P-P models results in a higher WSS on the lumen wall than that of the V-P/V-NS0 models. As shown in Figure 5(2b) and Figure 8(3a&3b), the differences of the WSS between the V-P/V-NS0 and the P-V/P-P models decrease as flow moves from the inlet into the first stenosis. From the neck of the first stenosis to downstream, the WSS resulted from the V-P/V-NS0 models becomes slightly higher than that of the P-V/P-P models, with the maximum and the minimum mean WSS for all models listed in Table 1.
BCs have notable influence on the pressure and the PD
Figure 7(2&3) and Figure 8(1)) indicate that the boundary conditions also have a non-negligible effect on the PD and the WPG, especially near the inlet. The V-P/V-NS0 models underestimate the pressure drop in upper stream segment of the artery. At the third lumen slice near the inlet of the artery, the pressure drop resulted from the V-P/V-NS0 models is less by about 36% at the peak flow (t = 0.425), compared with that resulted from the P-V model, which applied the in vivo data for both inlet and the outlet boundary conditions. This difference is gradually decreased till the neck of the first stenosis. The plots of the mean PD along the axial length of the artery in Figure 8(1a&1b) can further illustrate the above observations. Figure 7(3a&3b) show the similar pattern of the WPG on the difference between the V-P/V-NS0 models and the P-P/P-V models. Compared with the P-V model, the V-P/V-NS0 models underestimate the mean PD and the mean WPG by 43% and10%, respectively, on average over the lumen wall of the segment proximal to the first stenosis, while near the outlet of the artery, the V-P/V-NS0 models slightly overestimate the mean PD and the mean WPG.
Effect of BCs on temporal distribution patterns of blood flows
The ways the boundary conditions are imposed impact the temporal distribution patterns of the blood flows. It is a popular approach in blood flow simulation that the outlet boundary condition is imposed with traction free/zero normal stress condition. When this outlet boundary condition is combined with a velocity or a pressure inlet boundary condition, the simulation may result in certain misleading temporal distribution patterns of blood flow, such as the pressure or the velocity.
As shown in Figure 7(1a&1b), the shape of the pressure phasic profile resulted from the V-NS0 model is very different from that of the pressure waveform extracted from the patient and the total variation of the pressure during a cardiac cycle is relatively small. The difference between the maximum and the minimum values of the pressure in a cycle is only approximately 20 mmHg at the inlet of the artery. The difference decreases along the axial length and the phasic waveform of the pressure becomes almost flat near the outlet of the boundary. It is obvious that the V-NS0 model can't provide correct information regarding the temporal distribution pattern of the pressure.
From Figure 5(1a&1b) we can see that the shape of the phasic profile of the velocity magnitude resulted from the P-NS0 model is not same as these of the other models. Instead, it is similar to that of the pressure waveform imposed at the inlet. The shape of the phasic profile of the pressure, the WPG, and the WSS for P-NS0 model are also all similar to that of the inlet pressure waveform. This is expected since the shape of the phasic blood flow profile follows the shape of the only pressure pulse waveform imposed at the boundary. The phasic velocity waveform of the P-NS0 model is clearly wrong compared with the velocity data extracted from the patient.
The shape of the phasic profiles of the PD, the WPG, the velocity, and the WSS for V-P model, P-P model, V-NS0 model, and P-V model are all similar to that of the velocity pulse waveform (Figure 5 and 7). This may suggest that when both the velocity and the pressure pulse waveforms are imposed at the inlet/outlet boundaries of the artery, the shape of the velocity pulse waveform dominates the temporal distribution patterns of both the velocity and the pressure difference related quantities, such as the PD, the WPG, velocity and the WSS. Even though there is no velocity waveform imposed directly on the boundaries for the P-P model, the reason for this observation still true is that the Pout applied at the outlet boundary was resulted from the computation of the P-V model, where the pulse waveform of Pout was already influenced by the velocity waveform.
BCs have no significant effect on spatial distribution patterns of blood flows
Plots in Figure 8 can further quantitatively confirm the observations on the spatial distribution patterns of the blood flow from Figure 9. The pressure has two sharper drops along the artery, each corresponding to the neck of a stenosis (see Figure 8(1a&1b)). These sharper drops of the wall pressure result in the elevated values of the spatial gradient of the wall pressure (see Figure 8(2a&2b)). The wall shear stress also peaks at the neck of each stenosis (see Figure 8(3a&3b)). Compared with the plots along the inner curve, the plots along the outer curve have some oscillations around the middle segment of the artery between the inlet and the outlet. It is probably due to the fact that the outer wall around the bend of the artery is less smooth. For the P-NS0 model (plots not included), the time averaged means of the pressure drop, the spatial wall pressure gradient, and the wall shear stress along the inner curve and along the outer curve have the similar patterns as those of other models, but with much larger values.
When and where the WSS is less than 1 N/m 2
Low and oscillating WSS have been associated with atherosclerotic lesions within the coronary arteries . Intimal thickening likely occurs when the mean WSS is below 1 N/m 2 , which presents an inverse hyperplasia with respect to the shear stress . It is of special interest to know the size of the region on the surface of the lumen where the WSS is lower than 1 N/m 2 and to examine the duration of the low WSS area in a cardiac cycle. At each point on the lumen surface, we calculated the total length of the time when the location experiences a low WSS < 1 N/m 2 in a cardiac cycle. We also recorded the distribution of the regions where the WSS is less than 1 N/m 2 at different time in the cardiac cycle and where the mean WSS is less than 1 N/m 2 .
The size of each region with the WSS under 1 N/m 2 varies depending on the time stage in a cardiac cycle. The maximum size of the regions is assumed at certain time in a cycle depending on the type of boundary conditions of the model. For the models imposed with velocity profile boundary condition either at the inlet or at the outlet (V-NS0, V-P, P-V), the maximum size of regions with the WSS under 1 N/m 2 occurs at the minimum flow rate (t = 0.3) (see plots in Figure 10(c) for V-P model). For the P-P model, the maximum size of the regions with the WSS under 1 N/m 2 occurs at the beginning of the cardiac cycle (t = 0) (see Figure 10(d)). For the P-NS0 model, there is no significant difference for the size of regions with the WSS under 1 N/m 2 at different time during the cycle.
Correlation of lumen cross-section area with the PD, the WPG, and the WSS
Correlations of the lumen cross section area with the slice averaged WSSmean, PDmean, and WPGmean.
Discussion and conclusion
In literature, the effect of the inflow boundary condition was investigated by examining the impact of the waveform and the shape of the spatial profile of the inlet velocity on the cardiac hemodynamics. Waveform dependence has been cited for flow in the aortic arch [18, 23] and in an end-to-side anastomosis . The study conducted by Myers et al. implies that waveform shape is not a significant factor in the right coronary artery time-averaged WSS features. Changes in the inlet velocity profiles did not produce significant changes in the arterial velocity and wall shear stress patterns . Our study was performed to examine the influence of the type of boundary conditions imposed at the inlet and the outlet on the phasic variation and the spatial distribution patterns of important right coronary artery hemodynamic features (velocity, the WSS, pressure, the PD, and the WPG) in a patient specific stenotic right coronary artery model.
Our results from this study indicate that the type of different boundary conditions imposed at the inlet and the outlet does not have significant effect on the spatial distribution patterns of the PD, the WPG and the WSS on the lumen surface of the artery in terms of where each quantity assumes maximum and where the minimum occurs. However, the temporal distribution patterns and the magnitudes of the velocity, the WSS, the pressure, the PD, and the WPG can be impacted notably depending on the type of the boundary conditions imposed at the inlet and the outlet. More in detail, it was observed that
The type of boundary conditions examined in this study does not have a significant effect on the time averaged mean PD, the mean WPG and the mean WSS on the lumen surface in the post-stenosis and the downstream regions. However, the boundary conditions have notable effect in the segment of the artery between the inlet and the first stenosis on the magnitude of the velocity at the center line, on the pressure, the PD, the WPG, and the WSS. The fully developed velocity profile assumption maybe the cause of the difference. Therefore, if the local hemodynamics is the focus of the study, the inlet spatial velocity profile must be handled with caution for V-P and V-NS0 models.
When the inlet boundary is imposed by a velocity waveform, the V-P and V-NS0 models with different outlet boundary conditions result in indistinguishable velocity, the PD, the WPG and the WSS. However, the temporal distribution patterns and the magnitudes of the pressure resulted from these two models are very different. The phasic variation of pressure in a cardiac cycle resulted from V-NS0 model is very different from that of the waveform extracted from the patient and thus may cause a misleading.
When both of the velocity and the pressure waveforms are imposed at inlet/outlet boundaries, the shapes of the phasic profiles of the PD, the WPG, and the WSS only follow the shape of the velocity pulse waveform no matter whether it is imposed at the inlet or the outlet boundary. This may suggest that the velocity boundary condition influents more than the pressure boundary condition on the phasic shape of the PD, the WPG, and the WSS.
P-NS0 model is not a practical mathematic model for the simulation of the blood flow in stenotic right coronary artery and it results in extremely large values of the PD and the WPD and a large variation of the WSS, which obviously are not realistic. It also results in a wrong shape of the pulse velocity waveform.
The analysis of the sensitivity of the blood flow simulation on the pressure outlet boundary condition was also carried out. For the model with a pulse velocity waveform imposed at the inlet boundary (V-P), a 1% mean perturbation of the pulse pressure waveform (3% at the peak flow time) imposed at the outlet boundary will only cause a change of 0.5% in the maximum mean pressure drop and no change in the maximum mean WPG and the maximum mean WSS. On the other hand, for a model with a pulse pressure waveform imposed at the inlet boundary (P-P), the same amount of perturbation on the pulse pressure waveform imposed at the outlet boundary will result in a 7% change of the maximum mean PD, a 7% change of the maximum mean WPG, and a 5% change of the maximum mean WSS and a 6% change of the minimum mean WSS, respectively. This may suggest that the magnitudes of these hemodynamic features are more sensitive to the outlet pressure boundary condition in the P-P model than to that in the V-P model. Therefore, the accuracy of the measurement of the outlet pressure data from the patient is more important to the blood simulation in a patient specific P-P model.
Comparing the plots of the inlet velocity and pressure waveforms with the plot of the pressure waveform at the outlet boundary, we can clearly see the influence of the flow rate on the pressure drop in a cardiac cycle. From the inlet to the outlet, the total pressure drop is larger at the time associated with each local maximum point of the inlet velocity waveform during a cardiac cycle, and the maximum pressure drop occurs at the peak flow (t = 0.425). A linear regression analysis on the correlation between the inlet pulse velocity magnitude and the percentage change of the pressure ((Pin(t)-Pout(t))/Pin(t)) reveals a correlation coefficient r = 0.9103 with p-value of 0.0000, and a linear fitting function: the percentage change of the pressure = 0.2134*velocity-0.0356 with R-square = 0.8286.
In conclusion, the observations based on the computer simulations in this study indicated that the ways how pressure and velocity boundary conditions are imposed in computational models have considerable impact on flow velocity and shear stress predictions, especially in upstream region. Computer simulations with improperly imposed type of boundary conditions may result in some misleading hemodynamic information, especially on the temporal distribution patterns and the magnitudes of the blood flow. Accuracy of in vivo measurements of blood pressure and velocity is of great importance for reliable model predictions. It must be recalled that we have examined only one artery in this work, and therefore any conclusions that we can draw are limited.
This work was partially supported by a grant from the Simons Foundation (#210082 to Biyue Liu) and by US NIH/NIBIB R01 EB004759.
Publication of this article was paid with funding from NIH/NIBIB R01 EB004759.
This article has been published as part of BioMedical Engineering OnLine Volume 14 Supplement 1, 2015: Cardiovascular Disease and Vulnerable Plaque Biomechanics. The full contents of the supplement are available online at http://www.biomedical-engineering-online.com/supplements/14/S1
- Myers JG, Moore JA, Ojha M, Johnston KW, Ethier CR: Factors influencing blood flow patterns in the human right coronary artery. Annals of Biomedical Engineering 2001, 29: 109–120.View ArticleGoogle Scholar
- Glagov S, Zarins CK, Giddens DP, Ku DN: Mechanical factors in the pathogenesis, localization and evolution of atherosclerotic plaques. In Diseases of the Arterial wall. Edited by: Camilleri B, Fiessinger B. Berlin: Springer-Verlag; 1989:217–39.View ArticleGoogle Scholar
- Samady H, Eshtehardi P, McDaniel MC, Suo J, Dhawan SS, Maynard C, Timmins LH, Quyyumi AA, Giddens DP: Coronary artery wall shear stress is associated with progression and transformation of atherosclerotic plaque and arterial remodeling in patients with coronary artery disease. Circulation 2011, 124: 779–788. 10.1161/CIRCULATIONAHA.111.021824View ArticleGoogle Scholar
- Torii R, Wood NB, Hadjiloizou N, Dowsey AW, Wright AR, et al.: Fluid-structure interaction analysis of a patient-specific right coronary artery with physiological velocity and pressure waveforms. Commun Numer Meth Engng 2009, 25: 565–580. 10.1002/cnm.1231MATHMathSciNetView ArticleGoogle Scholar
- Johnston BM, Johnston PR, Corney S, Kilpatrick D: Non-Newtonian blood flow in human right coronary arteries: transient simulations. J Biomech 2006, 39: 1116–1128. 10.1016/j.jbiomech.2005.01.034View ArticleGoogle Scholar
- Tang D, Kamm RD, Yang C, Zheng J, Canton G, Bach R, Huang X, Hatsukami TS, Zhu J, Ma G, Maehara A, Mintz GS, Yuan C: Image-based modeling for better understanding and assessment of atherosclerotic plaque progression and vulnerability: Data, modeling, validation, uncertainty and predictions. Journal of Biomechanics 2014, 47: 834–846. 10.1016/j.jbiomech.2014.01.012View ArticleGoogle Scholar
- Taylor CA, Humphrey JD: Open problems in computational vascular biomechanics: Hemodynamics and arterial wall mechanics. Comput Methods Appl Mech Engrg 2009, 198: 3514–3523. 10.1016/j.cma.2009.02.004MATHMathSciNetView ArticleGoogle Scholar
- Taylor CA, Steinman DA: Image-based modeling of blood flow and vessel wall dynamics: Applications, Methods and Future Directions. Annals of Biomedical Engineering 2010, 38: 1188–1203. 10.1007/s10439-010-9901-0View ArticleGoogle Scholar
- Zeng D, Ding Z, Friedman MH, Ethier CR: Effects of cardiac motion on right coronary artery hemodynamics. Annals of Biomedical Engineering 2003, 31: 420–429.View ArticleGoogle Scholar
- Konala BC, Ashish D, Rupak KB: Influence of arterial wall-stenosis compliance on the coronary diagnostic parameters. Journal of Biomechanics 2011, 44: 842–847. 10.1016/j.jbiomech.2010.12.011View ArticleGoogle Scholar
- Lee KW, Wood NB, Xu XY: Ultrasound image-based computer model of a common carotid artery with a plaque. Medical Engineering & Physics 2004, 26: 823–840. 10.1016/j.medengphy.2004.08.013View ArticleGoogle Scholar
- Liu B, Zheng J, Bach R, Tang D: Correlations of coronary plaque wall thickness with wall pressure and wall pressure gradient: a representative case study. BioMedical Engineering Online 2012, 11: 43. 10.1186/1475-925X-11-43View ArticleGoogle Scholar
- Ng EYK, Siauw WL: Unsteady viscous flow model on moving the domain through a stenotic artery. Proceedings of the Institution of Mechanical Engineers. Part H, Journal of engineering in medicine 2001, 215(2):237–49. 10.1243/0954411011533634View ArticleGoogle Scholar
- Nosovitsky VA, Ilegbusi OJ, Jiang J, Stone PH, Feldman CL: Effects of curvature and stenosis-like narrowing on wall shear stress in a coronary artery model with phasic flow. Computers and Biomedical Research 1997, 30: 61–82. 10.1006/cbmr.1997.1434View ArticleGoogle Scholar
- Perktold KR, Nerem M, Peter RO: A numerical calculation of flow in curved tube model of the left main coronary artery. J Biomech 1991, 24: 175–189. 10.1016/0021-9290(91)90176-NView ArticleGoogle Scholar
- Ethier CR, Steinman DA, Zhang X, Karpik SR, Ojha M: Flow waveform effects on end-to-side anastomotic flow patterns. Journal of Biomechanics 1998, 31: 609–617. 10.1016/S0021-9290(98)00059-1View ArticleGoogle Scholar
- Mark FF, Bargeron CB, Deters OJ, Friedman MH: Nonquasi-steady character of pulsatile flow in human coronary arteries. J Biomech Eng 1985, 107: 24–28. 10.1115/1.3138514View ArticleGoogle Scholar
- Morbiducci U, Ponzini R, Gallo D, Bignardi C, Rizzo G: Inflow boundary conditions for image-based computational hemodynamics: Impact of idealized versus measured velocity profiles in the human aorta. Journal of Biomechanics 2013, 46: 102–109. 10.1016/j.jbiomech.2012.10.012View ArticleGoogle Scholar
- Tang D, Yang C, Zheng J, Bach R, Wang L, Muccigrosso D, Billiar K, Zhu J, Ma G, Maehara A, Mintz GS, Fan R: Human coronary plaque wall thickness correlated positively with flow shear stress and negatively with plaque wall stress: an IVUS-based fluid-structure interaction multi-patient study. BioMedical Engineering OnLine 2014, 13: 32. 10.1186/1475-925X-13-32View ArticleGoogle Scholar
- Cho YI, Kensey KR: Effects of the non-Newtonian viscosity of blood on flows in a diseased arterial vessel. Part 1: steady flows. Biorheology 1991, 28: 241–262.Google Scholar
- Friedman MH, Deters OJ, Mark FF, Bargeron CB, Hutchins GM: Arterial geometry affects hemodynamics. A potential risk factor for atherosclerosis. Athero 1983, 46: 225–231. 10.1016/0021-9150(83)90113-2View ArticleGoogle Scholar
- Ku DN: Blood flow in arteries. J Ann Rev Fluid Mech 1997, 9: 399–434.MathSciNetView ArticleGoogle Scholar
- Chang L, Tarbell JM: Numerical simulation of fully developed sinusoidal and pulsatile (physiological) flow in curved tubes. J Fluid Mech 1985, 161: 175–198. 10.1017/S0022112085002877MATHMathSciNetView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.