... | ... | @@ -12,12 +12,6 @@ Box of size 12cm. Things go crazy when the particles drift off the edge of the b |
|
|
|
|
|
![energy_change_l_1_drift_speed_0](/uploads/5707ac2534646767c3d35dcfb586bf4/energy_change_l_1_drift_speed_0.png)
|
|
|
|
|
|
![energy_change_l_1_drift_speed_1](/uploads/fa522d5256cc000ae8847afac9075e8c/energy_change_l_1_drift_speed_1.png)
|
|
|
|
|
|
![energy_change_l_1_drift_speed_2](/uploads/ce4a4ced1d4a1c7a7bb1ebd9703c9916/energy_change_l_1_drift_speed_2.png)
|
|
|
|
|
|
![energy_change_l_1_drift_speed_3](/uploads/ea22369a82bdea43e0ed9d6b223e64c8/energy_change_l_1_drift_speed_3.png)
|
|
|
|
|
|
A line of particles with a total mass of 1 solar mass is frozen in place. Another particle is put into a circular orbit at 1AU in a plane normal to the line, with the plane bisecting the line. Below I show the change in kinetic, potential and total energy, which should all stay the same.
|
|
|
eta = 0.01 and epsilon = 0.01
|
|
|
First plot is simply the two-body problem again
|
... | ... | @@ -31,11 +25,21 @@ First plot is simply the two-body problem again |
|
|
|
|
|
Here I show the effect of decreasing eta on conservation of total energy (epsilon = 0.01)
|
|
|
|
|
|
![line_of_10_particles_eta_test](/uploads/32397720099e5df21c60e6e6f2b67860/line_of_10_particles_eta_test.png)
|
|
|
![eta_test_perp_orbit](/uploads/59a33b1e2a88403c425b2b545ab039ea/eta_test_perp_orbit.png)
|
|
|
|
|
|
Same test but instead I put the orbiting particle in the same plane as the line
|
|
|
|
|
|
![eta_test_plane_orbit](/uploads/14f07f4973aeb85501d3b826adb6620c/eta_test_plane_orbit.png)
|
|
|
|
|
|
Changing the mass of the orbiting particle has no effect:
|
|
|
|
|
|
![mass_test_plane_orbit](/uploads/03b7d3e7d7530729fd5bd43ad9fa6dbd/mass_test_plane_orbit.png)
|
|
|
|
|
|
Another set-up.
|
|
|
|
|
|
3 frozen particles, each with mass of 1/3 solar mass, in an equilateral triangle of side-length 1 AU. One of the sides is on the x-axis. An orbiting particle is given an initial velocity equal to that in the previous example.
|
|
|
|
|
|
![line_of_10_particles_plane_eta_test](/uploads/efce48af686059fb8b702d8ccffca39d/line_of_10_particles_plane_eta_test.png)
|
|
|
![eta_test_plane_orbit_triangle](/uploads/3dc370df4ce18be20da2c1898d139584/eta_test_plane_orbit_triangle.png)
|
|
|
is in-line with the x-axis.
|
|
|
|
|
|
|