Drop impact experiments allow the modelling of a wide variety of natural processes, from raindrop impacts to planetary impact craters. In particular, interpreting the consequences of planetary impacts requires an accurate description of the flow associated with the cratering process. In our experiments, we release a liquid drop above a deep liquid pool to investigate simultaneously the dynamics of the cavity and the velocity field produced around the air–liquid interface. Using particle image velocimetry, we analyse quantitatively the velocity field using a shifted Legendre polynomial decomposition. We show that the velocity field is more complex than considered in previous models, in relation to the non-hemispherical shape of the crater. In particular, the velocity field is dominated by degrees 0 and 1, with contributions from degree 2, and is independent of the Froude and the Weber numbers when these numbers are large enough. We then derive a semi-analytical model based on the Legendre polynomial expansion of an unsteady Bernoulli equation coupled with a kinematic boundary condition at the crater boundary. This model explains the experimental observations and can predict the time evolution of both the velocity field and the shape of the crater, including the initiation of the central jet.
more »
« less
Rayleigh–Taylor instability in impact cratering experiments
When a liquid drop strikes a deep pool of a target liquid, an impact crater opens while the liquid of the drop decelerates and spreads on the surface of the crater. When the density of the drop is larger than the target liquid, we observe mushroom-shaped instabilities growing at the interface between the two liquids. We interpret this instability as a spherical Rayleigh–Taylor instability due to the deceleration of the interface, which exceeds the ambient gravity. We investigate experimentally the effect of the density contrast and the impact Froude number, which measures the importance of the impactor kinetic energy to gravitational energy, on the instability and the resulting mixing layer. Using backlighting and planar laser-induced fluorescence methods, we obtain the position of the air–liquid interface, an estimate of the instability wavelength, and the thickness of the mixing layer. We derive a model for the evolution of the crater radius from an energy conservation. We then show that the observed dynamics of the mixing layer results from a competition between the geometrical expansion of the crater, which tends to thin the layer, and entrainment related to the instability, which increases the layer thickness. The mixing caused by this instability has geophysical implications for the impacts that formed terrestrial planets. Extrapolating our scalings to planets, we estimate the mass of silicates that equilibrates with the metallic core of the impacting bodies.
more »
« less
- Award ID(s):
- 2020249
- NSF-PAR ID:
- 10332387
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 937
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract We integrate high‐resolution full‐waveform velocity models with seismic reflection images to map the peak ring and impactite stratigraphy at the Chicxulub structure. International Ocean Discovery Program/International Continental scientific Drilling Program Site M0077 provides ground truth for our interpretations. The peak ring is narrower (∼10 km width) where it is high relief (600–700 m below seafloor) and wider (∼15 km width) where it is lower relief (1,000–1,200 m below seafloor). Both target asymmetry and angle of impact could have contributed to observed differences in peak ring morphology. We interpret a layer of lowered velocities as a resurge layer formed from the ocean resurge, seiche, and returning tsunami flowing into the newly formed impact basin. This graded suevite layer has an average thickness of 187 ± 58 m with only local thickness differences within the annular trough, peak ring, and central basin. These observations suggest that the returning ocean was of substantial height and energetic enough to carry debris across the entire topographic peak ring. We map impact melt rock throughout the crater, with a thick impact melt sheet in the central basin (>500 m), thin intermittent melt rock capping the peak ring, and a ∼500‐m thick layer of melt rock in the annular trough near the peak ring that thins toward the crater rim. We estimate that ∼70%–75% of the melt rock volume is in the central basin. We image features above and adjacent to the central basin melt sheet that we interpret as upflow zones associated with a long‐lasting hydrothermal system.
-
David Weitz (Ed.)more » « less
Water striders are abundant in areas with high humidity and rainfall. Raindrops can weigh more than 40 times the adult water strider and some pelagic species spend their entire lives at sea, never contacting ground. Until now, researchers have not systematically investigated the survival of water striders when impacted by raindrops. In this experimental study, we use high-speed videography to film drop impacts on water striders. Drops force the insects subsurface upon direct contact. As the ensuing crater rebounds upward, the water strider is propelled airborne by a Worthington jet, herein called the first jet. We show the water strider’s locomotive responses, low density, resistance to wetting when briefly submerged, and ability to regain a super-surface rest state, rendering it impervious to the initial impact. When pulled subsurface during a second crater formation caused by the collapsing first jet, water striders face the possibility of ejection above the surface or submersion below the surface, a fate determined by their position in the second crater. We identify a critical crater collapse acceleration threshold ∼ 5.7 gravities for the collapsing second crater which determines the ejection and submersion of passive water striders. Entrapment by submersion makes the water strider poised to penetrate the air–water interface from below, which appears impossible without the aid of a plastron and proper locomotive techniques. Our study is likely the first to consider second crater dynamics and our results translate to the submersion dynamics of other passively floating particles such as millimetric microplastics atop the world’s oceans.
-
Understanding the development and breakup of interfacial waves in a two-phase mixing layer between the gas and liquid streams is paramount to atomization. Due to the velocity difference between the two streams, the shear on the interface triggers a longitudinal instability, which develops to interfacial waves that propagate downstream. As the interfacial waves grow spatially, transverse modulations arise, turning the interfacial waves from quasi-two-dimensional to fully three-dimensional. The inlet gas turbulence intensity has a strong impact on the interfacial instability. Therefore, parametric direct numerical simulations are performed in the present study to systematically investigate the effect of the inlet gas turbulence on the formation, development and breakup of the interfacial waves. The open-source multiphase flow solver, PARIS, is used for the simulations and the mass–momentum consistent volume-of-fluid method is used to capture the sharp gas–liquid interfaces. Two computational domain widths are considered and the wide domain will allow a detailed study of the transverse development of the interfacial waves. The dominant frequency and spatial growth rate of the longitudinal instability are found to increase with the inlet gas turbulence intensity. The dominant transverse wavenumber, determined by the Rayleigh–Taylor instability, scales with the longitudinal frequency, so it also increases with the inlet gas turbulence intensity. The holes formed in the liquid sheet are important to the disintegration of the interfacial waves. The hole formation is influenced by the inlet gas turbulence. As a result, the sheet breakup dynamics and the statistics of the droplets formed also change accordingly.more » « less
-
Abstract Diffusive convection can occur when two constituents of a stratified fluid have opposing effects on its stratification and different molecular diffusivities. This form of convection arises for the particular temperature and salinity stratification in the Arctic Ocean and is relevant to heat fluxes. Previous studies have suggested that planetary rotation may influence diffusive-convective heat fluxes, although the precise physical mechanisms and regime of rotational influence are not well understood. A linear stability analysis of a temperature and salinity interface bounded by two mixed layers is performed here to understand the stability properties of a diffusive-convective system, and in particular the transition from non-rotating to rotationally-controlled heat transfer. Rotation is shown to stabilize diffusive convection by increasing the critical Rayleigh number to initiate instability. In the rotationally-controlled regime, a −4/3 power law is found between the critical Rayleigh number and the Ekman number, similar to the scaling for rotating thermal convection. The transition from non-rotating to rotationally-controlled convection, and associated drop in heat fluxes, is predicted to occur when the thermal interfacial thickness exceeds about 4 times the Ekman layer thickness. A vorticity budget analysis indicates how baroclinic vorticity production is counteracted by the tilting of planetary vorticity by vertical shear, which accounts for the stabilization effect of rotation. Finally, direct numerical simulations yield generally good agreement with the linear stability analysis. This study, therefore, provides a theoretical framework for classifying regimes of rotationally-controlled diffusive-convective heat fluxes, such as may arise in some regions of the Arctic Ocean.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2022.111
