Stability analysis of time-periodic shear flow generated by an oscillating density interface

This paper theoretically and numerically investigates the stability of time-periodic shear flows generated by oscillating density interfaces, demonstrating that instability occurs when a nondimensional parameter exceeds a critical threshold, leading to exponential growth of perturbations and the formation of Kelvin-Helmholtz billows, with implications for diapycnal mixing in real-world water bodies like Lake Geneva and Chesapeake Bay.

The Gist

Imagine a large, rectangular swimming pool filled with two layers of water: a lighter, warmer layer on top and a heavier, colder layer on the bottom. Usually, these layers sit quietly on top of each other, like oil on water. But what happens if you gently tilt the pool back and forth, like a seesaw?

This paper investigates exactly that scenario. It asks: If you rock a two-layered fluid back and forth, does the boundary between the layers stay smooth, or does it eventually get chaotic and break apart?

Here is the story of the research, broken down into simple concepts:

1. The Setup: The "Seesaw" Tank

The researchers imagine a tank with two fluids. They tilt the tank slightly and oscillate it (rock it) back and forth.

• The Physics: When the tank tilts, gravity pulls the heavy bottom layer "downhill" and the light top layer "uphill." Because the tank is moving, this creates a shear flow—the top layer slides one way, and the bottom layer slides the other.
• The Twist: Unlike a steady river flow, this shear is time-periodic. It speeds up, slows down, reverses, and changes direction in a rhythmic cycle, just like the tides or the sloshing of a lake during a storm.

2. The Discovery: The "Tunnel" to Chaos

The team found that the boundary between the two layers doesn't become unstable immediately. It's like a car waiting at a red light that turns green only at a specific moment.

• The Waiting Game: At the start of the rocking cycle, the boundary is stable. It wiggles a bit but holds together.
• The Turning Point: As the tank continues to rock, there comes a specific moment (a "turning point") where the physics changes. The stability "tunnels" through a barrier and suddenly becomes unstable.
• The Explosion: Once this threshold is crossed, tiny ripples on the boundary start to grow exponentially. They don't just get bigger; they roll up into giant, swirling clouds known as Kelvin-Helmholtz billows. You've likely seen this in nature: the way clouds roll up in the sky when wind blows over a layer of air, or how cream swirls into coffee.

3. The "Magic Number" ($\beta$)

The researchers developed a "magic number" (called $\beta$) to predict when this chaos will happen. Think of $\beta$ as a measure of how hard you are rocking the tank relative to how heavy the layers are.

• The Rule: If you rock the tank gently (low $\beta$), the layers stay calm forever.
• The Threshold: If you rock it hard enough (specifically, if $\beta$ is greater than 1/4 for equal layers, or slightly less for unequal layers), the layers will eventually break.
• The Correction: The paper includes a "corrigendum" (a correction note). The authors realized they made a small math error when the layers are of unequal depth. They fixed the formula, which slightly changes the threshold for when instability starts in real-world scenarios like lakes, but it doesn't change the main conclusion: rocking the tank causes the layers to mix.

4. How They Solved It

The math behind this is tricky because the forces are constantly changing. The authors used three different tools to understand it:

1. The "Steady" Guess: They tried pretending the tank was just tilted at its maximum angle and not moving. Surprisingly, this simple guess gave them the right answer for when the instability starts, even though it couldn't explain the timing.
2. The "WKB" Method (Modified Airy Function): This is a sophisticated mathematical technique used to track waves through changing environments. It's like using a high-tech GPS to track a car driving through a foggy, winding road. This method perfectly predicted the exact moment the waves would start growing.
3. The "Vortex Blob" Simulation: They built a computer model where they treated the boundary as a string of tiny, invisible spinning tops (vortices). As the tank rocked, these tops interacted, and the simulation showed the boundary rolling up into those famous billow clouds, just like in real life.

5. Real-World Application: Lakes and Bays

The authors didn't just stop at the math; they applied their findings to two real places:

• Lake Geneva: A deep lake in Europe.
• Chesapeake Bay: A large estuary in the US.

In these places, the "tank" is the lake itself, and the "rocking" is caused by tides or wind. The study suggests that even if the water looks calm, the internal waves caused by the tides can create enough shear to trigger these mixing events. This is important because this mixing helps distribute oxygen, nutrients, and heat throughout the water, which is vital for the ecosystem.

Summary

In short, this paper explains that rocking a two-layered fluid creates a rhythmic shear that eventually causes the layers to mix violently. It provides the precise mathematical rules for when this happens, corrects a small error in the math for uneven layers, and shows that this mechanism is likely a key driver of mixing in our oceans and lakes. The boundary between the layers acts like a dam that holds back chaos until the rocking rhythm hits a specific beat, at which point the dam breaks, and the water swirls into beautiful, turbulent clouds.

Technical Summary

Technical Summary: Stability Analysis of Time-Periodic Shear Flow Generated by an Oscillating Density Interface

Problem Statement
This study investigates the stability of a two-layered, inviscid fluid system subjected to time-periodic shear flow generated by the oscillation of a density interface (pycnocline). The physical model is a conceptual two-layered oscillating tank, analogous to low-frequency, basin-scale oscillations (e.g., seiches or internal tides) observed in stratified lakes, estuaries, and shelf seas. Unlike traditional stability analyses that assume steady background shear, this work addresses the specific case where the background shear is explicitly time-dependent and periodic. The primary objective is to understand the alternate pathway to turbulence and diapycnal mixing triggered by the self-induced shear of an oscillating pycnocline, distinct from wind-induced surface shear mechanisms.

Methodology
The analysis proceeds through a multi-stage theoretical and numerical framework:

1. Governing Equations and Linearization: The authors derive the governing equations for a Boussinesq, two-layered fluid in a tank oscillating at a small angle $\alpha(t)$. By introducing infinitesimal perturbations to the background state, they derive a second-order ordinary differential equation (ODE) for the interfacial displacement. This equation is identified as a Schrödinger-type ODE with a periodic potential, characterized by a small parameter $\epsilon = \omega_f / \omega \ll 1$, representing the separation between the low-frequency background oscillation ($\omega_f$) and the high-frequency interfacial waves ($\omega$).

2. Stability Analysis:

   • Necessary and Sufficient Conditions: The paper establishes that instability occurs when the periodic potential $F(\tau)$ becomes negative. This leads to a necessary and sufficient condition for instability involving a non-dimensional parameter $\beta$ (resembling the inverse Richardson number). For equal depth layers, instability arises when $\beta > 1/4$. For unequal depths, the condition is generalized to $\beta > \beta_{min}$, where $\beta_{min}$ depends on the depth ratio $r$.
   • Turning Point Dynamics: The analysis reveals that the flow is initially stable but "tunnels" into an unstable region after reaching a specific turning point time $\tau_c$. This delayed onset is a critical feature distinguishing this unsteady system from steady-state analyses.
   • Asymptotic Solution (MAF): To obtain a uniformly valid composite approximate solution across the turning points, the authors employ the Modified Airy Function (MAF) method. This technique improves upon the standard Wentzel–Kramers–Brillouin (WKB) theory by providing accurate solutions both near and far from the turning points.
   • Floquet and Numerical Analysis: The stability is further verified using Floquet theory and direct numerical integration of the linear ODE. These methods confirm the instability thresholds and growth rates predicted by the asymptotic analysis.

3. Nonlinear Evolution: The study extends the linear analysis to the fully nonlinear regime using a modified vortex blob method. The circulation evolution equation is updated to account for the unsteady background flow. This numerical simulation captures the roll-up of the interface into Kelvin-Helmholtz (KH) billows.

4. Case Studies: The theoretical framework is applied to real-world scenarios: Lake Geneva (summer 1950 conditions) and Chesapeake Bay. These case studies utilize observed stratification and oscillation parameters to predict instability thresholds and cut-off wavelengths.

Key Contributions and Results

• Instability Mechanism: The paper demonstrates that the instability is of the Kelvin-Helmholtz type, driven by the time-periodic shear. It clarifies that parametric resonance (as seen in Mathieu-type equations) is not the driving mechanism due to the large separation of time scales between the forcing and the natural frequency.
• Correction of Generalized Analysis (Corrigendum): A significant technical contribution is the correction of the generalized analysis for arbitrary layer depth ratios ($r \in (0,1)$). The corrected instability condition is $\beta > \beta_{min} = \frac{1}{4}[r^2 + (1-r)^2 + \frac{1}{2}]$. For the specific case of equal depths ($r=1/2$), the condition remains $\beta > 1/4$.
• Quantitative Predictions:
  • For Lake Geneva, the corrected analysis predicts a long-wave cut-off wavelength of $\lambda_{long-wave} \approx 17.03$ m (previously 12.77 m in the draft, corrected in the final text to reflect the specific depth ratio). The instability onset is predicted to occur approximately 40 hours after the inception of mode-1 oscillations.
  • For Chesapeake Bay, the instability condition is corrected to $\beta > 0.248$ (vs. 0.246), with a cut-off wavelength of 4.17 m.
• Validation: The MAF solutions show near-indistinguishable agreement with numerical solutions for the linear amplitude evolution. The nonlinear vortex blob simulations successfully reproduce the formation of KH billows, qualitatively matching Direct Numerical Simulation (DNS) results from previous literature (Inoue & Smyth, 2009; Lewin & Caulfield, 2022).
• Steady vs. Unsteady Comparison: The study highlights that while steady-state analysis using maximum background velocities can predict the correct instability threshold and growth rate magnitude, it fails to capture the critical delayed onset of instability and the specific temporal evolution of the amplitude envelope.

Significance
The paper claims significance in providing a bespoke linear stability analysis for non-autonomous systems where standard steady-state assumptions are insufficient. By identifying the "tunneling" behavior and the specific turning point time $\tau_c$, the study offers a more accurate physical perspective on how low-frequency basin oscillations can trigger high-frequency interfacial instabilities. The results suggest that even in marginally stable stratified environments, the addition of shear from pycnocline oscillations can be sufficient to trigger a cascade of energy into turbulence, thereby enhancing diapycnal mixing. The application to Lake Geneva and Chesapeake Bay provides a physical context, demonstrating the potential relevance of this mechanism in natural aquatic systems. The work bridges the gap between theoretical stability analysis and the observed late-time dynamics of KH billows in oscillating stratified flows.