Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSimulations of the Early Post-Bounce Phase of Core-Collapse Supernovae in Three-Dimensional Space with Full Boltzmann Neutrino Transport
76
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We report on the core-collapse supernova simulation we conducted for a 11.2 M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashis nuclear equation of state. We focus on the prompt convection at 10 ms after bounce and investigate how neutrinos are transported in the convective matter. We apply a new analysis based on the eigenvalues and eigenvectors of the Eddington tensor and make a comparison between the Boltzmann transport results and the M1 closure approximation in the transition regime between the optically thick and thin limits. We visualize the eigenvalues and eigenvectors using an ellipsoid, in which each principal axis is parallel to one of the eigenvectors and has a length proportional to the corresponding eigenvalue. This approach enables us to understand the difference between the Eddington tensor derived directly from the Boltzmann simulation and the one given by the M1 prescription from a new perspective. We find that the longest principal axis of the ellipsoid is almost always nearly parallel to the energy flux in the M1 closure approximation whereas in the Boltzmann simulation it becomes perpendicular in some transition regions, where the mean free path is 0.1 times the radius. In three spatial dimensions, the convective motions make it difficult to predict where this happens and to possibly improve the closure relation there.
rate research
Read More
We develop a neutrino transfer code for core-collapse simulations, that directly solves the multidimensional Boltzmann equations in full general relativity. We employ the discrete ordinate method, which discretizes the six dimensional phase space. The code is an extension of our special relativistic code coupled to a Newtonian hydrodynamics code, which is currently employed for core-collapse supernova simulations. In order to demonstrate our codes capability to treat general relativistic effects, we conduct some tests: we first compute the free streaming of neutrinos in the Schwarzschild and Kerr spacetimes and compare the results with the geodesic curves; in the Schwarzschild case we deploy not only a 1-dimensional grid in space under spherical symmetry but also a 2-dimensional spatial mesh under axisymmetry in order to assess the capability of the code to compute the spatial advection of neutrinos; secondly, we calculate the neutrino transport in a fixed matter background, which is taken from a core-collapse supernova simulation with our general relativistic but spherically symmetric Boltzmann-hydrodynamics code, to obtain a steady neutrino distribution; the results are compared with those given by the latter code.
One- (1D) and two-dimensional (2D) core-collapse supernova simulations using full Boltzmann neutrino transport for 11.2M and 15.0M progenitor models have been performed to verify the closure relation for the moment method used in the approximate radiation transfer. This study finds areas where the results of the closure relation are inconsistent with those of Boltzmann transport, even for rotational models. In 1D simulations, the Eddington factors p defined in the fluid rest frame (FR) are compared to evaluate the maximum entropy closure for the Fermi-Dirac distribution (MEFD), confirming that MEFD closure performs better than other closures if p < 1/3 and phase space occupancy e > 0.5. In 2D simulations for non-rotating progenitor models, similar results are obtained from the principal-axis analysis of the Eddington tensor kij measured in FR. However, for rotating progenitor models, the principal axes of kij for Boltzmann transport tilt toward oblique directions where matter and neutrinos move relatively fast in azimuthal directions, while the principal axes of kij for MEFD closure are always parallel or perpendicular to the neutrino flux. Thus, the assumption of axisymmetric angular distribution to the flux direction in the closure relation does not hold in the strongly rotating supernova core in the early post-bounce phase. It is also shown that the deviation of the principal axes of kij from the flux direction increases when evaluated in a laboratory frame (LB). The optically thin and thick terms of the pressure tensor in LB negatively impact results in optically thicker and thinner regions, respectively.
We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27-solar-mass star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a 3-species neutrino leakage scheme. The leakage scheme captures the essential aspects of neutrino cooling, heating, and lepton number exchange as predicted by radiation-hydrodynamics simulations. The 27-solar-mass progenitor was studied in 2D by B. Mueller et al. (ApJ 761:72, 2012), who observed strong growth of the SASI while neutrino-driven convection was suppressed. In our 3D simulations, neutrino-driven convection grows from numerical perturbations imposed by our Cartesian grid. It becomes the dominant instability and leads to large-scale non-oscillatory deformations of the shock front. These will result in strongly aspherical explosions without the need for large-scale SASI shock oscillations. Low-l-mode SASI oscillations are present in our models, but saturate at small amplitudes that decrease with increasing neutrino heating and vigor of convection. Our results, in agreement with simpler 3D Newtonian simulations, suggest that once neutrino-driven convection is started, it is likely to become the dominant instability in 3D. Whether it is the primary instability after bounce will ultimately depend on the physical seed perturbations present in the cores of massive stars. The gravitational wave signal, which we extract and analyze for the first time from 3D general-relativistic models, will serve as an observational probe of the postbounce dynamics and, in combination with neutrinos, may allow us to determine the primary hydrodynamic instability.
We present multi-dimensional core-collapse supernova simulations using the Isotropic Diffusion Source Approximation (IDSA) for the neutrino transport and a modified potential for general relativity in two different supernova codes: FLASH and ELEPHANT. Due to the complexity of the core-collapse supernova explosion mechanism, simulations require not only high-performance computers and the exploitation of GPUs, but also sophisticated approximations to capture the essential microphysics. We demonstrate that the IDSA is an elegant and efficient neutrino radiation transfer scheme, which is portable to multiple hydrodynamics codes and fast enough to investigate long-term evolutions in two and three dimensions. Simulations with a 40 solar mass progenitor are presented in both FLASH (1D and 2D) and ELEPHANT (3D) as an extreme test condition. It is found that the black hole formation time is delayed in multiple dimensions and we argue that the strong standing accretion shock instability before black hole formation will lead to strong gravitational waves.
We report on a set of long-term general-relativistic three-dimensional (3D) multi-group (energy-dependent) neutrino-radiation hydrodynamics simulations of core-collapse supernovae. We employ a full 3D two-moment scheme with the local M1 closure, three neutrino species, and 12 energy groups per species. With this, we follow the post-core-bounce evolution of the core of a nonrotating $27$-$M_odot$ progenitor in full unconstrained 3D and in octant symmetry for $gtrsim$$ 380,mathrm{ms}$. We find the development of an asymmetric runaway explosion in our unconstrained simulation. We test the resolution dependence of our results and, in agreement with previous work, find that low resolution artificially aids explosion and leads to an earlier runaway expansion of the shock. At low resolution, the octant and full 3D dynamics are qualitatively very similar, but at high resolution, only the full 3D simulation exhibits the onset of explosion.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
