Computational code for "Generalized Elasticity Principle for Optimal Branching in Tree-Like Duct Networks

This repository contains the computational code that reproduces all numerical results, figures, and validation tables reported in: H. Mishra, "Generalized Elasticity Principle for Optimal Branching in Tree-Like Duct Networks," submitted to Physics of Fluids (2026). The code implements: The iterated-map formulation for optimal branching in constant-height duct networks (rectangular, elliptical, rhombic, and trapezoidal cross-sections) The body-fitted Poisson solver for trapezoidal ducts with validation against Shah and London (1978) Computation of the generalized elasticity, effective branching exponent, and optimal width ratio across all duct families Direct multi-variable constrained optimization for independent validation of the iterated map Murray's-law penalty calculations (pressure-drop comparison at matched volume) Power-law fluid extension for flow indices n = 0.5, 0.75, 1, 1.5, 2 All figures (Figs. 1–5) and tables (Tables I–VII) in the manuscript