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.
تحميل البحث