## New Research In

### Physical Sciences

### Social Sciences

#### Featured Portals

#### Articles by Topic

### Biological Sciences

#### Featured Portals

#### Articles by Topic

- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology

# Supernova, nuclear synthesis, fluid instabilities, and interfacial mixing

Edited by William A. Goddard III, California Institute of Technology, Pasadena, CA, and approved October 3, 2018 (received for review December 31, 2017)

## Abstract

Supernovae and their remnants are a central problem in astrophysics due to their role in the stellar evolution and nuclear synthesis. A supernova’s explosion is driven by a blast wave causing the development of Rayleigh–Taylor and Richtmyer–Meshkov instabilities and leading to intensive interfacial mixing of materials of a progenitor star. Rayleigh–Taylor and Richtmyer–Meshkov mixing breaks spherical symmetry of a star and provides conditions for synthesis of heavy mass elements in addition to light mass elements synthesized in the star before its explosion. By focusing on hydrodynamic aspects of the problem, we apply group theory analysis to identify the properties of Rayleigh–Taylor and Richtmyer–Meshkov dynamics with variable acceleration, discover subdiffusive character of the blast wave-induced interfacial mixing, and reveal the mechanism of energy accumulation and transport at small scales in supernovae.

- supernovae
- nuclear synthesis
- blast waves
- Rayleigh–Taylor instabilities
- Rayleigh–Taylor interfacial mixing

Supernovae are violent, disruptive explosions of stars (1). They have been a central problem in astrophysics since their discovery and identification in the 1930s. The debris ejected from a supernova mixes with the interstellar medium, forming a supernova remnant (SNR) (2⇓–4). Young SNRs still retain information concerning the explosion process. Explosions are initial value problems; solution requires the details of what explodes (1⇓⇓⇓–5).

## Astrophysics and Fluid Dynamics

Our understanding of stellar evolution is based on the convenient assumption that stars may be treated as spherically symmetric objects, at least on average (1). This approximation, which has proven to be surprisingly successful, allows us to evolve stars up to a presumed presupernova condition (1⇓–3). We may then invent an explosion model and compare the predictions with observed SNR. Unfortunately, this is a long extrapolation, with many untested choices along the way (1). The conditions in stars imply very high Reynolds numbers, and therefore, plasma motion in stars is expected to be turbulent (or at least, disordered) (4⇓⇓–7). This may break the spherical symmetry in detail but might be approximately valid on the global scale. However, as explosion begins, spherical symmetry is firmly broken on the global scale too, so that conventional stellar evolution theory must fail (1, 5).

How this happens is of considerable interest (1⇓⇓⇓⇓⇓–7). Some supernovae produce neutron stars (pulsars), and others produce black holes. They are thought to be the major source of galactic cosmic rays. Supernovae are the dominant source of elements not produced by the Big Bang (those being hydrogen and helium). The calcium in our bones and the iron in our blood were synthesized in a supernova as were silver, gold, uranium, and thorium. The latter elements can also be produced in neutron star mergers in the “r-process,” evidenced in the mergers’ observations and the abundances in dwarf galaxies (1⇓⇓⇓⇓⇓⇓⇓⇓⇓–11). Even carbon and nitrogen are partially produced in supernovae but more so in stars that die less violently as are the “s-process” heavy elements. The details are still in debate. Supernovae are a major source of radioactive nuclei for meteorites and for the early solar system (1⇓⇓⇓⇓⇓–7).

Fig. 1 provides a detailed look at the Cassiopeia A (Cas A) SNRs, the youngest nearby SNR known in the Milky Way. The SNRs of the Cas A have been produced by the explosion of a massive star. The image in Fig. 1 is combined from 18 images taken by NASA’s Hubble Space Telescope in 2004. It shows the Cas A remnants as a broken shell of filamentary and clumpy stellar ejecta glowing with the heat generated by the passage of a shock wave from the supernova blast. The various colors of the gas indicate differences in chemical composition. Bright green filaments are rich in the oxygen, red and purple ones are rich in the sulfur, and blue ones are composed mostly of the hydrogen and the nitrogen. The oxygen and the sulfur are produced by thermonuclear burning during and just before the explosion. The sulfur is formed by the oxygen nuclear burning. Some of the sulfur lies farther from the center than some of the oxygen, which is interpreted as due to Rayleigh–Taylor instabilities (RTIs) during the explosion (12⇓–14).

Hence, the key question is: can we work backward from the SNR toward the underlying explosion to provide insight into the event that is independent of conventional stellar evolution theory?

To do so, we have to better understand fluid dynamic aspects of the multiphysics problem of supernovae, particularly Rayleigh-Taylor instability, Richtmyer–Meshkov instability (RMI), and Rayleigh–Taylor (RT)/Richtmyer–Meshkov (RM) interfacial mixing that are caused by supernova’s blast (1, 12⇓⇓⇓–16). RTI develops at the interface of the fluids with distinct densities that are accelerated against their density gradients (12⇓–14). RMI develops when the acceleration is induced by a shock refracting a perturbed interface of the fluids with distinct acoustic impedances (17, 18). Intense interfacial RT/RM mixing of the fluids ensues with time (12⇓⇓⇓⇓⇓–18).

RTI/RMI and RT/RM mixing occur at a broad range of astrophysical phenomena in low- and high-energy density regimes (1, 14). Examples include the appearance of stiff light years-long structures in molecular hydrogen clouds, the formation of accretion disks and black holes, and the processes of stellar evolution (1⇓⇓⇓⇓⇓–7). The latter ranges from a birth of a star due to the interstellar gas collapse to life of a star with the extensive material mixing in the stellar interior and to death of a star in the supernova (1⇓⇓⇓⇓⇓–7, 14, 15). In supernovae, the blast wave-driven RT/RM mixing of the outer and inner layers of the progenitor star creates conditions for synthesis of heavy and intermediate mass elements in addition to light mass elements synthesized in the star before its explosion (1⇓⇓⇓⇓⇓–7, 14, 15, 19⇓⇓⇓⇓⇓–25).

In everyday life and in extreme astrophysics environments, RT flows are observed to have similar qualitative features of their evolution (12⇓⇓⇓–16). RTI starts to develop when the fluid interface is slightly perturbed near its equilibrium state (12⇓–14). The flow transits from an initial stage, where the perturbation amplitude grows quickly, to a nonlinear stage, where the growth rate slows and the interface is transformed into a composition of small-scale shear-driven vortical structures and a large-scale coherent structure of bubbles and spikes [with the bubble (spike) being the portion of the light (heavy) fluid penetrating the heavy (light) fluid]. RT flows are usually three-dimensional (3D) and have two macroscopic-length scales—the amplitude in the acceleration direction and the spatial period or the wavelength in the normal plane. The flow eventual stage is the self-similar interfacial mixing (12⇓⇓–15).

While RT and RM flows are similar in many regards, there are also important distinctions (12⇓⇓⇓⇓⇓–18, 26⇓⇓–29). For instance, postshock RM dynamics is a superposition of two motions (17, 18, 26⇓⇓–29). These are the background motion of the fluid bulk and the growth of interface perturbations (26⇓⇓–29). In the background motion, both fluids and their interface move as whole unit in the transmitted shock direction; this motion occurs even for an ideally planar interface, and is supersonic for strong shocks. The growth of the interface perturbations is due to impulsive acceleration by the shock; it develops only for a perturbed interface (26⇓⇓–29). The rate of this growth is subsonic, and the associated motion is incompressible. The growth rate is constant initially and decays with time later. RM unstable interface is transformed to a composition of a large-scale structure of bubbles and spikes and small-scale shear-driven vortical structures. Small-scale nonuniform structures also appear in the bulk (26⇓⇓–29). Self-similar RM mixing develops, and the energy supplied initially by the shock gradually dissipates (17, 18, 26⇓⇓–29).

In RT mixing with constant acceleration, the length scale, the velocity scale, and the Reynolds number increase with time (14⇓–16). At a first glance, such flow should quickly proceed to a fully disordered turbulent state (30). Turbulence is a stochastic process insensitive to deterministic conditions with intense energy transport from the large to the small scales (30⇓–32). In supernova environments, such transport may supply to small scales some activation energy required for the synthesis of heavy mass elements (1). Recent advances in the theory and experiment of RTI and RT mixing have found, however, that, in the high- and low-energy density regimes, the properties of heterogeneous, anisotropic, nonlocal, and statistically unsteady RT mixing depart from those of homogeneous, isotropic, local, and statistically steady canonical turbulence (14⇓–16, 30⇓⇓⇓–34). High Reynolds number, while necessary, is not a sufficient condition for turbulence to occur. RT mixing exhibits more order and has stronger correlations, weaker fluctuations, and stronger sensitivity to deterministic conditions compared with canonical turbulence (14⇓–16). For RM mixing, strong sensitivity to deterministic conditions has also been found along with its nominally large Reynolds number and small-scale interfacial vortical structures (26⇓⇓–29).

In supernovae in a hydrodynamic approximation, accelerations are induced by blast waves (1). Blast waves can be viewed as strong variable shocks (1, 19⇓⇓⇓⇓⇓–25, 33). Blast wave dynamics is self-similar, and the blast wave-induced acceleration is a power law function of time (spatial coordinate) (1, 19⇓⇓⇓⇓⇓–25, 33). Several questions thus appear. What are the properties of RTI with variable acceleration, and how do they differ from those of RTI with constant acceleration and from those of shock-driven RMI? How can these properties be applied for interpretation of astrophysical data? What are their potential outcomes for stellar evolution and nuclear synthesis?

In this work, we consider hydrodynamic aspects on the multiphysics problem of supernovae and their remnants. We focus on the dynamics of RTI and RT mixing with variable acceleration in a broad parameter regime (Figs. 2–7). The acceleration is a power law function of time (spatial coordinate). We find that, for variable acceleration, RT/RM dynamics is multiscale and has two macroscopic-length scales—the amplitude in the acceleration direction and the spatial period in the normal plane (14⇓–16, 34). Depending on the exponent of the acceleration power law, the dynamics can be RT type or RM type. For RT type, the acceleration sets the timescale at early stage and defines the nonlinear dynamics and the interfacial mixing at later stages. For RM type, the initial growth rate sets the timescale at early stage; at late stages, the drag defines the nonlinear dynamics and the interfacial mixing, and the initially supplied energy gradually dissipates (33, 35). The critical values of the exponent at which the transition occurs from RT- to RM-type dynamics are distinct for the linear, nonlinear, and mixing regimes. Particularly for blast wave-induced accelerations, the linear and nonlinear dynamics are RT type, and the mixing is RM type, but RM-type mixing develops quicker than the acceleration prescribes (33, 35). While for subdiffusive mixing dynamics, superdiffusive canonical turbulence may be a challenge to occur, other mechanisms are possible for energy accumulation at small scales. They are due to small-scale nonuniform structures developing in the fluid bulk and including cumulative jets, hot and cold spots, high- and low-pressure regions, and localizations as illustrated by numerical simulations of strong shock-driven RMI (17, 18, 26⇓⇓–29, 36). Such effects should be considered in the interpretation of observational data (1). Technical details of this work are given in *SI Appendix* for corresponding sections.

## Fluid Instabilities and Interfacial Mixing

### Governing Equations.

Dynamics of ideal fluids is governed by the conservation of mass, momentum, and energy:

To rigorously describe RT dynamics, one has to solve the problem of extreme complexity: solve a system of nonlinear partial differential equations in 4D space–time, solve the boundary value problems for a subset of nonlinear partial differential equations at a nonlinear freely evolving interface and at the outside boundaries and also solve the ill-posed initial value problem, with account for singularities and secondary instabilities developing in a finite time (14, 34). A complete theory of RTI applicable at all scales and all times has yet to be developed. Rigorous theories have successfully handled the problem in well-defined approximations; empirical models have repeatedly described a broad set of data with nearly the same set of parameters. The reader is referred to review and research papers (12⇓⇓⇓⇓⇓–18, 26⇓⇓⇓–30, 33⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓⇓–54) for details of theoretical and numerical studies of RT dynamics.

It is worth noting that, despite complexity and noisiness resulting from interactions of all of the scales, RT dynamics is observed to have certain features of universality and order, and is thus eligible to first principle considerations, such as group theory (14, 15, 34, 55, 56). For linear and nonlinear RTI, group theory analysis uses theory of discrete groups to solve the boundary value and initial value problems (34, 48⇓–50). For RT mixing, group theory is implemented in the momentum model with equations that have the same symmetries and scaling transformations as the governing equations (14, 16, 34, 40, 54). Some principal results of group theory analysis—the multiscale character of nonlinear dynamics, to which both the spatial period λ and the amplitude h contribute; the tendency to keep isotropy in the plane normal to the acceleration and the discontinuous dimensional cross-over; the order in RT mixing that may be coexistent with a quasiturbulent state—self-consistently explain the observations (14, 15, 26⇓⇓–29, 34, 48⇓–50, 54⇓⇓⇓–58).

### Group Theory Approach.

Group theory for linear and nonlinear dynamics studies RT flows that consist of large-scale coherent structures in which dynamics is potential,

Group theory for the mixing dynamics finds that, in RT mixing, the momentum and energy are gained and lost at any scale; the dynamics of a parcel of fluid is governed by a balance per unit mass the rates of momentum gain,

### Smoothed Particle Hydrodynamics Simulations.

Supernova environments are characterized by the conditions of high energy density, strong shocks, sharp changes of flow fields, large perturbations, and small effects of dissipation and diffusion. Numerical modeling of these extreme regimes is a challenge. This is because numerical methods should satisfy numerous competing requirements to capture shocks, track interfaces, and accurately account for dissipative processes. An efficient approach for modeling strong shock-driven dynamics is the smoothed particle hydrodynamics (SPH) implemented in the smoothed particle hydrodynamics code (SPHC) (36).

The SPH is a Lagrangian method representing a continuous fluid by means of fixed mass SPH particles and thus, reducing the governing partial differential equations to ordinary differential equations (26⇓⇓–29, 36). SPHC has been originated in astrophysics, and in addition to astrophysical gravitational problems, it has been successfully applied, tested, and validated in multiphysics problems in fluids, plasmas, materials for modeling the strong shock-driven RMI, the Noh problem, and other flow problems as well as reactive and supercritical fluids, material transformation under impact, ablation process in hypersonic flows, and charge imbalance in plasmas (36). When applied to strong shock-driven RMI, it has achieved an excellent agreement with the experiments and with the rigorous zero-order, linear, and group theories. The latter includes the multiscale character of the nonlinear dynamics and is evinced, for instance, in the flattening of RM bubbles, flow fields structure, and sensitivity to deterministic conditions (26⇓⇓–29, 36).

Here, we apply group theory to study RTI/RMI and RT/RM mixing with variable acceleration in a broad parameter regime. We further conduct SPHC simulations of strong shocks-driven RM flows to study their sensitivity to deterministic conditions and their small-scale dynamics at the interface and in the bulk. We find the properties of blast wave-induced RT/RM dynamics and propose the mechanisms for energy accumulation and transport at small scales in supernovae.

## RTI with Variable Acceleration

We consider RT dynamics with time-varying acceleration *Gk* parameterizes the acceleration strength. Time is

At early stage, for

At late stage, for

For

The multiplicity of the nonlinear solutions is associated with the nonlocal and singular character of the interfacial dynamics (34, 45, 46, 48⇓–50). The solutions exist and converge with increase in approximation order. The number of the family parameters is identified by symmetry of the global flow (34). For group p6mm, the dynamics is highly isotropic,

The multiplicity is also due to the presence of shear at the interface (58). Defining the shear as the spatial derivative of the jump of tangential velocity at the interface,

The fastest stable solution is the physically significant solution (Fig. 2) (14, 34, 48⇓–50, 58). For the fastest stable solutions, the dependence of curvature and velocity

Regular asymptotic solutions have important physics properties. For these solutions, there is effectively no motion of the fluids in the bulk and away from the interface, there is intense motion near the interface, and shear is present at the interface leading to formation of interfacial vortical structures, in agreement with observations (Fig. 3) (34, 48⇓–50).

Regular asymptotic solutions have important global properties. The flow tends to conserve isotropy in the plane (34, 48⇓–50). The 3D highly symmetric dynamics is universal (34, 48). That is, on the substitution

At

RT and RM families have similar physical, mathematical, and global properties (Figs. 2 and 3). These include the structure of the velocity fields, the existence of the family of solutions, their dependence on the flow symmetry and the interfacial shear, the tendency to keep isotropy in the plane normal to the acceleration, and the discontinuity of the 3D–2D cross-over. However, their local properties are distinct (Fig. 2) (34, 50). Particularly, for RM bubbles, the velocity is a monotone function on the curvature, and the fastest stable solution corresponds to a flat bubble with

## Mixing Dynamics

### RT/RM Dynamics for Time-Varying Acceleration.

In the nonlinear regime with

### RT-Type Mixing.

Properties of the asymptotic solutions indicate that, for

For

### RM-Type Mixing.

RM-type mixing develops for

### Space-Varying Acceleration.

Similar analyses can be conducted when the acceleration is a power law function on the spatial coordinate,

The mixing dynamics is RT type for

Fig. 4 illustrates RT/RM-type mixing for various values of the exponents a and n. Solutions are derived for C being a stochastic process with log-normal distribution with the mean

### Blast Wave-Induced Mixing.

Consider now RT/RM mixing induced by blast waves with the first kind (Sedov–Taylor) self-similarity and the second kind (Guderley–Stanyukovich–Landau) self-similarity (19⇓⇓⇓⇓⇓–25). Note that exponents

For Sedov–Taylor self-similar dynamics, the scaling dependence of the solution is set by the energy release Ε,

For Guderley–Stanyukovich–Landau self-similar dynamics, the load history and momentum transport should be accounted for. The corresponding invariant

### Small-Scale Nonuniform Structures of the Flow Fields.

Hence, depending on the drag value, the blast-wave induced mixing can be RT type with *a* ∼ *a*_{cr} or RM type with *a* < *a*_{cx}, *a*_{cr} ∈ (−2,−1); in either case, larger velocities (velocity fluctuations) correspond to small length scales, and the dynamics is sensitive to deterministic conditions. In canonical turbulence, large velocities (velocity fluctuations) correspond to large length scales, and the dynamics is independent of deterministic conditions (31, 32). In the blast wave-induced (subdiffusive) mixing, the canonical (superdiffusive) turbulence may be a challenge to develop (unless there is a source, other than gravity, supplying turbulent energy to the fluid system). If so, what are other possible mechanisms for energy accumulation at small scales, which is necessary for nuclear synthesis in supernovae (1)?

Such mechanisms may exist due to small-scale nonuniform structures of the flow fields (26⇓⇓–29, 59). Particularly, subdiffusive processes are often characterized by localizations. If a spot with high energy fluctuation appears in the flow, it can be trapped due to the slow dynamics. If the spot is hot enough, it may initiate a nuclear reaction accompanied by some energy release and may further induce a chain of other reactions and processes (1⇓–3). Blast waves are the special strong shocks (24). In strong shock-driven RMI, the shock–interface interaction may lead to intense production of small-scale nonuniform structures in the bulk in addition to small-scale shear-driven vortical structure at the interface (26⇓⇓–29). These nonuniform structures may include cumulative jets, hot and cold spots, high- and low-pressure regions, and may enable strong energy fluctuations at small scales.

## Strong Shock-Driven RM Flows

### SPH Simulations.

To illustrate scale coupling in strong shock-driven RMI, we use the SPHC in a hydrodynamic approximation for ideal monoatomic gases with adiabatic indexes *M* = {3,5,7,10} defined relatively to the light fluid with the speed of sound *c*_{l} (26⇓⇓–29). The interface is normal to the shock; the initial perturbation wavelength and amplitude are λ and

Fig. 5*A* illustrates SPHC simulations of the postshock dynamics of strong shock-driven RMI. The superposition of the growth of the interface perturbation with the background motion of the fluids, the formation of large-scale coherent structure of bubbles and spikes, the bubble flattening at late times, and the occurrence at small scales of interfacial vortical structures and bulk-immersed cumulative jets are seen in Fig. 5*A*. The flow regions are shown, with red (blue) for the light (heavy) fluid particles and green for the light fluid interfacial particles in Fig. 5*A*.

### Sensitivity of RM Dynamics to Deterministic Conditions.

In RMI with a single-wave initial perturbation, the initial growth rate *B* shows the dependence of the initial growth rate on initial perturbation amplitude for given *B* and *C*), suggesting that, in addition to wavelength λ, RM dynamics has the characteristic amplitude scale

For a multimode initial perturbation, the order and disorder in RM flow are highly sensitive to deterministic conditions, including the wavelengths, the amplitudes, and the relative phases of the initial perturbation waves (29). Fig. 6 show snapshots of the flow regions of late stages of the RMI, with the same corresponding values of the Mach and Atwood numbers and the observational time. The initial perturbations have two waves with the same amplitudes and wavelengths *Left*, the waves are antiphase, *Right*, the waves have random phase,

### Small-Scale Structures in RM Flow.

SPHC simulations accurately capture small-scale dynamics in strong shock-driven RMI. At the interface, we observe small-scale shear-driven vortical structures due to the Kelvin–Helmholtz instability (26⇓⇓–29). In the bulk, between the transmitted shock and the interface, we observe small-scale nonuniform localized structures of the flow fields (Fig. 7) (26). These include the “reverse” cumulative jets, which are short and energetic and develop due to collisions of converging fluid flows in the bulk (56, 60); the hot (cold) spots, which are the localized regions with the temperature much higher (lower) than that of the ambient; and the regions with high (low) pressure (26⇓⇓–29). These well-pronounced small-scale structures are volumetric in nature (Fig. 7). They develop, since the flow dynamics is adjusted to the motion of the interface.

Hence, the strong shock-driven RM dynamics is sensitive to deterministic conditions, may keep order at large and small scales, and may have localized nonuniform structures at small scales in the bulk in addition to shear-driven vortical structures at the interface.

## Discussion

Supernovae and their remnants are a central problem in astrophysics due to their role in the formation of neutron stars and black holes in the processes of stellar evolution and nuclear synthesis (Fig. 1) (1). We have considered this multiphysics problem in a hydrodynamic approximation, with supernova blast causing the development of RTI, RMI, and RT/RM interfacial mixing with variable acceleration (1⇓⇓⇓⇓⇓–7). We have applied group theory (14, 34) to study RTI/RMI and RT/RM mixing with variable acceleration, and identified properties of linear, nonlinear, and mixing RT/RM dynamics that have not been discussed before (Figs. 2–7).

We have found that, for

The mixing is RT type for

For RT-type mixing with

In blast wave-driven dynamics, the acceleration exponent can be

What are potential outcomes of this hydrodynamics for supernova and nuclear synthesis (1)?

In blast wave-driven RT/RM mixing in supernovae, canonical turbulence is traditionally considered as the mechanism for energy accumulation at small scales. According to our results, for the acceleration parameters typical for blast waves (19⇓⇓⇓⇓⇓–25), superdiffusive turbulence (31, 32) may be a challenge to implement in subdiffusive RT/RM mixing. However, the conditions of heterogeneity, nonlocality, anisotropy, and statistical unsteadiness that are common for RT/RM flows (14, 15, 33, 34) may lead to appearance of small-scale nonuniform structures in the bulk of the blast wave-driven mixing flows, whereas the slow subdiffusive transport may result in energy localization and trapping at small scales (26⇓⇓–29, 59). These effects are consistent with and may explain the richness of structures observed in supernovae, including the Cas A (Fig. 1) (1). We suggest that such effects be considered in interpretation of observational data.

Our work focuses on hydrodynamic aspects of the multiphysics problem of supernovae and their remnants and serves for systematic studies of the problem in perspective.

## Acknowledgments

We thank the University of Western Australia, the National Science Foundation, and the Theoretical Astrophysics Program and the Steward Observatory at the University of Arizona for support.

## Footnotes

- ↵
^{1}To whom correspondence may be addressed. Email: snezhana.abarzhi{at}gmail.com or wdarnett{at}gmail.com.

Author contributions: S.I.A. and W.D.A. designed research; S.I.A., A.K.B., A.N., A.P., N.C.S., R.F.S., and W.D.A. performed research; R.F.S. contributed new reagents/analytic tools; S.I.A., A.K.B., A.N., A.P., R.F.S., and W.D.A. analyzed data; S.I.A. guided research on fluid instabilities and mixing, including scientific and organizational aspects; A.K.B., A.N., A.P., N.C.S. studied fluid instabilities and mixing analytically and numerically; R.F.S. developed and supported the smoothed particle hydrodynamics code and provided advice on the code use; W.D.A. guided the project connecting the observational data, supernova, and nuclear synthesis to fluid instabilities and mixing; and S.I.A., A.K.B., A.N., A.P., R.F.S., and W.D.A. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1714502115/-/DCSupplemental.

Published under the PNAS license.

## References

- ↵
- Arnett D

- ↵
- Arnett WD

- ↵
- ↵
- ↵
- Bazan G,
- Arnett WD

- ↵
- Arnett WD,
- Meakin C

- ↵
- Arnett WD,
- Meakin C,
- Young PA

- ↵
- ↵
- Rosswog S, et al.

- ↵
- Rosswog S, et al.

- ↵
- Hansen TT,
- Simon JD,
- Marshall JL, et al.

- ↵
- Lord R

- ↵
- Davies RM,
- Taylor GI

- ↵
- ↵
- ↵
- ↵
- ↵
- Meshkov EE

- ↵
- Sedov LI

- ↵
- Taylor GI

- ↵
- Guderley G

- ↵
- Stanyukovich KP

- ↵
- Stanyukovich KP

- ↵
- Sedov LI

- ↵
- Landau LD,
- Lifshitz EM

- ↵
- ↵
- Dell Z,
- Stellingwerf RF,
- Abarzhi SI

- ↵
- Dell Z, et al.

- ↵
- Pandian A,
- Stellingwerf RF,
- Abarzhi SI

- ↵
- Neuvazhaev VE

- ↵
- Kolmogorov AN

- ↵
- Kolmogorov AN

- ↵
- Swisher N, et al.

- ↵
- Abarzhi SI

- ↵
- Pandian A,
- Swisher NC,
- Abarzhi SI

- ↵
- Stellingwerf RF

*Advances in the Free-Lagrange Method*, Lecture Notes in Physics, eds Trease HE, Fritts MF, Crowley WP (Springer, Berlin), pp 239–247. - ↵
- Chandrasekhar S

- ↵
- Kull HJ

- ↵
- ↵
- Abarzhi SI,
- Rosner R

- ↵
- Haan SW

- ↵
- ↵
- ↵
- ↵
- Garabedian PR

- ↵
- Inogamov NA

- ↵
- ↵
- ↵
- Abarzhi SI,
- Nishihara K,
- Glimm J

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Meshkov EE

*Studies of Hydrodynamic Instabilities in Laboratory Experiments*(FGUC-VNIIEF, Sarov, Russia). Russian. - ↵
- ↵
- Bhowmick AK,
- Abarzhi SI

- ↵
- Metzler R,
- Klafter J

- ↵
- Sokolov IV

## Citation Manager Formats

## Sign up for Article Alerts

## Article Classifications

- Physical Sciences
- Astronomy