Took a 4th stab (or perhaps 5th, I've lost count) at simulating n-body physics in a #bevy plugin that supports time acceleration. This one finally seems to have solved all the orbital stability and energy conservation issues that have plagued my past attempts.
Check out the post here, and see the link at the bottom for an interactive example.