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

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 }_{\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 }^{\ast }=2.72$, within $0.1\sigma$ of morphometric data. Sensitivity analysis confirms $|\Delta {\alpha }^{\ast }|<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.

26 pages, 4 images, https://zenodo.org/records/19027393 and supplement material available