[论文解读] A strongly hyperbolic first-order CCZ4 formulation of the Einstein equations and its solution with discontinuous Galerkin schemes

Michael Dumbser, Federico Maria Guercilena+3 more|arXiv (Cornell University)|Jul 31, 2017

Numerical methods for differential equations被引用 8

一句话总结

本文提出了一种强双曲型的一阶爱因斯坦方程形式化，即FO-CCZ4，该形式化源自带有约束阻尼的CCZ4系统，可实现高阶间断伽辽金（DG）模拟的稳定性。该研究首次实现了使用一阶形式化进行三维、高阶DG模拟的移动穿刺双黑洞系统，通过多维真空测试进行了验证，并借助自适应网格细化和亚单元限制器得到增强。

ABSTRACT

We present a strongly hyperbolic first-order formulation of the Einstein equations based on the conformal and covariant Z4 system (CCZ4) with constraint-violation damping, which we refer to as FO-CCZ4. As CCZ4, this formulation combines the advantages of a conformal and traceless formulation, with the suppression of constraint violations given by the damping terms, but being first order in time and space, it is particularly suited for a discontinuous Galerkin (DG) implementation. The strongly hyperbolic first-order formulation has been obtained by making careful use of first and second-order ordering constraints. A proof of strong hyperbolicity is given for a selected choice of gauges via an analytical computation of the entire eigenstructure of the FO-CCZ4 system. The resulting governing partial differential equations system is written in non-conservative form and requires the evolution of 58 unknowns. A key feature of our formulation is that the first-order CCZ4 system decouples into a set of pure ordinary differential equations and a reduced hyperbolic system of partial differential equations that contains only linearly degenerate fields. We implement FO-CCZ4 in a high-order path-conservative arbitrary-high-order-method-using-derivatives (ADER)-DG scheme with adaptive mesh refinement and local time-stepping, supplemented with a third-order ADER-WENO subcell finite-volume limiter in order to deal with singularities arising with black holes. We validate the correctness of the formulation through a series of standard tests in vacuum, performed in one, two and three spatial dimensions, and also present preliminary results on the evolution of binary black-hole systems. To the best of our knowledge, these are the first successful three-dimensional simulations of moving punctures carried out with high-order DG schemes using a first-order formulation of the Einstein equations.

研究动机与目标

• 开发一种适用于高阶间断伽辽金（DG）方法的一阶、强双曲型爱因斯坦方程形式化。
• 结合共形与迹消去形式化的优势，并引入约束违反阻尼，以提升数值稳定性。
• 利用DG方法实现高阶、稳定的黑洞时空模拟，特别是双黑洞系统。
• 通过先进的限制技术和自适应网格细化，解决黑洞演化中奇点带来的数值挑战。

提出的方法

• 通过仔细设计的一阶和二阶排序约束，从CCZ4系统推导出一阶形式化，以确保强双曲性。
• 构建一个包含58个演化变量的非守恒双曲型PDE系统，该系统可分解为纯常微分方程和仅含线性退化场的简化双曲系统。
• 在高阶任意高阶方法使用导数（ADER）-DG格式中实现该系统，采用局部时间推进和自适应网格细化。
• 集成三阶ADER-WENO亚单元有限体积限制器，以处理黑洞模拟中出现的奇点。
• 通过一维、二维和三维标准真空测试验证该形式化。
• 将该格式应用于双黑洞系统的演化，首次实现了使用高阶DG和一阶形式化的三维移动穿刺模拟。

实验结果

研究问题

• RQ1 能否构建一种一阶、强双曲型的爱因斯坦方程形式化，以支持高阶DG离散化？
• RQ2 在标准规范选择下，FO-CCZ4形式化如何确保强双曲性？其特征结构为何？
• RQ3 结合自适应网格细化和局部时间推进的高阶DG格式能否成功演化具有移动穿刺的双黑洞系统？
• RQ4 ADER-WENO亚单元限制器在稳定黑洞奇点附近的模拟方面效果如何？
• RQ5 FO-CCZ4形式化在多维真空测试和双黑洞演化中的性能与精度如何？

主要发现

• 通过解析计算完整特征结构，确认FO-CCZ4形式化在选定规范选择下为强双曲型。
• 该系统可分解为纯常微分方程和仅含线性退化场的简化双曲系统，简化了数值处理。
• 该形式化首次实现了使用高阶间断伽辽金格式对移动穿刺双黑洞系统进行成功的三维模拟。
• 该格式在一维、二维和三维真空测试中均实现了高精度和高稳定性。
• 自适应网格细化、局部时间推进与ADER-WENO亚单元限制器的结合，有效管理了黑洞奇点附近的数值不稳定性。
• 双黑洞演化的初步结果表明，高阶DG方法在使用一阶形式化时，对复杂时空动力学具有可行性。