[论文解读] A Way to 3D Numerical Relativity --- Coalescing Binary Neutron Stars ---

Takashi Nakamura, Ken-ichi Oohara|ArXiv.org|Dec 16, 1998

Pulsars and Gravitational Waves Research被引用 4

一句话总结

本文提出了一种用于合并双中子星的3D数值相对论的基础方法，采用爱因斯坦方程的(3+1)形式，结合共形时间切片和伪最小形变规范条件。该研究展示了初步的测试模拟，验证了该框架在使用向量并行超级计算机资源应对完全相对论性、非轴对称动力学的计算需求下，对中子星并合最后几毫秒的建模可行性，这是引力波的重要来源。

ABSTRACT

One of the most promising sources of gravitational radiation is coalescence of binary neutron stars or black holes. In order to study gravitational radiation at the merging phase of coalescing binary neutron stars which is called the last three-milliseconds, full general relativistic simulations are required. Since coalescence of binary stars is a completely non-axisymmetric and 3 dimensional (in space) event, which needs great powers of computers, we have to develop a method of 3D numerical relativity using a vector-parallel supercomputer. As a first step of this final goal, we started simulations using post-Newtonian hydrodynamics including radiation reaction of gravitational waves from 1989. They gave a lot of perspectives on coalescing events and gravitational radiation from them. Next we started to attack 3D, fully relativistic simulations. First, basic equations on the (3+1)-formalism of the Einstein equations are shown. As for gauge conditions, we use conformal time slicing and psuedo minimal distortion condtions at present. For these conditions as well as to solve initial value equations, we should solve some elliptic partial differential equations. It consumes the greatest part of CPU time. We show recent results of test simulations on coalescing binary neutron stars.

研究动机与目标

• 开发一种3D数值相对论框架，能够模拟双中子星并合的最后几毫秒，这是引力波探测的关键阶段。
• 解决非轴对称、完全相对论性动力学带来的计算挑战，需要高性能向量并行超级计算机。
• 通过首先测试关键组件（如规范条件和椭圆方程求解器）来为未来完整的广义相对论模拟奠定基础。
• 实现对并合阶段引力辐射的精确建模，这是地面引力波探测器的关键目标。
• 在进行全尺度动力学演化之前，通过测试模拟验证该方法的可行性。

提出的方法

• 采用爱因斯坦方程的(3+1)形式作为3D数值相对论的核心数学框架。
• 实施共形时间切片和伪最小形变规范条件，以控制坐标奇点并保持数值稳定性。
• 求解椭圆偏微分方程以获得初始数据和规范条件，这是模拟中最耗计算资源的部分。
• 使用后牛顿流体动力学结合引力波辐射反作用作为初步步骤，以深入了解并合动力学。
• 设计模拟以可移植至向量并行超级计算机，以应对3D相对论流体动力学的高计算负载。
• 进行测试模拟，以在全动力学演化前验证数值方法和规范条件。

实验结果

研究问题

• RQ1 如何有效实现广义相对论的(3+1)形式，以用于双中子星并合的3D模拟？
• RQ2 在高度动态的并合阶段，何种规范条件最适于保持数值稳定性？
• RQ3 在3D中，控制初始数据和规范条件的椭圆方程求解过程在计算上有多吃力？
• RQ4 使用后牛顿流体动力学和辐射反作用的测试模拟能否提供对最终并合动力学的可靠见解？
• RQ5 在实现中子星并合的全3D数值相对论模拟中，主要计算瓶颈是什么？

主要发现

• 采用共形时间切片和伪最小形变规范条件的(3+1)形式被证明是3D数值相对论模拟双中子星的可行框架。
• 求解控制初始数据和规范条件的椭圆偏微分方程消耗了大部分CPU时间，凸显了主要的计算瓶颈。
• 测试模拟成功证明了所提出的数值方法在模拟非轴对称、完全相对论性动力学方面的可行性。
• 使用后牛顿流体动力学结合引力波辐射反作用的方法，为并合动力学和引力波辐射提供了宝贵的初步见解。
• 该框架具有可扩展性，并设计用于利用向量并行超级计算机架构，以支持未来全动力学模拟。
• 该方法为未来模拟双中子星并合最后几毫秒的完整广义相对论模拟奠定了坚实基础。