Related papers: A parallel-GPU code for asteroid aggregation probl…
Geometrical cross sections of dust aggregates determine their coupling with disk gas, which governs their motions in protoplanetary disks. Collisional outcomes also depend on geometrical cross sections of initial aggregates. In the previous…
I begin with a brief history of N-body simulation and visualization and then go on to describe various methods for creating images and animations of modern simulations in cosmology and galactic dynamics. These techniques are incorporated…
In a previous paper we introduced a new method for simulating collisional gravitational $N$-body systems with linear time scaling on $N$, based on the Multi-Particle Collision (MPC) approach. This allows us to simulate globular clusters…
This paper presents novel approaches to parallelizing particle interactions on a GPU when there are few particles per cell and the interactions are limited by a cutoff distance. The paper surveys classical algorithms and then introduces two…
We developed a new direct-tree hybrid N-body algorithm for fully self-consistent N-body simulations of star clusters in their parent galaxies. In such simulations, star clusters need high accuracy, while galaxies need a fast scheme because…
The gravitational $N$-body problem, which is fundamentally important in astrophysics to predict the motion of $N$ celestial bodies under the mutual gravity of each other, is usually solved numerically because there is no known general…
We present the first detailed comparison between million-body globular cluster simulations computed with a H\'enon-type Monte Carlo code, CMC, and a direct $N$-body code, NBODY6++GPU. Both simulations start from an identical cluster model…
Preliminary results are presented about a fully self-consistent N-body simulation of a sample of four massive globular clusters in close interaction within the central region of a galaxy. The N-body representation (with N=1.5x10^6 particles…
Simulating large-scale articulated assemblies poses a significant challenge due to the numerical stiffness and geometric complexity of jointed structures. Conventional rigid body solvers struggle with the high nonlinearity induced by…
We study the motion of charged particles constrained to arbitrary two-dimensional curved surfaces but interacting in three-dimensional space via the Coulomb potential. To speed-up the interaction calculations, we use the parallel compute…
We present an alternative GPU acceleration for plane waves pseudopotentials electronic structure codes designed for systems that have small unit cells but require a large number of k points to sample the Brillouin zone as happens, for…
We investigate the performance gains from hyper-systolic implementations of n^2-loop problems on the massively parallel computer Quadrics, exploiting its 3-dimensional interprocessor connectivity. For illustration we study the communication…
We present our experience of porting the code used in the wave-packet convergent-close-coupling (WP-CCC) approach to run on NVIDIA V100 and AMD MI250X GPUs. The WP-CCC approach is a method used in the field of ion-atom collision physics to…
The stability of a galaxy model is most easily assessed through N-body simulation. Particle-mesh codes have been widely used for this purpose, since they enable the largest numbers of particles to be employed. We show that the functional…
Preliminary results are presented about a fully self-consistent N-body simulation of a sample of four massive globular clusters in close interaction within the central region of a galaxy. The N-body representation (with N=1.5x10^6 particles…
The kd-tree is a fundamental tool in computer science. Among other applications, the application of kd-tree search (by the tree method) to the fast evaluation of particle interactions and neighbor search is highly important, since the…
There is increasing evidence that many km-sized bodies in the Solar System are piles of rubble bound together by gravity. We present results from a project to map the parameter space of collisions between km-sized spherical rubble piles.…
A way to probe alternative theories of gravitation is to study if they could account for the structures of the universe. We then modified the well-known Gadget-2 code to probe alternative theories of gravitation through galactic dynamics.…
The Tree-Particle-Mesh (TPM) N-body algorithm couples the tree algorithm for directly computing forces on particles in an hierarchical grouping scheme with the extremely efficient mesh based PM structured approach. The combined TPM…
The most difficult aspect of the realistic modeling of granular materials is how to capture the real shape of the particles. Here we present a method to simulate granular materials with complex-shaped particles. The particle shape is…