Our first simulation uses the simplest model for a young SNR: a self-similar driven wave. Even in this case of spherically symmetric ejecta and a smooth circumstellar medium, there is significant structure in the SNR as a result of the instability of the contact discontinuity separating the shocked ejecta and the shocked CSM.


This image shows the gas density from a 2D simulation. The two black lines mark the reverse shock (left) and forward shock (right) as the SNR propagates to the right. As the SNR decelerates, the dense shell of shocked ejecta shown in purple is subject to the Raleigh-Taylor instability. In 2D this produces "fingers" of dense ejecta that push forward roughly half way toward the forward shock.

In this next image we show a volume rendering of the gas density from a 3D simulation, as viewed looking directly down onto the shocked layer. The low density of the shocked CSM is virtually transparent, so this image is showing the structure of the unstable layer of shocked ejecta. The Raleigh-Taylor fingers seen in 2D simulations are now seen as a network of filaments.



The impact of this instability on the observed morphology of a young SNR is illustrated here in a map of the emission measure of the ejecta integrated along the line of sight through the 3D model.

The color map corresponds to the line of sight emission measure of the ejecta only. The black circle is the location of the forward shock.

These animations show the evolution of the structure of the shocked ejecta in the expanding reference frame of the SNR:

SG301.mpg, SG301.qt




Close-up of a region containing the edge-brightened reverse shock. Links to animations:

SG301K.mpg, SG301K.qt


Here we see images generated from 3D ssdw simulations, the top row is a model running into a relic stellar wind (s=2) and the bottom row is a model driven wave running into a constant density csm (s=0). Both models have an ejecta density powerlaw with n=6.

The left images show the emission measure (density squared) of shocked ejecta integrated along a line of sight through the 3D model. The right images show the radial column density (integrating the density of shocked ejecta along radial rays) on a theta-phi plane.

For reference we show a 2D slice of the gas density for the s=2 (top) and s=0 (bottom) models.