复制(N/R)+

SENR/NRPy+:奇异曲线坐标系中的数值相对论。我们报告了一个新的开源、用户友好的数值相对论代码包SENR/NRPy+。我们的代码将BSSN参考度量公式的先前实现扩展到更广泛的一类曲线坐标系中,使其非常适合于对具有近似或精确对称性的物理配置进行建模。在模拟黑洞动力学的背景下,它比其他广泛使用的开源数值相对论代码效率高出一个数量级。NRPy+提供了一个基于Python的接口,其中方程以自然张量形式编写,并以任意有限差分顺序输出为高效的C代码,将复杂的张量方程置于科学家的指尖,而不需要昂贵的软件许可证。SENR提供了一个算法框架,它将NRPy+生成的C代码组合成一个有效的数值相关性代码。我们根据另外两个已建立的、最先进的代码进行验证,并取得了极好的一致性。我们第一次证明——在穿孔、喇叭和双黑洞演化的背景下——当空间有限差分导数的阶数增加时,约束违反和引力波形误差几乎指数收敛到零,同时保持球面坐标网格在中等分辨率下固定。这种在视界外的行为是显著的,因为数值误差在视界内不会收敛到零,而且沿极轴的所有点都是坐标奇点。该公式通过以单元为中心的网格和一个简单的基变换来处理这种坐标奇点,该变换将张量分量相对于坐标进行解析正则化。未来的计划包括扩展这一公式,允许动态坐标网格和点的双球面分布来有效捕捉轨道上紧凑的双星动力学。