Radius exponent in elastic and rigid arterial models optimized by the least energy principle

Abstract

It was analyzed in normal physiological arteries whether the least energy principle would suffice to account for the radius exponent x. The mammalian arterial system was modeled as two types, the elastic or the rigid, to which Bernoulli's and Hagen-Poiseuille's equations were applied, respectively. We minimized the total energy function E, which was defined as the sum of kinetic, pressure, metabolic and thermal energies, and loss of each per unit time in a single artery transporting viscous incompressible blood. Assuming a scaling exponent α between the vessel radius (r) and length (l) to be 1.0, x resulted in 2.33 in the elastic model. The rigid model provided a continuously changing x from 2.33 to 3.0, which corresponded to Uylings' and Murray's theories, respectively, through a function combining Reynolds number with a proportional coefficient of the l - r relationship. These results were expanded to an asymmetric arterial fractal tree with the blood flow preservation rule. While x in the optimal elastic model accounted for around 2.3 in proximal systemic (r >1 mm) and whole pulmonary arteries (r ≥0.004 mm), optimal x in the rigid model explained 2.7 in elastic-muscular (0.1 < r ≤1 mm) and 3.0 in peripheral resistive systemic arteries (0.004 ≤ r ≤0.1 mm), in agreement with data obtained from angiographic, cast-morphometric, and in vivo experimental studies in the literature. The least energy principle on the total energy basis provides an alternate concept of optimality relating to mammalian arterial fractal dimensions under α = 1.0.