A windy sea surface with Stokes waves

Technical Summary: A Windy Sea Surface with Stokes Waves

Problem Statement
The transition of wind-forced, gravito-capillary (GC) surface waves from smooth to corrugated states remains a longstanding challenge in nonlinear wave mechanics. While the theoretical threshold for wind-generated waves was established by Kelvin 155 years ago, the specific mechanisms governing the onset of surface roughness (corrugations) on short, progressive waves (4–13 cm) are not fully resolved. Existing models often rely on empirical estimates or numerical simulations initialized without wind-forcing, lacking first-principles theoretical results that account for the balance between wind energy input and viscous dissipation in maintaining permanent waveforms. Furthermore, previous regime maps distinguishing smooth from corrugated waves often assumed a sharp boundary, potentially overlooking the complexity of the transition region.

Methodology
The authors approach the problem by modeling the wind-driven sea surface as a collection of nonlinear, permanent-form Stokes waves. The study employs a mathematical framework based on driven-dissipative, nonlinear potential flow equations, incorporating viscous dissipation and wind forcing.

• Linear Analysis: The authors begin with a linearized analysis of a steady, sinusoidal waveform in a co-moving frame. They utilize the Helmholtz decomposition for perturbation velocity and apply Miles' form for wind perturbation pressure. This yields a viscous dispersion relation linking the Froude number ($F_\lambda$), wind stress ($P_\lambda$), Bond number ($\alpha_\lambda$), and Reynolds number ($Re_\lambda$).
• Nonlinear Analysis: The study extends to finite-amplitude waves by solving the full nonlinear potential flow equations in a co-moving frame. The formulation includes a vortical contribution to the vertical velocity component and utilizes a dynamic boundary condition derived from the nonlinear Bernoulli equation with a linearized viscous contribution.
• Dimensionality Reduction: A critical methodological constraint is the fixation of the Morton number ($M_o$) to the specific value for air-water interfaces ($2.29 \times 10^{-3}$). This reduces the dimensionality of the parameter space, allowing the Bond number to be expressed in terms of $Re_\lambda$ and $F_\lambda$.
• Stability Analysis: The authors perform a subharmonic instability analysis on the computed waveforms to determine how linear instabilities affect the evolution of corrugations.

Key Contributions and Results

1. Regime Maps in $(Re_\lambda, E)$ Space:
   The primary result is the construction of regime maps in the space of the wave-based Reynolds number ($Re_\lambda$) and wave energy ($E$). These maps distinguish between:

   • Region I: Smooth nonlinear carrier waves.
   • Region II: Waves with corrugations (ripples) predominantly on the leeward (downwind) face.
   • Region IV: Waves with corrugations on both windward and leeward faces.
   • Region III: A transition band linking Regions II and IV.

2. The "Band-Like" Transition:
   Contrary to previous studies that suggested a sharp boundary between smooth and corrugated states, the authors identify a transition band of finite thickness. Within this band, the wave steepness ($\epsilon$) exhibits a non-monotonic dependence on the Reynolds number. This implies that the transition is not a single critical point but a region where multiple wave states (smooth vs. corrugated) can coexist or where the steepness varies complexly.

3. Revised Critical Conditions:
   The study significantly revises earlier estimates (e.g., from [8, 21]) regarding the onset of corrugations.

   • The critical Bond number ($\alpha_\lambda^{-1}$) beyond which only smooth waves appear is found to be $>2500$ at $Re_\lambda \approx 6 \times 10^4$.
   • Unlike fixed empirical estimates of steepness cutoffs, the critical steepness in this study varies significantly with $Re_\lambda$.
   • The inclusion of wind stress, a crucial ingredient absent in prior DNS studies initialized without wind, is shown to be essential for accurately identifying these regimes.

4. Symmetry Trends:
   Within the corrugated regime, decreasing the Reynolds number transitions asymmetric, wind-forced GC waves toward symmetric waveforms. This trend is noted as opposite to observations in previous literature (e.g., [31]), which the authors attribute to the specific constraint of fixing the Morton number for air-water interfaces in their calculations.

5. Subharmonic Instability:
   Linear stability analysis of Region II (one-sided ripples) reveals that the fastest-growing subharmonic mode intensifies corrugations on alternate faces of the carrier wave while smoothing others. This behavior is reminiscent of the Benjamin-Feir instability but occurs in the presence of wind forcing.

Significance and Claims
The paper claims to provide the first first-principles regime maps for wind-forced GC waves that explicitly account for the balance of energy input and viscous dissipation. By reducing the dimensionality of the parameter space through the air-water constraint, the authors offer a more precise delineation of smooth versus corrugated states than previously available.

The findings have broad implications for:

• Oceanic Remote Sensing: Predicting the probability and nature of bound-wave components (corrugations) on carrier waves is crucial for interpreting remote-sensed ocean images. The study provides a theoretical basis for estimating these probabilities without relying solely on fixed empirical cutoffs.
• Nonlinear Wave Mechanics: The identification of a finite-thickness transition band and the non-monotonic behavior of wave steepness challenge the notion of a sharp phase transition in wave mechanics.
• Interfacial Pattern Formation: The insights into how wind stress and viscosity interact to generate specific ripple patterns on wave faces contribute to the fundamental understanding of interfacial dynamics.

The authors maintain a modest tone, noting that while their results revise early estimates and offer new regime distinctions, further work (including long-term evolution and two-sided corrugation stability) will be presented in forthcoming publications.