- 1PKU-HKUST Shenzhen-Hongkong Institution, Shenzhen, China
- 2Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing, China
- 3California Medical Innovations Institute, San Diego, CA, United States
A pulsatile hemodynamic analysis was carried out in the right coronary arterial (RCA) tree of control and RV hypertrophy (RVH) hearts. The shape of flow and wall shear stress (WSS) waves was hypothesized to be maintained throughout the RCA tree in RVH (i.e., similar patterns of normalized flow and WSS waves in vessels of various sizes). Consequently, we reconstructed the entire RCA tree down to the first capillary bifurcation of control and RVH hearts based on measured morphometric data. A Womersley-type model was used to compute the flow and WSS waves in the tree. The hemodynamic parameters obtained from experimental measurements were incorporated into the numerical model. Given an increased number of arterioles, the mean and amplitude of flow waves at the inlet of RCA tree in RVH was found to be two times larger than that in control, but no significant differences (p > 0.05) were found in precapillary arterioles. The increase of stiffness in RCA of RVH preserved the shape of normalized flow and WSS waves, but increased the PWV in coronary arteries and reduced the phase angle difference for the waves between the most proximal RCA and the most distal precapillary arteriole. The study is important for understanding pulsatile coronary blood flow in ventricular hypertrophy.
Introduction
There is compensatory vascular remodeling that accompanies RV hypertrophy (RVH) (Cooper et al., 1981; Manohar et al., 1981; Botham et al., 1984; Manohar, 1985; White et al., 1992; Kassab et al., 1993). In porcine model, systemic blood pressure in right coronary artery (RCA) was unchanged, but an increase of blood flow occurred in large epicardial branches at 5 weeks after pulmonary hypertension (Lu et al., 2011). Based on the measured morphometric data in arrested, vasodilated porcine heart of RVH, we carried out a steady-state flow analysis of arterial tree down to first capillary segments (Huo et al., 2007). The increase of blood flow was found to be caused by the compensatory growth of small vessels, which resulted in restoration of blood flow and wall shear stress (WSS) to normal level in the perfusion arterioles (diameter <100 μm) of RCA tree in the diastolic state of RVH hearts (Huo et al., 2007). The pulsatility of coronary blood flow is a significant hemodynamic feature (Fung, 1997; Nichols and McDonald, 2011). However, the pulsatile pressure-flow relationship has not been investigated in the coronary circulation during the progression of RVH. We have simulated the pulsatile blood flow in diastole in the absence of vessel tone in the entire coronary arterial tree of normal porcine hearts using the Womersley-type mathematical model (Huo and Kassab, 2006). Hence, the objective of present study is to carry out a complete pulsatile flow analysis in the coronary arterial tree of control and RVH hearts and to determine the effect of RVH on the flow and WSS waves in diastole in the absence of vessel tone.
Here, we hypothesized that the flow and WSS waves was preserved in RVH (i.e., similar patterns of normalized flow and WSS waves in vessels of various sizes). To test the hypothesis, pulsatile blood flows were computed by a Womersley-type numerical model (Huo and Kassab, 2006) in each vessel of the entire RCA tree down to the first capillary bifurcation (excluding the sub-tree distal to the posterior descending artery) of control and RVH hearts in diastole, which was reconstructed from morphometric data (i.e., vessel diameters, lengths and numbers) (Kassab et al., 1993). The experimental measurements of coronary wall thickness (Guo and Kassab, 2004) and stiffness (Garcia and Kassab, 2009) were also incorporated into the numerical model. The constitutive equation, based on experimental measurements, was similar to a previous study (Huo and Kassab, 2006). The predictions of the mathematical model showed good agreement with the experimental measurements in control and RVH hearts. A detailed comparison of pulsatile blood flows was made in vessels of various sizes throughout the entire RCA arterial trees between control and RVH.
Materials and Methods
Anatomical Model
Previously, Kassab et al. have carried out morphometric measurements of the RCA-posterior descending arterial (PDA) trees of control and RVH hearts (Kassab et al., 1993). Briefly, the morphometric data on the coronary arterial vessels of diameters <40 μm were obtained from histological specimens and the morphometric data on the coronary arterial vessels of diameters >40 μm were obtained from cast studies. The entire RCA-PDA tree down to the first capillary bifurcations was reconstructed in control and RVH using a growth algorithm (Mittal et al., 2005), based on the experimental measurements of the RV branches excluding the distal tree to the PDA. In summary, the present anatomical mathematical model has exact data (diameters, lengths and connectivity) for the larger vessels and statistical data for the microvessels different from the previous models (Kaimovitz et al., 2005, 2010), which were based on the statistical data (tables with means and standard deviations for diameters, lengths, connectivity) for the entire tree (Kassab et al., 1993).
Flow Simulation
After the branching pattern and vascular geometry of RV branches were generated, a pulsatile flow analysis was performed similar to a previous study (Huo and Kassab, 2006, 2007). Briefly, in the frequency domain, the governing equations (transformed from the conversion of mass and momentum) for flow (Q) and pressure (P) in a vessel are written as:
Where a and b are arbitrary constants of integration, x the axial coordinate along the vessel, ω the angular frequency, () is the wave velocity, h/R the ratio of wall thickness to radius, E the Young's modulus, ρ the density, and F10(α)= (, μ is the dynamic viscosity, J0 the Bessel function of zero order and first kind, and J1 the Bessel function of first order and first kind). (A(n) is the cross-sectional area in a vessel) is defined as the characteristic admittance, Z0 = 1/Y0 the characteristic impedance, , and . The impedance and admittance in a vessel is:
In a given vessel segment, at x = 0 and x = L, we have the following inlet and outlet impedance:
A combination of Equations (5) and (6) yields:
Since there are two or more vessels that emanate from the junction points of the entire RCA tree, the junction boundary condition (determined from the continuous pressure and mass conservation at the junction) is written as:
Equations (7) and (8) were used to calculate the impedance/admittance in the entire coronary tree from inlet to the capillary vessels.
The terminal impedance/admittance of the first capillary was assumed to be equal to the steady value as (g·sec/cm4), from which we proceeded backwards to iteratively calculate the impedance/admittance in the entire RCA tree of control and RVH hearts using Equations (7) and (8). The pulsatile pressure was the same as the previous study (Huo and Kassab, 2006) and discretized by a Fourier transformation as the inlet boundary condition. The flow and pressure in each vessel were then calculated by using Equations (1) and (2) coupled with the continuous pressure at junctions. The blood flow density (ρ) was assumed to be 1.06 g/cm3. The variation of viscosity with vessel diameter and hematocrit was based on Pries' viscosity model (Pries et al., 1992). The coronary wall thickness for every order was adopted from the previous measurements (Guo and Kassab, 2004; Choy and Kassab, 2009). The in situ static Young's modulus in control was made as ~8.0 × 106 (dynes/cm2) based on the experimental data (Kassab and Molloi, 2001), which was doubled in RVH (Garcia and Kassab, 2009). The dynamic Young's modulus was also considered for various frequencies (Bergel, 1961; Douglas and Greenfield, 1970; Gow et al., 1974), i.e., the Young's modulus increases with the increase of frequency ω (see Figure 2 in Bergel, 1961). Once the flow wave was determined in each vessel, the WSS waves, τ, was calculated as Zheng et al. (2010):
where Λ2 = i3α2. Unless otherwise stated, all computations used the previously measured physical properties and parameters as described above.
Statistical Analysis
ANOVA (SigmaStat 3.5) was used to detect statistical differences between control and RVH. A p < 0.05 was indicative of a significant difference between the two populations.
Results
The pulsatile model has been previously validated experimentally in normal hearts. Here, good agreement was found between experimental measurements and computational results for the RCA of RVH hearts, as shown in Figure 1. The steady-state flows were measured under loading and unloading of pressures as the RCA tree was perfused by cardioplegic solution, which were consistent with the computed pulsatile flows with the frequency approaching to zero since a steady-state flow can be mimicked by a pulsatile flow as ω → 0.
Figure 1. Experimental and computed pressure-flow relationship of the RCA of control and RVH hearts. The experimental results were measured under loading and unloading of pressures as the RCA tree was perfused by cardioplegic solution (μ = 1.1 cp and ρ = 1 g/cm3). The pulsatile model was used to compute the pressure and flow at very low frequency (ω → 0).
The flow waves were calculated in each vessel of RCA tree (excluding the distal tree to the PDA), which has RV branches with a mean (averaged over five anatomic reconstructions) of 0.36 and 1.3 million vessels for control and RVH, respectively. Table 1 summarizes morphometric and hemodynamic parameters in diameter-defined Strahler orders from precapillary arterioles (order 1) to the epicardial RCA tree (the highest order) in control and RVH hearts. The flow waves at the inlet and primary branches of RCA tree of control heart were compared with those in the RVH heart (Figure 2 vs. Figure 3). Given such an increase of vessel numbers, the mean (i.e., time-averaged value over a cardiac cycle) and amplitude (i.e., the change between peak and trough) of flow wave at the inlet of RCA tree in RVH is much larger than that in control. Figures 4A,B show the relationship between the time-averaged flow in a vessel and the cumulative length of the vessel from the RCA to the precapillary arteriole through similar primary branches in control and RVH hearts, respectively. Accordingly, Figures 4C,D show the relationship between the time-averaged pressure in a vessel and the cumulative length of the vessel. Figures 5A,B show the amplitude and phase angle of the impedance in the most proximal RCA, a distal vessel (6 cm from the RCA), and the most distal precapillary arteriole in control and RVH hearts, respectively. Accordingly, Figures 5C,D show a decrease of flow waves sequentially along the path from the RCA to the precapillary arteriole.
Table 1. Morphometric and hemodynamic parameters in orders according to the diameter-defined Strahler system.
Figure 2. Flow waves at the inlet and primary branches (P1–P6) of RCA tree of a control heart. The dimensions in the parenthesis refer to diameter and length (μm).
Figure 3. Flow waves at the inlet and primary branches (P1–P6) of RCA tree of a RVH heart. The dimensions in the parenthesis refer to diameter and length (μm).
Figure 4. (A,B) Relationship between the mean flow (averaged over a cardiac cycle) in a vessel and the cumulative length of the vessel from the root to the precapillary arteriole through similar primary branches in: (A) control and (B) RVH pig hearts (i.e., P2 in Figure 2 and P1 in Figure 3, respectively). (C,D) Relationship between the mean pressure (averaged over a cardiac cycle) in a vessel and the cumulative length of the vessel corresponding to (A,B).
Figure 5. (A,B) Impedance |Z(0,ω)| vs. angular frequency ω and phase Z(0,ω) vs. angular frequency ω sequentially along the corresponding paths to Figures 4A,B. (C,D) Flow waves sequentially along the corresponding paths to Figures 4A,B. The arrows indicate the one-to-one correspondence between (A–D).
Despite the large changes of flow waves from the root to the precapillary, we previously reported similar pattern of normalized waves in control hearts. Figures 6A,B show the flow and WSS waves normalized by the mean values at the most proximal RCA and the most distal precapillary arteriole of RCA tree in control. Figures 6C,D show the normalized flow and WSS waves in RVH as compared with those in control. The normalized waves at the RCA and precapillary arteriole are similar to each other (p >> 0.05) in both control and RVH. In comparison with the control, the difference of phase angles for flow waves between RCA and precapillary arteriole decreased by about 50% in various frequencies and the pulse wave velocity (PWV) increased by about 40% in the hypertrophic RCA tree.
Figure 6. (A) Flow and (B) WSS waves normalized by the mean values (averaged over a cardiac cycle) at the RCA and precapillary arteriole of control heart. (C,D) Normalized flow and WSS waves at the RCA and precapillary arteriole of RVH heart in correspondence with (A,B), respectively.
The effects of vessel compliance on the pattern of flow and WSS waves were examined. Figures 7A,B show a sensitivity analysis of normalized flow and WSS waves, respectively, at the inlet of RCA tree of control hearts with a 50% increase/decrease of Young's modulus of each vessel wall (including both static and dynamic Young's moduli). Figure 7C shows the change of the phase angle of the inlet impedance with a 50% increase/decrease of Young's modulus of each vessel wall. The change of vessel compliance in physiological range had negligible effects on the normalized waves (p >> 0.05), but significantly affected the phase angle of the impedance.
Figure 7. (A,B) A sensitivity analysis of normalized (A) flow and (B) WSS waves at the inlet of RCA tree with a 50% increase/decrease of Young's modulus of each vessel wall, where the baseline (i.e., the unity) refers to the static Young's modulus of 8 × 106 dynes/cm2 with adjustment for frequency. (C) Phase Z(0,ω) vs. angular frequency ω at the inlet of RCA tree with a 50% increase/decrease of Young's modulus of each vessel wall.
Discussion
The major finding was that RVH maintained similar patterns of normalized flow and WSS waves in different size vessels. RVH, however, did significant increase in the amplitude and mean (time-averaged over a cardiac cycle) of waves in large epicardial coronary branches due to the increase in number of small arterioles.
Steady-State Flow vs. Pulsatile Flow
Based on morphometric measurements and steady-state hemodynamic analysis, the remodeling of structure and function of the entire RCA tree was determined in RVH after 5 weeks of pulmonary banding (Kassab et al., 1993; Huo et al., 2007). The number of vessel segments in RVH was about four times larger than that in control hearts. The steady-state flow (which approximately equaled to the time-averaged value of pulsatile flow) was significantly increased in the large epicardial vessels in RVH, but distributed through substantially more vessels such that the flow and WSS in the small perfusion arterioles were not statistically different from the control (Huo et al., 2007). Here, we compared pulsatile flows between control and RVH using the Womersley-type model (Huo and Kassab, 2006) and morphometric measurements (Kassab et al., 1993; Kassab and Molloi, 2001; Garcia and Kassab, 2009). Corresponding to the steady-state flow, the amplitudes of flow and WSS waves in RVH were significantly increased in the epicardial vessels, but similar to those in the perfusion arterioles of control, as shown in Table 1 and Figures 2, 3. The amplitude and mean of flow and pressure waves decreased along the path from the root down to the precapillary arteriole in the RVH, which has the same trend as the control (see Figure 4) and also agrees with a previous study (Huo and Kassab, 2006). The increase of stiffness in RVH reduces the phase angle of vessel impedance (see Figure 5) and leads to a significant decrease (~50%) of phase angle difference for flow waves between the most proximal RCA and the most distal precapillary arteriole in various frequencies (see Figure 6), which reflects an about 40% increase of PWV in the RCA tree of RVH. Although an increase of aortic PWV is widely known as a marker of cardiovascular risk in hypertension (Mohiaddin et al., 1993; Butlin et al., 2013), it still requires further investigation in the interaction of large and small arteries, given the unknown causal relation between stiffening and hypertension in vessels of various sizes.
Flow and WSS Waves
The flow waves at various coronary vessels of diastolic control hearts were found to have a tendency of scaling to a single curve, except for a small phrase angle difference, when they were normalized by the mean values (Huo and Kassab, 2007). In the beating hearts, Ashikawa et al. measured the flow velocities in small arterioles (diameter of 12.8 ± 4.1 μm), capillaries, and small venules (diameter of 16.5 ± 6.5 μm; Ashikawa et al., 1986) and Toyota et al. measured the flow velocities in arterioles (diameter of ~100 μm; Toyota et al., 2005), which also showed similar normalized flow waves. Here, the normalized flow and WSS waves in the RCA tree of RVH maintained the same scaling characteristic, as shown in Figures 5, 6. Moreover, the sensitivity analysis in Figures 7A,B showed the change of stiffness of RCA tree in physiological range, which had negligible effect on the normalized flow and WSS waves. This can be explained from the constitutive equation [i.e., , in which A is the CSA, p the pressure, h/R the ratio of wall thickness to radius, E the Young's modulus]. From the equation, Young's modulus (the denominator) is almost twenty times larger than (p−pref)·(R/h) (the numerator) in normal coronary arterial tree such that is much less than one, where the static Young's modulus in control is 8 × 106 dynes/cm2 and increases as frequency increases, (p−pref) varies from −20 to +20 mmHg, and (R/h) approximately equals to ten. A 50% perturbation for the stiffness had negligible effect on the waves.
Physiological Implications
Coronary blood vessel wall is subjected to various types of hemodynamic forces (e.g., hydrostatic pressure, cyclic stretch, and fluid shear stress) caused by the pulsatile blood pressure and flow (Chiu and Chien, 2011). Intraspecific scaling power laws of vascular trees, derived from the steady-state analysis (Huo and Kassab, 2012), characterize coronary vasculature with remarkable simplicity. The compensatory remodeling in RVH was found to maintain the structure-function hierarchy (preserved scaling exponents) of fractal-like coronary arterial tree (Gong et al., 2016). Furthermore, the present study shows the preserved shape of normalized wave of pulsatile blood flows in the entire coronary arterial tree, which is unchanged during the remodeling in RVH. This implies that intraspecific scaling power laws of vascular trees can be extended from the steady to dynamic states. On the other hand, the altered phase angle of vessel impedance and flow and WSS waves shows the effects of the remodeling on the pulsatile blood pressure and flow in RVH, which requires further investigations.
Critique of Study
The present study carried out the pulsatile flow analysis based on the morphometric data and physiological measurements, but the constitutive equation was determined in a vasodilated state in the absence of vasomotor tone where coronary flow reserve was substantially reduced. The stiffness of large arteries was assumed to equal to that of small arterioles. Future studies are needed to identify the contribution of vascular tone to the vascular remodeling and flow waves in RVH with consideration of the stiffness in various size vessels.
Significance of Study
RVH is a compensatory response to pulmonary hypertension and the compensatory adaptations of coronary circulation serve to maintain perfusion of the increased myocardial mass. Here, we showed that the shape of normalized flow and WSS waves was preserved despite a significant increase of the amplitude and mean of flow and WSS waves in the epicardial vessels. The increase of stiffness in RVH increased the PWV in coronary arteries and reduced the phase angle difference for flow waves between the most proximal RCA and the most distal precapillary arteriole, but had negligible effect on the pattern of flow and WSS waves. The precise prediction of flow and WSS waves in coronary microcirculation is physiologically and clinically important to understand the ventricular hypertrophy.
Author Contributions
Data analysis were performed by YH at the college of engineering, Peking University. Paper was drafted and revised by YH and GK at the California Medical Innovations Institute.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This research is supported in part by the National Institute of Health-National Heart, Lung, and Blood Institute Grant HL134841 and HL118738 (GK) and the National Natural Science Foundation of China Grant 11672006 and Shenzhen Science and Technology R&D Grant JCYJ20160427170536358 (YH).
References
Ashikawa, K., Kanatsuka, H., Suzuki, T., and Takishima, T. (1986). Phasic blood flow velocity pattern in epimyocardial microvessels in the beating canine left ventricle. Circ. Res. 59, 704–711. doi: 10.1161/01.RES.59.6.704
Bergel, D. H. (1961). The dynamic elastic properties of the arterial wall. J. Physiol. (Lond). 156, 458–469. doi: 10.1113/jphysiol.1961.sp006687
Botham, M. J., Lemmer, J. H., Gerren, R. A., Long, R. W., Behrendt, D. M., and Gallagher, K. P. (1984). Coronary vasodilator reserve in young dogs with moderate right ventricular hypertrophy. Ann. Thorac. Surg. 38, 101–107. doi: 10.1016/S0003-4975(10)62214-X
Butlin, M., Qasem, A., Battista, F., Bozec, E., McEniery, C. M., Millet-Amaury, E., et al. (2013). Carotid-femoral pulse wave velocity assessment using novel cuff-based techniques: comparison with tonometric measurement. J Hypertens. 31, 2237–2243, discussion 2243. doi: 10.1097/HJH.0b013e328363c789
Chiu, J. J., and Chien, S. (2011). Effects of disturbed flow on vascular endothelium: pathophysiological basis and clinical perspectives. Physiol. Rev. 91, 327–387. doi: 10.1152/physrev.00047.2009
Choy, J. S., and Kassab, G. S. (2009). Wall thickness of coronary vessels varies transmurally in the LV but not the RV: implications for local stress distribution. Am. J. Physiol. Heart Circ. Physiol. 297, H750–H758. doi: 10.1152/ajpheart.01136.2008
Cooper, G. T., Tomanek, R. J., Ehrhardt, J. C., and Marcus, M. L. (1981). Chronic progressive pressure overload of the cat right ventricle. Circ. Res. 48, 488–497. doi: 10.1161/01.RES.48.4.488
Douglas, J. E., and Greenfield, J. C. Jr. (1970). Epicardial coronary artery compliance in the dog. Circ. Res. 27, 921–929. doi: 10.1161/01.RES.27.6.921
Garcia, M., and Kassab, G. S. (2009). Right coronary artery becomes stiffer with increase in elastin and collagen in right ventricular hypertrophy. J. Appl. Physiol. (1985). 106, 1338–1346. doi: 10.1152/japplphysiol.90592.2008
Gong, Y., Feng, Y., Chen, X., Tan, W., Huo, Y., and Kassab, G. S. (2016). Intraspecific scaling laws are preserved in ventricular hypertrophy but not in heart failure. Am. J. Physiol. Heart Circ. Physiol. 311, H1108–H1117. doi: 10.1152/ajpheart.00084.2016
Gow, B. S., Schonfeld, D., and Patel, D. J. (1974). The dynamic elastic properties of the canine left circumflex coronary artery. J. Biomech. 7, 389–395. doi: 10.1016/0021-9290(74)90001-3
Guo, X., and Kassab, G. S. (2004). Distribution of stress and strain along the porcine aorta and coronary arterial tree. Am. J. Physiol. Heart Circ. Physiol. 286, H2361–H2368. doi: 10.1152/ajpheart.01079.2003
Huo, Y., and Kassab, G. S. (2006). Pulsatile blood flow in the entire coronary arterial tree: theory and experiment. Am. J. Physiol. Heart Circ. Physiol. 291, H1074–H1087. doi: 10.1152/ajpheart.00200.2006
Huo, Y., and Kassab, G. S. (2007). A hybrid one-dimensional/Womersley model of pulsatile blood flow in the entire coronary arterial tree. Am. J. Physiol. Heart Circ. Physiol. 292, H2623–H2633. doi: 10.1152/ajpheart.00987.2006
Huo, Y., and Kassab, G. S. (2012). Intraspecific scaling laws of vascular trees. J. R. Soc. Interface 9, 190–200. doi: 10.1098/rsif.2011.0270
Huo, Y., Linares, C. O., and Kassab, G. S. (2007). Capillary perfusion and wall shear stress are restored in the coronary circulation of hypertrophic right ventricle. Circ. Res. 100, 273–283. doi: 10.1161/01.RES.0000257777.83431.13
Kaimovitz, B., Lanir, Y., and Kassab, G. S. (2005). Large-scale 3-D geometric reconstruction of the porcine coronary arterial vasculature based on detailed anatomical data. Ann. Biomed. Eng. 33, 1517–1535. doi: 10.1007/s10439-005-7544-3
Kaimovitz, B., Lanir, Y., and Kassab, G. S. (2010). A full 3-D reconstruction of the entire porcine coronary vasculature. Am. J. Physiol. Heart Circ. Physiol. 299, H1064–H1076. doi: 10.1152/ajpheart.00151.2010
Kassab, G. S., Imoto, K., White, F. C., Rider, C. A., Fung, Y. C., and Bloor, C. M. (1993). Coronary arterial tree remodeling in right ventricular hypertrophy. Am. J. Physiol. 265, H366–H375. doi: 10.1152/ajpheart.1993.265.1.H366
Kassab, G. S., and Molloi, S. (2001). Cross-sectional area and volume compliance of porcine left coronary arteries. Am. J. Physiol. Heart Circ. Physiol. 281, H623–H628. doi: 10.1152/ajpheart.2001.281.2.H623
Lu, X., Dang, C. Q., Guo, X., Molloi, S., Wassall, C. D., Kemple, M. D., et al. (2011). Elevated oxidative stress and endothelial dysfunction in right coronary artery of right ventricular hypertrophy. J. Appl. Physiol. (1985). 110, 1674–1681. doi: 10.1152/japplphysiol.00744.2009
Manohar, M. (1985). Transmural coronary vasodilator reserve, and flow distribution during tachycardia in conscious young swine with right ventricular hypertrophy. Cardiovasc. Res. 19, 104–112. doi: 10.1093/cvr/19.2.104
Manohar, M., Bisgard, G. E., Bullard, V., and Rankin, J. H. (1981). Blood flow in the hypertrophied right ventricular myocardium of unanesthetized ponies. Am. J. Physiol. 240, H881–H888. doi: 10.1152/ajpheart.1981.240.6.H881
Mittal, N., Zhou, Y., Ung, S., Linares, C., Molloi, S., and Kassab, G. S. (2005). A computer reconstruction of the entire coronary arterial tree based on detailed morphometric data. Ann. Biomed. Eng. 33, 1015–1026. doi: 10.1007/s10439-005-5758-z
Mohiaddin, R. H., Firmin, D. N., and Longmore, D. B. (1993). Age-related changes of human aortic flow wave velocity measured noninvasively by magnetic resonance imaging. J. Appl. Physiol. (1985). 74, 492–497. doi: 10.1152/jappl.1993.74.1.492
Nichols, W. W., and McDonald, D. A. (2011). McDonald's Blood Flow In Arteries: Theoretic, Experimental, and Clinical Principles, 6th Edn. London: Hodder Arnold.
Pries, A. R., Neuhaus, D., and Gaehtgens, P. (1992). Blood viscosity in tube flow: dependence on diameter and hematocrit. Am. J. Physiol. 263, H1770–H1778. doi: 10.1152/ajpheart.1992.263.6.H1770
Toyota, E., Ogasawara, Y., Hiramatsu, O., Tachibana, H., Kajiya, F., Yamamori, S., et al. (2005). Dynamics of flow velocities in endocardial and epicardial coronary arterioles. Am. J. Physiol. Heart Circ. Physiol. 288, H1598–H1603. doi: 10.1152/ajpheart.01103.2003
White, F. C., Nakatani, Y., Nimmo, L., and Bloor, C. M. (1992). Compensatory angiogenesis during progressive right ventricular hypertrophy. Am. J. Cardiovasc. Pathol. 4, 51–68.
Keywords: pulsatile flow, right ventricular hypertrophy, right coronary arterial tree, Womersley-type model, Pulsatile wall shear stress
Citation: Huo Y and Kassab GS (2018) Coronary Blood Flow Is Increased in RV Hypertrophy, but the Shape of Normalized Waves Is Preserved Throughout the Arterial Tree. Front. Physiol. 9:675. doi: 10.3389/fphys.2018.00675
Received: 08 December 2017; Accepted: 15 May 2018;
Published: 31 May 2018.
Edited by:
Timothy W. Secomb, University of Arizona, United StatesReviewed by:
C. Alberto Figueroa, University of Michigan, United StatesDaniel Goldman, University of Western Ontario, Canada
Copyright © 2018 Huo and Kassab. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ghassan S. Kassab, Z2thc3NhYkBjYWxtaTIub3Jn