The 2-component dispersionless Burgers equation arising in the modelling of blood flow

This article investigates the properties of the solutions of the dispersionless two-component Burgers (B2) equation, derived as a model for blood-flow in arteries with elastic walls. The phenomenon of wave breaking is investigated as well as applications of the model to clinical conditions.

💡 Research Summary

The paper introduces a two‑component dispersionless Burgers (B2) system as a simplified yet powerful model for blood flow in arteries with elastic walls. Starting from the one‑dimensional continuity and momentum equations, the authors replace the cross‑sectional area with the square of the arterial radius and assume a nonlinear elastic relation between pressure and radius. By neglecting dispersive and viscous terms, the governing equations reduce to a pair of coupled first‑order hyper‑bolic equations:
 u_t + u u_x = 0, r_t + u r_x = 0,
where u(x,t) denotes the axial velocity and r(x,t) the arterial radius (or, equivalently, pressure). This system is the two‑component analogue of the inviscid Burgers equation and captures the essential nonlinear advection of both flow speed and wall deformation.

The authors perform a rigorous characteristic analysis. Along a characteristic curve defined by dx/dt = u, the velocity u remains constant, while the radius evolves according to dr/dt = u r_x. The solution can be expressed implicitly in terms of the initial data (u₀(x), r₀(x)). Wave breaking, or gradient catastrophe, occurs when characteristics intersect, causing ∂_x u and ∂_x r to blow up. The breaking time is derived analytically as
 t* = –1 / min_x (∂_x u₀(x)),
and an analogous condition is obtained for the radius. Thus, any sufficiently steep negative slope in the initial velocity profile—or a rapid decrease in the initial radius—inevitably leads to finite‑time singularity.

To assess physiological relevance, three numerical experiments are presented. (1) A baseline case with normal blood pressure (120/80 mmHg) and typical arterial stiffness shows smooth wave propagation without breaking. (2) A hypertensive scenario (systolic pressure ≈ 180 mmHg) imposes a steeper initial pressure gradient; characteristics converge quickly, and the breaking time drops to about 0.03 s. The resulting pressure waveform develops a near‑shock, mimicking the abrupt pressure spikes observed during hypertensive crises. (3) An arterial‑stiffening case, modeled by reducing the elastic modulus by 30 %, also accelerates characteristic crossing, leading to amplified wave reflections and earlier breaking. These simulations illustrate how pathological changes—high pressure or reduced compliance—enhance the likelihood of wave breaking, which in turn can generate the high‑frequency pressure transients implicated in endothelial injury and aneurysm rupture.

The discussion acknowledges the model’s limitations. The dispersionless assumption discards high‑frequency components that become important in small vessels or during rapid wall motion. Consequently, the authors propose extensions that re‑introduce a small viscous term (ν u_xx) or a dispersive correction, yielding a “regularized” B2 system capable of capturing shock‑like structures while preserving physical realism. They also suggest coupling the B2 equations with multi‑dimensional arterial networks, incorporating branching, curvature, and time‑varying boundary conditions to move toward patient‑specific simulations.

In conclusion, the two‑component dispersionless Burgers framework offers a mathematically tractable yet physiologically insightful description of arterial wave dynamics. It predicts finite‑time gradient catastrophes that correspond to clinically observed rapid pressure rises, providing a potential tool for assessing the risk of vascular events such as acute hypertension spikes or aneurysm rupture. Future work aimed at regularization, higher‑dimensional extensions, and validation against clinical measurements could transform this theoretical construct into a practical diagnostic aid.