# a) Newtonian Gravity alileo studied bodies falling to Earth under gravity and concluded that all bodies fell with the same acceleration independent of size and material. Tycho Brahe made extensive astronomical observations which led Kepler to formulate his three famous laws of planetary motion relative to the Sun. All of these observations were of relative motion but the mass of one body was, in each case, much greater than that of the other. These led Newton to propose his theory of gravity using the concept of force and yielding an equation which gives the acceleration of a body relative to the centre of mass. He could just as well have presented it in the form ( ) 2 / / A B B A A B r m m G a + ? = (1) without invoking the concept of force and only requiring one definition of mass. This means that the principal of equivalence does not appear. That is, the acceleration of body B relative to A, in the radial direction, is proportional to the sum of their masses and inversely proportional to the square of their separation. G is the gravitational constant. # b) Gravitomagnetics It is now proposed that equation (1) be extended to include the relative velocity. The axioms are. a) It is assumed that in mass-free space light travels in straight lines. This defines a non-rotating frame of reference. b) Because all motion is relative there are no other restrictions on the frame of reference. c) Gravity propagates at the same speed as light. d) Mass, or rest mass, is simply the quantity of matter and is regarded as constant. It could be a count of the number of basic particles. The initial proposed equation is based on comparisons with electromagnetics. This equation gives results which agree with the measured results of the precession of the perihelion of Mercury and with the deflection of light grazing the Sun. Also it gives the correct definition for the Schwarzschild Radius. However, it suggests that the speed of light is constant. As a result it does not predict the Shapiro Time Delay. An extra term is then added which gives agreement with the time delay and also generates the accepted value for the Last Stable Orbit. See equations (2a), (3a), (4a) and (5a). The proposed equation is ( ) where a = acceleration of body B relative to body A, v = the relative velocity, r = the separation and ?? ?? = the unit vector from body A to body B. Also c = speed of light, K = G (m A + m B ) and v r is the radial component of velocity. Note that G is a constant which could be incorporated into the definition of the quantity of matter. These equations reduce to equation (1) when v << c. 2 12 2 2 ? 2 r V N e e a a a ? ? ? ? ? ? + ? = + = c v r K r K (4) where ? is the angle between the velocity and the radius. It should be noted that the velocities of the individual bodies do not appear in these equations, only the relative velocity. For two isolated bodies the relative motion is the only measureable value. A convenient definition of force is ( ) 2 1 2 2 2 2 2 r r e v v e a P × × + ? ? ? ? ? ? ? ? ? ? = = c r m Gm c v r m Gm B A B A µ (5) where ) /( B A B A m m m m + = µ , the reduced mass. By definition of the centre of mass (or the centre of momentum) the total momentum is zero with reference to the centre of mass. It is now proposed that the motion of the centre of mass of two bodies is not affected by collision. From this it follows that for a group of particles the motion of the centre of mass is unaffected by internal impacts. The relative acceleration is only radial when the relative velocity is either radial or tangential. In general the moment of momentum can be shown to be a function of the relative position. So, for an elliptic orbit it remains within bounds. # General inferences from equation (2). Reverts to Newtonian form when v << c. The second term of (2) is normal to the velocity. Moment of velocity (or moment of momentum per total quantity of matter) is shown to be a function of r. The equivalence of inertial mass to gravitational mass does not arise. # c) Modified Equations An extra term is added in the direction of the relative velocity. This will affect the speed of light but not its deflection. As the term is a function of c 4 it only has a very small effect on the motion of large bodies in Solar orbits. The new equation is 2 3 2 2 2 2 t e e a a a a 2 r V N ? ? ? ? ? ? + ? ? ? ? ? ? + ? = + + = c v c v r K c v r K r K r T ? (4a) Where v is the relative velocity and v r is the radial component. Also t = v/|v| is the vector in the direction of the velocity. The angle between the velocity and the radius is ?. For the third term the sign of the acceleration depends only on the sign of the radial velocity. K = G(m A + m B ) and c is the speed of light in a gravity free vacuum. This equation will be considered to be the basic for Post Newtonian Gravity. Justification will come from agreement with verified experimental data. Figure 1 Equation (4a) may be re-written as ( ) t e v v e a r r 3 2 2 2 2 2 2 2 2 1 ? ? ? ? ? ? + × × + ? ? ? ? ? ? ? ? ? ? = c v c v r K c r K c v r K r (2a) or 2 2 1 3 2 2 2 2 2 2 t v e a ? ? ? ? ? ? + + ? ? ? ? ? ? ? ? + ? = c v c v r K c r Kv c v r K r r r (3a) From which it is seen that the additional term is negligible when (v/c) 4 is small compared to unity Again, noting that v / v t = means that (3a) may be written as v e a r ? ? ? ? ? ? ? ? ? ? ? ? ? ? + + ? ? ? ? ? ? ? ? + ? = 2 2 2 2 2 2 1 2 1 c v c r Kv c v r K r . or v e a r Q c r Kv c v r K r 2 2 2 2 2 2 1 + ? ? ? ? ? ? ? ? + ? = . where ? ? ? ? ? ? ? ? ? ? ? ? ? ? + = 2 1 c v Q And ( ) 2 2 1 3 2 2 2 2 2 2 2 t e v v e a P r r ? ? ? ? ? ? + × × + ? ? ? ? ? ? ? ? ? ? = = c v c v c r m Gm c r m Gm c v r m Gm r B A B A B A µ (5a) where ) /( B A B A m m m m + = µ , the reduced mass. Equations (2 -5) will account for the precession of the perihelion of Mercury and will predict the observed value for the deflection of light grazing the Sun. These results were heralded as confirmation of Einstein's General Theory of Relativity. They also give the accepted value for the Schwarzschild Radius. Equations (2a -5a) will generate the same results as mentioned above but will also give the accepted value for the Shapiro Time Delay because the speed of light is now affected by gravity. Light is slowed down by gravity which is the opposite of mass particles. The value generated for the Last Stable Orbit is the accepted value. # d) Lorentz force We shall look at the standard theory of electromagnetism, but this is only to obtain some guidance as to the possible form of a gravitomagnetic theory. The electromagnetic notation is based on reference [7] for SI units. The force on a charge q 2 in an electromagnetic field is B v E F × + = 2 2 2 2 q q (1.4.1) B is the magnetic field and E is the electrostatic field and v 2 is the velocity of the charge. This equation defines the Lorentz force. The magnetic field, due to a length of conductor dl 1 carrying a current i 1 is e r v a V r a N ? 2? A B a T Global Journal of Researches in Engineering (A ) Volume Xx XII Issue I V ersion I 1 / 2 1 1 2 1 / 2 4 e l B × = d i r o ? µ . (1.4.2) but i d q 1 1 1 1 l v = so, for point charges, we can write for the force on charge 2 due to charge 1 is F e v v e 2 1 0 1 2 2 1 2 2 1 1 2 2 1 2 2 1 2 1 1 4 4 / / / / / ( ) = + × × ? ? µ ? q q r q q r o (see also page 256 of ref. [7]) More detailed analysis of the electrodynamics shows that the right hand side of the above equations should be multiplied by ?, but is usually omitted for small velocities. The speed of light c o o = 1 ? µ so we have ? ? ? ? ? ? ? ? ? ? ? ? ? ? × × ? ? ? ? ? ? + = 1 / 2 1 2 1 / 2 2 1 / 2 2 1 1 / 2 4 e v v e F c c r q q o ?? ? (1.4.3) For mass elements we assume that ? ? ? ? ? ? ? ? ? ? ? ? ? ? × × ? ? ? ? ? ? + ? = 1 / 2 1 2 1 / 2 2 1 / 2 2 1 1 / 2 e v v e F c c r m Gm ? (1.4.4) Here c, the speed of gravity waves, is assumed to be the same as that of light. The form of this equation will be taken as a guide; the problem is to choose the most appropriate values for the velocities. The assumptions that are made in the following development are: a) That the field generated by body (1) depends on the velocity of body (1) relative to body (2) i.e. (v 1/2 ) and b) The force on body (2) depends on the velocity of body (2) relative to the field i.e. (v 2/1 ). Let v 2/1 = v, therefore v 1/2 = -v. It must be emphasised that force and force fields are inventions solely for the purpose of visualisation. Substituting these values into equation (1.4.4) gives ?? 2/1 = ? ???? ?? 1 ?? 2 ?? 2/1 2 ??? ??/?? ? ? ?? ?? ? × ? ?? ?? × ?? ??/?? ?? = ??? 1/2 (1.4.5) For the Newtonian case the acceleration of body (2) relative to body ( 1) is 1 1 ) ( / 2 1 2 1 / 2 1 / 1 1 / 2 2 1 / 2 1 2 µ 1 2 1 2 1 2 F F F F F r r a = ? ? ? ? ? ? ? ? + = ? ? ? ? ? ? ? ? + = ? ? = ? = m m m m m m m m ? ? ? ? where 2 1 2 1 m m m m + = µ is known as the reduced mass. For the relativistic case we shall define force as ( ) ( ) 1 1 where 2 2 c v c dt d ? = ? ? ? ? ? ? ? ? ? + = = ? µ ? ? ?µ v a v a v F 2 3 Using this definition of force we obtain, using equation (1.4.5), ( ) ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? × × ? ? ? ? ? ? ? + ? = ? + 1 / 2 1 / 2 2 1 / 2 2 1 2 3 ) ( e v v e v a v a c c r m m G c ? ? ? or ( ) ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? + + ? = ? + 1 / 2 2 1 / 2 1 / 2 2 1 / 2 2 1 2 ) ( e v v e e v a v a 2 c c c v r m m G c ? (1.4.6) Where m 1 and m 2 , are the respective rest masses. Change the suffix of the unit vector to r, let K = G(m 1 + m 2 ) and let r = |r 2/1 | be the separation. Equation (1.4.6) now becomes ( ) v e v a v a 2 2 2 2 2 2 2 1 c r Kv c v r K c r r + ? ? ? ? ? ? + ? = ? + ? (1.4.7) Note that this equation is symmetric with regard to the two bodies. The scalar product of equation (1.4.7) with v yields 2 2 2 2 1 ) ( r Kv c v r ? = ? ? ? ? ? ? + ? ? v a By definition, ( ) ? 2 2 2 1 1 = ? / / v c therefore © 2022 Global Journals Global Journal of Researches in Engineering (A ) Volume Xx XII Issue I V ersion I 2 ) ( r Kv r ? = ? a v ? (1.4.8) Substituting equation (1.4.8) into equation (1.4.7) gives ( ) 2 1 2 2 2 2 2 r r e v v e a × × + ? ? ? ? ? ? ? ? ? ? = c r K c v r K as (2) 2 1 2 2 2 2 2 v e a c r Kv c v r K r r + ? ? ? ? ? ? ? ? + ? = as (3) # II. Application to Two Mass Problem a) Polar Coordinates The following development is based on the conventional treatment of the two body gravitational problem. For the dynamics of bodies in Solar orbits the modified equations are not required. Here, r is the separation and e r is the unit vector in the direction of body 2 as seen from body 1. ? is the orientation of the unit vector with respect to the 'fixed' stars and e ? is the unit vector normal to e r in the plane of the motion. Now ? ? ? ? e e a ) 2 ( ) ( 2 ? ? ? ? ? ? ? r r r r r + + ? = and ( ) ? ? e e v ? r r r + = Equation (3a) can now be expressed in component form Q c r Kv c v r K r r r 2 2 2 2 2 2 2 1 + ? ? ? ? ? ? ? ? ? ? ? ? ? ? + ? = ? ? ? ? ? (6) ( ) Q c r v Kv r dt d r r r ? ? ? d du h dt d d du u r ? = ? = 2 1 ? , ? ? d du h d u d h u r ? ? ? ? ? = 2 2 2 2 , and h u d du h u r h r r h ? ? ? ? ? 2 2 3 2 3 2 2 + = + ? = ? ? Equations (6, 7) may now be written 2 2 2 2 2 2 2 2 2 ? ? ? ? ? ? ? + = + + ? ? ? d du c K c Ku h K u h u h d du d u d ?(8) and 2 2 3 2 ? d du h u c K h u ? = ? (9) Since ? ? h dh d dh d u h = = ? ? ? 2 combining with equation (9) gives,(10) Integrating equation ( 10) leads to 2 0 / ) ( 2 0 c u u K e h h ? ? = (11) Therefore, for small variations ) / ) ( 4 1 ( 2 0 2 0 2 c u u K h h ? ? ? (12) where the suffix 0 refers, in this case, to the position ? = ?/2 measured from the periapsis. Substituting in equation ( 8) for h, using equation ( 12), we obtain ) ( 4 2 2 2 0 0 2 2 0 2 2 ? ? ? ? ? ? ? ? + ? ? ? ? ? ? + ? + = + u d du h u u K c K h K u d u d ? ? (13) # b) Precession of the Periapsis Equation ( 4) is very much easier to apply. This equation is equally applicable to the prediction of satellite trajectories. Because in these aces the relative speeds are not close to the speed of light. The equation which was developed in reference [19] for calculating the precession of the perihelion of Mercury per orbit is ( ) ( ) 2 2 2 1 1 6 e a c m m G P ? + = ? ? (14) where a is the semi-major axis and e is the eccentricity. This generates 42.89 arcsec/century. For the binary pulsar PSR 1913+16, which was discovered by Hulse and Taylor in 1974, (see reference [22]), the accepted data is that the masses of the two stars are 1.441 and 1.387 times the mass of the Sun, the semi-major axis is 1,950,100 km, the eccentricity is 0.617131 and the orbital period is 7.751939106 hr. Using equation ( 14) we obtain the result 4.22 deg/yr, which is in agreement with the measured value and that predicted by General Relativity. The orbital decay, or inward spiralling, of binary pulsars is said to be simply due to energy loss caused by gravitational wave emission. This may be the case but energy loss alone will not account for the phenomenon. The loss of mass # Global # c) Moment of Momentum If the additional term is negligible then it can be shown that the moment of momentum is ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? = 2 1 2 1 2 1 1 2 exp r r c K h h which depends on separation but is constant when c = infinity. # d) Schwarzschild Radius For a constant radius v r = 0 and v = v ? so equation ( 3) or (3a) becomes ? ? ? = ? ? ? ? ? ? ? ? + ? so if v ? = c then g r c M G c K r = + = = 2 2 ) 0 ( 2 2 which is known as the Schwarzschild Radius. # e) Last Stable Orbit Numerical integration of equation (3a) shows that the Last Stable Orbit occurs when the radius of the orbit is 3 times the Schwarzschild Radius, which is the accepted result based on General Relativity. If equation ( 3) is used then a value of 2.62 r g may be calculated algebraically. However if Q is not unity, as shown in equation (3a), then equation (12), with Q included is, ) ) ( 2 1 ( 0 2 0 2 Q u u r h h g ? ? ? (12a) where the suffix 0 refers to circular motion when u = u o Substituting in equation ( 8) for h, using equation (12a), equation ( 13) becomes ( ) 1 2 ) ( 4 2 2 2 0 0 2 2 0 2 2 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? + ? ? ? ? ? ? + ? + = + Q d du u d du h Q u u K c K h K u d u d ? ? ? (13a) If ? + = o u u then, for small variations ( ) ( ) ( ) 0 1 2 1 2 = + + ? + ? ? Q u r Q u r o g o g ? ? (13b) For circular motion it can be shown that g o g o r u r u c v ? = ? ? ? ? ? ? 2 2 so ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? + = g o g o r u r u Q 2 1 For a stable near circular orbit then ( ) ( ) ( ) 0 1 2 1 2 > + + ? Q u r Q u r o g o g so when the factor of ? in equation (13b) is zero, algebraic manipulation of (13b) gives r o /r g = 3, which is the accepted value. It also gives a value of 0.5. # f) Deflection of Light In equation (3a) terms 2 and 3 are parallel to the velocity so the component normal to the velocity is With K =M sun G and R s being the radius of the Sun the deflection is 1.75 arcsec. This value agrees with the measured value and with General Relativity. This confirms the assumption that the deflection of light grazing the Sun is small. For small variation of the speed of light assume that v = c . Also, for small deflections the scalar product of e r and n can be seen from Figure (2) to be R s / r. dt d v c v r K ? = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? + ? = ? n e n a r 2 2 1 r R r K s 2 2 ? = ? n a now # Therefore Integration leads to [ ] ? ? ? ? ? ? ? ? ? = i i rc K c 2 , c r K c c i i 2 0 + = ? where c is now the speed of light where r tends to infinity. Or ? ? ? ? ? ? ? ? ? = 2 2 1 rc GM c c i (16) Consider the case of light grazing the Sun at a radius R s and calculate the journey time. As an approximation assume the path to be a straight line. Then, since dt dx c i = ( ) 2 / 2 1 rc GM dx dt c ? = where 2 2 s R x r + = also as 2 2 rc GM << 1 dx R x c GM dt c s ? ? ? ? ? ? ? ? ? ? + + = 2 2 2 / 2 1 ? ? ? ? ? ? ? ? ? ? ? + + + = x R x c GM x s 2 2 2 # ln 2 Between the limits x = 0 to x and t = 0 to t we have the total time ( ) speed change is so small the additional time is less than 1%. ( ) s s R c GM x R x c GM c x t 2 ln 2 2 2 ln 2 # h) Gravity and the refraction of light The form of the extended equation is based on the known observations or deductions. The extra term is required to be in the direction of the relative velocity. Also, because the speed of light is at its maximum then passing though empty space it must reduce when moving into a gravitational field so assume that it depends on the magnitude of the radial speed. From this it follows that the acceleration will be repulsive. Also, for circular orbits the speed remains constant, which is true when v r = 0. The additional term could be ?? ?? = ?? ?? ?? 2 ? ?? ?? ?? ? ?? ?? , where ? is a constant depending on application. When applied to the Last Stable Orbit with ?? = 2 ? ?? ?? ? 2 the value agrees with the generally accepted value. The constant seems reasonable, so ?? ?? = 2?? ?? 2 ? ?? ?? ? 2 ? ?? ?? ?? ? ?? ?? This is the form shown in equation (2a) which leads to predicting the speed of light due to gravity. Consider the case when speed is very close to the speed of light in a gravity free vacuum. # Equation (2a) gives the acceleration parallel to the velocity ?? ? ?? = ? 2?? ?? 2 ?? ?? ?? ? From equation (3a) the acceleration normal to the velocity ?? ? ?? = ? 2?? ?? 2 ?? ?? ? ?? With ? being the angle between the radius from the gravity source and the velocity ?? ?? = ?? cos ? and ?? ð??"ð??" ? ?? = ? sin ? therefore ?? ?? = ?? ? ?? = ? 2?? ?? 2 ? cos ?and?? ?? = ?? ? ?? = ? 2?? ?? 2 ? sin ? So ?? ?? = ?? ?? tan ? (a) The refractive index n is defined as speed of light in a gravity free vacuum/ speed of light in a transparent media. For light passing through different media Snell's Law states that ?? 1 ?? 2 = ??/?? 1 ??/?? 2 = sin ? 2 sin ? 1 = ?? 2 ?? 1 , or ?? 1 / sin ? 1 = ?? 2 / sin ? 2 = ???????????????? # Thus ?? ?? = ?? ?? tan ? as equation (a). This is applicable to the passage of light through the Earth's atmosphere. Hence light passing through a strong gravity field will be affected in the same way as light passing through the atmosphere. This result gives more confirmation of the applicability of the basic equation of the paper. # III. Gravitomagnetics Applied to Rotating Bodies a) Basic Equations When equation ( 2) is applied to two body systems the equation generated is identical to the de Sitter form and agrees with the measurement of precession of the perihelion of Mercury and of the Binary Pulsar PSR 1913+16. The equation is equally applicable if one spherical body is large and non-rotating. Note that the additional term, which is a function of c 3 , is negligible for Solar dynamics. ( ) 2 1 2 2 2 2 2 r r e v v e a × × + ? ? ? ? ? ? ? ? ? ? = c r K c v r K (2rpt) (18a) where c is the speed of light, a is relative acceleration, v is relative velocity and r is relative position. Also e r is the unit vector from body A to body B. where G is the gravitational constant The calculations are made easier for multi-body systems by the use of a defined force as shown in equation ( 5). ( ) 2 1 2 2 2 2 2 r r e v v e a P × × + ? ? ? ? ? ? ? ? ? ? = = c r m Gm c v r m Gm B A B A µ (5rpt) (18b) µ is the reduced mass ( ) B A B A m m m m + / b) Gravity Probe B Gravity Probe B is the study of the precession of a gyroscope in a polar orbit about the Earth. The spin axis of the gyroscopes is perpendicular to the spin axis of the Earth. After over 3 decades of study and design at Stanford University the final results of the experiment were published in 2011 [27]. K = G(m A + m B ) When the new theory is applied to Gravity Probe B the following equations are derived algebraically using equation (5). The modified equation is not required because the relative velocities are not close to the speed of light. # c) Precession of the Periapsis of a small body orbiting a large rotating mass This problem is similar to the discussion of the precession of the perihelion of Mercury except that now the rotation of the Earth is taken into account. For Mercury the rotation of the Sun has negligible effect. The LAGEOS satellites yield results for the so called frame dragging, or Lense Thirring effect, which results from the rotation of the Earth. The Earth is regarded as a uniform spherical body which can be regarded as a set of uniform spherical shells. The sphere, of mass M and moment of inertia I, rotates about the Z axis at a constant angular speed ?. Consider a test body in orbit around the Earth performing an elliptical orbit where e is the eccentricity and a is the semi-major axis and a period of T. The plane has an inclination (inc) relative to the equatorial XY plane of the Earth. Again, based on equation (5).the rate of precession of the periapsis, as seen from the plane of the orbit, in radians per orbit, is ( ) ( ) T inc e a c GI e a c GM P ) cos( 1 2 1 6 2 / 3 2 3 2 2 2 ? ? ? ? = ? ? ?(19) The first term, the de Sitter precession, has been derived algebraically from equation (2). It agrees exactly with the generally accepted form and agrees with the measured results for the precession of the perihelion of Mercury and for the binary pulsar PSR 1913+16. However, the second term, the Lense-Thirring term, justified by numerical integration, is only half of the generally accepted value. # d) Anti-Gravity Weight loss of a rotating ring. Consider a ring rotating about a horizontal axis above a large body, such as the Earth. Using equation (3a) evaluate the component of the acceleration in the radial direction. 2 2 1 3 2 2 2 2 2 2 t v e a ? ? ? ? ? ? + + ? ? ? ? ? ? ? ? + ? = c v c v r K c r Kv c v r K r r r (3a) Noting that v / v t = means that (3a) may be written as v e a r ? ? ? ? ? ? ? ? ? ? ? ? ? ? + + ? ? ? ? ? ? ? ? + ? = 2 2 2 2 2 2 1 2 1 c v c r Kv c v r K r or ?? ??/?? 2 = ? ?1 + ?? 2 ?? 2 ? ?? ð??"ð??" + ?1 + ?? 2 ?? 2 ? 2?? ?? ?? ?? ?? Also as K is G times the total mass therefore as the mass of the ring is very small compared to the large body its mass is negligible. So, K = G (mass of the large body). Now ?? ?? = ?????????? also ?? ? ?? ð??"ð??" = ??????????, where ? is the angle between the velocity and the radius from the large mass. ?? ? ?? ð??"ð??" ??/?? 2 = ?1 + ?? 2 ?? 2 ? ??1 + 2 ? ?? ?? ? ?? ?????? 2 ??? This equation is the radial acceleration divided by the magnitude of the Newtonian gravity acceleration. If v/c is negligible then the right hand side is -1, as Global Journal of Researches in Engineering (A ) Volume Xx XII Issue I V ersion I expected. This equation is applicable to a point mass so for the whole ring we need to integrate for point mass md/2? divided by total mass of the ring m. ? ?1 + ?? 2 ?? 2 ? ??1 + 2 ? ?? ?? ? ?? ?????? 2 ??? ???? ?? ???? ???? Note that the integral of ?????? 2 ?? between 0 and 2? equals ?. The second term in the square bracket is always positive therefore this term of the acceleration will be radially outwards for all parts of the ring. Integration for all parts of the ring leads to ????? ð??"ð??" ??/?? 2 = (?1 + ?? 4 ), where ?? = ?? ?? ? . The value runs from -1, which is the Newtonian, to 0, which means weightless. When the spinning speed is small compared that of light then the chance of a measurable change in gravitational acceleration, or weight, is very unlikely since it depends of the forth power of v/c. negative gravity ever existed then start with a pool of particles with mixed gravity and electrical charge. Electrical charge is stronger than gravity so the positive electrical charged particles will attract the negative ones and will have zero electrical charge. The particles with gravitational positive charge will repel the negatively charges ones so the groups will separate. # IV. Discussion Equation (4a) is easier to apply than the theory of General Relativity (GR) and therefore leaves less room for misinterpretation. That force is a secondary quantity was strongly advocated by H. R. Hertz who regarded force as "a sleeping partner". Force is to dynamics as money is to commerce. Once force has been demoted to a defined quantity then force fields and inertia are also defined quantities, similarly for work and energy. Equation ( 2) is loosely modelled on the Lorentz force but this relationship is for guidance only in the same way that Maxwell used a mechanical model to form his equations. However, he abandoned the reference in his final paper on the subject once he had established that his equations predicted the then known observations. As shown above, when the new approach is applied to two body systems it agrees with the well verified observations of the precession of the perihelion of Mercury, deflection of light passing the Sun and the definition of the Schwarzschild Radius. All agree with the results obtained from the General Theory of Relativity. The third term in equation (4a) was added as it agrees with the measurements of the Shapiro Time Delay and generates a value equal to the accepted value for the Last Stable Orbit. The Gravity Probe B experiment testing the precession of gyroscopes in Earth orbit displays two equations, one for the geodesic term and one for the frame-dragging effect. The geodesic term does not involve the rotation of the Earth but the frame-dragging term does. The same form of equations have been generated algebraically using equation (5). The framedragging term is half of the published value, however, the geodesic term is about two thirds of the published value. The de Sitter effect agrees with the accepted results of analysis whether algebraically or by numerical integration for two body systems or large non-rotating bodies. This is true whether using equation (2) or equation (5). However, for the Lense-Thirring terms there is an unresolved factor which affects the periapsis precession. The published nodal precession test on the Earth satellites LAGEOS I & II, see reference [28], appear to agree with the accepted theory. The inclination of the satellites is approximately 90 0 +/-20 0 . The reason for this is that the accepted Lense-Thirring term does not depend on the inclination but all other effects do and therefore can be cancelled out. See also references [10] and [29]. The gravitational effect on the speed of light is still discussed but apart from the Shapiro Time Delay the effect is negligible when dealing with the motion of bodies. The decrease of the speed of light grazing the Sun is only 4 parts per million. Gravitational Redshift it is sometimes regarded as a proof of GR, however, it can be derived from other fundamental theories. As shown, light passing through a gravitational field refracts in accordance with Snell's Law. It has proved to be impossible, so far, to find any modification to equation ( 4) such that it gives the generally accepted value for the Lense-Thirring effect without changing the de Sitter effect applications. The de Sitter results have been obtained by several observations but the Lense-Thirring effect is very small compared to other effects. In the LAGEOS experiments for the precession of the periapsis the Lense-Thirring effect is less than 1% of the de Sitter effect, which makes it more difficult to evaluate. The GP-B test results have recently been published, reference [27]. There are four gyroscopes, two of which have original framedragging results which are close to that predicted by the new theory. The geodesic results are, on average, close to those of the accepted value. Nevertheless, over a one month period two of the gyroscopes precess at a rate close to the new theory predictions. It is widely stated that the inward spiralling of a binary star system is due to gravitational radiation. The loss of energy alone is not the cause of this effect. Energy loss can be related to outward spiralling, as is the case for the Earth Moon system. However, radiation pressure could be the cause. When general relativity is applied to multiple body systems several authors have produced slightly different results. Some results even do not return to the © 2022 Global JournalsGravitomagnetics a Simpler Approach Applied to Dynamics within the Solar System Year 2022 © 2022 Global JournalsGravitomagnetics a Simpler Approach Applied to Dynamics within the Solar System Year 2022 © 2022 Global JournalsGravitomagnetics a Simpler Approach Applied to Dynamics within the Solar System Newtonian form when the velocities are zero but only if the speed of light is taken to be infinite. This new approach does not undermine the General Theory of Relativity but because it is a simpler method it leaves less room for misinterpretation. Many of the extensions of GR are very complex mathematically, making errors more likely. * Philosophiae naturalis principia mathematica INewton 1687 London * The Meaning of Relativity AEinstein 1967 Chapman and Hall 6th ed. * A History of Mechanics RDugas 1988 Dover * The Principles of Mechanics Presented in a New Form HHertz 1899 11 MacMillan * JAlembert Traité De Dynamique 1758 Paris * Theory and experiment in gravitational physics Cambridge 1993 333 * HCOhanian Electrodynamics 1988 256 Allyn & Bacon * HGoldstein Mechanics 1980 Addison Wesley 2 * Stanford Relativity Gyro Experiment CWEveritt Francis 1986 * WheelerCiufolini Gravitation and Inertia. Princeton Series in Physics 1995 * Gravitational radiation RainerWeiss Rev. Mod. Phys 71 2 1999 * Applications of classical physics KipThorne Caltech. Ch 136 2003 Ph * ThorneMisner WheelerGravitation 1973 Freeman * HRHarrison TNettleton Advanced Engineering Dynamics. Arnold 1997 Wiley * Space-based Gravity Tests CMWill Nature 347 * MSLongair Theoretical Concepts in Physics Cambridge 1984 * A Gravitational and Electromagnetic Analogy. The Electrician OHeaviside 31 1893 * KRSymon Mechanics 1971 Addison Wesley 3 rd ed. * HRHarrison Gravity Galileo to Einstein and Back. (Revised) U.P 2016 * The Search for Gravity Waves PC WDavies 1980 Cambridge U. P * Space-based Gravity Tests CMWill Nature 347 * FowlerTaylor Mcculloch Binary Pulsar PSR1913+16. Nature 277 8 th Feb. 1979 * Six High-Precision Transits of OGLE-TR-113b ERAdams Astrophysical Journal 721 2010 Oct 1 * Instant Expert, General Relativity PedroFerreira New Scientist 3 July 2010 * MathewsPeters Physical Review 131 1963 * TaylorWeisberg arXiv:astro-ph/0407149v1 2004 * Gravity Probe B: Final Results CW FEveritt Phys.Rev. Lett 106 221101 May 2011 * Detection of the Lense-Thirring Effect due to Earth/s Spin ICiufolini arXiv:gr-qc/9704065v1 23 Apr. 1997 * Phenomenology of the Lense-Thirring effecting the Solar System LIorio arXiv:1009.3225v2[gr-qc]30 Sept 2010 * Gravitomagnetics Simpler Approach Applied to Rotating Bodies. H. Ron Harrison June 2011