Methodology
This simulation models the Lunar Gateway's Near-Rectilinear Halo Orbit (NRHO) using a three-body gravitational framework (Earth, Moon, Gateway) with verified NASA and NIST constants. Below are the hard-coded values used by the physics engine.
NASA / NIST Physics Constants
Independent Variable
Initial orbital insertion velocity, tested at three levels: 0.855 km/s (−5%), 0.900 km/s (baseline), 0.945 km/s (+5%).
Dependent Variable
Perilune altitude deviation (km) from the ideal 3,000 km perilune, measured over 30 simulated days.
Control Variables
- Gateway satellite mass: 40,000 kg (locked)
- Orbital insertion angle: 90° (locked)
- Lunar mass: 7.342 × 10²² kg (constant)
- Earth-Moon distance: 384,400 km (constant)
- Simulation duration: 30 days per trial (constant)
Trial Noise Model
Each simulation run introduces ±0.01° insertion-angle noise plus ±0.2% velocity-magnitude jitter, modelling combined sensor margin-of-error and thruster impulse precision. These produce the trial-to-trial variance that emerges in the recorded perilune deviations.
Three-Body Physics Note
The engine models Moon gravity directly plus a scaled Earth perturbation (5% of true Earth gravity) representing third-body effects. The reduced scale is a deliberate simplification that preserves the qualitative three-body character (non-Keplerian orbit shape, slow drift) while keeping the velocity → perilune-deviation signal visible across the 30-day window. A full N-body simulation would require station-keeping logic beyond the scope of this engine.