In the majority of applications steady flows are better than unsteady flows. Steady flows are usually associated with smaller fuel consumption, less fatigue, and less noise. Theoretically, a steady flow always exists as a solution of the governing equations. However, in engineering applications fluid flows are usually unsteady, or even turbulent, because the corresponding steady flow is unstable, that is, if disturbed it will never return to the steady state. In the bulk of the work on flow stability the flow velocity is represented as a sum of the steady solution and a small perturbation, and the non-linear terms are neglected. The resulting linear problem is much easier to solve. However, steady flows are often stable with respect to infinitesimally small perturbations but unstable with respect to finite perturbations. Moreover, the finite amplitude required to destabilise the flow is often small. Hence, it is the stability with respect to finite disturbances, which represents major practical interest.
Stability of a dynamical system with respect to finite perturbations can be established by finding a Lyapunov function. However, there is no general systematic method for constructing Lyapunov functions -- the discovery of such a function is dependent on the ingenuity and creativity of the investigator. Fortunately, in 2000 a breakthrough in control theory, made in a PhD thesis by P.Parrilo, has provided a constructive method for generating Lyapunov functions for systems whose dynamics can be described by polynomial functions. This method is based on sum-of-squares (SOS) optimization, which reduces the problem of finding a Lyapunov function for a polynomial system to one of constructing a polynomial function that satisfies a selection of algebraic conditions. Using a number of key results from semialgebraic geometry, in particular the Positivestellensatz theorem, the resulting problem can be reformulated as an optimization problem in the form of a Semi-Definite Programme (SDP). This is particularly promising since SDP optimization problems are convex and tractable (e.g. solvable in a number of operations that is a polynomial function of the problem size). A variety of well-supported software codes for solving such problems are freely available.
It is well known that the Galerkin approximation of the incompressible Navier-Stokes equations leads to an ODE with quadratic non-linearity. The SOS approach can be applied to such a system since the dynamics are represented by polynomials. SOS was applied to a hydrodynamic-type ODE system of the 9th order. A method of exploiting in the SOS approach the energy-invariance property of the bilinear component of the systems of hydrodynamic type was developed. This allowed the construction of a Lyapunov function for the value of the Reynolds number about seven times larger than the largest value for which non-linear stability of this system can be proven by the standard energy stability approach. This demonstrates the strong potential of applying SOS in fluid dynamics.
Accurate approximation of a fluid flow requires considering systems with a large number of Galerkin modes. However, application of conventional SOS methods to very large systems is problematic, since the SDP problems to be solved in such cases are prohibitively large. A method of reducing the order of the system in SOS studies, based on representing the higher modes with their energy only, has been developed. The method is guaranteed always to give results at least as good as the standard energy stability analysis. Moreover, it has also been proven mathematically that if the flow remains globally stable for the Reynolds number range beyond that given by the energy stability analysis, then there exists a polynomial Lyapunov function giving a better stability bound than the energy stability theory.
More information:
Try it, it is easy and it is fun!
Oberwolfach Workshop 2431 - Polynomial Optimization for Nonlinear Dynamics: Theory, Algorithms, and Applicationsat the Mathematisches Forschungsinstitut Oberwolfach, Germany 28 July - 2 August 2024.
Studying fluid flows with auxiliary functions and LMIsat the IFAC World Congress, held in Yokohama, Japan on 8-14th July 2023.
Bounding time averages: a road to solving the problem of turbulenceat Institut de Mathématiques de Bordeaux, Bordeaux, May 4, 2023.
Bounding time averagesand
How quasi-steady is the modulation of near-wall turbulence by large-scale structures?(with Yunjiu Yang).
Auxiliary functionals: a path to solving the problem of turbulenceat The Seminar in the Analysis and Methods of PDE (SIAM PDE) on March 4, 2021. Links to the abstract and the video.
Accelerating time averagingat 73rd Annual Meeting of the APS Division of Fluid Dynamics, November 22, 2020: abstract and video.
Accelerating time averaging using auxiliary functionsat the Aerodynamics and Flight Mechanics group seminar, University of Southampton, on 6 February 2019
Coherent structures in wall-bounded turbulence: new directions in a classic problem, London, August 29-31, 2018, with a talk
Large-scale motions for the QSQH theory(with Chi Zhang).
Questions concerning quasi-steady mechanism of the Reynolds number, pressure gradient, and geometry effect on drag reductionat the Workshop on Active Drag Reduction, Aachen, Germany, 15-16 March 2018.
The problem of turbulence: bounding solutions to equations of fluid mechanics & other dynamical systems, with Giovanni Fantuzzi providing exercise sessions, at The 6th Bremen Winter School
Dynamical systems and turbulence, March 12-16, 2018.
Sergei Chernyshenko