Computational Models of Stellar Collapse and Core-Collapse Supernovae

C. D. Ott,E. Schnetter,A. Burrows,E. Livne,E. O'Connor,F. Loeffler
DOI: https://doi.org/10.1088/1742-6596/180/1/012022
2009-07-23
Abstract:Core-collapse supernovae are among Nature's most energetic events. They mark the end of massive star evolution and pollute the interstellar medium with the life-enabling ashes of thermonuclear burning. Despite their importance for the evolution of galaxies and life in the universe, the details of the core-collapse supernova explosion mechanism remain in the dark and pose a daunting computational challenge. We outline the multi-dimensional, multi-scale, and multi-physics nature of the core-collapse supernova problem and discuss computational strategies and requirements for its solution. Specifically, we highlight the axisymmetric (2D) radiation-MHD code VULCAN/2D and present results obtained from the first full-2D angle-dependent neutrino radiation-hydrodynamics simulations of the post-core-bounce supernova evolution. We then go on to discuss the new code Zelmani which is based on the open-source HPC Cactus framework and provides a scalable AMR approach for 3D fully general-relativistic modeling of stellar collapse, core-collapse supernovae and black hole formation on current and future massively-parallel HPC systems. We show Zelmani's scaling properties to more than 16,000 compute cores and discuss first 3D general-relativistic core-collapse results.
High Energy Astrophysical Phenomena,Solar and Stellar Astrophysics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the detailed understanding of the explosion mechanism of Core - Collapse Supernovae (CCSN). Specifically, the paper focuses on the following aspects: 1. **Multi - dimensional, multi - scale, and multi - physical properties**: The CCSN problem is essentially a problem of multi - dimensional, multi - scale, and multi - physical processes. This means that it is necessary to simulate various physical phenomena from microscopic particle interactions to macroscopic hydrodynamics at different spatial scales. 2. **Computational strategies and requirements**: The paper discusses the computational strategies and requirements for solving the CCSN problem, especially how to implement these strategies on existing high - performance computing (HPC) systems. This includes using the axisymmetric (2D) radiation magnetohydrodynamics (MHD) code VULCAN/2D for full 2D angle - dependent neutrino radiation - hydrodynamics simulations, and a new code Zelmani based on the open - source HPCCactus framework, which provides an extensible adaptive mesh refinement (AMR) method for 3D fully general - relativistic modeling. 3. **The influence of neutrino transport**: The paper explores the influence of different neutrino transport techniques (such as the multi - component flux - limited diffusion method MGFLD and the discrete - ordinate method Sn) on the dynamics of the post - bounce stage of CCSN. By comparing the results of MGFLD and Sn, the researchers found that the Sn method can more accurately capture the asymmetry of the radiation field, thus affecting the neutrino heating efficiency and the evolution of the shock wave. 4. **Rotation effects**: The paper also studies the influence of rapid rotation on the dynamics of CCSN. By comparing the results of non - rotating models and rapid rotating models, it is found that rapid rotation will lead to enhanced neutrino heating in the polar regions and reduced convective overturning in the equatorial region, which further affects the shape and stability of the shock wave. 5. **General - relativistic effects**: The paper introduces a 3D fully general - relativistic method based on the Cactus computational toolkit, which can not only more accurately simulate the hydrodynamics of CCSN, but also dynamically form black holes, which cannot be achieved by Newtonian or pseudo - general - relativistic methods. In summary, this paper aims to deeply understand the explosion mechanism of CCSN through multi - dimensional, multi - scale, and multi - physical computational models, especially the influence in terms of neutrino transport, rotation effects, and general - relativistic effects.