[论文解读] Neutrino transport in type II supernovae: Boltzmann solver vs. Monte Carlo method

Shoichi Yamada, Han-Thomas Janka+1 more|arXiv (Cornell University)|Sep 2, 1998

Neutrino Physics Research被引用 7

一句话总结

本文比较了在静态爆发后核心模型中，针对II型超新星中中微子输运的有限差分玻尔兹曼求解器与蒙特卡洛方法。结果总体上显示两种方法具有良好的一致性，但玻尔兹曼求解器由于角向分辨率有限，在高度各向异性的区域低估了通量和Eddington因子，可能导致在热泡区域对中微子加热的高估——这凸显了采用自适应角向网格以提高精度的必要性。

ABSTRACT

We have coded a Boltzmann solver based on a finite difference scheme (S_N method) aiming at calculations of neutrino transport in type II supernovae. Close comparison between the Boltzmann solver and a Monte Carlo transport code has been made for realistic atmospheres of post bounce core models under the assumption of a static background. We have also investigated in detail the dependence of the results on the numbers of radial, angular, and energy grid points and the way to discretize the spatial advection term which is used in the Boltzmann solver. A general relativistic calculation has been done for one of the models. We find overall good agreement between the two methods. However, because of a relatively small number of angular grid points (which is inevitable due to limitations of the computation time) the Boltzmann solver tends to underestimate the flux factor and the Eddington factor outside the (mean) ``neutrinosphere'' where the angular distribution of the neutrinos becomes highly anisotropic. This fact suggests that one has to be cautious in applying the Boltzmann solver to a calculation of the neutrino heating in the hot-bubble region because it might tend to overestimate the local energy deposition rate. A comparison shows that this trend is opposite to the results obtained with a multi-group flux-limited diffusion approximation of neutrino transport. The accuracy of the Boltzmann solver can be considerably improved by using a variable angular mesh to increase the angular resolution in the semi-transparent regime.

研究动机与目标

• 开发并验证用于II型超新星中微子输运的有限差分玻尔兹曼求解器。
• 在真实的超新星条件下，将玻尔兹曼求解器的精度与蒙特卡洛输运代码进行对比。
• 评估网格分辨率（径向、角向、能量）以及空间对流离散化方案对玻尔兹曼求解器性能的影响。
• 评估求解器局限性对热泡区域中微子加热的影响。
• 探讨通过在半透明区域采用可变角向网格化以提升性能的潜力。

提出的方法

• 采用有限差分S_N方法求解相空间（位置、能量、角度）中的中微子输运方程。
• 玻尔兹曼求解器在角向和能量方向采用固定和可变的网格间距，包括在半透明区域采用自适应角向分辨率。
• 空间对流采用不同的离散化方案，其对解精度的影响已进行敏感性测试。
• 求解器应用于来自超新星模拟的静态爆发后核心模型，假设背景场固定。
• 在一种模型中引入广义相对论修正，以评估其对输运结果的影响。
• 在相同的物理条件下，将结果与高精度蒙特卡洛输运代码进行基准测试。

实验结果

研究问题

• RQ1 在真实的超新星大气中，玻尔兹曼求解器的结果与蒙特卡洛方法相比在定量上如何？
• RQ2 角向网格点数量在多大程度上影响玻尔兹曼求解器的精度，特别是在中微子各向异性较强的区域？
• RQ3 空间对流离散化方案的选择如何影响玻尔兹曼求解器的解稳定性与精度？
• RQ4 采用可变角向网格在多大程度上能改善半透明区域的分辨率？
• RQ5 由于通量和Eddington因子的不准确，玻尔兹曼求解器是否在热泡区域高估了局部能量沉积率？

主要发现

• 玻尔兹曼求解器在超新星模型的大多数区域与蒙特卡洛方法总体上表现出良好的一致性。
• 由于角向分辨率不足，求解器在半透明区和外区低估了通量因子和Eddington因子。
• 这种低估意味着在各向异性较强的热泡区域，可能导致局部中微子加热率被高估。
• 该误差趋势与多群通量限制扩散近似相反，后者倾向于高估加热。
• 采用可变角向网格显著提高了半透明区域的角向分辨率，降低了通量和Eddington因子的误差。
• 本研究证实，角向分辨率是玻尔兹曼求解器在超新星中微子输运中面临的关键限制。