[论文解读] Core-Collapse Supernovae as Supercomputing Science: a status report toward 6D simulations with exact Boltzmann neutrino transport in full general relativity

Kei Kotake, Kohsuke Sumiyoshi+5 more|arXiv (Cornell University)|May 29, 2012

Gamma-ray bursts and supernovae参考文献 33被引用 5

一句话总结

本文报告了在全般相对论框架下，利用六维中微子输运进行核心坍缩超新星模拟的进展，结合了多维辐射流体动力学与精确的玻尔兹曼中微子输运。研究显示，当前的超级计算机（如K计算机）可实现高分辨率的二维模拟，但包含完整中微子输运的三维模拟因内存和计算需求超过3 PB和每时间步8×10^16次操作，必须依赖百亿亿次（exa-scale）计算能力才能实现。

ABSTRACT

This is a status report on our endeavor to reveal the mechanism of core-collapse supernovae (CCSNe) by large-scale numerical simulations. Multi-dimensionality of the supernova engine, general relativistic magnetohydrodynamics, energy and lepton number transport by neutrinos emitted from the forming neutron star as well as nuclear interactions there, are all believed to play crucial roles in repelling infalling matter and producing energetic explosions. These ingredients are nonlinearly coupled with one another in the dynamics of core-collapse, bounce, and shock expansion. Serious quantitative studies of CCSNe hence make extensive numerical computations mandatory. Since neutrinos are neither in thermal nor in chemical equilibrium in general, their distributions in the phase space should be computed. This is a six dimensional (6D) neutrino transport problem and quite a challenge even for those with an access to the most advanced numerical resources such as the "K computer". To tackle this problem, we have embarked on multi-front efforts. In particular we report in this paper our recent progresses in the treatments of multi-dimensional (multi-D) radiation-hydrodynamics. We are currently proceeding on two different paths to the ultimate goal; in one approach we employ an approximate but highly efficient scheme for neutrino transport and treat 3D hydrodynamics and/or general relativity rigorously; some neutrino-driven explosions will be presented and comparisons will be made between 2D and 3D models quantitatively; in the second approach, on the other hand, exact but so far Newtonian Boltzmann equations are solved in two and three spatial dimensions; we will show some demonstrative test simulations. We will also address the perspectives of exa-scale computations on the next generation supercomputers.

研究动机与目标

• 通过大规模、多维数值模拟，推进对核心坍缩超新星爆发机制的理解。
• 解决在全般相对论下使用精确玻尔兹曼方程模拟六维中微子输运所面临的计算挑战。
• 评估在当前及下一代超级计算机上，长期、高分辨率三维模拟中完整中微子输运的可行性。
• 通过近似与精确中微子输运方案的对比，比较二维与三维模型，评估维度效应对爆发动力学的影响。
• 识别限制三维玻尔兹曼模拟的计算瓶颈——内存、运算量及矩阵求逆，并预测百亿亿次系统的需求。

提出的方法

• 采用多路径方法：一条路径使用近似但高效的中微子输运，结合完整的三维流体动力学与一般相对论；另一条路径则在二维与三维中使用精确的牛顿玻尔兹曼输运。
• 对中微子实施六维相空间离散化，分别对空间坐标、能量和角度方向建立独立的网格，以捕捉非平衡态中微子分布。
• 使用稀疏矩阵求解器处理中微子输运产生的大型线性系统，其中包含耦合不同角度与能量网格点的密集子矩阵。
• 估算未来百亿亿次平台的计算成本，以内存（三维模拟最高达3 PB）和浮点运算量（每时间步最高达8×10^16次）为指标。
• 考虑非平衡效应，如非弹性散射、多普勒频移和引力红移，这些效应分别使矩阵规模扩大Nε和Nε³倍。
• 将一般相对论整合进玻尔兹曼求解器作为最终目标，但因计算需求极端庞大，目前暂予推迟。

实验结果

研究问题

• RQ1 在全般相对论下，使用精确六维玻尔兹曼中微子输运进行三维核心坍缩超新星模拟，其计算需求是什么？
• RQ2 使用精确中微子输运的二维与三维模拟在爆发动力学与能量沉积方面有何差异？
• RQ3 在当前的超级计算机（如K计算机）上，三维中微子输运模拟在空间、能量和角度分辨率方面能达到何种水平？
• RQ4 非弹性散射与相对论效应（多普勒频移、引力红移）如何影响中微子输运中线性系统的规模与复杂度？
• RQ5 为实现长期三维模拟并具备完整中微子输运能力，百亿亿次超级计算机在内存与计算吞吐量方面需要达到何种水平？

主要发现

• K计算机可实现高分辨率的二维玻尔兹曼模拟，包含完整中微子输运，但长期三维模拟因内存与计算限制仍不可行。
• 当前使用精确玻尔兹曼输运的三维模拟，稀疏矩阵需约2TB内存，中微子分布函数需20GB内存，每时间步每中微子种类约有6×10^12次浮点运算。
• 若在K计算机上将三维模拟的空间、能量与角度分辨率提高一倍，矩阵内存需求将增至200TB，分布函数内存增至600GB，运算量将提升至每时间步2×10^15次。
• 为实现长期、高分辨率的三维模拟，需约3 PB内存与每时间步8×10^16次运算能力，表明必须依赖百亿亿次超级计算机。
• 引入非弹性散射与相对论效应会使矩阵规模扩大Nε²倍，运算量扩大Nε³倍，此类模拟在无百亿亿次系统支持下将难以实现。
• 在玻尔兹曼求解器中实现全般相对论仍是未来关键挑战，因其将进一步显著提升计算需求。