当前位置: X-MOL 学术Commun. Nonlinear Sci. Numer. Simul. › 论文详情
Our official English website, www.x-mol.net, welcomes your feedback! (Note: you will need to create a separate account there.)
Non-intrusive hierarchical coupling strategies for multi-scale simulations in gravitational dynamics
Communications in Nonlinear Science and Numerical Simulation ( IF 3.9 ) Pub Date : 2020-03-04 , DOI: 10.1016/j.cnsns.2020.105240
Simon Portegies Zwart , Inti Pelupessy , Carmen Martínez-Barbosa , Arjen van Elteren , Steve McMillan

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 of Bridge include 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 Newtonian N-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 patterns

We 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).



中文翻译:

用于重力动力学多尺度模拟的非侵入式分层耦合策略

分层代码耦合策略可以将单个数值求解器的结果组合为一个自洽的辛解。我们探索了允许这种耦合策略非侵入性的可能性。在这种情况下,底层的数值实现不受耦合本身的影响,但其功能会在接口中继续。这种方法对于求解自重系统的运动方程是很有效的。我们采用专门的集成商来解决问题的每个特定部分,并将结果组合为一个自洽的解决方案。特别是,我们探索了将嵌入宏观系统中的一个或多个微观系统的演化结合起来的可能性。这里介绍了桥接器包括更高阶的耦合策略(从经典的二阶到十阶),但是我们还演示了如何嵌套多个桥接器,以及如何在桥接器时步处引入其他过程以丰富物理,例如通过合并耗散过程或 这种扩充允许在经典的牛顿N体积分器中包括其他过程,而无需更改基础代码。这些附加过程包括例如Yarkovsky效应,动态摩擦或相对论动力学。这些过程中的一些对所有粒子起作用,而其他过程仅应用于子集。

在不改变代码的基础上,底层方法保持可操作性的意义上,所提出的方法是非侵入性的(除了添加get和set函数以启用桥运算符之外)。结果,基本集成商继续按照其内部时间步长运行,并保留其本地优化和并行性。多个桥可以嵌套并分层地耦合,从而允许构建多个嵌套的增强桥的复杂环境。尽管耦合拓扑可能变得相当复杂,但我们引入了层次耦合语言(HCL),这是一种可以描述复杂桥拓扑的元语言。元语言旨在激发关于更复杂的层次结构的讨论,在该层次结构中,桥运算符作为模式被引入。

我们提供了其中一些情况的示例应用程序,并讨论了可以应用这些集成器的条件。当我们将方法应用于模拟在银河系中运行的恒星团(100 kpc-时空和10 Gyr-时空)中的行星系统(时空和年-时尺度)时,典型的应用在时空尺度上的变化范围超过10个数量级。规模)。

更新日期:2020-03-04
down
wechat
bug