Jump to content

Analytical equation for Space-Time geodesics and relativistic orbit equation


pengkuan

Recommended Posts

Orbit equation and orbital precession

General Relativity explains gravity as Space-Time curvature and orbits of space objects as Space-Time geodesics. This concept is extremely hard to understand and geodesics hard to compute. However, if we can express Space-Time geodesics in analytical form, that is, as orbit equation like Newtonian orbit equation for planets, relativistic gravity will be much simpler and most people can understand.

From gravitational and inertial accelerations, I have derived the orbit equation for relativistic gravity, see equation (57), of which I will explain the derivation below. I will refer this equation as the relativistic orbit equation. Albert Einstein had correctly predicted the orbital precession of planet Mercury which had definitively validated General Relativity. The orbital precession that this orbit equation gives is identical to the one Albert Einstein had given [1][2], see equation (65).

Because the relativistic orbit equation gives the same result than Space-Time geodesics, it is an analytical expression of Space-Time geodesics which is so simple that personal computer can use to compute geodesics rather than big or super computer. The direct derivation from gravitational and inertial accelerations is a new insight into General Relativity. 

 Relativistic dynamics

  • Velocity in local frame

Figure 1 shows a body of mass m orbiting around a body of mass M which sits at the origin of the polar coordinate system. The position of m with respect to M is specified by the radial position vector r, of which the magnitude is r and the polar angle is q

Let frame_m be an inertial frame of reference that instantaneously moves with m. Newton’s laws apply in frame_m because m’s velocity is 0, so the inertial acceleration vector of m am equals the gravitational acceleration ag, see (1).  

Let frame_l be the local frame of reference in which M is stationary, m move at the velocity vl and the acceleration of m is a l. As frame_m moves with m, its velocity is vl in frame_l. We transform al into am using the transformation of accelerations between relatively moving frames, see equation (18) in «Relativistic kinematics and gravity»[3][4]. This transformation is in (2) where we have replaced a1 with al, a2 with am and u with vl

Equating (1) with (2) we get (3), both sides of which are dotted with the vector vl·dtl in (4) to give  on the left hand side, see (5), and  on the right hand side, see (6). Plugging (5) and (6) into (4) gives (7), which we integrate to give (8), with K being the integration constant. The value of eK is given in (9). Rearranging (8) gives (10), which expresses the local orbital velocity vl due to the gravitational field of M. 

  • Velocity in the frame at infinitely far

The orbit of a planet is the curve drawn by its position with respect to a frame in which the gravitational effect of M is 0. This frame is infinitely far from M and is labeled frame_¥. In this frame the time is t¥, the radial position vector is r and the velocity vector v¥ equals the derivative of r with respect to t¥, see (11)

The velocity vector vl in frame_l equals the derivative of the same r with respect to the time in frame_l, tl. We transform vl into the product of v¥ and  in (12).  is the time dilation factor of the Schwarzschild metric around M, see (13). By plugging (13) into (12) we obtain the relation between vl and v¥ in (14). Substituting (14) for vl in (10) gives (15), which expresses v¥ with the gravitational field of M. In the case of planets,  is very small and the term  is simplified into the Taylor’s series in (16)

Notice that the right hand side of (15) is the product of two parentheses, the first represents gravitational acceleration, see (10), the second represents the time dilation of Schwarzschild metric, see (13)

Substituting (16) for  in (15) gives (17), from which we drop out the third order term  to get  in (18), where u¥ is the velocity vector of m in frame_¥ correct to the 2nd order of  

  • Differential equation for orbit

The general expression of velocity in polar coordinates is given in (19) and the velocity squared u¥2 in (20). For expressing u¥2 with h the constant of Kepler’s second law given in (21) [5][6], we transform  in (22), which is plugged with (21) into (20) to give (23), which expresses u¥2 in terms of    instead of r

We divide (23) with c2 to express  in terms of , see (24). For transforming (18), we plug (24) and  into (18) to get (25) and (26), which is the differential equation for orbit in frame_¥. Solving (26) will give the relativistic orbit equation. 

 

Figures and equations are in the pdf here:Analytical equation for Space-Time geodesics and relativistic orbit equation.pdf

Link to comment
Share on other sites

Create an account or sign in to comment

You need to be a member in order to leave a comment

Create an account

Sign up for a new account in our community. It's easy!

Register a new account

Sign in

Already have an account? Sign in here.

Sign In Now
×
×
  • Create New...

Important Information

We have placed cookies on your device to help make this website better. You can adjust your cookie settings, otherwise we'll assume you're okay to continue.