Complete position function of time for multibody systems and gravity (Newtonian)

Click For Summary
SUMMARY

The discussion focuses on deriving a position function of time for multiple bodies in a gravitational system, particularly in a solar system context. The user seeks to extend basic Newtonian mechanics, specifically the equation Y(t) = Yo + Vo(t) + [a/2](t)^2, to account for variable gravitational acceleration as bodies move further from Earth. The conversation highlights the complexity of the N-body problem, emphasizing that a closed-form solution does not exist for more than two bodies, while a solution does exist for the two-body problem involving point masses, leading to orbital elements.

PREREQUISITES
  • Basic understanding of Newton's laws of motion
  • Familiarity with calculus concepts, particularly derivatives and integrals
  • Knowledge of gravitational force equations, specifically Fg = GMm/y^2
  • Understanding of the N-body problem in classical mechanics
NEXT STEPS
  • Research the closed-form solution for the Newtonian two-body problem and its implications for orbital mechanics
  • Study numerical methods for solving the N-body problem, such as the Runge-Kutta method
  • Explore the concept of orbital elements and their calculation for celestial bodies
  • Investigate simulation tools for modeling gravitational interactions in multi-body systems, such as MATLAB or Python libraries
USEFUL FOR

Students and professionals in physics, astrophysics, and engineering who are interested in gravitational dynamics, orbital mechanics, and the complexities of multi-body systems.

bleist88
Messages
6
Reaction score
0
Hi,

{First off I should share my ultimate goal, which is to have a position function of time of several bodies, let's say in a solar system, so that the position of each body and the acceleration due to gravity on all the other bodies in the system are given at any time t. Although Kepler's formulas are useful, I am not wanting only to understand how bodies in orbits work, but also bodies that are randomly thrown through the system such as the Voyager stations or meteors or what have you.

So to take this problem into steps and build up to this:

There is first being able to describe a rock dropped close to Earth. Here, acceleration of gravity is the constant 9.8m/s. This is simple and takes less thought.

Next, Imagine dropping this body from much higher where the gravity is weaker and so the acceleration is smaller, and as the body falls, the acceleration is not constant. This is where I get stuck.

Next, I think, would be to imagine now that body given some x component of velocity and to determine the path it would take in the ( x, y, z, t ) as well as the ( r ) it would have.

This would give the path the moon takes but also the path a meteor close to Earth would take if it fell into Earth.

Next would be an even greater challenge, to Imagine a multi-body system so that the mass of both bodies would be considerable and be constantly changing all the others position. This would involve not only a function based on changin acceleration, but also changing position of the major gravity players.}

I understand basic Calculus, which will obtain a position function of time based on Newton's F=ma and the constant acceleration of gravity close to the Earth's surface: Y(t) = Yo + Vo(t) + [a/2](t)^2

...where Yo is original height, Vo is original dy/dt, and a is the constant acceleration of gravity.

Where I am stuck is the next ste, where you would drop a ball close to Earth, yet far enough away so that the acceleration due to gravity is no longer constant but smaller, and a function on the distance from the center of the Earth the object is...

Fg = GMm/y^2 = ma ... so... a = GM/y^2

But the problem lies where you would apply this to the former equation I used when gravity was constant and unchanging. (Now Y is not just the height but the distance from the two objects centers)

It gives: Y(t) = Yo + Vo(t) + {GM/2[Y(t)]^2}(t)^2...

This is a bit ugly in computer type, but basically the same as before except that the acceleration is not a constant but another function based on Y. But Y(t) is dependent on a(t) which is really a function of Y(t). This means the function is a function of itself, which at my misunderstanding of Calculus makes no sense. What did Newton do from here?



This is the question I have and I would love for it to be answered. I have searched through my Dynamics book, Calculus, and gone to professors and still have yet to get further.
 
Physics news on Phys.org
You are looking for a closed form solution to the N-body problem. It doesn't exist. Google "N-body problem" and you will see.

There is a closed form solution to the Newtonian two body problem -- if the two bodies are point masses. This will lead to things like orbital elements (I suggest you google that phrase as well).
 

Similar threads

  • · Replies 9 ·
Replies
9
Views
3K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 6 ·
Replies
6
Views
1K
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
6
Views
2K
  • · Replies 14 ·
Replies
14
Views
22K
  • · Replies 10 ·
Replies
10
Views
2K