Hurdle Theorem for Parallel Flows

• The paper reveals the main contribution: an analytical method that sharply bounds hitting probabilities and instability criteria in coupled parallel systems.
• It employs techniques such as Rayleigh quotient minimization and Sturm–Liouville analysis to derive precise thresholds for instability and rare events.
• Applications span fluid dynamics and queueing theory, offering concrete predictions for flow instability and overflow probabilities in multidimensional models.

The Hurdle Theorem for Parallel Flows provides a sharp analytical framework for understanding instability and rare-event probabilities in systems where two or more parallel processes interact under drift, boundary, and constraint conditions. It encompasses results in hydrodynamic stability, queueing theory, and optical transformations, unifying diverse problems under the common structure of barrier-crossing in multidimensional Markov models and Sturm–Liouville eigenvalue problems. The theorem yields both sufficient criteria for instability in fluid flows and precise asymptotic characterizations of rare event probabilities in parallel queueing models.

1. Mathematical Formulation in Parallel Queue Models

Consider the constrained random walk $X = \{X_k: k\ge 0\}$ on $\mathbb{Z}_+^2$ with increments $I_k \in \{(1,0), (-1,0), (0,1), (0,-1)\}$. The transitions are determined by arrival rates $\lambda_i$ and service rates $\mu_i$ for $i=1,2$, representing the dynamics of two parallel $M/M/1$ queues, with traffic intensities $\rho_i = \lambda_i/\mu_i < 1$.

Define the domain $$A_n = \{ x \in \mathbb{Z}_+^2 : x_1 + x_2 \leq n \}$$ and its boundary $\partial A_n$. The stopping times

$\mathbb{Z}_+^2$0

are used to study the overflow probability $\mathbb{Z}_+^2$1.

Analysis of $\mathbb{Z}_+^2$2 proceeds by mapping the problem to an auxiliary random walk $\mathbb{Z}_+^2$3 on $\mathbb{Z}_+^2$4 with reversed first-coordinate jump probabilities and reflection at $\mathbb{Z}_+^2$5. The hitting time of the diagonal $\mathbb{Z}_+^2$6, denoted $\mathbb{Z}_+^2$7, becomes the effective "hurdle" to cross.

The core result—the Hurdle Theorem for Parallel Flows—states that, for scaled starting states $\mathbb{Z}_+^2$8 ($\mathbb{Z}_+^2$9),

$I_k \in \{(1,0), (-1,0), (0,1), (0,-1)\}$0

for some $I_k \in \{(1,0), (-1,0), (0,1), (0,-1)\}$1 and all large $I_k \in \{(1,0), (-1,0), (0,1), (0,-1)\}$2, with $I_k \in \{(1,0), (-1,0), (0,1), (0,-1)\}$3 (Ünlü et al., 2018).

In the "critical regime" $I_k \in \{(1,0), (-1,0), (0,1), (0,-1)\}$4, $I_k \in \{(1,0), (-1,0), (0,1), (0,-1)\}$5, the diagonal-hitting probability is given explicitly by

$I_k \in \{(1,0), (-1,0), (0,1), (0,-1)\}$6

In the generic regime $I_k \in \{(1,0), (-1,0), (0,1), (0,-1)\}$7, $I_k \in \{(1,0), (-1,0), (0,1), (0,-1)\}$8 can be approximated arbitrarily well via superpositions of log-linear harmonic functions parameterized on the characteristic surface $I_k \in \{(1,0), (-1,0), (0,1), (0,-1)\}$9.

2. Analytical Criteria in Hydrodynamic Parallel Flows

The original hydrodynamic context for the hurdle theorem is the stability analysis of inviscid, incompressible, 2D parallel shear flows $\lambda_i$0 on a finite interval $\lambda_i$1. The Rayleigh eigenvalue problem

$\lambda_i$2

is studied for criteria under which instability arises.

Define the curvature function

$\lambda_i$3

at a critical level $\lambda_i$4 where $\lambda_i$5, and set $\lambda_i$6. The domain "hurdle"

$\lambda_i$7

corresponds to the square of the half-wavenumber for Dirichlet conditions. The Hurdle Theorem for parallel flows states: if $\lambda_i$8 for all $\lambda_i$9, a neutral mode exists at $\mu_i$0 for some $\mu_i$1, and a pair of unstable complex-conjugate modes emerges under perturbation, establishing inviscid instability (Deguchi et al., 25 Jan 2026, Deguchi et al., 2024).

The Rayleigh-quotient criterion asserts that finding a test function $\mu_i$2 with Rayleigh quotient $\mu_i$3 guarantees the existence of a neutral root; a Sturm–Liouville argument and perturbation in $\mu_i$4 demonstrate transition to instability.

3. Extension to Axisymmetric Annular and Pipe Flows

The hurdle theorem is generalized to axisymmetric flows in cylindrical coordinates. The base flow $\mu_i$5 and disturbances with axial wavenumber $\mu_i$6 and azimuthal number $\mu_i$7 lead to the inviscid stability equation

$\mu_i$8

where $\mu_i$9, $i=1,2$0. The generalized curvature

$i=1,2$1

serves as the analog of $i=1,2$2.

For annular domains $i=1,2$3, the theorem asserts: if $i=1,2$4 for all $i=1,2$5 in $i=1,2$6, with

$i=1,2$7

and criticality conditions on $i=1,2$8, then a neutral axisymmetric mode exists, and a nearby branch is unstable for perturbed $i=1,2$9 (Deguchi et al., 25 Jan 2026).

An analog applies to full pipe geometry, where the hurdle $M/M/1$0 is replaced with $M/M/1$1 using domain constants derived from geometric and spectral parameters.

4. Sturm–Liouville and Rayleigh-Quotient Methods

Underlying the hurdle theorem in all settings is the conversion of the linear stability or hitting-time problem to a Sturm–Liouville eigenvalue problem. The Rayleigh quotient

$M/M/1$2

or its axisymmetric generalization is minimized over admissible test functions. Instability or positive hitting probability corresponds to $M/M/1$3, i.e., a lowest eigenvalue below zero. The "hurdle" emerges from comparing with the first eigenvalue of a problem with constant potential over a subdomain, yielding the explicit $M/M/1$4 barrier. The existence of a test function (often the first Dirichlet eigenfunction) for which the Rayleigh quotient crosses this hurdle suffices to guarantee instability or exceeding rare-event probabilities.

In probabilistic parallel-queue systems, similar spectral logic governs the construction of harmonic functions—built from single or conjugate points on the process's characteristic surface—used to approximate hitting probabilities.

5. Physical and Probabilistic Interpretation

In fluid mechanical applications, $M/M/1$5 plays the role of a reciprocal local Rossby–Mach number; subsonic regions (everywhere $M/M/1$6) correspond to instability, while local supersonic patches ($M/M/1$7 somewhere) can yield stability, recovering the Kelvin–Arnol'd theorems in various limits (Deguchi et al., 25 Jan 2026, Deguchi et al., 2024). In the annular and pipe extensions, $M/M/1$8 encodes inviscid centrifugal effects and geometric confinement.

In queueing theory, the diagonal $M/M/1$9 acts as a "hurdle," with the hitting probability encoding the chance of excessive backlog. The exponential bound for the relative error in the main theorem demonstrates that in the many-server or large-buffer limit, the rare-event behavior is controlled by the much-simpler reflected process hitting a linear barrier (Ünlü et al., 2018).

6. Relations to Broader Stability and Instability Criteria

The hurdle theorem both sharpens and unifies existing sufficient conditions for instability, notably the Kelvin–Arnol'd I/II theorems and the classical Rayleigh and Tollmien criteria. It applies in a broader class of settings, including non-monotonic base flows, stratified and quasi-geostrophic fluids, magnetohydrodynamics (when the eigenproblem reduces to appropriate Sturm–Liouville form), and even alternating jet flows in planetary atmospheres (Deguchi et al., 2024).

Its predictions have been shown, in both queueing and hydrodynamic contexts, to closely bracket numerically computed or exact boundaries for instability and rare-event onset. In practical computations, it provides nearly sharp, analytically tractable regimes separating stable from unstable (or high-probability from rare-event) zones.

7. Illustrative Examples and Applications

• Hydrodynamic flows: For the sinusoidal base profile $\rho_i = \lambda_i/\mu_i < 1$0 with constant $\rho_i = \lambda_i/\mu_i < 1$1, the domain-wide hurdle $\rho_i = \lambda_i/\mu_i < 1$2 recovers exact agreement between the threshold for instability in $\rho_i = \lambda_i/\mu_i < 1$3 and numerically computed critical values (Deguchi et al., 2024).
• Parallel queue overflow: The critical regime closed-form for diagonal-hitting probability, and robust harmonic-function constructions in non-degenerate regimes, yield sharp asymptotics for overflow events in high-traffic parallel servers (Ünlü et al., 2018).
• Annular and pipe flows: Numerical comparison of the hurdle theorem's thresholds with full eigenvalue computations in model flows confirms that the theorem predicts the instability region with high accuracy; in particular, it outperforms older criteria under confinement or strong curvature (Deguchi et al., 25 Jan 2026).

Application Area	Instability Criterion	Reference
2D parallel shear flow	$\rho_i = \lambda_i/\mu_i < 1$4 for all $\rho_i = \lambda_i/\mu_i < 1$5	(Deguchi et al., 25 Jan 2026)
Annular/pipe flows	$\rho_i = \lambda_i/\mu_i < 1$6 for all $\rho_i = \lambda_i/\mu_i < 1$7	(Deguchi et al., 25 Jan 2026)
Parallel queues	Hitting $\rho_i = \lambda_i/\mu_i < 1$8 probability approximates $\rho_i = \lambda_i/\mu_i < 1$9	(Ünlü et al., 2018)
Alternating jets	$$A_n = \{ x \in \mathbb{Z}_+^2 : x_1 + x_2 \leq n \}$$0 on a subdomain	(Deguchi et al., 2024)

The hurdle theorem thus serves as a cornerstone for modern analysis of instability and rare events in parallel flows of both physical and stochastic origin, connecting spectral theory, probability, and applied mechanics.