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
Drop impact on solids: contact-angle hysteresis filters impact energy into modal vibrations
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
- Award ID(s):
- 1935590
- NSF-PAR ID:
- 10311622
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 923
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We perform direct numerical simulations of a gas bubble dissolving in a surrounding liquid. The bubble volume is reduced due to dissolution of the gas, with the numerical implementation of an immersed boundary method, coupling the gas diffusion and the Navier–Stokes equations. The methods are validated against planar and spherical geometries’ analytical moving boundary problems, including the classic Epstein–Plesset problem. Considering a bubble rising in a quiescent liquid, we show that the mass transfer coefficient $k_L$ can be described by the classic Levich formula $k_L = (2/\sqrt {{\rm \pi} })\sqrt {\mathscr {D}_l\,U(t)/d(t)}$ , with $d(t)$ and $U(t)$ the time-varying bubble size and rise velocity, and $\mathscr {D}_l$ the gas diffusivity in the liquid. Next, we investigate the dissolution and gas transfer of a bubble in homogeneous and isotropic turbulence flow, extending Farsoiya et al. ( J. Fluid Mech. , vol. 920, 2021, A34). We show that with a bubble size initially within the turbulent inertial subrange, the mass transfer coefficient in turbulence $k_L$ is controlled by the smallest scales of the flow, the Kolmogorov $\eta$ and Batchelor $\eta _B$ microscales, and is independent of the bubble size. This leads to the non-dimensional transfer rate ${Sh}=k_L L^\star /\mathscr {D}_l$ scaling as ${Sh}/{Sc}^{1/2} \propto {Re}^{3/4}$ , where ${Re}$ is the macroscale Reynolds number ${Re} = u_{rms}L^\star /\nu _l$ , with $u_{rms}$ the velocity fluctuations, $L^*$ the integral length scale, $\nu _l$ the liquid viscosity, and ${Sc}=\nu _l/\mathscr {D}_l$ the Schmidt number. This scaling can be expressed in terms of the turbulence dissipation rate $\epsilon$ as ${k_L}\propto {Sc}^{-1/2} (\epsilon \nu _l)^{1/4}$ , in agreement with the model proposed by Lamont & Scott ( AIChE J. , vol. 16, issue 4, 1970, pp. 513–519) and corresponding to the high $Re$ regime from Theofanous et al. ( Intl J. Heat Mass Transfer , vol. 19, issue 6, 1976, pp. 613–624).more » « less
-
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
-
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 $rThe YBJ equation (Young & Ben Jelloul, J. Marine Res. , vol. 55, 1997, pp. 735–766) provides a phase-averaged description of the propagation of near-inertial waves (NIWs) through a geostrophic flow. YBJ is obtained via an asymptotic expansion based on the limit $\mathit{Bu}\rightarrow 0$ , where $\mathit{Bu}$ is the Burger number of the NIWs. Here we develop an improved version, the YBJ + equation. In common with an earlier improvement proposed by Thomas, Smith & Bühler ( J. Fluid Mech. , vol. 817, 2017, pp. 406–438), YBJ + has a dispersion relation that is second-order accurate in $\mathit{Bu}$ . (YBJ is first-order accurate.) Thus both improvements have the same formal justification. But the dispersion relation of YBJ + is a Padé approximant to the exact dispersion relation and with $\mathit{Bu}$ of order unity this is significantly more accurate than the power-series approximation of Thomas et al. (2017). Moreover, in the limit of high horizontal wavenumber $k\rightarrow \infty$ , the wave frequency of YBJ + asymptotes to twice the inertial frequency $2f$ . This enables solution of YBJ + with explicit time-stepping schemes using a time step determined by stable integration of oscillations with frequency $2f$ . Other phase-averaged equations have dispersion relations with frequency increasing as $k^{2}$ (YBJ) or $k^{4}$ (Thomas et al. 2017): in these cases stable integration with an explicit scheme becomes impractical with increasing horizontal resolution. The YBJ + equation is tested by comparing its numerical solutions with those of the Boussinesq and YBJ equations. In virtually all cases, YBJ + is more accurate than YBJ. The error, however, does not go rapidly to zero as the Burger number characterizing the initial condition is reduced: advection and refraction by geostrophic eddies reduces in the initial length scale of NIWs so that $\mathit{Bu}$ increases with time. This increase, if unchecked, would destroy the approximation. We show, however, that dispersion limits the damage by confining most of the wave energy to low $\mathit{Bu}$ . In other words, advection and refraction by geostrophic flows does not result in a strong transfer of initially near-inertial energy out of the near-inertial frequency band.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2021.547
Warning: Leaving National Science Foundation Website
You are now leaving the National Science Foundation website to go to a non-government website.
Website:
NSF takes no responsibility for and exercises no control over the views expressed or the accuracy of the information contained on this site. Also be aware that NSF's privacy policy does not apply to this site.
