Review Article

A Practical Review on Linear and Nonlinear Global Approaches to Flow Instabilities

[+] Author and Article Information
D. Fabre, D. Ferreira Sabino, P. Bonnefis, J. Sierra

Institut de Mécanique des fluides
de Toulouse (IMFT),
Université de Toulouse,
Toulouse, France

V. Citro

Dipartimento di Ingegneria (DIIN),
Universitá degli Studi di Salerno,
Via Giovanni Paolo II, 132,
Fisciano 84084, Italy;
Institut de Mécanique des fluides
de Toulouse (IMFT),
Université de Toulouse,
Toulouse, France

F. Giannetti

Dipartimento di Ingegneria (DIIN),
Universitá degli Studi di Salerno,
Via Giovanni Paolo II, 132,
Fisciano 84084, Italy

M. Pigou

Institut de Mécanique
des fluides de Toulouse (IMFT),
Université de Toulouse,
Toulouse, France

Manuscript received June 15, 2018; final manuscript received January 23, 2019; published online February 13, 2019. Editor: Harry Dankowicz.

Appl. Mech. Rev 70(6), 060802 (Feb 13, 2019) (16 pages) Paper No: AMR-18-1070; doi: 10.1115/1.4042737 History: Received June 15, 2018; Revised January 23, 2019

This paper aims at reviewing linear and nonlinear approaches to study the stability of fluid flows. We provide a concise but self-contained exposition of the main concepts and specific numerical methods designed for global stability studies, including the classical linear stability analysis, the adjoint-based sensitivity, and the most recent nonlinear developments. Regarding numerical implementation, a number of ideas making resolution particularly efficient are discussed, including mesh adaptation, simple shift-invert strategy instead of the classical Arnoldi algorithm, and a simplification of the recent nonlinear self-consistent (SC) approach proposed by Mantič-Lugo et al. (2014, “Self-Consistent Mean Flow Description of the Nonlinear Saturation of the Vortex Shedding in the Cylinder Wake,” Phys. Rev. Lett., 113(8), p. 084501). An open-source software implementing all the concepts discussed in this paper is provided. The software is demonstrated for the reference case of the two-dimensional (2D) flow around a circular cylinder, in both incompressible and compressible cases, but is easily customizable to a variety of other flow configurations or flow equations.

Copyright © 2018 by ASME
Your Session has timed out. Please sign back in to continue.


Drazin, P. , and Reid, W. , 2004, Hydrodynamic Stability, Cambridge University Press, Cambridge, UK.
Giannetti, F. , and Luchini, P. , 2007, “ Structural Sensitivity of the First Instability of the Cylinder Wake,” J. Fluid Mech., 581, pp. 167–197.
Marquet, O. , Sipp, D. , and Jacquin, L. , 2008, “ Sensitivity Analysis and Passive Control of Cylinder Flow,” J. Fluid Mech., 615, pp. 221–252.
Sipp, D. , and Lebedev, A. , 2007, “ Global Stability of Base and Mean Flows: A General Approach and Its Applications to Cylinder and Open Cavity Flows,” J. Fluid Mech., 593, pp. 333–358.
Mantič-Lugo, V. , Arratia, C. , and Gallaire, F. , 2014, “ Self-Consistent Mean Flow Description of the Nonlinear Saturation of the Vortex Shedding in the Cylinder Wake,” Phys. Rev. Lett., 113(8), p. 084501. [PubMed]
Sipp, D. , Marquet, O. , Meliga, P. , and Barbagallo, A. , 2010, “ Dynamics and Control of Global Instabilities in Open-Flows: A Linearized Approach,” ASME Appl. Mech. Rev., 63(3), p. 030801.
Luchini, P. , and Bottaro, A. , 2014, “ Adjoint Equations in Stability Analysis,” Annu. Rev. Fluid Mech., 46(1), pp. 493–517.
Schmid, P. J. , and Brandt, L. , 2014, “ Analysis of Fluid Systems: Stability, Receptivity, SensitivityLecture Notes From the FLOW-NORDITA Summer School on Advanced Instability Methods for Complex Flows, Stockholm, Sweden, 2013,” ASME Appl. Mech. Rev., 66(2), p. 024803.
Hecht, F. , 2012, “ New Development in FreeFem++,” J. Numer. Math., 20(3–4), pp. 251–265. https://pdfs.semanticscholar.org/549b/bc057248bfbc67b15d65d8a7f22432ed7c43.pdf
Schmid, P. J. , and Henningson, D. S. , 2001, Stability and Transition in Shear Flows, Springer-Verlag, New York.
Chomaz, J.-M. , 2005, “ Global Instabilities in Spatially Developing Flows: Non-Normality and Nonlinearity,” Annu. Rev. Fluid Mech., 37(1), pp. 357–392.
Schmid, P. J. , 2007, “ Nonmodal Stability Theory,” Annu. Rev. Fluid Mech., 39(1), pp. 129–162.
Hill, D. , 1992, “ A Theoretical Approach for Analyzing the Restabilization of Wakes,” 30th Aerospace Sciences Meeting and Exhibit, Reno, NV, p. 67.
Arnoldi, W. E. , 1951, “ The Principle of Minimized Iterations in the Solution of the Matrix Eigenvalue Problem,” Q. Appl. Math., 9(1), pp. 17–29.
Hecht, F. , 1998, “ The Mesh Adapting Software: Bamg,” INRIA, Rocquencourt, France.
Kouhi, M. , Onate, E. , and Mavriplis, D. , 2015, “ Adjoint-Based Adaptive Finite Element Method for the Compressible Euler Equations Using Finite Calculus,” Aerosp. Sci. Technol., 46, pp. 422–435.
Marquet, O. , and Lesshafft, L. , 2015, “ Identifying the Active Flow Regions That Drive Linear and Nonlinear Instabilities,” arXiv:1508.07620. https://arxiv.org/abs/1508.07620
Luchini, P. , Giannetti, F. , and Citro, V. , 2017, “ Error Sensitivity to Refinement: A Criterion for Optimal Grid Adaptation,” Theor. Comput. Fluid Dyn., 31(5–6), pp. 595–605.
Fani, A. , Citro, V. , Giannetti, F. , and Auteri, F. , 2018, “ Computation of the Bluff-Body Sound Generation by a Self-Consistent Mean Flow Formulation,” Phys. Fluids, 30(3), p. 036102.
Gallaire, F. , Boujo, E. , Mantic-Lugo, V. , Arratia, C. , Thiria, B. , and Meliga, P. , 2016, “ Pushing Amplitude Equations Far From Threshold: Application to the Supercritical Hopf Bifurcation in the Cylinder Wake,” Fluid. Dyn. Res., 48(6), p. 061401.
Tchoufag, J. , Fabre, D. , and Magnaudet, J. , 2015, “ Weakly Nonlinear Model With Exact Coefficients for the Fluttering and Spiraling Motion of Buoyancy-Driven Bodies,” Phys. Rev. Lett., 115, p. 114501. [PubMed]
Carte, G. , Dušek, J. , and Fraunié, P. , 1995, “ A Spectral Time Discretization for Flows With Dominant Periodicity,” J. Comput. Phys., 120(2), pp. 171–183.
Carte, G. , Dušek, J. , and Fraunié, P. , 1996, “ Numerical Simulation of the Mechanisms Governing the Onset of the Bénard-Von Kármán Instability,” Int. J. Numer. Methods Fluids, 23(8), pp. 753–785.
Hall, K. C. , Thomas, J. P. , and Clark, W. S. , 2002, “ Computation of Unsteady Nonlinear Flows in Cascades Using a Harmonic Balance Technique,” AIAA J., 40(5), pp. 879–886.
Moulin, J. , Jolivet, P. , and Marquet, O. , 2017, “ Augmented Lagrangian Preconditioner for Linear Stability Analysis of Incompressible Fluid Flows on Large Configurations,” FreeFem++ Days, Paris, France, Dec. 14–15.
Williamson, C. H. , 1988, “ Defining a Universal and Continuous Strouhal–Reynolds Number Relationship for the Laminar Vortex Shedding of a Circular Cylinder,” Phys. Fluids, 31(10), pp. 2742–2744.
Navrose, N. , and Mittal, S. , 2016, “ Lock-In in Vortex-Induced Vibration,” J. Fluid Mech., 794, pp. 565–594.
Tchoufag, J. , Fabre, D. , and Magnaudet, J. , 2014, “ Global Linear Stability Analysis of the Wake and Path of Buoyancy-Driven Disks and Thin Cylinders,” J. Fluid Mech., 740, pp. 278–311.
Fabre, D. , Bonnefis, P. , Charru, F. , Russo, S. , Citro, V. , Giannetti, F. , and Luchini, P. , 2014, “ Application of Global Stability Approaches to Whistling Jets and Wind Instruments,” International Symposium on Musical Acoustics (ISMA'14), Le Mans, France, July, pp. 7–12.
Longobardi, R. , Fabre, D. , Bonnefis, P. , Citro, V. , Giannetti, F. , and Luchini, P. , 2018, “ Studying Sound Production in the Hole-Tone Configuration Using Compressible and Incompressible Global Stability Analyses,” IUTAM Symposium on Critical Flow Dynamics Involving Moving/Deformable Structures With Design Applications, Santorini, Greece, June 18–22. https://www.researchgate.net/publication/327894084_Studying_sound_production_in_the_hole-tone_configuration_using_compressible_and_incompressible_global_stability_analyses
Meliga, P. , Sipp, D. , and Chomaz, J.-M. , 2010, “ Effect of Compressibility on the Global Stability of Axisymmetric Wake Flows,” J. Fluid Mech., 660, pp. 499–526.
Chireux, V. , Fabre, D. , Risso, F. , and Tordjeman, P. , 2015, “ Oscillations of a Liquid Bridge Resulting From the Coalescence of Two Droplets,” Phys. Fluids, 27(6), p. 062103.
Mougel, J. , Fabre, D. , Lacaze, L. , and Bohr, T. , 2017, “ On the Instabilities of a Potential Vortex With a Free Surface,” J. Fluid Mech., 824, pp. 230–264.
Fabre, D. , Auguste, F. , and Magnaudet, J. , 2008, “ Bifurcations and Symmetry Breaking in the Wake of Axisymmetric Bodies,” Phys. Fluids, 20(5), p. 051702.


Grahic Jump Location
Fig. 1

Illustration of the usage of the stabfem software to produce an adapted mesh and study the base flowand the linear stability properties of the wake flow around a cylinder (extracted from script SCRIPT_CYLINDER_ALLFIGURES.m)

Grahic Jump Location
Fig. 2

Illustration of the stucture of mesh M2 (adapted to both the base flow and structural sensitivity)

Grahic Jump Location
Fig. 3

Recirculation length Lx (a) and nondimensional drag Fx (b) of the base flow over a cylinder as function of Re

Grahic Jump Location
Fig. 4

Base flow for the flow over a cylinder at Re = 60. Pressure field (color or grayscale levels) and streamlines (iso-levels of the streamfunction ψ).

Grahic Jump Location
Fig. 5

Growth rate σ (a) and Strouhal number St = ω/2π (b) as function of the Reynolds number

Grahic Jump Location
Fig. 6

Contour plot of the streamwise velocity component: (a) (direct) eigenmode and (b) adjoint mode at Re = 60

Grahic Jump Location
Fig. 7

Structural sensitivity Sw for the cylinder's wake at Re = 60. The red line represents the streamline bounding the recirculation region.

Grahic Jump Location
Fig. 8

Structure of the unstable eigenmode (streamwise velocity component) for the compressible flow around a cylinder (Re = 150, M =0.2)

Grahic Jump Location
Fig. 9

Structure of the unstable eigenmode (pressure component) for the compressible flow around a cylinder (Re = 150, M =0.2). Upper half: real part and lower half: imaginary part.

Grahic Jump Location
Fig. 10

Illustration of the procedure for nonlinear calculations using stabfem (extract from script SCRIPT_CYLINDER_ALLFIGURES.m)

Grahic Jump Location
Fig. 11

Comparison between the WNL results, the harmonic balance data (HB1), and baseflow/linear results: Strouhal number (a), mean drag (b), amplitude of oscillating lift (c), energy-amplitude of the nonlinear perturbation (d), and recirculation length of mean/base flows (e)

Grahic Jump Location
Fig. 12

Structure of the mean flow over a cylinder for Re = 60, as computed by the HB1 model. Color levels: pressure; streamlines.

Grahic Jump Location
Fig. 13

Illustration of the stucture of mesh M4 (adapted to both the base flow and direct eigenmode)

Grahic Jump Location
Fig. 14

Illustration of the implementation of the Newton algorithm for base-flow computation (extract from freefem++ program Newton2D.edp)

Grahic Jump Location
Fig. 15

Illustration of the implementation of the shift-invert algorithm for single eigenmode computation (extract from freefem++ program Stab2D.edp)



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In