Search

We present a new formulation for numerically obtaining axisymmetric equilibrium structures of rotating stars in two spatial dimensions. With a view to apply it to the secular evolution of rotating stars, we base it on the Lagrangian description, i.e., we solve the force-balance equations to find the spatial positions of fluid elements endowed individually with a mass, specific entropy and angular momentum. The system of nonlinear equations obtained by finite-differencing the basic equations are solved with the W4 method, which is a new multi-dimensional root-finding scheme of our own devising. We augment it with a remapping scheme to avoid distortions of the Lagrangian coordinates. In this first one of a series of papers, we will give a detailed description of these methods initially. We then present the results of some test calculations, which include the construction of both rapidly rotating barotropic and baroclinic equilibrium states. We gauge their accuracies quantitatively with some diagnostic quantities as well as via comparisons with the counterparts obtained with an Eulerian code. For a demonstrative purpose, we apply the code to a toy-model cooling calculation of a rotating white dwarf., Comment: 15pages, 14 figures, submitted to MNRAS

We present a new formulation to construct numerically equilibrium configurations of rotating stars in general relativity. Having in mind the application to their quasi static evolutions, we adopt a Lagrangian formulation of our own devising, in which we solve force balance equations to seek for the positions of fluid elements assigned to the grid points, instead of the ordinary Eulerian formulation. Unlike previous works in the literature, we do not employ the first integral of the Euler equation, which is not obtained by an analytic integration in general. We assign a mass, specific angular momentum and entropy to each fluid element in contrast to the previous methods, in which the spatial distribution of the angular velocity or angular momentum is specified. Those distributions are determined after the positions of all fluid elements (or grid points) are derived in our formulation. We solve the large system of algebraic nonlinear equations that are obtained by discretizing the time-independent Euler and Einstein equations in the finite-elements method by using our new multi-dimensional root-finding scheme, named the W4 method. To demonstrate the capability of our new formulation, we construct some rotational configurations both barotropic and baroclinic. We also solve three evolutionary sequences that mimic the cooling, mass-loss, and mass-accretion as simple toy models., Comment: 19 pages, 13 figures

We proposed in this paper a new method, which we named the W4 method, to solve nonlinear equation systems. It may be regarded as an extension of the Newton-Raphson~(NR) method to be used when the method fails. Indeed our method can be applied not only to ordinary problems with non-singular Jacobian matrices but also to problems with singular Jacobians, which essentially all previous methods that employ the inversion of the Jacobian matrix have failed to solve. In this article, we demonstrate that (i) our new scheme can define a non-singular iteration map even for those problems by utilizing the singular value decomposition, (ii) a series of vectors in the new iteration map converges to the right solution under a certain condition, (iii) the standard two-dimensional problems in the literature that no single method proposed so far has been able to solve completely are all solved by our new method., Comment: 13 pages, 4 figures

The progenitor stars of core-collapse supernovae (CCSNe) are asymmetrically fluctuating due to turbulent convections in the late stages of their lives. The progenitor asymmetry at the pre-supernova stage has recently caught the attention as a new ingredient to facilitate shock revival in the delayed neutrino-heating mechanism. In this paper, we investigate the importance of the progenitor asymmetries to shock revival with a semi-analytical approach. Free parameters were chosen such that the time evolution of shock radii and mass accretion rates are compatible with the results of detailed numerical simulations of CCSNe in spherical symmetry. We first estimate the amplitude of asymmetries required for the shock revival by the impulsive change of pre-shock flows in the context of neutrino heating mechanism, and then convert the amplitude to the corresponding amplitude in the pre-supernova phase by taking into account the growth of asymmetries during infall. We apply our model to various types of progenitors and find that the requisite amplitude of pre-supernova asymmetry is roughly three times larger than the prediction by current stellar evolution models unless other additional physical ingredients such as multi-dimensional fluid instabilities and turbulent convections in post-shock flows aid shock revival. We thus conclude that progenitor asymmetries can not trigger the shock revival by the impulsive way but rather play a supplementary role in reality., Comment: Accepted to MNRAS

We propose a new class of method for solving nonlinear systems of equations, which, among other things,has four nice features: (i) it is inspired by the mathematical property of damped oscillators, (ii) it can be regarded as a simple extention to the Newton-Raphson(NR) method, (iii) it has the same local convergence as the NR method does, (iv) it has a significantly wider convergence region or the global convergence than that of the NR method. In this article, we present the evidence of these properties, applying our new method to some examples and comparing it with the NR method., Comment: 17 pages, 7 figures

We investigate axisymmetric steady solutions of (magneto)hydrodynamics equations that describe approximately accretion flows through a standing shock wave and discuss the effects of rotation and magnetic field on the revival of the stalled shock wave in supernova explosions. We develop a new powerful numerical method to calculate the 2-dimensional (2D) steady accretion flows self-consistently. We first confirm the results of preceding papers that there is a critical luminosity of irradiating neutrinos, above which there exists no steady solution in spherical models. If a collapsing star has rotation and/or magnetic field, the accretion flows are no longer spherical owing to the centrifugal force and/or Lorentz force and the critical luminosity is modified.In fact we find that the critical luminosity is reduced by about 50% - 70% for rapid rotations and about 20% - 50% for strong toroidal magnetic fields, depending on the mass accretion rate. These results may be also interpreted as an existence of the critical specific angular momentum or critical magnetic field, above which there exists no steady solution and the standing shock wave will revive for a given combination of mass accretion rate and neutrino luminosity., Comment: 22 pages, 9 figures, to be submitted

The explosion mechanism of core-collapse supernovae has not been fully understood yet but multi-dimensional fluid instabilities such as standing accretion shock instability (SASI) and convection are now believed to be crucial for shock revival. Another multi-dimensional effect that has been recently argued is the asymmetric structures in progenitors, which are induced by violent convections in silicon/oxygen layers that occur before the onset of collapse, as revealed by recent numerical simulations of the last stage of massive star evolutions. Furthermore, it has been also demonstrated numerically that accretions of such non-spherical envelopes could facilitate shock revival. These two multi-dimensional may hence hold a key to successful explosions. In this paper, we performed a linear stability analysis of the standing accretion shock in core-collapse supernovae, taking into account non-spherical, unsteady accretion flows onto the shock to clarify the possible links between the two effects. We found that such pre-shock perturbations can excite the fluid instabilities efficiently and hence help the shock revive in core-collapse supernovae., Comment: 31 pages, 19 figures, 1 table, accepted for publication in The Astrophysical Journal