A Critical Reassessment of the Hess–murray Law

The Hess–Murray law is a correlation between the radii of successive branchings in bi/trifurcated vessels in biological tissues. First proposed by the Swiss physiologist and Nobel laureate Walter Rudolf Hess in his 1914 doctoral thesis and published in 1917, the law was " rediscovered " by the American physiologist Cecil Dunmore Murray in 1926. The law is based on the assumption that blood or lymph circulation in living organisms is governed by a " work minimization " principle that—under a certain set of specified conditions—leads to an " optimal branching ratio " of r i`1 r i " 1 3 ? 2 " 0.7937. This " cubic root of 2 " correlation underwent extensive theoretical and experimental reassessment in the second half of the 20th century, and the results indicate that—under a well-defined series of conditions—the law is sufficiently accurate for the smallest vessels (r of the order of fractions of millimeter) but fails for the larger ones; moreover, it cannot be successfully extended to turbulent flows. Recent comparisons with numerical investigations of branched flows led to similar conclusions. More recently, the Hess–Murray law came back into the limelight when it was taken as a founding paradigm of the Constructal Law, a theory that employs physical intuition and mathematical reasoning to derive " optimal paths " for the transport of matter and energy between a source and a sink, regardless of the mode of transportation (continuous, like in convection and conduction, or discrete, like in the transportation of goods and people). This paper examines the foundation of the law and argues that both for natural flows and for engineering designs, a minimization of the irreversibility under physically sound boundary conditions leads to somewhat different results. It is also shown that, in the light of an exergy-based resource analysis, an amended version of the Hess–Murray law may still hold an important position in engineering and biological sciences.


Introduction
In this paper, the history of the Hess-Murray law, a correlation that relates the diameters of successive bifurcations in blood flow, is briefly revisited, its major consequences analyzed and its applicability to real (natural or artificial) artefacts examined.It is then argued that an exergy-based reinterpretation of the law may lead to more rigorously based and more general applications.

The Original Formulation by Hess
The Hess-Murray's law is a formula for relating the radii of daughter branches to the radii of the parent branch of a bi-or multi-furcated circulatory or respiratory system.The original formulation was published by Hess (1949 Nobel laureate in Physiology) [1,2].Hess stated his law in its final form his 1914 thesis: as a historical remark, it must be mentioned though that in his 1808 Croonian Lecture [3], Thomas Young, though, in a rather obscure wording and without providing any theoretical justification, had already stated a "cubic root of 2" law.Hess' intent was to determine the vessel radius that required minimum expenditure of energy by the organism.This approach, novel for his times, led him to adopt an "energy minimization" procedure to derive the optimal ratio between the radii of two successive branchings in vascular tissue.In Hess' words, (The optimal ratio) between the main and the secondary branch is reached when the tendency to control the blood flow is such that the regulatory action leads to the smallest possible reduction of the cross-section.This ratio is attained when the main-and secondary branch have the same effect on the resistance to the flow.
-Walter Rudolf Hess [2] (p.473) Hess provided an explicit proof of the "2 ´1 3 " 0.7937 ratio", which he obtained in three quite simple steps: (a) In a symmetric branching (Figure 1) the mass flowrate is equally split between the secondary branches: m i`1 " 0.5m i , and, for an incompressible fluid, Q i`1 " 0.5Q i , i.e., the two volumetric flowrates in the bifurcated sections are each equal to half the total volumetric flowrate in the main branch.Hess first relates the pressure drop ∆p in both the main and the secondary branches to the respective flow rates through a "viscous resistance" R, and then calculates a "total work expenditure" as the sum of the metabolic work related to the blood flow and the pumping work needed to overcome the frictional pressure drop; (b) The resistance to the flow in a circular tube of length L and cross-sectional area A is given by Poiseuille's formula, which Hess rearranges to obtain R " ∆p Q " 8πµL A 2 kg/(s¨m 4 ).He argues that, for a fixed flowrate and for the laminar isothermal flow of an incompressible newtonian fluid with constant properties, the "optimal control" is established when the variations of the hydraulic resistance induced by an infinitesimal variation in the flow area attain the same value in the main and in the daughter branches, i.e., when: ¯2, it follows that: It is noteworthy that Hess was the first scientist to make explicit mention of the teleology of blood transport: he wrote " . . .(the bifurcation) splits the blood flow delivered from the central system in such a way that the supply requested by the internal organs is ensured with the minimum possible flow volume of blood" [2].
Entropy 2016, 18, 283 2 of 17 the vessel radius that required minimum expenditure of energy by the organism.This approach, novel for his times, led him to adopt an "energy minimization" procedure to derive the optimal ratio between the radii of two successive branchings in vascular tissue.In Hess' words, (The optimal ratio) between the main and the secondary branch is reached when the tendency to control the blood flow is such that the regulatory action leads to the smallest possible reduction of the cross-section.This ratio is attained when the main-and secondary branch have the same effect on the resistance to the flow.
-Walter Rudolf Hess [2] (p.473) Hess provided an explicit proof of the "2 = 0.7937 ratio", which he obtained in three quite simple steps: (a) In a symmetric branching (Figure 1) the mass flowrate is equally split between the secondary branches: = 0.5 , and, for an incompressible fluid, = 0.5 , i.e., the two volumetric flowrates in the bifurcated sections are each equal to half the total volumetric flowrate in the main branch.Hess first relates the pressure drop ∆p in both the main and the secondary branches to the respective flow rates through a "viscous resistance" R, and then calculates a "total work expenditure" as the sum of the metabolic work related to the blood flow and the pumping work needed to overcome the frictional pressure drop; (b) The resistance to the flow in a circular tube of length L and cross-sectional area A is given by Poiseuille's formula, which Hess rearranges to obtain = ∆ = kg/(s•m 4 ).He argues that, for a fixed flowrate and for the laminar isothermal flow of an incompressible newtonian fluid with constant properties, the "optimal control" is established when the variations of the hydraulic resistance induced by an infinitesimal variation in the flow area attain the same value in the main and in the daughter branches, i.e., when: (c) Since = , it follows that: It is noteworthy that Hess was the first scientist to make explicit mention of the teleology of blood transport: he wrote "…(the bifurcation) splits the blood flow delivered from the central system in such a way that the supply requested by the internal organs is ensured with the minimum possible flow volume of blood" [2].

Murray's Rediscovery of the Hess' Law
Murray's analysis [4] was originally aimed at the identification of the most efficient oxygen transport in man.Since oxygen is supplied to the body organs by the blood through the vascular

Murray's Rediscovery of the Hess' Law
Murray's analysis [4] was originally aimed at the identification of the most efficient oxygen transport in man.Since oxygen is supplied to the body organs by the blood through the vascular system, he was faced with the problem of, as he put it, "the cost of blood volume".In line with the logical empiricism of his time, he explicitly remarks that the most suitable model for the analysis of metabolic processes is one based on an energy (power) balance, "all other factors remaining constant".Following Hess' ideas, Murray states "If the vessels are too small, the work required to drive the blood through them becomes too great; if the volume of the vessels is too large, the volume of blood, being equally large, becomes a burden to the whole body" [4] (p. 210).By a procedure slightly different from that proposed by Hess, he re-derived the same "cubic root of 2" rule: (a) Murray [4] starts by expressing the total power involved in the circulation of blood as the sum of two terms: the first needed to overcome the frictional resistance through the vessel and the second corresponding to the metabolic effort needed to "create and maintain" the blood volume therein: (b) In line with Hess, the pressure drop is expressed by the Hagen-Poiseuille's law: (c) Murray then uses Equation (3) to compute the metabolic cost c b , and then equates the total resistances R 0 and R 1+2 w.r.t. the flow area to derive his "cubic root of 2" law.Notice that the appearance of the cubic root does not depend on the value of the "specific metabolic cost of blood" c b , which is assumed constant; (d) By differentiating Equation (4) w.r.t.r and equating the result to 0, Murray argues that, for given c b and L, Q is proportional to r 3 ; (e) If the above point is correct, and since the flow obeys Poiseuille's law, the pressure gradient in successive branchings is proportional to r ´1: (f) Under the same stipulation, the flow velocity correspondingly decreases: (g) Correspondingly, the Reynolds number Re scales with the radius ratio as: and also decreases along successive branchings.
Murray was aware of Hess' work, but considered it preliminary and incomplete, and made only one reference in passing to it [4].Although his derivation does not make use of an "optimal control" assumption (he reasons purely in terms of metabolic energy needed to maintain the blood flow and energy dissipation by friction generated by the bifurcation), his equations indeed derive, contrary to what he wrote, from Hess' formulae, because both derivations are equivalent to the assumption of the constancy of the wall stress along successive branches (see Section 1.5 below).In spite of this slip, Murray's papers contain genuine novelty and are quite propositive: he uses the "cubic root of 2" paradigm to draw some inference about the correspondence between the results of his law and the experimental oxygen rates in capillary and pulmonary vessels available at that time.He states that, while the experimental results are in satisfactory agreement with the law for capillary flow, for larger vessels the agreement is weak, and attributes this discrepancy to "side effects" like turbulence and vessel wall structure.
In [5] he also provides an original formula for the branching angle, given the splitting point and the two endpoints of the branchings (Figure 2), extends the treatment to non-symmetric branchings and to trifurcations, and remarks that the results ought to apply equally well (and with the same limitations) to the branchings of xylems, the sap transport system in plants.
An interesting finding in Murray's analysis is the correlation between the branching angle and the radius ratio, which he derived using a constrained optimization of the pumping work: while referring readers to [5] for the details, the formulae he derives are: which in the case of symmetric splitting become: which, upon substitution of δ 1 " 1 3 ?

2
, provides α 1 = α 2 = 37.5 ˝.Note that, though Murray's original drawing shows different lengths for the branches, nowhere in his derivation is the length L related to the vessel diameter: he posits no "slenderness ratio" for the vessels.
Entropy 2016, 18, 283 4 of 17 and to trifurcations, and remarks that the results ought to apply equally well (and with the same limitations) to the branchings of xylems, the sap transport system in plants.
An interesting finding in Murray's analysis is the correlation between the branching angle and the radius ratio, which he derived using a constrained optimization of the pumping work: while referring readers to [5] for the details, the formulae he derives are: which in the case of symmetric splitting become: which, upon substitution of = √ , provides α1 = α2 = 37.5°.Note that, though Murray's original drawing shows different lengths for the branches, nowhere in his derivation is the length L related to the vessel diameter: he posits no "slenderness ratio" for the vessels.

A brief Summary of Experimental Results and of Other Allometric Laws
As remarked by Sherman [6] in a very well-documented review, the Hess-Murray law (H-M in the following) is only approximately valid for natural branched systems (animal tissues and tree structures): the exponent in the root is found experimentally to range between 2 and 3 (though some experimental findings seem to substantially widen this range, [7]) and the introduction of a correction for turbulent effects leads to a 7/3 = 2.3 value [8].Confirming the previous findings by Murray, the predictions of the H-M law are closer to reality for the smaller vessels (arterioles and capillaries).Additional deductions from-or comments & modifications to-the law were made by different authors: (a) Assuming thin-walled vessels, the external radius does not play any significant role, because it is-at the considered scales-negligibly larger than the internal radius [6]; (b) For 3-D branchings, the H-M still applies, because over 90% of the bifurcations are experimentally found to be coplanar [9]; (c) A different "optimal branching ratio" can be derived by using fractals [10,11]: however, even for the smallest pulmonary alveoli (for which values 0.4 < δ < 0.85 are found experimentally) its

A brief Summary of Experimental Results and of Other Allometric Laws
As remarked by Sherman [6] in a very well-documented review, the Hess-Murray law (H-M in the following) is only approximately valid for natural branched systems (animal tissues and tree structures): the exponent in the root is found experimentally to range between 2 and 3 (though some experimental findings seem to substantially widen this range, [7]) and the introduction of a correction for turbulent effects leads to a 7/3 = 2.3 value [8].Confirming the previous findings by Murray, the predictions of the H-M law are closer to reality for the smaller vessels (arterioles and capillaries).Additional deductions from-or comments & modifications to-the law were made by different authors: (a) Assuming thin-walled vessels, the external radius does not play any significant role, because it is-at the considered scales-negligibly larger than the internal radius [6]; (b) For 3-D branchings, the H-M still applies, because over 90% of the bifurcations are experimentally found to be coplanar [9]; (c) A different "optimal branching ratio" can be derived by using fractals [10,11]: however, even for the smallest pulmonary alveoli (for which values 0.4 < δ < 0.85 are found experimentally) its value is not in agreement with the H-M law; (d) Both experimental and numerical findings show that while the H-M law leads-under the above stated assumptions-to an optimal value of δ = 0.7937 for laminar flow, for turbulent flow in veins and arteries the result becomes 0.74 and for the largest arteries, 0.77.Moreover, contrary to Hess & Murray's result of an inverse scaling of the pressure gradient in successive branchings with δ (∆p/L = δ ´1, see Equation ( 4)), radius-dependent slenderness factors are found experimentally: L/r = 15.75r 1.10 for arterioles with r ě 50 µm, 1.79r 0.47 for arterioles with r ď 50 µm, 14.54r for veins with r ď 2000 µm, [12,13].Notice that the fundamental deviation of Huang's results [12] from the H-M law depend also on the different ordering of "successive branches" adopted therein, whose classification cannot be directly extended to capillary networks.Nevertheless, their findings contradict the H-M predictions in arterioles and veins (see also Table 3).Notice that [14] the H-M law is independent of the slenderness ratio r/L.Very well documented critiques of the H-M law are presented in [7,15]; (e) The H-M law assumes impervious vessel walls: this is bound to cause a discrepancy with real lymph and blood flows, in which a major transport mechanism is through selectively permeable walls [14,[16][17][18]; (f) Nothing in the derivation of the H-M.law makes explicit reference to biological vessels: in principle, the law should hold also for human artefacts (pipes and any other kind of conduit) [6]: if the cross-sectional area is not circular, the law of the area ¯2 still approximately holds if the hydraulic radius r h is substituted for r and a proper correction is introduced [19]; (g) The law has important consequences when the growth of a biological organism is taken into account: if it holds, then the radii of the individual vessels in the vascular system must "grow" according to cubic root of 2 because the volume flow rate of the organism grows as well, and so does r 0 , the radius of the primary stem [6,20].This applies also to stress situations: when a higher rate of blood (or sap) flow is needed, the existence (or the emergence) of a self-regulatory response which suitably adjusts the vessels radii, as originally proposed by Hess, seems reasonable [15]; (h) The H-M law also allows for an "allometric scaling" between different species [20]: if the relative sizes of the respective blood (or sap) vessels scale according to this law, then the physical dimensions depend on the radius of the main branch, which in turn depends on the metabolic rate of the species [21]; (i) If blood is modelled as a non-Newtonian fluid, the 2 ´1/3 law does not apply anymore, because-as expected-the exponent is found to have a non negligible dependence on the viscosity law (from 1 to 2.5 in [22,23]).Miguel [17,18] analyzed the irreversibility in branched flows of a power fluid and showed that (a) the "optimal" radius ratio depends on the "degree of non-newtonianity" of the fluid, ω; (b) that in the limit of Newtonian behavior, ω tends to 1 and the H-M law is recovered; (j) It is also possible to derive the optimal radius ratio in a dendritic structure by reformulating the problem as a constrained multi-variable optimization [11,24,25]: this approach, useful as it may be in engineering applications, loses somewhat of its significance in biological investigations.
On the basis of the above remarks, we can state that although the H-M law represents an innovative and theoretically interesting approach to the scaling in dendritic structures, its experimental confirmation and its generality are not completely satisfactory.We shall return to this point in Section 3.

Bejan's Constructal Law
In 1997 Adrian Bejan [26] presented a reinterpretation and a very clever extension of the H-M law: expanding one of Murray's statements (see comment to Figure 2 above) he reasoned that, if the law holds, then every bifurcated stream that carries some material (sap, blood, particles) or immaterial (energy under any form) flow from one point to two (or more) others in a given domain can-in an evolutionary sense-only develop along a global structure exactly specified by the "cubic root of 2" rule.This intuition forms the basis for his Constructal Theory (CT in the following), which has since seen a plethora of applications in a wide range of fields (see for instance [20,27]).Though CT is surely not as omni-comprehensive as it is being represented, its greatest merit is to have shown that-in the words of Bejan-"Nature is not fractal": real structures follow from goals (Hess' and Murray's teleology!)rather than from fractal exponents, the goals being determined by external factors related, in nature, to evolutional pressures, and in engineering, to design specifications.Contrary to previous applications of the H-M law, Bejan's findings are in a real sense "engineering" findings, not only qualitative but quantitative: they lead on the one side to a better interpretation and classification of natural (abiotic and biological) structures, and on the other side to measurable and repeatable design procedures for engineered artefacts.For the purpose of this paper, it suffices here to report one of the clearest definitions of the CT "thought procedure": The flow path was constructed in a sequence of steps that started with the smallest building block (elemental area) and continued in time with larger building blocks (assemblies, or constructs).The mode of transport with the highest resistivity (slow flow, diffusions, walking, and high cost) was placed on the smallest scales, completely filling the smallest elements.Modes of transport with successively lower resistivities (fast flow, streams, vehicles and low cost) were placed in the larger constructs, where they were used to connect the area-point or volume-point flows integrated over the constituents.The geometry of each building block was optimised for area-point access.
-Wishsanuruk Wechsatol, et al. [28] When applied to bifurcating flows, this translates into "higher pressure gradients in the smallest vessels", in complete agreement with the H-M law.In the next sections, it will be though argued that: (a) This very same "construction" derives in fact from a single principle, namely that of minimizing irreversibility in the transport of a certain mass flow rate of fluid from one point to a pair of a triplet of "users"; (b) The objective function whose minimum identifies the optimal configuration is the specific irreversible entropy generation rate; (c) That both natural and artificial dendritic flow structures do not admit of a single "optimal radius ratio" or of an "optimal splitting angle", because both of these quantities are determined by the imposed boundary conditions.Therefore, in an evolutionary sense, it is the random variation of the latter that determines the "shape" of the winning configuration.

A fluiddynamic Interpretation
For the laminar flow of a Newtonian, single-phase, incompressible fluid in a circular tube with a rigid non-permeable skin, the Hagen-Poiseuille law (Equation (4) above) relates the pressure drop over a length L of tube to the volumetric flowrate: Under the same stipulations, the wall shear stress is given by: In a symmetric bifurcation as those shown in Figure 1, Q 1 = Q 2 = 0.5Q 0 , and the ratio of the wall stresses in the main and in the daughter branches depends on the radius ratio: The required pumping power is proportional to the wall stress and to the flowrate: Now, if we denote as θ the wall stress ratio τ w,1 /τ w,0 , by δ the radius ratio r 1 /r 0 and by λ the length ratio L 1 /L 0 , the total pumping power required to deliver the flow from the bifurcation inlet (point O in Figure 1) to its two outlets (points A & B) is: For a given bifurcation (L 0 , r 0 , Q 0 assigned, so that τ w,0 is known), one obtains: For the purpose of the present paper, it is useful to examine Equations ( 11) and ( 14) jointly: the H-M law is equivalent to the assumption of a constant τ w along successive branches.But Equation (15) shows that τ w is not the only factor influencing the total power expenditure (what Hess and Murray call "the pumping work"), because there are several possible cases, based on different assumptions, that are relevant to the problem: Other cases may be of importance in physiological flows, like for instance a lower velocity in smaller branches, to facilitate nutrients absorption from the permeable wall [14,16,17], or wall thickness proportional to radius, which may be convenient for sap flow on plants, but they are outside of the goals of the present study and will not be considered here.

Scope of This Work
In conclusion, the scope of this work is to prove that: 1.The Hess-Murray law is only one of several possible instantiations of a class of rules that can be derived from irreversibility minimization; 2. That there are therefore different structures that can develop if one or more of the assumptions are modified; 3.That if an exergy-based resource cost paradigm is employed, the Hess-Murray law is amenable to a more fundamental and interesting interpretation that complements and extends the henceforth published results.
Since our goal is to demonstrate that the "configurator" of branched structures is the specific entropy generation rate, which depends on the shape of the domain, on the fluid properties and on the physical characteristics of the material of the vessel walls and of the surrounding tissues, only symmetric bifurcations will be considered here.Including non-symmetric cases would introduce additional geometric parameters that would complicate the treatment while adding nothing to the phenomenological significance of the model and to the interpretation of the results.For the same reason, the bifurcation losses, i.e., the viscous loss generated at the cusp, was omitted as well.For a more complete discussion on the effects of localized entropy generation, see [29].

Irreversibility in Viscous Isothermal Flows
Consider the flow of a viscous Newtonian fluid through a circular vessel: the irreversibility of the flow can be measured by its specific entropy generation rate (per unit mass).Following [30] we can separate the viscous and thermal effects, and in cylindrical coordinates (r,x) we obtain: where U and V are the axial and radial velocity respectively and tangential periodicity is assumed.U, V and T are local and time-dependent quantities, and therefore in the context of the present analysis it at the moment is irrelevant whether the flow is turbulent or laminar.The derivative in the axial direction has been maintained to include entrance length effects, which will though be neglected for the moment.
Let us consider, in line with Murray's original suggestion, that the scope of the branching is to deliver the fluid from an origin or source (point O in Figure 1) to two endpoints or sinks (A & B).Our goal is to show that, whether the "basin of influence" of the branching is a rectangle or a circular sector, the total entropy generation within the region where the bifurcation appears is an explicit function of, among other parameters, the aspect ratio of the surrounding domain [14,29].The first step is to relate the values of L 0 and L 1 (a similar derivation applies to any L i+1 ) to the splitting angle α 1 and to the aspect ratio of the area "served" by the flow, which may be defined as a = H/L or a = 2α 0 for rectangular and triangular domains respectively.A first relevant finding is that for triangular domains (circular sectors, a geometry relevant to heat exchanger design), Murray's angle law does not apply, because α 0 ď α 1 .
Figure 1 shows the respective formulae, derived using elementary trigonometry.We shall henceforth present the calculations for a single bifurcation in rectangular domains, as a function of the aspect ratio a = H/L that defines the domain in which the bifurcation takes place.A remark is in order here: all calculations presented in this paper are related to steady flow, and pulsating effects and compliant walls are excluded from this model: the emphasis is placed here on rigorous fundamental principles rather than on a higher degree of complexity.
Let us first assume, in accordance with the premises of the Hess-Murray law, that the flow structure "evolves" so as to minimize the irreversibility: this is equivalent to Hess' statement of minimizing the pressure loss per unit length, except that now, by means of Equations ( 16) and (17) we could also reformulate it in terms of total enthalpy.This paper deals though only with isothermal, fully developed laminar flows, for which Equation (17) vanishes identically, and so do the terms B/Bx in Equation (16).It is convenient to introduce dimensionless quantities: u = U/U 0 , v = V/U 0 , r r = r/r 0 , r x = x/L, etc., with r 0 and U 0 respectively the radius of the main stem (the first branch) and the average velocity therein, and L the length of the domain shown in Figure 1.Then, substituting the Poiseuille velocity profile and integrating Equation ( 16) from r = 0 to r = 1, and then again from x = 0 to x = L 0 /L, we obtain for the entropy generation rate for a circular vessel of length L 0 : Notice that in Equation ( 18) entrance length effects are neglected: a detailed investigation [29] shows that, as long as the flow is isothermal, they do not significantly modify the results.Under the first stipulation of the H-M law, let us now consider only symmetric bifurcations in which each daughter branch carries m i+1 = 0.5m i kg/s.Denoting as δ i the ratio r i+1 /r i and as λ i = L i+1 /L i (δ 0 = 1, λ 0 = L 0 /L), the viscous entropy generation rates in successive branches are shown in Table 1: S 0 denotes the entropy generation rate in a tube of length L 0 = L and of radius r 0 (no splitting).
In the line of thought explained in Section 1.5 above, four cases are of interest: 1.The velocity remains the same before and after the bifurcation, U i+1 = U i : thence δ i " 2. The Reynolds number remains the same before and after the bifurcation, Re 1 = Re 0 : δ i " r i`1 r i " 1 2 ; 3. The sum of the volumes of the daughter branches equals the volume of the main stem: 4. The pressure drop in the daughter branches equals that in the main stem and the radius ratio is dictated by the Hess-Murray law (δ i " λ i " 1 3 ?
To better understand the premises of this study, it is important to remark once more that there is no a priori physical reason for either one of the four above choices: their occurrence may be in fact facilitated by evolutionary pressures in nature and imposed by designer's choice in engineered artefacts.For instance, a constant velocity may be preferred if the volume to be "fed" by the fluid varies with the cross section of the vessel (as in a leaf) or when a constant residence time in the vessel is more efficient in an evolutionary sense; a constant Reynolds number leads to a constant heat exchange rate per unit length; a "constant volume" law implies that the daughter branches use roughly the same biological material as the main branch.This remark is in line with the argument presented in [19], in which a branching parameter X is introduced and the radius ratio is defined as δ i " 1 3 ?2X : the authors show that values of X different from 1 (X = 1 recovers the H-M limit) give origin to different flow configurations, some of which may be "optimal" in vivo.
Notice that in Equation (19) the length ratio between successive branches appears as well, and since-except for case 4-Murray's "optimal angle" relations do not apply, we need to use some different physical reasoning (cases 1,2,3 above) to establish a direct link between the radius-and the length ratio.In practice, we shall use the general geometric derivation shown in Figure 1 and select α 1 as the running variable.
The results for the first bifurcation in a rectangular domain are shown in Figure 3: for slender domains (a = 0.2) the constant volume case (case 3 above) is the one that generates the least irreversibility for bifurcation angles lower than about 20 ˝, and thereafter the CT case 4 is the least irreversible one.For a = 0.5 the constant volume case 3 is again the most efficient one up to α 1 = 35 ˝, beyond which the CT case is again less dissipative.For a = 0.75, the constant volume and the CT are in practice equivalent beyond α 1 = 40 ˝, while for lower angles the constant volume configuration is the least dissipative.This simple example demonstrates that in practical (albeit quite idealized) cases the failure of the Hess-Murray law to consider the "teleological" aspects of the structure (i.e., the goal for which a bifurcation should occur in nature or by design) makes its application questionable: a radius ratio equal to 1 3 ? 2 is indeed a feasible option, but other values are equally feasible on physical grounds and, with the possible exception of the constant velocity case (not included in Figure 3), all are more or less compatible with the experimental measurements.radius ratio equal to √ is indeed a feasible option, but other values are equally feasible on physical grounds and, with the possible exception of the constant velocity case (not included in Figure 3), all are more or less compatible with the experimental measurements.Along the same line of reasoning it is possible to calculate, with a larger computational effort but with no substantial methodological change, the entropy generation in the successive n > 1 branchings (Figure 4 and Table 1).

Is the Hess-Murray Energy Balance the Correct One?
Hess and Murray derived their "cubic root of 2" formula on the basis of an energy (power) balance applied to two different situations: (a) Hess assumed that the variation of the hydraulic viscous resistance w.r.t. a small change in the blood flow is the same in the parent and daughter branches.Today we know that blood supply pressure is regulated also in the arterioles, upstream of the capillaries, and that the latter may function also as temporary "blood reservoirs", and therefore Hess' balance cannot be applied to derive the radius ratio of capillary bifurcations; (b) Murray starts by calculating the total power (which he calls cost) of blood circulation as the sum of the power dissipated by friction losses and of a "metabolic cost" of the blood, calculated as the pumping power requested of the heart.This calculation is though inaccurate, because it neglects for instance the work needed to deliver the blood through the veins back to the heart (in fact, Murray [4,31] obtains three different values for his metabolic cost depending on the different assumptions he posits).Thus, his balance between the blood pumping work (proportional to the vessels radius) and the friction work is completely contingent on the rather arbitrary assumptions he posits to calculate the "pumping power" of the heart.
S/S0 S/S0 S/S0 Figure 3. Minimum irreversibility ratio S/S 0 as a function of the aspect ratio H/L and of the splitting angle α 1 for different physical assumptions [29].
Along the same line of reasoning it is possible to calculate, with a larger computational effort but with no substantial methodological change, the entropy generation in the successive n > 1 branchings (Figure 4 and Table 1).radius ratio equal to √ is indeed a feasible option, but other values are equally feasible on physical grounds and, with the possible exception of the constant velocity case (not included in Figure 3), all are more or less compatible with the experimental measurements.Along the same line of reasoning it is possible to calculate, with a larger computational effort but with no substantial methodological change, the entropy generation in the successive n > 1 branchings (Figure 4 and Table 1).

Is the Hess-Murray Energy Balance the Correct One?
Hess and Murray derived their "cubic root of 2" formula on the basis of an energy (power) balance applied to two different situations: (a) Hess assumed that the variation of the hydraulic viscous resistance w.r.t. a small change in the blood flow is the same in the parent and daughter branches.Today we know that blood supply pressure is regulated also in the arterioles, upstream of the capillaries, and that the latter may function also as temporary "blood reservoirs", and therefore Hess' balance cannot be applied to derive the radius ratio of capillary bifurcations; (b) Murray starts by calculating the total power (which he calls cost) of blood circulation as the sum of the power dissipated by friction losses and of a "metabolic cost" of the blood, calculated as the pumping power requested of the heart.This calculation is though inaccurate, because it neglects for instance the work needed to deliver the blood through the veins back to the heart (in fact, Murray [4,31] obtains three different values for his metabolic cost depending on the different assumptions he posits).Thus, his balance between the blood pumping work (proportional to the vessels radius) and the friction work is completely contingent on the rather S/S0 S/S0 S/S0

Is the Hess-Murray Energy Balance the Correct One?
Hess and Murray derived their "cubic root of 2" formula on the basis of an energy (power) balance applied to two different situations: (a) Hess assumed that the variation of the hydraulic viscous resistance w.r.t. a small change in the blood flow is the same in the parent and daughter branches.Today we know that blood supply pressure is regulated also in the arterioles, upstream of the capillaries, and that the latter may function also as temporary "blood reservoirs", and therefore Hess' balance cannot be applied to derive the radius ratio of capillary bifurcations; (b) Murray starts by calculating the total power (which he calls cost) of blood circulation as the sum of the power dissipated by friction losses and of a "metabolic cost" of the blood, calculated as the pumping power requested of the heart.This calculation is though inaccurate, because it neglects for instance the work needed to deliver the blood through the veins back to the heart (in fact, Murray [4,31] obtains three different values for his metabolic cost depending on the different assumptions he posits).Thus, his balance between the blood pumping work (proportional to the vessels radius) and the friction work is completely contingent on the rather arbitrary assumptions he posits to calculate the "pumping power" of the heart.
Neither Hess nor Murray explicitly formulated their law in its simplest form (which is likely to be the sense of Young's statement mentioned in Section 1.1): for laminar flow, and assuming a Poiseuille profile, the wall shear stress is τ w " ´du dr ¯r"R " 4µQ πR 3 , and if one imposes for τ w to remain constant across a symmetric bifurcation, a radius ratio 2 ´1/3 results (in Murray's derivation, the statement that Q is proportional to R 3 -point d in Section 1.2-is in practice equivalent to assuming a constant wall stress; see Section 1.5 above).Their "energy based" approach, though, is completely sound and cannot be discarded: considering the evolutionary pressures that favor selective competition, it is very likely that the circulatory system has evolved in such a way as to increase its "efficiency", however defined.Let us then reconsider the original scenario under which the H-M law was derived: a net of capillaries (Figure 5) carries a nutrients-laden fluid to a portion of a biological structure.Such a net performs two functions: first, by selective permeability on the outer wall of the vessels, it provides the surrounding organic material with the chemicals it needs; and second, possibly in a different region, absorbs from the surrounding tissue organic detritus and other discharges and delivers them to a different portion of the organism.Globally, in the assumption of steady-state conditions, the in-and out-flows from the walls compensate, and thus the mass flowrate through the individual vessels may be considered constant over their respective lengths.Therefore, in a very real sense, each bifurcation "serves" a small region of its surroundings, the one Al-Shammari [16] correctly identifies as "the supply region".Neither Hess nor Murray explicitly formulated their law in its simplest form (which is likely to be the sense of Young's statement mentioned in Section 1.1): for laminar flow, and assuming a Poiseuille profile, the wall shear stress is = = , and if one imposes for τw to remain constant across a symmetric bifurcation, a radius ratio 2 −1/3 results (in Murray's derivation, the statement that Q is proportional to R 3 -point d in Section 1.2-is in practice equivalent to assuming a constant wall stress; see Section 1.5 above).Their "energy based" approach, though, is completely sound and cannot be discarded: considering the evolutionary pressures that favor selective competition, it is very likely that the circulatory system has evolved in such a way as to increase its "efficiency", however defined.Let us then reconsider the original scenario under which the H-M law was derived: a net of capillaries (Figure 5) carries a nutrients-laden fluid to a portion of a biological structure.Such a net performs two functions: first, by selective permeability on the outer wall of the vessels, it provides the surrounding organic material with the chemicals it needs; and second, possibly in a different region, absorbs from the surrounding tissue organic detritus and other discharges and delivers them to a different portion of the organism.Globally, in the assumption of steady-state conditions, the in-and out-flows from the walls compensate, and thus the mass flowrate through the individual vessels may be considered constant over their respective lengths.Therefore, in a very real sense, each bifurcation "serves" a small region of its surroundings, the one Al-Shammari [16] correctly identifies as "the supply region".
(a) (b) If we observe the detailed pictures of some natural occurrences of such nets (Figure 6), the most striking character of each structure is without any doubt its connectivity: all images resemble "networks" rather than "dendrites".Therefore, it seems reasonable to assume a simplified rectangular geometry for the basin of influence of a single bifurcation, and since the permeability of the wall tissues remains constant, the evolution-driven efficiency of natural structures suggests that each branch drains or provides a mass flow rate proportional to the immediate surrounding portion of the domain.If we observe the detailed pictures of some natural occurrences of such nets (Figure 6), the most striking character of each structure is without any doubt its connectivity: all images resemble "networks" rather than "dendrites".Therefore, it seems reasonable to assume a simplified rectangular geometry for the basin of influence of a single bifurcation, and since the permeability of the wall tissues remains constant, the evolution-driven efficiency of natural structures suggests that each branch drains or provides a mass flow rate proportional to the immediate surrounding portion of the domain.With reference to schematic representation of Figure 5a, the process may be thus described as follows (Figure 7): the mass flow rate of the permeate from (or to) the vessel wall area 2πL0r0 is proportional to the domain portion A0, while that of the permeate from (or to) area 4πr1L1 is proportional to A1 (similar considerations are made in the detailed numerical simulation presented in [31]).This simple apportioning is based on the assumption of a uniform "specific absorbance" by the surrounding cells, and leads to a calculation of the radius ratio based on the function of the net rather than on the "energy balance" proposed by Hess and Murray:  (b) Hydrangea leaf [34]; (c) Capillary network in the human retina [35]; (d) Nettle leaves [36].
With reference to schematic representation of Figure 5a, the process may be thus described as follows (Figure 7): the mass flow rate of the permeate from (or to) the vessel wall area 2πL 0 r 0 is proportional to the domain portion A 0 , while that of the permeate from (or to) area 4πr 1 L 1 is proportional to A 1 (similar considerations are made in the detailed numerical simulation presented in [31]).This simple apportioning is based on the assumption of a uniform "specific absorbance" by the surrounding cells, and leads to a calculation of the radius ratio based on the function of the net rather than on the "energy balance" proposed by Hess and Murray: With reference to schematic representation of Figure 5a, the process may be thus described as follows (Figure 7): the mass flow rate of the permeate from (or to) the vessel wall area 2πL0r0 is proportional to the domain portion A0, while that of the permeate from (or to) area 4πr1L1 is proportional to A1 (similar considerations are made in the detailed numerical simulation presented in [31]).This simple apportioning is based on the assumption of a uniform "specific absorbance" by the surrounding cells, and leads to a calculation of the radius ratio based on the function of the net rather than on the "energy balance" proposed by Hess and Murray:  The rationale of the apportioning of the delivery of nutrients is the following: (a) The portion of the supply region denoted as A 0 receives all of (and only) the nutrients permeated through the main branch: . m A0 " 2πr 0 L 0 µ, where µ is the wall permeability (kg/(m 2 s)), assumed to remain constant throughout.By analogy for the remaining portion A 1 , .m A1 " 4πr 1 L 1 µ; (b) The "costs" of the irreversibility, S 0 and S 1 , must be attributed to the mass flow rates m 0 and m 1 respectively; (c) And the objective function becomes Equation ( 20), together with the formulae that link λ 0 and λ 1 to the splitting angle α 1 (half of bifurcation angle, Figure 1), sheds new light on the phenomenology: the structure of the capillary network is dictated by the instantaneous local balance between the drained or absorbed mass flowrate of nutrients and the permeability of the capillary walls (notice that the introduction of this "wall permeability" is in contrast with one of the initial assumptions, namely, that the walls are non-permeable.The relaxation of this constraint is acceptable for the case in study though, where the mass fraction of the permeate is of orders of magnitude smaller than the main flow rate of fluid in the vessel).In a very real sense, this is a "form from function" derivation!And since, in an evolutionary competition, conversion efficiency is a bonus, it is to be expected that structures will prevail that generate the minimum possible irreversibility per unit mass of supplied or absorbed fluid (same product, less resources used).The formula for the total entropy generation rate in the vessels is (second line of Equation ( 19)): where S L " 8µLQ 2 0 πr 4 0 T 0 is the irreversible entropy generation rate, in W/K, in a vessel of radius r 0 bisecting the domain for its entire length L. The total drained/leaked mass flowrate is: And the most efficient configuration, i.e., the one that requires the least pumping power to ensure the survival of the cells in the domain, is that for which the ratio ∆S/m dr is lowest.Table 2 displays the values obtained for the dimensionless ratio σ " ∆S {∆S L ∆m d{∆m d , L .The corresponding geometries are displayed in Figure 8.Some consequences of the above derivation are noteworthy: (1) The bifurcation angle (2α1) varies between 54° and 68° and is a function of the domain aspect ratio.This is not in agreement with Murray's cosine rule for the "optimal angle", that predicts 2α1 = 75°; (2) The radius ratio is a function of the splitting angle only, which is a remarkable result, because it was derived from a wall suction/drainage balance which in fact minimizes irreversible effects.Though much lower than that (0.7937) predicted by the H-M law, the values of δ1 shown in Table 2 are within the range of the available experimental data: a series of indicative measurements on the images displayed in see Figure 6a,b and d show that the radius ratio of the first bifurcation is in the range 0.3 ÷ 0.75; (3) The more slender the domain (lower aspect ratio), the shorter λ1 and the higher δ1.This implies that the pressure gradient in the main and in the daughter branches is also a function of the splitting angle, and (Table 2) it is higher in the latter, as predicted also by the H-M law; (4) The pressure drop is also a function of the aspect ratio and it does not assume equal values in the parent and daughter branches, its value being always higher in the latter: while for lower aspect ratios the total pressure drop is almost equally allocated, for square domains the ∆p1 is 25 times higher than the ∆p0.Incidentally, this also negates the assumption of "equal pumping work" that is at the basis of the H-M law, but agrees with the dictates of Constructal Theory; (5) Considering that entropy generation rate per unit mass of transported fluid is proportional to the "lost pumping power", the present model predicts that non-slender domains (aspect ratio a > 0.5) generate a much higher "cost" of the blood or sap distribution.Some consequences of the above derivation are noteworthy: (1) The bifurcation angle (2α 1 ) varies between 54 ˝and 68 ˝and is a function of the domain aspect ratio.This is not in agreement with Murray's cosine rule for the "optimal angle", that predicts 2α 1 = 75 ˝; (2) The radius ratio is a function of the splitting angle only, which is a remarkable result, because it was derived from a wall suction/drainage balance which in fact minimizes irreversible effects.Though much lower than that (0.7937) predicted by the H-M law, the values of δ 1 shown in Table 2 are within the range of the available experimental data: a series of indicative measurements on the images displayed in see Figure 6a,b,d show that the radius ratio of the first bifurcation is in the range 0.3 ˜0.75; (3) The more slender the domain (lower aspect ratio), the shorter λ 1 and the higher δ 1 .This implies that the pressure gradient in the main and in the daughter branches is also a function of the splitting angle, and (Table 2) it is higher in the latter, as predicted also by the H-M law; (4) The pressure drop is also a function of the aspect ratio and it does not assume equal values in the parent and daughter branches, its value being always higher in the latter: while for lower aspect ratios the total pressure drop is almost equally allocated, for square domains the ∆p 1 is 25 times higher than the ∆p 0 .Incidentally, this also negates the assumption of "equal pumping work" that is at the basis of the H-M law, but agrees with the dictates of Constructal Theory; (5) Considering that entropy generation rate per unit mass of transported fluid is proportional to the "lost pumping power", the present model predicts that non-slender domains (aspect ratio a > 0.5) generate a much higher "cost" of the blood or sap distribution.
Table 3 shows the "optimal" ratios and angles in the present and in some previous studies.

Conclusions
The Hess/Murray law was revisited on the basis of the most modern representations of capillary nets.If the structures are represented as bifurcations, it is possible to show that the "cubic root of 2" rule derived for the radius ratio from the H-M law is not the unique optimum that can be obtained under physically sound assumptions.In particular, there is no a priori reason for a natural structure to bifurcate under constant volume, constant Reynolds number, or constant velocity constraints (though this latter case is almost surely at odds with experimental observations).The original energy balance proposed by Hess & Murray is then scrutinized, and an alternative form based on the functionality of the capillary network is developed: this new energy balance shows that both the "optimal splitting angle" and the "optimal radius ratio" depend on the shape of the domain "served" by the bifurcation.In this new representation, the general features of the capillary structure depend thus on its function and are in complete agreement with an evolution-driven adaptive model based on a minimum dissipation principle: in nature, organisms show a tendency to minimize the energy dissipation by friction for a given delivered mass flowrate of nutrients; for artificial artifacts, engineers may follow the same guideline to design optimal distribution networks for material and energy flows.A more rational "economic" approach than Murray's "cost of blood" is certainly the amount of entropy generation per unit mass of delivered blood or sap.If the goal is to carry fluid to a "user", and the transportation is affected by irreversible entropy generation, then the most efficient transportation network is the one that displays the lowest entropy generation rate per unit mass of transported fluid.In thermodynamics, the degree of irreversibility is best measured by a parameter called exergy destruction, Ex δ , proportional to ∆S irr .Since exergy-based optimization procedures are well known in engineering and thermal sciences, it is likely that a systematic exergy analysis of fluid-carrying networks may shed more light on the evolution of living organisms, as it has on the design of energy systems and heat exchangers.
As a final remark, it is necessary to recall that the study presented herein was conducted assuming stationary flow and domain conditions: pulsating and compliant wall effects are therefore neglected.While this obviously detracts of the direct applicability of the theory to real flows, the existence of a rigorously based steady model may help future investigations.
(a) Rectangular domain (b) Triangular domain

Figure 2 .
Figure 2. Murray's original drawing for the derivation of his branching angles rule [5] (in the notation of this paper, x = α1, y = α2).

Figure 3 .
Figure 3. Minimum irreversibility ratio S/S0 as a function of the aspect ratio H/L and of the splitting angle α1 for different physical assumptions [29].

Figure 3 .
Figure 3. Minimum irreversibility ratio S/S0 as a function of the aspect ratio H/L and of the splitting angle α1 for different physical assumptions [29].

Figure 5 .
Figure 5. Schematic representation of the function of capillary networks in biological tissue.(a) functional scheme of a capillary network [32]; (b) its realistic representation (adapted from [16]).

Figure 5 .
Figure 5. Schematic representation of the function of capillary networks in biological tissue.(a) functional scheme of a capillary network [32]; (b) its realistic representation (adapted from [16]).

Figure 7 .
Figure 7.The two regions "served" by the main and by the daughter branches.

Figure 7 .
Figure 7.The two regions "served" by the main and by the daughter branches.

Figure 7 .
Figure 7.The two regions "served" by the main and by the daughter branches.

Figure 8 .
Figure 8. Optimal geometries for equal area drain.
Wall stress, N/m 2

Table 1 .
Viscous entropy generation rate in successive symmetrical branchings.

Table 2 .
Radius ratio δ 1 and optimal splitting angle for different domain aspect ratios.

Table 3 .
"Optimal" ratios and angles in the present and in some previous studies.Optimal geometries for equal area drain.

Table 3 .
"Optimal" ratios and angles in the present and in some previous studies.