An improved model of an actively contracting lymphatic vessel composed of several lymphangions: pumping characteristics

Using essentially our 2011 numerical model of a multi-lymphangion segment of a collecting lymphatic vessel, but augmented by inclusion of a refractory period and definition of a mid lymphangion pressure, we explore the effect of several parameters on the form of pump function curves. Pump function is sensitively dependent on the shape of the passive constitutive relation between lymphangion diameter and transmural pressure. Maximum flow-rate increases with the diameter scale applied to the constitutive relation and decreases with the pressure scale. Both maximum flow-rate and maximum pressure difference which can be overcome increase as the excess of lymphangion chain inlet pressure over external pressure is reduced, until inlet pressure is low enough that lymphangion collapse intervenes. The results are discussed in comparison with findings from biological experiments.

💡 Research Summary

This paper presents an enhanced numerical model of a actively contracting lymphatic vessel composed of multiple lymphangions, building directly on the authors’ 2011 multi‑lymphangion framework. Two key physiological refinements are introduced: (1) a refractory period following each active contraction, during which a lymphangion cannot be re‑stimulated, and (2) an explicit definition of a “mid‑lymphangion pressure” that represents the pressure in the segment between the upstream and downstream valves. These additions allow the model to capture the asymmetric pressure transmission, valve timing, and non‑linear contraction‑relaxation dynamics observed in real lymphatic vessels.

The core of the model remains the passive constitutive relation linking lymphangion diameter (D) to transmural pressure (ΔP). The authors systematically vary the shape of this D‑ΔP curve and introduce two scaling parameters: a diameter scale (α) and a pressure scale (β). Simulations reveal that increasing α (larger effective diameters for a given pressure) raises the maximal flow rate (Qmax) that the vessel chain can sustain, whereas increasing β (greater pressure required for the same diameter change) depresses Qmax and reduces overall pump efficiency. Thus, the passive mechanical properties of the vessel wall dominate the shape of the pump function curve (ΔPin versus Q).

A second set of investigations examines the effect of the excess inlet pressure (Pin – Pext). When this excess is reduced, each lymphangion operates closer to its optimal pressure range, allowing both Qmax and the maximal pressure difference the chain can overcome (ΔPmax) to increase simultaneously. However, if Pin falls below a critical threshold, lymphangion collapse occurs, leading to a precipitous drop in flow and eventual pump failure. This collapse phenomenon mirrors experimental observations of low‑pressure lymphatic failure.

Model outputs are compared with published experimental data on lymphatic pumping, particularly regarding valve opening/closing timing and the duration of the refractory period. The agreement supports the physiological relevance of the added features. Moreover, the study emphasizes that accurate knowledge of the passive D‑ΔP relationship is essential for predicting pump performance; small changes in its curvature or scaling can dramatically alter flow and pressure capabilities.

The authors discuss several practical implications. Clinically, manipulating the excess inlet pressure or augmenting lymphangion diameter (e.g., via pharmacological smooth‑muscle relaxation) could improve lymph transport in conditions such as lymphedema. From an engineering perspective, the findings guide the design of biomimetic lymphatic pumps: incorporating a refractory interval and a mid‑segment pressure sensor into control algorithms could enhance efficiency and stability. Finally, the work calls for refined experimental protocols to measure the passive D‑ΔP curve in isolated lymphangions, which would enable more precise parameterization of future models.

In summary, by integrating a refractory period and a mid‑lymphangion pressure definition into an established multi‑lymphangion framework, the paper provides a comprehensive analysis of how vessel wall mechanics, pressure scaling, and inlet‑external pressure differentials shape lymphatic pump function. The results deepen our mechanistic understanding of lymph propulsion, offer testable predictions for physiological experiments, and lay a solid theoretical foundation for therapeutic and bio‑engineering applications targeting the lymphatic system.