Theory of the N-Body Problem
June 9, 1996
23
1.5 Types of Deviations In Star Movement
There are several different types of deviations from the ideal star movement that
show up from the discretization and the round off errors. Most of these errors come from
having one or more of the constants of motion changing over time.
1.5.1 Gaining or Losing Energy.
When a star gains or loses energy, instead of the star making an elliptical orbit
around a collapsar, the orbits decay and either spiral in or spiral out. So, instead of a
perfect ellipse, the star movement looks like FIGURE 15.
This type of error tends to make a very uninteresting star system because a lot of
stars will be lost either due to collisions or by stars moving off the screen. Euler's method
tends to gain energy over time, but the fourth order Adam-Bashford method tends to very
slowly lose energy.
1.5.2 Distortion
Some methods, such as the taylor3 method, both gains and loses energy, depending
on the situation. When a star orbits a collapsar, it keeps it's elliptical orbit, but the orbit is
more elongated than it should be. It is hard to even notice this unless you monitor the total
energy level as the system progresses.
It is hard to predict what this type of error will do to a system. With only two stars
it is hard to tell that this type of error even exists, but with many stars it is clear that the
taylor3 method does not give as accurate of results as other methods.
1.5.3 Perihelion Shift
Some methods maintain the elliptical orbit around a collapsar, but the perihelion
(the spot on the ellipse closest to the collapsar) rotates around the collapsar over time.
Energy is neither gained nor lost over time, but the angular momentum changes.
FIGURE 15. Results of a method that gains/loses energy.
Previous page    Home    Table of Contents    Next page
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57
This document is best viewed as n-body.pdf because the translation to html was buggy.