An experimental study of the dynamics and droplet production in three mechanically generated plunging breaking waves is presented in this two-part paper. In the present paper (Part 2), in-line cinematic holography is used to measure the positions, diameters (
Droplet dynamics under an impinging air jet
Partially wetting droplets under an airflow can exhibit complex behaviours that arise from the coupling of surface tension, inertia of the external flow and contact-line dynamics. Recent experiments by Hooshanginejad et al. ( J. Fluid Mech. , vol. 901, 2020) revealed that a millimetric partially wetting water droplet under an impinging jet can oscillate in place, split or depin away from the jet, depending on the magnitude (i.e. $5\unicode{x2013}20\ {\rm m}\ {\rm s}^{-1}$ ) and position of the jet. To rationalise the experimental observations, we develop a two-dimensional lubrication model of the droplet that incorporates the external pressure of the impinging high-Reynolds-number jet, in addition to the capillary and hydrostatic pressures of the droplet. Distinct from the previous model by Hooshanginejad et al. ( J. Fluid Mech. , vol. 901, 2020), we simulate the motion of the contact line using precursor film and disjoining pressure, which allows us to capture a wider range of droplet behaviours, including the droplet dislodging to one side. Our simulations exhibit a comparable time-scale of droplet deformations and similar outcomes as the experimental observations. We also obtain the analytical steady-state solutions of the droplet shapes and construct the minimum criteria for splitting and depinning.
more »
« less
- Award ID(s):
- 2042194
- PAR ID:
- 10403990
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 943
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
), times and velocities of droplets generated by the three plunging breaking waves studied in Part 1 (Erinin$d\geq 100\ \mathrm {\mu }{\rm m}$ et al. ,J. Fluid Mech. , vol. 967, 2023, A35) as the droplets move up across a horizontal measurement plane located just above the wave crests. It is found that there are four major mechanisms for droplet production: closure of the indentation between the top surface of the plunging jet and the splash that it creates, the bursting of large bubbles that were entrapped under the plunging jet at impact, splashing and bubble bursting in the turbulent zone on the front face of the wave and the bursting of small bubbles that reach the water surface at the crest of the non-breaking wave following the breaker. The droplet diameter distributions for the entire droplet set for each breaker are fitted with power-law functions in separate small- and large-diameter regions. The droplet diameter where these power-law functions cross increases monotonically from 820 to 1480 from the weak to the strong breaker, respectively. The droplet diameter and velocity characteristics and the number of the droplets generated by the four mechanisms are found to vary significantly and the processes that create these differences are discussed.$\mathrm {\mu }{\rm m}$ -
Fritts et al. (J. Fluid Mech., vol. xx, 2022, xx) describe a direct numerical simulation of interacting Kelvin–Helmholtz instability (KHI) billows arising due to initial billow cores that exhibit variable phases along their axes. Such KHI exhibit strong ‘tube and knot’ dynamics identified in early laboratory studies by Thorpe ( Geophys. Astrophys. Fluid Dyn. , vol. 34, 1985, pp. 175–199). Thorpe ( Q.J.R. Meteorol. Soc. , vol. 128, 2002, pp. 1529–1542) noted that these dynamics may be prevalent in the atmosphere, and they were recently identified in atmospheric observations at high altitudes. Tube and knot dynamics were found by Fritts et al. ( J. Fluid. Mech. , 2022) to drive stronger and faster turbulence transitions than secondary instabilities of individual KH billows. Results presented here reveal that KHI tube and knot dynamics also yield energy dissipation rates $\sim$ 2–4 times larger as turbulence arises and that remain $\sim$ 2–3 times larger to later stages of the flow evolution, compared with those of secondary convective instabilities (CI) and secondary KHI accompanying KH billows without tube and knot influences. Elevated energy dissipation rates occur due to turbulence transitions by tube and knot dynamics arising on much larger scales than secondary CI and KHI where initial KH billows are misaligned. Tube and knot dynamics also excite large-scale Kelvin ‘twist waves’ that cause vortex tube and billow core fragmentation, more energetic cascades of similar interactions to smaller scales and account for the strongest energy dissipation events accompanying such KH billow evolutions.more » « less
-
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
-
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
-
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