N-Body Gravitational Simulator

The simulation is governed by Newton's Law of Universal Gravitation and his Second Law of Motion. The gravitational force between any two bodies i and j is:

Here, $G$ is the gravitational constant, $m_i$ and $m_j$ are the masses, $r_{ij}$ is the distance between them, and $\hat{r}_{ij}$ is the unit vector pointing from body i to body j. The force is always attractive.

To find the motion of a specific body, we sum the gravitational forces from ALL other (N-1) bodies to get the net force. Then, we use Newton's Second Law to find the acceleration:

This gives us $\vec{a}_j$. The simulator then numerically integrates this acceleration over tiny time steps to update the body's velocity and position, a process repeated for all N bodies every frame. This is the core computational challenge of the N-body problem.

Solar System & Exoplanet Dynamics: Astronomers use N-body simulations to model the long-term stability of our solar system and to predict the orbital configurations of newly discovered exoplanet systems. They can test if a proposed system of planets would survive for billions of years or if one would be ejected.

Galaxy Formation & Evolution: On a colossal scale, cosmologists simulate the gravitational interaction of billions of stars and dark matter particles to understand how galaxies like the Milky Way form, merge, and evolve over cosmic time.

Space Mission Trajectory Planning: When planning a mission like the Voyager probes, engineers use precise N-body calculations (including the Sun, planets, and major moons) to plot gravity-assist flybys, saving immense amounts of fuel by using planetary gravity as a slingshot.

Star Cluster Dynamics: In dense globular clusters, stars can undergo close encounters, gravitational slingshots, and even collisions. N-body simulations help astrophysicists understand phenomena like "core collapse" and the creation of binary star systems through dynamical interactions.

First, you might think "the results are accurate because it's running in real-time," but that is a major misconception. This simulator prioritizes calculation speed and clarity of visualization. For instance, if you set the time step Δt too large, "numerical instability" can occur, where energy is not conserved and orbits continuously expand or conversely fall into the center. In practical applications, careful settings suited to the phenomenon are necessary, such as using a time step smaller than 1/1000th of the planet's orbital period.

Next, do not underestimate the "chaotic" nature where tiny differences in initial conditions can drastically change the results. For example, if you run a simulation with Earth's position shifted by just 1 meter, its orbit could be completely different hundreds of years later. This is not a bug in the simulator but an inherent property of the three-body problem. For experiments requiring reproducibility, it is a golden rule to precisely record and manage the values of initial conditions.

Finally, regarding the handling of "collisions". In this tool, stars merge when they overlap, but real celestial collisions are far more complex. They involve not just simple mergers, but also fragmentation, vaporization, and scattering of debris. Interpreting the simulator's results directly as reality is dangerous. In fields specializing in collision phenomena like "Spaceguard," other simulations using more detailed material models are required.