My God, It’s Full Of Stars: The N-Body-Problem Series

I thought it would be a very nice idea to start my blog with a small project I wanted to do since a long time. A small tool/solver of the n-body-problem. It’s not just only the interest of solving a challenging mathematical problem with modern C++. I also love clean, concise and readable code, which is also something I would like to train and share with this project. You can find the empty repository setup at GitHub. I will always try to tag the end state of the repository after each post of this blog series.

So lets start with the bare minimum of theory we need. We want to start with a simple example of two still standing masses (or stars if you like) which get attracted to each other. For the sake of simplicity we consider them as point masses which mean they are infinitesimal small, or in other words, they have no volume.

From the image above we can see that a particle is moving the path , with its velocity and accelerated by . From this point we get the equations of motion of one such point as

with as the resulting position, as the starting position, as the velocity, as the time step and as the acceleration. We also know that the force is defined by the product of mass and the acceleration

We know as well that the force due to gravity between two masses can be calculated as

with the gravitational constant the two masses of them and their distance . With the equations of the force and the force due to gravity, we can derive the acceleration a mass point is experiencing while he is attracted by another mass point.

As a result we get our final equation which is deriving the change of position of a mass point by its starting position, its velocity and its acceleration.

With this I will close this post and we will proceed next time with setting up a basic CMAKE example project which is building our solver library and a separate test library.

Yes and no, yes they have no volume, but they have mass. These objects are called point masses and their property is to be infinite small (singularity) but have full mass of the object. Therefore m > 0

Is there any reason you don’t simplify your expression for “a” so that is “a = G*(m_2 / r^2)” ? (sorry, I don’t know if your comment boxes allow LaTeX or similar)! Excited to see the rest of these articles

If the objects have no volume, doesn’t that mean that your last equation has a divide by zero?

LikeLike

Yes and no, yes they have no volume, but they have mass. These objects are called point masses and their property is to be infinite small (singularity) but have full mass of the object. Therefore m > 0

LikeLike

Is there any reason you don’t simplify your expression for “a” so that is “a = G*(m_2 / r^2)” ? (sorry, I don’t know if your comment boxes allow LaTeX or similar)! Excited to see the rest of these articles

LikeLike

Just for sake of readability. You’re right, clearly m_1 has no meaning in the equation anymore, and it would have been also better to name a into a_1.

LikeLike

Oh good, I was worried I greatly missed something. I do like the subscript idea, but I’m joining your party late so don’t let me get in the way.

LikeLike