Derivation of planetary orbit equation with Lagrangian

Click For Summary
SUMMARY

The discussion centers on the derivation of the planetary orbit equation using Lagrangian mechanics, specifically addressing the constancy of angular momentum represented by the equation $$\dot{φ}=\frac{l}{\mu{r^2}}$$. The participants clarify that the constancy arises from the Lagrange equation, where $$\frac{\partial L}{\partial \varphi}=0$$ leads to the conclusion that $$\frac{d}{dt}\frac{\partial L}{\partial \dot\varphi}=0$$. Additionally, they explore the general solution of the differential equation $$\frac{\mathrm{d}^2 w}{\mathrm{d} \varphi^2}=-w$$, confirming that it can be expressed as $$w(\varphi)=A \cos(\varphi+\delta)$$, where $$A$$ and $$\delta$$ are constants derived from integration constants $$A_1$$ and $$A_2$$.

PREREQUISITES
  • Understanding of Lagrangian mechanics and the Lagrange equations.
  • Familiarity with differential equations, particularly second-order linear equations.
  • Knowledge of trigonometric identities and their applications in physics.
  • Basic concepts of angular momentum in classical mechanics.
NEXT STEPS
  • Study the derivation of Lagrangian equations in classical mechanics.
  • Learn about the properties of second-order linear differential equations.
  • Explore the relationship between angular momentum and orbital mechanics.
  • Investigate the implications of conic sections in planetary motion.
USEFUL FOR

Students and professionals in physics, particularly those focusing on classical mechanics, orbital dynamics, and mathematical physics. This discussion is beneficial for anyone looking to deepen their understanding of Lagrangian mechanics and its applications in celestial mechanics.

TimeRip496
Messages
249
Reaction score
5
I am stuck at this part page 1,
$$\frac{\partial{L}}{\partial{\dot{φ}}}=\mu{r^2}\dot{φ}=const=l------->\dot{φ}=\frac{l}{\mu{r^2}}......(8)$$
Why is this a constant? Isn't r and dφ/dt variables of time?

Source: http://www.pha.jhu.edu/~kknizhni/Mechanics/Derivation_of_Planetary_Orbit_Equation.pdf
 
Last edited by a moderator:
Physics news on Phys.org
it is because ##\frac{\partial L}{\partial \varphi}=0## and due to the Lagrange equation: ##\frac{d}{dt}\frac{\partial L}{\partial \dot\varphi}=\frac{\partial L}{\partial \varphi}## one gets ##\frac{d}{dt}\frac{\partial L}{\partial \dot\varphi}=0##
 
  • Like
Likes   Reactions: vanhees71 and TimeRip496
wrobel said:
it is because ##\frac{\partial L}{\partial \varphi}=0## and due to the Lagrange equation: ##\frac{d}{dt}\frac{\partial L}{\partial \dot\varphi}=\frac{\partial L}{\partial \varphi}## one gets ##\frac{d}{dt}\frac{\partial L}{\partial \dot\varphi}=0##
ok Thanks a lot!
I have one more question as to how did the source get this w(φ) = Acos(φ+d)? The equation is below eqn(14) which is stated in here " which has the well known solution w(φ) = Acos(φ+d), where both A and are constants. We can always choose d= 0 by a convenient choice of φ"
 
Last edited:
The most general solution of the differential equation
$$\frac{\mathrm{d}^2 w}{\mathrm{d} \varphi^2}=-w$$
is given by
$$w(\varphi)=A_1 \cos \varphi+A_2 \sin \varphi,$$
where ##A_1## and ##A_2## are constants of integration. This follows, because for a homogeneous 2nd-order linear differential equation the solution is the linear combination of two arbitrary linearly independent solutions, and cos and sin fulfill these conditions.

Now you can write this general solution also in another form,
$$w(\varphi)=A \cos(\varphi+\delta),$$
where ##A## and ##\delta## are constants. To see this just use the addition theorem for cos
$$w(\varphi)=A \cos \delta \cos \varphi - A \sin \delta \sin \varphi.$$
Comparing with the general solution above you have
$$A_1=A \cos \delta, \quad A_2=-A \sin \delta,$$
and you always find ##A>0## and ##\delta \in (-\pi,\pi]## to fulfill these equations. First we have
$$A=\sqrt{A_1^2+A_2^2}, \quad \delta =-\mathrm{sign} A_2 \arccos \left (\frac{A_1}{\sqrt{A_1^2+A_2^2}} \right).$$
So indeed the given solution is the most general one for the differential equation for ##w##.
 
vanhees71 said:
The most general solution of the differential equation
$$\frac{\mathrm{d}^2 w}{\mathrm{d} \varphi^2}=-w$$
is given by
$$w(\varphi)=A_1 \cos \varphi+A_2 \sin \varphi,$$
where ##A_1## and ##A_2## are constants of integration. This follows, because for a homogeneous 2nd-order linear differential equation the solution is the linear combination of two arbitrary linearly independent solutions, and cos and sin fulfill these conditions.

Now you can write this general solution also in another form,
$$w(\varphi)=A \cos(\varphi+\delta),$$
where ##A## and ##\delta## are constants. To see this just use the addition theorem for cos
$$w(\varphi)=A \cos \delta \cos \varphi - A \sin \delta \sin \varphi.$$
Comparing with the general solution above you have
$$A_1=A \cos \delta, \quad A_2=-A \sin \delta,$$
and you always find ##A>0## and ##\delta \in (-\pi,\pi]## to fulfill these equations. First we have
$$A=\sqrt{A_1^2+A_2^2}, \quad \delta =-\mathrm{sign} A_2 \arccos \left (\frac{A_1}{\sqrt{A_1^2+A_2^2}} \right).$$
So indeed the given solution is the most general one for the differential equation for ##w##.
Thanks a lot for the reply! But i am stuck at this part
and you always find ##A>0## and ##\delta \in (-\pi,\pi]## to fulfill these equations. First we have
$$A=\sqrt{A_1^2+A_2^2}, \quad \delta =-\mathrm{sign} A_2 \arccos \left (\frac{A_1}{\sqrt{A_1^2+A_2^2}} \right).$$
So indeed the given solution is the most general one for the differential equation for ##w##.
Isn't $$\quad \delta=arccos(\frac{A_1}{A})=-arcsin(\frac{A_2}{A})$$?
Besides do you mean that $$w(\varphi)=A \cos(\varphi+\delta),$$ is more general than $$w(\varphi)=A_1 \cos \varphi+A_2 \sin \varphi,$$? If so, I can't see the reason behind it.
 
No, both equations are equivalent, as demonstrated in the previous posting. You can always map from the pair ##(A_1,A_2)## of integration constants to the pair ##(A,\delta)##. It just depends on what you consider more convenient. In this case, it's easier to see that the orbit is a conic section (ellipse, parabola, hyperbola), using the ##(A,\delta)## notation.
 
  • Like
Likes   Reactions: TimeRip496
vanhees71 said:
No, both equations are equivalent, as demonstrated in the previous posting. You can always map from the pair ##(A_1,A_2)## of integration constants to the pair ##(A,\delta)##. It just depends on what you consider more convenient. In this case, it's easier to see that the orbit is a conic section (ellipse, parabola, hyperbola), using the ##(A,\delta)## notation.
Oh I see. But how did you obtain this $$\quad \delta =-\mathrm{sign} A_2 \arccos \left (\frac{A_1}{\sqrt{A_1^2+A_2^2}} \right).$$?
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 11 ·
Replies
11
Views
2K
  • · Replies 19 ·
Replies
19
Views
4K
  • · Replies 1 ·
Replies
1
Views
1K
Replies
6
Views
2K
  • · Replies 1 ·
Replies
1
Views
601
  • · Replies 21 ·
Replies
21
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K