# Murray's law

Murray's law predicts the thickness of branches in transport networks, such that the cost for transport and maintenance of the transport medium is minimized. This law is observed in the vascular and respiratory systems of animals, xylem in plants, and the respiratory system of insects. Its simplest version states that if a branch of radius $r$ splits into two branches of radii $r_{1}$ and $r_{2}$ , then $r^{3}=r_{1}^{3}+r_{2}^{3}$ . Like Hagen–Poiseuille equation and Fick's laws, which were also formulated from a biological context, Murray's law is a basic physical principle for transfer networks.

Murray's law is also a powerful biomimetics design tool in engineering. It has been applied in the design of self-healing materials, batteries, photocatalysts, and gas sensors.[citation needed] However, since its discovery, little attention has been paid to exploit this law for designing advanced materials, reactors and industrial processes for maximizing mass or energy transfer to improve material performance and process efficiency.

The original version of Murray's law is only applicable to mass-conservative transport in the network. There are generalizations for non-conservative networks, which describe effects such as chemical reactions and diffusion through the walls.

## Murray's law in mass-conservative networks

Murray's original analysis is based on the assumption that the radii inside a lumen-based system are such that the work for transport and upkeep is minimized. Larger vessels lower the expended energy for transport, but increase the overall volume of blood in the system; blood being a living fluid and hence requiring metabolic support. Murray's law is therefore an optimisation exercise to balance these factors.

For $n$  child branches splitting from a common parent branch, the law states that:

$r^{3}=r_{1}^{3}+r_{2}^{3}+r_{3}^{3}+...+r_{n}^{3}$

where $r$  is the radius of the parent branch, and $r_{1},\ldots ,r_{n}$  are the radii of the child branches. This law is only valid for laminar flow, since its derivation uses the Hagen–Poiseuille equation as a measure for the transport work (see below). Williams et al. deduced the formula for turbulent flow:

$r^{(7/3)}=r_{1}^{(7/3)}+r_{2}^{(7/3)}+r_{3}^{(7/3)}+...+r_{n}^{(7/3)}$

### Derivation

As stated above, the underlying assumption of Murray's law is that the power (energy per time) for transport and upkeep is minimal in a natural transport system. Hence we seek to minimize $P=P_{t}+P_{m}$ , where $P_{t}$  is the power required for transport and $P_{m}$  the power required to maintain the transport medium (e.g. blood).

#### Laminar flow

We first minimize the power of transport and upkeep in a single channel of the system, i.e. ignore bifurcations. Together with the assumption of mass-conservation this will yield the law. Let $Q$  be the laminar flow rate in this channel, which is assumed to be fixed. The power for transport in a laminar flow is $P_{t}=Q\Delta p$ , where $\Delta p$  is the pressure difference between the entry and exit of a tube of radius $r$  and length $l$ . The Hagen–Poiseuille equation for laminar flow states that ${\textstyle Q={\frac {\pi \,r^{4}}{8\,\mu \,l}}\Delta p}$  and hence ${\textstyle \Delta p={\frac {8\,\mu \,l}{\pi \,r^{4}}}Q}$  , where $\mu$  is the dynamic viscosity of the fluid. Substituting this into the equation for $P_{t}$ , we get

$P_{t}={\frac {8\,\mu \,l}{\pi \,r^{4}}}Q^{2}.$

We furthermore make the assumption that the power necessary for maintenance is proportional to the volume of the cylinder $V=\pi \,l\,r^{2}$  :
$P_{m}=\lambda \,V=\lambda \,\pi \,l\,r^{2}$

where $\lambda >0$  is a constant number (metabolic factor). This yields
$P=P_{t}+P_{m}={\frac {8\,\mu \,l}{\pi \,r^{4}}}Q^{2}+\lambda \,\pi \,l\,r^{2}.$

Since we seek the radius with minimal power, we calculate the stationary point with respect to r
${\operatorname {d} \over \operatorname {d} \!r}P=0\;\Leftrightarrow \;-{\frac {32\,\mu \,l}{\pi \,r^{5}}}Q^{2}+2\,\lambda \,\pi \,l\,r=0\;\Leftrightarrow \;Q^{2}={\frac {\pi ^{2}}{16}}{\frac {\lambda }{\mu }}r^{6}\;\Leftarrow \;Q={\frac {\pi }{4}}{\sqrt {\frac {\lambda }{\mu }}}r^{3}.$

With this relationship between the flow rate $Q$  and radius $r$  on an arbitrary channel, we return to the bifurcation: Since mass is conserved, the flow rate of the parent branch $Q$  is the sum of the flow rates in the children branches $Q_{1},\ldots ,Q_{n}$ . Since in each of these branches the power is minimized, the relationship above holds and hence, as claimed:
$Q=\sum _{i=1}^{N}Q_{i}\;\Rightarrow \;{\frac {\pi }{4}}{\sqrt {\frac {\lambda }{\mu }}}\,r^{3}=\sum _{i=1}^{N}{\frac {\pi }{4}}{\sqrt {\frac {\lambda }{\mu }}}\,r_{i}^{3}\;\Rightarrow \;r^{3}=\sum _{i=1}^{N}r_{i}^{3}.$

#### Diffusion

The power for transport spent by diffusion is given by

$P_{t}=Q\,\Delta C=\left(\pi \,D\,r^{2}{\frac {\Delta C}{l}}\right)\Delta C$

where the flow rate is given by Fick's law, whose $D$  is the diffusivity constant and $\Delta C$  is the difference of concentration between the ends of the cylinder. Similarly to the case of laminar flow, the minimisation of the objective function results in

${\frac {\Delta C}{l}}={\sqrt {\frac {\lambda }{D}}}\Rightarrow Q=\pi \,D\,r^{2}{\sqrt {\frac {\lambda }{D}}}\Leftrightarrow Q=k\,r^{2}$

Hence,

$Q_{p}=\sum _{i=1}^{N}Q_{i}\Rightarrow k\,r_{p}^{2}=\sum _{i=1}^{N}k\,r_{i}^{2}\Leftrightarrow r_{p}^{2}=\sum _{i=1}^{N}r_{i}^{2}$

## The generalized Murray’s law

However, the special Murray’s law is only applicable to flow processes involving no mass variations. Significant theoretical advances need to be made for more broadly applicable in the fields of chemistry, applied materials, and industrial reactions.

The generalized Murray's law deduced by Zheng et al. can be applicable for optimizing mass transfer involving mass variations and chemical reactions involving flow processes, molecule or ion diffusion, etc.

For connecting a parent pipe with radius of r0 to many children pipes with radius of ri , the formula of generalized Murray's law is: $r_{0}^{a}={1 \over 1-X}\sum _{i=1}^{N}r_{i}^{a}$ , where the X is the ratio of mass variation during mass transfer in the parent pore, the exponent α is dependent on the type of the transfer. For laminar flow α = 3; for turbulent flow α = 7/3; for molecule or ionic diffusion α = 2; etc.

## Murray materials

The generalized Murray’s law defines the basic geometric features for porous materials with optimum transfer properties. The generalized Murray’s law can be used to design and optimize the structures of an enormous range of porous materials. This concept has led to materials, termed as the Murray materials, whose pore-sizes are multiscale and are designed with diameter-ratios obeying the generalized Murray’s law.

A diagram of Murray materials with macro-meso-micropores built by nanopaticles as building blocks. (Zheng, CC-BY-SA)

As lithium battery electrodes, the Murray materials can reduce the stresses in these electrodes during the charge/discharge processes, improving their structural stability and resulting in a longer life time for energy storage devices . This material could also be used for boosting the performance of a gas sensor and a photocatalysis process that broke down a dye using light.

To achieve substances or energy transfer with extremely high efficiency, evolution by natural selection has endowed many classes of organisms with Murray materials, in which the pore-sizes regularly decrease across multiple scales and finally terminate in size-invariant units. For example, in plant stems and leaf veins, the sum of the radii cubed remains constant across every branch point to maximize the flow conductance, which is proportional to the rate of photosynthesis. For insects relying upon gas diffusion for breathing, the sum of radii squared of tracheal pores remains constant along the diffusion pathway, to maximize gases diffusion. From plants, animals and materials to industrial processes, the introduction of Murray material concept to industrial reactions can revolutionize the design of reactors with highly enhanced efficiency, minimum energy, time, and raw material consumption for a sustainable future.