The toroidal iron atmosphere of a protoneutron star: Numerical solution |
| |
Authors: | V S Imshennik K V Manukovskii M S Popov |
| |
Institution: | (1) Institute for Theoretical and Experimental Physics, ul. Bol’shaya Cheremushkinskaya 25, Moscow, 117259, Russia |
| |
Abstract: | 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⊙ is used as the initial data for the distribution of thermodynamic quantities in the outer shells of a high-mass star. Our computational region includes the outer part of the iron core (without its central part with a mass of 1M⊙ that forms the embryo of a protoneutron star at the preceding stage of the collapse) and the silicon and carbon-oxygen shells with a total mass of (1.8–2.5)M⊙. 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, Mtot=(0.117–0.122)M⊙, and total angular momentum, Jtot=(0.445–0.472)×1050 erg 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 (~10 s). |
| |
Keywords: | plasma astrophysics hydrodynamics shock waves |
本文献已被 SpringerLink 等数据库收录! |
|