---
title: Transition Flow Dynamics
url: https://www.emergentmind.com/topics/transition-flow
type: topic
---

# Transition Flow Dynamics

Transition flow denotes the flow state or dynamical evolution through which a system passes between qualitatively distinct regimes. In fluid mechanics, the term most commonly refers to finite-amplitude, often localized dynamics that connect laminar and turbulent states, such as puffs, slugs, oblique bands, laminar gaps, or time-periodic post-bifurcation motions. In other settings it denotes regime-changing flows produced by pulsation, vibration, interfacial breakup, topological rearrangements, or imposed parameter jumps. A distinct modern usage appears in generative modeling, where “transition flow” is defined as a learned mapping from an intermediate state at time \(t\) to a later state at time \(r\) [1706.10075], [1709.03738], [2603.15689].

## 1. Canonical meaning in wall-bounded shear flows

In wall-bounded shear flows, transition flow is classically associated with a globally subcritical scenario in which laminar and turbulent domains coexist over a finite Reynolds-number interval. Manneville’s review characterizes this interval by an upper threshold \(R_t\), above which turbulence is uniform, and a lower threshold \(R_g\), below which all turbulence eventually decays. For \(R_g<R<R_t\), the statistically steady state consists of laminar and turbulent coexistence, typically in the form of oblique bands. Near \(R_g\), the decay of these bands is described in terms of directed percolation, with turbulent fraction \(F_t \sim (R-R_g)^\beta\), spatial correlation length \(\xi \sim |R-R_g|^{-\nu_\perp}\), and temporal correlation \(\tau \sim |R-R_g|^{-\nu_\parallel}\), using the two-dimensional directed-percolation exponents \(\beta \approx 0.58\), \(\nu_\perp \approx 0.73\), and \(\nu_\parallel \approx 1.29\) [1706.10075].

Barkley’s reduced description places this phenomenology in a two-variable excitable-medium framework with turbulence intensity \(q(x,t)\) and mean shear \(u(x,t)\). In the continuous model, the laminar state is excitable for \(r<r_c=\epsilon_2/(\epsilon_1+\epsilon_2)\), which yields localized puff-like pulses, while for \(r>r_c\) the system becomes bistable and supports expanding slugs. The discrete extension embeds a chaotic repeller and reproduces finite puff lifetimes, puff splitting, sustained spatiotemporal intermittency, and a directed-percolation-like onset of sustained turbulence at \(R_c \simeq 2046.2\), with order-parameter scaling \(F_t \sim (R-R_c)^{\beta_{\rm DP}}\) and \(\beta \approx 0.28\) [1101.4125].

The dynamical-state inventory has been extended further by Frishman and Grafke, who introduced the antipuff and the gap-edge as counterparts to the puff and decay edge state. In their Barkley-model landscape, laminar gaps inside homogeneous turbulence are interpreted as antipuffs that nucleate and decay through the gap-edge, producing a symmetric picture in which localized turbulent states and localized laminar states each organize part of the transitional dynamics [2111.00233].

A complementary unification was proposed by Tao and Xiong through a locally defined Reynolds number,
\[
Re_L \equiv \frac{\max_{x_2}\bigl[U_0(x_2)\,\partial_{x_2}U_0(x_2)\bigr]\,R_h}{\tau_{w,0}},
\]
which is interpreted as the maximum local energy-input rate relative to local dissipation. Using \(Re_L\), they identified three transition stages common to linearly stable shear flows: equilibrium localized turbulence, temporally persistent turbulence, and uniform turbulence. In plane-Poiseuille flow, these occur at approximately \(Re_{L,1} \simeq 330.7\), \(Re_{L,2} \simeq 385 \pm 19\), and \(Re_{L,3} \simeq 519.6 \pm 19\), and the corresponding threshold bands collapse closely for pipe, channel, and Couette flows [1710.02258]. A related comparison of fluid-element-based criteria showed that Hanks’s \(K\), Tao’s \(Re_m\), and \(Re_L\) reduce to the same algebraic forms in Hagen–Poiseuille and plane–Poiseuille flows, while \(Re_m\) and \(Re_L\) remain well defined in plane Couette flow [2311.01067].

## 2. Thresholds and time-scale competition in pulsating pipe flow

In pulsating pipe flow, transition flow is controlled not only by Reynolds number but also by the interaction between turbulence dynamics and the pulsation period. The relevant dimensionless groups are the Womersley number
\[
\alpha=\frac{D}{2}\sqrt{\frac{2\pi f}{\nu}},
\]
the mean Reynolds number \(Re_s=U_sD/\nu\), the oscillatory Reynolds number \(Re_o=U_oD/\nu=ARe_s\), and the instantaneous Reynolds number
\[
Re(t)=Re_s[1+A\sin(2\pi f t)].
\]
For moderate pulsation amplitudes, the first instability remains subcritical and produces localized turbulent puffs analogous to those in steady pipe flow [1709.03738].

The onset was mapped through puff survival rather than immediate breakdown. Using the instantaneous escape rate \(\kappa(Re)=1/\tau(Re)\), with steady-flow puff lifetime
\[
\tau(Re)=\exp\{\exp[a\,Re+b]\},
\]
the survival probability over one cycle is
\[
P(Re_s,A,\alpha)=\exp\!\left[-\int_0^T \kappa\!\left(Re_s[1+A\sin(2\pi t/T)]\right)\,dt\right],
\]
and the threshold \(Re_{s,c}(\alpha;A)\) is defined by \(P=0.5\) [1709.03738].

This construction yields three regimes. For \(\alpha \lesssim 2\!-\!3\), the flow is quasi-steady over each cycle, turbulence decays during the slow phase, and transition is significantly delayed. In the \(\alpha \to 0\) limit,
\[
Re_{s,c} \simeq \frac{Re_{c,0}}{1-A},
\]
because the minimum instantaneous Reynolds number must stay above the steady critical value. For \(A=0.4\) and \(Re_{c,0}=2040\), this gives \(Re_{s,c}\approx 3400\). For \(\alpha \gtrsim 10\!-\!12\), turbulence cannot adjust to the instantaneous \(Re(t)\), puff lifetimes revert to their steady-flow values at \(Re_s\), and the threshold becomes essentially unaffected by pulsation, \(Re_{s,c}\simeq Re_{c,0}\approx 1860\!-\!2040\), depending on the survival criterion. Between these limits, roughly \(2.5\lesssim \alpha \lesssim 12\), the threshold drops steeply from the quasi-steady asymptote toward the steady-pipe value [1709.03738].

A different threshold proposal, due to Trinh, uses the normalized viscous-layer thickness \(\delta_\nu^+=\delta_\nu u_\tau/\nu\). In that construction, the asymptotic value \(\delta_{\nu,\infty}^+ \simeq 65\!-\!67\) is equated with the normalized pipe radius \(R^+=Ru_\tau/\nu\), leading to
\[
Re_c=\frac{(\delta_{\nu,\infty}^+)^2}{2}\approx 2.1\times 10^3
\]
for Newtonian pipe flow [1007.0810]. This is a source-specific criterion rather than a consensus one, but it is representative of attempts to express transition in terms of near-wall elemental dynamics.

## 3. Geometry, perturbation structure, and constitutive effects

Outside the straight-pipe setting, transition flow is strongly shaped by geometry, perturbation form, curvature, interfacial physics, and constitutive response. The studies below show that the transitional state is often not fixed by a single global Reynolds number, but by a coupled structure of local instability, nonlinear interaction, and finite-amplitude pathway selection [1807.10665], [1007.0081], [1901.08931], [2508.21362], [1005.1380].

| System | Transition sequence or threshold | Dominant mechanism |
|---|---|---|
| Sudden-expansion pipe | \(A_c \propto Re_t^{-3}\); LS, US2, US1; hysteresis | Symmetry-broken recirculation and convective shear instability |
| Square duct | \(\alpha=0\) optimals fail nonlinearly; \(\alpha=1,2\) sub-optimals trigger | Quadratic generation of streamwise-mean distortion and edge dynamics |
| Wide-gap Taylor–Couette | ATVF for \(68\le Re\le425\); transitional non-axisymmetry for \(425<Re<575\) | Outer-wall separation-region viscous-layer thickening |
| Liquid–liquid TCP flow | SF \(\to\) DSF \(\to\) transitional droplet flow \(\to\) UBF \(\to\) SBF | Interplay of \(Re\), \(We\), and droplet inertia \(St\) |
| Wavy viscoelastic channel | Convective \(\to\) transition \(\to\) elastic regimes | Competition between convective and elastic forcing |

In the sudden-expansion pipe with expansion ratio \(1:2\), the minimum inlet-vortex amplitude required for sustained unsteady flow obeys the power law \(A_c \propto Re_t^{-3}\), with least-squares fit \(A_c \simeq k Re_t^{-3}\), \(k \simeq 5.3\times10^8\). A weak inlet vortex first generates a weak unsteady pattern in the recirculation bubble, and a new instability, denoted US2, emerges from the high-shear region produced by symmetry-broken recirculation. The disturbance peak convects downstream according to \(z \simeq V_0 t_p + {\rm const}\) with \(V_0 \simeq 0.67\bar U\). The same system exhibits hysteresis: at \(A=0.2\), one reported case gives \(Re_t^u \simeq 1400\), \(Re_t^l \simeq 1325\), \(\Delta Re_h=75\), and \(H \approx 28\%\) [1807.10665].

In square-duct flow, linear optimal perturbations take the form of longitudinal vortices, and at \(Re_\tau=150\) the global optimum at \(\alpha=0\) reaches \(G_{\max}\approx873\) at \(T\approx1.3\). Yet these optimals decay in nonlinear DNS even for initial energy up to \(E_0=10^{-1}\). By contrast, streamwise-modulated sub-optimal disturbances at \(\alpha=1\) and \(\alpha=2\) trigger transition for \(E_0\) of order \(10^{-3}\), because quadratic interactions rapidly generate a streamwise-homogeneous distorted mean mode that becomes unstable. At the laminar–turbulent boundary the edge state is a self-sustained orbit with alternating pairs of large-scale vortices, and the reported threshold is \(Re_\tau\approx150\) with minimal perturbation energy \(E_0^\star\approx7.8\times10^{-3}\) for \(\alpha=1\) [1007.0081].

In wide-gap Taylor–Couette flow at radius ratio \(\eta=0.5\), the axisymmetric Taylor-vortex flow occupies \(68\le Re\le425\), a transitional non-axisymmetric regime appears for \(425<Re<575\), and wavy-vortex flow emerges beyond \(Re\approx575\). The proposed source of non-axisymmetric disturbance is a sudden increase in viscous-layer thickness in the outer-wall separation region. This localized thickening produces periodic secondary flow, raises the natural wavelength from \(NNW\approx2.00\) to \(NNW\approx2.67\) in the transition window, and reduces the dimensionless torque by as much as \(11\%\) relative to axisymmetric Taylor-vortex flow [1901.08931].

In liquid–liquid Taylor–Couette–Poiseuille flow with axial through-flow, the experimentally observed cascade is stratified flow, disturbed stratified flow, transitional droplet flow, unstable banded flow, and stable banded flow. The two principal transition lines in the \(Re\)–\(We\) plane obey power laws \(We_I=C_I Re^{A_I}\) and \(We_{II}=C_{II}Re^{A_{II}}\). For TR I, \(A_I\) remains in the range \(0.40\!-\!0.57\) across glycerol concentration, while for TR II, \(A_{II}\) decreases from \(0.57\) to \(0.03\) with increasing viscosity. Stable banded flow is linked to loss of droplet inertia, and the UBF \(\to\) SBF transition coincides with a critical \(St \simeq 1\) in the high-viscosity case [2508.21362].

In a two-dimensional wavy-walled viscoelastic channel governed by the FENE-P model, Yatou identified three steady regimes—convective, transition, and elastic—in the \(Re\)–\(Wi\) plane. The first transition occurs when streamwise elastic forcing becomes comparable to convection. During this transition a separation vortex disappears, a jet induced by viscoelasticity moves into the bulk, viscous wall friction drops abruptly, and elastic wall friction rises sharply. Near the first transition boundary the stress scales as \(\tau \sim Re\), and the critical line obeys \(Wi_{c1}-A_3 \propto (Re-Re_c)^{1/2}\) [1005.1380].

## 4. Coherent structures, front dynamics, and reduced diagnostics

Transition flow is often diagnosed through coherent-structure analysis rather than solely through bulk statistics. In hypersonic axisymmetrical compression-ramp flow at Mach \(5\), the transition scenario was decomposed using Spectral Proper Orthogonal Decomposition and resolvent analysis about the mean flow. The fluctuation field satisfies
\[
\partial_t q' = J q' + f',
\qquad
H(i\omega,m)=(i\omega I-J_m)^{-1},
\]
and the SPOD and resolvent modes were compared across the attached boundary layer, mixing layer, and reattachment region. The observed transition starts with the linear amplification of oblique first modes over a broad frequency band, followed by nonlinear interaction of oblique modes into streamwise streaks, and then linear amplification and breakdown of those streaks near reattachment. At \((m=72,f=51\,{\rm kHz})\), the leading SPOD mode captures approximately \(92\%\) of the energy, while the corresponding resolvent and SPOD structures align with \(|\langle \Psi,\psi\rangle| \approx 0.98\). The paper explicitly concludes that early nonlinear interaction is essential to the transition process [2009.08359].

A very different but equally structured example appears in quasi-two-dimensional bounded granular heap flow after a step change in feed rate. There, the transition between steady states is modeled through local mass conservation,
\[
\partial_t h + \partial_x q = 0,
\]
and a linear flux–slope law,
\[
q(x,t) = -A\,\partial_x h + B.
\]
In deviation variables, the problem reduces to a diffusion equation,
\[
\partial_t \tilde q = A\,\partial_{xx}\tilde q,
\]
together with a moving boundary at the wedge front. The front position obeys
\[
x_f(t)=\sqrt{Ct},
\]
the front speed decays as \(v_f(t)\propto t^{-1/2}\), and a second, slower transient follows once the front reaches the endwall. The transition flow is therefore represented as a diffusion-controlled moving-boundary problem rather than an instability cascade [1708.07270].

These two cases use different mathematics—one spectral and non-normal, the other diffusive and kinematic—but both treat the transition flow as a finite-time organization of localized structures: oblique modes and streaks in the hypersonic boundary-layer interaction, and a propagating wedge in the granular heap [2009.08359], [1708.07270].

## 5. Pattern-forming, topological, and oscillatory transition flows

Transition flow also denotes regime change in systems where the relevant order parameter is not turbulence intensity alone. In thermal vibrational convection, direct numerical simulation in a square box revealed three regimes: periodic circulation, columnar, and columnar-broken. With dimensionless vibration amplitude \(a\), angular frequency \(\omega\), and fixed \(Pr=4.38\), the critical control parameter for the first transition is the vibrational Rayleigh number
\[
Ra_{\rm vib}=\frac{a^2\omega^2 Pr}{2}.
\]
The onset of the columnar regime is reported at \(Ra_{\rm vib}^* \simeq 2.19\times10^6\). Below this threshold, the flow is an oscillatory circulation; above it, thermal plumes merge into nearly stationary columns and mid-height velocity fluctuations decrease; at still higher forcing the columns break and the flow reorganizes into a large-scale circulation with a sudden increase in fluctuations [2303.16752].

In quasi-geostrophic channel flow, the phrase “transition flow” refers to time-periodic solutions born when a zonal jet loses stability through a Hopf bifurcation. The normal form on the center manifold is
\[
\frac{dz}{dt}=\overline{\beta(R)}\,z+b\,z|z|^2+o(|z|^3),
\]
and the real part of the cubic coefficient,
\[
\gamma=\Re b,
\]
determines whether the transition is continuous (Type I, \(\gamma<0\)) or catastrophic (Type II, \(\gamma>0\)). Numerical evaluation for \(0.01\le E\le0.03\) and \(0.1\le a\le0.6\) yielded \(\gamma>0\), typically \(0.13\!-\!0.20\), suggesting that catastrophic transition is preferred in the parameter region examined [1502.04974].

In two-dimensional passive nematics, flow coupling changes the topology of the defect-mediated transition itself. Without flow alignment, the transition is consistent with the Berezinskii–Kosterlitz–Thouless scenario, with the standard defect-pair free-energy balance leading to \(\hat T_c^Q \approx \pi/8\). With hydrodynamics and \(\lambda=0\), the defect-creation threshold is approximately \(T_{\rm dc}^*\approx0.042\); without flow coupling it is approximately \(0.13\). For strain-rate-aligning nematics with \(\lambda=1\), bend–splay walls appear already at \(T^*\approx0.004\), the defect-creation threshold drops to \(T_{\rm dc}^*\approx0.004\), and once defects are created they remain unbound across the entire fluctuation range in both forward and backward protocols. In this sense, flow alignment suppresses the reversible binding–unbinding mechanism characteristic of BKT [2510.08282].

In confluent epithelial monolayers, the flow generated by a T1 transition is a localized saddle. After rotating and averaging many events, Jain et al. obtained an ensemble-averaged velocity field with
\[
v_x \simeq +\lambda x,\qquad v_y \simeq -\lambda y
\]
near the four-cell vertex, together with quadrupolar vorticity around the rearranging cells. The corresponding tissue fluidization is quantified through the mean separation \(\Pi(t)\) of initially neighboring cells, which satisfies
\[
\Pi(t)=\Pi_0+\alpha N_{T1}(t),
\]
with \(\alpha\sim0.004\) cell-diameters per T1 for the chosen parameters, and the same linear law remains nearly unchanged under variations in deformability, rotational noise, adhesion, and active fraction [2403.20100].

## 6. Transition flow as a learned map in generative modeling

A conceptually distinct use of the term appears in diffusion- and flow-based generative modeling. In “Transition Flow Matching,” the transition flow is defined on a standard flow-matching coupling \(X_0\sim p_0\), \(X_1\sim p_1\), with interpolation
\[
X_t=\alpha(t)X_0+\beta(t)X_1,\qquad 0\le t\le1.
\]
Given a later time \(r\ge t\), the conditional transition state is
\[
x_{t\to r}^Z=x_t+\int_t^r v(x_\tau,\tau\mid Z)\,d\tau,
\]
and the transition flow is the conditional expectation
\[
T(x_t,t,r)\equiv \mathbb E[X_{t\to r}^Z\mid X_t=x_t].
\]
By definition, \(T(x_t,t,r)=x_{t\to r}\), and the associated average velocity is
\[
u(x_t,t,r)=\frac{T(x_t,t,r)-x_t}{r-t}.
\]
This makes transition flow a global transport quantity rather than a local velocity field [2603.15689].

Training uses a neural approximation \(X^\theta(x_t,t,r)\) and a conditional loss based on the computable linear-interpolation target \(x_{t\to r}^Z=(1-r)X_0+rX_1\). The total derivative is obtained by a Jacobian–vector product,
\[
\frac{d}{dt}X^\theta=\partial_{x_t}X^\theta\cdot v(X_t,t\mid Z)+\partial_t X^\theta,
\]
and the standard Transition Flow Matching objective is
\[
\mathcal L_{\rm TFM}(\theta)=\mathbb E\left\|
\left[(1-r)X_0+rX_1+(r-t)\frac{d}{dt}X^\theta(X_t,t,r)\right]
-X^\theta(X_t,t,r)
\right\|^2.
\]
When \(r=1\), the formulation reduces to mean-velocity flow, since
\[
T(x_t,t,1)=x_t+(1-t)\bar v(x_t,t).
\]
The method supports one-step sampling \(x_1=X^\theta(x_0,0,1)\) and few-step sampling on arbitrary time grids [2603.15689].

The empirical results reported in the paper place this abstract notion of transition flow in direct algorithmic use. On CIFAR-10 with a \(\sim55\)M-parameter UNet, the reported FIDs are \(2.77\) for \(NFE=1\), \(2.08\) for \(NFE=2\), \(1.96\) for \(NFE=5\), and \(1.91\) for \(NFE=10\). On latent ImageNet \(256\times256\) with a \(\sim676\)M-parameter transformer, the one-step FID is \(3.02\) and the two-step FID is \(2.77\) [2603.15689].

Taken together, these literatures suggest that transition flow is not a single mechanism but a class of structured intermediate dynamics. A plausible implication is that the most useful descriptions are rarely purely global: they are localized in space, scale, or state space, and they emphasize puffs, bands, streaks, wedges, columns, defect walls, or explicit maps \(T(x_t,t,r)\) as the natural objects that organize the transition.

Source: https://www.emergentmind.com/topics/transition-flow