Review Articles

Instability in Stratified Shear Flow: Review of a Physical Interpretation Based on Interacting Waves

[+] Author and Article Information
Jeffrey R. Carpenter

Eawag, Swiss Federal Institute of
Aquatic Science and Technology,
Surface Waters – Research and Management,
Kastanienbaum 6047, Switzerland;
Department of Geology and Geophysics,
Yale University,
New Haven, CT 06520-8109
e-mail: jeffcarp@gmail.com

Edmund W. Tedford

Marine Science Institute,
University of California, Santa Barbara,
Santa Barbara, CA 93106-9620
e-mail: ttedford@eos.ubc.ca

Eyal Heifetz

Department of Geophysical,
Atmospheric and Planetary Sciences,
Tel Aviv University,
Tel Aviv 69978, Israel;
Department of Meteorology (MISU),
Stockholm University,
Stockholm SE-106 91, Sweden
e-mail: eyalh@cyclone.tau.ac.il

Gregory A. Lawrence

Department of Civil Engineering,
University of British Columbia,
Vancouver, BC, V6T 1Z4, Canada
e-mail: lawrence@civil.ubc.ca

Here we are neglecting a part of the solution referred to as the continuous spectrum, which is not essential for the purposes of this paper. For more information the reader is referred to Schmid and Henningson [4].

The Rayleigh and Fjørtoft theorems can also be interpreted in terms of the conservation of pseudomomentum and pseudoenergy, respectively, and the interested reader is referred to Heifetz et al. [11] for further details.

Although the work of Garcia [23] was published 6 years previous to Holmboe's [14] paper, Garcia directly acknowledges Holmboe with the wave interaction concept in an acknowledgments section that appears (unusually) in the main body of the paper.

The stability of this particular flow may also be interpreted in terms of an energetics perspective. The stable flow is due to the lack of a shear growth source—the waves must either be embedded within the shear or at its sides in order to be able to extract energy from it. See Harnik and Heifetz [37] and Rabinovich et al. [29] for further details.

It is more appropriate to refer to this vorticity-vorticity wave interaction as the Rayleigh [22] mechanism, however, it has become conventional in the stratified shear flow literature to use KH after Kelvin [40] and Helmholtz [41], and we shall stick with this notation in the remainder of the paper.

1Corresponding author.

Manuscript received May 20, 2011; final manuscript received October 17, 2012; published online January 23, 2013. Editor: Harry Dankowicz.

Appl. Mech. Rev 64(6), 060801 (Jan 23, 2013) (17 pages) doi:10.1115/1.4007909 History: Received May 20, 2011; Revised October 17, 2012

Instability in homogeneous and density stratified shear flows may be interpreted in terms of the interaction of two (or more) otherwise free waves in the velocity and density profiles. These waves exist on gradients of vorticity and density, and instability results when two fundamental conditions are satisfied: (I) the phase speeds of the waves are stationary with respect to each other (“phase-locking“), and (II) the relative phase of the waves is such that a mutual growth occurs. The advantage of the wave interaction approach is that it provides a physical interpretation to shear flow instability. This paper is largely intended to purvey the basics of this physical interpretation to the reader, while both reviewing and consolidating previous work on the topic. The interpretation is shown to provide a framework for understanding many classical and nonintuitive results from the stability of stratified shear flows, such as the Rayleigh and Fjørtoft theorems, and the destabilizing effect of an otherwise stable density stratification. Finally, we describe an application of the theory to a geophysical-scale flow in the Fraser River estuary.

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


Lindzen, R., 1988, “Instability of Plane Parallel Shear Flow (Toward a Mechanistic Picture of How it Works),” PAGEOPH, 126, pp. 103–121. [CrossRef]
Drazin, P., and Reid, W., 1982, Hydrodynamic Stability, Cambridge University Press, Cambridge, UK.
Betchov, R., and Criminale, W., 1967, Stability of Parallel Flows, Academic Press, New York.
Schmid, P., and Henningson, D., 2001, Stability and Transition in Shear Flows, Springer, New York.
Koppel, D., 1964, “On the Stability of a Thermally Stratified Fluid Under the Action of Gravity,” J. Math. Phys., 5, pp. 963–982. [CrossRef]
Smyth, W., and Peltier, W., 1991, “Instability and Transition in Finite-Amplitude Kelvin-Helmholtz and Holmboe Waves,” J. Fluid Mech., 228, pp. 387–415. [CrossRef]
Chagelishvili, G. D., Khujadze, G. R., Lominadze, J. G., and Rogava, A. D., 1997, “Acoustic Waves in Unbounded Shear Flows,” Phys. Fluids, 9(7), pp. 1955–1962. [CrossRef]
Alexakis, A., Young, Y., and Rosner, R., 2002, “Shear Instability of Fluid Interfaces: Stability Analysis,” Phys. Rev. E, 65, p. 026313. [CrossRef]
Umurhan, O., and Heifetz, E., 2007, “Holmboe Modes Revisited,” Phys. Fluids, 19, p. 064102. [CrossRef]
Bretherton, F., 1966, “Baroclinic Instability and the Short Wavelength Cut-Off in Terms of Potential Vorticity,” Q. J. R. Meteorol. Soc., 92, pp. 335–345. [CrossRef]
Heifetz, E., Bishop, C., Hoskins, B., and Methven, J., 2004, “The Counter-Propagating Rossby-Wave Perspective on Baroclinic Instability. I: Mathematical Basis,” Q. J. R. Meteorol. Soc., 130, pp. 211–231. [CrossRef]
Smyth, W., and Peltier, W., 1990, “Three-Dimensional Primary Instabilities of a Stratified, Dissipative, Parallel Flow,” Geophys. Astrophys. Fluid Dyn., 52, pp. 249–261. [CrossRef]
Gelfgat, A. Y., and Kit, E., 2006, “Spatial Versus Temporal Instabilities in a Parametrically Forced Stratified Mixing Layer,” J. Fluid Mech., 552, pp. 189–227. [CrossRef]
Holmboe, J., 1962, “On the Behavior of Symmetric Waves in Stratified Shear Layers,” Geofys. Publ., 24, pp. 67–112.
Baines, P., and Mitsudera, H., 1994, “On the Mechanism of Shear Flow Instabilities,” J. Fluid Mech., 276, pp. 327–342. [CrossRef]
Kundu, P., Cohen, I., and Hu, H., 2004, Fluid Mechanics, 2nd ed., Elsevier, Amsterdam.
Caulfield, C., 1994, “Multiple Linear Instability of Layered Stratified Shear Flow,” J. Fluid Mech., 258, pp. 255–285. [CrossRef]
Heifetz, E., Bishop, C. H., and Alpert, P., 1999, “Counter-Propagating Rossby Waves in the Barotropic Rayleigh Model of Shear Instability,” Q. J. R. Meteorol. Soc., 125, pp. 2835–2853. [CrossRef]
Heifetz, E., and Methven, J., 2005, “Relating Optimal Growth to Counterpropagating Rossby Waves in Shear Instability,” Phys. Fluids, 17, p. 064107. [CrossRef]
Heifetz, E., Reuveni, Y., Gelfgat, A., Kit, E., and Methven, J., 2006, “Counterpropagating Rossby Wave Perspective on Kelvin Helmholtz Instability as a Limiting Case of a Rayleigh Shear Layer With Zero Width,” Phys. Fluids, 18, p. 018101. [CrossRef]
Farrell, B., 1988, “Optimal Excitation of Perturbations in Viscous Shear Flow,” Phys. Fluids, 31, pp. 2093–2101. [CrossRef]
Rayleigh, J., 1880, “On the Stability, or Instability, of Certain Fluid Motions,” Proc. Lond. Math. Soc., 12, pp. 57–70.
Garcia, R., 1956, “Barotropic Waves in Straight Parallel Flow With Curved Velocity Profile,” Tellus, 8, pp. 82–93. [CrossRef]
Baines, P., 1995, Topographic Effects in Stratified Flows, Cambridge University Press, Cambridge, UK.
Michalke, A., 1964, “On the Inviscid Instability of Hyperbolic-Tangent Velocity Profile,” J. Fluid Mech., 19, pp. 543–556. [CrossRef]
Hazel, P., 1972, “Numerical Studies of the Stability of Inviscid Stratified Shear Flows,” J. Fluid Mech., 51, pp. 39–61. [CrossRef]
Carpenter, J., Balmforth, N., and Lawrence, G., 2010, “Identifying Unstable Modes in Stratified Shear Layers,” Phys. Fluids, 22, p. 054104. [CrossRef]
Harnik, N., Heifetz, E., Umurhan, O., and Lott, F., 2008, “A Buoyancy-Vorticity Wave Interaction Approach to Stratified Shear Flow,” J. Atmos. Sci., 65, pp. 2615–2630. [CrossRef]
Rabinovich, A., Umurhan, O., Harnik, N., Lott, F., and Heifetz, E., 2011, “Vorticity Inversion and Action-at-a-Distance Instability in Stably Stratified Shear Flow,” J. Fluid Mech., 670, pp. 301–325. [CrossRef]
Balmforth, N. J., Roy, A., and Caulfield, C. P., 2012, “Dynamics of Vorticity Defects in Stratified Shear Flow,” J. Fluid Mech., 694, pp. 292–331. [CrossRef]
Craik, A., 1985, Wave Interactions and Fluid Flows, Cambridge University Press, Cambridge, UK.
Cairns, R., 1979, “The Role of Negative Energy Waves in Some Instabilities of Parallel Flows,” J. Fluid Mech., 92, pp. 1–14. [CrossRef]
Fjørtoft, R., 1950, “Application of Integral Theorems in Deriving Criteria of Stability of Laminar Flows and for the Baroclinic Circular Vortex,” Geofys. Publ., 17, pp. 1–52.
Drazin, P., 2002, Introduction to Hydrodynamic Stability, 1st ed., Cambridge University Press, Cambridge, UK.
Howard, L., and Maslowe, S., 1973, “Stability of Stratified Shear Flows,” Boundary-Layer Meteorol., 4, pp. 511–523. [CrossRef]
Redekopp, L., 2001, “Elements of Instability Theory for Environmental Flows,” Environmental Stratified Flows, Kluwer, Boston.
Harnik, N., and Heifetz, E., 2007, “Relating Overreflection and Wave Geometry to the Counterpropagating Rossby Wave Perspective: Toward a Deeper Understanding of Shear Instability,” J. Atmos. Sci., 64, pp. 2238–2261. [CrossRef]
Taylor, G., 1931, “Effect of Variation in Density on the Stability of Superposed Streams of Fluid,” Proc. R. Soc. Lond. A, 132, pp. 499–523. [CrossRef]
Caulfield, C., Peltier, W., Yoshida, S., and Ohtani, M., 1995, “An Experimental Investigation of the Instability of a Shear Flow With Multilayered Density Stratification,” Phys. Fluids, 7, pp. 3028–3041. [CrossRef]
Kelvin, W., 1871, “Hydrokinetic Solutions and Observations,” Philos. Mag., 42, pp. 362–377.
Helmholtz, H., 1868, “On Discontinuous Movements of Fluids,” Philos. Mag., 36, pp. 337–346.
Lawrence, G., Browand, F., and Redekopp, L., 1991, “The Stability of a Sheared Density Interface,” Phys. Fluids, 3(10), pp. 2360–2370. [CrossRef]
Alexakis, A., 2005, “On Holmboe's Instability for Smooth Shear and Density Profiles,” Phys. Fluids, 17, p. 084103. [CrossRef]
Haigh, S., and Lawrence, G., 1999, “Symmetric and Nonsymmetric Holmboe Instabilities in an Inviscid Flow,” Phys. Fluids, 11(6), pp. 1459–1468. [CrossRef]
Tedford, E., Carpenter, J., Pawlowicz, R., Pieters, R., and Lawrence, G., 2009, “Observation and Analysis of Shear Instability in the Fraser River Estuary,” J. Geophys. Res., 114, p. C11006. [CrossRef]
Smyth, W., Klaassen, G., and Peltier, W., 1988, “Finite Amplitude Holmboe Waves,” Geophys. Astrophys. Fluid Dyn., 43, pp. 181–222. [CrossRef]
Zhu, D., and Lawrence, G., 2001, “Holmboe's Instability in Exchange Flows,” J. Fluid Mech., 429, pp. 391–409. [CrossRef]
Tedford, E., Pieters, R., and Lawrence, G., 2009, “Symmetric Holmboe Instabilities in a Laboratory Exchange Flow,” J. Fluid Mech., 636, pp. 137–153. [CrossRef]
Miles, J., 1961, “On the Stability of Heterogeneous Shear Flows,” J. Fluid Mech., 10, pp. 496–508. [CrossRef]
Howard, L., 1961, “Note on a Paper of John W. Miles,” J. Fluid Mech., 10, pp. 509–512. [CrossRef]
Thorpe, S., 1973, “Experiments on Instability and Turbulence in a Stratified Shear Flow,” J. Fluid Mech., 61, pp. 731–751. [CrossRef]
Caulfield, C., and Peltier, W., 2000, “The Anatomy of the Mixing Transition in Homogenous and Stratified Free Shear Layers,” J. Fluid Mech., 413, pp. 1–47. [CrossRef]
Smyth, W., and Winters, K., 2003, “Turbulence and Mixing in Holmboe Waves,” J. Phys. Oceanogr., 33, pp. 694–711. [CrossRef]
Smyth, W., Carpenter, J., and Lawrence, G., 2007, “Mixing in Symmetric Holmboe Waves,” J. Phys. Oceanogr., 37, pp. 1566–1583. [CrossRef]
Carpenter, J., Lawrence, G., and Smyth, W., 2007, “Evolution and Mixing of Asymmetric Holmboe Instabilities,” J. Fluid Mech., 582, pp. 103–132. [CrossRef]


Grahic Jump Location
Fig. 1

(a) Piecewise profiles used to demonstrate vorticity and internal gravity waves. The interface of density and vorticity is located at zℓ. (b) The vertical velocity eigenfunction w∧(z) for vorticity and internal gravity waves. Since w∧(z) is generally complex it has both an amplitude |w∧(z)|=|A|e-k|z-zℓ| and a phase given by tan(θw)=Im(w∧)/Re(w∧)=Im(A)/Re(A).

Grahic Jump Location
Fig. 2

Dispersion relation for the profiles in Fig. 1 showing (a) the phase speed cr(k), and (b) the frequency ωr(k) for vorticity and internal gravity waves. The vorticity wave mode shown here has a leftward intrinsic propagation and is denoted by V- (thin line), whereas the gravity wave modes are indicated by G± (thick lines) and have intrinsic propagation in both directions.

Grahic Jump Location
Fig. 3

Structure of the vorticity wave broken down into three important fields: interface displacement η˜, vorticity perturbation q˜, and the vertical velocity w˜. These fields are all sinusoidal with wavenumber k and illustrated on the right. Throughout the remainder of the paper we will use the more compact notation shown on the left, with circular arrows denoting crests (clockwise) and troughs (counterclockwise) in q˜ and similarly with the vertical arrows denoting crests and troughs in the w˜-field. These diagrams will be referred to as wave field diagrams.

Grahic Jump Location
Fig. 4

Structure of the rightward-propagating internal gravity wave broken down into its important fields. The notation is the same as in Fig. 3 except we have also plotted the generation of baroclinic vorticity, described by Eq. (16), as the dashed line. In addition to the rotational arrows, the q˜-field is indicated by the gray shading. This is because once wave interactions are accounted for the q˜-field is not generally directly related to the η˜-field, as is the case in the vorticity wave (see Sec. 5).

Grahic Jump Location
Fig. 5

Illustration of the relationship between the η˜ and w˜ fields in a growing and rightward-propagating interfacial wave

Grahic Jump Location
Fig. 6

The piecewise shear layer profile on left along with the eigenfunction amplitude of the two interface contributions to the total eigenfunctions |w∧(z)| and |ψ∧(z)|, which are related by w∧=-ikψ∧. The wave field diagram illustrating the unstable interaction of vorticity waves is shown on the right for the particular value of k = 0.4, corresponding to the maximum growth rate. As in Fig. 3, the vertical arrows indicate the peaks and troughs in the w˜-fields at each interface, due only to the displacement of that interface. In addition, we also show the vorticity q˜ and displacement η˜ of each interface.

Grahic Jump Location
Fig. 7

Dispersion relation for the piecewise shear layer. Upper panel shows the real frequency ωr(k) of the dispersion relation with the dashed line indicating unstable (ℜ) modes, solid lines denote the two stable (S) modes, and the gray lines correspond to the isolated vorticity waves of the upper (V1-) and lower (V2+) interfaces. The bottom panel shows the growth rate ωi(k) of the unstable mode.

Grahic Jump Location
Fig. 8

(a) The smooth `tanh’ shear layer velocity profile U*(z*) (gray) and the associated vorticity gradient distribution U*''(z*) (black). (b) Comparison between the growth rates of the smooth tanh-profile and the piecewise linear profile. Growth rates for both the traditional scaling of U*(z*) (dashed line) and the vorticity interface scaling (thick solid line) are shown, along with the piecewise result (thin line) for comparison.

Grahic Jump Location
Fig. 9

(a) Dispersion relation ωr(k) for two stable interacting vorticity waves used to illustrate the Rayleigh inflection point condition. As in Fig. 2 the gray lines indicate the dispersion of waves in isolation from each other. (b) Wave field diagram showing that it is not possible for the waves to cause mutual growth in one another—regardless of the relative phase difference. Here an arbitrary phase difference has been chosen for illustration purposes.

Grahic Jump Location
Fig. 10

(a) Dispersion relation for two stable vorticity waves that demonstrate Fjørtoft's extension (notation as in Fig. 2). (b) Wave field diagram showing the reinforcement of the intrinsic wave propagation by the advection of the mean flow profile that is stable by Fjørtoft's extension.

Grahic Jump Location
Fig. 11

Profiles of U(z) and Q(z) for the triangular jet flow, together with the wave field diagram at k = 0.5

Grahic Jump Location
Fig. 12

Dispersion relation for the triangular jet. The frequencies of the isolated vorticity waves (gray curves, V1,2,3), the stable waves (solid curves, S1,2), and the unstable jet mode (dashed curve, J) are shown in (a), with the growth rate in (b).

Grahic Jump Location
Fig. 13

(a) Profiles of U(z) and ρ¯(z) that demonstrate the unstable Holmboe interaction between a vorticity wave (at z = 1) and an internal gravity wave (located at z = 0). (b) The wave field diagram showing the structure of the unstable mode with J = 2, k = 1.87. The gray shading is used to indicate the q˜2 field, which may have a phase difference from the η˜2 field. (c) Stability diagram showing the Holmboe unstable region. Contours are of growth rate (ωi) with an interval of 0.03, and the thick solid lines indicate the stability boundary with stable regions in gray. The resonance condition is denoted by the dot-dashed line and follows very closely the unstable region. An example of the dispersion relation ωr(k) shown in (d) is taken at J = 2.

Grahic Jump Location
Fig. 14

(a) Profiles that demonstrate a stable interaction between a gravity wave and a vorticity wave. An example of the dispersion relation ωr(k) is shown in (b) for J = 0.25. See Fig. 2 for further notation.

Grahic Jump Location
Fig. 15

(a) Dimensionless profile of background velocity U(z) that demonstrates the unstable Taylor–Caulfield interaction between two internal gravity waves (at z=±1). (b) The wave field diagram showing the structure of the unstable mode with J = 2, k = 1. (c) Stability diagram showing the Taylor–Caulfield unstable region (see Fig. 13(c) for notation). An example of the dispersion relation ωr(k) shown in (d) is taken at J = 2.

Grahic Jump Location
Fig. 16

(a) Profiles for Holmboe's [14] symmetric (ε = 0) stratified shear layer. (b) The associated symmetric stability diagram where the gray area represents stable regions, with the thick dark line denoting the stability boundary, as well as the transition from propagating to stationary modes. The thin dark lines are contours of growth rate, ωi, with spacing of 0.03, and the thin gray contours represent the phase speed cr, with spacing 0.2. (c), (d) Profiles and stability diagram for the asymmetric stratified shear layer with ε = 0.25. In both (b), (d) the dashed-dotted line represents the resonance condition between vorticity and gravity waves.

Grahic Jump Location
Fig. 17

Profiles and stability diagram for the smooth profiles of Holmboe's [14] symmetric (ε = 0) stratified shear layer with R = 5. All notation the same as in Fig. 16. Oscillations appear on the high-k stability boundary due to issues in numerical resolution (see Ref. [43] for more details).

Grahic Jump Location
Fig. 18

Profiles measured in the Fraser River estuary. (a) Density and N*2, (b) velocity and vorticity gradients, along with the displacement eigenfunction from the most amplified wavenumber of the linear stability predictions in (c). An echosounding image is shown in (d) where instability waves can be seen at the depth where the eigenfunction predicts a maximum displacement (indicated by the horizontal line in all panels). The wavelength of maximum growth rate is shown by the horizontal bar in (d). (Figure is modified from Tedford et al. [45].)



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.

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