数值相对论

恒星坍缩为黑洞并释放出引力波的过程中的一个时刻。图中白色球体是黑洞的表观视界，周围的彩色波表示多极引力波（l=2，m=2）。右下角的图是远处的观测者采集到的波形。

数值相对论（英語：numerical relativity）是广义相对论的一个分支，旨在通过数值方法求解爱因斯坦场方程，以模拟强引力场中的物理过程。相对论天文学中的物理系统，如引力坍缩、中子星、黑洞及引力波等等，以及其他不能利用弱场低速情形中结论[a]近似的现象都可以利用数值相对论模拟。[1]

由于爱因斯坦方程的复杂性与非线性，[b]这一领域的模拟需要特定的数值方法。[2]而计算量巨大的三维问题则需要借助超级计算机解决。一些数学与天体物理学问题，比如密接联星及其引力波的数值模拟，目前可以利用数值相对论求解。[1]

概述

首次被人类直接探测到的引力波，GW150914。图中展示了实验探测结果与利用数值相对论模拟的结果。

数值相对论的要旨在于，利用数值计算方法研究目前不能解析地描述的引力场。这些引力场可动可静，既可能是真空场，也可能包含物质场。[c]静止场问题主要研究场的稳定性。而动态引力场问题则分为两个主要方面，初值问题与演化问题。两种问题各自有不同的求解方法。[3]

数值相对论可用于研究宇宙学模型，比如引力坍缩的临界现象及有黑洞或中子星参与的过程（特别是这两种天体的并合与受到的扰动）。这些情形都涉及到对时空演化的追踪。描述这些情形的爱因斯坦场方程可以用不同形式表述。常用方法包含柯西问题方法、特征线法[4]以及雷奇理论（英语：Regge theory）法等。[5]这些方法都是从某个超曲面上引力场某一时刻的情形开始推导的，也就是说基于初始值按照时间的流向追踪它向邻近超曲面的演化。[6]

与其他数值分析问题一样，数值相对论也需要注意数值解是否稳定并收敛，以及它们适用的初始与边界条件。而度规与坐标条件的存在更使得数值相对论显得尤其复杂。能否求出合适的数值解还会受到爱因斯坦方程表述方式的影响。

由于所涉及的场有所不同，经典场论的许多方法在数值相对论中并不可用。不过，在应对大尺度问题时，数值相对论则与计算流体力学、电动力学以及固体力学等计算科学有很多相似之处。数值相对论还涉及到数值分析、并发计算、偏微分方程以及几何学等数学领域。[7]

发展史

理论基础

阿尔伯特·爱因斯坦在1915年系统地阐释了广义相对论。[8]与狭义相对论中的情形类似，时间与空间在广义相对论中统一为时空，而演化则服从爱因斯坦方程。爱因斯坦方程是一组耦合非线性偏微分方程组。在这个方程组诞生后的近一个世纪的时间里，物理学家仅仅解出了少数几个精确解。并且，其中许多的精确解是在高度对称的条件下，也就是说在方程被很大程度简化时，推导出来的，比如均一且各向同性宇宙的弗里德曼解。[9]

为了在更普遍的条件下研究爱因斯坦方程，数值相对论应运而生。数值求解爱因斯坦方程的必要前提是将四维时空还原为各自独立的三维空间与一维时间，也就是时空的“3＋1解构”。这个过程可以通过几种或简或繁的途径实现。这个问题的首个突破是由理查德·阿诺维特、斯坦利·德塞尔以及查尔斯·W·米斯纳于1950年代后期利用哈密顿力学沿着保罗·狄拉克给出的途径实现的。他们对于爱因斯坦方程的推演结果称作ADM形式。[10]出于技术原因，他们的结果并不适用于数值计算。它的双曲性很弱，因而在实际计算中很难使用。不过数值相对论中许多用到3＋1解构的求解方法还是与ADM形式紧密相关。这种解构会将爱因斯坦方程整理为具有一定初始条件的柯西问题。这种形式已经是计算机数值求解惯常解决的问题。[11]

时空的坐标也不能唯一确定。即使固定了起始超曲面上的坐标，但在穿越到邻近的超曲面时，时空坐标也会因不同的取法而不同。[d]这是数值相对论的一个特征。这种度规上的自由度虽然并不会影响实际上的物理过程，但对过程的描述可能会因其改变。假设从起始的超曲面到邻近的超曲面，时间与空间坐标上分别移动了与。位移描述方式上的不同，也就是说与的取法不同，会导致推导出的方程组形式发生变化。适当的取法会为数值方程的数值求解带来潜在的好处，但坐标和度规许多看来“自然”的取法会导致数值解不稳定，并导致模拟失败。[12]

茨维·皮兰

在ADM形式发表时，计算机的算力并不足以应对这一计算。1964年，哈恩和林德奎斯特首次尝试数值求解爱因斯坦方程。[13]随后在1970年代，斯马[14][15]与埃普利[16]继续这方面的尝试。这些早期的尝试主要针对轴对称（“2＋1维”）时空的情况。茨维·皮兰写出了首个可以追踪辐射引力波的柱对称系统演化的程序。[17]皮兰在此项工作中为数值相对论中许多概念奠基，比如自由演化与约束演化。它们是推导起始数据随时间演化的两种途径。[18][19]对称性的引入降低了内存与运算速度的门槛，使得科学家可以通过当时的超级计算机解决一些问题。[17]

早期结果

理查德·斯塔克与茨维·皮兰在1980年代早期首次对于真实的天体物理学问题做了数值计算。[20]他们首次计算了旋转黑洞在产生阶段辐射出的引力波。在这个结果发表后的近二十年的时间里，数值相对论少有结果发表。这可能还是因为计算机算力不够。1990年代，美国“大挑战（英语：Grand Challenges）”中的双黑洞项目组成功地模拟了两个黑洞迎头相撞的过程。他们利用了问题轴对称的特点简化了求解过程。项目组还在后期处理阶段通过这个数值解计算了事件视界。[21]

而在三维时空中数值求解爱因斯坦方程方面最早的尝试专注于不自转的史瓦西黑洞。这种黑洞静止且具有球对称性。对于验证数值相对论，这一问题十分理想：一是由于科学家非常了解这个问题的解析解，二是由于不自转的黑洞的数值解必然对时间收敛，三是由于这个问题还包含了数值模拟中一种非常复杂的模拟对象——黑洞中心的引力奇点。1995年，安尼诺斯等人发表了这一问题最早的结果之一。[22]在文章中，他们写道：

目前的计算机缺少足够的内存与算力以模拟解析的三维时空。这一点阻碍了三维数值相对论的发展。

理论的发展

在之后的几年里，随着计算机功能更为强大，几组研究人员发展了几种可以提升运算效率的技术。拉撒路组使用了早前通过非线性ADM方程组对黑洞融合模拟的结果，基于单黑洞线性摄动理论，得到了更为稳定的程序。[23]而从黑洞模型入手，又有两种可以绕过奇点问题的技术：切除法与穿刺法。[24]2005年，比勒陀利乌斯综合上述方法，采用了适当的坐标条件在模拟双黑洞方面取得了重大突破。[25]这几种新方法的稳定性在之后的几年里进一步提升，使得人们可以模拟两个黑洞在并合前彼此旋转几十以至几百转的情形。除此之外，计算流体力学中的自适应计算网格细化方法也被引入到数值相对论中来。[26]

拉撒路计划

继“大挑战”后，科学家又在拉撒路计划（1998－2005）中利用天体物理学结果对双黑洞融合做短期数值模拟。不过当时，在超级计算机上模拟这个系统演变的过程中，对于爱因斯坦方程的数值求解无论如何都得不到稳定的结果。研究者前后结合了多种近似方法（比如后牛顿方法和单黑洞扰动）来对全过程做数值模拟。[23]

拉撒路计划标志着双黑洞问题当时最为前沿的进展，给出了许多可以用于天体物理学研究的相当精准结果，比如引力波带走的能量与角动量，[27][28]不同质量的黑洞融合时发出的脉冲[29]以及融合后的黑洞的最终质量、动量和角动量。[30]科学家还在这项计划中计算了融合过程中释出的引力波的具体形式。这项结果对于引力波探测器的研制非常重要。同时，他们还预言黑洞撞击会瞬间产生宇宙中最为强大的能暴。在仅仅一秒内以引力波形式放出的能量足以大过银河系自产生以来其中恒星放出的能量总和。而引力波辐射只占系统总质量减量的很小一部分。[31]

切除法

切除法在1990年代后期提出，[32]是指在模拟演化过程中将奇点周围的事件视界“切除”掉。理论上，由于因果律以及视界本身的性质，[e]这种方法不会影响到对于切除部分之外部分的演化的模拟。也就是说，即使没有解出黑洞内的情况也不会妨碍到外部的精确求解。人们可以通过设置适当的边界条件实现“切除”的过程。[33]

尽管切除法在实际应用中非常成功，但它也存在两个小缺陷。首先，在使用切除法时，人们必须审慎地选用坐标条件。尽管物理效应可能不会传播到视界外，但坐标效应可以。比如，在使用椭圆坐标条件时，黑洞内的网格变化会立即传播到视界外。[34]这意味着，人们需要采用特征速度小于光速的双曲型坐标条件（比如谐和坐标条件）。[35]其次，在黑洞运动时，所切除的部分也必须随黑洞一起运动。[33]

切除法在随后几年里随着校准条件的改进稳定性逐步提高。并且科学家也解决了切除部分随运算区域运动的问题。[36][37][38][39][40][35]借由切除法，双黑洞运动轨迹及融合过程的首个长期稳定的计算结果于2005年发表。[25]

穿刺法

在“穿刺法”中，解分为包含黑洞奇点的解析部分（也就是“刺”）[41]和奇点之外数值构造的部分。这种方法是由布里尔－林德奎斯特算法（针对起始时静止的黑洞）推广而得。[42]而如果再做一般化处理则可得到鲍恩－约克算法（针对起始时旋转且运动的黑洞）。[43]直到2005年，所有采用穿刺法的结果都需要刺在模拟过程中坐标位置不变。不过，彼此密联的黑洞自然会受到引力的作用运动。这意味着坐标系本身或会“伸展”或会“扭曲”。这会导致模拟的部分阶段出现数值不稳定。另外那种避开奇点的方法在模拟物质坍缩形成黑洞的过程中，由于也会采用这种坐标条件，同样也会造成类似效应。[44]

2005年，研究人员首次展示了刺在坐标系中运动的可行性，从而解决了这种方法早前存在的部分问题，使其可以用于精确追踪黑洞的长期演化。[25][45][46]通过选择合适的坐标条件并对奇点周围的物理场做粗略的估计，[f]人们可以得到绕着彼此旋转的双黑洞的数值解，并能精确算出它们释放的引力辐射。[47]

自适应网格细化

自适应网格细化（英语：Adaptive mesh refinement）是一种在数值相对论萌芽前即已发展的数值算法。它是在1980年代由乔普提克在研究标量场坍缩的临界现象时引入的。在那种情形中，场参数处于黑洞形成与空间扩张的临界点。[48][49]由于使用球对称，这项工作针对的是一维情形。之后这种算法被推广到二维情形中。[50]二维网格法也用于研究非均匀宇宙[51][52]和史瓦西黑洞。[53]目前，自适应网格已经成为数值相对论的标准工具，用于研究黑洞与其他稠密系统的融合，以及这种事件产生的引力波。[54][55]

近期进展

对于两个中子星并合生成黑洞过程的可视化模拟（NASA，2014）

目前，数值相对论领域每年都会有数十以至数百篇的相关论文发表，内容涵盖广义相对论数学、引力波以及涉及旋转黑洞的天体物理学问题。所使用的方法在对中子星、黑洞等星体组成的双星系统和多黑洞系统的研究中也有运用。[56][57]从这些工作中，科学家预测两个旋转黑洞并合产生的黑洞的速度能达到4000以至10000 km/s，可以让它穿越所有已知的星系。[58][59]而释出的能量占到总体静质量的8％。黑洞旋转轴所可能发出的剧变则可以解释从射电星系中观测到的相对论性喷流。[60]这一领域另一个重要的研究方向是分类LIGO及VIRGO不能探测到的黑洞并合产生的引力辐射。[61]

新近算法的精确程度足以用来验证引力波的探测结果，比如人类首次直接探测到的引力波GW150914，实验数据与数值模拟数据的误差在4％以内。[62]

脚注

注释

  1. ^ 比如基于爱因斯坦方程精确解的后牛顿解构与摄动理论。
  2. ^ 比如在不同的方程表述以及初始和边界条件下，时间演化的双曲性与柯西问题的适定性。
  3. ^ 在广义相对论中，引力场外的其他场都可称为物质场。
  4. ^ 类似地，在狭义相对论中，时间流向与速度在不同的惯性系中也不同。
  5. ^ 视界内部的相互作用不会影响到外部的物理性质
  6. ^ 由于没有物理效应可从黑洞逃逸，这种估计并不会影响到整体结果。

引注

  1. ^ 1.0 1.1 Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010，Preface. Sec. What is numerical relativity?.
  2. ^ Choptuik M. W. Computational Methods in General Relativity: The Theory 2006.
  3. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010，Sec. 2.6 The constraint and evolution equations and Introduction of Ch. 3 Constructing initial data.
  4. ^ Winicour, J., Characteristic Evolution and Matching, Living Reviews in Relativity, 2012, 15 (2), Bibcode:2012LRR....15....2W, doi:10.12942/lrr-2012-2, （原始内容存档于2015年8月15日） 
  5. ^ Gentle, Adrian P., Regge Calculus: A Unique Tool for Numerical Relativity, General Relativity and Gravitation, 2002, 34 (10): 1701—1718 [2018-03-14], ISSN 0001-7701, doi:10.1023/A:1020128425143, （原始内容存档于2018-03-14） 
  6. ^ Cook, G., Initial Data for Numerical Relativity, Living Reviews in Relativity, 2000, 3 (5), Bibcode:2000LRR.....3....5C, arXiv:gr-qc/0007085 可免费查阅, doi:10.12942/lrr-2000-5, （原始内容存档于2006年7月12日） 
  7. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010，Ch. 6 Numerical methods.
  8. ^ Einstein, Albert. Der Feldgleichungen der Gravitation. Sitzungsberiche der Deutschen Akademie der Wissenschaften zu Berlin, Klasse fur Mathematik, Physik, und Technik. 
  9. ^ Stephani, H.; Kramer, D.; MacCallum, M.; Hoenselaers, C.; Herlt, E., Sec. 1.1 What are exact solutions, and why study them?, Exact Solutions of Einstein's Field Equations, Cambridge Monographs on Mathematical Physics, Cambridge University Press, 2003 [2018-03-14], ISBN 9781139435024, （原始内容存档于2016-05-29） 
  10. ^ Arnowitt, R.; Deser, S.; Misner, C. W., The dynamics of general relativity, L. Witten (编), Gravitation: An Introduction to Current Research, New York: Wiley: 227–265, 1962 
  11. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010，Ch. 2 The 3+1 decompostion of Einstein’s equations.
  12. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010，Ch. 4 Choosing coordinates: the lapse and shift.
  13. ^ Hahn, S. G.; Lindquist, R. W., The two-body problem in geometrodynamics, Annals of Physics, 1964, 29: 304–331, Bibcode:1964AnPhy..29..304H, doi:10.1016/0003-4916(64)90223-4 
  14. ^ Smarr, L. The Structure of General Relativity with a Numerical Example (Ph.D.论文). University of Texas, Austin. 1975. 
  15. ^ Smarr, L., Spacetimes generated by computers: Black holes with gravitational radiation, Annals of the New York Academy of Sciences, 1977, 302: 569, doi:10.1111/j.1749-6632.1977.tb37076.x 
  16. ^ Eppley, K. The numerical evolution of the collision of two black holes (Ph.D. Dissertation论文). Princeton University. 1975. 
  17. ^ 17.0 17.1 Piran, T., Cylindrical general relativistic collapse, Physical Review Letters, 1978, 41 (16): 1085, Bibcode:1978PhRvL..41.1085P, doi:10.1103/PhysRevLett.41.1085 
  18. ^ Alcubierre M. Introduction to 3+1 Numerical Relativity 2008，Sec. 2.6 Free versus constrained evolution.
  19. ^ Bona C., Palenzuela-Luque C., Bona-Casas C. Elements of Numerical Relativity and Relativistic Hydrodynamics 2009，Subsec. 2.3.3 Evolution strategies.
  20. ^ Stark, R. F.; Piran, T., Gravitational-wave emission from rotating gravitational collapse, Physical Review Letters, 1985, 55: 891–894, Bibcode:1985PhRvL..55..891S, doi:10.1103/PhysRevLett.55.891 
  21. ^ Matzner, R. A.; Seidel, H. E.; Shapiro, S. L.; Smarr, L.; Suen, W.-M.; Teukolsky, S. A.; Winicour, J., Geometry of a Black Hole Collision, Science, 1995, 270 (5238): 941–947, Bibcode:1995Sci...270..941M, doi:10.1126/science.270.5238.941 
  22. ^ Anninos, P.; Massó, J.; Seidel, E.; Suen, W.-M.; Towns, J., Three-dimensional numerical relativity: The evolution of black holes, Physical Review D, 1995, 52 (4): 2059–2082, Bibcode:1995PhRvD..52.2059A, arXiv:gr-qc/9503025 可免费查阅, doi:10.1103/PhysRevD.52.2059 
  23. ^ 23.0 23.1 Baker, J.; Campanelli, M.; Lousto, C. O., The Lazarus project: A pragmatic approach to binary black hole evolutions, Physical Review D, 2002, 65 (4): 044001, Bibcode:2002PhRvD..65d4001B, arXiv:gr-qc/0104063 可免费查阅, doi:10.1103/PhysRevD.65.044001 
  24. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010，Ch. 13 Binary black hole evolution.
  25. ^ 25.0 25.1 25.2 Pretorius, F., Evolution of Binary Black-Hole Spacetimes, Physical Review Letters, 2005, 95 (12): 121101, Bibcode:2005PhRvL..95l1101P, arXiv:gr-qc/0507014 可免费查阅, doi:10.1103/PhysRevLett.95.121101 
  26. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010，Subsec. 6.2.5 Mesh refinement.
  27. ^ Baker, J.; Brügmann, B.; Campanelli, M.; Lousto, C. O.; Takahashi, R., Plunge Waveforms from Inspiralling Binary Black Holes, Physical Review Letters, 2001, 87 (12): 121103, Bibcode:2001PhRvL..87l1103B, arXiv:gr-qc/0102037 可免费查阅, doi:10.1103/PhysRevLett.87.121103 
  28. ^ Baker, J.; Campanelli, M.; Lousto, C. O.; Takahashi, R., Modeling gravitational radiation from coalescing binary black holes, Physical Review D, 2002, 65 (12): 124012, Bibcode:2002PhRvD..65l4012B, arXiv:astro-ph/0202469 可免费查阅, doi:10.1103/PhysRevD.65.124012 
  29. ^ Campanelli, M., Understanding the fate of merging supermassive black holes, Classical and Quantum Gravity, 2005, 22: 387, Bibcode:2005CQGra..22S.387C, arXiv:astro-ph/0411744 可免费查阅, doi:10.1088/0264-9381/22/10/034 
  30. ^ Baker, J.; Campanelli, M.; Lousto, C. O.; Takahashi, R., Coalescence remnant of spinning binary black holes, Physical Review D, 2004, 69 (2): 027505, Bibcode:2004PhRvD..69b7505B, arXiv:astro-ph/0305287 可免费查阅, doi:10.1103/PhysRevD.69.027505 
  31. ^ Berti, E.; Cardoso, V.; Gonzalez, J. A.; Sperhake, U.; Hannam, M.; Husa, S.; Brügmann, B., Inspiral, merger, and ringdown of unequal mass black hole binaries: A multipolar analysis, Physical Review D, 2007, 76 (6): 064034 (1—40), Bibcode:2007PhRvD..76f4034B, arXiv:gr-qc/0703053 可免费查阅, doi:10.1103/PhysRevD.76.064034 
  32. ^ Alcubierre, M.; Brügmann, B., Simple excision of a black hole in 3+1 numerical relativity, Physical Review D, 2001, 63 (10): 104006, Bibcode:2001PhRvD..63j4006A, arXiv:gr-qc/0008067 可免费查阅, doi:10.1103/PhysRevD.63.104006 
  33. ^ 33.0 33.1 Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010，Subsec. 13.1.2 Black hole excision.
  34. ^ Bona C., Palenzuela-Luque C., Bona-Casas C. Elements of Numerical Relativity and Relativistic Hydrodynamics 2009，Subsec. 6.1.3 Regular initial data.
  35. ^ 35.0 35.1 Shoemaker, D.; et al., Moving black holes via singularity excision, Classical and Quantum Gravity, 2003, 20 (16): 3729-3743, Bibcode:2003CQGra..20.3729S, arXiv:gr-qc/0301111 可免费查阅, doi:10.1088/0264-9381/20/16/313 
  36. ^ Bona, C.; Masso, J.; Seidel, E.; Stela, J., New formalism for numerical relativity, Physical Review Letters, 1995, 75 (4): 600–603, Bibcode:1995PhRvL..75..600B, arXiv:gr-qc/9412071 可免费查阅, doi:10.1103/PhysRevLett.75.600 
  37. ^ Cook, G. B.; et al., Boosted Three-Dimensional Black-Hole Evolutions with Singularity Excision, Physical Review Letters, 1998, 80: 2512–2516, Bibcode:1998PhRvL..80.2512C, arXiv:gr-qc/9711078 可免费查阅, doi:10.1103/PhysRevLett.80.2512 
  38. ^ Alcubierre, M., Hyperbolic slicings of spacetime: singularity avoidance and gauge shocks, Classical and Quantum Gravity, 2003, (20): 607–623, Bibcode:2003CQGra..20..607A, arXiv:gr-qc/0210050 可免费查阅 
  39. ^ Alcubierre, M.; et al., Gauge conditions for long-term numerical black hole evolutions without excision, Physical Review D, 2003, 67 (8): 084023, Bibcode:2003PhRvD..67h4023A, arXiv:gr-qc/0206072 可免费查阅, doi:10.1103/PhysRevD.67.084023 
  40. ^ Brügmann, B.; Tichy, W.; Jansen, N., Numerical Simulation of Orbiting Black Holes, Physical Review Letters, 2004, 92 (21): 211101, Bibcode:2004PhRvL..92u1101B, arXiv:gr-qc/0312112 可免费查阅, doi:10.1103/PhysRevLett.92.211101 
  41. ^ Brandt, S.; Brügmann, B., A Simple Construction of Initial Data for Multiple Black Holes, Physical Review Letters, 1997, 78 (19): 3606-3609, Bibcode:1997PhRvL..78.3606B, arXiv:gr-qc/9703066 可免费查阅, doi:10.1103/PhysRevLett.78.3606 
  42. ^ Brill, D. R.; Lindquist, R. W., Interaction Energy in Geometrostatics, Physical Review, 1963, 131 (1): 471–476, Bibcode:1963PhRv..131..471B, doi:10.1103/PhysRev.131.471 
  43. ^ Bowen, J. M.; York Jr., J. W., Time-asymmetric initial data for black holes and black-hole collisions, Physical Review D, 1980, 21 (8): 2047–2056, Bibcode:1980PhRvD..21.2047B, doi:10.1103/PhysRevD.21.2047 
  44. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010，Subsec. 13.1.1 Singularity avoiding coordinates.
  45. ^ Baker, J. G.; et al., Gravitational-Wave Extraction from an Inspiraling Configuration of Merging Black Holes, Physical Review Letters, 2006, 96 (11): 111102, Bibcode:2006PhRvL..96k1102B, arXiv:gr-qc/0511103 可免费查阅, doi:10.1103/PhysRevLett.96.111102 
  46. ^ Campanelli, M.; et al., Accurate Evolutions of Orbiting Black-Hole Binaries without Excision, Physical Review Letters, 2006, 96 (11): 111101, Bibcode:2006PhRvL..96k1101C, arXiv:gr-qc/0511048 可免费查阅, doi:10.1103/PhysRevLett.96.111101 
  47. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010，Subsec. 13.1.3 The moving puncture method.
  48. ^ Choptuik, M. W. Experiences with an adaptive mesh refinement algorithm in numerical relativity. Evans, C.; Finn, L.; Hobill, D. (编). Frontiers in numerical relativity. Cambridge: Cambridge University Press. 1989. ISBN 0521366666. 
  49. ^ Choptuik, M. W., Universality and scaling in gravitational collapse of a massless scalar field, Physical Review Letters, 1993, 70 (1): 9–12, Bibcode:1993PhRvL..70....9C, doi:10.1103/PhysRevLett.70.9 
  50. ^ Choptuik, M. W.; et al., Critical collapse of the massless scalar field in axisymmetry, Physical Review D, 2003, 68 (4): 044007, Bibcode:2003PhRvD..68d4007C, arXiv:gr-qc/0305003 可免费查阅, doi:10.1103/PhysRevD.68.044007 
  51. ^ Hern, S. D. Numerical relativity and inhomogeneous cosmologies (Ph.D.论文). Cambridge University. 1999. 
  52. ^ Belanger, Z. B. Adaptive mesh refinement in the T2 symmetric spacetime (Master论文). Oakland University. 2001. 
  53. ^ Schnetter, E.; et al., Evolutions in 3D numerical relativity using fixed mesh refinement, Classical and Quantum Gravity, 2004, 21: 1465–1488, Bibcode:2004CQGra..21.1465S, arXiv:gr-qc/0310042 可免费查阅, doi:10.1088/0264-9381/21/6/014 
  54. ^ Imbiriba, B.; et al., Evolving a puncture black hole with fixed mesh refinement, Physical Review D, 2004, 70 (12): 124025, Bibcode:2004PhRvD..70l4025I, arXiv:gr-qc/0403048 可免费查阅, doi:10.1103/PhysRevD.70.124025 
  55. ^ Fiske, D. R.; et al., Wave zone extraction of gravitational radiation in three-dimensional numerical relativity, Physical Review D, 2005, 71 (10): 104036, Bibcode:2005PhRvD..71j4036F, arXiv:gr-qc/0503100 可免费查阅, doi:10.1103/PhysRevD.71.104036 
  56. ^ Etienne, Zachariah B.; Liu, Yuk Tung; Shapiro, Stuart L.; Baumgarte, Thomas W. Relativistic Simulations of Black Hole-Neutron Star Mergers: Effects of black-hole spin. Phys. Rev. D. 2009, 76 (4). Bibcode:2009PhRvD..79d4024E. arXiv:0812.2245 可免费查阅. doi:10.1103/PhysRevD.79.044024. 
  57. ^ Lousto, C. O.; Zlochower, Y., Foundations of multiple-black-hole evolutions, Physical Review D, 2008, 77 (2): 024034, Bibcode:2008PhRvD..77b4034L, arXiv:0711.1165 可免费查阅, doi:10.1103/PhysRevD.77.024034 
  58. ^ Campanelli M.; et al., Maximum Gravitational Recoil, Physical Review Letters, 2007, 98 (23): 231102, Bibcode:2007PhRvL..98w1102C, arXiv:gr-qc/0702133 可免费查阅, doi:10.1103/PhysRevLett.98.231102 
  59. ^ Healy, J.; et al., Superkicks in Hyperbolic Encounters of Binary Black Holes, Physical Review Letters, 2009, 102 (4): 041101, Bibcode:2009PhRvL.102d1101H, arXiv:0807.3292 可免费查阅, doi:10.1103/PhysRevLett.102.041101 
  60. ^ Campanelli, M.; et al., Spin flips and precession in black-hole-binary mergers, Physical Review D, 2007, 75 (6): 064030, Bibcode:2007PhRvD..75f4030C, arXiv:gr-qc/0612076 可免费查阅, doi:10.1103/PhysRevD.75.064030 
  61. ^ Hinder, I.; et al., Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration, Classical and Quantum Gravity, 2013, 31 (2): 025012, Bibcode:2013CQGra..31b5012H, arXiv:1307.5307 可免费查阅, doi:10.1088/0264-9381/31/2/025012 
  62. ^ The LIGO Scientific Collaboration; The Virgo Collaboration, Tests of general relativity with GW150914 (PDF) ArXiv e-prints, 2016 [2018-03-15], Bibcode:2016arXiv160203841T, arXiv:1602.03841 可免费查阅, （原始内容 (PDF)存档于2016-02-15） 

引用书目

  • Choptuik M. W. Computational Methods in General Relativity: The Theory. Françoise, J. P.; Naber, G. L.; Tsun, T. S. (编). Encyclopedia of mathematical physics. Encyclopedia of Mathematical Physics 1. Elsevier. 2006: 604—610. ISBN 9780125126601. 
  • Baumgarte, T. W.; Shapiro S. L. Numerical Relativity: Solving Einstein's Equations on the Computer. Cambridge University Press. 2010. ISBN 9780521514071. 
  • Alcubierre, M. Introduction to 3+1 Numerical Relativity. Oxford University Press. 2008. ISBN 9780199205677. 
  • Bona, C.; Palenzuela-Luque, C.; Bona-Casas, C. Elements of Numerical Relativity and Relativistic Hydrodynamics: From Einstein's Equations to Astrophysical Simulations. Lecture Notes in Physics, Vol. 783 2. Springer-Verlag. 2009. ISBN 9783642011634. 
  • Gourgoulhon, E. 3+1 Formalism in General Relativity: Bases of Numerical Relativity. Lecture Notes in Physics, Vol. 846. Springer-Verlag. 2012. ISBN 9783642245244. 

延伸阅读

外部链接

基礎概念
现象
方程
進階理論
精确解
近似解与数值模拟
科學家