[论文解读] Variational Principles in General Relativity

Brian D. Baker|ArXiv.org|May 26, 2002

Pulsars and Gravitational Waves Research参考文献 35被引用 4

一句话总结

本文利用3+1形式化和曲面积分，为广义相对论中的静止、轴对称、无旋双中子星构建了一个变分原理，以确保质量-能量的一致性。通过非线性自适应多重网格方法求解，实现了引力波和多极矩的精确预测，并通过后牛顿理论进行了数值验证。

ABSTRACT

This thesis details an effort to generate astrophysically interesting solutions to the two-body problem in General Relativity. The thesis consists of two main parts. The first part presents an analytical variational principle for describing binary neutron stars undergoing irrotational fluid flow. The variational principle is a powerful tool for generating accurate estimates of orbital parameters for neutron stars in quasi-equilibrium orbits. The second part of the thesis details the numerical application of the variational principle by solving the initial value problem for binary black holes in quasi- equilibrium circular orbits. The analysis draws from the novel ``puncture'' method of describing the black holes, and relies on nonlinear adaptive multigrid techniques for generating numerical results. We arrive at two important conclusions. First, the analytical variational principle describing binary neutron stars in irrotational motion provides a road map for future numerical simulations, and also lends credence to previous simulations by other authors. Second, the numerical application and description of binary black holes in quasi-equilibrium circular orbits simplifies the analyses of previous authors, and allows for the imposition of realistic boundary data in simulations with relatively high grid densities. Both the variational principle and its application are used to generate accurate estimates of the orbital parameters.

研究动机与目标

• 为广义相对论中的静止、轴对称、无旋双中子星建立一个变分原理。
• 引入曲面积分和边界条件，以确保正确计算ADM质量与能量守恒。
• 实现非线性自适应多重网格求解器，以数值求解耦合的爱因斯坦方程。
• 将数值结果与后牛顿近似进行对比，验证多极矩的收敛性。
• 为建模致密双星系统提供一个自洽框架，实现引力波和质量矩的精确预测。

提出的方法

• 采用时空的3+1分解，将爱因斯坦方程划分为初值问题和演化问题。
• 应用对称迹零张量和托尔内的广义解，以描述流体流动和度规扰动。
• 实施非线性自适应多重网格算法，求解由变分原理导出的椭圆型方程。
• 通过网格加密与插值、限制算子（如 symrstrct、rstrct）处理多分辨率解。
• 采用Galerkin有限差分法结合罗宾边界条件，稳定偏微分方程的数值解。
• 通过计算网格上的数值积分计算ADM质量与多极矩，以验证全局守恒定律。

实验结果

研究问题

• RQ1 如何在广义相对论中为静止、轴对称、无旋双中子星一致地构建一个变分原理？
• RQ2 曲面积分与边界项在确保正确ADM质量与能量守恒中起什么作用？
• RQ3 非线性自适应多重网格方法如何提升求解致密双星爱因斯坦方程的收敛性与精度？
• RQ4 数值结果在多大程度上与后牛顿理论对 lapse 函数和引力波内容的预测一致？
• RQ5 该方法能否准确计算平衡态配置中的多极矩与质量-能量分布？

主要发现

• 变分原理成功生成了满足正确表面与边界条件的双中子星自洽平衡构型。
• 数值结果表明，lapse 函数与引力波振幅与后牛顿理论预测高度一致，验证了方法的准确性。
• 自适应多重网格求解器收敛高效，残差减少与误差控制在多级网格上均得到验证。
• 多极矩与ADM质量以高精度计算，证实了全局能量与动量守恒。
• 对称迹零张量与曲面积分项的实现确保了与广义相对论约束的一致性。
• 即使在强场区域，代码也能以极少的数值伪影准确计算引力势与度规分量。