skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Predicting the breaking strength of gravity water waves in deep and intermediate depth
We revisit the classical but as yet unresolved problem of predicting the strength of breaking 2-D and 3-D gravity water waves, as quantified by the amount of wave energy dissipated per breaking event. Following Duncan ( J. Fluid Mech. , vol. 126, 1983, pp. 507–520), the wave energy dissipation rate per unit length of breaking crest may be related to the fifth moment of the wave speed and the non-dimensional breaking strength parameter  $$b$$ . We use a finite-volume Navier–Stokes solver with large-eddy simulation resolution and volume-of-fluid surface reconstruction (Derakhti & Kirby, J. Fluid Mech. , vol. 761, 2014 a , pp. 464–506; J. Fluid Mech. , vol. 790, 2016, pp. 553–581) to simulate nonlinear wave evolution, with a strong focus on breaking onset and postbreaking behaviour for representative cases of wave packets with breaking due to dispersive focusing and modulational instability. The present study uses these results to investigate the relationship between the breaking strength parameter $$b$$ and the breaking onset parameter $$B$$ proposed recently by Barthelemy et al. ( J. Fluid Mech. , vol. 841, 2018, pp. 463–488). The latter, formed from the local energy flux normalized by the local energy density and the local crest speed, simplifies, on the wave surface, to the ratio of fluid speed to crest speed. Following a wave crest, when $$B$$ exceeds a generic threshold value at the wave crest (Barthelemy et al. 2018), breaking is imminent. We find a robust relationship between the breaking strength parameter $$b$$ and the rate of change of breaking onset parameter $$\text{d}B/\text{d}t$$ at the wave crest, as it transitions through the generic breaking onset threshold ( $$B\sim 0.85$$ ), scaled by the local period of the breaking wave. This result significantly refines previous efforts to express $$b$$ in terms of a wave packet steepness parameter, which is difficult to define robustly and which does not provide a generically accurate forecast of the energy dissipated by breaking.  more » « less
Award ID(s):
1756355
PAR ID:
10136849
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
848
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. While it has long been recognized that Lagrangian drift at the ocean surface plays a critical role in the kinematics and dynamics of upper ocean processes, only recently has the contribution of wave breaking to this drift begun to be investigated through direct numerical simulations (Deike et al. ,  J. Fluid Mech. , vol. 829, 2017, pp. 364–391; Pizzo et al. ,  J. Phys. Oceanogr. , vol. 49(4), 2019, pp. 983–992). In this work, laboratory measurements of the surface Lagrangian transport due to focusing deep-water non-breaking and breaking waves are presented. It is found that wave breaking greatly enhances mass transport, compared to non-breaking focusing wave packets. These results are in agreement with the direct numerical simulations of Deike  et al. ( J. Fluid Mech. , vol. 829, 2017, pp. 364–391), and the increased transport due to breaking agrees with their scaling argument. In particular, the transport at the surface scales with $$S$$ , the linear prediction of the maximum slope at focusing, while the surface transport due to non-breaking waves scales with $$S^{2}$$ , in agreement with the classical Stokes prediction. 
    more » « less
  2. null (Ed.)
    In this work, we study the nonlinear travelling waves in density stratified fluids with piecewise-linear shear currents. Beginning with the formulation of the water-wave problem due to Ablowitz et al.  ( J. Fluid Mech. , vol. 562, 2006, pp. 313–343), we extend the work of Ashton & Fokas ( J. Fluid Mech. , vol. 689, 2011, pp. 129–148) and Haut & Ablowitz ( J. Fluid Mech. , vol. 631, 2009, pp. 375–396) to examine the interface between two fluids of differing densities and varying linear shear. We derive a systems of equations depending only on variables at the interface, and numerically solve for periodic travelling wave solutions using numerical continuation. Here, we consider only branches which bifurcate from solutions where there is no slip in the tangential velocity at the interface for the trivial flow. The spectral stability of these solutions is then determined using a numerical Fourier–Floquet technique. We find that the strength of the linear shear in each fluid impacts the stability of the corresponding travelling wave solutions. Specifically, opposing shears may amplify or suppress instabilities. 
    more » « less
  3. We present high-resolution three-dimensional (3-D) direct numerical simulations of breaking waves solving for the two-phase Navier–Stokes equations. We investigate the role of the Reynolds number ( Re , wave inertia relative to viscous effects) and Bond number ( Bo , wave scale over the capillary length) on the energy, bubble and droplet statistics of strong plunging breakers. We explore the asymptotic regimes at high Re and Bo , and compare with laboratory breaking waves. Energetically, the breaking wave transitions from laminar to 3-D turbulent flow on a time scale that depends on the turbulent Re up to a limiting value $$Re_\lambda \sim 100$$ , consistent with the mixing transition in other canonical turbulent flows. We characterize the role of capillary effects on the impacting jet and ingested main cavity shape and subsequent fragmentation process, and extend the buoyant-energetic scaling from Deike et al. ( J. Fluid Mech. , vol. 801, 2016, pp. 91–129) to account for the cavity shape and its scale separation from the Hinze scale, $$r_H$$ . We confirm two regimes in the bubble size distribution, $$N(r/r_H)\propto (r/r_H)^{-10/3}$$ for $$r>r_H$$ , and $$\propto (r/r_H)^{-3/2}$$ for $$r 
    more » « less
  4. When liquid metal batteries are charged or discharged, strong electrical currents are passing through the three liquid layers that we find in their interior. This may result in the metal pad roll instability that drives gravity waves on the interfaces between the layers. In this paper, we investigate theoretically metal pad roll instability in idealised cylindrical liquid metal batteries that were simulated previously by Weber et al. ( Phys. Fluids , vol. 29, no. 5, 2017 b , 054101) and Horstmann et al. ( J. Fluid Mech. , vol. 845, 2018, pp. 1–35). Near the instability threshold, we expect weakly destabilised gravity waves, and in this parameter regime, we can use perturbation methods to find explicit formulas for the growth rate of all possible waves. This perturbative approach also allows us to include dissipative effects, hence we can locate the instability threshold with good precision. We show that our theoretical growth rates are in quantitative agreement with previous and new direct numerical simulations. We explain how our theory can be used to estimate a lower bound on cell size beneath which metal pad roll instability is unlikely. 
    more » « less
  5. The energetics of drop deposition are considered in the capillary-ballistic regime characterized by high Reynolds number and moderate Weber number. Experiments are performed impacting water/glycol drops onto substrates with varying wettability and contact-angle hysteresis. The impacting event is decomposed into three regimes: (i) pre-impact, (ii) inertial spreading and (iii) post contact-line (CL) pinning, conveniently framed using the theory of Dussan & Davis ( J. Fluid Mech. , vol. 173, 1986, pp. 115–130). During fast-time-scale inertial spreading, the only form of dissipation is CL dissipation ( $$\mathcal {D}_{CL}$$ ). High-speed imaging is used to resolve the stick-slip dynamics of the CL with $$\mathcal {D}_{CL}$$ measured directly from experiment using the $$\Delta \alpha$$ - $$R$$ cyclic diagram of Xia & Steen ( J. Fluid Mech. , vol. 841, 2018, pp. 767–783), representing the contact-angle deviation against the CL radius. Energy loss occurs on slip legs, and this observation is used to derive a closed-form expression for the kinetic K and interfacial $$\mathcal{A}$$ post-pinning energy $$\{K+\mathcal {A}\}_p/\mathcal {A}_o$$ independent of viscosity, only depending on the rest angle $$\alpha _p$$ , equilibrium angle $$\bar {\alpha }$$ and hysteresis $$\Delta \alpha$$ , which agrees well with experimental observation over a large range of parameters, and can be used to evaluate contact-line dissipation during inertial spreading. The post-pinning energy is found to be independent of the pre-impact energy, and it is broken into modal components with corresponding energy partitioning approximately constant for low-hysteresis surfaces with fixed pinning angle $$\alpha _p$$ . During slow-time-scale post-pinning, the liquid/gas ( $lg$ ) interface is found to vibrate with the frequencies and mode shapes predicted by Bostwick & Steen ( J. Fluid Mech. , vol. 760, 2014, pp. 5–38), irrespective of the pre-impact energy. Resonant mode decay rates are determined experimentally from fast Fourier transforms of the interface dynamics. 
    more » « less