[论文解读] A parallel implementation of a new fast algorithm for N-body simulations

P. Londrillo, Carlo Nipoti+1 more|arXiv (Cornell University)|Dec 5, 2002

Scientific Research and Discoveries被引用 3

一句话总结

本文提出了一种并行化的 Fortran-90 实现（FVFPS），用于 Dehnen 的 O(N) 动量守恒快速泊松求解器，适用于 N体模拟，实现了有效的 O(N) 缩放和高精度。星系并合模拟表明，等质量并合保持了基本平面，而吸积情景则不然，这是由于结构非同构性所致，尽管两种情景均未能重现 Faber-Jackson 关系或 Kormendy 关系。

ABSTRACT

A new, momentum preserving fast Poisson solver for N-body systems sharing effective O(N) computational complexity, has been recently developed by Dehnen (2000, 2002). We have implemented the proposed algorithms in a Fortran-90 code, and parallelized it by a domain decomposition using the MPI routines. The code has been applied to intensive numerical investigations of galaxy mergers, in particular focusing on the possible origin of some of the observed scaling relations of elliptical galaxies. We found that the Fundamental Plane is preserved by an equal mass merging hierarchy, while it is not in a scenario where galaxies grow by accretion of smaller stellar systems. In addition, both the Faber-Jackson and Kormendy relations are not reproduced by our simulations.

研究动机与目标

• 将 Dehnen 的 O(N) 快速泊松求解器实现为可移植的、并行化的 Fortran-90 代码，用于天体物理 N体模拟。
• 在保持高精度和动量守恒的前提下，实现有效的 O(N) 计算复杂度。
• 通过无耗散并合层次结构研究椭圆星系的结构演化。
• 检验模拟的并合残骸是否能重现如基本平面、Faber-Jackson 和 Kormendy 关系等观测到的标度关系。

提出的方法

• 采用 Dehnen 的算法，使用固定阶数 p=3 的笛卡尔多极展开以实现最佳性能。
• 实现基于质量的开口角准则 θ(M)，以提高精度和效率。
• 使用 MPI 进行域分解，实现跨多个处理器的并行化。
• 将粒子组织为八叉树结构，以实现分层单元分组和相互作用近似。
• 将该方法应用于模拟具有不同并合历史的无耗散星系并合。
• 通过与 GADGET 代码对比验证结果，确保数值可靠性。

实验结果

研究问题

• RQ1 基于等质量碰撞的并合层次结构是否能保持椭圆星系观测到的基本平面？
• RQ2 与基本平面相比，由吸积驱动的并合序列在结构同构性方面如何？
• RQ3 模拟的并合残骸在多大程度上能重现 Faber-Jackson 和 Kormendy 标度关系？
• RQ4 结构非同构性在吸积残骸偏离基本平面的过程中起什么作用？
• RQ5 具有动量守恒的 O(N) 算法能否产生可靠且高分辨率的星系形成模拟？

主要发现

• FVFPS 代码在粒子数 N ≥ 10^4 时实现了操作次数的有效 O(N) 缩放，证实了理论预测。
• 并行版本在 N ≈ 10^6 粒子时，16 个处理器内的强缩放性能良好。
• 等质量并合模拟保持了基本平面，而基于吸积的层次结构则未能保持。
• 吸积情景中的偏差源于结构非同构性，这通过 Sersic 拟合得到证实，即 m 随质量增加而减小。
• 无论是等质量并合还是吸积层次结构，均未能重现 Faber-Jackson 或 Kormendy 关系，尽管前者由于速度弥散和有效半径的补偿效应而保持了基本平面。
• FVFPS 代码的结果与 GADGET 代码结果高度一致，验证了其数值精度和可靠性。