Two-layer sharply stratified Euler fluids in three dimensions: a Hamiltonian setting

This paper derives a Hamiltonian structure for three-dimensional, two-layer stratified Euler fluids through a reduction process, demonstrating that the resulting effective 2D models correspond to the Kaup-Broer-Kupershmidt-Boussinesq and Kadomtsev-Petviashvili equations under specific approximations.

The Gist

Imagine you are looking at a giant, clear tank of water. Now, imagine that instead of just plain water, the tank is filled with two different liquids that don't mix—like oil sitting on top of water. This is a "two-layer stratified fluid."

If you poke the boundary where the two liquids meet, you create waves. These aren't just surface ripples; they are "internal waves" that travel deep inside the liquid. Scientists want to predict exactly how these waves move, how they grow, and how they eventually die out.

This paper is essentially a high-level "instruction manual" for a mathematical engine that predicts those waves. Here is the breakdown of how they did it:

1. The "Master Blueprint" (The Hamiltonian Setting)

In physics, there is a concept called a Hamiltonian. Think of it as the "DNA" or the "Master Blueprint" of a system. If you know the Hamiltonian, you know everything: how much energy is in the system, how much it wants to move, and how it will react to being pushed.

The authors started with a massive, incredibly complex 3D blueprint that describes every single molecule of fluid in the tank. However, that blueprint is too heavy to use—it’s like trying to fly a Boeing 747 just to go to the grocery store. It’s "overkill."

2. The "Simplification Trick" (Hamiltonian Reduction)

The researchers performed what they call a "reduction." Instead of tracking every drop of liquid in the whole tank, they realized they only really need to track what is happening at the interface—the "border patrol" line where the two liquids meet.

They mathematically "squashed" the 3D complexity down into a 2D model. It’s like instead of tracking every single person in a crowded stadium, you only track the movement of the people sitting on the very front row. If you know what the front row is doing, you can get a very good idea of the energy of the whole crowd.

3. The "Zoom Lens" (Asymptotic Expansion)

The paper then uses a technique called "asymptotics." This is like using a zoom lens on a camera.

• The Wide Shot: They look at the big, slow, long waves (the "Long Wave Limit").
• The Close-Up: They look at the tiny, subtle details of how the waves bend and wiggle (the "Weakly Non-Linear" part).

By zooming in on these specific scales, they derived two famous mathematical "characters":

• The KBK-B Model: This is a complex, 2D model that describes how waves dance and interact in a wide area.
• The KP Equation: This is a "specialized" version. Imagine the waves are mostly traveling in one direction (like a train on a track) but occasionally drifting slightly to the side. The KP equation is the perfect mathematical tool for that specific "one-way" motion.

4. Why does this matter? (The "Hardware" of the Ocean)

The most interesting part of their discovery is that the behavior of the waves depends on the "hardware" of the system—the thickness of the layers and how heavy they are.

They found a "tipping point." Depending on the density and depth of the liquids, the waves can change their personality:

• They can be "Bright" solitons: Waves that look like a sudden, tall hump of energy.
• They can be "Dark" solitons: Waves that look like a sudden, deep dip or a "hole" moving through the liquid.

Summary: The Big Picture

If the ocean or the atmosphere is a giant, complex musical instrument, this paper is providing the sheet music. By using advanced geometry and "mathematical squashing," the authors have created a way to write down the rules for how internal waves move, allowing us to predict the "rhythm" of the fluids without needing to track every single atom.

Technical Summary

This paper, titled "Two-layer sharply stratified Euler fluids in three dimensions: a Hamiltonian setting," presents a rigorous mathematical study of the dynamics of internal waves in a three-dimensional, two-layer incompressible Euler fluid. The authors employ the framework of Hamiltonian geometry to derive simplified, physically relevant models from the fundamental governing equations.

1. The Problem

The study focuses on a 3D fluid system consisting of two homogeneous layers of different densities ($\rho_1$ and $\rho_2$) separated by a sharp interface $z = \zeta(x, y, t)$, confined between two infinite horizontal plates.

While 2D models of such systems are well-documented, the 3D case is mathematically more complex due to the increased degrees of freedom and the potential for non-local effects. The primary challenge addressed is how to systematically reduce the full 3D Euler equations—which are difficult to solve directly—into simplified, low-dimensional models (like the Boussinesq or KP equations) while preserving the underlying Hamiltonian structure and its associated conservation laws.

2. Methodology

The authors utilize a multi-step geometric and asymptotic approach:

• Hamiltonian Reduction: They start with the Benjamin-Bowman Hamiltonian structure for 3D stratified fluids. Using the Marsden-Ratiu (MR) reduction theorem, they project the infinite-dimensional 3D Poisson structure onto a submanifold representing the two-layer configuration. This reduces the phase space to a triple of effective variables: the interface displacement ($\zeta$) and the components of the tangential momentum jump (vorticity sheet) across the interface ($e_\sigma, e_\tau$).
• Long-Wave Asymptotics: To derive tractable models, they apply a long-wave expansion using two small parameters: $\epsilon$ (the ratio of vertical to horizontal scales, representing dispersion) and $\alpha$ (the ratio of wave amplitude to depth, representing nonlinearity).
• Weakly Non-Linear (WNL) Approximation: By linking these parameters ($\alpha \sim \epsilon^2$), they obtain a local Hamiltonian density, avoiding the non-localities that typically plague 3D fluid models.
• Dirac Reduction: To derive the unidirectional model (the KP equation), they use a Dirac-type reduction process. This allows them to impose constraints (such as irrotationality and unidirectional propagation) directly onto the Hamiltonian structure.

3. Key Contributions and Results

A. Derivation of the 2D KBK-B Model

The authors derive the two-dimensional Kaup-Broer-Kupershmidt Boussinesq (KBK-B) model. This model governs the weakly nonlinear dispersive waves at the interface.

• Critical Parameters: They identify a "critical" parameter $B$ (dependent on layer densities and depths). A significant finding is that the sign of $B$ determines whether waves are of elevation (bright solitons) or depression (dark solitons).
• Stationary Solutions: They show that for stationary solutions, the sign of $B$ dictates whether the interface is attracted to the upper or lower plate, a phenomenon they note as a potential area for further study.

B. Derivation of the KP Equation

By breaking the rotational symmetry and assuming slow variation in one horizontal direction, the authors derive the Kadomtsev-Petviashvili (KP-II) equation.

• Hamiltonian Consistency: They prove that the KP equation's bi-Hamiltonian structure is a direct consequence of the parent 3D Euler system through the Dirac reduction of the KBK-B model.
• Soliton Solutions: They provide explicit formulas for "line" traveling solitary waves (line solitons) and establish the physical constraints (speed and amplitude limits) required to ensure the wave does not contact the confining plates.

C. Well-posedness and Regularization

The authors address the mathematical stability of the models. They demonstrate that the linearized WNL system is initially ill-posed, but can be regularized through a specific change of variables (moving from interfacial velocities to layer-averaged momentum densities), which ensures all Fourier modes are stable.

4. Significance

The significance of this work lies in its mathematical elegance and physical consistency:

1. Structural Integrity: Unlike many asymptotic derivations that treat equations as mere approximations, this paper ensures that the simplified models (KBK-B and KP) inherit the Hamiltonian nature of the original Euler equations. This guarantees that the conservation of mass, momentum, and energy is mathematically "built-in."
2. Unified Framework: It provides a unified geometric pipeline—from 3D Euler $\to$ 2D KBK-B $\to$ 1D/Unidirectional KP—using a consistent set of reduction techniques.
3. Physical Insight: The study clarifies how "hardware parameters" (density and depth ratios) fundamentally change the qualitative behavior of the fluid (e.g., the transition between elevation and depression waves), providing a bridge between abstract Hamiltonian geometry and geophysical fluid dynamics.