Suppression of Richtmyer-Meshkov Instability via Special Pairs of Shocks and Phase Transitions

The classical Richtmyer-Meshkov instability (RMI) is a hydrodynamic instability characterizing the evolution of an interface following shock loading. In contrast to other hydrodynamic instabilities such as Rayleigh-Taylor, it is known for being unconditionally unstable: regardless of the direction of shock passage, any deviations from a flat interface will be amplified. In this article, we show that for negative Atwood numbers, there exist special sequences of shocks which result in a nearly perfectly suppressed instability growth. We demonstrate this principle computationally and experimentally with stepped fliers and phase transition materials. A fascinating immediate corollary is that in specific instances, a phase-transitioning material may self-suppress RMI.

Scaling Law for the Onset of Solidification at Extreme Undercooling

Quasi-isentropic compression enables one to study the solidification of metastable liquid states that are inaccessible through other experimental means. The onset of this nonequilibrium solidification is known to depend on the compression rate and material-specific factors, but this complex interdependence has not been well characterized. In this study, we use a combination of experiments, theory, and computational simulations to derive a general scaling law that quantifies this dependence. One of its applications is a novel means to elucidate melt temperatures at high pressures.

Model for the Solid-Liquid Interfacial Free Energy at High Pressures

The free energy involved in the formation of an interface between two phases (e.g., a solid-liquid interface) is referred to as the interfacial free energy. For the case of solidification, the interfacial free energy dictates the height of the energy barrier required to nucleate stable clusters of the newly forming solid phase and is essential for producing an accurate solidification kinetics model using classical nucleation theory (CNT)-based methods. While various methods have been proposed for modeling the interfacial free energy for solid-liquid interfaces in prior literature, many of these formulations involve making restrictive assumptions or approximations, such as the system being at or near equilibrium (i.e., the system temperature is approximately equal to the melt temperature) or that the system is at pressures close to atmospheric. However, these approximations and assumptions may break down in highly non-equilibrium situations, such as in dynamic-compression experiments where metastable liquids that are undercooled by hundreds of kelvin or overpressurized by several gigapascals or more are formed before eventually solidifying. We derive a solid-liquid interfacial free-energy model for such high-pressure conditions by considering the enthalpies of interactions between pairs of atoms or molecules. We also consider the contribution of interface roughness (disordering) by incorporating a multilayer interface model known as the Temkin n-layer model. Our formulation is applicable to a diverse variety of materials, and we demonstrate it by developing models specifically for two different materials: water and gallium. We apply our interfacial free-energy formulation to CNT-based kinetics simulations of several suites of dynamic-compression experiments that cause liquid water to solidify to the high-pressure solid polymorph ice VII and have found good agreement to the observed kinetics with only minor empirical fitting.

Measuring the melting curve of iron at super-Earth core conditions

The discovery of more than 4500 extrasolar planets has created a need for modeling their interior structure and dynamics. Given the prominence of iron in planetary interiors, we require accurate and precise physical properties at extreme pressure and temperature. A first-order property of iron is its melting point, which is still debated for the conditions of Earth’s interior. We used high-energy lasers at the National Ignition Facility and in situ x-ray diffraction to determine the melting point of iron up to 1000 gigapascals, three times the pressure of Earth’s inner core. We used this melting curve to determine the length of dynamo action during core solidification to the hexagonal close-packed (hcp) structure. We find that terrestrial exoplanets with four to six times Earth’s mass have the longest dynamos, which provide important shielding against cosmic radiation.

Metastability of Liquid Water Freezing into Ice VII under Dynamic Compression

The ubiquitous nature and unusual properties of water have motivated many studies on its metastability under temperature- or pressure-induced phase transformations. Here, nanosecond compression by a high-power laser is used to create the nonequilibrium conditions where liquid water persists well into the stable region of ice VII. Through our experiments, as well as a complementary theoretical-computational analysis based on classical nucleation theory, we report that the metastability limit of liquid water under nearly isentropic compression from ambient conditions is at least 8 GPa, higher than the 7 GPa previously reported for lower loading rates.

Numerical modeling of solid-cluster evolution applied to the nanosecond solidification of water near the metastable limit

Classical nucleation theory (CNT) is a promising way to predictively model the submicrosecond kinetics of phase transitions that occur under dynamic compression, such as the suite of experiments performed over the past two decades on the solidification of liquid water to the high-pressure ice VII phase. Myint et al. [Phys. Rev. Lett. 121, 155701 (2018)] presented the first CNT-based model for these types of rapid phase transitions, but relied on an empirical scaling parameter in their transient induction model to simulate the lag time that occurs prior to the onset of significant formation of ice VII clusters in the system. To build on that study, we model the liquid water-ice VII phase transformation using a numerical discretization scheme to solve the Zel'dovich-Frenkel partial differential equation, which is a fundamental CNT-based kinetic equation that describes the statistical time-dependent behavior of solid cluster formation. The Zel'dovich-Frenkel equation inherently accounts for transience in the nucleation kinetics and eliminates the need for the empirical scaling factor used by Myint et al. One major result of this research is that transience is found to play a relatively small role in the nucleation process for the dynamic-compression time scales of the liquid water-ice VII experiments being simulated. Instead, we show that it is possible to properly model the lag time using steady-state CNT by making small refinements to the interfacial free energy value. We have also developed a new dimensionless parameter that may be applied a priori to predict whether or not transient nucleation will be important in a given dynamic-compression experiment.