Waves in a shear flow: transition between the KH, Holmboe and Miles instability

This paper investigates shear-driven wave generation in a two-fluid system with an exponential velocity profile, revealing a novel transition from Kelvin-Helmholtz to Holmboe and finally to Miles critical layer instabilities as the density ratio decreases, a finding validated by both linear theory and nonlinear simulations that demonstrate the first unified observation of all three canonical instabilities within a single background state.

The Gist

Imagine two layers of fluid sitting on top of each other, like oil floating on water, or air blowing over the ocean. The top layer is moving fast, while the bottom layer is sitting still. Because they are moving at different speeds, they rub against each other. This friction creates a "shear" that can make the flat surface between them wobble and eventually break into waves.

This paper is a deep dive into what kind of waves form and why, depending on how heavy the top fluid is compared to the bottom one.

The authors studied a specific setup: a smooth, exponential wind blowing over a sharp boundary between two fluids. They discovered that as you change the density ratio (how heavy the top fluid is relative to the bottom), the waves don't just get bigger; they fundamentally change their personality.

Here is the story of the three distinct "personalities" these waves can take, explained through simple analogies:

1. The "Miles" Wave: The Invisible Thief (Low Density Ratio)

Scenario: Think of wind blowing over the ocean (Air is very light, Water is very heavy). The density ratio is tiny (about 0.001).

• The Analogy: Imagine a pickpocket (the wind) walking past a crowd. The pickpocket doesn't need to push the whole crowd; they just need to brush against one specific person to steal a wallet.
• What happens: The instability happens at a very specific, invisible "critical layer" in the air, just above the water. The wind steals energy only from this thin slice of air.
• The Result: The waves grow slowly but steadily. They look like classic ocean swells. The paper found that even if the top fluid is 10 times heavier than air (but still much lighter than water), this "Miles" behavior persists. It's a very efficient, localized energy theft.

2. The "Holmboe" Wave: The Spiky Cusp (Medium Density Ratio)

Scenario: Think of fresh water meeting salt water, or two different oils mixing. The top fluid is now about half as heavy as the bottom one (Density ratio ~0.5).

• The Analogy: Imagine two people running side-by-side on a track, but one is slightly faster. Instead of a smooth hand-off, they start to trip over each other's feet, creating a jagged, chaotic tangle.
• What happens: The "invisible thief" mechanism breaks down. Now, the wind (or top fluid) interacts with the entire layer of moving fluid, not just a thin slice. The waves become asymmetric.
• The Result: The waves develop a sharp, jagged peak (a "cusp") at the top. Instead of rolling smoothly, the tip of the wave gets sheared off, shooting out tiny droplets (spume). It looks like a wave that has been "sliced" by a knife. This is the Holmboe instability.

3. The "Kelvin-Helmholtz" (KH) Wave: The Rolling Spiral (High Density Ratio)

Scenario: Think of two layers of water (or two very similar liquids) sliding past each other. The top fluid is almost as heavy as the bottom one (Density ratio ~0.9).

• The Analogy: Imagine two sheets of paper sliding past each other. Because they are so similar in weight, they don't just trip; they curl up into tight, perfect spirals.
• What happens: The distinction between the "critical layer" and the rest of the fluid disappears. The whole interface becomes unstable.
• The Result: The waves curl up into the classic, iconic spirals you see in clouds (mackerel sky) or in a cup of coffee when you stir milk in. These are the famous Kelvin-Helmholtz spirals. They distort rapidly and mix the two fluids violently.

The Big Discovery: A Smooth Transition

The most exciting part of this paper is that the authors didn't just find these three types; they found the smooth transition between them.

Usually, scientists study these three as separate, unrelated phenomena. This paper shows that if you start with air over water (Miles) and slowly make the air heavier and heavier (simulating different fluids), the wave smoothly morphs:

1. It starts as a Miles wave (stealing energy from a thin layer).
2. It morphs into a Holmboe wave (developing sharp, spiky crests that spray droplets).
3. It finally becomes a Kelvin-Helmholtz wave (curling into spirals).

Why Does This Matter?

• For the Ocean: It helps us understand how wind creates waves and how much spray (sea foam) is generated, which affects climate models and weather prediction.
• For the Stars: These same physics apply to the boundary layers of stars and accretion disks in space, where gases of different densities swirl around each other.
• For Industry: It helps engineers design better pipelines for oil and gas, or understand how liquids mix in chemical reactors.

The "Secret Sauce" of the Paper

The authors used a mix of math theory (solving complex equations to predict the waves) and computer simulations (watching virtual waves crash and break). They proved that their computer models matched the math perfectly in the beginning, and then showed us the beautiful, chaotic shapes the waves take when they get big and nonlinear.

In a nutshell: The paper tells us that the "personality" of a wave isn't fixed. It depends entirely on how heavy the top fluid is compared to the bottom. Change the weight, and you change the wave from a smooth swell, to a spiky spray, to a curling spiral.

Technical Summary

Here is a detailed technical summary of the paper "Waves in a shear flow: transition between the KH, Holmboe and Miles instability" by Anil Kumar, S. Ravichandran, and Ratul Dasgupta.

1. Problem Statement

The study investigates the generation of waves at the interface of two immiscible fluids driven by shear instabilities. The specific configuration involves a sharp density interface (stable stratification) located within a shear layer of the upper fluid, which follows a classic exponential velocity profile ($U(z) = U_\infty[1 - \exp(-z/\Delta)]$). The lower fluid is initially quiescent.

The central problem is to understand how the nature of the fastest growing instability mode transitions as the density ratio ($\delta = \rho_{upper}/\rho_{lower}$) varies from near unity (e.g., fresh/salt water, $\delta \approx 0.9$) down to the air-water limit ($\delta \approx 0.001$). While previous studies have examined transitions between specific pairs of instabilities (e.g., Kelvin-Helmholtz to Holmboe, or Miles to Kelvin-Helmholtz) for specific profiles, this paper aims to demonstrate a continuous, smooth transition between all three canonical inviscid instabilities (Kelvin-Helmholtz, Holmboe, and Miles) within a single background state model without using the Boussinesq approximation.

2. Methodology

The authors employ a dual approach combining linear stability analysis and nonlinear numerical simulations.

• Linear Stability Analysis:

  • Governing Equations: The inviscid, incompressible Euler equations with gravity and surface tension are linearized. The perturbation streamfunctions satisfy the Rayleigh equation.
  • Analytical Solution: For the exponential velocity profile, the Rayleigh equation is transformed into a Gauss hypergeometric equation, allowing for an analytical expression of the eigenfunctions.
  • Dispersion Relation: A dispersion relation is derived (Eq. 2.12) linking the complex phase speed ($c$) to the wavenumber ($k$), Froude number ($Fr$), Bond number ($Bo$), and density ratio ($\delta$).
  • Comparative Models: The results are benchmarked against:
    1. A Piecewise Linear (PL) velocity profile (to isolate the effects of curvature).
    2. A Discontinuous KH profile (vortex sheet).
    3. Existing literature (Morland & Saffman 1993; Young & Wolfe 2014).

• Nonlinear Numerical Simulations:

  • Solver: The open-source code Basilisk is used to solve the full incompressible Euler equations with surface tension (Volume-of-Fluid method).
  • Initial Conditions: Simulations are initialized with the base state plus the fastest growing eigenmode from the linear theory.
  • Parameters: The study focuses on high Bond numbers ($Bo \sim O(100)$, relevant to geophysical/astrophysical scales) and varying Froude numbers ($Fr \sim O(1)$ to $O(1000)$).

3. Key Contributions

1. Unified Transition Demonstration: The paper provides the first demonstration of a smooth transition between three distinct instability mechanisms (Miles $\to$ Holmboe $\to$ Kelvin-Helmholtz) within a single exponential background profile as the density ratio $\delta$ increases from 0.001 to 0.9.
2. Reynolds Stress Behavior: A novel theoretical explanation is provided for the qualitative change in the Reynolds stress ($\tau$) profile. The authors show that the sharp jump in Reynolds stress characteristic of the Miles instability (at the critical layer) smoothens out as $\delta$ increases, transitioning to a smooth variation characteristic of Holmboe and KH instabilities.
3. Persistence of Miles Instability: The study demonstrates that the Miles critical layer instability persists up to $\delta = 0.01$ (ten times the air-water value), challenging the notion that it is exclusive to extreme density contrasts.
4. Nonlinear Regime Characterization: The authors map the nonlinear saturation of these waves, identifying distinct morphological signatures:
   • Miles/Holmboe regime: Formation of sheared cusps and spume droplets.
   • KH regime: Formation of classic spiral vortices.

4. Key Results

A. Linear Stability Regimes

• Low $\delta$ ($\delta \le 0.01$): Miles Instability

  • The fastest growing mode travels at the surface-gravity wave speed ($c_g$).
  • Critical Layer Dynamics: There is a near-discontinuous jump in the inviscid Reynolds stress at the critical location ($z_c$ where $U(z_c) = c_r$).
  • Energy Extraction: Energy is extracted from the shear primarily in the thin region between the critical layer and the interface.
  • Phase Shift: A characteristic phase jump occurs in the streamfunction across $z_c$.

• Intermediate $\delta$ ($0.1 \lesssim \delta \lesssim 0.5$): Transition to Holmboe

  • As $\delta$ increases, the ratio of growth rate to phase speed ($c_i/c_r$) increases significantly (from $\sim 0.001$ to $\sim 0.4$).
  • The sharp jump in Reynolds stress at $z_c$ disappears, becoming a smooth variation.
  • Energy extraction spreads across the entire shear layer.
  • The growth rates and phase speeds align closely with the Piecewise Linear (PL) model, indicating a resonance between vorticity and density interfaces (Holmboe mechanism).

• High $\delta$ ($\delta \ge 0.9$): Kelvin-Helmholtz (KH) Instability

  • The instability is dominated by the interaction of vorticity interfaces.
  • The mode travels at speeds closer to the mean shear velocity ($U_\infty/2$).
  • The Reynolds stress and kinetic energy profiles show no distinct signature at the critical layer.

B. Nonlinear Simulations

• $\delta = 0.01$ (Miles regime): Waves grow exponentially and match linear theory for several periods. Upon saturation, tiny surface ripples appear.
• $\delta = 0.1$: Finite-amplitude waves exhibit larger ripples reminiscent of capillary-gravity Stokes waves.
• $\delta = 0.5$ (Holmboe regime): Waves develop a sheared cusp at the crest. These cusps emit spume droplets, resembling asymmetric Holmboe waves observed in experiments (Lawrence et al., 1998).
• $\delta = 0.9$ (KH regime): Interfacial waves rapidly distort into classic KH spirals (vortex roll-ups) within five wave periods.

C. Theoretical Explanation of Reynolds Stress Transition

The authors utilize Lin's relation for the vertical variation of Reynolds stress. They show that the term responsible for the jump in stress behaves like a Dirac delta function when the growth rate is very small relative to the phase speed ($c_i/c_r \ll 1$, the Miles case). As $\delta$ increases, $c_i/c_r$ increases, causing the "delta function" to broaden into a smooth distribution, thereby eliminating the sharp jump in stress.

5. Significance

• Geophysical and Astrophysical Relevance: The findings bridge the gap between wind-wave generation (air-water, low $\delta$) and internal wave dynamics in oceans (fresh/salt water, high $\delta$) and astrophysical accretion flows.
• Experimental Implications: The discovery that Miles instability persists at $\delta=0.01$ suggests that experimental measurement of this instability is feasible using fluid pairs with higher density ratios than air-water, where growth rates are larger and easier to detect.
• Modeling Accuracy: The study highlights the critical role of background profile curvature. Comparisons with Piecewise Linear models show that while growth rates may match, the phase speeds and physical mechanisms (critical layer vs. interface resonance) differ significantly, emphasizing the need for smooth profiles in accurate modeling.
• Wave Breaking and Spray: The identification of the "sheared cusp" and spume emission in the Holmboe regime provides a theoretical basis for understanding wave breaking and spray generation in intermediate density ratio scenarios, distinct from the classic KH spiral breaking.

In conclusion, this work unifies the understanding of shear-driven interfacial instabilities, demonstrating that the nature of the instability is not fixed but evolves continuously with the density ratio, governed by the interplay between critical layer dynamics and interfacial wave resonance.