On the Neutrino Distributions in Phase Space for the Rotating Core-collapse Supernova Simulated with a Boltzmann-neutrino-radiation-hydrodynamics Code

• Published in: The Astrophysical journal 2019-02, Vol.872 (2), p.181
• Main Authors: Harada, Akira, Nagakura, Hiroki, Iwakami, Wakana, Okawa, Hirotada, Furusawa, Shun, Matsufuru, Hideo, Sumiyoshi, Kohsuke, Yamada, Shoichi
• Format: Article
• Language: English
• Subjects: Angular distribution; Angular momentum; astronomy & astrophysics; ASTRONOMY AND ASTROPHYSICS; Astrophysics; Collapse; Computational fluid dynamics; Computer simulation; Distribution functions; Fluctuations; Fluid flow; Fluid mechanics; Flux; Gravitational collapse; Hydrodynamic equations; Hydrodynamics; Mathematical analysis; methods: numerical; Neutrinos; Radiation; radiative transfer; Rotation; shock waves; Supernova; supernovae: general; Tensors
• ISSN: 0004-637X; 1538-4357; 1538-4357
• DOI: 10.3847/1538-4357/ab0203

With the Boltzmann-radiation-hydrodynamics code, which we have developed to solve numerically the Boltzmann equations for neutrino transfer, the Newtonian hydrodynamics equations, and the Newtonian self-gravity simultaneously and consistently, we simulate the collapse of a rotating core of the progenitor with a zero-age- main-sequence mass of 11.2 M and a shellular rotation of at the center. We pay particular attention in this paper to the neutrino distribution in phase space, which is affected by the rotation. By solving the Boltzmann equations directly, we can assess the rotation-induced distortion of the angular distribution in momentum space, which gives rise to the rotational component of the neutrino flux. We compare the Eddington tensors calculated both from the raw data and from the M1-closure approximation. We demonstrate that the Eddington tensor is determined by complicated interplays of the fluid velocity and the neutrino interactions and that the M1-closure, which assumes that the Eddington factor is determined by the flux factor, fails to fully capture this aspect, especially in the vicinity of the shock. We find that the error in the Eddington factor reaches ∼20% in our simulation. This is due not to the resolution but to the different dependence of the Eddington and flux factors on the angular profile of the neutrino distribution function, and hence modification to the closure relation is needed.