... | ... | @@ -4,20 +4,19 @@ Compare SWIFT solution with analytic solution, and the solution from the scipy O |
|
|
|
|
|
Show the effect of decreasing 'eta'.
|
|
|
|
|
|
![point_mass_eta_test_r_3_m_1](/uploads/b0024b951379bc78d08ed6cda3aba4ba/point_mass_eta_test_r_3_m_1.png)
|
|
|
![point_mass_eta_test_r_3_m_1](/uploads/fb61ce85bc42eb07dce412813c7c3b6b/point_mass_eta_test_r_3_m_1.png)
|
|
|
|
|
|
Now solve a similiar problem but replace the central particle with a line of 10 particles, from -0.5 AU to 0.5 AU along the x-axis, with total mass of one solar mass. Initial position and velocity of the orbiting test particle is still the same.
|
|
|
|
|
|
Can test the effect of changing both the multipole order and 'eta'
|
|
|
|
|
|
![particle_line_m_order_test_r_3_minus_log_10_eta_3](/uploads/bad8f189532a5580bfb2e32a1dbb2de4/particle_line_m_order_test_r_3_minus_log_10_eta_3.png)
|
|
|
![particle_line_m_test_r_3_minus_log_eta_3](/uploads/e033aac6cb4d87f7fddde909feb608d8/particle_line_m_test_r_3_minus_log_eta_3.png)
|
|
|
|
|
|
![particle_line_eta_test_r_3_m_1](/uploads/b2d8b5907a9cd20b468bc5201d9f579b/particle_line_eta_test_r_3_m_1.png)
|
|
|
|
|
|
Now do the same this again but with a triangle of particles, centred on [0.5,0.5,0.5], evenly spaced around a circle of radius 0.3 AU in the x-y plane, this was done to ensure the particles were all in the same cell.
|
|
|
|
|
|
|
|
|
![triangle_m_order_test_r_3_m_2](/uploads/3c09f6b9e5550b9aad6618a22f3d8d11/triangle_m_order_test_r_3_m_2.png)
|
|
|
![triangle_m_test_r_3_minus_log_eta_3](/uploads/5ccba8f7367cafe7b817962d5a83947a/triangle_m_test_r_3_minus_log_eta_3.png)
|
|
|
|
|
|
![triangle_eta_test_r_3_m_1](/uploads/f9b0c781abe92d74baba943d0d99b58b/triangle_eta_test_r_3_m_1.png)
|
|
|
|
... | ... | |