A Unified Variational Principle for Branching Transport Networks: Wave Impedance, Viscous Flow, and Tissue Metabolism

TL;DR

This paper develops a unified variational framework for biological transport networks, explaining observed diameter scaling and structural features by balancing wave impedance, viscous flow, and tissue metabolism.

Contribution

It introduces a new network-level optimization model that predicts the observed diameter scaling exponent and structural properties of branching networks.

Findings

Predicted diameter scaling exponent $oldsymbol{eta ext{ around } 2.72}$ matches empirical data.

Derived binary branching as an optimal structural feature.

Quantified wave dissipation consistent with clinical estimates.

Abstract

The branching geometry of biological transport networks is characterized by a diameter scaling exponent $\alpha$. Two structural attractors compete: impedance matching ($\alpha \sim 2$) for pulsatile flow and viscous-metabolic minimization ($\alpha = 3$) for steady flow. Neither predicts the empirically observed $\alpha_{\mathrm{exp}} = 2.70 \pm 0.20$ in mammalian arterial trees. Incorporating sub-linear vessel-wall scaling $h(r) \propto r^p$ ($p = 0.77$) into a three-term metabolic cost rigorously breaks Murray's cubic law -- via Cauchy's functional equation -- bounding the static optimum to $\alpha_t \in [2.90, 2.94]$. We formulate a unified network-level Lagrangian balancing wave-reflection penalties against transport-metabolic costs. Because the operational duty cycle $\eta$ is uncertain over developmental timescales, we cast the optimization as a zero-sum game between network architecture and environment. Von Neumann's minimax theorem -- proved via strict monotonicity of the cost curves -- yields a unique saddle point $(\alpha^, \eta^)$ satisfying an exact equal-cost condition. We further prove $N = 2$ uniquely maximizes the network stiffness ratio $\kappa_{\mathrm{eff}}(N)$, deriving binary branching as a structural consequence of the framework. For the porcine coronary tree ($G = 11$ generations), $\alpha^* = 2.72$, within $0.1\sigma$ of morphometric data. Sensitivity analysis confirms $|\Delta\alpha^*| < 0.01$ across physiological metabolic ranges; the prediction depends critically only on the histological exponent $p$ -- a zero-parameter derivation from fundamental scaling principles that simultaneously recovers a cumulative wave dissipation of 6.3%, consistent with independent clinical estimates.