NcorpiON : A O ( N ) software for N-body integration in collisional and fragmenting systems

被引:0
|
作者
Couturier, Jeremy [1 ,2 ]
Quillen, Alice C. [1 ,2 ]
Nakajima, Miki [1 ,2 ]
机构
[1] Univ Rochester, Dept Earth & Environm Sci, 227 Hutchison Hall, Rochester, NY 14627 USA
[2] Univ Rochester, Dept Phys & Astron, 227 Hutchison Hall, Rochester, NY 14627 USA
基金
美国国家科学基金会;
关键词
N-body; Fast multipole method; Mesh; Fragmentation; Collision; FalcON; DISRUPTION; ALGORITHM; EVOLUTION; PARTICLE; PLANET; RING;
D O I
10.1016/j.newast.2024.102313
中图分类号
P1 [天文学];
学科分类号
0704 ;
摘要
NcorpiON is a general purpose N-body software initially developed for the time-efficient integration of collisional and fragmenting systems of planetesimals or moonlets orbiting a central mass. It features a fragmentation model, based on crater scaling and ejecta models, able to realistically simulate a violent impact. The user of NcorpiON can choose between four different built-in modules to compute self-gravity and detect collisions. One of these makes use of a mesh-based algorithm to treat mutual interactions in O ( N ) time. Another module, much more efficient than the standard Barnes-Hut tree code, is a O ( N ) tree-based algorithm called FalcON. It relies on fast multipole expansion for gravity computation and we adapted it to collision detection as well. Computational time is reduced by building the tree structure using a three-dimensional Hilbert curve. For the same precision in mutual gravity computation, NcorpiON is found to be up to 25 times faster than the famous software REBOUND. NcorpiON is written entirely in the C language and only needs a C compiler to run. A python addon, that requires only basic python libraries, produces animations of the simulations from the output files. NcorpiON can communicate with REBOUND's webGL viewer via MPI for 3D visualization. The name NcorpiON, reminding of a scorpion, comes from the French N-corps, meaning N-body, and from the mathematical notation O(N), due to the running time of the software being almost linear in the total number N of bodies. NcorpiON detects collisions and computes mutual gravity faster than REBOUND, and unlike other N-body integrators, it can resolve a collision by fragmentation. The fast multipole expansions are implemented up to order eight to allow for a high precision in mutual gravity computation.
引用
收藏
页数:20
相关论文
共 50 条
  • [42] VIBRATIONS OF RAPIDLY ROTATING N-BODY SYSTEMS
    JELLINEK, J
    LI, DH
    CHEMICAL PHYSICS LETTERS, 1990, 169 (05) : 380 - 386
  • [43] Application of cellular automata to N-body systems
    Lejeune, A
    Perdang, J
    Richert, J
    PHYSICAL REVIEW E, 1999, 60 (03): : 2601 - 2611
  • [44] MEASURING DYNAMICAL AGE OF N-BODY SYSTEMS
    VONHOERNER, S
    ASTRONOMY & ASTROPHYSICS, 1976, 46 (02) : 293 - 302
  • [45] Made-to-measure N-body systems
    Syer, D
    Tremaine, S
    MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, 1996, 282 (01) : 223 - 233
  • [46] N-BODY MODELS
    MILLER, RH
    ANNALS OF THE NEW YORK ACADEMY OF SCIENCES, 1990, 596 : 120 - 122
  • [47] Two-body correlations in N-body boson systems
    Sørensen, O.
    Fedorov, D.V.
    Jensen, A.S.
    Physical Review A - Atomic, Molecular, and Optical Physics, 2002, 66 (03): : 325071 - 325071
  • [48] THE N-BODY PROBLEM
    ADOMIAN, G
    FOUNDATIONS OF PHYSICS LETTERS, 1993, 6 (06) : 597 - 602
  • [49] Collisional evolution of galaxy clusters:: Fokker-Planck and N-body models
    Takahashi, K
    Sensui, T
    Funato, Y
    Makino, K
    ASTROPHYSICAL SUPERCOMPUTING USING PARTICLE SIMULATIONS, 2003, (208): : 449 - 450
  • [50] N-BODY PROBLEM
    VETROVEC, M
    ACTA CIENTIFICA VENEZOLANA, 1968, 19 (01): : 24 - &