🌌 Solar System Physics Simulation

This simulation serves as a comprehensive demonstration of classical gravitational physics and celestial mechanics, implementing Newton's law of universal gravitation as the core governing principle. Every celestial body in the system exerts an attractive force on every other body, with the force magnitude calculated using the inverse square law formula F = G × m₁ × m₂ / r², where G represents the gravitational constant, m₁ and m₂ are the masses of the interacting bodies, and r is the distance between their centers. The simulation accurately models the complex n-body problem, showing how multiple gravitational interactions create stable orbital configurations, resonant patterns, and dynamic instabilities. By allowing users to manipulate parameters like gravitational strength and time acceleration, the simulation reveals how small changes in fundamental constants can dramatically alter the stability and evolution of entire planetary systems, providing insights into why our solar system maintains its current configuration while others might collapse into chaos or eject planets entirely.

The simulation implements Kepler's three laws of planetary motion with mathematical precision, demonstrating how these empirical laws emerge from Newton's gravitational theory. Kepler's first law is visualized through elliptical orbits where the Sun occupies one focus, with orbital eccentricity determining the degree of elongation from perfect circles. The second law manifests in the varying orbital speeds of planets, where bodies move faster when closer to the Sun and slower when farther away, ensuring that the line connecting a planet to the Sun sweeps equal areas in equal time intervals. Kepler's third law is demonstrated through the relationship between orbital period and semi-major axis, where P² ∝ a³ for objects orbiting the same central mass. The simulation also shows perturbations from Keplerian orbits caused by gravitational interactions between planets, asteroids, and other bodies, illustrating how these small deviations can lead to orbital resonances, chaotic behavior, and long-term orbital evolution that shapes the architecture of planetary systems over astronomical timescales.

The n-body problem represents one of the most challenging computational challenges in physics, involving the prediction of motions for multiple celestial bodies under mutual gravitational attraction. Unlike the two-body problem which has analytical solutions, the n-body problem generally requires numerical integration methods for accurate solutions. This simulation employs the Velocity Verlet integration algorithm, a symplectic integrator that preserves the fundamental conservation laws of energy and momentum over long simulation times. The algorithm calculates gravitational forces between all pairs of bodies in each time step, then updates positions and velocities using finite difference approximations. While exact analytical solutions exist only for special cases like hierarchical systems, the numerical approach provides sufficiently accurate results for most astronomical scenarios. The simulation demonstrates how small perturbations can accumulate over time, potentially leading to chaotic behavior where tiny changes in initial conditions result in dramatically different outcomes, a phenomenon known as sensitive dependence on initial conditions.

Gravitational force vectors are calculated in real-time for every pair of celestial bodies using vector mathematics and Newton's gravitational law. For each body, the net gravitational force is computed as the vector sum of individual forces from all other bodies, with each force component calculated using F = G × m₁ × m₂ / r² × û, where û represents the unit vector pointing from one body to another. Orbital trajectories emerge naturally from the integration of these forces over time, with the simulation using numerical methods to solve the differential equations of motion. The velocity of each body is updated based on the net acceleration caused by gravitational forces, while positions are updated based on current velocities. The simulation includes collision detection algorithms that handle close encounters between bodies, either through elastic collisions that conserve momentum and energy, or inelastic mergers that combine masses and momenta. Orbital trails are rendered by storing historical position data and connecting these points with smooth curves, allowing users to visualize the complex paths that celestial bodies follow under the influence of multiple gravitational sources.

Conservation laws form the mathematical foundation that ensures the physical accuracy and stability of the simulation over extended periods. The law of conservation of momentum is strictly enforced, meaning the total momentum vector of the system remains constant unless external forces are applied, which is crucial for maintaining realistic orbital dynamics. Energy conservation is implemented through the use of symplectic integration methods that minimize numerical dissipation, ensuring that the total mechanical energy of the system (kinetic plus potential) remains nearly constant over time. Angular momentum conservation is particularly important for rotational systems and orbital mechanics, preventing unphysical changes in the rotational state of the system. The simulation demonstrates how these conservation laws constrain the possible motions of celestial bodies, explaining why certain orbital configurations are stable while others lead to catastrophic disruptions. By monitoring these conserved quantities in real-time, users can assess the numerical accuracy of the simulation and detect when approximations might be breaking down due to computational limitations.

The simulation distinguishes between different types of celestial bodies through their physical properties and behavioral characteristics, each following appropriate physical models. Stars are modeled as massive, luminous bodies that dominate the gravitational potential of the system, providing the central anchor around which planets orbit. Planets are represented as large, spherical bodies with significant mass and gravitational influence, capable of clearing their orbital zones of smaller debris. Moons and satellites orbit planets rather than the central star, creating complex hierarchical systems with multiple centers of gravity. Asteroids are small, rocky bodies with minimal gravitational influence but significant numbers, often forming belts or scattered distributions. Comets are modeled with highly elliptical orbits and volatile compositions that affect their long-term stability. Dwarf planets and Kuiper belt objects demonstrate the transition between planetary and minor body characteristics. Each body type has unique visual representations, orbital behaviors, and physical interactions that reflect their astronomical properties and dynamical roles within the solar system architecture.

Orbital resonances occur when celestial bodies have orbital periods that are related by simple integer ratios, leading to regular gravitational perturbations that can stabilize or destabilize orbits. The simulation models these resonances by tracking the mean motion ratios between bodies and applying resonant perturbation terms to their orbital elements. First-order resonances like the 2:1 ratio between Jupiter and Mars create stable configurations that protect certain orbits from chaotic diffusion. Higher-order resonances can lead to complex dynamical behavior, including the formation of Kirkwood gaps in the asteroid belt where resonant perturbations eject bodies from stable orbits. The simulation demonstrates how resonances affect orbital eccentricity and inclination over time, showing how Jupiter's resonances sculpt the architecture of the asteroid belt and protect the inner planets from cometary impacts. Users can observe how changing the mass distribution or orbital parameters affects the establishment and maintenance of resonant configurations, providing insights into the dynamical history of our solar system and the formation of planetary systems around other stars.

Collision detection and gravitational interactions are implemented using sophisticated algorithms that balance computational efficiency with physical accuracy. The simulation uses a multi-stage collision detection system that first employs broad-phase algorithms to identify potentially colliding pairs, followed by narrow-phase precise collision detection using geometric intersection tests. Gravitational interactions are calculated using a Barnes-Hut tree algorithm for large numbers of bodies, approximating distant gravitational effects while computing exact forces for nearby interactions. Collision outcomes depend on the selected interaction mode: elastic collisions conserve both momentum and kinetic energy, inelastic collisions result in mergers that combine masses and redistribute momentum, while fragmentation modes can break bodies into smaller components. The simulation also models tidal interactions, where differential gravitational forces cause tidal heating, orbital circularization, and rotational synchronization. These complex interactions demonstrate how gravitational dynamics shape the evolution of planetary systems, from the formation of planets through accretion to the long-term orbital stability that allows life to develop on planetary surfaces.

The simulation offers valuable insights into the formation and evolution of planetary systems through its demonstration of accretion processes, orbital migration, and dynamical instabilities. By allowing users to add bodies with different masses and orbital parameters, it shows how gravitational interactions can lead to runaway growth where massive bodies accrete smaller ones more rapidly. The simulation illustrates the concept of orbital migration, where planets can spiral inward or outward due to gravitational interactions with protoplanetary disks or other planets. Dynamical instabilities are demonstrated through the sudden ejection of planets or the onset of chaotic orbital behavior when certain mass ratios or orbital configurations are reached. The formation of resonant chains and orbital gaps provides insights into how planetary systems self-organize into stable configurations. By experimenting with different initial conditions, users can explore how the formation history of a planetary system influences its final architecture, explaining why some systems have closely packed inner planets while others have wide orbital gaps or eccentric orbits that challenge theories of gentle planetary formation.

This simulation serves as a powerful educational tool for teaching physics, astronomy, and computational science across multiple levels. For undergraduate physics students, it provides hands-on experience with gravitational physics, numerical methods, and the challenges of computational modeling. Graduate students can use it to explore advanced topics like chaotic dynamics, stellar system stability, and the n-body problem. The simulation's interactive nature allows researchers to test hypotheses about planetary system formation, orbital resonances, and dynamical evolution without the computational complexity of professional astronomical software. Educational applications include demonstrating Kepler's laws, explaining tidal locking, and illustrating the concept of gravitational potential wells. The ability to manipulate parameters in real-time makes abstract mathematical concepts concrete and accessible. For research purposes, the simulation can model specific astronomical scenarios, test the stability of exoplanetary systems, or explore the dynamics of star clusters and galactic structures. Its open architecture allows for extensions to include additional physical effects like general relativity corrections or magnetic field interactions, making it a versatile platform for both teaching and scientific investigation.