Celestial Body Simulator

Your browser does not support WebGL.

Start

Settings

$\mathrm{G}=$

About

This is a Newtonian gravitational simulator written using Three.js. The bodies orbit using the following formula: $$m_i\frac{d^2{\mathrm{<b>q</b>}}_i}{d{t}^2}=\sum _{j=0,j\ne i}^{n-1}\frac{\mathrm{G}{m}_i{m}_j({\mathrm{<b>q</b>}}_j-{\mathrm{<b>q</b>}}_i)}{||{\mathrm{<b>q</b>}}_j-{\mathrm{<b>q</b>}}_i{||}^3}$$ where:

• $m_k$ is the mass of body $k$
• ${\mathrm{<b>q</b>}}_k$ is the position of body $k$
• $\mathrm{G}$ is the gravitational constant
• $||{\mathrm{<b>q</b>}}_l-{\mathrm{<b>q</b>}}_k||$ is the magnitude of the Euclidean distance between bodies $k$ and $l$

The simulation often behaves differently between runs. My assumption is that this is due to miniscule differences in timing. An $n$-body system is inherrently chaotic, meaning that small changes to input parameters can cause huge changes in overall behaviour, so these small difference in timing may have a large effect.

You can see the source code hosted on Github, as well as planned features.