I I made an N-body simulation program

Click For Summary
A project focused on solving the 3-body gravity problem evolved into a comprehensive N-body simulator that incorporates both gravitational and electromagnetic interactions, utilizing special relativistic mechanics. The program allows for unique simulations, such as modeling a hydrogen atom alongside a microscopic black hole, demonstrating advanced numerical integration techniques. Users are encouraged to explore the simulation's capabilities, with feedback welcomed for further improvement. The discussion also highlights the complexity of N-body simulations and shares resources for additional learning. Overall, the project showcases innovative programming and physics integration.
Random Physics Fan
Messages
1
Reaction score
7
My little project - 3body simulator started as a programming exercise, which I attempted to make a program that solves the famous 3-body gravity problem using numerical integration.
Later I gained interest in the underlying physics and expanded the program to solve N-body problems involving both electromagnetic and gravitational interactions. I also replaced Newtonian mechanics with special relativistic mechanics (as for the model of gravity, general relativity is too challenging for me, so I used an approximation under the weak field limit which was Lorentz-invariant).
You can do all sorts of weird simulation such as the following:

Object 1 (e-) and object 2 (p+) form a hydrogen atom under the classical Bohr's model, with n = 10; object 3 is a hypothetical microscopic blackhole with the mass of the Moon, placed 10 mm under the hydrogen atom.

Object 1:
mass: 9.11e-31 / charge: -1.60e-19
position: 5.29e-9,0,0 / velocity: 0,218500,0

Object 2:
mass: 1.67e-27 / charge: 1.60e-19
position: 0,0,0 / velocity: 0,0,0

Object 3:
mass: 7.35e22 / charge: 0
position: 0,0,-0.01 / velocity: 0,0,0

Duration of simulation: 1e-12 / number of steps: 1000
-> Use advanced option and turn on scientific notation (to show the position of the electron)

(Note: the Schwarzschild radius of the blackhole is about 0.11mm (<<10mm), so general relativistic effects like gravitational time dilation are insignificant.)

Please try out my little program and you comments are appreciated!
 
  • Like
Likes davyzhu, Bandersnatch, vanhees71 and 4 others
Physics news on Phys.org
Hi,

Congratulation. I did an N-body simulation (yours is N=3) some time ago. It is really not easy nor it is dead, in my eyes. Perhaps you are interested in the following material: N-Body

Have fun!
 
  • Like
Likes aaroman and vanhees71
Great work @Random Physics Fan and @Omega0

While I can't comment on the accuracy of your simulations, from a programming standpoint, I'm sure others can benefit from what you folks did.

As an aside, there are also some good simulations at the Open Source Physics website done in Java using the OSP library of simulation routines.

www.compadre.org/osp

And with that said, its a good time to close this thread.

Jedi
 
Thread 'Why higher speeds need more power if backward force is the same?'
Power = Force v Speed Power of my horse = 104kgx9.81m/s^2 x 0.732m/s = 1HP =746W Force/tension in rope stay the same if horse run at 0.73m/s or at 15m/s, so why then horse need to be more powerfull to pull at higher speed even if backward force at him(rope tension) stay the same? I understand that if I increase weight, it is hrader for horse to pull at higher speed because now is backward force increased, but don't understand why is harder to pull at higher speed if weight(backward force)...

Similar threads

  • · Replies 15 ·
Replies
15
Views
3K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 19 ·
Replies
19
Views
4K
Replies
3
Views
1K
  • · Replies 19 ·
Replies
19
Views
3K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K