Orbital velocities in the Schwartzschild geometry

In summary: Yes, it is very wrong.From the correct equation \frac{d^2\phi}{ds^2}=0 you should obtain (no surprise):\frac{d\phi}{ds}=constant=\omegaThe trajectory is completed by the other obvious equationr=R=constantYou get one more interesting equation, that gives u the time dilation. Start with:ds^2=(1-r_s/R)dt^2-(Rd\phi)^2 and you get:\frac{dt}{ds}=\sqrt{\frac{1+(R\omega)^2}{1-r_s/R}}or:\frac{ds}{dt}=\sqrt{1-r_s/R}\sqrt{1-\frac{(R
  • #456
starthaus said:
[tex]\frac{d}{dt}(K/\alpha)=\frac{d}{ds}(K/\alpha)\frac{ds}{dt}=-\frac{K}{\alpha^2}\frac{2m}{r^2}\frac{dr}{dt}[/tex]

See if you can figure the intermediate steps.
I deleted the question because 1) I am perfectly capable of carrying out the calculation myself and 2) no one is disputing that K is a constant with respect to s. The fact that I obtain the same result by treating K as a constant with respect to r is (to me) ample proof that K is a constant with respect to r.

I don't claim to derive everything from base principles. I take results from various authors and documents and extrapolate or extend them a little and draw my own conclusions and cross check my conclusions against known results. Sometimes I get it wrong and hope to learn from those times. The advantage of a forum is that you get feedback before you go too far off track. I am not a physicist or a mathematician. This is just a hobby for me and I like to pick the brains of people like Dalespam, Altabeh, Espen, etc, etc that are more advanced than me and see what I can learn from them and cross relate their results to my work and other references and generally build up "the big picture". That is just my way of doing things. Sorry if you do not like it.

Finally a question. Do you withdraw your claim that K is a function of r?

Yes, no or not sure?
 
Last edited:
Physics news on Phys.org
  • #457
starthaus said:
I wasn't going to put in all the intermediate calculations.
Of course you weren't going to. That would be too helpful.
 
  • #458
espen180 said:
I'm getting closer to a solution, but the expressions are ridicculously big, so doing the calculus/algebra/proofreading is a pain.

Here's how far I've gotten:
http://5554229043997450163-a-1802744773732722657-s-sites.googlegroups.com/site/espen180files/Schwartzschild.pdf?attachauth=ANoY7crvgkJaf2YSb9wipuKnvEoacyFDok-MWepVPK4tZrlLjK8If_HxjS_lzaZVPS-0StjbvJLQU3Bvx-t_oZcneyRG3Orj_V5UXTUFarMpmWYhgtEJNgxcRpypcyted9W0Qa7Re2gn-2ib6A0f-WJV5jukKhp_IfJhNoCnhscy-m2uBhjTO2cHlLX45JEfEI8kNFatSDfwlnBKGeus9Ufx5OdK13A-sw%3D%3D&attredirects=1" (Section 5, page 7)

I'm somewhat worried that I have an r-derivative raised to the 4th power. I can't see how it would disappear later on, either, so maybe I have an error somewhere.

Hi Espen,

From my earlier post:
kev said:
Re-insert the full form of H back into the equation:

[tex] \frac{d^2r}{ds^2} = -\frac{m}{r^2} + \frac{d\phi^2}{ds^2}(r-3m)[/tex]
(which is the same as (49) in your document) it would seem that (57) in your document should be:

[tex] \left(\frac{dr}{dt}\right)^2 \cdot \frac{d^2r}{ds^2} = \left(\frac{dr}{dt}\right)^2 \left(-\frac{m}{r^2} + (r-3m)\frac{d\phi^2}{ds^2}\right)[/tex]

using units of G=c=1.

You seem to have an extra (dt/ds)^2 expression on the right.

Equation (30) is fine in the context it is given in and as long as it not used to derive the general case (and you don't seem to have done so) it is not a problem.
 
Last edited by a moderator:
  • #459
kev said:
Hi Espen,

From my earlier post:
(which is the same as (49) in your document) it would seem that (57) in your document should be:

[tex] \left(\frac{dr}{dt}\right)^2 \cdot \frac{d^2r}{ds^2} = \left(\frac{dr}{dt}\right)^2 \left(-\frac{m}{r^2} + (r-3m)\frac{d\phi^2}{ds^2}\right)[/tex]

using units of G=c=1.

You seem to have an extra (dt/ds)^2 expression on the right.

I just have one question regarding this. Why are you multiplying d2r/ds2 by (dr/dt)2?

I used (dt/ds)2 to convert (dø/ds)2 to (dø/dt)2.
 
  • #460
espen180 said:
I just have one question regarding this. Why are you multiplying d2r/ds2 by (dr/dt)2?

I used (dt/ds)2 to convert (dø/ds)2 to (dø/dt)2.
My mistake. I misread the tau as an r. (That is why prefer the s symbol). Taking that into account and your conversion factor then it would seem that the first line of (57) is indeed correct.

I will go and polish my glasses and dig deeper. :wink: Sorry bout that :redface:
 
Last edited:
  • #461
kev said:
I deleted the question because 1) I am perfectly capable of carrying out the calculation myself and 2) no one is disputing that K is a constant with respect to s. The fact that I obtain the same result by treating K as a constant with respect to r is (to me) ample proof that K is a constant with respect to r.

It is ample proof that you are an incorigible hacker. You demanded to see how the equation was derived, I showed you the rigorous derivation , yet you persist in the hack.

BTW:

1. The equation I derived is not identical to your original one from post 277, it includes your post 277 as a particular case.

2. The equation in itself is pretty worthless since it is just a transformation from the simpler equation written in terms of proper time to uglier one written in terms of coordinate time.

3. The equation written in coordinate time can, through a lot of mathematical skill, lead to the equation of motion. I have challenged you , espen180 and Altabeh to find the expressions of [tex]r(s)[/tex], [tex]\frac{dr}{ds}[/tex], [tex]\frac{d^2r}{ds^2}[/tex], in other words to see if you can integrate the differential equation, if not for general orbits (way outside your skills) at least for radial motion. So far, nothing.
I don't claim to derive everything from base principles. I take results from various authors and documents and extrapolate or extend them a little and draw my own conclusions and cross check my conclusions against known results. Sometimes I get it wrong and hope to learn from those times. The advantage of a forum is that you get feedback before you go too far off track. I am not a physicist or a mathematician. This is just a hobby for me and I like to pick the brains of people like Dalespam, Altabeh, Espen, etc, etc that are more advanced than me and see what I can learn from them and cross relate their results to my work and other references and generally build up "the big picture". That is just my way of doing things. Sorry if you do not like it.

This is fine. What I do not understand is your reluctance to drop the hacks and to pick up rigorous approaches.
Finally a question. Do you withdraw your claim that K is a function of r?

Yes, no or not sure?

Why do you have to persist in horrible hacks? Can you ever learn how to derive things cleanly? I just showed you how to do it in the previous post.
 
  • #462
kev said:
My mistake. I misread the tau as an r. (That is why prefer the s symbol). Taking that into account and your conversion factor then it would seem that (57) is indeed correct. I think I might have been led on a "bum steer" by you know who. Should have known better. :P

Err, no. (57) is definitely wrong. It is amusing that you fail to see the error especially after a few minutes after I showed you the solution. So, I think that the "bum steer" is coming from you.
 
Last edited:
  • #463
starthaus said:
1. The equation I derived is not identical to your original one from post 277, it includes your post 277 as a particular case.
Dead wrong. my general equation, (derived long before your equation) contains all the information contained in your equation.
starthaus said:
2. The equation in itself is pretty worthless since it is just a transformation from the simpler equation written in terms of proper time to uglier one written in terms of coordinate time.
Equations in terms of proper time can get pretty ugly. The orbit of a photon for example, contains division by zero and infinities. It is much easier to measure the coordinate velocity of a particle and use the measurement directly in calculations, rather than have to also calculate the proper time time of the particle as well.
starthaus said:
3. The equation written in coordinate time can, through a lot of mathematical skill, lead to the equation of motion. I have challenged you , espen180 and Altabeh to find the expressions of [tex]r(s)[/tex], [tex]\frac{dr}{ds}[/tex], [tex]\frac{d^2r}{ds^2}[/tex], in other words to see if you can integrate the differential equation, if not for general orbits (way outside your skills) at least for radial motion. So far, nothing.
I think you are losing the plot here. Espen (the OP) has requested assistance with his calculation on orbital velocities. This is not the homework forum. It is not about setting a series of challenges with a sprinkling of dead-end diversions and red herrings. It is not just about your ego. We like to be collaborate and be cooperative in this forum. Your style would probably be better appreciated in the homework forum. Have you considered moving there?
 
  • #464
kev said:
Dead wrong. my general equation, (derived long before your equation) contains all the information contained in your equation.

Err, no. The general solution contains [tex]K[/tex] as a parameter, yours doesn't. Only if you set [tex]K[/tex] to particular values one recovers your equation. So, your solution is a subset of mine, like it or not.

Equations in terms of proper time can get pretty ugly.

Not the one that I am talking about. Compare it to the much uglier one expressed in coordinate time.
I think you are losing the plot here. Espen (the OP) has requested assistance with his calculation on orbital velocities.
Of which, despite yoour help he hasn't calculated any. It would be nice if you concentrated on helping espen180 in getting eqs (30) and (57) right.
 
Last edited:
  • #465
starthaus said:
It would be nice if you concentrated on helping espen180 in getting eqs (30) and (57) right.

Time and time again! You keep saying this! Yet, when asked to substansiate your claims, you back away, meaning you are either not in possesion of an answer (in other words, a groundless assertion) or purposefully withholding it. In both cases, such a claim has no merit.

Claim <-> Proof

These go hand in hand. You have yet to provide the latter.
 
  • #466
starthaus said:
Err, no. The general solution contains [tex]K[/tex] as a parameter, yours doesn't. Only if you set [tex]K[/tex] to particular values one recovers your equation. So, your solution is a subset of mine, like it or not.

So you substituted [itex] (\alpha/K) [/itex] for [itex] (ds/dt)[/itex]. Big deal!

I gave both versions, with and without constants:
kev said:
The metric is independent of [itex]\phi[/itex] and t, so there is a constant associated with coordinate angular velocity [itex](H_c)[/itex] which is obtained by finding the partial derivative of L with respect to [itex]d\phi/dt[/itex]

[tex]\frac{\delta{L}}{\delta(d\phi/dt)} = \frac{r^2}{\alpha} \frac{d\phi}{dt} = H_c [/tex]

The metric is independent of s and t, so there is a constant associated with time dilation [itex](K_c)[/itex] which is obtained by finding the partial derivative of L with respect to [itex]ds/dt[/itex]

[tex]\frac{\delta{L}}{\delta(ds/dt)} = \frac{1}{\alpha} \frac{ds}{dt} = K_c [/tex]
.
.
.
[tex]\frac{d^2r}{dt^2}= \frac{\alpha^2 H_c^2}{r^4}(r-5M) +\frac{M}{r^2}(2\alpha-3K_c^2\alpha^2)[/tex]

Re-inserting the full forms of [itex]H_c[/itex] and [itex]K_c[/itex] back in gives:

[tex]\frac{d^2r}{dt^2}=
\frac{d\phi^2}{dt^2}
(r-5M) +\frac{M}{r^2}(2\alpha-3
\frac{ds^2}{dt^2}
)[/tex]

Now the (ds/dt)^2 term is a little inconvenient, but we can find an alternative form by solving the Schwarzschild metric to directly obtain:

[tex] (ds/dt)^2 = (\alpha - dr^2/(\alpha dt^2) - r^2d\phi^2/dt^2)[/tex]

and substituting this form into the equation above it to obtain:

[tex]\frac{d^2r}{dt^2}=
\frac{d\phi^2}{dt^2}
(r-5M) +\frac{M}{r^2}(2\alpha-3
(\alpha - dr^2/(\alpha dt^2) - r^2d\phi^2/dt^2)
)[/tex]

which after a bit of algebra simplifies to:

[tex]\Rightarrow \frac{d^2r}{dt^2}= \alpha\left(
r \frac{d\phi^2}{dt^2}-\frac{M}{r^2}+
\frac{3M}{\alpha^2 r^2} \frac{dr^2}{dt^2}\right)[/tex]

This is the fully general form of acceleration in Schwarzschild coordinates in terms of coordinate time and gives Espen something to compare his results against in his latest document.
The third equation down contains only the single variable r and the constant parameters of the trajectory, [itex]K_c[/itex] and [itex]H_c[/itex]. In fact it is only if I set [itex]H_c[/itex] to particular values, that I recover your solution, so your solution is a subset of mine by your definition.
 
Last edited:
  • #467
espen180 said:
Time and time again! You keep saying this! Yet, when asked to substansiate your claims, you back away, meaning you are either not in possesion of an answer (in other words, a groundless assertion) or purposefully withholding it. In both cases, such a claim has no merit.

Claim <-> Proof

These go hand in hand. You have yet to provide the latter.

You need to look at your equations and compare with the ones already derived. Maybe kev can help you find the error in (57), it is quite obvious.
 
  • #468
kev said:
So you substituted [itex] (\alpha/K) [/itex] for [itex] (ds/dt)[/itex]. Big deal!

Err, no. Look again.
 
Last edited:
  • #469
starthaus said:
Maybe kev can help you find the error in (57), it is quite obvious.

It's exactly this I'm talking about. Learn some pedagogy!
 
  • #470
espen180 said:
It's exactly this I'm talking about. Learn some pedagogy!

There is value in working towards discovery, even if it is in order to discover your own mistakes. Demanding (in a very impolite manner) that others do your own work is not the right approach.
 
  • #471
starthaus said:
There is value in working towards discovery, even if it is in order to discover your own mistakes.

I can certainly appreciate such a philosophy, but saying "It's wrong" without further clarification is, and I hope you agree, insufficient when correcting someone. It simply makes it frustrating for the receiver when an answer is withheld, and I am certain that you'd feel the same had you been exposed to it.
 
  • #472
starthaus said:
You need to look at your equations and compare with the ones already derived. Maybe kev can help you find the error in (57), it is quite obvious.

Hi Espen, I have double checked and triple checked all 3 lines of (57) now and can assure you it is perfect. Ignore Starthaus. It is another red herring.
 
  • #473
starthaus said:
Err, no. The general solution contains [tex]K[/tex] as a parameter, yours doesn't. Only if you set [tex]K[/tex] to particular values one recovers your equation. So, your solution is a subset of mine, like it or not.

You have shot yourself in the foot there by declaring that a general solution is incomplete, if it does not contain the constant parameters of the trajectory. By that definition, your solution is incomplete. I will complete your solution for you:
starthaus said:
[tex]\frac{d^2r}{ds^2}=-\frac{m}{r^2}+(r-3m)(\frac{d\phi}{ds})^2[/tex]

Combine the above with the first Euler-Lagrange equation

[tex]\alpha\frac{dt}{ds}=K[/tex]

and, without any hacky assumptions about [tex]K[/tex] you will obtain that :

[tex]\frac{d^2r}{dt^2}=(r-3m)(\frac{d\phi}{dt})^2-(\frac{\alpha}{K})^2\frac{m}{r^2}+\frac{2m}{\alpha*r^2}(\frac{dr}{dt})^2[/tex]

The following can be directly obtained from the Schwarzschild metric:

[tex]\frac{dr^2}{dt^2} = \alpha \left(\alpha - \frac{ds^2}{dt^2} - r^2\frac{d\phi^2}{dt^2} \right) [/tex]

The second Euler-Langrange equation is

[tex]r^2 \frac{d\phi}{ds}=H[/tex]

Substitute the above expressions (along with the first Euler-Lagrange constant) into your incomplete solution and you obtain, after simplification;

[tex]\frac{d^2r}{dt^2}=\frac{\alpha^2}{r^4}\frac{H^2}{K^2} (r-5m) +\frac{m}{r^2}(2\alpha - 3\alpha^2 /K^2 ) [/tex]

There you go. You now have a fully complete solution just like the one I derived ages ago. It is not as nice, because you have a mixture of K and H in the first term while in my solution the constant parameters are cleanly isolated into separate terms.

The Lagrange constants in your solution, are related to my constants by the relations:

[tex]K_c = \frac{1}{K} = \frac{1}{\alpha} \frac{ds}{dt} [/tex]

and

[tex]H_c = \frac{H}{K} = \frac{r^2}{\alpha} \frac{d\phi}{dt} [/tex]
 
  • #474
kev said:
You have shot yourself in the foot there by declaring that a general solution is incomplete, if it does not contain the constant parameters of the trajectory. By that definition, your solution is incomplete. I will complete your solution for you:

LOL, give it a rest.
The following can be directly obtained from the Schwarzschild metric:

[tex]\frac{dr^2}{dt^2} = \alpha \left(\alpha - \frac{ds^2}{dt^2} - r^2\frac{d\phi^2}{dt^2} \right) [/tex]

The second Euler-Langrange equation is

[tex]r^2 \frac{d\phi}{ds}=H[/tex]

You are learning. Notice that [tex]\frac{dH}{dr}=2r\frac{d\phi}{ds}[/tex]. You need to give pause and think about this.
Substitute the above expressions (along with the first Euler-Lagrange constant) into your incomplete solution and you obtain, after simplification;

[tex]\frac{d^2r}{dt^2}=\frac{\alpha^2}{r^4}\frac{H^2}{K^2} (r-5m) +\frac{m}{r^2}(2\alpha - 3\alpha^2 /K^2 ) [/tex]

This is totally uunecessary, you did not understand an iota from the method of getting the general solution, you must NOT use all three of the Euler-Lagrange equations because only two are independent. I am quite sure I tried to teach you this long ago.
Do yourself a favor, learn the formalism, ok?
 
Last edited:
  • #475
kev said:
Hi Espen, I have double checked and triple checked all 3 lines of (57) now and can assure you it is perfect. Ignore Starthaus. It is another red herring.

Err, wrong. I'll give you a hint, compare the coefficients for [tex](\frac{d\phi}{dt})^2[/tex]
 
Last edited:
  • #476
starthaus said:
LOL, give it a rest.




You are learning. Notice that [tex]\frac{dH}{dr}=2r\frac{d\phi}{ds}[/tex]




This is totally uunecessary, you did not understand an iota from the method of getting the general solution, you must NOT use all three of the Euler-Lagrange equations because only two are independent. I am quite sure I tried to teach you this long ago.
Do yourself a favor, learn the formalism, ok?

After a day being off, I feel now refreshed because the fallacies you've been feeding all of us were just a few hours looked overshadowed by the beauties of nature.

Well, now back to our own discussion, give me the proof of proposition I gave you in post #389. If you couldn't make out proving it, then here is a hint: Use the geodesic equation [tex]\nabla_b u^a u_b=0[/tex].

AB
 
  • #477
kev said:
The following can be directly obtained from the Schwarzschild metric:

[tex]\frac{dr^2}{dt^2} = \alpha \left(\alpha - \frac{ds^2}{dt^2} - r^2\frac{d\phi^2}{dt^2} \right) [/tex]

The second Euler-Langrange equation is

[tex]r^2 \frac{d\phi}{ds}=H[/tex]
starthaus said:
You are learning. Notice that [tex]\frac{dH}{dr}=2r\frac{d\phi}{ds}[/tex]

LOL. That is an elementary calculus blunder. Shame on you!

Using the power rule:

[tex]\frac{dH}{dr}=2r\frac{d\phi}{ds} + r^2 \frac{d}{dr}\left(\frac{d\phi}{ds}\right) [/tex]

You need to prove

[tex]r^2 \frac{d}{dr}\left(\frac{d\phi}{ds} \right)= 0 [/tex]

That is only true if [itex](d\phi/ds)=0[/itex] or r=0 which it generally not the case for orbital motion, or if [itex](d\phi/ds)=0[/itex] is not a function of (r). We can easily show that [itex](d\phi/ds)[/itex] IS a function of r in the Newtonian limit when [itex]ds \approx t [/itex] because it is well known that angular momentum is conserved in Newtonian physics and [itex]d\phi/dt[/itex] changes with radius for a non circular orbit. You do realize GR has to agree with the Newtonian limit right? Try again.

kev said:
Hi Espen, I have double checked and triple checked all 3 lines of (57) now and can assure you it is perfect. Ignore Starthaus. It is another red herring.
starthaus said:
Err, wrong. I'll give you a hint, look at the coefficient for [tex]\frac{d\phi}{dt}[/tex]. Not your day today.

The equations in section (57) of Espen's document using G=c=1 are:

------------------------------------

[tex] \frac{ds^2}{dt^2} \cdot \frac{d^2r}{ds^2} = \frac{ds^2}{dt^2} \left(-\frac{M}{r^2} + (r-3M)\frac{dt^2}{ds^2}\frac{d\phi^2}{dt^2}\right)[/tex]

[tex] = -\frac{M}{r^2}\frac{ds^2}{dt^2} + (r-3M)\frac{d\phi^2}{dt^2}[/tex]

[tex] = -\frac{M}{r^2}
\alpha + \frac{M}{r^2\alpha}\frac{dr^2}{dt^2}
+ (r-2M)\frac{d\phi^2}{dt^2}[/tex]

------------------------------------

What you are missing is that from the Schwarzschild metric:

[tex]\frac{ds^2}{dt^2} = \alpha - \frac{1}{\alpha}\frac{dr^2}{dt^2} -{\color{red} r^2\frac{d\phi^2}{dt^2}} [/tex]


Using this information I can expand on Espen's work to make it clear that his work in section (57) contains no errors:

------------------------------------

[tex] \frac{ds^2}{dt^2} \cdot \frac{d^2r}{ds^2} = \frac{ds^2}{dt^2} \left(-\frac{M}{r^2} + (r-3M)\frac{dt^2}{ds^2}\frac{d\phi^2}{dt^2}\right)[/tex]

[tex] = -\frac{M}{r^2}\frac{ds^2}{dt^2} + (r-3M)\frac{d\phi^2}{dt^2}[/tex]

[tex] = -\frac{M}{r^2}\left(
\alpha - \frac{1}{\alpha}\frac{dr^2}{dt^2} {\color{red} -r^2\frac{d\phi^2}{dt^2}} \right)
+ (r-3M)\frac{d\phi^2}{dt^2}[/tex]

[tex] = -\frac{M}{r^2}
\alpha + \frac{M}{r^2\alpha}\frac{dr^2}{dt^2} {\color{red} + M\frac{d\phi^2}{dt^2}}
+ (r-3M)\frac{d\phi^2}{dt^2}[/tex]

[tex] = -\frac{M}{r^2}
\alpha + \frac{M}{r^2\alpha}\frac{dr^2}{dt^2}
+ (r-3M {\color{red} +M})\frac{d\phi^2}{dt^2}[/tex]

[tex] = -\frac{M}{r^2}
\alpha + \frac{M}{r^2\alpha}\frac{dr^2}{dt^2}
+ (r-2M)\frac{d\phi^2}{dt^2}[/tex]

------------------------------------

Get it now? Will you now stop telling Espen that (57) is wrong?

So in one day you have posted a elementary calculus blunder, an elementary physics blunder and an elementary algebra blunder and you were trying to tell me I was having a bad day? LOL.
 
  • #478
kev said:
LOL. That is an elementary calculus blunder. Shame on you!

Using the power rule:

[tex]\frac{dH}{dr}=2r\frac{d\phi}{ds} + r^2 \frac{d}{dr}\left(\frac{d\phi}{ds}\right) [/tex]

Err, [tex]r[/tex] and [tex]\phi[/tex] are independent variables, so your [tex]\frac{d}{dr}\left(\frac{d\phi}{ds}\right) [/tex] is pure nonsense. You really need to stop winging the calculus rules. Especially when multivariate functions are not covered by the sw package you are using in attempting your differentiations.
 
Last edited:
  • #479
starthaus said:
Err, [tex]r[/tex] and [tex]\phi[/tex] are independent variables...

No, not really. Sure, you can set the initial conditions freely, but after that, the geodesic equation, in which every coordinate, as well as its first and second derivatives are related to the other coodinates, dictates their behavior, so no coordinates are free from the influence of other coordinates after you set the initial conditions.
 
  • #480
starthaus said:
You really need to stop winging the calculus rules.

yowsa. pot-to-kettle.

There's only ONE independent variable in this problem. The worldine
coordinate s. So what we have are a bunch of functions of 1-variable.
r = r(s), t = t(s)... Now listen up if H were a function it would be a
function of s! (and only s).

When we write
[tex]r^2\frac{d\phi}{ds} = H[/tex]
What we mean is
[tex]r(s)^2 \frac{d\phi}{ds}(s) = H(s).[/tex]

When we're solving a trajectory problem in mechanics we
are solving a system of ODE(<- these are NOT PDE's).
For this problem, the system of ODE's that we need to solve is
[tex] \frac{d^2r}{ds^2} = \ldots,
\frac{d^2 t}{ds^2} = \ldots,
\frac{d}{ds} \left( r^2 \frac{d\phi}{ds} \right)= 0. [/tex]

This means H is independent of r, and [itex]\phi[/itex]. it's a CONSTANT.
 
  • #481
qbert said:
even though it's late in the game --

i think a point has come up which hasn't been answered very well. So I'll
take a stab at it. You need to be careful when using the calculus of
variations to keep in mind which variables are dependent and which are
independent. In the current case, the only independent variable is the
worldline coordinate (t or s). Everything else is a function of (say t) t.
That means H can be at most a function of t. period. since we know
dH/dt = 0 it is the constant function. meaning a real, honest to goodness,
constant.

If I've expressed myself clearly -- this is enough to prove the result.
The rest is for the formula lovers out there.
S'pose [itex]L = L(q_i(t), \dot{q}_i(t))[/itex] then by the E-L eqns we know
[tex] \frac{d}{dt}\left( \frac{\partial L}{\partial \dot{q}_j} \right)
= \frac{\partial L}{\partial q_j}[/itex]

Now suppose that L is cyclic in one variable say [itex]q_k[/itex]. That is
[tex] \frac{\partial L}{\partial q_k} = 0
= \frac{d}{dt} \left( \frac{\partial L}{\partial \dot{q}_k} \right)[/tex]

So that [itex] \frac{\partial L}{\partial \dot{q}_k} = C[/itex]
Now suppose to the contrary that C were really a function of the "variables"
ie C = f(qi).

Then we get
[tex] \frac{\partial f}{\partial q_j} = \frac{\partial^2 L}{\partial q_j \partial \dot{q}_k}
= \frac{\partial^2 L}{\partial \dot{q}_k \partial q_j}[/tex]

[tex] = \frac{\partial }{\partial \dot{q}_k} \left( \frac{\partial L}{\partial q_j} \right)
= \frac{\partial}{\partial \dot{q}_k}\left( \frac{d}{dt} \frac{\partial L}{\partial \dot{q}_j} \right) = \frac{\partial}{\partial \dot{q}_j}\left( \frac{d}{dt} \frac{\partial L}{\partial \dot{q}_k} \right) = \frac{\partial}{\partial \dot{q}_k}\left( 0 \right) = 0.[/tex]

This is one of the bestr answers I have seen in this thread. Unfortunately, you haven't taken the time to understand the dispute between me and kev. The derivation that you show above is based on the Euler-lagrange equations, so, your results are valid only only along the arc defined by:

[tex]ds^2=\alpha dt^2-dr^2/\alpha-(rd\phi)^2[/tex]

kev is trying to hack his derivation by differentiating as if [tex]H[/tex] and [tex]K[/tex] are constant everywhere. This is obviously not true. To wit, in your example [tex]f(r)=r^2\frac{d\phi}{ds}[/tex] so, [tex]\frac{df}{dr}[/tex] isn't zero.
 
  • #482
qbert said:
yowsa. pot-to-kettle.

There's only ONE independent variable in this problem. The worldine
coordinate s. So what we have are a bunch of functions of 1-variable.
r = r(s), t = t(s)... Now listen up if H were a function it would be a
function of s! (and only s).

When we write
[tex]r^2\frac{d\phi}{ds} = H[/tex]
What we mean is
[tex]r(s)^2 \frac{d\phi}{ds}(s) = H(s).[/tex]

Good, so you just confirmed that kev's

[tex]\frac{dH}{dr}=2r\frac{d\phi}{ds} + r^2 \frac{d}{dr}\left(\frac{d\phi}{ds}\right) [/tex]

is nonsense. Thank you
 
Last edited:
  • #483
kev said:
[tex] = -\frac{M}{r^2}
\alpha + \frac{M}{r^2\alpha}\frac{dr^2}{dt^2}
+ (r-2M)\frac{d\phi^2}{dt^2}[/tex]

------------------------------------

Get it now? Will you now stop telling Espen that (57) is wrong?

.

Err, wrong, post 430 gives you the correct (unhacked) solution. As you can see, the correct coefficient is [tex]r-3M[/tex] not [tex]r-2M[/tex].
 
  • #484
kev said:
LOL. That is an elementary calculus blunder. Shame on you!

Using the power rule:

[tex]\frac{dH}{dr}=2r\frac{d\phi}{ds} + r^2 \frac{d}{dr}\left(\frac{d\phi}{ds}\right) [/tex]

You need to prove

[tex]r^2 \frac{d}{dr}\left(\frac{d\phi}{ds} \right)= 0 [/tex]
starthaus said:
Err, [tex]r[/tex] and [tex]\phi[/tex] are independent variables, so your [tex]\frac{d}{dr}\left(\frac{d\phi}{ds}\right) [/tex] is pure nonsense. You really need to stop winging the calculus rules. Especially when multivariate functions are not covered by the sw package you are using in attempting your differentiations.

Note: I have completely re-edited this post to reflect the reminder by Espen that the anti-derivative of zero is a constant.

If as you claim, r and [itex]\phi[/itex] are independent variables, then it follows that

[tex] \frac{d}{dr}\left(\frac{d\phi}{ds}\right) =0 \qquad \qquad (1)[/tex]

Taking the anti-derivative of both sides with respect to s gives:

[tex]\frac{d\phi}{ds} = C_{(r)} \qquad \qquad (2)[/tex]

where [itex]C_{(r)}[/itex] is a constant with respect to r.

We agree that H is constant with respect to s so we obtain:

[tex]\frac{dH}{ds}= \frac{d}{ds}\left(r^2\frac{d\phi}{ds}\right) =0 [/tex]

[tex]\Rightarrow \frac{d\phi}{ds}\frac{d}{ds}(r^2) +r^2 \frac{d}{ds}\left(\frac{d\phi}{ds}\right) =0 [/tex]

[tex]\Rightarrow \frac{d\phi}{ds}\frac{d}{ds}(r^2) +r^2 \frac{dr}{ds}\frac{d}{dr}\left(\frac{d\phi}{ds}\right) =0 \qquad \qquad (3)[/tex]

Substituting (2) into (3) gives:

[tex] \frac{d\phi}{ds}\frac{d}{ds}(r^2) +r^2 \frac{dr}{ds}\frac{d}{dr}\left(C_{(r)}\right) =0 [/tex]

[tex] \Rightarrow \frac{d\phi}{ds}\frac{d}{ds}(r^2) =0 [/tex]

[tex] \Rightarrow \frac{d}{ds}(r^2) =0 [/tex]

Taking the anti-derivative of both sides with respect to s gives:

[tex] \Rightarrow r^2 = C_{(s)} [/tex]

where [itex]C_{(s)} [/itex] is a constant with respect to proper time s.

This shows that Starthaus's claims are only valid for circular orbits where r is constant.
 
Last edited:
  • #485
The antiderivative is ...=constant, but that shouldn't invalidate your proof. In fact, you can use it to prove that H is a constant.
 
  • #486
starthaus said:
Err, wrong, post 430 gives you the correct (unhacked) solution. As you can see, the correct coefficient is [tex]r-3M[/tex] not [tex]r-2M[/tex].


In post #430 you give:
starthaus said:
[tex]\frac{d^2r}{ds^2}=-\frac{m}{r^2}+(r-3m)(\frac{d\phi}{ds})^2[/tex]


Espen gives:
[tex] \frac{ds^2}{dt^2} \cdot \frac{d^2r}{ds^2} = -\frac{M}{r^2}
\alpha + \frac{M}{r^2\alpha}\frac{dr^2}{dt^2}
+ (r-2M)\frac{d\phi^2}{dt^2}[/tex]


Do you not see that the left hand sides of the equations are different and so by elementary algebra we expect the right hand sides to be different?

You are talking about two different quantities. Wrong again.
 
  • #487
kev said:
In post #430 you give:



Espen gives:



Do you not see that the left hand sides of the equations are different and so by elementary algebra we expect the right hand sides to be different?

You are talking about two different quantities. Wrong again.

Err, no. Bottom of post 430 says:

[tex]\frac{d^2r}{dt^2}=(r-3m)(\frac{d\phi}{dt})^2-(\frac{\alpha}{K})^2\frac{m}{r^2}+\frac{2m}{\alpha*r^2}(\frac{dr}{dt})^2[/tex]
 
  • #488
starthaus said:
Err, no. Bottom of post 430 says:

[tex]\frac{d^2r}{dt^2}=(r-3m)(\frac{d\phi}{dt})^2-(\frac{\alpha}{K})^2\frac{m}{r^2}+\frac{2m}{\alpha*r^2}(\frac{dr}{dt})^2[/tex]

Writing out the [tex]\frac{\alpha}{K}[/tex] term, what do you get?
 
  • #489
kev said:
and I have proven:

[tex] \frac{d}{dr}\left(\frac{d\phi}{ds}\right) \ne 0 [/tex]

kev,

I understand that you did not understand calculus very well and you are desperately trying to cover your mathematical blunder. So, let's try an "intuitive" explanation, what you "proved" above is that the angular speed depends on radius. Do you grasp the enormity of your error?
 
  • #490
espen180 said:
Writing out the [tex]\frac{\alpha}{K}[/tex] term, what do you get?

What is the puropse of this nonsensical request? The derivation is based on the fact that (see post 430)

[tex]\frac{ds}{dt}=\frac{\alpha}{K}[/tex]
 

Similar threads

Replies
4
Views
959
Replies
9
Views
1K
Replies
3
Views
2K
Replies
13
Views
4K
Replies
28
Views
1K
Replies
21
Views
1K
Back
Top