Search

A family of compressible exact periodic solutions is reported for the first time in an open cavity flow setup. These are found using a novel framework which permits the computation of such solutions in an arbitrary complex geometry. The periodic orbits arise from a synchronised concatenation of convective and acoustic events which strongly depend on the Mach number. This flow-acoustic interaction furnishes the periodic solutions with a remarkable stability and it is found to completely dominate the system's dynamics and the sound directivity. The periodic orbits, which could be called `exact Rossiter modes', collapse with a family of equilibrium solutions at a subcritical Hopf bifurcation, occurring in the quasi-incompressible regime. This shows compressibility has a destabilising effect in cavity flows, which we analyse in detail. By establishing a connection with previous 2D and 3D stability studies of cavity flows, we are able to isolate the effect of purely compressible two-dimensional flow phenomena across Mach number. A linear stability analysis of the equilibria provides insight into the compressible flow mechanisms responsible for the instability. A close look at the adjoint modes suggests that an eigenvalue merge occurs at a Mach number between 0.35 and 0.4, which boosts the receptivity of the leading mode and determines the onset of the unstable character of the system. The effect of the choice of base flow over the transition dynamics is also discussed, where in the present case, the frequencies associated to the leading eigenmodes show a strong connection with the frequencies of the periodic orbits at the same Mach numbers.

The relationship between Koopman mode decomposition, resolvent mode decomposition and exact invariant solutions of the Navier-Stokes equations is clarified. The correspondence rests upon the invariance of the system operators under symmetry operations such as spatial translation. The usual interpretation of the Koopman operator is generalised to permit combinations of such operations, in addition to translation in time. This invariance is related to the spectrum of a spatio-temporal Koopman operator, which has a travelling wave interpretation. The relationship leads to a generalisation of dynamic mode decomposition, in which symmetry operations are applied to restrict the dynamic modes to span a subspace subject to those symmetries. The resolvent is interpreted as the mapping between the Koopman modes of the Reynolds stress divergence and the velocity field. It is shown that the singular vectors of the resolvent (the resolvent modes) are the optimal basis in which to express the velocity field Koopman modes where the latter are not a priori known.

A novel adjoint-based framework oriented to optimal flow control in compressible direct numerical simulations is presented. Also, a new formulation of the adjoint characteristic boundary conditions is introduced, which enhances the stability of the adjoint simulations. The flow configuration chosen as a case study consists of a two dimensional open cavity flow with aspect ratio $L/H=3$ and Reynolds number $Re=5000$. This flow configuration is of particular interest, as the turbulent and chaotic nature of separated flows pushes the adjoint approach to its limit. The target of the flow actuation, defined as cost, is the reduction of the pressure fluctuations at the sensor location. To exploit the advantages of the adjoint method, a large number of control parameters is used. The control consists of an actuating sub-domain where a two-dimensional body force is applied at every point within the sub-volume. This results in a total of $2.256 \cdot 10^6$ control parameters. The final actuation achieved a successful reduction of the cost of $79.6\%$, by altering the directivity of the sound radiated by the trailing edge of the cavity and breaking up the large shear layer instabilities into smaller structures.

This paper describes a robust linear time-invariant output-feedback control strategy to reduce turbulent fluctuations, and therefore skin-friction drag, in wall-bounded turbulent fluid flows, that nonetheless gives performance guarantees in the nonlinear turbulent regime. The novel strategy is effective in reducing the supply of available energy to feed the turbulent fluctuations, expressed as reducing a bound on the supply rate to a quadratic storage function. The nonlinearity present in the equations that govern the dynamics of the flow is known to be passive and can be considered as a feedback forcing to the linearised dynamics (a Lur'e decomposition). Therefore, one is only required to control the linear dynamics in order to make the system close to passive. The ten most energy-producing spatial modes of a turbulent channel flow were identified. Passivity-based controllers were then generated to control these modes. The controllers require measurements of streamwise and spanwise wall-shear stress, and they actuate via wall transpiration. Nonlinear direct numerical simulations demonstrated that these controllers were capable of significantly reducing the turbulent energy and skin-friction drag of the flow.

We report that many exact invariant solutions of the Navier-Stokes equations for both pipe and channel flows are well represented by just few modes of the model of McKeon & Sharma J. Fl. Mech. 658, 356 (2010). This model provides modes that act as a basis to decompose the velocity field, ordered by their amplitude of response to forcing arising from the interaction between scales. The model was originally derived from the Navier-Stokes equations to represent turbulent flows and has been used to explain coherent structure and to predict turbulent statistics. This establishes a surprising new link between the two distinct approaches to understanding turbulence., Comment: 21 figures

This paper addresses the problem of obtaining low-order models of fluid flows for the purpose of designing robust feedback controllers. This is challenging since whilst many flows are governed by a set of nonlinear, partial differential-algebraic equations (the Navier-Stokes equations), the majority of established control theory assumes models of much greater simplicity, in that they are firstly: linear, secondly: described by ordinary differential equations, and thirdly: finite-dimensional. Linearisation, where appropriate, overcomes the first disparity, but attempts to reconcile the remaining two have proved difficult. This paper addresses these two problems as follows. Firstly, a numerical approach is used to project the governing equations onto a divergence-free basis, thus converting a system of differential-algebraic equations into one of ordinary differential equations. This dispenses with the need for analytical velocity-vorticity transformations, and thus simplifies the modelling of boundary sensing and actuation. Secondly, this paper presents a novel and straightforward approach for obtaining suitable low-order models of fluid flows, from which robust feedback controllers can be synthesised that provide a priori guarantees of robust performance when connected to the (infinite-dimensional) linearised flow system. This approach overcomes many of the problems inherent in approaches that rely upon model-reduction. To illustrate these methods, a perturbation shear stress controller is designed and applied to plane channel flow, assuming arrays of wall mounted shear-stress sensors and transpiration actuators. DNS results demonstrate robust attenuation of the perturbation shear-stresses across a wide range of Reynolds numbers with a single, linear controller., Comment: Accepted for publication in the Journal of Fluid Mechanics

An analytical framework for studying the logarithmic region of turbulent channels is formulated. We build on recent findings (Moarref et al., J. Fluid Mech., 734, 2013) that the velocity fluctuations in the logarithmic region can be decomposed into a weighted sum of geometrically self-similar resolvent modes. The resolvent modes and the weights represent the linear amplification mechanisms and the scaling influence of the nonlinear interactions in the Navier-Stokes equations (NSE), respectively (McKeon & Sharma, J. Fluid Mech., 658, 2010). Originating from the NSE, this framework provides an analytical support for Townsend's attached-eddy model. Our main result is that self-similarity enables order reduction in modeling the logarithmic region by establishing a quantitative link between the self-similar structures and the velocity spectra. Specifically, the energy intensities, the Reynolds stresses, and the energy budget are expressed in terms of the resolvent modes with speeds corresponding to the top of the logarithmic region. The weights of the triad modes -the modes that directly interact via the quadratic nonlinearity in the NSE- are coupled via the interaction coefficients that depend solely on the resolvent modes (McKeon et al., Phys. Fluids, 25, 2013). We use the hierarchies of self-similar modes in the logarithmic region to extend the notion of triad modes to triad hierarchies. It is shown that the interaction coefficients for the triad modes that belong to a triad hierarchy follow an exponential function. The combination of these findings can be used to better understand the dynamics and interaction of flow structures in the logarithmic region. The compatibility of the proposed model with theoretical and experimental results is further discussed., Comment: Submitted to J. Fluid Mech

We study the Reynolds number scaling and the geometric self-similarity of a gain-based, low-rank approximation to turbulent channel flows, determined by the resolvent formulation of McKeon & Sharma (2010), in order to obtain a description of the streamwise turbulence intensity from direct consideration of the Navier-Stokes equations. Under this formulation, the velocity field is decomposed into propagating waves (with single streamwise and spanwise wavelengths and wave speed) whose wall-normal shapes are determined from the principal singular function of the corresponding resolvent operator. Using the accepted scalings of the mean velocity in wall-bounded turbulent flows, we establish that the resolvent operator admits three classes of wave parameters that induce universal behavior with Reynolds number on the low-rank model, and which are consistent with scalings proposed throughout the wall turbulence literature. In addition, it was shown that a necessary condition for geometrically self-similar resolvent modes is the presence of a logarithmic turbulent mean velocity. We identify the scalings that constitute hierarchies of self-similar modes that are parameterized by the critical wall-normal location where the speed of the mode equals the local turbulent mean velocity. For the rank-1 model subject to broadband forcing, the integrated streamwise energy density takes a universal form which is consistent with the dominant near-wall turbulent motions. When the shape of the forcing is optimized to enforce matching with results from direct numerical simulations at low turbulent Reynolds numbers, further similarity appears. Representation of these weight functions using similarity laws enables prediction of the Reynolds number and wall-normal variations of the streamwise energy intensity at high Reynolds numbers (${Re}_\tau \approx 10^3 - 10^{10}$)., Comment: The paper is to appear in the Journal of Fluid Mechanics