Two-body equation of motion resolution

Click For Summary
SUMMARY

The discussion centers on the two-body vector differential equation of motion, represented as $$\ddot{\textbf{r}}+\frac{\mu}{r^3}\textbf{r}=0$$. This equation is classified as a non-linear differential equation, which complicates its analytical solution. While a partial solution exists that describes the size and shape of the orbit as a conic section, a complete solution, which includes time as a variable, is not straightforward and requires advanced techniques such as iterative calculations or geometric methods. The conversation emphasizes the importance of understanding both the mechanics of planar particle motion and Kepler's laws for a comprehensive grasp of the topic.

PREREQUISITES
  • Understanding of ordinary differential equations (ODEs)
  • Familiarity with Kepler's laws of planetary motion
  • Knowledge of mechanics of planar particle motion
  • Basic calculus and iterative methods for solving equations
NEXT STEPS
  • Study the mechanics of planar particle motion in detail
  • Learn about Kepler's equation and its applications
  • Explore numerical methods for solving transcendental equations
  • Investigate the Euler method for simulating orbits
USEFUL FOR

This discussion is beneficial for physics students, mathematicians, and anyone interested in celestial mechanics or the mathematical modeling of orbital dynamics.

greg_rack
Gold Member
Messages
361
Reaction score
79
Hi guys,

I have started studying differential equations on my own, sticking to my last high-school year's textbook, along with a few physics applications of ODEs.
Online I came across the n-body problem, which then took me to the basic two-body problem!

I'm here to ask you a few things about the two-body vector differential equation of motion:
$$\ddot{\textbf{r}}+\frac{\mu}{r^3}\textbf{r}=0$$

First: which type of DE is it? My textbook covers(approximately) linear DEs, explaining just a few resolution methods for first and second order homogeneous and non-homogeneous ones.
Is this non-linear, and so not solvable with "usual" and easier methods? If yes, why is it so?

Second: the textbook I found says: "though this equation looks easy, its complete solution is not".
Why?

Lastly, after the aforementioned consideration, it provides the steps to find a "partial solution that will tell us the size and shape of the orbit", which is a function ##r##(rel. position of the two bodies) of the form of a conic section: this "partial solution" thus tells us -only- the shape and size of the orbit, but doesn't provide us with a function of time for the movement.
How is that? What is a "partial solution", and what would be needed to find a "complete" one?

Online resources look all very confusing and over-complicated, and I'm having a really hard time deciphering those... hope you can help me, and my questions aren't too silly :)
 
Physics news on Phys.org
Hi,

greg_rack said:
First: which type of DE is it?
It took me a bit of searching to find your textbook. Next time, perhaps you can insert the link in your post ? Saves others time :smile: !

And on top of page 8 you missed the answer:

1621374656829.png

The assumptions that follow and that lead to ##\ddot{\textbf{r}}+\frac{\mu}{r^3}\textbf{r}=0## don't change that.

Is this non-linear, and so not solvable with "usual" and easier methods? If yes, why is it so?

Second: the textbook I found says: "though this equation looks easy, its complete solution is not".
Why?
I could shrug this off with 'ours is not to reason why' or 'I really don't know ' . Why things are as they are is not physics. There are no tools in our toolkit to solve this analytically.

Suffice it to point to Kepler's equation ##M = E - e\sin E\ ## that remains after a whole lot of chewing.
One could argue that it's so difficult because the effective one-dimensional potential in which the motion takes place is such an unpleasant function of the radial distance ( ##-\displaystyle{a\over r} + \displaystyle{b\over r^2}\ ## ). Does that help you ?

The partial solution is 1-41 or 1-42 ( ##r = \displaystyle{p\over 1+e\cos\nu} \ ## ). The complete solution would be ##{\textbf{r}(t)}##, for which the book refers to Chapter 4, over a hundred pages further. You have your work cut out for you :wink: !

greg_rack said:
all very confusing and over-complicated, and I'm having a really hard time deciphering those... hope you can help me, and my questions aren't too silly :)
Your questions are not silly and I admire your courage and stamina to tackle something like this. But I think you are trying to bite off more than you can chew. I can't help in making it easier. The authors do a good job and I can't rewrite their book.

My advice is twofold:
If your angle is physics stick with Mechanics_of_planar_particle_motion and Kepler's laws
And if it's differential equations, find an undergraduate calculus book you like and do the exercises.

https://www.physicsforums.com/threads/free-math-books.796225/
https://www.physicsforums.com/threads/free-physics-books.796223/

##\ ##
 
Last edited:
  • Like
Likes   Reactions: greg_rack and Ibix
BvU said:
Next time, perhaps you can insert the link in your post ? Saves others time :smile: !
Definitely will do! Haven't really thought about it before posting 😅

BvU said:

The assumptions that follow and that lead to ##\ddot{\textbf{r}}+\frac{\mu}{r^3}\textbf{r}=0## don't change that.
Great, and what does "defied solution" mean? I came across, I got the sense, but wasn't really able to get the precise meaning... I'm not English, and no translator gave me a satisfying translation for "defied solution"

BvU said:
One could argue that it's so difficult because the effective one-dimensional potential in which the motion takes place is such an unpleasant function of the radial distance ( ##-\displaystyle{a\over r} + \displaystyle{b\over r^2}\ ## ). Does that help you ?
That does definitely help, thanks a lot!

BvU said:
The partial solution is 1-41 or 1-42 ( ##r = \displaystyle{p\over 1+e\cos\nu} \ ## ).
Got it, but what is exactly meant by "partial solution"? Is it merely a way of solving the equation which gives the sought radius as a scalar, and not as ##{\textbf{r}(t)}##?

Heartfelt thanks for your kindness, answer, and advice... it's been precious! :)
 
greg_rack said:
defied solution
It's a part of the verb "to defy". "It has resisted our attempts at finding a solution" is an expanded translation. Or "we can't solve it".
 
  • Like
Likes   Reactions: greg_rack
greg_rack said:
what is exactly meant by "partial solution"
Is is only a part of the solution, namely the radial distance as a function of an angle that varies in time.
So it gives the shape of the trajectory.
For a complete solution you would still need how that angle varies in time.

##\ ##
 
  • Like
Likes   Reactions: greg_rack
Because the equation of the two body problem is what is known as a transcendental equation, and those often do not have a closed form solution. https://en.wikipedia.org/wiki/Transcendental_equation.

It is solvable though, but only using iterative calculation or geometry. If you just want to create some neat orbits then i suggest you use the Euler method instead of Keplers. Since the latter is probably a bit too difficult to start with.
 
Last edited:
  • Like
Likes   Reactions: greg_rack
  • Like
Likes   Reactions: greg_rack
Thank you so much guys!
 

Similar threads

  • · Replies 7 ·
Replies
7
Views
1K
  • · Replies 1 ·
Replies
1
Views
992
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 21 ·
Replies
21
Views
4K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 2 ·
Replies
2
Views
10K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
3K
Replies
10
Views
2K