Undergrad 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
 
Topic about reference frames, center of rotation, postion of origin etc Comoving ref. frame is frame that is attached to moving object, does that mean, in that frame translation and rotation of object is zero, because origin and axes(x,y,z) are fixed to object? Is it same if you place origin of frame at object center of mass or at object tail? What type of comoving frame exist? What is lab frame? If we talk about center of rotation do we always need to specified from what frame we observe?

Similar threads

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