[论文解读] GR-Athena++: General-relativistic magnetohydrodynamics simulations of neutron star spacetimes

William M. Cook, Boris Daszuta+6 more|arXiv (Cornell University)|Nov 8, 2023

Pulsars and Gravitational Waves Research被引用 5

一句话总结

该论文提出GR-Athena++，一种新型广义相对论磁流体动力学（GRMHD）代码，将Athena++的约束传递方法与爱因斯坦方程的Z4c形式化相结合，用于在八叉树自适应网格细化（AMR）下模拟中子星时空。该代码实现了孤立和双中子星系统的稳定、收敛模拟，约束保持精度极高（磁场散度自由约束违反 ≤10⁻¹⁶），重力波波形准确，并在超过10⁵个CPU核心上实现了超过80%的强缩放效率，使强引力场中磁化中子星动力学的百亿亿次（exascale）模拟成为可能。

ABSTRACT

We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in non-rotating stars with total relative violation of the divergence-free constraint of $10^{-16}$. It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin-Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above $80\%$ in excess of $10^5$ CPU cores and excellent weak scaling is shown up to $\sim 5 imes 10^5$ CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exascale computers.

研究动机与目标

• 开发一种稳健、可扩展的GRMHD代码，用于模拟动力学强引力场中的中子星时空。
• 实现具有真实微观物理和磁场动力学的孤立和双中子星长期、高分辨率模拟。
• 在复杂、非对称时空中实现对物理约束（如磁场无散度）的高保真度和重力波形的高精度。
• 在百亿亿级架构上展示GRMHD大规模模拟的强缩放和弱缩放效率。
• 通过模拟中子星并合过程中的磁场增强、湍流和喷流，支持多信使天体物理学研究。

提出的方法

• 该代码在Athena++基础上扩展了约束传递方法，以在广义相对论中保持磁场的无散度条件。
• 通过与爱因斯坦方程的Z4c形式化耦合，实现时空的动态演化，支持强场、动力学引力的模拟。
• 采用八叉树自适应网格细化（AMR）技术，以在需要时高分辨率解析小尺度特征（如开尔文-赫姆霍兹不稳定性）。
• 使用单元中心的度规重构方法以及高阶插值/限制算子，提升精度与稳定性。
• 采用高阶格式计算重力波形，以确保收敛性和精度。
• 代码针对百亿亿级性能进行了优化，强缩放测试达到10⁵个CPU核心，弱缩放测试达到约5×10⁵个核心。

实验结果

研究问题

• RQ1 GR-Athena++能否在长时间模拟中稳定模拟初始具有轴向磁场的孤立中子星，同时保持磁场无散度约束？
• RQ2 GR-Athena++在动力学时空中对环向磁场分量的增长及其引发的磁不稳定性建模的精度如何？
• RQ3 自适应网格细化在双中子星并合过程中对开尔文-赫姆霍兹不稳定性解析的分辨率提升程度如何？
• RQ4 GR-Athena++在大规模超算系统上的强缩放和弱缩放效率如何？
• RQ5 GR-Athena++能否产生收敛的重力波形，并准确模拟包含微观方程态和中微子输运的并合后动力学？

主要发现

• 在孤立中子星模拟中，GR-Athena++的磁场无散度约束总相对违反量达到10⁻¹⁶量级，表明其具有极佳的约束保持能力。
• 该代码成功捕捉到孤立中子星中环向磁场分量增长至总磁能约10%的情况，与磁流体不稳定性理论预测一致。
• 孤立磁化中子星的长期模拟在长达68.7毫秒内保持稳定，内部能量守恒性能优于以往工作。
• 在双中子星并合模拟中，基于磁场强度的AMR在增加一级细化后实现8倍分辨率增益，与均匀网格细化相比，并合时刻的网格块数量减少约2倍。
• 在超过10⁵个CPU核心上，强缩放效率超过80%；对于双星系统，弱缩放效率在约5×10⁴个核心内超过90%；对于单星系统，弱缩放效率在约5×10⁵个核心内超过90%。
• 该代码展示了重力波形的收敛性，并在双星系统中正确演化了引力辐射反作用力，验证了其物理精度。