🔢 mathematics

A Unified Model for Blood and Lymph Flow with Coupled Nonsmooth Biochemical Dynamics

This paper presents a unified mathematical framework that couples partial differential equations for blood and lymph fluid dynamics with non-smooth ordinary differential equations modeling calcium and nitric oxide kinetics to rigorously derive and validate the stable oscillatory pumping behavior of lymphangions.

Original authors: Bogna Jaszczak-Dyka, Łukasz Płociniczak

Published 2026-03-03

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Bogna Jaszczak-Dyka, Łukasz Płociniczak

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine your body as a bustling city. You know about the arteries and veins, which are like the main highways carrying blood (the delivery trucks) pumped by the heart (the central power plant). But there's a second, quieter network: the lymphatic system. Think of this as the city's sewage and recycling system. It collects excess fluid and waste from the streets (your tissues) and returns it to the main water supply.

Unlike the blood system, which has a giant pump (the heart) pushing everything forward, the lymphatic system has no central engine. Instead, it relies on thousands of tiny, self-contained "pumps" called lymphangions. These are like little segments of a garden hose that squeeze themselves to push water forward.

This paper is a mathematical blueprint for how these tiny pumps work, specifically focusing on how they decide when to squeeze and when to relax.

The Problem: How does a pump know when to work?

In a mechanical pump, you might use a timer or a sensor. In the body, the lymphangion uses a chemical conversation between two key players:

Calcium (Ca²⁺): The "muscle builder." When levels are high, the vessel wall squeezes tight (contracts).
Nitric Oxide (NO): The "muscle relaxer." When levels are high, the vessel wall loosens up.

The paper explains that these two chemicals are locked in a dance.

The Squeeze: When the vessel squeezes, it pushes fluid out. This flow creates friction (shear stress) against the vessel walls.
The Signal: That friction tells the wall to produce Nitric Oxide.
The Relax: The Nitric Oxide kills the Calcium, causing the vessel to relax and open up.
The Cycle: As it relaxes, fluid rushes in, friction drops, Nitric Oxide production stops, Calcium builds up again, and the vessel squeezes once more.

The Innovation: A Unified Mathematical Model

The authors built a super-accurate simulation of this process. Here is how they did it, broken down simply:

1. The Fluid Mechanics (The "Hose" Part)

They started with the fundamental laws of physics (Navier-Stokes equations) to describe how fluid moves through a flexible tube.

Analogy: Imagine a flexible garden hose lying on the ground. If you squeeze it in the middle, the water shoots out the end. The authors wrote math to predict exactly how the hose deforms and how fast the water moves, accounting for the fact that lymph is a bit like water, but blood is thicker (like a slurry).
The Twist: They didn't just look at the fluid; they looked at how the walls of the vessel react to the pressure, treating the vessel like a stretchy balloon.

2. The Biochemistry (The "Brain" Part)

This is where the paper gets really clever. They realized that the valve opening and closing isn't just mechanical; it's chemical.

They created a set of equations for Calcium and Nitric Oxide.
The "Non-Smooth" Surprise: In real life, valves don't open and close smoothly like a dimmer switch; they snap open or shut like a light switch. The math reflects this "jagged" behavior (called non-smooth dynamics). It's like a traffic light that is either strictly Red or strictly Green, with no "yellow" in between for the valve state.

3. The Connection (The "Unified" Part)

Most previous studies looked at the fluid or the chemicals separately. This paper glues them together.

They showed how the fluid flow creates the chemical signal (friction -> NO).
They showed how the chemical signal changes the fluid flow (NO -> relaxation -> more flow).
The Result: A complete loop. The math proves that if you set the parameters right (like the right amount of Calcium sensitivity), this system naturally starts oscillating. It doesn't need a timer; the chemistry and physics create their own rhythm.

Why Does This Matter?

The authors found that under specific conditions, this system creates a Limit Cycle.

Analogy: Think of a swing. If you push it at the right time, it keeps swinging forever without you needing to push it again. The lymphangion is that swing. The paper mathematically proves why the lymphatic system keeps pumping rhythmically on its own.

The Big Picture:
This research is like finding the instruction manual for the body's "self-cleaning" system.

For Healthy People: It explains how our bodies prevent swelling (edema) and keep tissues healthy.
For Sick People: If this chemical dance gets out of sync (e.g., too much Calcium or not enough Nitric Oxide), the pumps stop working. This leads to lymphedema (chronic swelling), which is common in cancer patients or after surgery.

By understanding the exact math behind this rhythm, doctors and engineers might one day design better drugs or devices to fix broken lymphatic pumps, helping people with swelling issues get their "recycling system" back on track.

In short: The paper takes the complex physics of fluid flow and the chemistry of cell signals, mixes them into a single mathematical recipe, and proves that this recipe naturally creates the rhythmic pumping that keeps our bodies from swelling up.

1. Problem Statement

The lymphatic system is a critical but historically understudied component of human physiology, responsible for fluid homeostasis, immune function, and lipid transport. Unlike the cardiovascular system, which relies on a central pump (the heart), the lymphatic system operates as an open network driven by intrinsic rhythmic contractions of lymphangions (segmented contractile chambers).

The Challenge: Existing mathematical models often lack a unified derivation from first physical principles or fail to adequately couple the fluid dynamics of lymph flow with the complex, non-smooth biochemical feedback loops (involving Calcium ions and Nitric Oxide) that regulate valve opening and closing.
The Goal: To develop a rigorous, unified mathematical framework that derives the governing equations for lymph flow from the Navier-Stokes equations and couples them with a non-smooth dynamical system describing the biochemical kinetics of $Ca^{2+}$ and NO, ultimately proving the existence of self-sustained oscillations (limit cycles) that mimic physiological pumping.

2. Methodology

The authors employ a multi-stage analytical approach combining asymptotic analysis, fluid mechanics, and non-smooth dynamical systems theory.

A. Derivation of the Flow Model (Fluid Dynamics)

Starting Point: The authors begin with the full Navier-Stokes equations for viscous, axisymmetric, laminar flow in a cylindrical tube with a dynamic boundary $R(x,t)$ .
Assumptions: The fluid is treated as Newtonian (a valid approximation for lymph in larger vessels), and the flow is driven primarily by pressure gradients.
Scaling and Lubrication Approximation: Using perturbation methods, they introduce a small aspect ratio parameter $\epsilon = R_0/L \ll 1$ (where $R_0$ is the radius and $L$ is the length). They perform non-dimensionalization and apply the lubrication approximation (retaining leading-order terms).
Constitutive Law: A pressure-radius relationship is derived to account for vessel wall elasticity and tension, incorporating empirical data for lymphangions.
Averaging: The equations are averaged over the cross-sectional area to derive a closed system of two nonlinear Partial Differential Equations (PDEs) governing the mean velocity $U$ and the vessel radius $R$ (or cross-sectional area $A$ and flux $Q$ ). This system resembles the Saint-Venant shallow water equations but includes specific terms for lymphatic compliance and wall tension.

B. Biochemical Dynamics Model (Valve Regulation)

Mechanism: The model focuses on the mechano-chemical oscillator governing lymphangion contraction.
- Calcium ( $Ca^{2+}$ ): Drives contraction. Its dynamics include decay, voltage-dependent influx, stretch-activated channels, and inhibition by Nitric Oxide.
- Nitric Oxide (NO): Drives relaxation. Produced in response to wall shear stress (which depends on flow velocity) and $Ca^{2+}$ levels.
Non-Smooth Formulation: The interaction is modeled using a system of two Ordinary Differential Equations (ODEs) containing Heaviside step functions ( $H$ ). These functions represent threshold-based switching (e.g., channel opening only when concentration exceeds a specific value), making the system non-smooth.
Simplification: To analyze the core oscillatory mechanism, the authors decouple the biochemical system from the fluid PDEs by assuming a constant vessel radius ( $R = \text{const}$ ). This reduces the problem to analyzing the 2D non-smooth ODE system for $Ca^{2+}$ ( $C$ ) and NO ( $N$ ).

C. Mathematical Analysis

Filippov Systems Theory: The reduced ODE system is analyzed as a piecewise-smooth dynamical system. The state space is divided by a discontinuity boundary ( $C=1$ ).
Equilibrium Classification: The authors classify equilibria as admissible, virtual, boundary, or pseudo-equilibria based on their location relative to the switching boundary.
Bifurcation Analysis: They perform a rigorous bifurcation analysis to identify parameter regimes (specifically relating to production/degradation rates $\alpha, \beta, \gamma, \zeta$ ) where the system transitions from stable fixed points to periodic behavior.
Limit Cycle Proof: Using a Poincaré map constructed on the switching boundary, they prove the existence of a stable limit cycle under specific conditions where both equilibria are virtual and stable.

3. Key Contributions

Unified Derivation from First Principles: Unlike previous studies that often use ad-hoc assumptions, this paper rigorously derives the lymph flow equations from the Navier-Stokes equations using asymptotic perturbation theory, providing a unified framework applicable to both blood and lymph (distinguished only by rheological parameters).
Coupling Fluid Mechanics with Non-Smooth Biochemistry: The paper introduces a novel coupling between the fluid dynamics and a non-smooth biochemical model. It explicitly models the valve opening/closing as a threshold-driven process dependent on $Ca^{2+}$ and NO concentrations.
Rigorous Proof of Oscillations: The authors provide a mathematical proof (Theorem 1) that the simplified non-smooth system admits a stable limit cycle. This validates the hypothesis that intrinsic biochemical feedback loops are sufficient to generate the rhythmic pumping observed in lymphangions without external forcing.
Bifurcation Characterization: The study maps out the parameter space, identifying specific regimes (e.g., when $\beta/\alpha > 1$ and $(\beta+\gamma)/(\alpha(\zeta+1)) < 1$ ) that guarantee oscillatory behavior, and classifies the types of boundary equilibrium bifurcations (persistence vs. non-smooth fold) that occur at transitions.

4. Results

Fluid Model: The derived governing equations successfully describe the spatio-temporal evolution of lymph flow, pressure, and vessel radius, capturing the interplay between transmural pressure, wall tension, and viscous forces.
Oscillatory Dynamics: The analysis of the reduced ODE system confirms that the interplay between $Ca^{2+}$ $C a^{2 +}$ -induced contraction and NO-induced relaxation creates a self-sustaining relaxation oscillator.
- Phase Plane Analysis: The phase portraits show trajectories moving between two virtual equilibria, crossing the discontinuity boundary, and forming a closed loop (limit cycle).
- Analytical Solutions: The authors derive closed-form expressions for the flow trajectories in both regions of the phase space using incomplete gamma functions, which are used to construct the Poincaré map.
Parameter Sensitivity: The study identifies that the existence of the limit cycle is robust within a specific range of dimensionless parameters, aligning with physiological estimates for lymphatic vessels.

5. Significance

Theoretical Validation: The work provides a deterministic mathematical explanation for the "relaxation oscillations" observed in lymphatic pumping, bridging the gap between microscopic biochemical kinetics and macroscopic fluid transport.
Physiological Insight: By identifying the parameter regimes for stable oscillations, the model offers insights into how lymphatic dysfunction (e.g., lymphedema) might arise from perturbations in biochemical parameters (e.g., altered NO production or calcium sensitivity).
Foundation for Future Research: The unified framework serves as a baseline for future studies on pathological conditions, therapeutic interventions, and the fully coupled fluid-biochemistry system (where the radius $R$ is dynamic). It opens the door for investigating nonlocal equations and complex network-scale simulations of the lymphatic system.

In summary, this paper presents a mathematically rigorous and physiologically grounded model that successfully unifies fluid dynamics and non-smooth biochemical kinetics to explain the self-regulating pumping mechanism of the lymphatic system.