Large time behavior and asymptotic stability of the two-dimensional Euler and linearized Euler equations

We study the asymptotic behavior and the asymptotic stability of the two-dimensional Euler equations and of the two-dimensional linearized Euler equations close to parallel flows. We focus on spectrally stable jet profiles $U(y)$ with stationary streamlines $y_{0}$ such that $U'(y_{0})=0$, a case that has not been studied previously. We describe a new dynamical phenomenon: the depletion of the vorticity at the stationary streamlines. An unexpected consequence, is that the velocity decays for large times with power laws, similarly to what happens in the case of the Orr mechanism for base flows without stationary streamlines. The asymptotic behaviors of velocity and the asymptotic profiles of vorticity are theoretically predicted and compared with direct numerical simulations. We argue on the asymptotic stability of these flow velocities even in the absence of any dissipative mechanisms.Comment: To be published in Physica D, nonlinear phenomena (accepted January 2010

Pattern formation in Hamiltonian systems with continuous spectra; a normal-form single-wave model

Pattern formation in biological, chemical and physical problems has received considerable attention, with much attention paid to dissipative systems. For example, the Ginzburg--Landau equation is a normal form that describes pattern formation due to the appearance of a single mode of instability in a wide variety of dissipative problems. In a similar vein, a certain "single-wave model" arises in many physical contexts that share common pattern forming behavior. These systems have Hamiltonian structure, and the single-wave model is a kind of Hamiltonian mean-field theory describing the patterns that form in phase space. The single-wave model was originally derived in the context of nonlinear plasma theory, where it describes the behavior near threshold and subsequent nonlinear evolution of unstable plasma waves. However, the single-wave model also arises in fluid mechanics, specifically shear-flow and vortex dynamics, galactic dynamics, the XY and Potts models of condensed matter physics, and other Hamiltonian theories characterized by mean field interaction. We demonstrate, by a suitable asymptotic analysis, how the single-wave model emerges from a large class of nonlinear advection-transport theories. An essential ingredient for the reduction is that the Hamiltonian system has a continuous spectrum in the linear stability problem, arising not from an infinite spatial domain but from singular resonances along curves in phase space whereat wavespeeds match material speeds (wave-particle resonances in the plasma problem, or critical levels in fluid problems). The dynamics of the continuous spectrum is manifest as the phenomenon of Landau damping when the system is ... Such dynamical phenomena have been rediscovered in different contexts, which is unsurprising in view of the normal-form character of the single-wave model

Absolute and convective instabilities of an inviscid compressible mixing layer: Theory and applications

This study aims to examine the effect of compressibility on unbounded and parallel shear flow linear instabilities. This analysis is of interest for industrial, geophysical, and astrophysical flows. We focus on the stability of a wavepacket as opposed to previous single-mode stability studies. We consider the notions of absolute and convective instabilities first used to describe plasma instabilities. The compressible-flow modal theory predicts instability whatever the Mach number. Spatial and temporal growth rates and Reynolds stresses nevertheless become strongly reduced at high Mach numbers. The evolution of disturbances is driven by Kelvin -Helmholtz instability that weakens in supersonic flows. We wish to examine the occurrence of absolute instability, necessary for the appearance of turbulent motions in an inviscid and compressible two-dimensional mixing layer at an arbitrary Mach number subject to a three-dimensional disturbance. The mixing layer is defined by a parametric family of mean-velocity and temperature profiles. The eigenvalue problem is solved with the help of a spectral method. We ascertain the effects of the distribution of temperature and velocity in the mixing layer on the transition between convective and absolute instabilities. It appears that, in most cases, absolute instability is always possible at high Mach numbers provided that the ratio of slow-stream temperature over fast-stream temperature may be less than a critical maximal value but the temporal growth rate present in the absolutely unstable zone remains small and tends to zero at high Mach numbers. The transition toward a supersonic turbulent regime is therefore unlikely to be possible in the linear theory. Absolute instability can be also present among low-Mach-number coflowing mixing layers provided that this same temperature ratio may be small, but nevertheless, higher than a critical minimal value. Temperature distribution within the mixing layer also has an effect on the growth rate, this diminishes when the slow stream is heated. These results are applied to the dynamics of mixing layers in the interstellar medium and to the dynamics of the heliopause, frontier between the interstellar medium, and the solar wind. (C) 2009 American Institute of Physics

Structure and stability of the compressible Stuart vortex

The structure and two- and three-dimensional stability properties of a linear array of compressible Stuart vortices (CSV; Stuart 1967; Meiron et al. 2000) are investigated both analytically and numerically. The CSV is a family of steady, homentropic, two-dimensional solutions to the compressible Euler equations, parameterized by the free-stream Mach number M_∞, and the mass flux _ inside a single vortex core. Known solutions have 0 < M_∞ < 1. To investigate the normal-mode stability of the generally spatially non-uniform CSV solutions, the linear partial-differential equations describing the time evolution of small perturbations to the CSV base state are solved numerically using a normal-mode analysis in conjunction with a spectral method. The effect of increasing M_∞ on the two main classes of instabilities found by Pierrehumbert & Widnall (1982) for the incompressible limit M_∞ → 0 is studied. It is found that both two- and three-dimensional subharmonic instabilities cease to promote pairing events even at moderate M_∞. The fundamental mode becomes dominant at higher Mach numbers, although it ceases to peak strongly at a single spanwise wavenumber. We also find, over the range of ε investigated, a new instability corresponding to an instability on a parallel shear layer. The significance of these instabilities to experimental observations of growth in the compressible mixing layer is discussed. In an Appendix, we study the CSV equations when ε is small and M_∞ is finite using a perturbation expansion in powers of ε. An eigenvalue determining the structure of the perturbed vorticity and density fields is obtained from a singular Sturm–Liouville problem for the stream-function perturbation at O(ε). The resulting small-amplitude steady CSV solutions are shown to represent a bifurcation from the neutral point in the stability of a parallel shear layer with a tanh-velocity profile in a compressible inviscid perfect gas at uniform temperature

CAT-generating mechanisms

The development of instability configurations; the transition from unstable growth of these configurations into turbulence; a description of the nature of that turbulence; the question of decay of turbulence; and the existence of what is called fossil turbulence are discussed

Cluster-based reduced-order modelling of a mixing layer

We propose a novel cluster-based reduced-order modelling (CROM) strategy of unsteady flows. CROM combines the cluster analysis pioneered in Gunzburger's group (Burkardt et al. 2006) and and transition matrix models introduced in fluid dynamics in Eckhardt's group (Schneider et al. 2007). CROM constitutes a potential alternative to POD models and generalises the Ulam-Galerkin method classically used in dynamical systems to determine a finite-rank approximation of the Perron-Frobenius operator. The proposed strategy processes a time-resolved sequence of flow snapshots in two steps. First, the snapshot data are clustered into a small number of representative states, called centroids, in the state space. These centroids partition the state space in complementary non-overlapping regions (centroidal Voronoi cells). Departing from the standard algorithm, the probabilities of the clusters are determined, and the states are sorted by analysis of the transition matrix. Secondly, the transitions between the states are dynamically modelled using a Markov process. Physical mechanisms are then distilled by a refined analysis of the Markov process, e.g. using finite-time Lyapunov exponent and entropic methods. This CROM framework is applied to the Lorenz attractor (as illustrative example), to velocity fields of the spatially evolving incompressible mixing layer and the three-dimensional turbulent wake of a bluff body. For these examples, CROM is shown to identify non-trivial quasi-attractors and transition processes in an unsupervised manner. CROM has numerous potential applications for the systematic identification of physical mechanisms of complex dynamics, for comparison of flow evolution models, for the identification of precursors to desirable and undesirable events, and for flow control applications exploiting nonlinear actuation dynamics.Comment: 48 pages, 30 figures. Revised version with additional material. Accepted for publication in Journal of Fluid Mechanic

Nonlinear evolution of the magnetized Kelvin-Helmholtz instability: from fluid to kinetic modeling

The nonlinear evolution of collisionless plasmas is typically a multi-scale process where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understanding cross-scale processes in plasmas. We report here the first comparative study of the evolution of a magnetized shear flow, through a variety of different plasma models by using magnetohydrodynamic, Hall-MHD, two-fluid, hybrid kinetic and full kinetic codes. Kinetic relaxation effects are discussed to emphasize the need for kinetic equilibriums to study the dynamics of collisionless plasmas in non trivial configurations. Discrepancies between models are studied both in the linear and in the nonlinear regime of the magnetized Kelvin-Helmholtz instability, to highlight the effects of small scale processes on the nonlinear evolution of collisionless plasmas. We illustrate how the evolution of a magnetized shear flow depends on the relative orientation of the fluid vorticity with respect to the magnetic field direction during the linear evolution when kinetic effects are taken into account. Even if we found that small scale processes differ between the different models, we show that the feedback from small, kinetic scales to large, fluid scales is negligable in the nonlinear regime. This study show that the kinetic modeling validates the use of a fluid approach at large scales, which encourages the development and use of fluid codes to study the nonlinear evolution of magnetized fluid flows, even in the colisionless regime