Bifurcation analysis of Stokes waves with piecewise smooth vorticity in deep water

Technical Summary: Bifurcation Analysis of Stokes Waves with Piecewise Smooth Vorticity in Deep Water

Problem Statement
This paper addresses the rigorous existence of two-dimensional, periodic, steady water waves (Stokes waves) propagating in an infinitely deep fluid domain under gravity. The study focuses on rotational flows where the vorticity distribution is piecewise smooth, specifically allowing for a finite number of jump discontinuities within the fluid. This setting models stratification effects observed in real oceanographic phenomena, such as sudden changes in current shear.

The mathematical formulation involves a free boundary problem where the fluid domain is unbounded below ($y \to -\infty$) and the free surface ($y = \eta(x)$) is unknown. The governing equations are the steady Euler equations for an inviscid, incompressible fluid. The primary challenges identified are:

1. The presence of an internal interface due to vorticity discontinuities.
2. The unbounded spatial domain, which prevents the linearized operator from being Fredholm (a standard requirement for classical bifurcation theory).
3. The need to handle weak solutions due to the discontinuous vorticity, requiring a reformulation that accommodates transmission conditions across the interface.

Methodology
The authors employ a multi-step analytical approach combining hodograph transformations, height function formulations, and singular bifurcation theory:

1. Reformulation via Height Function:
   The problem is transformed from the physical $(x, y)$ domain to a fixed rectangular domain using the Dubreil-Jacotin transformation (hodograph transformation). The independent variables become $(q, p)$, where $q=x$ and $p=-\psi$ (pseudo-stream function). The dependent variable is the height function $h(q, p) = y$. This converts the free boundary problem into a quasilinear elliptic PDE with oblique boundary conditions on a fixed domain. The piecewise smooth vorticity $\gamma(\psi)$ introduces a transmission problem with specific interface conditions at the depth where the vorticity jumps.

2. Approximation and Regularization:
   Due to the unbounded domain, the linearized operator is not Fredholm. To overcome this, the authors introduce a sequence of approximating problems (parameterized by $\epsilon$) where the operator is modified to be Fredholm of index zero. This follows the strategy used in previous works on deep water waves but adapted for the transmission problem structure.

3. Global Bifurcation Theory:
   The existence of non-trivial solutions is established using analytic global bifurcation theory (specifically the framework of Buffoni and Toland).

   • Local Bifurcation: The authors identify a bifurcation point from the trivial (laminar) flow solution by analyzing the kernel of the linearized operator. They prove the existence of a simple eigenvalue and verify the transversality condition.
   • Nodal Properties: A crucial step involves preserving the "nodal pattern" (monotonicity of the wave profile) along the bifurcation branch. The authors define a set of solutions satisfying specific sign conditions on derivatives and prove that the local bifurcation curve stays within this set, preventing the formation of closed loops or return to the laminar branch.
   • Global Continuation: Using Whyburn's topological lemma, the authors pass to the limit as the approximation parameter $\epsilon \to 0$. They verify that the set of solutions for the approximating problems is compact and closed, allowing the construction of a global continuum of solutions for the original problem.

4. Estimates and Compactness:
   To handle the unbounded domain in the limit, the authors establish uniform exponential decay estimates for the derivatives of the height function as $p \to -\infty$. These estimates, derived using Phragmén-Lindelöf type arguments and Schauder-type estimates tailored for transmission problems (citing Ladyzhenskaya and Gilbarg-Trudinger), ensure the relative compactness of the solution set in the appropriate function spaces.

Key Results
The main result is Theorem 2.1, which establishes the existence of a connected global bifurcation branch $\mathcal{K}$ of solutions $(c, h)$ (wave speed and height function) in the space $\mathbb{R} \times C^{0,\alpha}(\mathcal{D}) \cap C^{3,\alpha}$ (piecewise smooth regularity).

The continuum $\mathcal{K}$ originates from a laminar flow solution and extends globally. The paper proves that along this branch, one of the following two alternatives must occur:

• (a) Unbounded Wave Speed: The wave speed $c$ tends to infinity ($\lim c_k = \infty$).
• (b) Horizontal Stagnation: The wave profile approaches a state of horizontal stagnation, where the horizontal fluid velocity $u$ approaches the wave speed $c$ ($\lim \sup \partial_p h_k = \infty$, implying $u \to c$). This corresponds to the formation of a stagnation point on the free surface or within the flow.

The results hold for vorticity functions $\gamma$ that are piecewise $C^{1,\alpha}$ with a finite number of jumps, provided the vorticity decays sufficiently fast at infinity ($O(s^{-2-r})$) and satisfies a specific bound related to gravity ($-\Gamma_{inf} < g/2^{3/4}$).

Significance and Claims
The paper claims to fill a specific gap in the mathematical theory of water waves by addressing two difficult aspects simultaneously:

1. Deep Water: Unlike many existing results restricted to finite depth, this work handles the infinite depth case where the velocity field and vorticity decay at the bottom.
2. Discontinuous Vorticity: While previous studies (e.g., Constantin, Hur) have treated smooth vorticity or discontinuous vorticity in bounded domains, this paper provides a comprehensive treatment for piecewise smooth vorticity in deep water.

The authors emphasize that their approach overcomes the non-Fredholm nature of the linearized operator in unbounded domains and the technical challenges of weak solutions arising from vorticity jumps. By combining the height function formulation with a singular bifurcation argument and Whyburn's lemma, they rigorously construct large-amplitude waves that exhibit stratification effects, a feature noted by Constantin as a promising direction for future research in water wave theory. The work extends the applicability of global bifurcation theory to transmission problems in unbounded domains.