[论文解读] A Multiple-Grid-Patch Evolution Scheme for 3-D Black Hole Excision

Jonathan Thornburg|ArXiv.org|Dec 4, 2000

Pulsars and Gravitational Waves Research被引用 4

一句话总结

本文提出一种基于 r×(角度坐标) 拓扑结构的多网格区域拼接方案，用于在三维黑洞消去法中替代标准的 xyz 网格，实现平滑、非阶梯状的消去边界，从而减少数值不稳定性。采用四阶有限差分法与 ADM 方程，代码将 Kerr 初始数据演化至约 60M，约束误差在区域拼接角落达到约 1%，目前正实施 BSSN 形式化以提升稳定性。

ABSTRACT

When using black hole excision to numerically evolve a fully generic black hole spacetime, most 3-D 3+1 codes use an $xyz$-topology (spatial) grid. In such a grid, an $r = \constant$ excision surface must be approximated by an irregular and non-smooth "staircase-shaped" excision grid boundary, which may introduce numerical instabilities into the evolution. In this paper I describe an alternate scheme, which uses multiple grid patches, each with topology $\{r imes ({ m angular coordinates})\}$, to cover the slice outside the $r = \constant$ excision surface. The excision grid boundary is now smooth, so the evolution should be less prone to instabilities. With 4th order finite differencing, this code evolves Kerr initial data to ${\sim} 60M$ using the ADM equations; I'm currently implementing the BSSN equations in it in the hope that this will improve the stability.

研究动机与目标

• 解决标准三维 $3+1$ 有限差分代码中因不规则、阶梯状消去边界导致的数值不稳定性问题。
• 开发一种适用于通用黑洞时空的稳定、高阶精度演化的网格系统，采用多区域拼接网格结构。
• 用 r×(角度坐标) 区域替代标准的 xyz 网格拓扑结构，以实现平滑的消去表面。
• 通过仅使用内部数据在重叠区域边界进行一维插值，提升数值稳定性。
• 在新方案中实现并测试 BSSN 形式化，以增强长期演化的稳定性，超越当前的 ADM 形式化。

提出的方法

• 该方案使用多个网格区域，每个区域具有 $r \times (\text{角度坐标})$ 拓扑结构，用于覆盖 $r = \text{常数}$ 消去面之外的空间切片。
• 每个区域周围设有角度方向的虚拟边界区域，这些区域中的场变量通过相邻区域的一维张量插值得到。
• 相邻区域共享 $r$ 和一个角度坐标，确保插值与网格连续性的一致性。
• 网格在径向方向上呈平滑非均匀分布，空间方向采用四阶有限差分法。
• 时间积分采用四阶龙格-库塔方法的线法，与先前球对称实现一致。
• 通过使区域重叠多个网格间距，确保插值使用来自区域内部的数据，从而增强稳定性。

实验结果

研究问题

• RQ1 与标准 $xyz$-拓扑结构网格相比，采用 $r \times (\text{角度坐标})$ 拓扑结构的多区域网格系统是否能减少三维黑洞消去法中的数值不稳定性？
• RQ2 在重叠区域边界上使用一维插值对有限差分演化的稳定性与精度有何影响？
• RQ3 ADM 形式化在新方案中的表现如何，特别是在多个区域交汇处？
• RQ4 与 ADM 形式化相比，BSSN 形式化是否能在该多区域框架中进一步提升稳定性？
• RQ5 平滑的消去边界在多大程度上可消除与阶梯状边界相关的不稳定性？

主要发现

• 多网格区域拼接方案成功地将标准 $xyz$ 网格代码中不规则、阶梯状的消去边界替换为平滑的 $r = \text{常数}$ 边界。
• 采用四阶有限差分法与 ADM 方程，代码可将 Kerr 初始数据演化至约 $60M$ 的坐标时间。
• 在三个区域交汇的角落处出现有限差分法不稳定性，表现为约束违反在 $t \sim 60M$ 时达到约 1%。
• 区域交汇处的不稳定性表明当前 ADM 形式化在该多区域设置下存在局限性。
• 作者正在积极地在代码中实现 BSSN 形式化，以研究其在长期数值稳定性方面的改进潜力。
• 通过使用仅包含内部数据的重叠虚拟区域，确保插值过程保持稳定，避免边界伪影。