Computational Relativistic Astrophysics via Space-Time Discontinuous Galerkin Finite Element Methods

计算相对论天体物理学领域已经进入了一个新时代，它对引力波信号的预测与LIGO观测进行了比较，并用于解释观测事件的特征。第一次观测到的引力波只有一年多一点的历史，已经增加了世界范围内对建模可能的引力波源的兴趣。这种模型包括其他类型的辐射，如电磁或中微子对应物，在所谓的“多信使”天体物理学中。在不久的将来，我们将有突破性的发现，这只能与X射线或无线电频谱的发现所带来的科学革命相比较。引力波天文学这一新领域的基础工作正在奠定之中。 引力波是由爱因斯坦方程控制的非常紧凑（密集）的天体物理物体发出的，例如涉及黑洞，中子星，这些双星系统的系统，或者黑洞或中子星形成的坍缩场景。真实地模拟这样的系统不仅需要求解爱因斯坦方程，还需要模拟物质和辐射。这些系统是高度动态的，详细，准确（忠实）的大规模数值计算是理解它们的唯一途径。控制方程太复杂了，无法解析求解，也无法在台式计算机上计算的简单模型中求解。 在当今的高性能计算（HPC）系统上使用计算方法的复杂性和困难严重阻碍了进展。加速器（例如GPU）是常见的，未来的系统预计将需要更多的并行性，同时提供更少的内存带宽，增加科学程序员的负担。随着硬件架构的发展，许多以前高效的算法不再有效，因为它们只利用了新硬件计算能力的一小部分。 在这里提出的工作中，我们将开发新的数值算法来解决这些问题。虽然广义相对论将时空视为一个非常优雅的表述中的单一结构，但当前的主流数值方法并不这样做：它们明确地将时空分为空间和时间，以获得为非相对论场景设计的大量数值方法，但在此过程中也放弃了相对论的许多优雅。 这里开发的数值方法将离散时空，而不是单独的空间和时间，奇怪的是，这样做将有可能使可扩展性和效率提高一个数量级。这将反过来允许模型在物理上更加准确和逼真，因为更多的物理学可以被纳入其描述中。

The field of Computational Relativistic Astrophysics has entered a new era where its predictions of gravitational wave signals are compared to LIGO observations, and are used to interpret signatures of observed events. The first observed gravitational wave is only a little more than one year old, and has already multiplied the world-wide interest in modelling possible gravitational wave sources. This modelling includes other kinds of radiation, such as electromagnetic or neutrino counterparts, in so-called "multi-messenger" astrophysics. The very near future will bring us breakthrough discoveries, to be compared only to the scientific revolution brought about by the discovery of X-rays or the radio spectrum. The groundwork in the new field of gravitational wave astronomy is being laid now. Gravitational waves are emitted by very compact (dense) astrophysical objects which are governed by the Einstein equations, such as systems involving black holes, neutron stars, binary systems of these, or collapse scenarios where black holes or neutron stars are formed. Faithfully modelling such systems requires not only solving the Einstein equations, but also modelling matter and radiation. These systems are highly dynamic, and detailed, accurate (faithful) large-scale numerical calculations are the only road towards understanding them. The governing equations are far too complex to be solved analytically, or in simple models that can be calculated on a desktop computer. Progress is severely hindered by the complexity and difficulties in using computational methods on today's high-performance computing (HPC) systems. Accelerators (e.g. GPUs) are commonplace, and future systems are expected to require even more parallelism while providing less memory bandwidth, increasing the burden of scientific programmers. As hardware architectures evolve, many formerly highly efficient algorithms are not efficient any more, as they only make use of a small fraction of the computing power of newer hardware. In the work proposed here, we will develop novel numerical algorithms to address these issues. While general relativity treats spacetime as a single construct in a very elegant formulation, current mainstream numerical methods do not: They explicitly split spacetime into space and time to gain access to large body of numerical methods designed for non-relativistic scenarios, but also foregoing much of the elegance of relativity in the process. The numerical methods developed here will discretize spacetime, not space and time separately, and by doing so, will curiously have the potential to be an order of magnitude more scalable and efficient. This will in turn allow models that are significantly more physically accurate and realistic, as more physics can be incorporated into their description.