Astrophysics

Orbital Mechanics

N-BODY ENGINE
POST-PROCESSING: HIGH
WebGPU Context Active

Parameters

Star Mass (M)1250
Inject New Body
Orbit Radius (r)15
Inclination (θ)0°
Planet Mass (m)1
Kinetic E.
0.0 J
Potential E.
0.0 J
Energy Distribution
Kinetic E.Total E. 0.0JPotential E.
Live Equation
F⃗g=−GMmr2r^\vec{F}_g = -G\frac{Mm}{r^2}\hat{r}
Mstar=1250 M⊕M_{\text{star}} = 1250\,M_\oplus

Awaiting Telemetry

Governing Dynamics

Gravitational N-body simulation computing pairwise forces between all massive bodies:

mid2r⃗idt2=∑j≠iGmimj∣r⃗j−r⃗i∣3(r⃗j−r⃗i)m_i\frac{d^2\vec{r}_i}{dt^2} = \sum_{j \neq i} G\frac{m_i m_j}{|\vec{r}_j - \vec{r}_i|^3}(\vec{r}_j - \vec{r}_i)

Uses a symplectic integrator to conserve energy over long periods.

vcirc=GMrv_{\text{circ}} = \sqrt{\frac{GM}{r}}