No abstract

With the goal of providing the first example of application of a recently proposed method, thus demonstrating its ability to give results in principle, global stability of a version of the rotating Couette flow is examined. The flow depends on the Reynolds number and a parameter characterizing the magnitude of the Coriolis force. By converting the original Navier–Stokes equations to a finite-dimensional uncertain dynamical system using a partial Galerkin expansion, high-degree polynomial Lyapunov functionals were found by sum-of-squares of polynomials optimization. It is demonstrated that the proposed method allows obtaining the exact global stability limit for this flow in a range of values of the parameter characterizing the Coriolis force. Outside this range a lower bound for the global stability limit was obtained, which is still better than the energy stability limit. In the course of the study, several results meaningful in the context of the method used were also obtained. Overall, the results obtained demonstrate the applicability of the recently proposed approach to global stability of the fluid flows. To the best of our knowledge, it is the first case in which global stability of a fluid flow has been proved by a generic method for the value of a Reynolds number greater than that which could be achieved with the energy stability approach.

The sensitivity of long-time averages of a hyperbolic chaotic system to parameter perturbations can be determined using the shadowing direction, the uniformly-bounded-in-time solution of the sensitivity equations. Although its existence is formally guaranteed for certain systems, methods to determine it are hardly available. One practical approach is the Least-Squares Shadowing (LSS) algorithm (Q Wang, SIAM J Numer Anal 52, 156, 2014), whereby the shadowing direction is approximated by the solution of the sensitivity equations with the least square average norm. Here, we present an alternative, potentially simpler shadowing-based algorithm, termed periodic shadowing. The key idea is to obtain a bounded solution of the sensitivity equations by complementing it with periodic boundary conditions in time. We show that this is not only justifiable when the reference trajectory is itself periodic, but also possible and effective for chaotic trajectories. Our error analysis shows that periodic shadowing has the same convergence rates as LSS when the time span T is increased: the sensitivity error first decays as 1{T and then, asymptotically as 1{?T . We demonstrate the approach on the Lorenz equations, and also show that, as T tends to infinity, periodic shadowing sensitivities converge to the same value obtained from long unstable periodic orbits (D Lasagna, SIAM J Appl Dyn Syst 17, 1, 2018) for which there is no shadowing error. Finally, finite-difference approximations of the sensitivity are also examined, and we show that subtle non-hyperbolicity features of the Lorenz system introduce a small, yet systematic, bias.

Abstract

No abstract

No abstract

A rapid predictive tool based on the linearised Reynolds-averaged Navier–Stokes equations is proposed in this work to investigate secondary currents generated by streamwise-independent surface topography modulations in turbulent channel flow. The tool is derived by coupling the Reynolds-averaged momentum equation to the Spalart–Allmaras transport equation for the turbulent eddy viscosity, using a nonlinear constitutive relation for the Reynolds stresses to capture correctly secondary motions. Linearised equations, describing the steady flow response to arbitrary surface modulations, are derived by assuming that surface modulations are shallow. Since the equations are linear, the superposition principle holds and the flow response induced by an arbitrary modulation can be obtained by combining appropriately the elementary responses obtained over sinusoidal modulations at multiple spanwise length scales. The tool permits a rapid exploration of large parameter spaces characterising structured surface topographies previously examined in the literature. Here, channels with sinusoidal walls and with longitudinal rectangular ridges are considered. For sinusoidal walls, a large response is observed at two spanwise wavelengths scaling in inner and outer units respectively, mirroring the amplification mechanisms in turbulent shear flows observed from transient growth analysis. For longitudinal rectangular ridges, the model suggests that the analysis of the response and the interpretation of the topology of secondary structures is facilitated when the ridge width and the gap between ridges are used instead of other combinations proposed in the literature.

In this paper a novel nonlinear feedback control design methodology for incompressible fluid flows aiming at the optimisation of long-time averages of flow quantities is presented. It applies to reduced-order finite-dimensional models of fluid flows, expressed as a set of first-order nonlinear ordinary differential equations with the right-hand side being a polynomial function in the state variables and in the controls. The key idea, first discussed in Chernyshenko et al. (Phil. Trans. R. Soc. Lond. A, vol. 372, 2014, 20130350), is that the difficulties of treating and optimising long-time averages of a cost are relaxed by using the upper/lower bounds of such averages as the objective function. In this setting, control design reduces to finding a feedback controller that optimises the bound, subject to a polynomial inequality constraint involving the cost function, the nonlinear system, the controller itself and a tunable polynomial function. A numerically tractable and efficient approach to the solution of such optimisation problems, based on sum-of-squares techniques and semidefinite programming, is proposed. To showcase the methodology, the mitigation of the fluctuation kinetic energy in the unsteady wake behind a circular cylinder in the laminar regime at Re = 100, via controlled angular motions of the surface, is numerically investigated. A compact reduced-order model that resolves the long-term behaviour of the fluid flow and the effects of actuation, is first derived using proper orthogonal decomposition and Galerkin projection. In a full-information setting, feedback controllers are then designed to reduce the long-time average of the resolved kinetic energy associated with the limit cycle. These controllers are then implemented in direct numerical simulations of the actuated flow. Control performance, total energy efficiency and the physical control mechanisms identified are analysed in detail. Key elements of the methodology, implications and future work are finally discussed.