[论文解读] New numerical framework for the generalized Baumgarte-Shapiro-Shibata-Nakamura formulation: The vacuum case for spherical symmetry

M. A. Alcoforado, R. F. Aranha+2 more|arXiv (Cornell University)|May 19, 2021

Pulsars and Gravitational Waves Research参考文献 62被引用 7

一句话总结

本论文提出RIO，一种基于Galerkin-Collocation方法的新型开源光谱代码，用于在球对称情形下求解广义BSSN形式。通过使用在原点和空间无穷远处天然满足正则性的基函数，该方法在无需显式正则化的情况下实现了约束条件的指数收敛，展示了在真空球对称测试中哈密顿量和动量约束的近机器精度准确性。

ABSTRACT

Here we report a developed high performance and simplified version of the code denominated RIO, which can be easily extended, for the generalized BSSN formulation. We implement a code which is regular at the center of symmetry, without use a special procedure for regularization, as usual. We get exponential convergence for constraints. The numerical algorithm is based on the Galerkin-Collocation method developed successfully for diverse physical scenarios by the Numerical Relativity Group at UERJ. For the sake of clarity in presentation, we consider here the most simple case to display the most salient features of the procedure. Thus, we focus on the definite tests of the new numerical framework. The timing and performance of the code show that we can reach a better accuracy close to the machine precision, for the Hamiltonian and momentum constraints. RIO will be an open source code; currently is under continuous development to consider more general and realistic problems.

研究动机与目标

• 开发一个高性能、可扩展的广义BSSN形式在球对称下的数值框架。
• 通过将正则性嵌入光谱基函数中，消除在原点和极轴处对临时正则化程序的需求。
• 利用Galerkin-Collocation光谱方法，在约束违反量上实现指数收敛。
• 为未来扩展至柱对称和三维问题（包括引力坍缩和Brill波坍缩）奠定基础。
• 在UERJ建立统一、可维护的代码库（RIO），用于数值相对论研究。

提出的方法

• 采用Galerkin-Collocation（GC）光谱方法，结合Gauss-Legendre-Lobatto点的配点法与Galerkin投影。
• 基函数选自L2[0, ∞)希尔伯特空间，以确保在r=0和空间无穷远处的内在正则性。
• 通过这些基函数的谱展开，将偏微分方程组转化为自治动力系统。
• 时间演化通过固定步长的四阶Runge-Kutta积分器应用于谱系数来实现。
• 配点分布被设计以避免奇点，且通过基函数选择自然处理边界条件。
• 该方法在具有已知初始数据和约束的真空球对称情形下得到验证。

实验结果

研究问题

• RQ1 光谱Galerkin-Collocation方法是否能在球对称情形下对广义BSSN形式的哈密顿量和动量约束实现指数收敛？
• RQ2 是否可通过基函数选择在原点和空间无穷远处实现内在正则性，从而消除对特殊正则化技术的需求？
• RQ3 与现有BSSN实现相比，该框架在约束违反量和收敛率方面的性能和准确性如何？
• RQ4 该框架是否可扩展至更复杂的场景，如引力坍缩和柱坐标系下的二维问题？
• RQ5 RIO代码在真空球对称测试中的可实现精度和计算效率如何？

主要发现

• RIO代码在哈密顿量和动量约束中实现了指数收敛，表明具有谱精度。
• 约束违反量被降低至接近机器精度水平，展示了极高的数值精度。
• 通过基函数选择，该方法自然地在原点和空间无穷远处实现正则性，无需使用特殊正则化技术。
• 该代码在使用固定步长四阶Runge-Kutta积分器进行时间演化时表现出稳定性和良好性能。
• 该框架具有可扩展性，正积极开发用于未来在二维和三维问题（包括引力坍缩和Brill波坍缩）中的应用。
• RIO代码为开源，专为持续开发和在数值相对论研究中重复使用而设计。