From Taylor 1st Edition, problem 8.25, Classical Mechancis: Consider a particle

Question

From Taylor 1st Edition, problem 8.25, Classical Mechancis:

Consider a particle with mass m and angular momentum l in the field of a central force F= ? k/(r^(5/2)). To simplify your equations, choose units for which m=l=k=1.

b) Assuming now that the particle has energy E=-0.1 find an accurate value of rmin, the particle's distance of closest approach to the force center.

c) Assuming that the particle is a r= rmin when ?=0 use a computer program to solve the transformed radial equation and find the orbit in the form r=r(?) for 0???7?. Plot the orbit. Does it appear to be closed?

Explanation / Answer

The transformed radial equation:

u''()=-u()- l 2 u ( ) 2 l2u()2F

comes from the radial equation of the orbit

r . . r..=F(r)+ l 2 r 3 l2r3

by making the substitution u=1/r.

The attempt at a solution

I have completed everything but the final plot of the orbit and I'm not sure why it is giving me trouble. To do part a) simply find where the derivative of the effective potential is equal to zero and solve for r. Using the suggested simplifications for the units it is easy to find that r 0 r0=1. I plotted the effective potential in Mathematica along with the energy of the particle by defining

Ueff (r)=- 2 /(3 r3/2) + 1/2r2 and then using

Plot[{Ueff[r],-0.1},{r,0,5}].

From the plot it is easy to see that the distance of closest approach is around 0.7 and if you find the roots you will see that it is r 0 r0=0.6671. Now, for part c) you need to write F(r) as F(u) using the substitution mentioned above. After simplification

u''()= u ( ) u()-2u().

If this can be solved for u() then you can easily find r(). I tried using Mathematica again to numerically solve and then plot the reciprocal of my solution. Part of the problem is that I wasn't sure what to use for my second initial condition, u'(0), but I tried several different numbers just to see if I could get a reasonable plot and I did not. Anyway, here is what I entered

s=NDSolve[{u''[]== u [ ] u[]-2u[], u[0]=1.499, u'[0]=1}, u, {, 0, 7}]

Here, u[0] comes from the fact that u=1/r and r0=0.667. Once again, I just picked a number for u'[0] but I have tried various selections. To plot the solution PolarPlot

[Evaluate[1/u[]/.s] ,{,0,7}]

The result is not the orbit (not closed) that is pictured in the back of the book, it is more or less a straight line. Not sure where I went wrong here and any help would be appreciated. To check that the technique was correct I did check my work with other examples that I knew the solution to (a free particle, and a particle under the influence of an inverse square force).

Related Questions

Navigate

• Browse (All)
• Browse F
• Subjects

Integrity-first tutoring: explanations and feedback only — we do not complete graded work. Learn more.