Orbital velocities in the Schwartzschild geometry

  • Thread starter Thread starter espen180
  • Start date Start date
  • Tags Tags
    Geometry Orbital
Click For Summary
The discussion revolves around calculating the velocity of a particle in a circular orbit around a black hole using the tensor formulation of General Relativity (GR) and the Schwarzschild geometry. The original poster encounters a zero velocity result when applying the geodesic equations, raising concerns about potential miscalculations or misunderstandings. Participants point out a sign error in the equations and emphasize the importance of correctly interpreting the assumptions made, particularly regarding fixed radial coordinates. The conversation highlights the need for clarity in defining orbital velocity and its relationship to angular velocity in the context of GR, ultimately leading to the conclusion that the equations derived must accurately reflect the dynamics of circular motion in a relativistic framework.
  • #91
espen180 said:
I was writing #85 when you posted #84, so pardon me for not noticing it.

The equation in #53, after imposing dø/ds=0:

-d/ds(1/\alpha*2dr/ds)-(2m/r^2)(dt/ds)^2-d/dr(1/\alpha)(dr/ds)^2=0

\frac{\text{d}}{\text{d}s}\frac{2}{\alpha}\frac{\text{d}r}{\text{d}s}=2\frac{\text{d}r}{\text{d}s}\frac{\text{d}}{\text{d}s}\frac{1}{\alpha}+\frac{2}{\alpha}\frac{\text{d}^2r}{\text{d}s^2}

It may be because I'm quite tired, but I think I need a push on solving this.

You are on the right path (pay attention to the last term, the derivative wrt d/dr, you may have made a mistake , if you do the calculation right there will be some terms that will cancel out, at the end you will be rewarded by a very interesting (and correct) equation of motion.

An even easier approach, as I mentioned before, start with:

L=\alpha(dt/ds)^2-(dr/ds)^2/\alpha

and form the Euler-Lagrange equation wrt dr/ds and r (it is the same as one of the geodesic equations). We can continue tomorrow.
 
Last edited:
Physics news on Phys.org
  • #92
starthaus said:
What in \frac{d^2r}{dt^2}=\frac{d}{dt}(\frac{dr}{dt}) do you still struggle with?
What is the value of \frac{dr}{dt} for a particle in freefall at the moment of apogee?
 
Last edited by a moderator:
  • #93
Al68 said:
What is the value of \frac{dr}{dt} for a particle in freefall at the moment of apogee?

If you don't understand ordinary differential equations, I can recommend a few good courses. In the meanwhile, try trolling other threads.
 
  • #94
starthaus said:
Al68 said:
What is the value of \frac{dr}{dt} for a particle in freefall at the moment of apogee?
If you don't understand ordinary differential equations, I can recommend a few good courses. In the meanwhile, try trolling other threads.
LOL. You're such a sweetheart. How can you be so helpful and so pleasant at the same time? Thank you so much for your brilliant, congenial, and non-condescending answer.
 
  • #95
Al68 said:
LOL. You're such a sweetheart. How can you be so helpful and so pleasant at the same time? Thank you so much for your brilliant, congenial, and non-condescending answer.

If you spent less time trolling and more time studying you would have known that, given the ODE:

\frac{d^2r}{dt^2}+A\frac{dr}{dt}+B=0 for any a<t<b

if you make \frac{dr}{dt}=0 for any a<t<b

this means

\frac{d^2r}{dt^2}=0 any a<t<b

meaning that:

B=0

So, in general, \frac{dr}{dt}=0 is not a solution for the ODE. This is why we try solving the ODE with known methods rather than doing silly things like setting \frac{dr}{dt}=0.

Even worse is the naive attempt by certain members of this forum at calculating \frac{d^2r}{dt^2} by inserting \frac{dr}{dt}=0 in the above ODE and declaring that \frac{d^2r}{dt^2}=-B.

Have you slept through your calculus classes or you never took any?
Now, if you could go troll other threads and leave me to help espen180 find the answer to his second question in this thread, that would be nice. We are almost done with the solution, if you have something to contribute, you are welcome but if you only plan to create trolling noise , than troll elsewhere.
 
Last edited:
  • #96
starthaus said:
If you spent less time trolling and more time studying you would have known that, given the ODE:

\frac{d^2r}{dt^2}+A\frac{dr}{dt}+B=0 for any a<t<b

if you make \frac{dr}{dt}=0 for any a<t<b

this means

\frac{d^2r}{dt^2}=0 any a<t<b

meaning that:

B=0

So, in general, \frac{dr}{dt}=0 is not a solution for the ODE. This is why we try solving the ODE with known methods rather than doing silly things like setting \frac{dr}{dt}=0.

Have you slept through your calculus classes or you never took any?
Now, if you could go troll other threads and leave me to help espen180 find the answer to his second question in this thread, that would be nice.
Now you're just being too much of a sweetheart. How can you be so non-condescending and provide such a direct and specific answer to my question? I'm sure espen180 will appreciate such non-condescending and straightforward answers as much as I did.
 
  • #97
Al68 said:
Now you're just being too much of a sweetheart. How can you be so non-condescending and provide such a direct and specific answer to my question? I'm sure espen180 will appreciate such non-condescending and straightforward answers as much as I did.

He does. Especiially since I answer his questions.
 
  • #98
kev said:
Dear starthaus,

Since you seem to have a very poor grasp of the most elementary physics, let us go back back to Newtonian physics and review the basics.

Consider a ball thrown vertically upwards. To a good aproximation, if the ball is not thrown too high, the acceleration is 9.8m/s. At the apogee (the maximum height of the ball's trajectory or even more basically when the ball stops going upwards and starts falling back down again) the average velocity is zero or in the infinitesimal limit the velocity is momentarily exactly zero. At this point dr/dt = 0 and the acceleration is non-zero (d^2r/dt^2 = 9.8m/s. If the acceleration was zero at the apogee the ball would remain at its maximum height and not fall down again. You can experimentally prove that this is not the case in your own back garden.

Now in the Schwarzschild metric dr/dt=0 means the vertical velocity of the test particle is momentarily zero at a given instant. (Recall that by definition dt means and infinitesimal interval of time - Please review an introductory textbook on calculus if you have forgotten this basic fact.) Setting dr/dt=0 says nothing about the motion of the particle in the next instant or in the previous instant. Therefore setting dr/dt=0 in the metric does not by itself, imply anything about (as you put it) whether dr/dt=0 is an "initial condition" or a "general condition". The only way you can determine if it is an initial condition (I prefer a momentary condition) or a general condition (I prefer to say a condition that does not vary over time) is by determining whether or not the acceleration of the particle is non-zero or zero. Setting dr/dt equal to zero, does not by itself imply d^2r/dt^2 must also be zero. If you think this is mathematically impossible, then you need to refresh your basic math skills as well as your basic physics knowledge.

Completely seconded. In particular we must add to this the fact that in https://www.physicsforums.com/blog_attachment.php?attachmentid=167&d=1275053682" the author bizarrly argues that d^2r/ds^2 is the proper acceleration of a radially freely falling particle and this has been explained to him several times by many people over other threads but instead of accepting "the error" he continued backing his wrong idea with no consequences.

Physics is not just "mathematical skills"; it requires you to have a skill of painting a physial picture first (the thing that you really sound to be unfamiliar with) and then coloring it with the mathematical hues.

starthaus said:
If you spent less time trolling and more time studying you would have known that, given the ODE:

\frac{d^2r}{dt^2}+A\frac{dr}{dt}+B=0 for any a<t<b

if you make \frac{dr}{dt}=0 for any a<t<b

this means

\frac{d^2r}{dt^2}=0 any a<t<b

meaning that:

B=0

So, in general, \frac{dr}{dt}=0 is not a solution for the ODE. This is why we try solving the ODE with known methods rather than doing silly things like setting \frac{dr}{dt}=0.

Uh, this is worse than I thought it would be. Of course \frac{dr}{dt}=0 is not a solution for the ODE and this is 100% true. But you're missing the fact that there are other equations (basically conditions) involved so that differential equation is not alone to be given a condition (not a solution as you may wrongly have reckoned it to be), i.e. \frac{dr}{dt}=0, and this is by considering the motion to be momentarily or instantaneously at rest along the geodesic that particle follows. This has nothing wrong with it and therefore \frac{d^2r}{dt^2}=-B is very correct with the condition given.

Even worse is the naive attempt by certain members of this forum at calculating \frac{d^2r}{dt^2} by inserting \frac{dr}{dt}=0 in the above ODE and declaring that \frac{d^2r}{dt^2}=-B.

So if you don't have any idea as to what probably the whole "physical picture" is, please do not throw stones across the threads.

AB
 
Last edited by a moderator:
  • #99
Altabeh said:
Physics is not just "mathematical skills"; it requires you to have a skill of painting a physial picture first (the thing that you really sound to be unfamiliar with) and then coloring it with the mathematical hues.

If you don't know how differential equations describe the equations of motion, that's ok.
If you want to learn, then it is not necessary to add your name to the list of trollers, wait a little for espen180 to use the tools I gave him and you'll learn.
Uh, this is worse than I thought it would be. Of course \frac{dr}{dt}=0 is not a solution for the ODE and this is 100% true. But you're missing the fact

No, I'm not missing anything, I am just pointing out that several of you are blissfully basking in the same elementary mistake. Instead of trolling, can you try deriving the equation of motion? It is really simple, you know.

i.e. \frac{dr}{dt}=0, and this is by considering the motion to be momentarily or instantaneously at rest along the geodesic that particle follows. This has nothing wrong with it and therefore \frac{d^2r}{dt^2}=-B is very correct with the condition given.

It is not necessary to resort to your hacks about "momentary" and "instantaneous" motion. If you knew how, you could have derived the general equation of motion, applicable for any t. How about you tried that instead on spending so much energy in ranting? Feel free to use the hints that I gave out in this thread.
 
Last edited:
  • #100
starthaus said:
If you don't know how differential equations describe the equations of motion, that's ok.
If you want to learn, then it is not necessary to add your name to the list of trollers, wait a little for espen180 to use the tools I gave him and you'll learn.




No, I'm not missing anything, I am just pointing out that several of you are blissfully basking in the same elementary mistake. Instead of trolling, can you try deriving the equation of motion? It is really simple, you know.



It is not necessary to resort to your hacks about "momentary" and "instantaneous" motion. If you knew how, you could have derived the general equation of motion, applicable for any t. How about you tried that instead on spending so much energy in ranting? Feel free to use the hints that I gave out in this thread.
Forum rules can be found here: https://www.physicsforums.com/showthread.php?t=5374.

Your repeated violations are too obvious and numerous to point out, and continue in this post.
 
  • #101
kev said:
You are making the assumption that d\tau=dt when the particle is stationary in the metric but it easy to show that is not correct.

True.

Start with the full Schwarzschild metric:

c^2 d\tau^2 = (1-r_s/r)c^2 dt^2 - \frac{1}{(1-r_s/r)} dr^2 - r^2 d\theta^2 -r^2 \sin^2 (\theta) d\phi^2

True, espen180 already knows that.
Impose dr=d\phi=d\theta=0 to obtain:
Wrong.
espen180 has switched to radial motion, so your condition dr=0 is clearly incorrect. The rest of the derivation fails because your initial premise is false.
Once again, the correct condition is d\phi=d\theta=0 (but espen180 already knows that). The correct metric is c^2 d\tau^2 = (1-r_s/r)c^2 dt^2 - \frac{1}{(1-r_s/r)} dr^2

Starting from the correct metric you can derive the general equation of motion. Once you get the equation of motion you can get the correct expressions for acceleration, speed, whatever you need. Don't hack in dr=0. No need for that.
 
Last edited:
  • #102
starthaus said:
Impose dr=d\phi=d\theta=0 to obtain:
Wrong.
espen180 has switched to radial motion, so your condition dr=0 is clearly incorrect.
That makes no logical sense. How does "espen180 has switched to radial motion" result in "your condition dr=0 is clearly incorrect", given that the latter obviously wasn't referring to radial motion?

On what basis is that condition wrong when referring to a stationary particle, as in kev's post you quoted?
 
  • #103
Al68 said:
That makes no logical sense. How does "espen180 has switched to radial motion" result in "your condition dr=0 is clearly incorrect", given that the latter obviously wasn't referring to radial motion?

On what basis is that condition wrong when referring to a stationary particle, as in kev's post you quoted?

All espen180's posts after post 56 refer to radial motion.Since kev attempts to answer post 56, his attempt to set dr=0 in the metric is clearly wrong since r is variable.

(espen180 makes a lesser mistake than kev by attempting to set \frac{dr}{d\tau}=0)
 
Last edited:
  • #104
starthaus said:
All espen180's posts after post 56 refer to radial motion.Since kev attempts to answer post 56, his attempt to set dr=0 is clearly wrong.
So kev is wrong that dr=0 for a stationary particle because espen180's post refers to radial motion? Still makes no logical sense.
 
  • #105
That post was about the acceleration of a particle immediately after being dropped from rest. The word "rest" should justify the dr=0 in the geodesic equations.
 
  • #106
espen180 said:
That post was about the acceleration of a particle immediately after being dropped from rest. The word "rest" should justify the dr=0 in the geodesic equations.

Not if you are unable to derive the equations of motion.I don't know why you persist in hacking the metric. I thought that your whole thread is about deriving the equations of motion, no? I gave you the model of how this is done for circular orbits, the same exact model applies to radial motion. Once you get the equation of motion you can get the correct expressions for velocities, acceleration, whatever you may need.
 
Last edited:
  • #107
Al68 said:
So kev is wrong that dr=0 for a stationary particle because espen180's post refers to radial motion? Still makes no logical sense.

I spent quite a few posts explaining to you why you don't hack the metric (nor the equation of motion) by hacking in dr=0 (or \frac{dr}{d\tau}=0) by hand. Do you think that you could spend your time in a more useful way helping espen180 get the correct equation of motion from the correct (unhacked) metric? The benefit would be obtaining the correct equation of motion, velocity and acceleration.
 
Last edited:
  • #108
espen180 said:
I'm trying to use the tensor formulation of GR to calculate the velocity of a particle in a circular orbit around a black hole.

Here is the work I have done so far.

What concerns me is that I end up getting zero velocity when applying the metric to the differential equations I get from the geodesic equation. I wonder if I have made a miscalculation, but I am unable to find any, so maybe there is a misunderstanding on my part.

Any help is appreciated.

Post 53 answers your question. Now, that you see the general solution, how would you apply it to radial (instead of circular) motion?
 
  • #109
espen180 said:
Okay, if I want to find an expression for the acceleration of a particle without assuming dr/ds=0 initially, I have to solve

\frac{\text{d}^2t}{\text{d}\tau^2}+\left(1-\frac{r_s}{r}\right)^{-1}\frac{r_s}{r^2}\frac{\text{d}t}{\text{d}\tau}\frac{\text{d}r}{\text{d}\tau}=0

\frac{\text{d}^2r}{\text{d}\tau^2}+c^2\left(\frac{r_s}{2r^2}-\frac{r_s^2}{2r^3}\right)\left(\frac{\text{d}t}{\text{d}\tau}\right)^2-\left(\frac{r_s}{2\left(1-\frac{r_s}{r}\right)r^2}\right)\left(\frac{\text{d}r}{\text{d}\tau}\right)^2=0

So I figure the first step is to substitute the dt/dτ. The problem is that

\text{d}t=\text{d}\tau \sqrt{\left(1-\frac{r_s}{r}\right)\left(1-\left(\frac{\frac{\text{d}r}{\text{d}t}}{c^2}\right)^2\right)}

which, if correct, would mean that I get a mix of derivatives of r wrt. t and τ. So I have no idea where to start, or if the equations are analytically solvable.

Hi espen,

The equations are analytically solvable and all the information you need is already contained in the Schwarzschild metric for radial motion only:

c^2 d\tau^2 = (1-r_s/r) c^2 dt^2 - \frac{1}{(1-r_s/r)} dr^2

Divide through by (1-r_s/r) d\tau^2 and rearrange:

c^2\left(\frac{dt}{d\tau}\right)^2 = \frac{c^2}{(1-r_s/r)} + \frac{1}{(1-r_s/r)^2} \left(\frac{dr}{d\tau} \right)^2

Substitute this expression for (c dt/d\tau)^2 into your second equation and you obtain:

\frac{\text{d}^2r}{\text{d}\tau^2}+\left(\frac{r_s}{2r^2}-\frac{r_s^2}{2r^3}\right)\left( \frac{c^2}{(1-r_s/r)} + \frac{1}{(1-r_s/r)^2} \left(\frac{dr}{d\tau} \right)^2 \right)-\left(\frac{r_s}{2\left(1-\frac{r_s}{r}\right)r^2}\right)\left(\frac{\text{d}r}{\text{d}\tau}\right)^2=0

Simplify the above:

\frac{\text{d}^2r}{\text{d}\tau^2} = -\frac{c^2 r_s}{2r^2} = -\frac{GM}{r^2}

This result means that the acceleration using these coordinates is independent of the falling velocity dr/ds because we have made no assumption of dr/ds=0. It also means that the acceleration is non-zero for any value dr/ds and so it is proved that the claim by Starthaus that the acceleration is zero when dr/ds=0, is false.

Now you may have heard that acceleration IS dependent on velocity and also that the time dilation factor is the product of time dilation due to gravitational potential and time dilation due to velocity. The above calculations seem to contradict this, but it all makes sense when you are careful about who makes the measurements. The following is an attempt to clarify various coordinate measurements and define a notation convention:

Schwarzschild coordinate measurement:

These are measurements made by a stationary observer at infinity and use the notation r, t, dr and dt.

Local measurements:

These are measurements made by a stationary observer at r of a free falling particle passing r and are indicated by the use of primed variables, e.g. r', t', dr' and dt'.

Co-moving measurements:

These are measurements made by an free-falling observer that is co-moving with respect to the free falling particle and local to the particle. In these coordinates the, free falling particle's velocity is always zero and the acceleration of the particle is always zero.
These measurements are indicated by the zero substript as r_o, t_o, dr_o and dt_o.
In the Schwarzschild metric dt_o means the same as the infinitesimal proper time interval d\tau or ds/c.

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

Using these definitions and considering a free falling particle with radial motion only, the Schwarzschild coordinate acceleration according to an observer at infinity is:

\frac{d^2 r}{dt^2} = -\frac{GM}{r^2}\left((1-r_s/r)-\frac{3(dr/dt)^2}{(1-r_s/r)}\right)

(Note that this measurement of acceleration is velocity dependent.)

This equation can be obtained from the general equation I derived earlier in post #48:
kev said:
... the general coordinate radial acceleration of a freefalling particle in the metric:

a = \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)
by setting the angular velocity constant H_c to zero and substituting the full form of K_c back in.

The acceleration of the free falling particle according to a local observer that is stationary at r is:

\frac{d^2r '}{dt' ^2} = - \frac{GM}{r^2}\left(\frac{1-(dr '/dt ')^2}{\sqrt{1-2GM/r}}\right)

The derivation of the above equation can found in post #https://www.physicsforums.com/showpost.php?p=2747788&postcount=345".

The acceleration of the free falling particle, using a mixture of distance measured by the Schwarzschild observer at infinity and time as measured by a co-free-falling observer is:

\frac{d^2r}{dt_o^2} = -\frac{GM}{r^2}

This is the expression I derived at the start of this post from your acceleration equation which is in turn obtained from the derivation in your document http://sites.google.com/site/espen180files/Schwartzschild.pdf?attredirects=0&d=1 (using equation (9))

The acceleration of the free falling particle, using distance and time measured by a co-free-falling observer is:

\frac{d^2r_o}{dt_o^2} = 0

This I suspect is something like the acceleration derived by Starhaus, but he is using co-moving in an inconsistent way (the observer is not co-rotating in his orbital calculations) and he using the symbol dr to mean distance measured by a free falling observer when everyone else is using dr to mean distance measured by the Schewarzschild observer. This is why Starthaus said the velocity dr/dt of a stationary particle at r is not zero, because by his definition of dr/dt being measured by a free falling observer, the stationary particle does not appear to be stationary. Starthaus never defines his variables in physical terms, never checks the physical implications or conclusions of his equations and never checks if the symbols he is using have a different physical meaning to the symbols being used by everyone else. He just blindly applies calculus to symbols (defined differently to everyone else) and when his results do not look the same as everyone elses (not surprising really) he declares everyone else to be wrong. (For someone who is supposed to be good at calculus, his basic algebra is surprisingly shaky too.)
 
Last edited by a moderator:
  • #110
kev said:
Setting dr/dt equal to zero, does not by itself imply d^2r/dt^2 must also be zero. If you think this is mathematically impossible, then you need to refresh your basic math skills as well as your basic physics knowledge.
starthaus said:
What in \frac{d^2r}{dt^2}=\frac{d}{dt}(\frac{dr}{dt}) do you still struggle with?

In a different thread https://www.physicsforums.com/showpost.php?p=2737116&postcount=215 you said:

\frac{d^2r}{dt^2}=-\frac{m}{r^2}(1-2m/r)(3\frac{1-2m/r}{1-2m/R}-2})

Now R is the apogee height or release point and when R=r the velocity dr/dt is zero. Setting R=r gives the acceleration at apogee as:

\frac{d^2r}{dt^2}= -\frac{m}{r^2}(1-2m/r)

This is non-zero, so you are contradicting yourself when you claim that dr/dt=0 means that d^2r/dt^2 must also be zero.
 
  • #111
espen180 said:
So I figure the first step is to substitute the dt/dτ. The problem is that

\text{d}t=\text{d}\tau \sqrt{\left(1-\frac{r_s}{r}\right)\left(1-\left(\frac{\frac{\text{d}r}{\text{d}t}}{c^2}\right)^2\right)}

which, if correct, would mean that I get a mix of derivatives of r wrt. t and τ.

The time dilation equation is not correct. Start with the Schwarzschild metric for radial motion only:

c^2 dt_o^2 = (1-r_s/r) c^2 dt^2 - \frac{1}{(1-r_s/r)} dr^2

\left(\frac{dt_o}{dt}\right)^2 = (1-r_s/r) - \frac{(dr/dt)^2}{c^2(1-r_s/r)}

\left(\frac{dt_o}{dt}\right)^2 = (1-r_s/r) \left(1 - \frac{(dr/dt)^2}{c^2(1-r_s/r)^2} \right)

\left(\frac{dt_o}{dt}\right)^2 = (1-r_s/r) \left(1 - \frac{(dr '/dt ')^2}{c^2} \right)

dt_o = dt \sqrt{\left(1-\frac{r_s}{r}\right) \left(1 - \frac{(dr '/dt ')^2}{c^2}\right) }

Of course you now have a mixture of dt, dt_o and dt' but this is not a problem. I demonstrated how to get around this in the earlier post.
 
Last edited:
  • #112
kev said:
Hi espen,

The equations are analytically solvable and all the information you need is already contained in the Schwarzschild metric for radial motion only:

c^2 d\tau^2 = (1-r_s/r) c^2 dt^2 - \frac{1}{(1-r_s/r)} dr^2

Divide through by (1-r_s/r) d\tau^2 and rearrange:

c^2\left(\frac{dt}{d\tau}\right)^2 = \frac{c^2}{(1-r_s/r)} + \frac{1}{(1-r_s/r)^2} \left(\frac{dr}{d\tau} \right)^2

Substitute this expression for (c dt/d\tau)^2 into your second equation and you obtain:

\frac{\text{d}^2r}{\text{d}\tau^2}+\left(\frac{r_s}{2r^2}-\frac{r_s^2}{2r^3}\right)\left( \frac{c^2}{(1-r_s/r)} + \frac{1}{(1-r_s/r)^2} \left(\frac{dr}{d\tau} \right)^2 \right)-\left(\frac{r_s}{2\left(1-\frac{r_s}{r}\right)r^2}\right)\left(\frac{\text{d}r}{\text{d}\tau}\right)^2=0

Simplify the above:

\frac{\text{d}^2r}{\text{d}\tau^2} = -\frac{c^2 r_s}{2r^2} = -\frac{GM}{r^2}

Very good, you finally managed to reproduce the equation of motion from my blog. If you knew how to use the Euler-Lagrange method you could have done it in two lines weeks ago.
This result means that the acceleration using these coordinates is independent of the falling velocity dr/ds because we have made no assumption of dr/ds=0.

Good, so you finally understood that the hack \frac{dr}{ds}=0 is not correct and that you could get the equation of motion without resorting to this hack. This is what I've been telling you repeatedly.
It also means that the acceleration is non-zero for any value dr/ds and so it is proved that the claim by Starthaus that the acceleration is zero when dr/ds=0, is false.

What I've been telling you is something totally different. If you insist on hacking the metric by putting \frac{dr}{ds}=0 by hand as you've been doing, then, by virtue of elementary calculus, you'd get \frac{d^2r}{ds^2}=0.

This I suspect is something like the acceleration derived by Starhaus, but he is using co-moving in an inconsistent way (the observer is not co-rotating in his orbital calculations)

You mean that you don't understand the general solution at post 56? Why don'y you say so? I can explain it to you.
and he using the symbol dr to mean distance measured by a free falling observer when everyone else is using dr to mean distance measured by the Schewarzschild observer.

Umm, no. The definition is given by the Schwarzschild metric itself.

This is why Starthaus said the velocity dr/dt of a stationary particle at r is not zero,
No, what I said is that you can't hack the metric by putting \frac{dr}{ds}=0 by hand as you have been doing. That's all.
 
  • #113
Your post #53 is very suspect. Let's go through it.

starthaus said:
You are repeating the same errors , you just made L=1=constant, when you differentiate a constant, you get ...zero.

You have added new errors as well. If you want to obtain the lagrangian, then you shout divide by ds, not by dt.

If you work out the equations of motion for a photon (which by definition always has ds=0) then when you divide both sides by ds^2 you end up with:

ds^2=\alpha dt^2-\frac{1}{\alpha}dr^2-r^2d\phi^2

\frac{0}{0} =\alpha \frac{dt^2}{0}-\frac{1}{\alpha}\frac{dr^2}{0}-r^2d\frac{\phi^2}{0}

which is undefined on the left and a bunch of infinite terms on the right. This is a very shaky foundation for a rigorous derivation.

starthaus said:
If you do this, you get the correct Lagrangian:

L=\alpha (dt/ds)^2-\frac{1}{\alpha}(dr/ds)^2-r^2(d\phi/ds)^2

Once you get the Lagrangian, you can get one of the Euler-Lagrange equations:

r^2d\phi/ds=H_c (compare against your incorrect expression).

The other Euler-Lagrange equation is:

\alpha (dt/ds)=K (compare to your other incorrect expression).
I defined:

\frac{\delta{L}}{\delta(ds/dt)} = \frac{1}{\alpha} \frac{ds}{dt} = K_c

where K_c is the constant with respect to cordinate time. It is not incorrect just because it is not the same as the constant with respect to proper time. I never intended K = K_c. They are two different things. In fact K = 1/K_c.

H and H_c are two different things too.

starthaus said:
I know that I have written this stuff for you before.

There is a third Euler-Lagrange equation, it is dependent on the other two, so it does not contain any extra information per se:

-d/ds(1/\alpha*2dr/ds)-(2m/r^2(dt/ds)^2-d/dr(1/\alpha)(dr/ds)^2-2r(d\phi/ds)^2)=0

If you make r=R in the above, this means the cancellation of the terms in dr/ds and if you giving you

(d\phi/ds)^2 = \frac{m}{R^3}

Let's break down your steps:

-d/ds(1/\alpha*2dr/ds)-(2m/r^2(dt/ds)^2-d/dr(1/\alpha)(dr/ds)^2-2r(d\phi/ds)^2)=0

You then introduce R without defining it, but in our previous discussions we have used R to mean the radial coordinate when the particle is at its apogee or perigee and dr/ds=0. You then claim that when dr/ds=0 that the terms containing dr/ds cancel out (this is not always true so it is a false assumption) so that you obtain:

0-(2m/r^2(dt/ds)^2-0-2r(d\phi/ds)^2)=0

which solves to:

(d\phi/ds)^2 = \frac{m}{R^3} (dt/ds)^2

and not as you claim:

(d\phi/ds)^2 = \frac{m}{R^3}

You can add this to the list of basic algebra blunders you have already made in this thread.
 
Last edited:
  • #114
kev said:
In a different thread https://www.physicsforums.com/showpost.php?p=2737116&postcount=215 you said:

\frac{d^2r}{dt^2}=-\frac{m}{r^2}(1-2m/r)(3\frac{1-2m/r}{1-2m/R}-2})

Now R is the apogee height or release point and when R=r the velocity dr/dt is zero. Setting R=r gives the acceleration at apogee as:

\frac{d^2r}{dt^2}= -\frac{m}{r^2}(1-2m/r)

Yes, I solved this in my blog weeks ago. Did you just get around to reading it?


This is non-zero, so you are contradicting yourself when you claim that dr/dt=0 means that d^2r/dt^2 must also be zero.

Well, what I've been telling you is something entirely different :
-that you can't hack the metric by putting \frac{dr}{dt}=0 by hand as you've been doing. Once you do this hack, it has the consequence that \frac{d^2r}{dt^2}=0 curtesy of basic calculus. This is not about the relationship between instantaneous speed and instantaneous acceleration but between the function \frac{dr}{dt} and the function \frac{d^2r}{dt^2}. Do you understand the difference?
 
  • #115
kev said:
Your post #53 is very suspect. Let's go through it.



If you work out the equations of motion for a photon (which by definition always has ds=0) then when you divide both sides by ds^2 you end up with:

ds^2=\alpha dt^2-\frac{1}{\alpha}dr^2-r^2d\phi^2

\frac{0}{0} =\alpha \frac{dt^2}{0}-\frac{1}{\alpha}\frac{dr^2}{0}-r^2d\frac{\phi^2}{0}

which is undefined on the left and a bunch of infinite terms on the right. This is a very shaky foundation for a rigorous derivation.

LOL. You really need to take a class in calculus. If you want the lightlike metric, you only need to set ds=0. If you manage to do this correctly, you will be rewarded by getting:

\alpha dt^2-dr^2/\alpha=0

I snipped the rest of your nonsense.
 
Last edited:
  • #116
starthaus said:
Very good, you finally managed to reproduce the equation of motion from my blog. If you knew how to use the Euler-Lagrange method you could have done it in two lines weeks ago.
I have known all along that the equations given by Espen in his document are correct. You have been telling him he has been doing it wrong and he has lost confidence in his own calculations and abilities and that is why I have gone to the trouble to explain again that there is nothing wrong with his latest document.
starthaus said:
Good, so you finally understood that the hack \frac{dr}{ds}=0 is not correct and that you could get the equation of motion without resorting to this hack. This is what I've been telling you repeatedly.
It is not incorrect. It is just a special case at apogee and the conclusions drawn from it are only valid at apogee.
starthaus said:
You mean that you don't understand the general solution at post 56? Why don'y you say so? I can explain it to you.
If I could be bothered I could show you that your equation leads to the acceleration being zero for all r. This contradicts the claim you are now making that the acceleration is GM/r^2. It is obvious from this result that the acceleration is non-zero even when the dr/dt=0 because the acceleration in this form is independent of dr/dt.
 
  • #117
Hi Espen,

using software like Maxima you can reduce errors ( or at least check your results ). I've listed a Maxima script that does something similar to your calculation. It might be of interest.

Code:
/*  
Schwarzschild : orbital motion

save as Scwarz-orbit.mac

Load from the file menu - > File|Batch file
*/
kill(all);
load(ctensor);
/* set some flags */
cframe_flag: false;
ratchristof: true;
ratriemann : true;
ratfac : true;
ctrgsimp: true;
/* define the dimension */
dim: 4;
/* list the coordinates */
ct_coords: [t,r,theta,phi];
/* set up the metric */
/* assign to lg a matrix of zeros  ':' means 'assign' */
lg:zeromatrix(4,4);
/* now add the Schwarzschild coefficients */
lg[1,1]:-c^2*(1-2*m/r);
lg[2,2]:1/(1-2*m/r);
lg[3,3]:r^2;
lg[4,4]:r^2*sin(theta)^2;
/* make the inverse matrix */
ug:invert(lg);
/* get Christoffels */
christof(mcs);
cgeodesic(true);

/** for orbital motion set r,phi derivs to 0 **/
geod[1]:subst(Tdotdot, 'diff(t,s,2), geod[1])$
geod[1]:subst(0, 'diff(r,s,1), geod[1])$

geod[2]:subst(0, 'diff(r,s,2), geod[2])$
geod[2]:subst(0, 'diff(r,s,1), geod[2])$
geod[2]:subst(Thetadot, 'diff(theta,s,1), geod[2])$
geod[2]:subst(0, 'diff(phi,s,1), geod[2])$
geod[2]:subst(tdot, 'diff(t,s,1), geod[2])$

geod[3]:subst(0, 'diff(r,s,2), geod[3])$
geod[3]:subst(0, 'diff(r,s,1), geod[3])$
geod[3]:subst(Thetadot, 'diff(theta,s,1), geod[3])$
geod[3]:subst(0, 'diff(phi,s,1), geod[3])$
geod[3]:subst(tdot, 'diff(t,s,1), geod[3])$

geod[4]:subst(0, 'diff(r,s,2), geod[4])$
geod[4]:subst(0, 'diff(r,s,1), geod[4])$
geod[4]:subst(Thetadot, 'diff(theta,s,1), geod[4])$
geod[4]:subst(0, 'diff(phi,s,1), geod[4])$
geod[4]:subst(tdot, 'diff(t,s,1), geod[4])$

/* these are all = 0 giving 4 diff. equations */
geod[1];
geod[2];
geod[3];
geod[4];

solve([geod[2] ], [Thetadot]);
 
  • #118
kev said:
I have known all along that the equations given by Espen in his document are correct.

..up to the point where you introduce the hack dr=0

If I could be bothered I could show you that your equation leeds to the acceleration being zero for all r.

You really need to take the time and to understand how the results were derived. It will require that you learn calculus, the theory of differential equations and the lagrangian approach.
 
  • #119
#109-#111

This is great! I'll take some time to study these posts closely.

Mentz114 said:
Hi Espen,

using software like Maxima you can reduce errors ( or at least check your results ). I've listed a Maxima script that does something similar to your calculation. It might be of interest.

Thank you very much! I will check it out.

starthaus said:
..up to the point where you introduce the hack dr=0

You really need to take the time and to understand how the results were derived. It will require that you learn calculus, the theory of differential equations and the lagrangian approach.

I'm not using the Lagrangian approach. In post #1 I said I am using the geodesic equation. I don't understand where your expression comes from, so I don't like using it.
 
  • #120
espen180 said:
I'm not using the Lagrangian approach. In post #1 I said I am using the geodesic equation. I don't understand where your expression comes from, so I don't like using it.

The two approaches produce the same exact equations. Whether you use geodesic equations or the Euler-Lagrange equations, they are the same exact equations. Same number of equations, same exact form. The Lagrangian method produces the equations quicker and with less chances of making mistakes.
 

Similar threads

  • · Replies 13 ·
Replies
13
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 9 ·
Replies
9
Views
2K
  • · Replies 62 ·
3
Replies
62
Views
6K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 13 ·
Replies
13
Views
5K
  • · Replies 28 ·
Replies
28
Views
2K
  • · Replies 20 ·
Replies
20
Views
2K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 1 ·
Replies
1
Views
1K