Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Past high-energy density laboratory experiments provided insights into the physics of supernovae, supernova remnants, and the destruction of interstellar clouds. In a typical experimental setting, a laser-driven planar blast wave interacts with a compositionally-homogeneous spherical or cylindrical target. In this work we propose a new laboratory platform that accounts for curvature of the impacting shock and density stratification of the target. Both characteristics reflect the conditions expected to exist shortly after a supernova explosion in a close binary system. We provide details of a proposed experimental design (laser drive, target configuration, diagnostic system), optimized to capture the key properties of recent ejecta-companion interaction models. Good qualitative agreement found between our experimental models and their astrophysical counterparts highlights strong potential of the proposed design to probe details of the ejecta-companion interaction for broad classes of objects by means of high energy density laboratory experiments.

The classic single-degenerate model for the progenitors of Type Ia Supernova (SN Ia) predicts that the supernova ejecta should be enriched with solar-like abundance material stripped from the companion star. Spectroscopic observations of normal SNe Ia at late times, however, have not resulted in definite detection of hydrogen. In this Letter, we study line formation in SNe Ia at nebular times using non-LTE spectral modeling. We present, for the first time, multidimensional radiative transfer calculations of SNe Ia with stripped material mixed in the ejecta core, based on hydrodynamical simulations of ejecta-companion interaction. We find that interaction models with main sequence companions produce significant H$\alpha$ emission at late times, ruling out this type of binaries being viable progenitors of SNe Ia. We also predict significant He I line emission at optical and near-infrared wavelengths for both hydrogen-rich or helium-rich material, providing an additional observational probe of stripped ejecta. We produce models with reduced stripped masses and find a more stringent mass limit of $M_{st} \lesssim 1\times 10^{-4} M_\odot$ of stripped companion material for SN 2011fe.

We study the conditions required to produce self-sustained detonations in turbulent, carbon-oxygen degenerate plasma at low densities. We perform a series of three-dimensional hydrodynamic simulations of turbulence driven with various degrees of compressibility. The average conditions in the simulations are representative of models of merging binary white dwarfs. We find that material with very short ignition times is abundant in the case that turbulence is driven compressively. This material forms contiguous structures that persist over many ignition times, and that we identify as prospective detonation kernels. Detailed analysis of prospective kernels reveals that these objects are centrally-condensed and their shape is characterized by low curvature, supportive of self-sustained detonations. The key characteristic of the newly proposed detonation mechanism is thus high degree of compressibility of turbulent drive. The simulated detonation kernels have sizes notably smaller than the spatial resolution of any white dwarf merger simulation performed to date. The resolution required to resolve kernels is 0.1 km. Our results indicate a high probability of detonations in such well-resolved simulations of carbon-oxygen white dwarf mergers. These simulations will likely produce detonations in systems of lower total mass, thus broadening the population of white dwarf binaries capable of producing Type Ia supernovae. Consequently, we expect a downward revision of the lower limit of the total merger mass that is capable of producing a prompt detonation. We review application of the new detonation mechanism to various explosion scenarios of single, Chandrasekhar-mass white dwarfs.

We study supernova ejecta-companion interactions in a sample of realistic semidetached binary systems representative of Type Ia supernova progenitor binaries in a single-degenerate scenario. We model the interaction process with the help of a high-resolution hydrodynamic code assuming cylindrical symmetry. We find that the ejecta hole has a half-opening angle of 40--50$^\circ$ with the density by a factor of 2-4 lower, in good agreement with the previous studies. Quantitative differences from the past results in the amounts and kinematics of the stripped companion material and levels of contamination of the companion with the ejecta material can be explained by different model assumptions and effects due to numerical diffusion.We analyse and, for the first time, provide simulation-based estimates of the amounts and of the thermal characteristics of the shock-heated material responsible for producing a prompt, soft X-ray emission. Besides the shocked ejecta material, considered in the original model by Kasen, we also account for the stripped, shock-heated envelope material of stellar companions, which we predict partially contributes to the prompt emission. The amount of the energy deposited in the envelope is comparable to the energy stored in the ejecta. The total energy budget available for the prompt emission is by a factor of about 2-4 smaller than originally predicted by Kasen. Although the shocked envelope has a higher characteristic temperature than the shocked ejecta, the temperature estimates of the shocked material are in good agreement with the Kasen's model. The hottest shocked plasma is produced in the subgiant companion case.

We study the violent phase of the merger of massive binary white dwarf systems. Our aim is to characterize the conditions for explosive burning to occur, and identify a possible explosion mechanism of Type Ia supernovae. The primary components of our model systems are carbon-oxygen (C/O) white dwarfs, while the secondaries are made either of C/O or of pure helium. We account for tidal effects in the initial conditions in a self-consistent way, and consider initially well-separated systems with slow inspiral rates. We study the merger evolution using an adaptive mesh refinement, reactive, Eulerian code in three dimensions, assuming symmetry across the orbital plane. We use a co-rotating reference frame to minimize the effects of numerical diffusion, and solve for self-gravity using a multi-grid approach. We find a novel detonation mechanism in C/O mergers with massive primaries. Here the detonation occurs in the primary's core and relies on the combined action of tidal heating, accretion heating, and self-heating due to nuclear burning. The exploding structure is compositionally stratified, with a reverse shock formed at the surface of the dense ejecta. The existence of such a shock has not been reported elsewhere. The explosion energy ($1.6\times 10^{51}$ erg) and $^{56}$Ni mass (0.86 M$_\odot$) are consistent with a SN Ia at the bright end of the luminosity distribution, with an approximated decline rate of $\Delta m_{15}(B)\approx 0.99$. Our study does not support double-detonation scenarios in the case of a system with a 0.6 M$_\odot$ helium secondary and a 0.9 M$_\odot$ primary. Although the accreted helium detonates, it fails to ignite carbon at the base of the boundary layer or in the primary's core.

We study the evolution of the collapsing core of a 15 Msun blue supergiant supernova progenitor from the core bounce until 1.5 seconds later. We present a sample of hydrodynamic models parameterized to match the explosion energetics of SN 1987A. We find the spatial model dimensionality to be an important contributing factor in the explosion process. Compared to two-dimensional simulations, our three-dimensional models require lower neutrino luminosities to produce equally energetic explosions. We estimate that the convective engine in our models is 4% more efficient in three dimensions than in two dimensions. We propose that the greater efficiency of the convective engine found in three-dimensional simulations might be due to the larger surface-to-volume ratio of convective plumes, which aids in distributing energy deposited by neutrinos. We do not find evidence of the standing accretion shock instability nor turbulence being a key factor in powering the explosion in our models. Instead, the analysis of the energy transport in the post-shock region reveals characteristics of penetrative convection. The explosion energy decreases dramatically once the resolution is inadequate to capture the morphology of convection on large scales. This shows that the role of dimensionality is secondary to correctly accounting for the basic physics of the explosion. We also analyze information provided by particle tracers embedded in the flow, and find that the unbound material has relatively long residency times in two-dimensional models, while in three dimensions a significant fraction of the explosion energy is carried by particles with relatively short residency times.

We investigate the possibility of generating and studying turbulence in plasma by means of high-energy density laser-driven experiments. Our focus is to create supersonic, self-magnetized turbulence with characteristics that resemble those found in the interstellar medium (ISM). We consider a target made of a spherical core surrounded by a shell made of denser material. The shell is irradiated by a sequence of laser pulses sending inward-propagating shocks that convert the inner core into plasma and create turbulence. In the context of the evolution of the ISM, the shocks play the role of supernova remnant shocks and the core represents the ionized interstellar medium. We consider the effects of both pre-existing and self-generating magnetic fields and study the evolution of the system by means of two-dimensional numerical simulations. We find that the evolution of the turbulent core is generally, subsonic with rms-Mach number $M_t\approx 0.2$. We observe an isotropic, turbulent velocity field with an inertial range power spectra of $P(k)\propto k^{-2.3}$. We account for the effects of self-magnetization and find that the resulting magnetic field has characteristic strength $\approx 3\times 10^{4}$ G. The corresponding plasma beta is $\approx 1\times 10^{4}$--$1\times 10^{5}$, indicating that the magnetic field does not play an important role in the dynamical evolution of the system. The natural extension of this work is to study the system evolution in three-dimensions, with various laser drive configurations, and targets with shells and cores of different masses. The latter modification may help to increase the turbulent intensity and possibly create transonic turbulence. One of the key challenges is to obtain transonic turbulent conditions in a quasi-steady state environment.

(abridged) There exists a substantial disagreement between computer simulation results and high-energy density laboratory experiments of the Rayleigh-Taylor instability Kuranz et al. (2010). We adopt the Braginskii formulation for transport in hot, dense plasma, implement and verify the additional physics modules, and conduct a computational study of a single-mode RTI in two dimensions with various combinations of the newly implemented modules. We find that magnetic fields reach levels on the order of 11 MG in the absence of thermal conduction. We observe denting of the RT spike tip and generation of additional higher order modes as a result of these fields. Contrary to interpretation presented in earlier work Nishiguchi (2002), the additional mode is not generated due to modified anisotropic heat transport effects but due to dynamical effect of self-generated magnetic fields. The main effects of thermal conduction are a reduction of the RT instability growth rate (by about 20% for conditions considered here) and inhibited mixing on small scales. In this case, the maximum self-generated magnetic fields are weaker (approximately 1.7 MG). These self-generated magnetic fields are of very similar strength compared to magnetic fields observed recently in HED laboratory experiments Manuel et al. (2012). We find that thermal conduction plays the dominant role in the evolution of the model RTI system considered. It smears out small-scale structure and reduces the RTI growth rate. This may account for the relatively featureless RT spikes seen in experiments, but does not explain mass extensions observed in experiments. Resistivity and related heat source terms were not included in the present work, but we estimate their impact on RTI as modest and not affecting our main conclusions. Resistive effects will be discussed in detail in the next paper in the series.

Aims: We present neutrino light curves and energy spectra for two representative type Ia supernova explosion models: a pure deflagration and a delayed detonation. Methods: We calculate the neutrino flux from $\beta$ processes using nuclear statistical equilibrium abundances convoluted with approximate neutrino spectra of the individual nuclei and the thermal neutrino spectrum (pair+plasma). Results: Although the two considered thermonuclear supernova explosion scenarios are expected to produce almost identical electromagnetic output, their neutrino signatures appear vastly different, which allow an unambiguous identification of the explosion mechanism: a pure deflagration produces a single peak in the neutrino light curve, while the addition of the second maximum characterizes a delayed-detonation. We identified the following main contributors to the neutrino signal: (1) weak electron neutrino emission from electron captures (in particular on the protons Co55 and Ni56) and numerous beta-active nuclei produced by the thermonuclear flame and/or detonation front, (2) electron antineutrinos from positron captures on neutrons, and (3) the thermal emission from pair annihilation. We estimate that a pure deflagration supernova explosion at a distance of 1 kpc would trigger about 14 events in the future 50 kt liquid scintillator detector and some 19 events in a 0.5 Mt water Cherenkov-type detector. Conclusions: While in contrast to core-collapse supernovae neutrinos carry only a very small fraction of the energy produced in the thermonuclear supernova explosion, the SN Ia neutrino signal provides information that allows us to unambiguously distinguish between different possible explosion scenarios. These studies will become feasible with the next generation of proposed neutrino observatories.

(abridged) We study the hydrodynamic evolution of a non-spherical core-collapse supernova in two spatial dimensions. We find that our model displays a strong tendency to expand toward the pole. We demonstrate that this expansion is a physical property of the low-mode, SASI instability. The SASI leaves behind a large lateral velocity gradient in the post shock layer which affects the evolution for minutes and hours later. This results in a prolate deformation of the ejecta and a fast advection of Ni-rich material from moderate latitudes to the polar regions. This effect might actually be responsible for the global asymmetry of the nickel lines in SN 1987A. The simulations demonstrate that significant radial and lateral motions in the post-shock region, produced by convective overturn and the SASI during the early explosion phase, contribute to the evolution for minutes and hours after shock revival. They lead to both later clump formation, and a significant prolate deformation of the ejecta which are observed even as late as one week after the explosion. As pointed out recently by Kjaer et al., such an ejecta morphology is in good agreement with the observational data of SN 1987A. Systematic future studies are needed to investigate how the SASI-induced late-time lateral expansion depends on the dominant mode of the SASI, and to which extent it is affected by the dimensionality of the simulations. The impact on and importance of the SASI for the distribution of iron group nuclei and the morphology of the young SNR argues for future three-dimensional explosion and post-explosion studies on singularity-free grids that cover the entire sphere. Given the results of our 2D resolution study, present 3D simulations must be regarded as underresolved, and their conclusions must be verified by a proper numerical convergence analysis in three dimensions.

We provide a brief review of the physical processes behind the radiative driving of the winds of OB stars and the Bondi-Hoyle-Lyttleton capture and accretion of a fraction of the stellar wind by a compact object, typically a neutron star, in detached high-mass X-ray binaries (HMXBs). In addition, we describe a program to develop global models of the radiatively-driven photoionized winds and accretion flows of HMXBs, with particular attention to the prototypical system Vela X-1. The models combine XSTAR photoionization calculations, HULLAC emission models appropriate to X-ray photoionized plasmas, improved models of the radiative driving of photoionized winds, FLASH time-dependent adaptive-mesh hydrodynamics calculations, and Monte Carlo radiation transport. We present two- and three-dimensional maps of the density, temperature, velocity, ionization parameter, and emissivity distributions of representative X-ray emission lines, as well as synthetic global Monte Carlo X-ray spectra. Such models help to better constrain the properties of the winds of HMXBs, which bear on such fundamental questions as the long-term evolution of these binaries and the chemical enrichment of the interstellar medium.

In order to determine the circumstances under which isolated SNRs are capable of rising into and enriching the thick disk and galactic halo, simulations of supernova remnants are performed with the FLASH magnetohydrodynamic code. We performed simulations in which the interstellar magnetic field is parallel to or perpendicular to the galactic plane as well as a simulation without a magnetic field. The ambient gas density distribution and gravitational potential are based on observations of our galaxy. We evolve the remnants to ages of roughly 10,000,000 years. For our simulation without a magnetic field, we compare the evolution of the hot bubble's velocity with the velocity evolution calculated from the buoyant and drag accelerations. We found surprisingly small vertical velocities of the hot gas, from which we estimated the drag coefficient to be ten for the non-magnetic simulation. Although we found little buoyant motion of the hot gas during the remnant's lifetime, we found rapid vertical motion of the associated cool dense gas near the end of the remnants life. This motion deformed the remnant into a mushroom cloud structure similar to those found in previous simulations. The simulation in which we have a 4 micro-Gauss magnetic field parallel to the galactic mid-plane shows a dramatically elongated bubble parallel to the magnetic field. The magnetic field pins the supernova remnant preventing it from rising. In the simulation with the 4 micro-Gauss magnetic field perpendicular to the midplane the hot bubble rises more, indicating that having the magnetic field in the same direction as the gravitational force enhances the rise of the bubble.

We describe preliminary results of a global model of the radiatively-driven photoionized wind and accretion flow of the high-mass X-ray binary Vela X-1. The full model combines FLASH hydrodynamic calculations, XSTAR photoionization calculations, HULLAC atomic data, and Monte Carlo radiation transport. We present maps of the density, temperature, velocity, and ionization parameter from a FLASH two-dimensional time-dependent simulation of Vela X-1, as well as maps of the emissivity distributions of the X-ray emission lines.

We develop and demonstrate the methodology of testing multi-dimensional supernova models against observations by studying the properties of one example of the detonation from failed deflagration (DFD) explosion model of thermonuclear supernovae. Using time-dependent multi-dimensional radiative transfer calculations, we generate the synthetic broadband optical light curves, near-infrared light curves, color evolution curves, full spectral time-series, and spectropolarization of the model, as seen from various viewing angles. All model observables are critically evaluated against examples of well-observed, standard Type Ia supernovae (SNe Ia). We explore the consequences of the intrinsic model asphericity by studying the dependence of the model emission on viewing angle, and by quantifying the resulting dispersion in (and internal correlations between) various model observables. These statistical properties of the model are also evaluated against those of the available observational sample of SNe Ia. On the whole, the DFD model shows good agreement with a broad range of SN Ia observations. Certain deficiencies are also apparent, and point to further developments within the basic theoretical framework. We also identify several intriguing orientation effects in the model which suggest ways in which the asphericity of SNe Ia may contribute to their photometric and spectroscopic diversity and, conversely, how the relative homogeneity of SNe Ia constrains the degree of asymmetry allowable in the models. The comprehensive methodology adopted in this work proves an essential component of developing and validating theoretical supernova models, and helps motivate and clearly define future directions in both the modeling and the observation of SNe Ia.

We present a detonating failed deflagration model of Type Ia supernovae. In this model, the thermonuclear explosion of a massive white dwarf follows an off-center deflagration. We conduct a survey of asymmetric ignition configurations initiated at various distances from the stellar center. In all cases studied, we find that only a small amount of stellar fuel is consumed during deflagration phase, no explosion is obtained, and the released energy is mostly wasted on expanding the progenitor. Products of the failed deflagration quickly reach the stellar surface, polluting and strongly disturbing it. These disturbances eventually evolve into small and isolated shock-dominated regions which are rich in fuel. We consider these regions as seeds capable of forming self-sustained detonations that, ultimately, result in the thermonuclear supernova explosion. Preliminary nucleosynthesis results indicate the model supernova ejecta are typically composed of about 0.1-0.25 Msun of silicon group elements, 0.9-1.2 Msun of iron group elements, and are essentially carbon-free. The ejecta have a composite morphology, are chemically stratified, and display a modest amount of intrinsic asymmetry. The innermost layers are slightly egg-shaped with the axis ratio ~1.2-1.3 and dominated by the products of silicon burning. This central region is surrounded by a shell of silicon-group elements. The outermost layers of ejecta are highly inhomogeneous and contain products of incomplete oxygen burning with only small admixture of unburned stellar material. The explosion energies are ~1.3-1.5 10^51 erg.

The thermonuclear explosion of a massive white dwarf in a Type Ia supernova explosion is characterized by vastly disparate spatial and temporal scales. The extreme dynamic range inherent to the problem prevents the use of direct numerical simulation and forces modelers to resort to subgrid models to describe physical processes taking place on unresolved scales. We consider the evolution of a model thermonuclear flame in a constant gravitational field on a periodic domain. The gravitational acceleration is aligned with the overall direction of the flame propagation, making the flame surface subject to the Rayleigh-Taylor instability. The flame evolution is followed through an extended initial transient phase well into the steady-state regime. The properties of the evolution of flame surface are examined. We confirm the form of the governing equation of the evolution suggested by Khokhlov (1995). The mechanism of vorticity production and the interaction between vortices and the flame surface are discussed. The results of our investigation provide the bases for revising and extending previous subgrid-scale model.

AIMS. We study and discuss the time-dependent X-ray emission predicted by hydrodynamic modeling of the interaction of a SNR shock wave with an interstellar gas cloud. The scope includes: 1) to study the correspondence between modeled and X-ray emitting structures, 2) to explore two different physical regimes in which either thermal conduction or radiative cooling plays a dominant role, and 3) to investigate the effects of the physical processes at work on the emission of the shocked cloud in the two different regimes. METHODS. We use a detailed hydrodynamic model, including thermal conduction and radiation, and explore two cases characterized by different Mach numbers of the primary shock: M = 30 in which the cloud dynamics is dominated by radiative cooling and M = 50 dominated by thermal conduction. From the simulations, we synthesize the expected X-ray emission, using available spectral codes. RESULTS. The morphology of the X-ray emitting structures is significantly different from that of the flow structures originating from the shock-cloud interaction. The hydrodynamic instabilities are never clearly visible in the X-ray band. Shocked clouds are preferentially visible during the early phases of their evolution. Thermal conduction and radiative cooling lead to two different phases of the shocked cloud: a cold cooling dominated core emitting at low energies and a hot thermally conducting corona emitting in the X-ray band. The thermal conduction makes the X-ray image of the cloud smaller, more diffuse, and shorter-lived than that observed when thermal conduction is neglected.

We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in a protoplanetary disc during a few hundred orbital periods. We compare the surface density contours, potential vorticity and smoothed radial profiles at several times. The disc mass and gravitational torque time evolution are analyzed with high temporal resolution. There is overall consistency between the codes. The density profiles agree within about 5% for the Eulerian simulations while the SPH results predict the correct shape of the gap although have less resolution in the low density regions and weaker planetary wakes. The disc masses after 200 orbital periods agree within 10%. The spread is larger in the tidal torques acting on the planet which agree within a factor 2 at the end of the simulation. In the Neptune case the dispersion in the torques is greater than for Jupiter, possibly owing to the contribution from the not completely cleared region close to the planet.

Two-dimensional simulations of strongly anisotropic supernova explosions of a nonrotating 15 solar mass blue supergiant progenitor are presented, which follow the hydrodynamic evolution from times shortly after shock formation until hours later. It is shown that explosions which around the time of shock revival are dominated by low-order unstable modes (i.e. by a superposition of the l=2 and l=1 modes, in which the former is strongest), are consistent with all major observational features of SN 1987A, in contrast to models which show high-order mode perturbations only and were published in earlier work. Among other items, the low-mode models exhibit final iron-group velocities of up to 3300 km/s, strong mixing at the He/H composition interface, with hydrogen being mixed downward in velocity space to only 500 km/s, and a final prolate anisotropy of the ejecta with a major to minor axis ratio of about 1.6. The success of low-mode explosions with an energy of about 2x10**51 erg to reproduce these observed features is based on two effects: the (by 40%) larger initial maximum velocities of metal-rich clumps compared to our high-mode models, and the initial global deformation of the shock. The latter triggers the growth of a strong Richtmyer-Meshkov instability at the He/H interface that results in a global anisotropy of the inner ejecta at late times (i.e. t > 10000 s), although the shock itself has long become spherical by then. The simulations suggest a coherent picture, which explains the observational data of SN 1987A within the framework of the neutrino-driven explosion mechanism using a minimal set of assumptions. It is therefore argued that other paradigms, which are based on (more) controversial physics, may not be required to explain this event. (abbreviated)

We model the hydrodynamic interaction of a shock wave of an evolved supernova remnant with a small interstellar gas cloud like the ones observed in the Cygnus loop and in the Vela SNR. We investigate the interplay between radiative cooling and thermal conduction during cloud evolution and their effect on the mass and energy exchange between the cloud and the surrounding medium. Through the study of two cases characterized by different Mach numbers of the primary shock (M = 30 and 50, corresponding to a post-shock temperature $T\approx 1.7\times 10^6$ K and $\approx 4.7\times 10^6$ K, respectively), we explore two very different physical regimes: for M = 30, the radiative losses dominate the evolution of the shocked cloud which fragments into cold, dense, and compact filaments surrounded by a hot corona which is ablated by the thermal conduction; instead, for M = 50, the thermal conduction dominates the evolution of the shocked cloud, which evaporates in a few dynamical time-scales. In both cases we find that the thermal conduction is very effective in suppressing the hydrodynamic instabilities that would develop at the cloud boundaries.

We present our first nucleosynthesis results from a numerical simulation of the thermonuclear disruption of a static cold Chandrasekhar-mass C/O white dwarf. The two-dimensional simulation was performed with an adaptive-mesh Eulerian hydrodynamics code, FLASH, that uses as a flame capturing scheme the evolution of a passive scaler. To compute the isotopic yields and their velocity distribution, 10,000 massless tracer particles are embedded in the star. The particles are advected along streamlines and provide a Lagrangian description of the explosion. We briefly describe our verification tests and preliminary results from post-processing the particle trajectories with a modest (214 isotopes) reaction network.

We present an extension of the Piecewise Parabolic Method to special relativistic fluid dynamics in multidimensions. The scheme is conservative, dimensionally unsplit, and suitable for a general equation of state. Temporal evolution is second-order accurate and employs characteristic projection operators; spatial interpolation is piece-wise parabolic making the scheme third-order accurate in smooth regions of the flow away from discontinuities. The algorithm is written for a general system of orthogonal curvilinear coordinates and can be used for computations in non-cartesian geometries. A non-linear iterative Riemann solver based on the two-shock approximation is used in flux calculation. In this approximation, an initial discontinuity decays into a set of discontinuous waves only implying that, in particular, rarefaction waves are treated as flow discontinuities. We also present a new and simple equation of state which approximates the exact result for the relativistic perfect gas with high accuracy. The strength of the new method is demonstrated in a series of numerical tests and more complex simulations in one, two and three dimensions.

We have developed a three-dimensional radiative transfer method designed specifically for use with parallel adaptive mesh refinement hydrodynamics codes. This new algorithm, which we call hybrid characteristics, introduces a novel form of ray tracing that can neither be classified as long, nor as short characteristics, but which applies the underlying principles, i.e. efficient execution through interpolation and parallelizability, of both. Primary applications of the hybrid characteristics method are radiation hydrodynamics problems that take into account the effects of photoionization and heating due to point sources of radiation. The method is implemented in the hydrodynamics package FLASH. The ionization, heating, and cooling processes are modelled using the DORIC ionization package. Upon comparison with the long characteristics method, we find that our method calculates the column density with a similarly high accuracy and produces sharp and well defined shadows. We show the quality of the new algorithm in an application to the photoevaporation of multiple over-dense clumps. We present several test problems demonstrating the feasibility of our method for performing high resolution three-dimensional radiation hydrodynamics calculations that span a large range of scales. Initial performance tests show that the ray tracing part of our method takes less time to execute than other parts of the calculation (e.g. hydrodynamics and adaptive mesh refinement), and that a high degree of efficiency is obtained in parallel execution. Although the hybrid characteristics method is developed for problems involving photoionization due to point sources, the algorithm can be easily adapted to the case of more general radiation fields.

We consider some of the spectral and polarimetric signatures of the gravitational confined detonation scenario for Type Ia supernova explosions. In this model, material produced by an off-center deflagration (which itself fails to produce the explosion) forms a metal-rich atmosphere above the white dwarf surface. Using hydrodynamical simulations, we show that this atmosphere is compressed and accelerated during the subsequent interaction with the supernova ejecta. This leads ultimately to the formation of a high-velocity pancake of metal-rich material that is geometrically detached from the bulk of the ejecta. When observed at the epochs near maximum light, this absorbing pancake produces a highly blueshifted and polarized calcium IR triplet absorption feature similar to that observed in several Type~Ia supernovae. We discuss the orientation effects present in our model and contrast them to those expected in other supernova explosion models. We propose that a large sample of spectropolarimetric observations can be used to critically evaluate the different theoretical scenarios.

Two- and three-dimensional simulations demonstrate that hydrodynamic instabilities can lead to low-mode (l=1,2) asymmetries of the fluid flow in the neutrino-heated layer behind the supernova shock. This provides a natural explanation for aspherical mass ejection and for pulsar recoil velocities even in excess of 1000 km/s. We propose that the bimodality of the pulsar velocity distribution might be a consequence of a dominant l=1 mode in case of the fast component, while higher-mode anisotropy characterizes the postshock flow and SN ejecta during the birth of the slow neutron stars. We argue that the observed large asymmetries of supernovae and the measured high velocities of young pulsars therefore do not imply rapid rotation of the iron core of the progenitor star, nor do they require strong magnetic fields to play a crucial role in the explosion. Anisotropic neutrino emission from accretion contributes to the neutron star acceleration on a minor level, and pulsar kicks do not make a good case for non-standard neutrino physics in the nascent neutron star.

As a component of the Flash Center's validation program, we compare FLASH simulation results with experimental results from Los Alamos National Laboratory. The flow of interest involves the lateral interaction between a planar Ma=1.2 shock wave with a cylinder of gaseous sulfur hexafluoride (SF_6) in air, and in particular the development of primary and secondary instabilities after the passage of the shock. While the overall evolution of the flow is comparable in the simulations and experiments, small-scale features are difficult to match. We focus on the sensitivity of numerical results to simulation parameters.

We show by two-dimensional and first three-dimensional simulations of neutrino-driven supernova explosions that low (l=1,2) modes can dominate the flow pattern in the convective postshock region on timescales of hundreds of milliseconds after core bounce. This can lead to large global anisotropy of the supernova explosion and pulsar kicks in excess of 500 km/s.

We present a new mechanism for Type Ia supernova explosions in massive white dwarfs. The proposed scenario follows from relaxing the assumption of symmetry in the model and involves a detonation created in an unconfined environment. The explosion begins with an essentially central ignition of stellar material initiating a deflagration. This deflagration results in the formation of a buoyantly-driven bubble of hot material that reaches the stellar surface at supersonic speeds. The bubble breakout forms a strong pressure wave that laterally accelerates fuel-rich outer stellar layers. This material, confined by gravity to the white dwarf, races along the stellar surface and is focused at the location opposite to the point of the bubble breakout. These streams of nuclear fuel carry enough mass and energy to trigger a detonation just above the stellar surface. The flow conditions at that moment support a detonation that will incinerate the white dwarf and result in an energetic explosion. The stellar expansion following the deflagration redistributes stellar mass in a way that ensures production of intermediate mass and iron group elements consistent with observations. The ejecta will have a strongly layered structure with a mild amount of asymmetry following from the early deflagration phase. This asymmetry, combined with the amount of stellar expansion determined by details of the evolution (principally the energetics of deflagration, timing of detonation, and structure of the progenitor), can be expected to create a family of mildly diverse Type Ia supernova explosions.

We present the first high-resolution three-dimensional simulations of the deflagration phase of Type Ia supernovae that treat the entire massive white dwarf. We report the results of simulations in which ignition of the nuclear burning occurs slightly off-center. The subsequent evolution of the nuclear burning is surprisingly asymmetric with a growing bubble of hot ash rapidly rising to the stellar surface. Upon reaching the surface, the mass of burned material is $\approx 0.075 M_\sun$ and the kinetic energy is $4.3 \times 10^{49}$ ergs. The velocity of the top of the rising bubble approaches 8000 km s$^{-1}$. The amount of the asymmetry found in the model offers a natural explanation for the observed diversity in Type Ia supernovae. Our study strongly disfavors the classic central-ignition pure deflagration scenario by showing that the result is highly sensitive to details of the initial conditions.

We compare the results of two--dimensional simulations to experimental data obtained at Los Alamos National Laboratory in order to validate the FLASH code. FLASH is a multi--physics, block--structured adaptive mesh refinement code for studying compressible, reactive flows in various astrophysical environments. The experiment involves the lateral interaction between a planar Ma=1.2 shock wave with a cylinder of gaseous sulfur hexafluoride (SF$_6$) in air. The development of primary and secondary flow instabilities after the passage of the shock, as observed in the experiments and numerical simulations, are reviewed and compared. The sensitivity of numerical results to several simulation parameters are examined. Computed and experimentally measured velocity fields are compared. Motivation for experimental data in planes parallel to the cylinder axis is provided by a speculative three--dimensional simulation.

Assuming that the neutrino luminosity from the neutron star core is sufficiently high to drive supernova explosions by the neutrino-heating mechanism, we show that low-mode (l = 1, 2) convection can develop from random seed perturbations behind the shock. A slow onset of the explosion is crucial, requiring the core luminosity to vary slowly with time, in contrast to the burst-like exponential decay assumed in previous work. Gravitational and hydrodynamic forces by the globally asymmetric supernova ejecta were found to accelerate the remnant neutron star on a timescale of more than a second to velocities above 500 km/s, in agreement with observed pulsar proper motions.

Two-dimensional simulations of a Type II and a Type Ib-like supernova explosion are presented that encompass shock revival by neutrino heating, neutrino-driven convection, explosive nucleosynthesis, the growth of Rayleigh-Taylor instabilities, and the propagation of newly formed metal clumps through the exploding star. In both cases we find very high Ni56 velocities of 17000 km/s shortly after shock-revival, and a complete fragmentation of the progenitor's metal core within the first few minutes after core bounce, due to the growth of Rayleigh-Taylor instabilities at the Si/O and (C+O)/He composition interfaces. This leads to the formation of high-velocity, metal-rich clumps which eventually decouple from the flow and move ballistically through the ejecta. Maximum final metal velocities of 3500-5500 km/s and 1200 km/s are obtained for the Type Ib model and the Type II model, respectively. The low maximum metal velocities in the Type II model, which are significantly smaller than those observed in SN 1987A, are due to the massive hydrogen envelope of the progenitor. The envelope forces the supernova shock to decelerate strongly, leaving behind a reverse shock below the He/H interface, which interacts with the clumps and slows them down significantly. This reverse shock is absent in the Type Ib-like model. The latter is in fairly good agreement with observations of Type Ib supernovae. In addition, in this case the pattern of clump formation in the ejecta is correlated with the convective pattern prevailing during the shock-revival phase. This might be used to deduce observational constraints for the dynamics during this early phase of the evolution, and the role of neutrino heating in initiating the explosion.

One of the central problems in supernova theory is the question how massive stars explode. Understanding the physical processes that drive the explosion is crucial for linking the stellar progenitors to the final remnants and for predicting observable properties like explosion energies, neutron star and black hole masses, nucleosynthetic yields, explosion anisotropies, and pulsar kicks. Here we review different suggestions for the explosion mechanism and discuss the constraints that can or cannot be deduced from observations. The most refined and detailed hydrodynamical models are available for neutrino-driven supernova explosions, although the viability of this mechanism has yet to be demonstrated convincingly. Since spherical models do not explode, the hope rests on the helpful effects of convection inside the nascent neutron star and in the neutrino-heated region behind the stalled shock. We present the first two-dimensional simulations of these processes which have been performed with a Boltzmann solver for the neutrino transport and a state-of-the-art description of neutrino-matter interactions. Our most complete models fail, but convection brings them encouragingly close to a success. An explosion could be obtained by just a minor modification of the neutrino transport, in which case the exploding model fulfills important requirements from observations (abridged).

The rapidly growing base of observational data for supernova explosions of massive stars demands theoretical explanations. Central of these is a self-consistent model for the physical mechanism that provides the energy to start and drive the disruption of the star. We give arguments why the delayed neutrino-heating mechanism should still be regarded as the standard paradigm to explain most explosions of massive stars and show how large-scale and even global asymmetries can result as a natural consequence of convective overturn in the neutrino-heating region behind the supernova shock. Since the explosion is a threshold phenomenon and depends sensitively on the efficiency of the energy transfer by neutrinos, even relatively minor differences in numerical simulations can matter on the secular timescale of the delayed mechanism. To enhance this point, we present some results of recent one- and two-dimensional computations, which we have performed with a Boltzmann solver for the neutrino transport and a state-of-the-art description of neutrino-matter interactions. Although our most complete models fail to explode, the simulations demonstrate that one is encouragingly close to the critical threshold because a modest variation of the neutrino transport in combination with postshock convection leads to a weak neutrino-driven explosion with properties that fulfill important requirements from observations.

We look in detail at the process of mapping an astrophysical initial model from a stellar evolution code onto the computational grid of an explicit, Godunov type code while maintaining hydrostatic equilibrium. This mapping process is common in astrophysical simulations, when it is necessary to follow short-timescale dynamics after a period of long timescale buildup. We look at the effects of spatial resolution, boundary conditions, the treatment of the gravitational source terms in the hydrodynamics solver, and the initialization process itself. We conclude with a summary detailing the mapping process that yields the lowest ambient velocities in the mapped model.

We present a case study of validating an astrophysical simulation code. Our study focuses on validating FLASH, a parallel, adaptive-mesh hydrodynamics code for studying the compressible, reactive flows found in many astrophysical environments. We describe the astrophysics problems of interest and the challenges associated with simulating these problems. We describe methodology and discuss solutions to difficulties encountered in verification and validation. We describe verification tests regularly administered to the code, present the results of new verification tests, and outline a method for testing general equations of state. We present the results of two validation tests in which we compared simulations to experimental data. The first is of a laser-driven shock propagating through a multi-layer target, a configuration subject to both Rayleigh-Taylor and Richtmyer-Meshkov instabilities. The second test is a classic Rayleigh-Taylor instability, where a heavy fluid is supported against the force of gravity by a light fluid. Our simulations of the multi-layer target experiments showed good agreement with the experimental results, but our simulations of the Rayleigh-Taylor instability did not agree well with the experimental results. We discuss our findings and present results of additional simulations undertaken to further investigate the Rayleigh-Taylor instability.

A hydrodynamic interaction between an old supernova remnant and a massive black hole is studied with the help of a high-resolution adaptive mesh refinement code. A series of multidimensional numerical models is obtained for a fixed distance between supernova explosion site and the black hole but different densities of the ambient medium. The interaction between the shell of the supernova remnant and the black hole results in a rapid increase of the emission from the vicinity of the black hole. Light curves of the gas heated by gravitational compression and hydrodynamical interactions are obtained assuming optically thin conditions. It is shown that the duration of the event and the peak luminosity critically depend on the density of the ambient medium. The model is applied to the data obtained by the Chandra satellite for the Galactic center.

We use hydrodynamic simulations with adaptive grid refinement to study the dependence of hot gas flows in X-ray luminous giant elliptical galaxies on the efficiency of heat supply to the gas. We consider a number of potential heating mechanisms including Type Ia supernovae and sporadic nuclear activity of a central supermassive black hole. As a starting point for this research we use an equilibrium hydrostatic recycling model (Kritsuk 1996). We show that a compact cooling inflow develops, if the heating is slightly insufficient to counterbalance radiative cooling of the hot gas in the central few kiloparsecs. An excessive heating in the centre, instead, drives a convectively unstable outflow. We model the onset of the instability and a quasi-steady convective regime in the core of the galaxy in two-dimensions assuming axial symmetry. Provided the power of net energy supply in the core is not too high, the convection remains subsonic. The convective pattern is dominated by buoyancy driven large-scale mushroom-like structures. Unlike in the case of a cooling inflow, the X-ray surface brightness of an (on average) isentropic convective core does not display a sharp maximum at the centre. A hybrid model, which combines a subsonic peripheral cooling inflow with an inner convective core, appears to be stable. We also discuss observational implications of these results.

After decades of one-dimensional nucleosynthesis calculations, the growth of computational resources has meanwhile reached a level, which for the first time allows astrophysicists to consider performing routinely realistic multidimensional nucleosynthesis calculations in explosive and, to some extent, also in non-explosive environments. In the present contribution we attempt to give a short overview of the physical and numerical problems which are encountered in these simulations. In addition, we assess the accuracy that can be currently achieved in the computation of nucleosynthetic yields, using multidimensional simulations of core collapse supernovae as an example.

Implementation details and test cases of a newly developed hydrodynamic code, AMRA, are presented. The numerical scheme exploits the adaptive mesh refinement technique coupled to modern high-resolution schemes which are suitable for relativistic and non-relativistic flows. Various physical processes are incorporated using the operator splitting approach, and include self-gravity, nuclear burning, physical viscosity, implicit and explicit schemes for conductive transport, simplified photoionization, and radiative losses from an optically thin plasma. Several aspects related to the accuracy and stability of the scheme are discussed in the context of hydrodynamic and astrophysical flows.

Motivated by observations of supernova SN 1987A, various authors have simulated Rayleigh-Taylor (RT) instabilities in the envelopes of core collapse supernovae (for a review, see Mueller 1998). The non-radial motion found in these simulations qualitatively agreed with observations in SN 1987A, but failed to explain the extent of mixing of newly synthesized 56Ni quantitatively. Here we present results of a 2D hydrodynamic simulation which re-addresses this failure and covers the entire evolution of the first 5 hours after core bounce.

The advantages and disadvantages of two approaches to astrophysical hydrodynamics, Smoothed Particle Hydrodynamics and Adaptive Mesh Refinement, are briefly discussed together with some current problems of computational hydrodynamics.

New two-dimensional, high-resolution calculations of a core collapse supernova in a 15 Msol blue supergiant are presented, which cover the entire evolution from shock revival until the first few hours of the explosion. Explosive nucleosynthesis, its dependence upon convective mixing during the first second of the evolution and the growth of Rayleigh-Taylor instabilities at the composition interfaces of the progenitor star are all modeled consistently and allow for a comparison with observational data. We confirm our earlier findings, that the perturbations induced by neutrino driven convection are sufficiently strong to seed large-scale Rayleigh-Taylor mixing and to destroy the onion-shell structure of the stellar He-core. As in our earlier calculations, the strong deceleration of the nickel clumps in the layers adjacent to the He/H interface suggests that the high velocities of iron-group elements observed in SN 1987A cannot be explained on the basis of currently favored progenitor models. Possible solutions to this dilemma and the implications of the mixing for type Ib explosions are briefly discussed.

A problem of mass flow in the immediate vicinity of a planet embedded in a protoplanetary disk is studied numerically in two dimensions. Large differences in temporal and spatial scales involved suggest that a specialized discretization method for solution of hydrodynamical equations may offer great savings in computational resources, and can make extensive parameter studies feasible. Preliminary results obtained with help of Adaptive Mesh Refinement technique and high-order explicit Eulerian solver are presented. This combination of numerical techniques appears to be an excellent tool which allows for direct simulations of mass flow in vicinity of the accretor at moderate computational cost. In particular, it is possible to resolve the surface of the planet and to model the process of planet growth with minimal set of assumptions. Some issues related to visualization of the results and future prospects are discussed briefly.

High-resolution two-dimensional simulations were performed for the first five minutes of the evolution of a core collapse supernova explosion in a 15 solar mass blue supergiant progenitor. The computations start shortly after bounce and include neutrino-matter interactions by using a light-bulb approximation for the neutrinos, and a treatment of the nucleosynthesis due to explosive silicon and oxygen burning. We find that newly formed iron-group elements are distributed throughout the inner half of the helium core by Rayleigh-Taylor instabilities at the Ni+Si/O and C+O/He interfaces, seeded by convective overturn during the early stages of the explosion. Fast moving nickel mushrooms with velocities up to about 4000 km/s are observed. This offers a natural explanation for the mixing required in light curve and spectral synthesis studies of Type Ib explosions. A continuation of the calculations to later times, however, indicates that the iron velocities observed in SN 1987 A cannot be reproduced because of a strong deceleration of the clumps in the dense shell left behind by the shock at the He/H interface.

We present results of high-resolution two-dimensional simulations which follow the first five minutes of a core collapse supernova explosion in a 15 solar mass blue supergiant progenitor. The computations start shortly after core bounce and include neutrino-matter interactions by using a light-bulb approximation for the neutrinos, and a treatment of the nucleosynthesis due to explosive silicon and oxygen burning. We find that newly formed iron-group elements are distributed throughout a significant fraction of the stellar helium core by the concerted action of convective and Rayleigh-Taylor instabilities. Fast moving nickel mushrooms with velocities up to 4000 km/s are observed. A continuation of the calculations to later times, however, indicates, that the iron velocities observed in SN 1987 A cannot be reproduced due to a strong deceleration of the clumps during their interaction with the dense shell left behind by the shock at the He/H interface. Therefore, we cannot confirm the claim that convective "premixing" of the nickel in the early phases of the explosion solves the problem of the high iron velocities.

We present a numerical study of the cosmic density vs. velocity divergence relation (DVDR) in the mildly non-linear regime. We approximate the dark matter as a non-relativistic pressureless fluid, and solve its equations of motion on a grid fixed in comoving coordinates. Unlike N-body schemes, this method yields directly the volume-averaged velocity field. The results of our simulations are compared with the predictions of the third-order perturbation theory (3PT) for the DVDR. We investigate both the mean `forward' relation (density in terms of velocity divergence) and the mean `inverse' relation (velocity divergence in terms of density), with emphasis on the latter. On scales larger than about 20 megaparsecs, our code recovers the predictions of 3PT remarkably well, significantly better than recent N-body simulations. On scales of a few megaparsecs, the DVDR predicted by 3PT differs slightly from the simulated one. In particular, approximating the inverse DVDR by a third-order polynomial turns out to be a poor fit. We propose a simple analytical description of the inverse relation, which works well for mildly non-linear scales.

Simple modifications for higher-order Godunov-type difference schemes are presented which allow for accurate advection of multi-fluid flows in hydrodynamic simulations. The constraint that the sum of all mass fractions has to be equal to one in every computational zone throughout the simulation is fulfilled by renormalizing the mass fractions during the advection step. The proposed modification is appropriate for any difference scheme written in conservation form. Unlike other commonly used methods it does not violate the conservative character of the advection method. A new steepening mechanism, which is based on modification of interpolation profiles, is used to reduce numerical diffusion across composition discontinuities. Additional procedures are described, which are necessary to enforce monotonicity. Several numerical experiments are presented which demonstrate the capability of our Consistent Multi-fluid Advection (CMA) method in case of smooth and discontinuous distributions of fluid phases and under different hydrodynamic conditions. It is shown that due to the reduced diffusivity of the proposed scheme the abundance of some heavy elements obtained from hydrodynamic simulations of type II supernova explosions can change by a factor of a few in the most extreme cases.

In recent years the technique of Adaptive Mesh Refinement (AMR) has been successfully applied to several astrophysical problems. This method should also allow a consistent multi-dimensional modelling of supernova explosions. In this contribution we address some of the computational difficulties encountered when trying to apply the AMR method to explosive nucleosynthesis and supernova envelope ejection.

(shortened version) We show that many observations of W44, a supernova remnant in the galactic plane at a distance of about 2500 pc, are remarkably consistent with the simplest realistic model. The model remnant is evolving in a smooth ambient medium of fairly high density, about 6 cm^-3 on average, with a substantial density gradient. At the observed time it has an age of about 20,000 years, consistent with the age of the associated pulsar, and a radius of 11 to 13 pc. Over most of the outer surface, radiative cooling has become important in the post shock gas; on the denser end there has been sufficient compression of the cooled gas to develop a very thin dense half shell of about 450 M_sun, supported against further compression by nonthermal pressure. The half shell has an expansion velocity of about 150 km s^-1, and is bounded on the outer surface by a radiative shock with that speed. We provide several analytic tools for the assembly of models of this type. We review the early evolution and shell formation analyses and their generalizations to evolution in a density gradient. We also calculate the density and temperature that should be present in the hot interior of a remnant with thermal conduction. We supply the van der Laan mechanism in a particularly useful form for the calculation of radio continuum from radiative remnants. Finally, we demonstrate a simple technique for estimating the optical emission expected. These tools are employed to choose parameters of models which we then explore with our 1d and 2d hydrocodes, providing, respectively, the detailed x-ray spectra and dynamical characteristics.

We present simulations of relativistic jets propagating in a nonuniform medium. Specifically, we study the bending of jets propagating obliquely to the vector of the density-gradient. Our results are applied to the NGC 4258, where such a medium is assumed to be provided by the atmosphere of the sub-parsec accretion disk tilted with respect to the original direction of the jet propagation. As a result, the jet is bent on a scale comparable to the density scaleheight of the disk atmosphere. The magnitude of the bending effect is found to be largest for light jets with low Lorentz factors. The predicted direction of bending is consistent with the observations.

Here we discuss the interaction of massive stars with their parental molecular clouds. A summary of the dynamical evolution of HII regions and wind-driven bubbles in high-pressure cloud cores is given. Both ultracompact HII regions and ultracompact wind-driven bubbles can reach pressure equilibrium with their surrounding medium. The structures stall their expansion and become static and, as long as the ionization sources and the ambient densities remain about constant, the resulting regions are stable and long lived. For cases with negative density gradients, and depending on the density distribution, some regions never reach the static equilibrium condition. For power-law density stratifications, $\rho \propto r^{-w}$, the properties of the evolution depend on a critical exponent, $w_{crit}$, above which the ionization front cannot be slowed down by recombinations or new ionizations, and the cloud becomes fully ionized. This critical exponent is $w_{crit}=3/2$ during the expansion phase. For $w>3/2 $ the gas expands supersonically into the surrounding ionized medium, and there are two regimes separated by $w=3$. For $3/2<w\leq 3$, the slow regime, the inner region drives a weak shock moving with almost constant velocity through the cloud. For $w>3$, the fast regime, the shock becomes strong and accelerates with time. Finally, the evolution of slow winds in highly pressurized region is described briefly.

We present the main features in the evolution of the gas pressure in star forming regions, from the formation of the parental cloud to the moment when the region is pressurized by interacting stellar winds. The main processes for cloud destruction and the self-limiting properties of star formation are described. The high star forming efficiency in nuclear starbursts is a consequence of the high gas pressures. The evolution of slow winds in highly pressurized region is also sketched.

We present a new Eulerian code able to follow the evolution of large-scale structures in the Universe in the weakly nonlinear regime. We compare test results with a N-body code and analytical results.

We discuss some typical problems related to numerical hydrodynamics of a dense interstellar medium. A newly developed hydrodynamical code based on adaptive mesh refinement technique is presented and applied to simulate the evolution of a supernova remnant in a high-density medium. Advantages of this new approach in comparison to commonly used hydrodynamical methods are briefly discussed.

We examine the evolution of inhomogeneities (fragments) of supernova ejecta in compact supernova remnants by means of hydrodynamical modeling and simplified analytical calculations. Under the influence of intense post-shock cooling the fragments become strongly compressed as they traverse the hot shocked region between the reverse and outer shocks of the remnant. We find that the most likely outcome of the interaction of fragments with the reverse shock and the hot shocked region is their disruption resulting in generation of secondary fragments. Secondary fragments arriving at the thin and dense outer shell of the remnant give rise to brief X-ray flashes. Under suitable conditions the primary fragments may traverse the hot shocked region without being completely destroyed, to eventually reach the outer shell as dense, elongated structures. Collisions of such fragments with the shell are likely to give rise to powerful X-ray flares.

Evolution of fast, radiative shocks in high density medium is presented. Ionizing spectra and approximate broad band light curves of the shocked gas are calculated. Emergent shock spectra, as seen by a distant observer, are obtained from photoionization models. The emergent spectra have a power-law shape $F_{\nu}\propto{\nu}^{-\alpha}$ with mean spectral index $\alpha\sim0.6-1.0$ in the energy range $0.01-10$ keV, and have a high-energy cutoff corresponding to the original shock velocity. It is shown that the models exhibit promising features that may account for some photometric and spectral properties of Active Galactic Nuclei.