PDA

View Full Version : Spherical Oscillator


T-7
Mar24-08, 04:41 PM
Hi,

I'd really appreciate some speedy advice on this question. I need to know whether the expressions for energy and the angular momentum are correct, as I seem unable to produce the final expression desired from what I have derived.

Here's the problem:

1. The problem statement, all variables and given/known data

The Lagrangian is

L = \mu \nu^{2}/2 - kr^{2}/2

i. Calculate the angular momentum vector, L, proving that L is constant and that, unless L vanishes, the oscillation occurs in the x,y plane.

ii. Using x = x_{0}cos(\omega t), y = y_{0}sin(\omega t), show the trajectory is an ellipse. Express the energy, E, and then both the eccentricity and the parameter in terms of E and Lz.

iii. Representing the trajectory in polar coordinates by x = rcos(\phi), y = rsin(\phi), show that you get the simple relation


cos (2\phi) = \frac{1 - \nu^{2} / (2\xi^{2})}{\sqrt{1 - \nu^{2}}} ,
r = \sqrt{ \frac{2E}{k} }\xi ,
\nu = \frac{ L_{z}\omega }{E}


3. The attempt at a solution

For i. I found that L = Lz = \mu \omega x_{0}y_{0}sin(\phi_{x}-\phi_{y}) (assuming L is not zero), by using three sine equations for the x,y and z components (with angle offsets) and taking L = m(r \times v); I chose my axis so that x and y could be rewritten as x = x_{0}cos(\omega t), y = y_{0}sin(\omega t), and then find I can write Lz = \mu \omega x_{0}y_{0}.

For ii., it seemed all that was required was to square both expressions and combine thus:

\frac{x}{x_{0}}^{2} + \frac{y}{y_{0}}^{2} = 1

Writing the energy E = \mu \nu^{2}/2 + kr^{2}/2 and substituting the expressions for x and y, I found the expression to boil down to E = \frac{(\mu \omega^{2})}{2} (x_{0}^{2}+y_{0}^{2}) where \omega^{2} = k/m.

I'm guessing I need to go ahead and solve this equation for x0 (assuming it's right) by using Lz to get rid of y0, and hope that something which looks like the eccentricity and the parameter will emerge (?).

Part iii. has me really stuck. I substitute x = rcos(\phi), y = rsin(\phi) into the ellipse equation, which looks like a promising start, use double angle results to turn the cos and the sine into cos (2\phi), and eventually obtain (after some simply rearranging):

cos (2\phi) = \frac{1-r^{2}/2(1/x_{0}^{2}+1/y_{0}^{2})}{r^{2}/2(1/x_{0}^{2}-1/y_{0}^{2})}

Multiplying numerator and denomiator by x_{0}y_{0} and using my expressions derived earlier for Energy and Lz I arrive at:

cos (2\phi) = \frac{L_{z}^{2}/\mu^{2}\omega^{2}-r^{2}/2(2E/\mu\omega^{2})}{r^{2}/2(y_{0}^{2}-x_{0}^{2})}

It looks promising, but once the first term on the numerator has been made one, it seems no manipulation can make this equation correspond to the one given. I can only presume my energy, angular moment or ellipse equations are wrong? But are they?

Many thanks for your insights!

joe8232
Mar26-08, 09:00 AM
Hey, Can I just ask how you got the expression for E=\frac{(\mu\omega^2)}{2}(x_0^2+y^2_0)

Ruddiger28
Mar26-08, 01:22 PM
I've been working on this problem as well, and I think that the expression for the energy should have a k term

E=\frac{(\mu\omega^2)}{2}(x_0^2+y^2_0)+\frac{k^2}{ 2}}(x_0^2+y^2_0)

Because if you put the terms for xo and yo in for r I can't see why the k term would disappear. Also for your expression for the ellipse the denominators should be squared.

I haven't gone any further with this yet, but the method seems like it should work from what you've explained.

joe8232
Mar26-08, 01:55 PM
Can I just ask if we are all handing this in tomorrow for a certain classical mechanics module at a certain Scottish University!?!

I tried using Maple to get the expression for E which comes out as:
E=\frac{k}{2}(x_0^2+y_0^2)

I can't seem to get very far with the next question and haven't even looked at the last one. Oh Classical Mechanics!

has anyone had any luck getting the eccentricity out?