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: 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
PAR ID:
10311622
Author(s) / Creator(s):
; ;
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
  1. 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
  2. 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
  3. 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
  4. During solid surface impact, a falling drop's energy is transformed into oscillations of its liquid/gas interface. We consider drop deposition during oblique impact in the capillary-ballistic regime characterized by high Reynolds number and moderate Weber number. We treat this as an inverse problem showing that post-impact observations of the frequency spectrum and modal partition of energy allow one to determine a drop's pre-impact characteristics and wetting properties. Our analysis is useful for quantifying contact-line dissipation during inertial spreading and can be used as a diagnostic technique for determining substrate wetting properties. 
    more » « less
  5. 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