Non-intrusive hierarchical coupling strategies for multi-scale simulations in gravitational dynamics
Non-intrusive hierarchical coupling strategies for multi-scale simulations in gravitational dynamics
复制标题
重力动力学多尺度模拟的非侵入式分层耦合策略
DOI：
10.1016/j.cnsns.2020.105240
复制
发表时间：
2020
影响因子：
3.9
通讯作者：
McMillan, Steve
中科院分区：
文献类型：
--
作者：
Portegies Zwart, Simon;Pelupessy, Inti;Martínez-Barbosa, Carmen;van Elteren, Arjen;McMillan, Steve
分层的代码耦合策略使得将各个数值求解器的结果组合成一个自洽的辛解成为可能。我们探讨了允许这样的耦合策略是非侵入性的可能性。在这种情况下，底层的数值实现不受耦合本身的影响，但其功能在接口中延续。这种方法对于求解大尺度范围内自引力系统的运动方程是有效的。我们采用专门的积分器来解决问题的每个特定部分，并将结果组合成一个自洽的解决方案。特别是，我们探索了将嵌入到宏观系统中的一个或多个微观系统的演化结合在一起的可能性。这里给出的桥的概括包括更高阶的耦合策略(从经典的2阶到10阶)，但我们也演示了如何嵌套多个桥，以及如何在桥的时间步长引入额外的过程来丰富物理，例如通过加入耗散过程或。这样的增强允许在经典的NewtonianN-Body积分器中包括额外的过程，而不需要改变底层代码。例如，这些附加过程包括雅可夫斯基效应、动力摩擦或相对论动力学。所提出的方法是非侵入式的，因为底层方法在不改变代码的情况下保持可操作(除了添加Get和Set函数以启用桥运算符之外)。因此，基本积分器继续以其内部时间步长运行，并保持其局部优化和并行性。多个桥梁可以嵌套并分层耦合，从而允许构建多个嵌套的加强型桥梁的复杂环境。虽然耦合拓扑可能会变得相当复杂，但我们引入了层次耦合语言(HCL)，这是一种可以描述复杂桥接拓扑的元语言。元语言旨在促进对更复杂的层次结构的讨论，在这种层次结构中，桥操作符被引入为模式。我们给出了其中几种情况的示例应用程序，并讨论了这些积分器可以应用的条件。当我们应用该方法模拟绕银河系运行的星团中的行星系统(Au空间和年-时间尺度)时，典型的应用范围在时间和空间尺度上超过10个数量级(100 kpc-空间尺度和10个Gyr-时间尺度)。
Hierarchical code coupling strategies make it possible to combine the results of individual numerical solvers into a self-consistent symplectic solution. We explore the possibility of allowing such a coupling strategy to be non-intrusive. In that case, the underlying numerical implementation is not affected by the coupling itself, but its functionality is carried over in the interface. This method is efficient for solving the equations of motion for a self-gravitating system over a wide range of scales. We adopt a dedicated integrator for solving each particular part of the problem and combine the results to a self-consistent solution. In particular, we explore the possibilities of combining the evolution of one or more microscopic systems that are embedded in a macroscopic system. The here presented generalizations ofBridgeinclude higher-order coupling strategies (from the classic 2nd order up to 10th-order), but we also demonstrate how multiple bridges can be nested and how additional processes can be introduced at the bridge time-step to enrich the physics, for example by incorporating dissipative processesor. Such augmentation allows for including additional processes in a classic NewtonianN-body integrator without alterations to the underlying code. These additional processes include for example the Yarkovsky effect, dynamical friction or relativistic dynamics. Some of these processes operate on all particles whereas others apply only to a subset.The presented method is non-intrusive in the sense that the underlying methods remain operational without changes to the code (apart from adding the get- and set-functions to enable the bridge operator). As a result, the fundamental integrators continue to operate with their internal time step and preserve their local optimizations and parallelism. Multiple bridges can be nested and coupled hierarchically, allowing for the construction of a complex environment of multiple nested augmented bridges. While the coupling topology may become rather complicated, we introduce the hierarchical coupling language (HCL), a meta language in which complex bridge topologies can be described. The meta language is meant for stimulating the discussion on even more complex hierarchies in which the bridge operators are introduced as patternsWe present example applications for several of these cases and discuss the conditions under which these integrators can be applied. Typical applications range over 10 orders of magnitude in temporal and spatial scales when we apply the method to simulating planetary systems (au spatial and year-temporal scale) in a star cluster that orbits in the Galaxy (100 kpc-spatial and 10 Gyr-temporal scale).
登录
查看更多内容
DOI：
10.1111/j.1365-2966.2011.20137.x
发表时间：
2011
影响因子：
4.8
作者：
F. I. Pelupessy;S. Zwart
通讯作者：
S. Zwart
使用 18600 个 GPU 模拟银河系的引力树代码需要 24.77 Pflops
DOI：
--
发表时间：
2014
期刊：
International Conference for High Performance Computing, Networking, Storage and Analysis
影响因子：
--
作者：
J. Bédorf;E. Gaburov;M. Fujii;Keigo Nitadori;T. Ishiyama;S. Zwart
通讯作者：
S. Zwart
galpy：银河动力学包
DOI：
--
发表时间：
2014
期刊：
影响因子：
--
作者：
J. Bovy
通讯作者：
J. Bovy
基于 Ahmad-Cohen 方案的 Hermite 积分器解决引力多体问题
DOI：
--
发表时间：
1992
期刊：
影响因子：
--
作者：
J. Makino;S. Aarseth
通讯作者：
S. Aarseth
影响因子：
2.2
作者：
M. Sofroniou;G. Spaletta
通讯作者：
M. Sofroniou;G. Spaletta