Phil Lucht Math & Physics Archive
Home / Math and Physics Files / Physics / Mechanics / frames doc stuff / Reviewed Docs

goldstein v equation attempt 2 on 7_12_12

DOCX · 24.8 KB
Open DOCX file

Dated 3.26.05 in its text and apparently revised 7/12/12, this is a Word document of Phil's own derivation with bracketed self-commentary. It computes the acceleration in rotating-frame terms using the rotation matrix R(t) and its time derivative, giving an extra term ω x v. It then splits the Goldstein equation into radial and angular components and checks them, verifying the first but ending with r = 0 and noting something is still wrong. Symbols are partly lost in the extraction.

AI-written summary; may contain errors. This description is approximate.

Extracted text (machine-read; may contain errors)
This is the Title PhL 3.26.05 Note that page numbering is turned on in this template. [ This entire document is an example of : garbage in, garbage out ] (a) A Pre-calculation. Recall from above that vbody(t) = R(t) v'body(t) [ what is vbody ? ] [ OK, let's imagine that v'body is the velocity of our ant as measured in the rotating frame S'. Later I will call this vector v'S' = v' = (dr'/dt)S' with components (v')'i in the rotating frame. What then is vbody in the above equation? It cannot be vS= v because we know the equation is not right for that. And v'S = (dr'/dt)S which from the G Rule we know is v' – ω x r' and that doesn't seem to match our equation either. The fourth is vS' = (dr/dt)S' and that seems pretty distant. So I think vbody is none of the actual four velocities in the problem. It is just an actively rotated version of v' . So why do a calculation on something that you don't know what is ? I think I am just continuing my erroneous ways. Recall how we keep getting that "extra factor of R" in the G Rule which I though you had to then set to R = 1] Then dt vbody(t) = (dtR) v'body(t) + R(t) dtv'body(t) or abody(t) = (dtR) v'body(t) + R(t) a'body(t) . And we must carefully take note that we don't just have abody(t) = R(t) a'body(t). In fact abody is not a "vector" due to the above transformation rule with its extra piece. Using our usual method, [(dtR) v'body(t)]i = [(dtexp(-iφ(t)Jz)) v'body(t)]i = {[-iφ'(t) Jz] [(exp(-iφ(t)Jz) ) v'body(t)] }i // i (Jk)ij = kij = {[-iφ'(t) Jz] [Rv'body(t)] }i = {[-iφ'(t) Jz] [vbody(t)] }i = -iφ'(t) (Jz)ij [vbody(t)]j = - φ'(t) ε3ij [vbody(t)]j = - ω3 ε3ij [vbody(t)]j = εi3j ω3[vbody(t)]j = εikj ωk[vbody(t)]j So we have shown that [(dtR) v'body(t)]i = εikj ωk[vbody(t)]j [(dtR) v'body(t)] = ω x vbody(t) Therefore, going way back we then have abody(t) = (dtR) v'body(t) + R(t) a'body(t) = R(t) a'body(t) + ω x vbody(t) [ and here is that fully factor of R mentioned above ] and so we have restarted our "extra term". (b) The Goldstein equation for v . That equation says this dv = dφ x v + dvbody dv/dt = ω x v + dvbody/dt [ OK, if correctly interpreted: dv/dt = ∂Sv and dvbody/dt = ∂S'v ] a = ω x v + abody Installing our pre-calculation result, this says a = ω x v + R(t) a'body(t) + ω x vbody(t) We may regard our first equation as a first order PDE in variable v(t) dv/dt – ω x v = R(t) a'body(t) + ω x vbody(t) We think we know everything on the RHS, but let's make sure. For the first term we get, a'body = dv'body /dt = dt[' ' + r' ' '] = five terms = ' ' + ' ' ' + ' ' ' + r'' ' + r' ' [- ' '] = (' - r''2) ' + (r'' + 2 ' ' ) ' and then (this removes the primes from the unit vectors) R(t) a'body(t) = (' - r''2) + (r'' + 2 ' ' ) For the second term we have, vbody = ' + r' ' // from Section 3 ω x vbody(t) = [ω] x [' + r' ' ] = ω' - ω r' ' The combined RHS is then (' - r''2) + (r'' + 2 ' ' ) + ω ' - ω r' ' = (' - r''2- ω r' ') + (r'' + 2 ' ' + ω ') Therefore, our Goldstein equation is this dv/dt – ω x v = (' - r''2- ω r' ') + (r'' + 2 ' ' + ω ') Now we want to get the components for the terms on the LHS: v = vr + vθ ω = ω => ω x v = [ω] x [vr + vθ] = ω vr – ω vθ dv/dt = r + vr + θ + vθ[- ] = (r - vθ) + (vr + θ) So here is our grand equation (r - vθ) + (vr + θ) – ω vr + ω vθ = (' - r''2- ω r' ') + (r'' + 2 ' ' + ω ') Separating by components, we get these two equations (r - vθ) + ω vθ = (' - r''2- ω r' ') (vr + θ) – ω vr = (r'' + 2 ' ' + ω ') or (r - vθ)+ ω vθ = (' - r''2- ω r' ') (vr + θ)– ω vr = (r'' + 2 ' ' + ω ') or r - vθ + ω vθ = ' - r''2- ω r' ' vr + θ – ω vr = r'' + 2 ' ' + ω ' or r - (-ω) vθ = ' - r''2 – ω r' ' θ + ( – ω)vr = r'' + 2 ' ' + ω ' And now the two equations are not decoupled! Things are much more complicated than before. No doubt we could solve these equations somehow, but I am not going to try. My main interest is to verify that these equations are true!! From the previous section we found that v(t) = vr + vθ where vr = vθ = r So we can say vr = r = vθ = r θ = r + Then our two coupled equations above become r - (-ω) vθ = ' - r''2 – ω r' ' θ + ( – ω)vr = r'' + 2 ' ' + ω ' or - (-ω) r = ' - r''2 – ω r' ' r + + ( – ω) = r'' + 2 ' ' + ω ' Let's whittle away on the first equation: - (-ω) r = ' - r''2 – ω r' ' - (-ω) r = - r'2 – ω r ' // since r = r' - (-ω) r = - r'2 – ω r ' - ' r = - r'2 – ω r ' - ' = - '2 – ω ' - = - ' – ω = ' + ω and this is true, so we have verified our first equation ! Now let's whittle away at the second equation r + + ( – ω) = r'' + 2 ' ' + ω ' r + + ( – ω) = r' + 2 ' + ω // since r = r' r + + – ω = r' + 2 ' + ω r + 2 – ω = r' + 2 ' + ω r ( –' ) + 2 ( - ') = 2ω r + 2 ω = 2ω r = 0 Close but something is still wrong.