skip to main content

Title: Hydrodynamic instability at impact interfaces and planetary implications

Impact-induced mixing between bolide and target is fundamental to the geochemical evolution of a growing planet, yet aside from local mixing due to jetting – associated with large angles of incidence between impacting surfaces – mixing during planetary impacts is poorly understood. Here we describe a dynamic instability of the surface between impacting materials, showing that a region of mixing grows between two media having even minimal initial topography. This additional cause of impact-induced mixing is related to Richtmyer-Meshkov instability (RMI), and results from pressure perturbations amplified by shock-wave refraction through the corrugated interface between impactor and target. However, unlike RMI, this new impact-induced instability appears even if the bodies are made of the same material. Hydrocode simulations illustrate the growth of this mixing zone for planetary impacts, and predict results suitable for experimental validation in the laboratory. This form of impact mixing may be relevant to the formation of stony-iron and other meteorites.

; ;
Award ID(s):
Publication Date:
Journal Name:
Nature Communications
Nature Publishing Group
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 geophysicalmore »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.« less
  2. Abstract

    Recent developments of burst-mode lasers and imaging systems have opened new realms of simultaneous diagnostics for velocity and density fields at a rate of 1 kHz–1 MHz. These enable the exploration of previously unimaginable shock-driven turbulent flow fields that are of significant importance to problems in high-energy density physics. The current work presents novel measurements using simultaneous measurements of velocity and scalar fields at 60 kHz to investigate Richtmyer-Meshkov instability (RMI) in a spatio-temporal approach. The evolution of scalar fields and the vorticity dynamics responsible for the same are shown, including the interaction of shock with the interface. This temporal information is used to validate two vorticity-deposition models commonly used for initiation of large scale simulations, and have been previously validated only via simulations or integral measures of circulation. Additionally, these measurements also enable tracking the evolution and mode merging of individual flow structures that were previously not possible owing to inherently random variations in the interface at the smallest scales. A temporal evolution of symmetric vortex merging and the induced mixing prevalent in these problems is presented, with implications for the vortex paradigms in accelerated inhomogenous flows.

  3. Abstract

    A combined analysis of smooth and non-smooth bifurcations captures the interplay of different qualitative transitions in a canonical model of an impact pair, a forced capsule in which a ball moves freely between impacts on either end of the capsule. The analysis, generic for the impact pair context, is also relevant for applications. It is applied to a model of an inclined vibro-impact energy harvester device, where the energy is generated via impacts of the ball with a dielectric polymer on the capsule ends. While sequences of bifurcations have been studied extensively in single- degree-of-freedom impacting models, there are limited results for two-degree-of-freedom impacting systems such as the impact pair. Using an analytical characterization of impacting solutions and their stability based on the maps between impacts, we obtain sequences of period doubling and fold bifurcations together with grazing bifurcations, a particular focus here. Grazing occurs when a sequence of impacts on either end of the capsule are augmented by a zero-velocity impact, a transition that is fundamentally different from the smooth bifurcations that are instead characterized by eigenvalues of the local behavior. The combined analyses allow identification of bifurcations also on unstable or unphysical solutions branches, which we termmore »ghost bifurcations. While these ghost bifurcations are not observed experimentally or via simple numerical integration of the model, nevertheless they can influence the birth or death of complex behaviors and additional grazing transitions, as confirmed by comparisons with the numerical results. The competition between the different bifurcations and their ghosts influences the parameter ranges for favorable energy output; thus, the analyses of bifurcation sequences yield important design information.

    « less

    Planet engulfment can be inferred from enhancement of refractory elements in the photosphere of the engulfing star following accretion of rocky planetary material. Such refractory enrichments are subject to stellar interior mixing processes, namely thermohaline mixing induced by an inverse mean-molecular-weight gradient between the convective envelope and radiative core. Using mesa stellar models, we quantified the strength and duration of engulfment signatures following planet engulfment. We found that thermohaline mixing dominates during the first ∼5–45 Myr post-engulfment, weakening signatures by a factor of ∼2 before giving way to depletion via gravitational settling on longer time-scales. Solar metallicity stars in the 0.5–1.2 M⊙ mass range have observable signature time-scales of ∼1 Myr–8 Gyr, depending on the engulfing star mass and amount of material engulfed. Early type stars exhibit larger initial refractory enhancements but more rapid depletion. Solar-like stars (M = 0.9–1.1 M⊙) maintain observable signatures (>0.05 dex) over time-scales of ∼20 Myr–1.7 Gyr for nominal 10 M⊕ engulfment events, with longer-lived signatures occurring for low-metallicity and/or hotter stars (1 M⊙, ∼2–3 Gyr). Engulfment events occurring well after the zero-age main sequence produce larger signals due to suppression of thermohaline mixing by gravitational settling of helium (1 M⊙, ∼1.5 Gyr). These results indicate that it may be difficult to observe engulfment signatures in solar-like starsmore »that are several Gyr old.

    « less
  5. Abstract. The realization of the difficulty of limiting global-meantemperatures to within 1.5 or 2.0 ∘C abovepre-industrial levels stipulated by the 21st Conference of Parties inParis has led to increased interest in solar radiation management (SRM)techniques. Proposed SRM schemes aim to increase planetary albedo to reflectmore sunlight back to space and induce a cooling that acts to partiallyoffset global warming. Under the auspices of the Geoengineering ModelIntercomparison Project, we have performed model experiments whereby globaltemperature under the high-forcing SSP5-8.5 scenario is reduced to followthat of the medium-forcing SSP2-4.5 scenario. Two different mechanisms toachieve this are employed: the first via a reduction in the solar constant(experiment G6solar) and the second via modelling injections of sulfurdioxide (experiment G6sulfur) which forms sulfate aerosol in thestratosphere. Results from two state-of-the-art coupled Earth system models(UKESM1 and CESM2-WACCM6) both show an impact on the North AtlanticOscillation (NAO) in G6sulfur but not in G6solar. Both models show apersistent positive anomaly in the NAO during the Northern Hemisphere winterseason in G6sulfur, suggesting an increase in zonal flow and an increase inNorth Atlantic storm track activity impacting the Eurasian continent and leadingto high-latitude warming over Europe and Asia. These results are broadlyconsistent with previous findings which show similar impacts fromstratospheric volcanicmore »aerosol on the NAO and emphasize that detailedmodelling of geoengineering processes is required if accurate impacts of SRMeffects are to be simulated. Differences remain between the two models inpredicting regional changes over the continental USA and Africa, suggestingthat more models need to perform such simulations before attempting to drawany conclusions regarding potential continental-scale climate change underSRM.« less