[论文解读] A Proof of Convergence for Numerical Approximations Generated by the Locally Inertial Godunov Method in General Relativity

Zeke Vogler, Blake Temple|arXiv (Cornell University)|Dec 21, 2011

Differential Equations and Boundary Problems参考文献 9被引用 4

一句话总结

本文证明了在广义相对论中求解爱因斯坦方程的局部惯性Godunov方法的收敛性，表明有界且非振荡的数值近似收敛于弱解。该方法使用动态时间步长和边界条件的Riemann问题解，且在总变差有界性和通量与源项光滑性假设下建立了收敛性。

ABSTRACT

We provide details for the proof of convergence of the locally inertial Godunov method with dynamical time dilation.

研究动机与目标

• 为数值相对论中局部惯性Godunov方法的收敛性提供严格的证明。
• 解决在有限模拟中，针对具有可变时间步长和空间边界条件的Godunov格式，证明收敛性的挑战。
• 将Groah和Temple对Glimm格式的论证扩展至Godunov框架，将其适应于平均化方法而非随机抽样。
• 验证数值收敛性和总变差有界性可推出对爱因斯坦方程弱解的收敛性。
• 为使用局部惯性Godunov方法模拟广义相对论中冲击波相互作用的数值模拟提供理论基础。

提出的方法

• 在标准史瓦西坐标系下，使用具有动态时间膨胀的局部惯性Godunov方法，计算爱因斯坦方程的近似解。
• 在每个空间单元内应用Riemann问题解来演化系统，通量和源项由弱等价的爱因斯坦-欧拉系统导出。
• 采用时间依赖的ODE步长，使用光滑函数G、f和h，度规A被视为Lipschitz连续。
• 施加反映有限模拟域边界的边界条件，包括左边界和右边界数据。
• 在每个时间步对解施加总变差有界性，以控制振荡并确保稳定性。
• 采用修改后的Groah和Temple论证，将其适应于Godunov格式的平均化方法，而非随机抽样。

实验结果

研究问题

• RQ1 在何种条件下，局部惯性Godunov方法收敛于爱因斯坦方程的弱解？
• RQ2 当假设总变差有界性而非推导出该有界性时，如何证明数值近似解的收敛性？
• RQ3 可变时间步长和空间边界在广义相对论中Godunov型格式的收敛性分析中起什么作用？
• RQ4 包含项A′·∇A f(Aij,û,x)如何处理Riemann单元边界处度规的不连续性？
• RQ5 即使未事先证明总变差有界性，数值收敛性和有界变差是否仍可推出收敛于弱解？

主要发现

• 该定理表明，若由局部惯性Godunov方法生成的近似解有界且无振荡，则其收敛于爱因斯坦方程的弱解。
• 证明显示，数值格式中的残差误差为O(Δx)，在给定假设下确认了相容性和收敛性。
• 单元平均值与Riemann问题解之间的差异由函数在单元上的总变差有界，如引理2.2所形式化。
• 通过利用G的光滑性和总变差有界性，对ODE步长误差进行有界处理，从而在整体格式中得到O(Δx)的误差。
• 通过在ODE更新中包含项A′·∇A f(Aij,û,x)，该方法考虑了单元边界处度规的不连续性。
• 该结果特别适用于模拟冲击波相互作用，因为在这些情况下，数值收敛性和有界变差的验证是可行的。