Breaking wave field statistics with a multi-layer model

The statistics of breaking wave fields are characterised within a novel multi-layer framework, which generalises the single-layer Saint-Venant system into a multi-layer and non-hydrostatic formulation of the Navier-Stokes equations. We simulate an ensemble of phase-resolved surface wave fields in physical space, where strong nonlinearities, including directional wave breaking and the subsequent highly rotational flow motion, are modelled, without surface overturning. We extract the kinematics of wave breaking by identifying breaking fronts and their speed, for freely evolving wave fields initialised with typical wind wave spectra. The ?(c) distribution, defined as the length of breaking fronts (per unit area) moving with speed c to c + dc following Phillips (J. Fluid Mech., vol. 156, 1985, pp. 505-531), is reported for a broad range of conditions. We recover the ?(c) ? c(-6) scaling without wind forcing for sufficiently steep wave fields. A scaling of ?(c) based solely on the root-mean-square slope and peak wave phase speed is shown to describe the modelled breaking distributions well. The modelled breaking distributions are in good agreement with field measurements and the proposed scaling can be applied successfully to the observational data sets. The present work paves the way for simulations of the turbulent upper ocean directly coupled to a realistic breaking wave dynamics, including Langmuir turbulence, and other sub-mesoscale processes.

Automated summary

The statistics of breaking wave fields are studied using a novel multi-layer framework, which extends the single-layer Saint-Venant system to a multi-layer and non-hydrostatic formulation of the Navier-Stokes equations. Phase-resolved surface wave fields with strong nonlinearities, including directional wave breaking and rotational flow motion, are simulated without surface overturning. The kinematics of wave breaking are extracted by identifying breaking fronts and their speed, and a scaling based on root-mean-square slope and peak wave phase speed is shown to accurately describe the modelled breaking distributions.