Global properties of different types of supernovae are overviewed. A special emphasis is made on the issues of interest both for the supernova and cosmic ray physics such as supernova shock waves and neutrino signal.

We numerically solved the two-dimensional axisymmetric hydrodynamic problem of the explosion of a low-mass neutron star in a circular orbit. In the initial conditions, we assumed a nonuniform density distribution in the space surrounding the collapsed iron core in the form of a stationary toroidal atmosphere that was previously predicted analytically and computed numerically. The configuration of the exploded neutron star itself was modeled by a torus with a circular cross section whose central line almost coincided with its circular orbit. Using an equation of state for the stellar matter and the toroidal atmosphere in which the nuclear statistical equilibrium conditions were satisfied, we performed a series of numerical calculations that showed the propagation of a strong divergent shock wave with a total energy of 0.2x10^51 erg at initial explosion energy release of 1.0x10^51 erg. In our calculations, we rigorously took into account the gravitational interaction, including the attraction from a higher-mass (1.9M_solar) neutron star located at the coordinate origin, in accordance with the rotational explosion mechanism for collapsing supernovae.W e compared in detail our results with previous similar results of asymmetric supernova explosion simulations and concluded that we found a lower limit for the total explosion energy.

A numerical method presented by Imshennik et al. (2002) is used to solve the two dimensional axisymmetric hydrodynamic problem on the formation of a toroidal atmosphere during the collapse of an iron stellar core and outer stellar layers. An evolutionary model from Boyes et al. (1999) with a total mass of $25M_{\odot}$ is used as the initial data for the distribution of thermodynamic quantities in the outer shells of a high-mass star. We analyze in detail the results of three calculations in which the difference mesh and the location of the inner boundary of the computational region are varied. In the initial data, we roughly specify an angular velocity distribution that is actually justified by the final result - the formation of a hydrostatic equilibrium toroidal atmosphere with reasonable total mass, $M^{tot} = (0.117 \div 0.122)M_{\odot}$, and total angular momentum, $J^{tot} = (0.445 \div 0.472) x 10^{50} erg \cdot s$, for the two main calculations. We compare the numerical solution with our previous analytical solution in the form of toroidal atmospheres (Imshennik and Manukovskii 2000). This comparison indicates that they are identical if we take into account the more general and complex equation of state with a nonzero temperature and self-gravitation effects in the atmosphere. Our numerical calculations, first, prove the stability of toroidal atmospheres on characteristic hydrodynamic time scales and, second, show the possibility of sporadic fragmentation of these atmospheres even after a hydrodynamic equilibrium is established. The calculations were carried out under the assumption of equatorial symmetry of the problem and up to relatively long time scales $(\approx 10s)$.

We analyze the behavior of the outer envelope in a massive star during and after the collapse of its iron core into a protoneutron star (PNS) in terms of the equations of one-dimensional spherically symmetric ideal hydrodynamics. The profiles obtained in the studies of the evolution of massive stars up to the final stages of their existence, immediately before a supernova explosion (Boyes et al. 1999), are used as the initial data for the distribution of thermodynamic quantities in the envelope.We use a complex equation of state for matter with allowances made for arbitrary electron degeneracy and relativity, the appearance of electron-positron pairs, the presence of radiation, and the possibility of iron nuclei dissociating into free nucleons and helium nuclei. We performed calculations with the help of a numerical scheme based on Godunov's method. These calculations allowed us to ascertain whether the emersion of the outer envelope in a massive star is possible through the following two mechanisms: first, the decrease in the gravitational mass of the central PNS through neutrino-signal emission and, second, the effect of hot nucleon bubbles, which are most likely formed in the PNS corona, on the envelope emersion. We show that the second mechanism is highly efficient in the range of acceptable masses of the nucleon bubbles ($\leq 0.01M_\odot$) simulated in our hydrodynamic calculations in a rough, spherically symmetric approximation.

An analytic solution has been found in the Roche approximation for the axially symmetric structure of a hydrostatically equilibrium atmosphere of a neutron star produced by collapse. A hydrodynamic (quasione-dimensional) model for the collapse of a rotating iron core in a massive star gives rise to a heterogeneous rotating protoneutron star with an extended atmosphere composed of matter from the outer part of the iron core with differential rotation (Imshennik and Nadyozhin, 1992). The equation of state of a completely degenerate iron gas with an arbitrary degree of relativity is taken for the atmospheric matter. We construct a family of toroidal model atmospheres with total masses $M \approx 0.1 \div 0.2 M_{\odot}$ and total angular momenta $J \approx (1 \div 5.5) \cdot 10^{49} erg \cdot s$, which are acceptable for the outer part of the collapsed iron core, in accordance with the hydrodynamic model, as a function of constant parameters $\omega_{0} and r_{0}$ of the specified differential rotation law $\Omega = \omega_{0}\exp{\Big[-\frac{(r\sin{\Theta})^{2}}{r_{0}^{2}}\Big]}$ in spherical coordinates. The assumed rotation law is also qualitatively consistent with the hydrodynamic model for the collapse of an iron core.

We consider an improved rotational mechanism of the explosion of a collapsing supernova. We show that this mechanism leads to two-stage collapse with a phase difference of \sim 5 h. Based on this model, we attempt a new interpretation of the events in underground neutrino detectors on February 23, 1987, related to the supernova SN 1987A.

A model of a carbon-oxygen (C--O) presupernova core with an initial mass 1.33 M_\odot, an initial carbon mass fraction 0.27, and with an average mass growth-rate 5 x 10^{-7} M_\odot/yr due to accretion in a binary system was evolved from initial central density 10^9 g/cm^3, and temperature 2.05 x 10^8 K through convective core formation and its subsequent expansion to the carbon runaway at the center. The only thermonuclear reaction contained in the equations of evolution and runaway was the carbon burning reaction 12C + 12C with an energy release corresponding to the full transition of carbon and oxygen (with the same rate as carbon) into 56Ni. As a parameter we take \alpha_c - a ratio of a mixing length to the size of the convective zone. In spite of the crude assumptions, we obtained a pattern of the runaway acceptable for the supernova theory with the strong dependence of its duration on \alpha_c. In the variants with large enough values of \alpha_c=4.0 x 10^{-3} and 3.0 x 10^{-3} the fuel combustion occurred from the very beginning as a prompt detonation. In the range of 2.0 x 10^{-3} >= \alpha_c >= 3.0 x 10^{-4} the burning started as a deflagration with excitation of stellar pulsations with growing amplitude. Eventually, the detonation set in, which was activated near the surface layers of the presupernova (with m about 1.33 M_\odot) and penetrated into the star down to the deflagration front. Excitation of model pulsations and formation of a detonation front are described in detail for the variant with \alpha_c=1.0 x 10^{-3}.