Gravitational Field of a Stationary Homogeneous Oblate Spheroidal Massive Body

Transcription

1 International Journal of Theoretical Mathematical Physics 08, 8(3): 6-70 DOI: 0.593/j.ijtmp Gravitational Field of a Stationary Homogeneous Oblate Spheroidal Massive Body Nura Yakubu,*, S. X. K. Huwusu, A. Adamu Department of Physics, University of Maiduguri, Maiduguri, Nigeria Department of Theoretical Physics, National Mathematical Center Abuja, Nigeria Abstract The application of Newton s laws of motion in gravitational fields has been restricted to the fields of perfectly spherical bodies. In this paper the gravitational fields due to a Spheroidal body is derived to open the way for the extension of Newtonian mechanics from the fields of spherical bodies to these Spheroidal bodies. The results showed that Newton s gravitational field equations for a Spheroidal body are linear separable can be solved in terms of the well known special functions of mathematical physics, the Legendre functions. We also obtained gravitational intensity field g of the body its corresponding unit vectors in components forms. We further obtained general acceleration vector a its associated acceleration vector components forms. We finally obtained the components of the equation of motion for a particle of non zero rest mass in the gravitational field exterior to an oblate spheroidal massive body. The Newton s universal gravitational field equations for a homogeneous oblate Spheroidal body with the exact complete results are available for mathematical analysis hence physical interpretation experimental investigation for all bodies in the universe. Keywords Gravitational field, Homogeneous, Static, Spheroidal. Introduction In the present paper some theoretical background of gravitational field of Spheroidal bodies is stated. Also in this paper the gravitational scalar potential for both the interior exterior regions are evaluated by means of Legendre differential equations. Then the gravitational intensities of oblate Spheroidal bodies for both the regions are determined. Prior to 950, theoretical study of gravitation was restricted almost exclusively to the fields of massive bodies of perfect spherical geometry. For example, in Newton s theory of universal gravitation (T. U. G.) the motion of particles (such as projectile, satellites penduli) in earth s atmosphere is treated under the assumption that the earth is a perfect sphere. Similarly, in the solar system the motion of bodies (such as planets comets asteroids) is treated under the assumption that the sun is perfect sphere. Also in Einstein theory of gravitation, called General Relativity Theory (G. R. T.) the motion of bodies (such as planet) particles (such as photons) is treated under the * Corresponding author: nuraphy@gmail.com (Nura Yakubu) Published online at Copyright 08 The Author(s). Published by Scientific & Academic Publishing This work is licensed under the Creative Commons Attribution International License (CC BY). assumption that the sun is a perfect sphere (Schwrchiild space time). But it is well known that the only reason for these restrictions is mathematical convenience simplicity. The real fact of nature is that all rotating planets stars galaxies in the universe are Spheroidal is obvious that the Spheroidal geometry will have corresponding consequences effects in the motion of all particles in their gravitational fields. These effects will exist in both Newton s T. U. G Einstein s G. R. T. consequently the way is prepared for the study of motion of all particles in the gravitational fields of Spheroidal bodies. As an example it is now well known that satellite orbits around the earth are not govern by just simple universe distance square gravitational field of perfect spherical geometry. They are also govern by second harmonics (pole of order 3) as well as fourth harmonics (pole of order 5) of gravitational scalar potential due to imperfect spherical geometry. Now for comparison with these approximations the exact analytical gravitational scalar potential for a perfect oblate Spheroidal body is derived (Howusu, 00).. Analysis Consider a homogeneous oblate Spheroidal body of rest mass M 0. Let S be the referenc frame whose origin 0 coincides with centre of the body as shown in figure.

2 International Journal of Theoretical Mathematical Physics 08, 8(3): z y x Figure. Oblate Spheroid Then the oblate Spheroidal coordinates (η, ξ ) are defined in terms of the Cartesian coordinates (x, y, z) by where a is a constant x = a[( + ξ )( η )] / cos φ () y = a[( + ξ )( η )] / sin φ () z = aηξ (3) 0 ξ η (4) 0 π ξ = ξ 0 (5) where ξ 0 = a constant. Now since the body is homogenous, the density of active mass l is given by ρ(r) = ρ 0: ξ ξ 0 0: ξ > ξ 0 (6) where ρ 0 is the constant density of rest mass. As is well known, Newton s gravitational field equation for the gravitational scalar potential Φ due a distribution of mass density ρ is given by Φ g (r) = 4πGρ(r) (7) where G is the universal gravitational constant. It follows from the explicit expression for the Laplacian operator in oblate Spheroidal coordinate that the interior exterior gravitational scalar potentials Φ respectively satisfy the equations (Nura, et. al., 07), 4πGρ 0 = 0 = ( + a (ξ +η ) ξ ξ ) + ( ξ η η ) + ξ η η (+ξ )(η ) ( + a (ξ +η ) ξ ξ ) + ( ξ η η ) + ξ η η (+ξ )(η ) φ φ g (ξ, η, φ) (8) φ φ g (ξ, η, φ) (9) To solve these equations it may be noted that by the symmetry of the distribution of mass about the polar axis the dynamical gravitation a potential will be independent of the azimuthal angle φ. Hence equations (8) (9) reduce to 4πGρ 0 a (ξ + η ) = ξ ( + ξ ) ξ φ g (ξ, η) + η ( η ) η φ g + (ξ, η) 0 = ξ ( + ξ ) ξ φ g (ξ, η) + η ( η ) η φ g + (ξ, η) (0a) (0b)

3 64 Nura Yakubu et al.: Gravitational Field of a Stationary Homogeneous Oblate Spheroidal Massive Body Next, let us seek a separable complementary solution of the interior, equation, Φ g e, as where λ is the separation constant. Consequently, for the choice ξ ( + ξ ) ξ U (ξ) λu (ξ) = 0 () η ( + η ) η V (η) λv (η) = 0 () Equation (0) becomes the Legendre s equation hence has solution of the form λ = l(l + ), l = 0,,, (3) V (η) = P l(η), l = 0,,, (4) Q l (η) where P l Q l are the linearly independent solutions of the Legendre s differential equation. Also by the transformation equation equation () becomes ξ = iτ (5) τ ( + τ ) ξ U (τ) l(l + )U (τ) = 0 (6) where equation (6) is the Legendre s differential equation hence has the solutions of the form Consequently U (ξ) = P l(g), l = 0,,, (7) Q l (iξ) Φ g e (ξ, η) = [A l P l (iξ) + B l Q l (iξ)][c l P l (η) + D l Q l (η)] l=0 where A l, B l, C l D l are arbitrary constants. Next let us seek a particular solution of the interior equation Φ g p Substituting (9) into equation (0) a it follows that Hence the general solution of the interior equation Φ g is given by (8) Φ g p = K(ξ η ) (9) 4πGρ 0 a = 6K (0) K = 3 πgρ 0a () Φ g (ξ, η) = 3 πgρ 0a (ξ η ) + [A l P l (iξ) + B l Q l (iξ)][c l P l (η) + D l Q l (η)] l=0 where A l, B l, C l D l are arbitrary constants P l Q l are the Legendre s function of order. Similarly, since the exterior equation is homogeneous its general solution Φ g + is given by Φ g + (ξ, η) = [A l + P l (iξ) + B l + Q l (iξ)][c l + P l (η) + D l + Q l (η)] l=0 where A + l, B l, C l D l are arbitrary constants. Now since the interior exterior regions both contain the coordinates η = 0 which is a singularity of Q l D + l = 0, D l = 0 l = 0,,, (4) Also since Q l (iξ) is not defined at the center of the body ξ = 0 in the interior region B + l 0, l = 0,,, (5) Also since Q l (iξ) is not defined at ξ = 0 in the ixterior region A + l 0, l = 0,,, (6) Hence from the continuity of the gravitational scalar potentials () (3) across the boundary ξ = ξ 0 the orthogonality of the Legendre function P l it follows that, the resulting equations are A 0 + πgρ 0a ξ = B Q 0 (iξ) (7) 3 as () (3)

4 International Journal of Theoretical Mathematical Physics 08, 8(3): Equations (7) (9) we get A P (iξ 0 ) = B + Q (iξ) (8) A P (iξ 0 ) 4πGρ 0a = B + Q (iξ 0 ) (9) 9 4πGρ 0 a ξ 0 3 = B 0 + d dξ Q 0(iξ)ξ=ξ 0 (30) Thus, solving (7) (3) a simultaneously: A 0 + ξ P l(iξ 0 )ξ=ξ 0 = B l + ξ Q l(iξ 0 )ξ=ξ 0 (3a) B 0 + = 4πGρ 0a ξ 0 3 ξ Q 0(iξ) ξ=ξ 0 (3b) A 0 + = 4πGρ 0a ξ 0 Q 0 (iξ 0 ) 3 ξ Q 0(iξ) ξ=ξ 0 πgρ 0a ξ (3) 3 3 Next, let us consider equation (9) above i.e. A P (iξ 0 ) 4πGρ 0a = B + 9 Q (iξ 0 ) Now differentiate both sides of this equation with respect to ξ, if follows that A 0 = 0 (33) B + = 0 (34) A ξ P (iξ)ξ=ξ 0 = B + ξ Q (iξ)ξ=ξ 0 (35) By multiplying both sides of equation (9) with ξ P (iξ)ξ=ξ 0 it follows that A P (iξ 0 ) P 4πGρ0a ξ (iξ) P = B ξ=ξ 9 ξ 0 (iξ)ξ=ξ Q (iξ 0 ) P (36) ξ 0 (iξ)ξ=ξ 0 This follows by multiplying both side of equation (36) with P (iξ 0 ). Thus where B + = 4πGρ 0 a ξ P (iξ) ξ=ξ 0 9[ ] (37) [ ] = Q (iξ 0 ) ξ P (iξ)ξ=ξ 0 P (iξ 0 ) ξ P (iξ)ξ=ξ 0 (38) 4πGρ 0 a A + ξ Q (iξ) ξ=ξ 0 = (39) 9[ ] Hence the general solution of the interior- exterior equation in equation () (3) reduced to Φ g (ξ, η) = πgρ 0a ξ η + A 3 P (iξ)p (η) + A 0 P 0 (iξ)p 0 (η) (40) Φ + g (ξ, η) = B + 0 Q 0 (iξ)p 0 (η) + B + Q (iξ)p (η) (4) These are the gravitational scalar potentials of the body in terms of rest mass density. It follows that gravitational scalar potential (40) (4) may be expressed in terms of the total rest mass of the body by the relation It may be noted that results from mathematical tables gives (Nura, 0) M 0 = 4 9 a3 ξ 0 πρ 0 (3 + ξ ) (4) Q 0 (t) = l n t+ (43) t

5 66 Nura Yakubu et al.: Gravitational Field of a Stationary Homogeneous Oblate Spheroidal Massive Body Consequently Q 0 (t) = 4 (3t )l n t+ + 3 t (44) t Q 0 (iξ) = i ξ + 3 ξ3 + 5 ξ5 + (45) Q (iξ) = ξ ξ ξ5 + (46) 3. Evaluation of Gravitational Field Intensity for Both the Interior Exterior Scalar Potentials The Del operator is given by equation below Hence putting It follows that = u h u u + v h v v + w h w w (47) u = η (48) v = ξ (49) w = Φ (50) (η, ξ, φ) = η h η η + ξ h ξ ξ + Φ h φ φ Therefore, it follows from definition that the gravitational intensity field g, of the body is given by Now, Putting the general solution g+ in the above equation (5) give the following results (5) g = Φ g = g η η + g ξ ξ + g φ φ (5) g η + = h η η Φ g + (ξ, η) (53) Or Equation (55) can take the form: Since Similarly Or g + η = η / [B a ξ +η η 0 + Q 0 (iξ)p 0 (η) + B + Q (iξ)p (η)] (54) g + η = η / + B0 Q a ξ +η 0 (iξ) P η 0(η) + B + Q (iξ) P η (η) (55) g η + = B + a η ξ +η / Q (iξ) η P (η) (56) P 0 =, (57) g ξ + = h ξ ξ Φ g + (ξ, η) (58) Or g + ξ = ξ / + B0 Q a ξ +η 0 (iξ) P ξ 0(η) + B + Q (iξ) P ξ (η) (59) g φ + = h φ φ Φ g + (ξ, η) (60) g φ + = 0 (6) To obtain a gravitational field intensity for the interior potential Φ g + (ξ, η), we substitute (44) into the equation

6 International Journal of Theoretical Mathematical Physics 08, 8(3): Or Since Similarly hence g + η = η a / 4πGρ 0 a ξ +η 3 g + ξ = η a g + ξ = η a / πgρ 0 a ξ +η 3 g η + = h η η Φ g + (ξ, η) (6) + A 0 P 0 (iξ) η P 0(η) + A P (iξ) η P (η) (63) / πgρ 0 a ξ +η 3 + A P (iξ)p (η) (64) P 0 =, (65) g ξ + = h ξ ξ Φ g + (ξ, η) (66) + A 0 P 0 P 0 (η) ξ (iξ) + A P (η) ξ P (iξ) (67) g φ + = 0 (68) The results (5), (57) (68) (6) are the gravitational intensity field interior to the body while (6), (670 (68) are those exterior to the body. Next, precisely as is well known cases of the spheroidal polar coordinates, the definition of the oblate spheroidal coordinates where The scalar factors are And the corresponding unit vector are given by η = η+ξ I = xı + yȷ + zk (9) x = a( η ) / ( ξ ) / cos φ (0) y = a( η ) / ( ξ ) / sin φ () z = aηξ () η, 0 ξ <, 0 φ π (3) cos φ (η +ξ ) / h η = a η +ξ η / h ξ = a η +ξ ξ / (4) (5) h φ = a[( η )( + ξ )] / (6) ı η+ξ sin φ (η +ξ ) / ȷ + ξξ k (7) (η +ξ ) / ξ = ξη cos φ (η +ξ ) / ı + ηη sin φ (η +ξ ) / ȷ + ηξ k (8) (η +ξ ) / The general position vector r is given by And the general velocity vector U is given by φ = { sin φ ı + cos φ ȷ } (9) r = aηη η + aξ+ξ ξ (0) (η +ξ ) / (η +ξ ) /

7 68 Nura Yakubu et al.: Gravitational Field of a Stationary Homogeneous Oblate Spheroidal Massive Body u = dr dt = aη η +ξ / η + aξ η +ξ / ξ + aφ ( η ) ( + ξ ) / φ () (η ) (η ) For oblate spheroidal coordinates, the general acceleration vector a is given by a = aη η +ξ (η ) + aξ η +ξ (ξ ) + aηη +ξ (η +ξ ) (η ) + aξξ η (η +ξ ) (+ξ ) aηξ η (η +ξ ) (+ξ ) + aξη +ξ (η +ξ ) (η ) aξηξ (η +ξ ) (η + aηφ η +ξ ) (η +ξ ) + η aηη ξ (η +ξ ) (+ξ aξφ η +ξ ) (η +ξ ) + a( η )( + ξ )φ + aξη ξ φ aη+ξ η φ φ (3) Equation (3) above reduced to the form where a η = aη +ξ η + (η ) a ξ = aη +ξ ξ + (ξ ) aη+ξ (η ) (η +ξ ) aξη (η +ξ ) (+ξ ) (+ξ ) η + ξ (η ) ξ () a = a η η + a ξ ξ + a φ φ (4) aηη (η +ξ ) (+ξ ) aξ+ξ (η +ξ ) (η ) ξ + aηη +ξ φ + (η +ξ ) η + aξη +ξ φ + (η +ξ ) aξ (η +ξ ) (η ) aη (η +ξ ) (+ξ ) η ξ (5) η ξ (6) a φ = a( η )( + ξ )φ + aξη ξ φ aη+ξ η φ (7) (+ξ ) (η ) Consequently, it follows from Newton s law of motion in gravitational fields a = g (8) And equation (7) (9) (5) (7) that the components of the equation of motion given by η+ξ η + η + ξη η ξ + ηη ξ + ηη +ξ φ (η )(η +ξ ) (η +ξ ) (η +ξ )(ξ ) (η +ξ ) = B + η a ξ +η Q (iξ) η P (η) (9) = +ξ a (η +ξ ) B 0 + P 0 (η) ξ Q 0(iξ) + B + P (η) ξ Q (iξ) (30) ξη ξ+ξ ξ + ξ η + η ξ + η + ξη +ξ φ (3) (+ξ )(η +ξ ) (η +ξ ) (η +ξ )(η ) These are the complete equations of motion for a particle of non zero rest mass in the gravitational field exterior to an oblate spheroidal massive body. 4. Summary Conclusions In this paper, Newton s Universal Gravitational Potential field equation for a homogeneous oblate Spheroidal body is (η +ξ ) φ = ξ +ξ ξ φ η η η φ (3) formulated solved with the exact complete results given by (40) (4). Consequently, these potentials are now available for application in physics. An application of the results obtained in this work is that the gravitational scalar potential (40) (4) imply the corresponding gravitational intensities (56), (59), (6), (64), (67) (68). A profound philosophical inference is suggested by the results obtained in the theses that Newton s gravitational field equations for a Spheroidal body are linear separable

8 International Journal of Theoretical Mathematical Physics 08, 8(3): hence solvable in terms of the well known special functions of mathematical physics, the Legendre function. This fact suggests that Einstein s geometrical theory of gravitation, GRT, whose gravitational field equation is non linear may not be the most natural generalization of Newton s dynamical theory of gravitation. In the present work, the gravitational field intensity or force fields for both the interior exterior gravitational scalar potentials were derived. These sets of equations pave the way for the equations of motion for a particle of non-zero rest mass using spheroidal coordinates. It is now known that most bodies in the universe are Spheroidal n nature. As an example, it is now well known that satellite orbits round the earth are not governed by just the simple inverse distance squared gravitational fields of perfect spherical geometry. They are also governed by second harmonics (pole of order 3), as well as fourth harmonics (pole of order 5) of gravitational scalar potential due to imperfect spherical geometry. In the first place this work opens the door for the physical interpretation of all the solutions obtained in this work hence experimental investigation in the motion of all bodies in the earth s atmosphere solar system as well as all other gravitating systems in the universe. In the second place the door in henceforth opened for the mathematical study of all the unsolved equations of motion for all particles in all types of motions in all systems in the universe, such as more accurate calculation of (i) Missile satellite space craft trajectory in the earth s atmosphere (ii) motions of moons around their planets, (iii) motion of planets (iv) comets asteroids around the sun (v) motions of stars around their galactic nuclei in the universe. In the third place it is most interesting instructive to note that by the definition of the function Q 0 Q the gravitational intensity fields due to oblate Spheroidal body are all even orders in the universe coordinates. Since this work pave the way for new physical investigations experimental verifications of the motion of bodies in the universe. It may be noted that precisely as in the case of the gravitational field of a stationary homogeneous spherical massive body the door is henceforth opened for the theoretical investigation of the following problems. Pure ξ motion of a particle of non zero rest mass in the gravitational field exterior to a stationary oblate spheroidal massive body. Plane motion of a particle of non zero rest mass in the gravitational field exterior to a stationary homogeneous oblate spheroidal massive body. Pure ξ motion of a photon in the gravitational field interior to a stationary homogeneous oblate spheroidal massive body. Gravitational frequency shift of a photon in the gravitational field exterior to a stationary homogeneous oblate spheroidal massive body. Gravitational time dilation in the gravitational field exterior to a stationary oblate spheroidal massive body. Gravitational time dilation in the gravitational field interior to stationary oblate massive body. Electrogravitational coupling in the gravitational field exterior to a homogeneous oblate spheroidal massive body. Electrogravitational coupling in the gravitational field interior to a homogeneous oblate spheroidal massive body. Gravitational length contraction in the gravitational field exterior to a stationary homogeneous oblate spheroidal massive body. Gravitational length contraction in the gravitational field interior to a stationary homogeneous oblate spheroidal massive body. Radar sounding in the gravitational field exterior to a stationary homogeneous oblate spheroidal massive body. Radar sounding in the gravitational field interior to a stationary homogeneous oblate spheroidal massive body. REFERENCES [] Albert, E. (979), Quanta cosmology: Johnson Corporation, London. [] Ashatague, H. Kshif, A. (994), Curvature Corlineatiolls of some Spherically Symmetry Space-Tme. American Journal of Mathematical Physics 37(7) [3] Atwater, H. E. (974), Introduction to General Relativity. Pergamen Press, New York. [4] Fock, V. (997), the Theory of Space-Time Gravitation. Pergamen Press, New York. [5] Forster, J. Nightingale, J. D. (979), A Short Course in General Relativity, Longman, New York. [6] G. Arfken, Mathematical Methods for Physics. (Academic Press, New York, 968). [7] Howusu, S. X. K. (996), Cosmological Redshift. Experimental Detection of Gravitational Radiation Aperion (Montreal) 3 () [8] Howusu, S. X. K. (00), Natural Philosophy of Gravitation, Jos University Press, Jos Nigeria. [9] Howusu, S. X. K. (003), Gravitational Field of Spherical Bodies Extension of Gravitational Fields of Spherical Bodies, Galilean Electrodynamics. [0] J. L. Anderson. Principles of Relativity (academic Press, New York, 967). [] Mark, K. (986), Gravitational Field Intensity. American Journal of Physics 46() [] N. Yakubu S. X. K. Howusu (0). Gravitational Field of a Stationary Homogeneous Spheroidal Massive Body. Journal of Science Technology Research, Vol. 0, No. 3, pp 0 9.

9 70 Nura Yakubu et al.: Gravitational Field of a Stationary Homogeneous Oblate Spheroidal Massive Body [3] N. Yakubu, S. X. K. Howusu, O. Nwagbara I. Nuhu (07). Solution of Newton s Gravitational Field Equation of a Static Homogeneous Prolate Spheroidal Massive Body. International Journal of Engeneering Science Innovative Technology (IJESIT). Vol. 6, Issue, pp [4] Rosser, W. G. (967), Introductory Relativity. Butter worth, London. [5] Slater, J. C. (933), Introductory Theoretical Physics, Mc Graw Hill Publishing Company Ltd, London. [6] Solelger, I. (98), Fundamentals of Theoretical Physics. Vol. Mir Publishers, Moscow. [7] Weinberge, S. (975), Gravitation Cosmology Principles applications of the General Theory Relativity. (J Willey, New York). [8] William, H. (974), H Book of Mathematical Sciences CRC Press International. Boca Roton. Florida.