[论文解读] Constraint Preserving Boundary Conditions for Hyperbolic Formulations of Einstein's Equations

Tarfulea, Nicolae|ArXiv.org|Aug 2, 2005

Black Holes and Theoretical Physics参考文献 40被引用 6

一句话总结

本文为有界区域中爱因斯坦方程的一阶对称双曲形式构造了适定且保持约束的边界条件。通过将约束系统转化为等价的无约束‘扩展系统’，并推导出最大非负边界条件，作者确保了在闵可夫斯基时空线性化下的双曲形式中长期数值稳定性和约束保持性。

ABSTRACT

Einstein's system of equations in the ADM decomposition involves two subsystems of equations: evolution equations and constraint equations. For numerical relativity, one typically solves the constraint equations only on the initial time slice, and then uses the evolution equations to advance the solution in time. Our interest is in the case when the spatial domain is bounded and appropriate boundary conditions are imposed. A key difficulty, which we address in this thesis, is what boundary conditions to place at the artificial boundary that lead to long time stable numerical solutions. We develop an effective technique for finding well-posed constraint preserving boundary conditions for constrained first order symmetric hyperbolic systems. By using this technique, we study the preservation of constraints by some first order symmetric hyperbolic formulations of Einstein's equations derived from the ADM decomposition linearized around Minkowski spacetime with arbitrary lapse and shift perturbations, and the closely related question of their equivalence with the linearized ADM system. Our main result is the finding of well-posed maximal nonnegative constraint preserving boundary conditions for each of the first order symmetric hyperbolic formulations under investigation, for which the unique solution of the corresponding initial boundary value problem provides a solution to the linearized ADM system on polyhedral domains.

研究动机与目标

• 解决在具有人工边界的人工有界空间区域中构造稳定、长时间数值解的爱因斯坦方程挑战。
• 确保边界条件保持系统的约束，防止约束违反，从而避免数值精度下降。
• 建立将约束一阶对称双曲系统系统性地转化为等价无约束系统（扩展系统）的方法。
• 在柯西问题和初值边界问题设定下，证明原始约束系统与其扩展对应物之间的等价性。
• 为线性化爱因斯坦方程的具体双曲形式（EC、AA、A）推导显式、适定、最大非负边界条件。

提出的方法

• 将线性化爱因斯坦方程表述为带有约束的一阶对称双曲系统。
• 应用一种通用技术，通过将约束嵌入演化方程中，构造等价的无约束‘扩展系统’。
• 利用对称双曲系统理论，推导出保持约束的最大非负边界条件。
• 验证在纯柯西问题和初值边界问题设定下，扩展系统与原始约束系统在数学上等价。
• 将该框架应用于特定形式：爱因斯坦-克里斯托费尔（EC）、阿列克谢坚科-阿诺德（AA）和阿诺德（A），为每种形式推导边界条件。
• 分析非齐次边界条件及其在扩展系统框架下与约束的一致性。

实验结果

研究问题

• RQ1 一阶对称双曲形式的爱因斯坦方程在有界区域中，什么样的边界条件能保持约束并确保适定性？
• RQ2 如何系统性地将约束双曲系统转化为等价的无约束系统（扩展系统），同时保持数学等价性？
• RQ3 在初值问题和初值边界问题设定下，扩展系统是否与原始约束系统等价？
• RQ4 能否为EC、AA和A形式构造最大非负边界条件，使其保持约束并确保长期数值稳定性？
• RQ5 规范选择（如谐和规范）在确保扩展系统与原始线性化ADM系统之间的一致性中起什么作用？

主要发现

• 成功推导出线性化爱因斯坦方程的爱因斯坦-克里斯托费尔（EC）、阿列克谢坚科-阿诺德（AA）和阿诺德（A）形式的适定、最大非负、保持约束的边界条件。
• 通过将约束嵌入演化方程中构造的扩展系统，已在柯西问题和初值边界问题设定下被证明与原始约束系统在数学上等价。
• 在所推导边界条件下，扩展系统初值边界问题的唯一解，可导出多面体区域上线性化ADM系统的解。
• 该方法确保约束偏差不会随时间增长，从而实现长期稳定的数值模拟。
• 该框架具有通用性，适用于任何带约束的一阶对称双曲系统，为数值相对论提供了一个稳健工具。
• 非齐次边界条件也得到处理，推导出的条件在扩展系统形式下保持与约束的一致性。