We present two-dimensional MHD simulations of the evolution of a young Type Ia supernova remnant during its interaction with an interstellar cloud of comparable size at impact. We include for the first time in such simulations explicit relativistic electron transport, including spectral information using a simple but effective scheme that follows their acceleration at shocks and subsequent transport. From this information we also model radio synchrotron emission, including spectra. The principal conclusions from these experiments are: 1) Independent of the cloud interaction, the SNR reverse shock can be an efficient site for particle acceleration in a young SNR. 2) At these early times the synchrotron spectral index due to electrons accelerated at the primary shocks should be close to 0.5 unless those shocks are modified by cosmic-ray pressures. However, interaction with the cloud generates regions of distinctly steeper spectra, which may complicate interpretation in terms of global dynamical models for SNR evolution. 3) The internal motions within the SNR become highly turbulent following the cloud interaction. 4) An initially uniform interstellar magnetic field is preferentially amplified along the magnetic equator of the SNR, primarily due to biased amplification by instabilities. Independent of the external field configuration, there is a net radial direction to this field inside the SNR. 5) Filamentary radio structures correlate well with magnetic filaments, while diffuse emission follows the electron distribution. 6) Interaction with the cloud enhances both the electron population and the radio emission.

Recent HST observations of the Crab Nebula show filamentary structures that appear to originate from the Rayleigh-Taylor (R-T) instability operating on the supernova ejecta accelerated by the pulsar-driven wind. In order to understand the origin and formation of the filaments in the Crab Nebula, we study the interaction of a pulsar wind with the uniformly expanding supernova remnant by means of numerical simulation. By performing two-dimensional numerical simulations, we find three independent instabilities in the interaction region between the pulsar wind and the expanding supernova remnant. The most important instability develops as the shock driven by the pulsar bubble becomes accelerated ($r \propto t^{6/5}$). The instability produces pronounced filamentary structures that resemble the observed filaments in the Crab Nebula. Our numerical simulations can reproduce important observational features of the Crab Nebula. The high density heads in the R-T finger tips are produced because of the compressibility of the gas. The density of these heads is found to be about 10 times higher than other regions in the fingers. The mass contained in the R-T fingers is found to be $60 % - 75 %$ of the total shocked mass and the kinetic energy within the R-T fingers is $55 % - 72 %$ of the total kinetic energy of the shocked flow. The R-T fingers are found to accelerate with a slower rate than the shock front, which is consistent with the observations. By comparing our simulations and the observations, we infer that the some finger-like filaments (region F or G in Hester's observation) started to develop about 657 years ago.

This paper reports on a workshop hosted by the University of Minnesota, March 23-26, 1997. It addressed fundamental dynamical issues associated with the evolution of shell supernova remnants and the relationships between supernova remnants and their environments. The workshop considered, in addition to classical shell SNRs, dynamical issues involving X-ray filled composite remnants and pulsar driven shells, such as that in the Crab Nebula. Approximately 75 participants with wide ranging interests attended the workshop. An even larger community helped through extensive on-line debates prior to the meeting. Each of the several sessions, organized mostly around chronological labels, also addressed some underlying, general physical themes: How are SNR dynamics and structures modified by the character of the CSM and the ISM and vice versa? How are magnetic fields generated in SNRs and how do magnetic fields influence SNRs? Where and how are cosmic-rays (electrons and ions) produced in SNRs and how does their presence influence or reveal SNR dynamics? How does SNR blast energy partition into various components over time and what controls conversion between components? In lieu of a proceedings volume, we present here a synopsis of the workshop in the form of brief summaries of the workshop sessions. The sharpest impressions from the workshop were the crucial and under-appreciated roles that environments have on SNR appearance and dynamics and the critical need for broad-based studies to understand these beautiful, but enigmatic objects. \\

We examine the formation and evolution of the density enhancement (density spike) that appears downstream of strong, cosmic-ray-modified shocks. This feature results from temporary overcompression of the flow by the combined cosmic-ray shock precursor/gas subshock. Formation of the density spike is expected whenever shock modification by cosmic-ray pressure increases strongly. That occurence may be anticipated for newly generated strong shocks or for cosmic-ray-modified shocks encountering a region of higher external density, for example. The predicted mass density within the spike increases with the shock Mach number and with shocks more dominated by cosmic-ray pressure. We find this spike to be linearly unstable under a modified Rayleigh-Taylor instability criterion at the early stage of its formation. We confirm this instability numerically using two independent codes based on the two-fluid model for cosmic-ray transport. These two-dimensional simulations show that the instability grows impulsively at early stages and then slows down as the gradients of total pressure and gas density decrease. Observational discovery of this unstable density spike behind shocks, possibly through radio emission enhanced by the amplified magnetic fields would provide evidence for the existence of strongly cosmic-ray modified shock structures.

We study the radio emission from young supernova remnants by means of 3D numerical MHD simulations of the Rayleigh-Taylor instability in the shell of the remnant. The computation is carried out in spherical polar coordinates ($r, \theta, \phi$) by using a moving grid technique which allows us to finely resolve the shell. Three-dimensional result shows more turbulent (complex) structures in the mixing region than the two-dimensional result, and the instability is found to deform the reverse shock front. Stokes parameters (I,Q, and U) are computed to study the radio properties of the remnant. The total intensity map shows two distinctive regions (inner and outer shells). The inner shell appears to be complex and turbulent exhibiting loop structures and plumes as a result of the Rayleigh-Taylor instability, while the outer shell is faint and laminar due to the shocked uniform ambient magnetic fields. The inner shell resembles the observed radio structure in the main shell of young SNRs, which is evidence that the Rayleigh-Taylor instability is an ongoing process in young SNRs. When only the peculiar components of the magnetic fields generated by the instability are considered, the polarization B-vector in the inner radio shell is preferentially radial with about $20 \sim 50\%$ of fractional polarization which is higher than the observed value. The fractional polarization is lowest in the turbulent inner shell and increases outward, which is attributed to the geometric effect. The polarized intensity is found to be correlated with the total intensity. We demonstrate that the polarized intensity from the turbulent region can dominate over the polarized intensity from the shocked uniform fields if the amplified field is sufficiently strong. Therefore, we conclude that the Rayleigh-Taylor instability can explain the dominant radial magnetic

We discover a new dynamical mechanism that significantly enhances the growth of Rayleigh-Taylor fingers developed near the contact interface between the supernova ejecta and swept-up ambient gas in young supernova remnants if the supernova remnant expands into a clumpy (cloudy) circumstellar medium. Our numerical simulation demonstrates that large Rayleigh-Taylor fingers can obtain a sufficient terminal velocity to protrude through the forward shock front by taking extra kinetic energy from vorticies generated by shock-cloud interactions. We suggest this mechanism as a means to generate the aspherical expansion of the supernova ejecta. Ambient magnetic fields are stretched and amplified as the Rayleigh-Taylor fingers protrude, possibly leading to strongly enhanced radio emission. The material in the protrusions originates from the ejected stellar material with greatly enhanced heavy elements. Therefore, it can be a strong X-ray emitter. The timescale for the Rayleigh-Taylor fingers to reach the forward shock depends on the size, mass density and distribution of clouds being engulfed by the supernova shock, although the details will require further numerical investigation.

Young supernova remnants such as Tycho generally exhibit a bright circular clumpy shell in both radio and X-ray emission. For several young remnants, various arguments suggest that the magnetic field is larger than can be explained by compression of a few $\mu G$ ambient magnetic field by the shock wave. Radio polarization studies reveal a net radial orientation of magnetic fields in the shell which cannot be explained by the simple compression either. We model Rayleigh-Taylor instability at the interface of the ejecta and the shocked ambient medium to explain these observations. We have performed multidimensional MHD simulations of the instability in the shell of a Type-I supernova remnant for the first time utilizing a moving grid technique which allows us to follow the growth of the instability and its effect on the local magnetic field in detail. We find that the evolution of the instability is very sensitive to the deceleration of the ejecta and the evolutionary stage of the remnant. As the reverse shock enters the inner uniform density region, Chevalier's self-similar stage ends and the thickness of radio shell increases and the instability weakens. Our simulation shows that Rayleigh-Taylor and Kelvin-Helmholtz instabilities amplify ambient magnetic fields locally by as much as a factor of 60 around dense fingers due to stretching, winding, and compression. Globally, the amount of magnetic-field amplification is low and the magnetic energy density reaches only about 0.3% of the turbulent energy density at the end of simulation. Strong magnetic field lines draped around the fingers produce the radial B-vector polarization, whereas thermal bremsstrahlung from the dense fingers themselves produce the clumpy X-ray emission. As a result, the X-ray emission peaks inside of the radio emission.