Couette-Poiseuille Flow: Analysis & Stability

• Couette-Poiseuille flow is a parallel shear flow characterized by a superposition of wall-driven and pressure-driven forces, resulting in a linear-plus-parabolic velocity profile.
• It plays a central role in energy stability analyses and the determination of critical Reynolds numbers, with various formulations adapting the base flow for different experimental and theoretical settings.
• Research applications extend to studying turbulent transitions, roll-streak interactions, and extensions to non-Newtonian, rotating, or multi-phase flow contexts.

Couette-Poiseuille flow denotes the family of parallel shear flows between two plates in which wall-driven shear and pressure-driven or body-force-driven transport are superposed. In its classical incompressible Newtonian form, the streamwise velocity is a linear-plus-parabolic profile, so the flow interpolates continuously between plane Couette flow and plane Poiseuille flow. Across the literature, the same structure appears in energy-stability theory, laboratory channels with zero net flux, rotating shearing boxes, and rigorous steady Navier-Stokes analysis, with changes only in scaling, parametrization, and boundary data (Lam, 2012, Ghosh et al., 2021, Jiang et al., 2022, Klotz et al., 2017).

1. Canonical formulations of the base flow

The defining feature of Couette-Poiseuille flow is a unidirectional basic state

$\mathbf{U}=(U,0,0),$

with $U$ quadratic in the wall-normal coordinate after nondimensionalization. Several parameterizations are standard, each adapted to a different analytical setting.

Setting	Base profile	Distinguished parameters
Energy-stability family	$U(z)=A(1-z^2)+z$ for $0\le A\le \tfrac12$, with an equivalent form for $\tfrac12<A\le1$	$A=0$ Couette, $A=1$ Poiseuille
Steady 2D channel flow	$u_b(y)=\bigl(U_c\,y/H+U_p[1-(y/H)^2],0\bigr)$	$U_c$ wall speed, $U_p$ pressure amplitude
Rotating shearing box	$U$0	$U$1
Zero-flux laboratory profile	$U$2	$U$3

In the Reynolds-Orr formulation, the Couette-Poiseuille family is written for $U$4 as

$U$5

with $U$6 yielding plane Couette flow $U$7 and $U$8 yielding plane Poiseuille flow $U$9 (Lam, 2012). In steady 2D Navier-Stokes theory, the same family is written as

$U(z)=A(1-z^2)+z$0

so that $U(z)=A(1-z^2)+z$1 gives pure Couette flow and $U(z)=A(1-z^2)+z$2 pure Poiseuille flow (Jiang et al., 2022).

In forced rotating shear flows, the steady $U(z)=A(1-z^2)+z$3-momentum balance produces

$U(z)=A(1-z^2)+z$4

which becomes the dimensionless Couette-Poiseuille profile

$U(z)=A(1-z^2)+z$5

This form is used to model local accretion-disk shear modified by residual forcing (Ghosh et al., 2021). In experiments designed to suppress mean advection, the pressure gradient is adjusted so that the net flux vanishes, producing a zero-mean profile with a low-shear negative Poiseuille region and a high-shear positive Couette region (Klotz et al., 2017).

A recurrent source of confusion is that these formulas are not competing definitions. They are equivalent realizations of the same linear-plus-parabolic structure under different scalings, wall speeds, and flux constraints.

2. Energy stability, modal reduction, and critical Reynolds numbers

A central formulation of Couette-Poiseuille stability is the Reynolds-Orr energy equation for a disturbance velocity $U(z)=A(1-z^2)+z$6. Its kinetic energy

$U(z)=A(1-z^2)+z$7

satisfies the exact identity

$U(z)=A(1-z^2)+z$8

where $U(z)=A(1-z^2)+z$9 is the strain-rate tensor of the basic flow and $0\le A\le \tfrac12$0 (Lam, 2012). For normal modes proportional to $0\le A\le \tfrac12$1, the Euler-Lagrange equations reduce to a coupled system for $0\le A\le \tfrac12$2 with

$0\le A\le \tfrac12$3

The minimum Reynolds number for global energy decay is obtained as an unconstrained minimization over divergence-free disturbances,

$0\le A\le \tfrac12$4

Numerically, the normal-mode problem is converted from a fourth- or sixth-order ODE system into a Fredholm integral equation of the second kind and then discretized with a Gauss-Legendre Nystrom rule (Lam, 2012).

The resulting thresholds distinguish sharply between two-dimensional and fully three-dimensional disturbances. For streamwise disturbances $0\le A\le \tfrac12$5, the reported minima are $0\le A\le \tfrac12$6 for plane Couette flow, $0\le A\le \tfrac12$7 for a Couette-Poiseuille profile at $0\le A\le \tfrac12$8, and $0\le A\le \tfrac12$9 for plane Poiseuille flow. For spanwise disturbances $\tfrac12<A\le1$0, the corresponding values are $\tfrac12<A\le1$1, $\tfrac12<A\le1$2 at $\tfrac12<A\le1$3, and $\tfrac12<A\le1$4. For fully three-dimensional disturbances, the minima become $\tfrac12<A\le1$5 at $\tfrac12<A\le1$6, $\tfrac12<A\le1$7 at $\tfrac12<A\le1$8, and $\tfrac12<A\le1$9 at $A=0$0, with the Poiseuille minimum attained near $A=0$1 (Lam, 2012).

The principal conclusion is that, except for plane Couette flow, the smallest energy-stability threshold is attained by genuinely three-dimensional disturbances. This establishes that in nonlinear energy stability, three-dimensional modes are more dangerous than purely two-dimensional ones for the Couette-Poiseuille family, whereas the linear-profile limit $A=0$2 preserves the special role of two-dimensional spanwise modes (Lam, 2012).

A distinct but related linear-stability setting is plane Couette-Poiseuille flow with uniform cross-flow. There the modified Orr-Sommerfeld equation contains a constant $A=0$3 term, so Squire’s transformation still applies and the first instability is two-dimensional. For each wall-velocity parameter $A=0$4, two unconditional-stability intervals in the cross-flow Reynolds number $A=0$5 are reported. The lower cutoff is associated with profile skewing, critical-layer displacement, and reduced energy production; an intermediate regime exhibits long-wave or resonant Tollmien-Schlichting instability; and at very large $A=0$6 the flow restabilizes as energy production decays like $A=0$7 (Guha et al., 2010).

3. Convective, absolute, and rotation-induced instabilities

For spatio-temporal instability, plane Couette-Poiseuille flow with one moving plate is often written as

$A=0$8

where $A=0$9 is the plate-speed ratio and $A=1$0. Two-dimensional perturbations

$A=1$1

satisfy the Orr-Sommerfeld equation

$A=1$2

with $A=1$3 (Srinivas et al., 2022).

Using the Briggs-Bers criterion, the flow exhibits a sequence

$A=1$4

where $A=1$5 is downstream convective instability, $A=1$6 absolute instability, and $A=1$7 upstream convective instability. For $A=1$8, the flow remains convectively unstable downstream. For negative $A=1$9, increasing reverse plate motion first opens an absolute-instability window $u_b(y)=\bigl(U_c\,y/H+U_p[1-(y/H)^2],0\bigr)$0, and then produces an upstream convective regime in which the unstable wave packet travels against the bulk flow (Srinivas et al., 2022).

This regime is notable because the flow is non-inflectional. The reported result identifies Couette-Poiseuille flow as the first instance of a non-inflectional absolute instability within constant-viscosity formulation. High-$u_b(y)=\bigl(U_c\,y/H+U_p[1-(y/H)^2],0\bigr)$1 analysis shows that for plane Poiseuille flow both the leading- and trailing-edge velocities approach zero as $u_b(y)=\bigl(U_c\,y/H+U_p[1-(y/H)^2],0\bigr)$2; consequently, even slight reverse plate motion can trigger absolute instability and then upstream convection. The same work interprets the instability pocket in a complex Ginzburg-Landau framework and attributes to viscosity a dual role: at moderate $u_b(y)=\bigl(U_c\,y/H+U_p[1-(y/H)^2],0\bigr)$3 dispersion is large enough to sustain a wide absolute-instability window, whereas at very low or very high $u_b(y)=\bigl(U_c\,y/H+U_p[1-(y/H)^2],0\bigr)$4 the window shrinks (Srinivas et al., 2022).

Rotation changes the picture again. In a rotating shearing box with extra forcing $u_b(y)=\bigl(U_c\,y/H+U_p[1-(y/H)^2],0\bigr)$5, the background linear shear becomes the Couette-Poiseuille profile

$u_b(y)=\bigl(U_c\,y/H+U_p[1-(y/H)^2],0\bigr)$6

and perturbations satisfy coupled Orr-Sommerfeld and Squire equations with Coriolis terms. For Keplerian rotation $u_b(y)=\bigl(U_c\,y/H+U_p[1-(y/H)^2],0\bigr)$7, pure vertical perturbations $u_b(y)=\bigl(U_c\,y/H+U_p[1-(y/H)^2],0\bigr)$8 become unstable at $u_b(y)=\bigl(U_c\,y/H+U_p[1-(y/H)^2],0\bigr)$9, while fully three-dimensional perturbations $U_c$0 become unstable at $U_c$1. For vertical modes, the reported critical values are $U_c$2 at $U_c$3 for $U_c$4, and $U_c$5 at $U_c$6 for $U_c$7 (Ghosh et al., 2021).

The contrast between Couette and Poiseuille curvature is explicit in that setting. Plane Couette flow with rotation is stable for $U_c$8, with eigenvalues corresponding to epicyclic oscillations plus viscous damping, whereas plane Poiseuille flow develops a growing branch under rotation. This distinction is used to argue that small forcing in a local accretion disk can convert linearly stable shear into a Poiseuille-type profile that becomes linearly unstable even for Keplerian rotation (Ghosh et al., 2021).

These statements do not conflict with the energy-stability thresholds discussed above. They concern different notions of stability: sufficient conditions for global energy decay, temporal or spatio-temporal modal growth, and rotating versus nonrotating base states.

4. Subcritical transition, turbulent spots, and the roll-streak cycle

Laboratory work has used Couette-Poiseuille flow to isolate subcritical transition in a configuration with nearly zero mean advection velocity. In the reported apparatus, the moving belt drives a Couette component while the induced pressure gradient generates a counterflowing Poiseuille component, so turbulent structures are nearly stationary in the laboratory frame. For the zero-mean profile, natural transition remains globally laminar for $U_c$9, becomes intermittent with localized spots for $U_p$0, produces oblique turbulent bands above $U_p$1, and becomes featureless turbulent for $U_p$2 (Klotz et al., 2017).

Under permanent finite-amplitude forcing by a sphere attached near the belt, a localized turbulent spot appears above $U_p$3. The spot has two nested regions: an active core characterized by wavy streaks similar to traveling waves, and a surrounding region that also contains weak undisturbed streaks and oblique waves at the laminar-turbulent interface. Both streamwise and spanwise extents increase with $U_p$4, and the centroid of the spot crosses from the Poiseuille-dominated side to the Couette-dominated side near $U_p$5, indicating that sustained activity shifts toward the stronger shear near the moving wall as $U_p$6 increases. Spatio-temporal filtering further shows that the wavy-streak core travels in the downstream Poiseuille direction, even though the spot as a whole is nearly stationary (Klotz et al., 2017).

Quench experiments examine the reverse process, namely decay after abrupt reduction of $U_p$7. In a large-aspect-ratio plane Couette-Poiseuille channel, the spanwise velocity component $U_p$8, which contains the roll signature, decays faster than the streamwise streak component $U_p$9. The turbulent fraction based on thresholding $U$00,

$U$01

decays linearly across the full range of final Reynolds numbers studied, while the corresponding spanwise energy

$U$02

decays exponentially. The slope $U$03 extrapolates to zero at $U$04, and the decay rate $U$05 extrapolates to zero at $U$06, both close to the independently determined self-sustainment threshold $U$07 (Liu et al., 2020).

The same experiments identify a two-stage decay of streaks at $U$08: a slow stage while rolls remain present and regenerate streaks by lift-up, followed by a faster decay once the rolls have vanished. This directly links the decay dynamics to the roll-streak feedback cycle rather than to a single viscous time scale (Liu et al., 2020).

Direct numerical simulations near transition make that feedback more explicit. In a periodic Couette-Poiseuille domain with $U$09, the optimal transient growth for infinitesimal disturbances is reported as

$U$10

with $U$11 and $U$12. The flow variables are organized into streak amplitude $U$13, roll amplitude $U$14, and waviness amplitude $U$15, and the Waleffe feedback closure is written as

$U$16

After a short transient, $U$17 becomes negligible in both sustained turbulence and large-amplitude decay, giving a quadratic roll-waviness law with $U$18, $U$19, and

$U$20

for sufficiently large waviness $U$21 (Etchevest et al., 5 Aug 2025).

That DNS study distinguishes three regimes: laminar decay, nonlinear decay, and turbulent steady state. Sustained turbulence is reported for $U$22 and sufficiently large initial waviness, whereas a quench analysis in the appendix places the critical Reynolds number for sustained turbulence between $U$23 and $U$24 (Etchevest et al., 5 Aug 2025). Taken together with the experiments, this suggests that reported transition thresholds are strongly geometry- and protocol-dependent, while the central roll-streak-waviness mechanism is robust across configurations.

5. Steady Navier-Stokes theory, uniqueness, and the inviscid limit

A separate line of work treats Couette-Poiseuille flow as a steady boundary-value problem for the 2D Navier-Stokes equations. In a finite channel $U$25, the steady system

$U$26

admits the parallel basic flow

$U$27

with no-slip boundary conditions $U$28 and $U$29 (Jiang et al., 2022).

For $U$30, smooth shear profiles $U$31 that are $U$32-close to $U$33 and satisfy $U$34 give rise, for sufficiently small $U$35, to unique nearby $U$36 steady Navier-Stokes solutions stable under infinitesimal perturbations. In the special Couette case $U$37, the stability statement is stronger: under the additional smallness and degeneracy assumptions stated in the paper, one can expand to arbitrarily high order and prove existence of a unique nearby steady solution for perturbations of $U$38 that are $U$39 in size but $U$40 as $U$41. The same framework extends to controlled forcing: for any smooth monotone shear $U$42 with $U$43 and bounded $U$44-norm, a small distributed forcing can be chosen so that the steady Navier-Stokes system has a unique solution in an $U$45-neighborhood of $U$46 (Jiang et al., 2022).

The same paper derives a zero-viscosity limit. From the uniform-in-$U$47 bounds, the steady Navier-Stokes solutions converge strongly to the corresponding Euler shear flow: $U$48 with

$U$49

In the Couette-only case, convergence holds in $U$50 for any $U$51, and also in $U$52 if $U$53 is slightly smoother (Jiang et al., 2022).

A complementary result addresses uniqueness and structural stability in an infinite straight channel $U$54 with prescribed wall velocities and arbitrary flux. The Couette-Poiseuille base flow is written as

$U$55

with pressure $U$56, $U$57, and flux

$U$58

Under the no-flow-reversal assumption

$U$59

the linearized steady operator is continuously invertible in both global and local Sobolev settings, yielding local uniqueness of the base solution and nonlinear structural stability under small external forces (Galdi et al., 18 Feb 2026).

Within a natural symmetric class, global uniqueness is also obtained. When $U$60 and $U$61, the only symmetric solution of the full steady problem with zero forcing, or with sufficiently small forcing, is the base flow itself. A stated corollary is that for any flux $U$62, the Poiseuille profile

$U$63

is the unique solution in the symmetric class, with no smallness restriction on $U$64 (Galdi et al., 18 Feb 2026).

The same analysis also shows the importance of excluding flow reversal. If $U$65 changes sign in the interior, the linearized operator can lose injectivity; the paper constructs such a counterexample by combining periodicity in $U$66 with Orr-Sommerfeld spectral-instability results (Galdi et al., 18 Feb 2026).

6. Extensions to suspensions, rarefied gases, and non-Newtonian constitutive laws

Couette-Poiseuille configurations remain analytically useful when the constitutive model is changed. In concentrated suspensions, a two-phase formulation introduces a local solid volume fraction $U$67, separate velocities $U$68 and $U$69, fluid pressure $U$70, and interphase drag

$U$71

For plane Couette flow, the base state has $U$72 and uniform $U$73. For pressure-driven Poiseuille flow, the base state can develop a jammed central plug $U$74 in which $U$75 and $U$76, with a sheared region outside. The model is well posed near maximum packing only if the friction parameter satisfies $U$77, equivalently $U$78. If $U$79, a collision-pressure branch becomes ill posed; when $U$80, that branch is stable but a convection-induced transient growth mechanism remains (Ahnert et al., 2018).

The plug region in the suspension model has a Bingham-type structure, and the paper explicitly compares it with single-phase Bingham flow. The conclusion is not that the two models are equivalent, but rather that they share the plug-shear geometry while differing in stability mechanisms: the two-phase model retains both collision-pressure-driven ill-posedness and convection-induced transient growth, neither of which appears in the single-phase Bingham fluid when physically correct boundary conditions are imposed (Ahnert et al., 2018).

In dilute gases, the steady Couette-Poiseuille problem can be treated kinetically rather than hydrodynamically. For a BGK model with plates in relative motion and a uniform body force $U$81, the distribution function satisfies

$U$82

and the solution is expanded in the scaled force $U$83 at arbitrary shear-rate Knudsen number

$U$84

At second order when $U$85, the midplane temperature curvature is

$U$86

Hence the temperature profile has a local minimum at the center when $U$87, corresponding to the bimodal Poiseuille-type state, and becomes Couette-like and parabolic for larger shear (Tij et al., 2010).

That kinetic solution also quantifies departures from Navier-Stokes-Fourier theory. While the shear stress agrees with the Newtonian constitutive prediction at the order considered, kinetic theory produces normal-stress differences, a finite streamwise heat flux $U$88, and a nonzero $U$89 even when $U$90. These are explicitly identified as non-Newtonian effects beyond the NSF description (Tij et al., 2010).

For generalized Newtonian liquids, exact steady Couette-Poiseuille solutions have also been obtained for power-law fluids with partial slip and uniform cross-flow. In nondimensional form the governing relations are

$U$91

with linear Navier slip at both walls. Depending on the pressure-gradient parameter $U$92, cross-flow parameter $U$93, and power-law index $U$94, the resulting flow can be monotonic Couette-type or nonmonotonic Poiseuille-type, and Couette-type profiles may be convex, linear, or concave. Closed-form solutions are given for $U$95, $U$96, and $U$97. The same analysis identifies a breakdown of the pure power-law model for the dilatant case $U$98, where beyond a finite critical cross-flow no real solution for the wall shear parameter exists (El-Mistikawy, 2018).

Across these extensions, Couette-Poiseuille flow functions less as a single formula than as a structural template: a wall-driven and pressure-driven superposition whose stability, transition, and constitutive response are then reshaped by rotation, concentration effects, rarefaction, slip, cross-flow, and non-Newtonian rheology.