Technical Brief

A Modification of Murray's Law for Shear-Thinning Rheology

[+] Author and Article Information
Patrick M. McGah

Department of Mechanical Engineering,
University of Washington,
Stevens Way, Box 352600,
Seattle, WA 98195
e-mail: pmcgah@u.washington.edu

Massimo Capobianchi

Professor Department of Mechanical Engineering, Gonzaga University,
502 East Boone Avenue,
Spokane, WA 99258
e-mail: capobianchi@gonzaga.edu

1Corresponding author.

Manuscript received September 5, 2014; final manuscript received December 29, 2014; published online March 10, 2015. Assoc. Editor: Hai-Chao Han.

J Biomech Eng 137(5), 054503 (May 01, 2015) (6 pages) Paper No: BIO-14-1440; doi: 10.1115/1.4029504 History: Received September 05, 2014; Revised December 29, 2014; Online March 10, 2015

This study reformulates Murray's well-known principle of minimum work as applied to the cardiovascular system to include the effects of the shear-thinning rheology of blood. The viscous behavior is described using the extended modified power law (EMPL), which is a time-independent, but shear-thinning rheological constitutive equation. The resulting minimization problem is solved numerically for typical parameter ranges. The non-Newtonian analysis still predicts the classical cubic diameter dependence of the volume flow rate and the cubic branching law. The current analysis also predicts a constant wall shear stress throughout the vascular tree, albeit with a numerical value about 15–25% higher than the Newtonian analysis. Thus, experimentally observed deviations from the cubic branching law or the predicted constant wall shear stress in the vasculature cannot likely be attributed to blood's shear-thinning behavior. Further differences between the predictions of the non-Newtonian and the Newtonian analyses are highlighted, and the limitations of the Newtonian analysis are discussed. Finally, the range and limits of applicability of the current results as applied to the human arterial tree are also discussed.

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.


Murray, C. D., 1926, “The Physiological Principle of Minimum Work: I. The Vascular System and the Cost of Blood Volume,” Proc. Natl. Acad. Sci. U.S.A., 12(3), pp. 207–214. [CrossRef] [PubMed]
Zamir, M., 1977, “Shear Forces and Blood Vessel Radii in the Cardiovascular System,” J. Gen. Physiol., 69(4), pp. 449–461. [CrossRef] [PubMed]
Zamir, M., and Medeiros, J., 1982, “Arterial Branching in Man and Monkey,” J. Gen. Physiol., 79(3), pp. 353–360. [CrossRef] [PubMed]
Zamir, M., Wrigley, S., and Langille, B., 1983, “Arterial Bifurcations in the Cardiovascular System of a Rat,” J. Gen. Physiol., 81(3), pp. 325–335. [CrossRef] [PubMed]
Rossitti, S., and Löfgren, J., 1993, “Vascular Dimensions of the Cerebral Arteries Follow the Principle of Minimum Work,” Stroke, 24(3), pp. 371–377. [CrossRef] [PubMed]
Kassab, G. S., and Fung, Y.-C. B., 1995, “The Pattern of Coronary Arteriolar Bifurcations and the Uniform Shear Hypothesis,” Ann. Biomed. Eng., 23(1), pp. 13–20. [CrossRef] [PubMed]
Taber, L. A., Ng, S., Quesnel, A. M., Whatman, J., and Carmen, C. J., 2001, “Investigating Murray's Law in the Chick Embryo,” J. Biomech., 34(1), pp. 121–124. [CrossRef] [PubMed]
Sather, B. A., Hageman, D., and Wagenseil, J. E., 2012, “Murray's Law in Elastin Haploinsufficient (Eln+/-) and Wild-Type (WT) Mice,” ASME J. Biomech. Eng., 134(12), p. 124504. [CrossRef]
Mayrovitz, H. N., and Roy, J., 1983, “Microvascular Blood Flow: Evidence Indicating a Cubic Dependence on Arteriolar Diameter,” Am. J. Physiol., 245(6), pp. H1031–H1038. [PubMed]
Riva, C. E., Grunwald, J. E., Sinclair, S. H., and Petrig, B., 1985, “Blood Velocity and Volumetric Flow Rate in Human Retinal Vessels,” Invest. Ophth. Vis. Sci., 26(8), pp. 1124–1132.
Taber, L. A., 1998, “An Optimization Principle for Vascular Radius Including the Effects of Smooth Muscle Tone,” Biophys. J., 74(1), pp. 109–114. [CrossRef] [PubMed]
Mayrovitz, H. N., 1987, “An Optimal Flow-Radius Equation for Microvessel Non-Newtonian Blood Flow,” Microvasc. Res., 34(3), pp. 380–384. [CrossRef] [PubMed]
LaBarbera, M., 1990, “Principles of Design of Fluid Transport Systems in Zoology,” Science, 249(4972), pp. 992–1000. [CrossRef] [PubMed]
Capobianchi, M., 2008, “Pressure Drop Predictions for Laminar Flows of Extended Modified Power Law Fluids in Rectangular Ducts,” Int. J. Heat Mass Transfer, 51(5), pp. 1393–1401. [CrossRef]
Revellin, R., Rousset, F., Baud, D., and Bonjour, J., 2009, “Extension of Murray's Law Using a Non-Newtonian Model of Blood Flow,” Theor. Biol. Med. Modell., 6(1), pp. 1–9. [CrossRef]
Chien, S., Usami, S., Taylor, H. M., Lundberg, J. L., and Gregersen, M. I., 1966, “Effects of Hematocrit and Plasma Proteins on Human Blood Rheology at Low Shear Rates,” J. Appl. Physiol., 21(1), pp. 81–87. [PubMed]
Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P., 2002, Numerical Recipes in C: The Art of Scientific Computing, 2nd ed., Cambridge University Press, New York.
Pries, A., Neuhaus, D., and Gaehtgens, P., 1992, “Blood Viscosity in Tube Flow: Dependence on Diameter and Hematocrit,” Am. J. Physiol., 263(6 Pt.2), pp. H1770–H1778. [PubMed]
Alarcón, T., Byrne, H. M., and Maini, P. K., 2005, “A Design Principle for Vascular Beds: The Effects of Complex Blood Rheology,” Microvasc. Res., 69(3), pp. 156–172. [CrossRef] [PubMed]
Gaenzer, H., Neumayr, G., Marschang, P., Sturm, W., Kirchmair, R., and Patsch, J. R., 2001, “Flow-Mediated Vasodilation of the Femoral and Brachial Artery Induced by Exercise in Healthy Nonsmoking and Smoking Men,” J. Am. Coll. Cardiol., 38(5), pp. 1313–1319. [CrossRef] [PubMed]
Mitchell, G. F., Parise, H., Vita, J. A., Larson, M. G., Warner, E., Keaney, J. F., Jr., Keyes, M. J., Levy, D., Vasan, R. S., and Benjamin, E. J., 2004, “Local Shear Stress and Brachial Artery Flow-Mediated Dilation: The Framingham Heart Study,” Hypertension, 44(2), pp. 134–139. [CrossRef] [PubMed]
Loth, F., Fischer, P. F., and Bassiouny, H. S., 2008, “Blood Flow in End-to-Side Anastomoses,” Ann. Rev. Fluid Mech., 40, pp. 367–393. [CrossRef]
Cheng, C., Helderman, F., Tempel, D., Segers, D., Hierck, B., Poelmann, R., van Tol, A., Duncker, D. J., Robbers-Visser, D., Ursem, N. T., van Haperen, R., Wentzel, J. J., Gijsen, F., van der Steen, A. F., de Crom, R., and Krams, R., 2007, “Large Variations in Absolute Wall Shear Stress Levels Within one Species and Between Species,” Atherosclerosis, 195(2), pp. 225–235. [CrossRef] [PubMed]
Dammers, R., Stifft, F., Tordoir, J. H. M., Hammeleers, J. M. M., Hoeks, A. P. G., and Kitslaar, P. J. E. H. M., 2003, “Shear Stress Depends on Vascular Territory: Comparison Between Common Carotid and Brachial Artery,” J. Appl. Physiol., 94(2), pp. 458–489. [CrossRef]
Wu, S. P., Ringgaard, S., Oyre, S., Hansen, M. S., Rasmus, S., and Pedersen, E. M., 2004, “Wall Shear Rates Differ Between the Normal Carotid, Femoral, and Brachial Arteries: An In Vivo MRI Study,” J. Magn. Reson. Imaging, 19(2), pp. 188–193. [CrossRef] [PubMed]
Zamir, M., Sinclair, P., and Wonnacott, T. H., 1992, “Relation Between Diameter and Flow in Major Branches of the Arch of the Aorta,” J. Biomech., 25(11), pp. 1303–1310. [CrossRef] [PubMed]
Liu, Y., and Kassab, G. S., 2007, “Vascular Metabolic Dissipation in Murray's Law,” Am. J. Physiol. Heart Circ. Physiol., 292(3), pp. H1336–H1339. [CrossRef] [PubMed]
Huo, Y., and Kassab, G. S., 2012, “Intraspecific Scaling Laws of Vascular Trees,” J. Roy Soc. Interface, 9(66), pp. 190–200. [CrossRef]
Brewster, R. A., 2013, “Pressure Drop Predictions for Laminar Fully-Developed Flows of Purely-Viscous Non-Newtonian Fluids in Circular Ducts,” ASME J. Fluid Eng., 135(10), p. 101106. [CrossRef]


Grahic Jump Location
Fig. 1

Viscosity versus shear rate for the least-squares fit of EMPL model (solid line) to the experimental data (symbols) of Chien et al. [16]

Grahic Jump Location
Fig. 2

Cost function versus radius for a flow rate of 6 mL/min using b = 778 erg/cm3/s for both Newtonian and non-Newtonian cases

Grahic Jump Location
Fig. 3

Volume flow rate versus optimal radius for a non-Newtonian fluid in the case where b = 778 erg/cm3/s. The slope of the curve is equal to 3 to within a margin 10−12, indicating that Murray's cubic relationship is valid even when non-Newtonian behavior is considered.

Grahic Jump Location
Fig. 4

Wall shear stress versus optimal radius using b = 778 erg/cm3/s for both Newtonian and non-Newtonian cases

Grahic Jump Location
Fig. 5

Velocity profiles at the optimal radius using b = 778 erg/cm3/s for both Newtonian and non-Newtonian cases



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In