Thread Closed

Orbital velocities in the Schwartzschild geometry

 
Share Thread Thread Tools
Jul8-10, 07:39 PM   #528
 
Recognitions:
Science Advisor Science Advisor

Orbital velocities in the Schwartzschild geometry


Quote by espen180 View Post
Alright! I have now finished the calculation, thanks to kev's help and support!

Here is the newest document:
Download

I have removed the sections 2-4 which dealt with special cases, focusing on the general, unrestricted case. I abandoned the calculation I was currently working on (with ridicculously large expressions etc) for a simpler approach. Please take a look.

Edit: Updated with angular accelerations.
At the start of section 2, should d(theta)/dt=0 be d(theta)/d(tau) instead?

Also, I don't understand the comment "turning to the metric", just before Eqn 15, since the metric is Eqn 1?
 
Jul9-10, 12:45 AM   #529
 
Recognitions:
Science Advisor Science Advisor
Quote by atyy View Post
Also, I don't understand the comment "turning to the metric", just before Eqn 15, since the metric is Eqn 1?
I see. It's not the metric alone that you use to get Eqn 15. It's the metric plus the geodesic equation which implies that gab.dxa/dtau.dxb/dtau is a constant on a geodesic (following say Eq 2.4 and 2.5 of http://www.blau.itp.unibe.ch/lecturesGR.pdf; also see Altabeh's post #342)
 
Jul9-10, 02:38 AM   #530
 
Quote by starthaus View Post
Err, wrong. The above doesn't even make sense.



Basic calculus says that it doesn't follow.
Tell us which basic calculus supprts your nonsense here! LOL! I see you're completely bogged down with basics of calculus. Work hard!

[tex]r[/tex] is a coordinate while contrary to your fallacious claims [tex]R=R(\phi)[/tex] is not a "plane polar coordinate" but rather a trajectory. You are in no position to "teach" since you don't even know the difference between coordinates and trajectories. Please stop trolling this thread.
Nonsense. Read the page 196 of D'inverno to see why you were blindly correcting the author.

"... Then (15.20) can be integrated directly to give

[tex]r^2\dot{\phi}=h,[/tex]

where h is a constant. This is conservation of angular momentum (compare with (15.6) and note that, in the equatorial plane, the spherical polar coordinate is the same as the plane polar coordinate R).
"

Eq. (15.6) from the book is:

"[itex]R^2\dot{\phi}=h,[/itex]"

and

the equation (15.8) exactly shows that

"[itex]R=R(\phi)[/itex]."
Is you labelling us a "troller" another dead-end escape route towards not standing corrected? If so, then I don't want to say "you stop trolling" because I see that you're double finished by now.

AB
 
Jul9-10, 04:24 AM   #531
 
Quote by atyy View Post
At the start of section 2, should d(theta)/dt=0 be d(theta)/d(tau) instead?

Also, I don't understand the comment "turning to the metric", just before Eqn 15, since the metric is Eqn 1?
Thanks for pointing out, I've fixed it now.

I said "turning to the metric" instead of "from (1) we obtain" since I used the matric in a different form. (spacetime element vs. matrix).
 
Jul9-10, 04:55 AM   #532
 
Quote by espen180 View Post
Alright! I have now finished the calculation, thanks to kev's help and support!

Here is the newest document:
Download

I have removed the sections 2-4 which dealt with special cases, focusing on the general, unrestricted case. I abandoned the calculation I was currently working on (with ridicculously large expressions etc) for a simpler approach. Please take a look.

Edit: Updated with angular accelerations.
Good job espen. Way to go!

AB
 
Jul9-10, 05:20 AM   #533
 
I think the OP might be making life difficult for himself. With any form of variational problem, you should try to use Noether's theorem to integrate your equations up, rather than try to go from the Euler-Lagrange equations themselves. This will invariably make your life easier. Your Lagrangian density reads:

[tex] \mathcal{L} = \left( 1- \frac{1}{r}\right) \left( \frac{\mathrm{d} t}{\mathrm{d} \lambda}\right)^2 - \left( 1- \frac{1}{r}\right)^{-1} \left( \frac{\mathrm{d} r}{\mathrm{d} \lambda}\right)^2 - r^2 \left( \frac{\mathrm{d} \theta}{\mathrm{d} \lambda}\right)^2 - r^2 \sin^2\theta \left( \frac{\mathrm{d} \phi}{\mathrm{d} \lambda}\right)^2 [/tex]

using natural units. It is clear that the vector fields:

[tex] \frac{\partial}{\partial t} \quad \textrm{and}\quad \frac{\partial}{\partial \phi}[/tex]

are Killing, so Noether's theorem integrates up two of the Euler-Lagrange equations for you and gives you two constants of motion:

[tex] \left(1-\frac{1}{r}\right) \frac{\mathrm{d} t}{\mathrm{d}\lambda} = \mathrm{const} \,(=E) \quad \textrm{and} \quad r^2 \sin^2\theta \frac{\mathrm{d} \phi}{\mathrm{d} \lambda} = \mathrm{const} \,(=h) \qquad (*) [/tex]

i.e. on a given geodesic, these quantities remain unchanged. Similarly, since [tex]\partial_\lambda \mathcal{L}=0[/tex], we know [tex] \mathcal{L}[/tex] remains constant, and we set it to {+1,-1,0} depending on whether you're interested in timelike, spacelike or nulll geodesics. Call this constant k. Note that all our ODEs are now 1st order. Setting [tex]\theta = \pi/2[/tex] (validity can be deduced from the [tex]\theta[/tex] E-L equation) and using (*) in [tex]\mathcal{L}=k[/tex] gives the ODE:

[tex] \left( \frac{\mathrm{d} r}{\mathrm{d}\lambda}\right)^2 = E^2 - \left( 1-\frac{1}{r}\right) \left(k + \frac{h^2}{r^2}\right)[/tex]

If you'd prefer to parameterise your geodesics using [tex]\phi[/tex], use the second of the constraints in (*) again and you get:

[tex] \frac{h^2}{r^4} \left( \frac{\mathrm{d} r}{\mathrm{d} \phi}\right)^2 = E^2 - \left( 1-\frac{1}{r}\right) \left( k + \frac{h^2}{r^2}\right) [/tex]

If you'd prefer to do all this using tensors, just apply Noether's theorem in the form: if [tex]L_V g=0[/tex] (i.e. [tex]V[/tex] is a Killing vector) then [tex]V^\mu \dot{x}^\nu g_{\mu\nu} = \mathrm{const}[/tex].
 
Jul9-10, 05:27 AM   #534
 
Quote by Altabeh View Post
Good job espen. Way to go!

AB
Thanks! :)

And thanks for your hard work!
 
Jul9-10, 09:47 AM   #535
 
Quote by Anthony View Post
I think the OP might be making life difficult for himself. With any form of variational problem, you should try to use Noether's theorem to integrate your equations up, rather than try to go from the Euler-Lagrange equations themselves. This will invariably make your life easier. Your Lagrangian density reads:

[tex] \mathcal{L} = \left( 1- \frac{1}{r}\right) \left( \frac{\mathrm{d} t}{\mathrm{d} \lambda}\right)^2 - \left( 1- \frac{1}{r}\right)^{-1} \left( \frac{\mathrm{d} r}{\mathrm{d} \lambda}\right)^2 - r^2 \left( \frac{\mathrm{d} \theta}{\mathrm{d} \lambda}\right)^2 - r^2 \sin^2\theta \left( \frac{\mathrm{d} \phi}{\mathrm{d} \lambda}\right)^2 [/tex]

using natural units. It is clear that the vector fields:

[tex] \frac{\partial}{\partial t} \quad \textrm{and}\quad \frac{\partial}{\partial \phi}[/tex]

are Killing, so Noether's theorem integrates up two of the Euler-Lagrange equations for you and gives you two constants of motion:

[tex] \left(1-\frac{1}{r}\right) \frac{\mathrm{d} t}{\mathrm{d}\lambda} = \mathrm{const} \,(=E) \quad \textrm{and} \quad r^2 \sin^2\theta \frac{\mathrm{d} \phi}{\mathrm{d} \lambda} = \mathrm{const} \,(=h) \qquad (*) [/tex]

i.e. on a given geodesic, these quantities remain unchanged. Similarly, since [tex]\partial_\lambda \mathcal{L}=0[/tex], we know [tex] \mathcal{L}[/tex] remains constant, and we set it to {+1,-1,0} depending on whether you're interested in timelike, spacelike or nulll geodesics. Call this constant k. Note that all our ODEs are now 1st order. Setting [tex]\theta = \pi/2[/tex] (validity can be deduced from the [tex]\theta[/tex] E-L equation) and using (*) in [tex]\mathcal{L}=k[/tex] gives the ODE:

[tex] \left( \frac{\mathrm{d} r}{\mathrm{d}\lambda}\right)^2 = E^2 - \left( 1-\frac{1}{r}\right) \left(k + \frac{h^2}{r^2}\right)[/tex]

If you'd prefer to parameterise your geodesics using [tex]\phi[/tex], use the second of the constraints in (*) again and you get:

[tex] \frac{h^2}{r^4} \left( \frac{\mathrm{d} r}{\mathrm{d} \phi}\right)^2 = E^2 - \left( 1-\frac{1}{r}\right) \left( k + \frac{h^2}{r^2}\right) [/tex]

If you'd prefer to do all this using tensors, just apply Noether's theorem in the form: if [tex]L_V g=0[/tex] (i.e. [tex]V[/tex] is a Killing vector) then [tex]V^\mu \dot{x}^\nu g_{\mu\nu} = \mathrm{const}[/tex].
Exactly! This simply is based on the proposition given in post #389 and by this perfect explanation we are done here.

AB
 
Jul9-10, 05:26 PM   #536
 
I'm glad it was of some use (I don't have the will power to find post #389).
 
Jul9-10, 06:52 PM   #537
 
Mentor
Blog Entries: 10
It's on page 25, if you are using the forum default of 16 posts per page.
 
Thread Closed
Thread Tools


Similar Threads for: Orbital velocities in the Schwartzschild geometry
Thread Forum Replies
Strange geodesic in Schwartzschild metric Special & General Relativity 7
Orbital Velocity vs Mass and Orbital Velocity vs Radius Introductory Physics Homework 2
How do Orbits behave as orbital velocities become significantly relatavistic Special & General Relativity 8
Differential Geometry And Difference Geometry? Differential Geometry 1
Spherical Geometry and Flat Geometry Space-Spacetime's Special & General Relativity 0