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

基于时空不连续伽辽金有限元方法的计算相对论天体物理学

• 批准号: RGPIN-2017-04581
• 负责人: Schnetter, Erik
• 金额: $ 4.37万
• 依托单位: University of Waterloo
• 依托单位国家: 加拿大
• 项目类别: Discovery Grants Program - Individual
• 财政年份: 2022
• 资助国家: 加拿大
• 起止时间: 2022-01-01 至 2023-12-31
• 项目状态: 已结题
• 来源: https://www.nserc-crsng.gc.ca/ase-oro/Details-Detailles_eng.asp?id=762579
• 关键词: Computational; Relativistic; Astrophysics; via; Space

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.

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