Extension of Murray's law using a non-Newtonian model of blood flow

1742-4682-6-7 1742-4682 Research Extension of Murray's law using a non-Newtonian model of blood flow Revellin Rémi remi.revellin@insa-lyon.fr Rousset François francois.rousset@insa-lyon.fr Baud David David.Baud@chuv.ch Bonjour Jocelyn jocelyn.bonjour@insa-lyon.fr

Université de Lyon, CNRS, INSA-Lyon, CETHIL, UMR5008, F-69621, Villeurbanne, France, Université Lyon 1, F-69622, France

Department of Obstetrics and Gynaecology, University Hospital of Lausanne, Maternity-CHUV, CH-1011 Lausanne, Switzerland

Theoretical Biology and Medical Modelling 1742-4682 2009 6 1 7 http://www.tbiomed.com/content/6/1/7 19445663 10.1186/1742-4682-6-7 09 4 2009 15 5 2009 15 5 2009 2009 Revellin et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

So far, none of the existing methods on Murray's law deal with the non-Newtonian behavior of blood flow although the non-Newtonian approach for blood flow modelling looks more accurate.

Modeling

In the present paper, Murray's law which is applicable to an arterial bifurcation, is generalized to a non-Newtonian blood flow model (power-law model). When the vessel size reaches the capillary limitation, blood can be modeled using a non-Newtonian constitutive equation. It is assumed two different constraints in addition to the pumping power: the volume constraint or the surface constraint (related to the internal surface of the vessel). For a seek of generality, the relationships are given for an arbitrary number of daughter vessels. It is shown that for a cost function including the volume constraint, classical Murray's law remains valid (i.e. ΣR^c= cste with c = 3 is verified and is independent of n, the dimensionless index in the viscosity equation; R being the radius of the vessel). On the contrary, for a cost function including the surface constraint, different values of c may be calculated depending on the value of n.

Results

We find that c varies for blood from 2.42 to 3 depending on the constraint and the fluid properties. For the Newtonian model, the surface constraint leads to c = 2.5. The cost function (based on the surface constraint) can be related to entropy generation, by dividing it by the temperature.

Conclusion

It is demonstrated that the entropy generated in all the daughter vessels is greater than the entropy generated in the parent vessel. Furthermore, it is shown that the difference of entropy generation between the parent and daughter vessels is smaller for a non-Newtonian fluid than for a Newtonian fluid.

Introduction

Since several decades, many studies have been carried out on the optimal branching pattern of a vascular system. Based on the simple assumption of a steady Poiseuille blood flow, the well known Murray's law 1 has been established. It links the radius of a parent vessel R₀(immediately upstream from a vessel bifurcation) to the radii of the daughter vessels R₁and R₂(immediately downstream after a vessel bifurcation) as R₀/R₁= R₀/R₂= 2^-1/3. From Murray's analysis, the required condition of minimum power occurs when Q ∝ R³where Q denotes the volumetric flow. This relation, called "cube law", is determined assuming that two energy terms contribute to the cost of maintaining blood flow in any section of any vessel: (i) the pumping power and (ii) the energy metabolically required to maintain the volume of blood which is referred to as "volume constraint". A generalization of this relation can be proposed as Q ∝ R^cwhere c is determined from the condition of minimum power by assuming other constraints (for instance surface constraint yields Q ∝ R^2.52). Under the condition c = 3, the shear stress on the vessel walls is uniform and independent of vessel diameter 3. Several studies have been carried out to determine the value of c 45678 which usually ranges between 2 and 3. The influence of the value of c from 2 to 4 has also been investigated 9. The in vivo wall shear stress in an arterial system has been measured 10. It was found that mean wall shear stress was far from constant along the arterial tree, which implied that Murray's cube law on flow diameter relations could not be applied to the whole arterial system. According to the authors, c likely varies along the arterial system, probably from 2 in large arteries near the heart to 3 in arterioles. A method allowing for estimation of wall shear rate in arteries using the flow waveforms has been developed 11. This work allowed to determine the time-dependent wall shear rates occurring in fully developed pulsatile flow using Womersley's theory. They found a non-uniform distribution of wall shear rates throughout the arterial system.

Following the cubic law, Murray 12 proposed the optimal branching angle. Optimally, the larger branch makes a smaller branching angle than the smaller branch. This work was extended to non-symmetrical bifurcations 13. The arterial bifurcations in the cardiovascular system of a rat have been investigated 14. The results were found to be consistent with those previously reported in humans and monkeys. Murray's optimization problem has also been reproduced computationally using a three dimensional vessel geometry and a time-dependent solution of the Navier-Stokes equations 15.

From Murray's law, some relationships have been proposed between the vessel radius and the volumetric flow, the average linear velocity flow, the velocity profile, the vessel-wall shear stress, the Reynolds number and the pressure gradient 9. In the same way, based on the Poiseuille assumptions, scaling relationships have been described between vascular length and volume of coronary arterial tree, diameter and length of coronary vessel branches and lumen diameter and blood flow rate in each vessel branch 1617.

It is also possible to determine Murray's law using other approaches. A model have been suggested based on a "delivering" artery system of an organ characterized, (i) by the space-filling fractal embedding into the tissue and (ii) by the uniform distribution of the blood pressure drop over the artery system 18. The minimalist principles were not used but the result remains the same. Murray's energy cost minimization have been extended to the pulsatile arterial system, by analysing a model of pulsatile flow in an elastic tube 19. It is found that for medium and small arteries with pulsatile flow, Murray's energy minimization leads to Murray's Law.

Surprisingly, so far, none of the existing methods on Murray's law deal with the non-Newtonian behavior of blood flow although, the non-Newtonian approach for blood flow modeling looks more accurate. Blood is a multi-component mixture with complex rheological characteristics. Experimental investigations showed that blood exhibits non-Newtonian properties such as shear-thinning, viscoelasticity, thixotropy and yield stress 202122. Blood rheology has been shown to be related to its microscopic structures (e.g. aggregation, deformation and alignment of red blood cells). The non-Newtonian steady flow in a carotid bifurcation model have been investigated 2324. The authors showed that in that case, viscoelastic properties may be ignored. The fact that blood exhibits a viscosity that decreases with increasing rate of deformation (shear-thining or so called pseudoplastic behavior) is thus the predominant non-Newtonian effect. There are several inelastic models in the literature to account for the non-Newtonian behavior of blood 2526. The most popular models are the power-law 2728, the Casson 29 and the Carreau 30 fluids. The power-law model is the most frequently used as it provides analytical results for many flow situations. On the usual log-log coordinates, this model results in a linear relation between the viscosity and the shear rate. Blood viscosity have been measured by using a falling-ball viscometer and a cone-plate viscometer for shear rate from 0.1 to 400 s^-131. For both techniques, the authors found that measured values are aligned on a straight line suggesting that the power-law model fits experimental data with sufficient precision.

From the literature review, it can be established that none of the existing studies deal with the minimalist principle along with non-Newtonian models. The combination of both aspects will be studied and are presented hereafter. For a seek of generality, the relationships will be given for an arbitrary number of daughter vessels.

Non-Newtonian model of blood flow

Consider the laminar and isothermal flow of an incompressible inelastic fluid in a straight rigid circular tube of radius R and length l as shown in Fig. (1). For steady fully-developed flow, we make the following hypotheses on the velocity components:

Figure 1

Definition sketch

Definition sketch.

where we have used standard cylindrical coordinates such that the z-axis is aligned with the pipe centerline. It means that the only nonzero velocity component is the axial component which is a function of the distance to the pipe centerline only.

The equation of motion may be written as

where p denotes the pressure and σ_rzis the only nonzero deviatoric stress tensor component. It follows that the pressure is independent from both r and θ. Moreover ∂p/∂z is a constant which we will denote as -Δp/l identified as the pressure drop over the length l. Integrating the third component of the equation of motion, we have:

where K is a constant. As the stress remains finite at r = 0, the constant K must be set to 0. We thus obtain:

where σ_wis the shear stress evaluated at the wall. For a purely viscous fluid, the shear stress σ_rzreads:

where is the generalized viscosity and is the effective deformation rate which is given here by .

Let us consider the power-law constitutive equation proposed by Ostwald 27 and De Waele 28 given by:

It features two parameters: dimensionless flow index n and consistency m with units Pa.sⁿ. On log-log coordinates, this model results in a linear relation between viscosity and shear rate. The fluid is shear-thinning like blood (i.e. viscosity decreases as shear rate increases) if n<1 and shear-thickening (i.e. viscosity increases as shear rate increases) if n > 1. When n = 1 the Newtonian fluid is recovered and in that case parameter m represents the constant viscosity of the fluid. This model is very popular in engineering work because a wide variety of flow problems have been solved analytically for it.

Combining Eqs. (2), (3) and (4) and using the condition of no slip at the wall, we obtain the following velocity field:

It can be noted that the Poiseuille parabolic velocity profile is recovered for n = 1. For a shear-thinning fluid, the velocity profile becomes blunter as n decreases. The flow rate is:

The pressure drop for the flow can be evaluated from Eqs. (2) and (6) to be:

The previous relation can be put in the form:

which reduces in the Newtonian case to the classical Hagen-Poiseuille relation .

Extension of Murray's Law

Let us express Eq. (7) for a vessel k included in a tree structure, it comes:

In a more general manner (also suggested in 2), Eq. (9) may be written as:

where Ψ_κis function of the length l_kof the vessel k and the properties of blood, Q_κis the mass flow rate of blood and R_κis the radius of the vessel k. The parameters a and b are only function of the fluid properties. From Eq. (8), these parameters can be identified for a power law fluid as:

Introduce a cost function Φ_k, as a linear combination of two quantities: the pumping power Q_k·Δp_kand the energy cost to maintain the blood volume π·l_k·R²_k, it yields:

where A_kis a cost factor for pumping and B_kis a sort of maintenance cost of the blood volume. In other words, the metabolic rate of energy required to maintain the volume of blood. The cost function can be written as:

where B_k' = B_k·π·l_kand A_k' = A_k·Ψ_k.

The derivative of this expression with respect to the radius at constant mass flow rate and channel length gives the value of an extremum:

The second derivative of the cost function is

which is found to be always positive because b>0, and so are A_k' and B_k'. As a result, the extremum is a minimum.

According to Eq. (13), the relation between the radius R_kand the mass flow rate Q_kis as follows:

It is noteworthy that if the constraint is not the energy cost to maintain the blood volume but that of the internal area of the vessel (2 π·l_k·R_k), we get:

This expression may be useful when one wants to include mass and/or heat transfer through the vessel wall. Actually, when the vessel diameter decreases, blood catch up with a non-Newtonian fluid and the heat and mass transfer through the vessel wall becomes more and more significant.

Now consider a parent vessel (0) divided into a finite number of vessels (1,..., j) as shown in Fig. (2). Conservation of mass yields the following relation:

Figure 2

Schematic of a bifurcation: parent vessel 0 divided into j daughter vessels

Schematic of a bifurcation: parent vessel 0 divided into j daughter vessels.

Combining Eqs. (17) and (15), we get:

Equation (18) is the generalization of the cube law. In case of laminar Newtonian flow (a = 1 and b = 4, thus c = 3), the classical cube law is recovered. Note that whatever the value of n in the Poiseuille case, c is equal to 3.

The ratio of two consecutive vessels (parent and daughter) is thus written as:

Expression (19) is the generalization of Murray's law and was also proposed by 32. For instance, if one assumes a Poiseuille flow and two daughter vessels one gets the well known result:

Relation (19) takes into account a general form of the pressure drop (not necessarily a Poiseuille flow), a finite but not fixed number of daughter vessels (but not necessarily equal to two) and a possible unequal distribution of the flow in each daughter vessel. Two interesting parameters may thus be calculated:

- The bifurcation index α_irepresents the relative caliber of the symmetry of the bifurcation:

- The area ratio (expansion parameter) β which is the ratio of the combined cross-sectional area of the daughters over that of the parent vessel. Values of β greater than unity produce expansion in the total cross-sectional area available to flow as it progresses from one of the tree to the next. It can be written as:

The following relations are then deduced:

Until now, all equations are formulated with a parent vessel that divid into a finite number of daughter vessels (1,..., j). However, since a parent vessel divide into two daughter vessels in animals and humans, further equations will be formulated with j = 2.

Meaning of the results

In this section, we will examine the meaning of the results obtained in the general case. Particularly, we will focus on the variation of each parameter with the radius of the vessel.

Volumetric flow

As established above, see Eq. (15), the volumetric flow rate is proportional to the radius to the power c when minimizing the cost function:

Velocity of flow

The volumetric flow is proportional to R^cand the cross-area of a vessel is proportional to R², thus the flow velocity (v) is expressed as:

Velocity profile

The maximum velocity, denoted by v_max, is attained at the center of the vessel. It can thus be obtained by setting r to 0 in Eq. (5).

The mean velocity <v> is defined as the ratio between the flow rate (Eq. (6)) and the vessel cross section area:

From these expressions, the velocity profile is independent of the radius of the vessel and is given by the following relation:

As a result, the velocity profile is only function of the fluid properties and remains the same whatever the radius.

Vessel wall shear stress

The vessel wall shear stress may be expressed as:

In the particular case of c = 3, the vessel wall shear stress remains unchanged all along the vascular system.

If c<3, when blood flows from the parent to the daughter vessels, the vessel wall shear stress increases because the vessel radii decrease in the flow direction. On the contrary, if c>3, the vessel wall shear stress decreases because the vessel radii increase in the flow direction.

Reynolds number

The Reynolds number (Re) is proportional to the radius R multiplied by the flow velocity v. According to the relation obtained above for the velocity of flow, it comes:

Since c is often greater than two, the Reynolds number will always increases in the direction of the blood path.

Pressure gradient

The pressure gradient is proportional to Q^a/R^b, i.e. R^ca/R^bsince Q ∝ R^c. The relation between the pressure drop and the vessel radius is then:

Conductance and resistance

In a Murray system, ΣR^cis constant. In addition, the resistance of the fluid is proportional to R^-b. We thus obtain:

and the reciprocal of resistance is the conductance defined as:

Cross sectional area

ΣR^cis constant in a Murray system and the cross sectional area of vessels is proportional to ΣR^-2. As a consequence, the cross sectional area is:

Entropy generation

Entropy generation has several origins: heat transfer, mass transfer, pressure drop... Entropy generation depends on the internal physical phenomena encountered in a process. In the case of an isothermal flow, entropy generation exists, can be quantified, and is related to mass transfer and pressure drop. In our case, entropy generation (S') may be obtained by dividing the cost function (based on the surface constraint) by temperature, which is assumed uniform and constant here (T = 310.15K) 2. The volume constraint cannot be used in that case because the cost of blood maintenance is a process which is external to the system. The minimum entropy generation is thus reached at the minimum of the cost function. In case of surface constraint, the expression of the entropy generation is defined as:

where A_k" = A_k·Ψ_k/T and B_k" = B_k·2·π·l_k/T. Combining with Eq. (10), Expression (25) reduces to

In case of two daughter vessels, the link between the entropy generation upward and downward the arterial bifurcation is defined as:

For the Newtonian case, c = 2.5 and = 1.52. For the non-Newtonian case, n = 0.74 for instance and = 1.50. Whatever the fluid, as >1, the entropy generated in all the daughter vessels is greater than the entropy generated in the parent vessel. Furthermore, this result means that the difference of entropy generation between the parent and daughter vessels is smaller for a non-Newtonian fluid than for a Newtonian fluid. This behaviour can be related to the velocity profile, which is blunter for a non-Newtonian fluid, as shown by Eq. (5).

Illustrating example

Egushi and Karino 31 measured blood viscosity as a function of shear rate using the classical cone-plate viscometer and obtained

Thanks to a falling-ball viscometer, they obtained:

Table 1 shows the effect of n on the c parameter for both constraints. When the cost function involves the volume constraint, parameter c equals three whether the fluid is described by the Newtonian law or by a power-law model. In that case, Murray's law remains valid for a shear-thinning fluid like blood. The optimal ratio between a parent vessel radius and daughter radii is the same whether the fluid is Newtonian or not. In contast, parameter c depends on the value of n when the cost function involves the surface constraint. The optimal ratio of parent vessel radius to daughter radii thus depends on the fluid properties. It is found that decreasing n leads to a drop of parameter c. In conclusion, the more shear-thinning the fluid is, the lower the optimal ratio. c thus ranges from 2.42 to 3, depending on the constraint, which corresponds to the typical range of c measured experimentally and reported in the literature 45678.

Table 1

Influence of n on the c parameter for the two different constraints.

Nominal value of n

Volume constraint

Surface constraint

2.5

0.81

0.78

2.44

0.81

2.45

0.84

2.46

0.74

0.72

2.42

0.74

2.43

0.76

2.43

Table 2 summarizes the influence of index n on different parameters. For example, when assuming a volume constraint, the Reynolds number varies as the radius to power two. It means that whatever the fluid model, the Reynolds number follows the same law. For a surface constraint, decreasing n leads to a decrease in the exponent of the radius. In other words, a shear-thinning fluid the Reynolds number varies less with the vessel radius, all other parameters being kept constant. In addition, the flow resistance is proportional to R to the power c' where c' is always negative. This result agrees with the fact that the greatest part of the resistance of the arterial tree is in the smallest vessels.

Table 2

Influence of n on different parameters

Parameters

Newtonian

n = 1

Non-Newtonian

n = 0.81

Non-Newtonian

n = 0.74

ΣR³

ΣR^2.5

ΣR³

ΣR^2.45

ΣR ³

ΣR ^2.43

Volumetric flow

R ³

R ^2.5

R ³

R ^2.45

R ³

R ^2.43

Velocity of flow

R ^0.5

R ^0.45

R ^0.43

Vessel wall shear stress

R ^-0.5

R ^-0.55

R ^-0.57

Reynolds number

R ²

R ^1.5

R ²

R ^1.45

R ²

R ^1.43

Pressure gradient

R ^-1

R ^-1.5

R ^-1

R ^-1.45

R ^-1

R ^-1.42

Conductance

R ^1.5

R ^0.43

R ^0.98

R ^0.22

R ^0.79

Resistance

R ^-1

R ^-1.5

R ^-0.43

R ^-0.98

R ^-0.22

R ^-0.79

Cross sectional area

R ^-1

R ^-0.5

R ^-1

R ^-0.45

R ^-1

R ^-0.43

Entropy generation

R ²

βand

1.26

1.52

1.26

1.51

1.26

1.50

Conclusion

Blood is a multi-component mixture with complex rheological characteristics. Experimental investigations have shown that blood exhibits non-Newtonian properties such as shear-thinning, viscoelasticity, thixotropy and yield stress. Blood rheology is shown to be related to its microscopic structures (e.g. aggregation, deformation and alignment of blood cells and plattelets). Shear-thinning is the predominant non-Newtonian effect in bifurcations of blood flows.

In this study, we have proposed for the first time an analytical expression of Murray's law using a non-Newtonian blood flow model (power law model), assuming two different constraints in addition to the pumping power: (i) the volume constraint and (ii) the surface constraint. Surface constraint may be useful if one wants to include heat and/or mass transfer in the cost function, specially in capillaries. For a seek of generality, the relationships have been given for an arbitrary number of daughter vessels. Note that there is an alternative formulation of the constrained optimization problem using the Lagrange multipliers, as discussed in 32. However, using this approach, the results presented in this paper would not have been modified.

It has been showed that for a cost function including the volume constraint, classical Murray's law remains valid (i.e. ΣR^c= cste with c = 3 is verified). In other words, the value of c is independent of the fluid properties. On the contrary, for a cost function including the surface constraint, different values of c may be calculated depending on the fluid properties, i.e. the value of n. The fluid is shear-thinning if n<1 and shear-thickening if n>1. When n = 1 the Newtonian fluid is recovered. In the present study, we have used two different blood values of n found in the literature, namely n = 0.81 and n = 0.74. In summary, it has been found that c varies from 2.42 to 3 depending on the constraint and the index n. For the particular Newtonian model, the surface constraint leads to c = 2.5.

Entropy generation has several origins: heat transfer, mass transfer, pressure drop, etc. The cost function (based on the surface constraint) can be related to entropy generation by dividing it by the temperature. It has been demonstrated that the entropy generated in all the daughter vessels is greater than the entropy generated in the parent vessel. Furthermore, it is shown that the difference of entropy generation between the parent and daughter vessels is smaller for a non-newtonian fluid than for a Newtonian fluid. This behaviour can be related to the velocity profile, which is blunter for a non-Newtonian fluid, as shown by Eq. (5).

Based on the literature review and on our work, we propose in the following, further possible investigations:

- The effect of singularities on the cost function has hardly ever been investigated. Few works exist on this aspect 33 and Tondeur et al. (Tondeur D, Fan Y, Luo L: Constructal optimization of arborescent structures with flow singularities. Chem. Eng. Sci. 2009, submitted.). However, the effect of the singularities on the entropy generation might not be negligible. This contribution should be added in the cost function in the future.

- In reality, heat and mass might be transfered through the vessel wall, leading to resistance that should be included in the cost function. Indeed, gases, nutrients and metabolic waste products are exchanged between blood and the underlying tissue. Substances pass through the vessels by active or passive transfer, i.e. diffusion, filtration or osmosis. Moreover, pathologic states such as edema and inflammation might increase such phenomena.

- It is accepted that pulsatile blood flow is more realistic than steady-state flow. The cost function should also include the effect of pulsatile flow in an elastic tube.

- Blood is essentially a two-phase fluid consisting of formed cellular elements suspended in a liquid medium, the plasma. The corpuscular nature of blood raises the question of whether it can be treated as a continuum, and the peculiar makeup of plasma makes it seem different from more common fluids. In particular, when the vessel radius decreases down to the smallest capillaries, the continuum approach diverges from the reality. Treating blood as a non-continuum fluid should be a possible next step.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

RR developed the general equations of Murray's law and carried out the entropy generation analysis.

FR developed the equations for the non-Newtonian model of blood flow.

DB participated in the design of the study and the manuscript.

JB participated in the design of the study and the manuscript.

All the authors read and approved the final manuscript.

The physiological principle of minimum work. I. The vascular system and the cost of blood volume Murray CD Proc Natl Acad Sci U S A. 1926 12 3 207 214 1084489 16576980 Design and scaling laws of ramified fluid distributors by the constructal approach Tondeur D Luo L Chem Eng Sci 2004 59 1799 1813 Branching exponent heterogeneity and wall shear stress distribution in vascular trees Karau KL Krenz GS Dawson CA Am J Physiol Heart Circ Physiol. 2001 280 3 H1256 H1263 11179071 Structure-function relationships in the pulmonary arterial tree Dawson CA Krenz GS Karau KL Haworth ST Hanger CC Linehan JH J Appl Physiol 1999 86 569 583 9931193 Basal EDRF activity helps to keep the geometrical configuration of arterial bifurcations close to the Murray optimum Griffith TM Edwards DH J Theor Biol 1990 146 545 573 2273900 Diameters and cross-sectional areas of branches in the human pulmonary arterial tree Horsfield K Woldenberg MJ Anat Rec 1989 223 245 251 2923275 Microvascular blood flow: evidence indicating a cubic dependence on arteriolar diameter Mayrovitz HN Roy J Am J Physiol. 1983 245 6 H1031 H1038 6660303 Estimation of intravascular blood pressure gradient by mathematical analysis of arterial casts Suwa N Niwa T Fukasuwa H Sasaki Y Tohoku J Exp Med. 1963 79 168 198 13979440 On connecting large vessels to small: the meaning of Murray's law Sherman TF J Gen Physiol 1981 78 431 453 2228620 7288393 Wall shear stress as measured in vivo: consequences for the design of the arterial system Reneman RS Hoeks APG Medical and Biological Engineering and Computing 2008 46 499 507 Distribution of wall shear rate throughout the arterial tree: A case study Stroev PV Hoskins PR Easson WJ Atherosclerosis 2007 191 276 280 16828101 The physiological principle of minimum work applied to the angle of branching of arteries Murray CD J Gen Physiol 1926 9 835 841 Nonsymmetrical bifurcations in arterial branching Zamir M J Gen Physiol 1978 72 837 845 2228494 731200 Arterial Bifurcations in the Cardiovascular System of a Rat Zamir M Wrigley SM Langille BL J Gen Physiol 1983 81 325 335 2215579 6842176 A computational framework for derivative-free optimization of cardiovascular geometries Marsden AL Feinstein JA Taylor CA Comput Methods Appl Mech Eng 2008 197 1890 1905 Scaling laws of vascular trees: of form and function Kassab GS Am J Physiol Heart Circ Physiol. 2006 290 2 H894 H903 16143652 Design of coronary circulation: A minimum energy hypothesis Kassab GS Comput Methods Appl Mech Eng 2007 196 3033 3042 On the Principles of the Vascular Network Branching Gafiychuk VV Lubashevsky IA J Theor Biol 2001 212 1 9 11527441 Pulsatile blood flow, shear force, energy dissipation and Murray's Law Painter PR Edén P Bengtsson HU Theor Biol Med Model 2006 3 31 1590016 16923189 Sheardependent deformation of erythrocytes in rheology of human blood Chien S Usami S Dellenback RJ Gregersen MI Am J Physiol. 1970 219 1 136 142 5424839 Frequency and shear rate dependence of viscoelasticity of human blood Thurston GB Biorheology 1973 10 375 381 4772010 Rheological parameters fort he viscosity, viscoelasticity and thixotropy of blood Thurston GB Biorheology 1979 16 3 149 162 508925 The influence of the non-Newtonian properties of blood on the flow in large arteries: steady flow in a carotid bifurcation model Gijsen FJH van de Vosse FN Janssen JD J Biomech. 1999 32 6 601 608 10332624 The influence of the non-Newtonian properties of blood on the flow in large arteries: unsteady flow in a 90 degrees curved tube Gijsen FJH Allanic E van de Vosse FN Janssen JD J Biomech. 1999 32 7 705 713 10400358 Evaluation of Continuum Models for Characterizing the Constitutive Behavior of Blood Yeleswarapu KK PhD Thesis University of Pittsburgh, Pittsburgh, PA 1996 An Overview of Some Mathematical Models of Blood Rheology Sequeira A Janela J A Portrait of State-of-the-Art Research at the Technical University of Lisbon Springer Netherlands Seabra Pereira M 2007 65 87 About the rate function of the viscosity of dispersed systems Ostwald W Kolloid Z 1925 36 99 117 Viscometry and plastometry De Waele A Oil Color Chem Assoc J 1923 6 33 88 A flow equation for pigment-oil suspensions of the printing ink type Casson N Rheology of Disperse Systems Pergamon Press, New York, NY Mill CC 1959 An Analysis of the Viscous Behavior of Polymeric Solutions Carreau PJ De Kee D Daroux M Can J Chem Eng 1979 57 135 141 Measurement of rheologic property of blood by a falling-ball blood viscometer Egushi Y Karino T Ann Biomed Eng. 2008 36 4 545 553 18259867 Optimal distribution of viscous dissipation in a multi-scale branched fluid distributor Luo L Tondeur D Int J Thermal Sci 2005 44 1331 1141 Tree-shaped flow structures with local junction losses Wechsatol W Lorente S Bejan A Int J Heat Mass Transfer 2006 49 2957 2964