Bounds on mean energy in the Kuramoto-Sivashinsky equation computed using semidefinite programming

Abstract

We present methods for bounding infinite-time averages in dynamical systems governed by nonlinear PDEs. The methods rely on auxiliary functionals, which are similar to Lyapunov functionals but satisfy different inequalities. The inequalities are enforced by requiring certain expressions to be sums of squares of polynomials, and the optimal choice of auxiliary functional is posed as a semidefinite program (SDP) that can be solved computationally. To formulate these SDPs we approximate the PDE by truncated systems of ODEs and proceed in one of two ways. The first approach is to compute bounds for the ODE systems, increasing the truncation order until bounds converge numerically. The second approach incorporates the ODE systems with analytical estimates on their deviation from the PDE, thereby using finite truncations to produce bounds for the full PDE. We apply both methods to the Kuramoto-Sivashinsky equation, where we compute upper bounds on the spatiotemporal average of energy by employing polynomial auxiliary functionals up to degree six. The first approach is used for most computations, but a subset of results are checked using the second approach, and the results agree to high precision. These bounds apply to all odd solutions of period $2\pi L$, where $L$ is varied. Sharp bounds are obtained for $L\le 10$, and trends suggest that more expensive computations would yield sharp bounds at larger $L$ also. The bounds are known to be sharp (to within 0.1% numerical error) because they are saturated by the simplest nonzero steady states, which apparently have the largest mean energy among all odd solutions. Prior authors have conjectured that mean energy remains $O(1)$ for $L\gg 1$ since no particular solutions with larger energy have been found. Our bounds constitute the first positive evidence for this conjecture, albeit up to finite $L$, and they offer some guidance for analytical proofs.

32 pages, 6 figures; v2: revisions after review, new proposition added

References in corpus (2)

• Symmetry groups, semidefinite programs, and sums of squares
• Optimal bounds and extremal trajectories for time averages in nonlinear dynamical systems

Cited by in corpus (26)

• Bounding extreme events in nonlinear dynamics using convex optimization
• Concurrent multi-parameter learning demonstrated on the Kuramoto-Sivashinsky equation
• The background method: Theory and computations
• Relationship between the methods of bounding time averages
• On the well-posedness of an anisotropically-reduced two-dimensional Kuramoto-Sivashinsky equation
• Bounds on heat transport for convection driven by internal heating
• A Framework for Input-Output Analysis of Wall-Bounded Shear Flows
• Bounding extrema over global attractors using polynomial optimisation
• Inferring symbolic dynamics of chaotic flows from persistence
• Heat transport bounds for a truncated model of Rayleigh-Bénard convection via polynomial optimization
• Finding extremal periodic orbits with polynomial optimisation, with application to a nine-mode model of shear flow
• Maximum Amplification of Enstrophy in 3D Navier-Stokes Flows
• Bounds for internally heated convection with fixed boundary heat flux
• Optimal cooling of an internally heated disc
• Rigorous bounds on the heat transport of rotating convection with Ekman pumping
• Regularity criteria for the Kuramoto-Sivashinsky equation in dimensions two and three
• Finding unstable periodic orbits: a hybrid approach with polynomial optimization
• A study of the double pendulum using polynomial optimization
• Convex relaxations of integral variational problems: pointwise dual relaxation and sum-of-squares optimization
• Minimum wave speeds in monostable reaction-diffusion equations: sharp bounds by polynomial optimization
• Data-driven Discovery of Invariant Measures
• Remarks on the stabilization of large-scale growth in the 2D Kuramoto-Sivashinsky equation
• A PIE Representation of Scalar Quadratic PDEs and Global Stability Analysis Using SDP
• Bounds on dissipation in three-dimensional planar shear flows: reduction to two-dimensional problems
• Stability of a Parametrically Driven, Coupled Oscillator System: An Auxillary Function Method Approach
• Sum-of-squares bounds on correlation functions in a minimal model of turbulence