Improved global stability bounds for two-dimensional plane Poiseuille flow

TL;DR

This paper improves the known bounds on the nonlinear stability limit of 2D plane Poiseuille flow by constructing Lyapunov functionals and using semidefinite programming, achieving a 22% higher stability Reynolds number.

Contribution

The authors develop a computational method combining mode decomposition and Lyapunov functionals to establish tighter nonlinear stability bounds for plane Poiseuille flow.

Findings

The stability limit exceeds the classical energy bound at certain streamwise lengths.

The flow is globally stable up to Re ≈ 106.8, a 22% increase over previous bounds.

A simple five-mode set suffices to produce improved stability bounds.

Abstract

This work provides new lower bounds on the global (nonlinear) stability limit of pressure-driven two-dimensional plane Poiseuille flow, improving on the energy stability limit, $Re_E$, originally computed by Orr in 1907. Using a computer we carefully construct quartic Lyapunov functionals of the velocity perturbations about the laminar profile, which certify the nonlinear stability of the flow to arbitrary perturbations. The formulation combines a decomposition of the velocity into finitely many energy eigenmodes, referred to as a 'mode set', and an infinite-dimensional 'tail', together with explicit bounds that recast the Lyapunov inequality conditions as semidefinite programs, whose feasibility is tested. Over the streamwise lengths considered, the certified stability limit exceeds the classical energy bound. In particular, at the critical energy-stable streamwise length, where $Re_E\approx 87.59$, the flow is found to be globally stable up to $Re \approx 106.8$ (representing a $22\%$ improvement). Various modestly-sized mode sets, capable of capturing sufficient features of the nonlinear dynamics of energy growth and subsequent decay, are proposed and found to be successful in producing improved bounds, with the simplest one involving only five modes.