Intermittency and transition to chaos in the cubical lid-driven cavity flow

Transition from steady state to intermittent chaos in the cubical lid-driven flow is investigated numerically. Fully three-dimensional stability analyses have revealed that the flow experiences an Andronov-Poincar\'e-Hopf bifurcation at a critical Reynolds number $Re_c$ = 1914. As for the 2D-periodic lid-driven cavity flows, the unstable mode originates from a centrifugal instability of the primary vortex core. A Reynolds-Orr analysis reveals that the unstable perturbation relies on a combination of the lift-up and anti lift-up mechanisms to extract its energy from the base flow. Once linearly unstable, direct numerical simulations show that the flow is driven toward a primary limit cycle before eventually exhibiting intermittent chaotic dynamics. Though only one eigenpair of the linearized Navier-Stokes operator is unstable, the dynamics during the intermittencies are surprisingly well characterized by one of the stable eigenpairs.Comment: Accepted for publication in Fluid Dynamics Researc

A domain decomposition matrix-free method for global linear stability

This work is dedicated to the presentation of a matrix-free method for global linear stability analysis in geometries composed of multi-connected rectangular subdomains. An Arnoldi technique using snapshots in subdomains of the entire geometry combined with a multidomain linearized Direct Numerical Finite difference simulations based on an influence matrix for partitioning are adopted. The method is illustrated by three benchmark problems: the lid-driven cavity, the square cylinder and the open cavity flow. The efficiency of the method to extract large-scale structures in a multidomain framework is emphasized. The possibility to use subset of the full domain to recover the perturbation associated with the entire flow field is also highlighted. Such a method appears thus a promising tool to deal with large computational domains and three-dimensionality within a parallel architecture

Nonlinear control of unsteady finite-amplitude perturbations in the Blasius boundary-layer flow

The present work provides an optimal control strategy, based on the nonlinear Navier–Stokes equations, aimed at hampering the rapid growth of unsteady finite-amplitude perturbations in a Blasius boundary-layer flow. A variational procedure is used to find the blowing and suction control law at the wall providing the maximum damping of the energy of a given perturbation at a given target time, with the final aim of leading the flow back to the laminar state. Two optimally growing finite-amplitude initial perturbations capable of leading very rapidly to transition have been used to initialize the flow. The nonlinear control procedure has been found able to drive such perturbations back to the laminar state, provided that the target time of the minimization and the region in which the blowing and suction is applied have been suitably chosen. On the other hand, an equivalent control procedure based on the linearized Navier–Stokes equations has been found much less effective, being not able to lead the flow to the laminar state when finite-amplitude disturbances are considered. Regions of strong sensitivity to blowing and suction have been also identified for the given initial perturbations: when the control is actuated in such regions, laminarization is also observed for a shorter extent of the actuation region. The nonlinear optimal blowing and suction law consists of alternating wall-normal velocity perturbations, which appear to modify the core flow structures by means of two distinct mechanisms: (i) a wall-normal velocity compensation at small times; (ii) a rotation-counterbalancing effect al larger times. Similar control laws have been observed for different target times, values of the cost parameter, and streamwise extents of the blowing and suction zone, meaning that these two mechanisms are robust features of the optimal control strategy, provided that the nonlinear effects are taken into account

Shock wave instability and the carbuncle phenomenon: same intrinsic origin ?

The theoretical linear stability of a shock wave moving in an unlimited homogeneous environment has been widely studied during the last fifty years. Important results have been obtained by Dyakov (1954), Landau & Lifchitz (1959) and then by Swan & Fowles (1975) where the fluctuating quantities are written as normal modes. More recently, numerical studies on upwind finite difference schemes have shown some instabilities in the case of the motion of an inviscid perfect gas in a rectangular channel. The purpose of this paper is first to specify a mathematical formulation for the eigenmodes and to exhibit a new mode which was not found by the previous stability analysis of shock waves. Then, this mode is confirmed by numerical simulations which may lead to a new understanding of the so-called carbuncle phenomenon

Numerical Study of the Effect of Freestream Turbulence on by-pass Transition in a Boundary Layer

We use direct numerical simulations in the presence of free-stream turbulence having different values of intensity, T u, and integral length scale, L, in order to determine which kind of structures are involved in the path to transition of a boundary-layer flow. The main aim is to determine under which conditions the path to transition involves structures similar to the linear or non-linear optimal perturbations. For high values of T u and L, we observe a large-amplitude path to transition characterized by localized vortical structures and patches of high- and low-momentum fluctuations. Such a scenario is found to correlate well with the L and hairpin structures resulting from the time evolution of non-linear optimal perturbations, whereas, for lower T u and L, a larger correlation is found with respect to linear optimal disturbances. This indicates that a large-amplitude path to transition exists, different from the one characterized by elongated streaks undergoing secondary instability. To distinguish between the two transition scenarios, a simple parameter linked to the streamwise localisation of high- and low-momentum zones is introduced. Finally, an accurate law to predict the transition location is provided, taking into account both T u and L, valid for both the transition scenarios

Global and Koopman modes analysis of sound generation in mixing layers

It is now well established that linear and nonlinear instability waves play a significant role in the noise generation process for a wide variety of shear flows such as jets or mixing layers. In that context, the problem of acoustic radiation generated by spatially growing instability waves of two-dimensional subsonic and supersonic mixing layers are revisited in a global point of view, i.e., without any assumption about the base flow, in both a linear and a nonlinear framework by using global and Koopman mode decompositions. In that respect, a timestepping technique based on disturbance equations is employed to extract the most dynamically relevant coherent structures for both linear and nonlinear regimes. The present analysis proposes thus a general strategy for analysing the near-field coherent structures which are responsible for the acoustic noise in these configurations. In particular, we illustrate the failure of linear global modes to describe the noise generation mechanism associated with the vortex pairing for the subsonic regime whereas they appropriately explain the Mach wave radiation of instability waves in the supersonic regime. By contrast, the Dynamic Mode Decomposition (DMD) analysis captures both the near-field dynamics and the far-field acoustics with a few number of modes for both configurations. In addition, the combination of DMD and linear global modes analyses provides new insight about the influence on the radiated noise of nonlinear interactions and saturation of instability waves as well as their interaction with the mean flow

Comparison of Subgrid-scale Viscosity Models and Selective Filtering Strategy for Large-eddy Simulations

Explicitly filtered large-eddy simulations (LES), combining high-accuracy schemes with the use of a selective filtering without adding an explicit subgrid-scales (SGS) model, are carried out for the Taylor-Green-vortex and the supersonic-boundary-layer cases. First, the present approach is validated against direct numerical simulation (DNS) results. Subsequently, several SGS models are implemented in order to investigate if they can improve the initial filter-based methodology. It is shown that the most accurate results are obtained when the filtering is used alone as an implicit model, and for a minimal cost. Moreover, the tests for the Taylor-Green vortex indicate that the discretization error from the numerical methods, notably the dissipation error from the high-order filtering, can have a greater influence than the SGS models

Rapid path to transition via nonlinear localized optimal perturbations in a boundary-layer flow

Recent studies have suggested that in some cases transition can be triggered by some purely nonlinear mechanisms. Here we aim at verifying such an hypothesis, looking for a localized perturbation able to lead a boundary-layer flow to a chaotic state, following a nonlinear route. Nonlinear optimal localized perturbations have been computed by means of an energy optimization which includes the nonlinear terms of the Navier- Stokes equations. Such perturbations lie on the turbulent side of the laminar-turbulent boundary, whereas, for the same value of the initial energy, their linear counterparts do not. The evolution of these perturbations toward a turbulent flow involves the presence of streamwise-inclined vortices at short times and of hairpin structures prior to breakdown

Time-stepping and Krylov methods for large-scale instability problems

With the ever increasing computational power available and the development of high-performances computing, investigating the properties of realistic very large-scale nonlinear dynamical systems has been become reachable. It must be noted however that the memory capabilities of computers increase at a slower rate than their computational capabilities. Consequently, the traditional matrix-forming approaches wherein the Jacobian matrix of the system considered is explicitly assembled become rapidly intractable. Over the past two decades, so-called matrix-free approaches have emerged as an efficient alternative. The aim of this chapter is thus to provide an overview of well-grounded matrix-free methods for fixed points computations and linear stability analyses of very large-scale nonlinear dynamical systems.Comment: To appear in "Computational Modeling of Bifurcations and Instabilities in Fluid Mechanics", eds. A. Gelfgat, Springe

Collective secondary instabilities: an application to three-dimensional boundary-layer flow

In some linearly unstable flows, secondary instability is found to have a much larger wavelength than that of the primary unstable modes, so that it cannot be recovered with a classical Floquet analysis. In this work, we apply a new formulation for capturing secondary instabilities coupling multiple length scales of the primary mode. This formulation, based on two-dimensional stability analysis coupled with a Bloch waves formalism originally described in Schmid et al. (2017), allows to consider high-dimensional systems resulting from several repetitions of a periodic unit, by solving an eigenproblem of much smaller size. Collective instabilities coupling multiple periodic units can be thus retrieved. The method is applied on the secondary stability of a swept boundary-layer flow subject to stationary cross-flow vortices, and compared with Floquet analysis. Two multi-modal instabilities are recovered: for streamwise wavenumber $\alpha_v$ close to zero, approximately twelve sub-units are involved in large-wavelength oscillations; whereas a staggered pattern, characteristic of subharmonic instabilities, is observed for $\alpha_v = 0.087$