Skip to main content

Section 10.2 Gravitational Interactions

Gravitational interactions are one of four fundamental action-at-a-distance interactions which are known to occur between particles. At the particle level, they are described by Newton’s law of gravitation. We are interested in a macroscopic description of gravitational interactions.

Subsection 10.2.1 Newton’s Law of Gravitation

The success of Newton in predicting the motion of a falling body on Earth or the motion of the planets around the Sun hinges on his law of universal attraction which can be stated as follows:

Definition 10.2.1. Law of Gravitation.

Two particles \(P_1\) and \(P_2\) of mass \(m_1\) and \(m_2\text{,}\) respectively, and at a distance \(r\) from each other, attract each other with a force proportional to their mass and inversely proportional to \(r^2\text{.}\) Hence, the force due to \(P_1\) on \(P_2\) can be written as
\begin{equation} \bF_{1 \to 2}^g = - G {m_1 m_2 \over |\br_{P_1 P_2}|^3 } \br_{P_1 P_2}\tag{10.2.1} \end{equation}
where \(G= 6.67 \times 10^{-11} \)N m\(^2\)kg\(^{-2}\) is a universal constant, called constant of gravitation.
Figure 10.2.2.
Note that \(\bF_{2 \to 1}^g = - \bF_{1 \to 2}^g\text{.}\) Every mass particle attracts any other mass particle. However the effect of gravitational forces is only appreciable if at least one of the “particles” has a very large mass. For instance, on Earth’s surface we neglect the mutual attraction of bodies in comparison with Earth’s gravitational pull. This law is applicable for applications ranging from terrestrial mechanics to celestial mechanics.

Subsection 10.2.2 Gravitational Field

Definition 10.2.3. Gravitational Field.

The effect of a discrete or continuous mass distribution \(\Sigma_1\) is equivalent to a gravitational field \(\bG_{\Sigma_1} (P)\) at any location \(P\) of space, such that, a particle of mass \(m\) located at point \(P\) at a particular instant would be subjected to the force
\begin{equation} \bF_{\Sigma_1 \to P}^g = m \bG_{\Sigma_1} (P)\tag{10.2.2} \end{equation}
For a discrete set of particles \(\Sigma_1 = \{ P_1 , P_2 , \cdots P_N \}\) of mass \(m_1\text{,}\) \(m_2\text{,}\) \(\cdots\) , \(m_N\text{,}\) we find
\begin{equation} \bG_{\Sigma_1} (P) = - G \sum_{i=1}^N {m_i \bu_i \over r_i^2} \tag{10.2.3} \end{equation}
with \(r_i = |\br_{P_i P}|\) the distance from \(P_i\) to \(P\text{,}\) and \(\bu_i = \br_{P_i P} /r_i\) the unit vector directed from \(P_i\) to \(P\text{.}\) Hence the gravitational field is obtained by a superposition of individual fields. For a continuous mass system \(\Sigma_1\text{,}\) we obtain
\begin{equation} \bG_{\Sigma_1} (P) = - G \int_{\Sigma_1} {\br_{P_1 P} \over | \br_{P_1 P}|^3} dm(P_1)\tag{10.2.4} \end{equation}
where \(dm (P_1) = \rho (P_1) dV\) is an infinitesimal mass of volume \(dV\) surrounding \(P_1\text{.}\) In general, the calculation of this integral is difficult, unless certain assumptions are made, as will be done in Subsection 10.2.3.
The macroscopic gravitational effect of a material system \(\Sigma_1\) on material system \(\Sigma\) can then be modeled in terms of the gravitational action screw which takes the following form
\begin{equation} \{\cA^g_{\Si_1 \to \Si}\} = \left\{ \begin{array}{c} \bR_{\Si_1 \to \Si}^g \\ \;\\ \bM_{A, \Si_1 \to \Si}^g \end{array} \right\} = \left\{ \begin{array}{c} \int_\Si \bG_{\Sigma_1}(P) \, dm(P) \\ \;\\ \int_\Si \br_{AP}\times \bG_{\Sigma_1} (P) \, dm(P) \end{array} \right\}_A\tag{10.2.5} \end{equation}
For most terrestrial applications, one considers only the effects of Earth’s gravitational forces, that is, one neglects the gravitational effects of all celestial bodies other than that of Earth, because such systems are too far, or their mass is too small relative to that of Earth. Furthermore, for applications which take place in domains of small dimensions compared to that of Earth, one may assume that the gravitational field \(\bG_{\Earth\to \Si} (P)\) to be approximately uniform.
However, in many applications of celestial mechanics, these assumptions are not adequate and one must then resort to use better approximations of Earth’s gravitational field \(\bG_{\Earth\to \Si} (P)\) to obtain useful expressions of \(\{ \cA^g_{\Earth \to \Si}\}\text{.}\)

Subsection 10.2.3 Gravitational Field of a Body with Spherical Symmetry

Most celestial bodies such as Earth are spherical to a first order of approximation. In addition, their mass density can be considered to be only a function of the distance \(r\) from their center. Such bodies are said to possess spherical material symmetry. Then, by symmetry, the resulting gravitational field at any point is only a function of the radial distance measured from this point to the center of the body. More specifically, we can state the following result:

Proof.

The proof is provided by the divergence theorem for the flux of gravitational field \(\bG_{\Sigma_1} (P)\) across a surface \(\partial \Omega\) enclosing a domain \(\Omega\)
\begin{equation*} \int_{\partial\Omega} \bG_{\Sigma_1} \cdot \bn \, dA = \int_\Omega \nabla \cdot \bG_{\Sigma_1} dV = -4\pi G \int_\Omega \rho dV \end{equation*}
where \(\bn\) being the outward unit normal to surface \(\partial \Omega\text{.}\) Recall that \(\nabla \cdot \bG_{\Sigma_1} = - 4\pi G \rho (P)\) where \(\rho\) is the local mass density. Let us choose for \(\partial\Omega\) the sphere of center \(G_1\) and radius \(r\text{.}\) Due to the spherical symmetry we expect \(\bG_{\Sigma_1}\) to be of the form \(\bG_{\Sigma_1}(P) = g(r) \be_r\) for all \(P\) on \(\partial \Omega\text{,}\) where \(\be_r = \bn\) is unit radial vector at \(P\text{.}\) This leads to
\begin{equation*} \int_{\partial\Omega} \bG_{\Sigma_1}(P) \cdot \be_r dA = 4\pi r^2 g(r)\text{.} \end{equation*}
With \(\int_\Omega\rho dV = M_1\text{,}\) we obtain \(g(r)=- G M_1/r^2\text{.}\) Note that the same approach allows for the determination of field \(\bG_{\Sigma_1}\) at a point \(P\) of the interior of \(\Sigma_1\text{.}\)
The Gravitational Field of Earth: We can treat Earth as a body of spherical material symmetry of center \(C\) and radius \(R\text{:}\) indeed, the difference between the equatorial radius \(R_E\) and the polar radius \(R_P\) leads to the value \(1/300\) of the ratio \((R_E - R_P)/R_E\text{.}\) Then, the gravitation field at a point \(P\) located at altitude \(z\) above Earth’s surface (see Figure 10.2.5) is given by (with \(r = R+z\)):
\begin{equation} \bG_{E}(P) = - G \frac{M_E}{(R+z)^2}\be_r \tag{10.2.7} \end{equation}
where \(M_E\) is Earth’s mass, and \(\be_r\) is the unit vector directed from \(C\) to \(P\text{.}\) If \(z \ll R\text{,}\) we can approximate \(1/(R+z)^2\) as \((1- 2z/R)/R^2\) to first order, and we obtain
\begin{equation} \bG_{E}(P) \approx - g_0 (1- \frac{2z}{R} )\be_r \tag{10.2.8} \end{equation}
where
\begin{equation} g_0 = {G M_E \over R^2}\tag{10.2.9} \end{equation}
is the gravitational field at ground level. Hence, it is seen that the magnitude of Earth’s gravitational field decreases with altitude, though the variation is very weak: with \(R = 6,400\) km, we find a decrease of \(1\%\) at elevation \(z = 32\) km. Its direction also varies weakly: two points on the surface of the Earth whose corresponding unit vectors \(\be_r\) make an angle of one degree are separated by a distance of \(110\) km. Hence, we see that, for most terrestrial applications, we may treat the gravitational field as a constant vector in both magnitude and direction:
\begin{equation} \bG_{E}(P) = - g_0 \be_r\tag{10.2.10} \end{equation}
With \(M_E = 6 \times 10^{24}\) kg, we find \(g_0 \approx 9.81\) m/s\(^2\text{.}\) In this case, the gravitational action screw \(\{\cA_{\Earth\to \Si}^g\}\) takes the form
\begin{equation} \{\cA_{\Earth\to \Si}^g \} = \left\{ \begin{array}{c} \int_\Si \bog \, dm(P) \\ \;\\ \int_\Si \br_{AP}\times \bog \, dm(P) \end{array} \right\}_A = \left\{ \begin{array}{ccc} m_\Si \bog \\ \;\\ \bze \end{array} \right\}_G\tag{10.2.11} \end{equation}
where \(m_\Si\) denotes the mass of system \(\Si\text{,}\) and \(G\) denotes the mass center of \(\Si\text{.}\)
Figure 10.2.5.
In the following problem, a correction term is obtained for the expression (10.2.11) of the gravitational action screw due to a body with spherical material symmetry on a “small” body.

Example 10.2.6. (Small body approximation).

Consider the gravitational effect of Earth \(\cE\) (of mass \(M\text{,}\) and mass center \(C\text{.}\) possessing spherical material symmetry on a satellite \(\cS\) (of mass \(m\text{,}\) and mass center \(G\)) considered small in relation to its distance to Earth. See Figure 10.2.7.
Show that the resulting gravitational action screw \(\{\cA^g_{\cE \to \cS}\}\) is given by
\begin{equation*} \{\cA^g_{\cE \to \cS}\} = - {G Mm \over \rho^2} \left\{ \begin{array}{c} \bx -\frac{3}{2m\rho^2} [ (2I_{Gx}-I_{Gy} - I_{Gz})\bx - 2 I_{Gxy} \by - 2I_{Gxz} \bz] \\ \\ {3 \over m\rho} (I_{Gxz} \by - I_{Gxy} \bz) \end{array} \right\}_G \end{equation*}
where unit vector \(\bx\) is defined by \(\br_{CG} = \rho \bx\text{,}\) and where \((\bx, \by, \bz)\) forms a right-handed basis. The scalars \(I_{Gx}= \bx \cdot \cI_G (\bx)\text{,}\) \(I_{Gy} =\by \cdot \cI_G (\by)\text{,}\) and \(I_{Gz}= \bz \cdot \cI_G (\bz)\) are moments of inertia about mass center \(G\text{.}\) The scalars \(I_{Gxy} = \bx \cdot \cI_G (\by)\text{,}\) \(I_{Gxz} = \bx \cdot \cI_G (\bz)\) are the corresponding products of inertia.
Hint.
To derive this result, assume that \(\br \equiv \br_{GP}\) is much smaller in magnitude than \(\rho = |\br_{CG}|\) for any point \(P\) of \(\cS\text{,}\) and expand the expression \((1+ 2 \bx\cdot \br /\rho+ \br^2 /\rho^2)^{-3/2}\) up to second-order in powers of \(|\br| /\rho \ll 1\text{.}\)
Figure 10.2.7.
Solution.
We first find the general expression of the resultant of the action screw by using \(\bG_{\cE} (P) = - G M(\rho\bx + \br)/|\rho\bx + \br|^3\) and denoting \(\br = \br_{GP}\text{:}\)
\begin{align*} \bR^g _{\cE \to \cS} \amp = -G M \int_\cS \frac{\rho\bx + \br}{(\rho^2 + 2 \rho \bx \cdot \br + \br^2)^{3/2}} dm\\ \amp = -\frac{G M}{\rho^2} \int_\cS \frac{\bx + {\br\over \rho}} {(1 + 2 \bx \cdot {\br \over \rho}+ {\br^2\over \rho^2})^{3/2}} dm \end{align*}
Assuming \(|\br /\rho| \ll 1\text{,}\) we then expand \((1 + 2 \bx \cdot {\br \over \rho}+ {\br^2\over \rho^2})^{-3/2}\) to second order:
\begin{equation*} \Big(1 + 2 \bx \cdot {\br \over \rho}+ {\br^2\over \rho^2}\Big)^{-3/2} = 1 - 3 \bx\cdot {\br\over \rho} - \frac{3}{2}{\br^2\over \rho^2} +\frac{15}{2}\big(\bx \cdot {\br\over \rho} \big)^2 + \cdots \end{equation*}
leading to
\begin{equation*} \bR^g _{\cE \to \cS} = -{G mM \over \rho^2}\bx + {3G M \over 2\rho^4} \int_\cS [(\br^2 -5(\bx\cdot\br)^2)\bx +2(\bx\cdot\br) \br] dm +\cdots \end{equation*}
where we have used the fact that \(\int_\cS \br \, dm = \bze\text{.}\) Each of the integrals in the brackets (the second-order terms) can be expressed in terms of the inertia characteristics of \(\cS\) about its mass center: by writing \(\br = x\bx + y\by + z\bz\) we find
\begin{equation*} \int_\cS [(\br^2 -5(\bx\cdot\br)^2)\bx +2(\bx\cdot\br) \br] dm = \int_{\cS} [(y^2+z^2-2x^2)\bx +2xy \by + 2 xz \bz ] dm \end{equation*}
This gives the final expression
\begin{equation*} \bR^g _{\cE \to \cS} = -{G mM \over \rho^2}\bx + {3G M \over 2\rho^4} \left[ (2I_{Gx}-I_{Gy} - I_{Gz})\bx - 2 I_{Gxy} \by - 2I_{Gxz} \bz \right]+ \cdots \end{equation*}
Note that, as expected, the second-order term vanishes for a body \(\cS\) of spherical material symmetry, that is, when \(I_{Gx}= I_{Gy}= I_{Gz}\) and \(I_{Gxy} = I_{Gxz} = 0\text{.}\) To find the moment about mass center \(G\text{,}\) we first write
\begin{equation*} \bM^g_{G, \cE \to \cS} = \int_{\cS} \br\times \bG_{\cE} (P) dm =-G M \int_{\cS} \br\times (\rho\bx + \br) (\rho^2 + 2 \rho \bx \cdot \br + \br^2)^{-3/2} dm \end{equation*}
We then expand to leading order of \(|\br|/\rho\text{:}\)
\begin{equation*} \bM^g_{G, \cE \to \cS} = - {G M \over \rho^2} \int_{\cS} \br \times \bx (1- 3\bx\cdot {\br \over \rho} + \cdots ) dm = {3 G M \over \rho^3} \int_{\cB} \br \times \bx (\bx\cdot \br) dm + \cdots \end{equation*}
With \(\int_{\cS} \br \times \bx (\bx\cdot \br) dm = \int_{\cS} x (-y \bz + z\by)dm = (I_{Gxy} \bz - I_{Gxz} \by)\text{,}\) we find the final results.

Subsection 10.2.4 Weight of a Particle and Gravitational Acceleration

Consider a plumb line experiment in a referential \(\cE\text{,}\) whereby a particle \(P\) of mass \(m\) is suspended by a string attached at a point \(Q\) fixed in \(\cE\text{.}\) See Figure 10.2.9. If \(P\) is at rest in \(\cE\text{,}\) we call vertical in \(\cE\) the direction taken by the plumb line.

Definition 10.2.8. Gravitational Acceleration.

The weight \(\bW\) of particle \(P\) of mass \(m\) in referential \(\cE\) is the force which balances the tension in the string of a plumb line experiment (in the absence of any other external perturbing influence). Then, the gravitational acceleration at the location of \(P\) in \(\cE\) is the vector field \(\bg (P)\) defined by
\begin{equation} \bW = m \bg (P)\tag{10.2.12} \end{equation}
Figure 10.2.9. Plumb line experiment.
Hence, the weight of a particle and the corresponding gravitational acceleration depend on the local value of the gravitational field \(\bG\) created by one or more celestial bodies and on the transport acceleration of referential \(\cE\) relative to a (Newtonian) referential \(\cF\) relative to which it is in motion:
\begin{equation} \bg(P) = \bG (P) - \ba_{P\in\cE /\cF}\tag{10.2.13} \end{equation}
The same particle of mass \(m\) will have a different weight whether the plumb line experiment is done on Earth, on the Moon, in the space station, or in an elevator in acceleration.
In terrestrial applications, \(\bg\) is mostly a function of Earth’s gravity field \(\bG = \bG_E\) (which is mostly a function of elevation) and to a much lesser extent a function of Earth’s rotation (thus dependent of latitude). The latter effect is obtained by evaluating \(\ba_{P\in\cE /\cF}\) relative to the geocentric referential \(\cF\text{.}\) The following problem quantifies the effect of Earth’s rotation.

Example 10.2.10.

Consider the plumb line experiment done in the terrestrial referential \(\cE\) attached to Earth at a colatitude \(\phi\text{.}\) Apply Newton’s second law relative to the geocentric referential \(\cF\) about which a terrestrial referential \(\cE\) rotates with angular velocity \(\bom_{\cE/\cF} = \omega_E \bk\text{,}\) that is, solve equation
\begin{equation*} -m g_0 \be_r + W \be_z = m \ba_{P/\cF} \end{equation*}
to find the weight \(W\) of a particle \(P\) of mass \(m\text{,}\) and the angle \(\theta\) taken by the vertical \(\be_z\text{.}\) Use the expression \(\ba_{P/\cF} = \om_E \bk \times (\om_E \bk \times \br_{CP})\approx \om_E^2 \br_{OQ}\text{.}\) See Figure 10.2.11 for the definitions of unit vectors \(\bk, \be_z, \be_r\text{.}\)
Figure 10.2.11.
Solution.
We treat Earth as a body of spherical material symmetry, and neglecting the effect of altitude (\(z=0\)), we approximate the gravitational field as \(\bG_{E}(P) = - g_0 \be_r\) with \(g_0 = {G M_E / R^2}\text{.}\) We denote by \(\be_z\) the unit vector along the string from \(P\) to the attachment point \(Q\text{,}\) by \(\be_r\) the unit vector directed from Earth’s center \(C\) to \(P\) or to any point “nearby” such as origin \(O\text{.}\) The string tension is given by \(W \be_z\text{,}\) denoting by \(W\) the weight of particle \(P\) in \(\cE\text{.}\) If we account for the rotation of the Earth, the equilibrium of \(P\) is now governed by the equation
\begin{equation*} m \ba_{P / \cF } = -m g_0 \be_r + W \be_z \end{equation*}
where \(\cF\) is the so-called geocentric referential defined by Earth’s axis of rotation \(SN\text{,}\) with two additional axes pointing to “fixed stars”, and about which referential \(\cE\) (Earth) has angular velocity \(\bom_{\cE / \cF}= \om_E \bk\text{,}\) with \(\om = 2 \pi\) \(\rm rad/day\text{.}\) \(\bk\) is the unit vector directed from the south pole \(S\) to the north pole \(N\text{.}\) The change of referential formula for acceleration gives us
\begin{align*} \ba_{P/\cF} \amp = \ba_{P / \cE }+ \ba_{P \in \cE / \cF} + 2 \bom_{\cE / \cF} \times \vel_{P/\cE} = \ba_{P \in \cE / \cF } \approx \om_E \bk \times (\om_E \bk \times \br_{CP}) \\ \amp = -R \om_E^2 \sin\phi \; \bu \end{align*}
since \(P\) is fixed in \(\cE\text{,}\) and the motion of \(\cE\) relative to \(\cF\) is a uniform rotation about axis \((C, \bk)\text{.}\) Unit vector \(\bu\) is defined as the vector perpendicular to \(\bk\) lying in the plane containing axis \((C, \bk)\) and point \(P\text{.}\) Hence, the equilibrium equation of \(P\) becomes
\begin{equation*} -m g_0 \be_r + W \be_z + m R \om_E^2 \sin\phi \; \bu = \bze \end{equation*}
Let us define the ratio \(\ep\) of accelerations by
\begin{equation*} \ep = {R \om_E^2 / g_0} \end{equation*}
With \(R = 6400\) km, \(\om_E = \frac{2\pi}{24\times 3600}\) \(\rm rad/s\text{,}\) and \(g_0 = 9.8\) \(\rm m/s^2\text{,}\) we obtain a very small value \(\ep = 0.0034\text{.}\) Now from the vectorial equation
\begin{equation*} -\be_r + {W \over mg_0} \be_z + \ep \sin\phi \; \bu = \bze \qquad (1) \end{equation*}
we can determine the direction of the vertical \(\be_z\) and the weight \(W\) in referential \(\cE\text{:}\) denote by \(\te\) the angle \((\be_r , \be_z)\text{.}\) To eliminate \(W\) from equation (1), take the cross-product with \(\be_z\) to obtain
\begin{equation*} \tan\te = \ep \frac{\sin\phi\cos\phi}{1 - \ep \sin^2\phi} \end{equation*}
Since \(\ep \ll 1\text{,}\) the angle \(\te\) is small and can be obtained from
\begin{equation*} \te = \ep \sin\phi\cos\phi . \end{equation*}
To obtain the weight \(W\) and the corresponding gravitational acceleration \(g\) in \(\cE\text{,}\) take the scalar product of equation (1) with \(\be_r\)
\begin{equation*} W = mg , \qquad g = g_0 (1 - \ep \sin^2 \phi) . \end{equation*}
The maximum angle \(\te\) leading to the minimum gravitational acceleration is obtained at colatitude \(\phi = \pi /4\text{:}\)
\begin{equation*} \te_{\max{}} = 0.0017, \qquad g_{\min{}} = 0.997 \; g_0 . \end{equation*}
These results show that considering \(\cE\) as a Newtonian referential for a body at rest is a valid approximation: the error incurred by neglecting the Earth’s rotation is very small.
A more detailed description of the gravitational acceleration \(\bg (P)\) at a point \(P\) of the surface of the Earth can be given by including the effect of the gravitational field generated by both the Moon \(\bG_M\) and the Sun \(\bG_S\text{,}\) and by accounting for the transport acceleration relative to a heliocentric referential \(\cF\text{:}\)
\begin{equation*} \bg (P) = \bG_E (P)+ \bG_M (P) + \bG_S (P) - \ba_{P\in\cE /\cF} \end{equation*}
where the transport acceleration \(\ba_{P\in\cE /\cF}\) includes the acceleration \(\ba_{C/\cF}\) of the mass center of the Earth around the Sun, in addition to the term \(-\omega_E^2 \br_{QP}\) accounting for the rotation of the Earth about its axis
\begin{equation*} \ba_{P\in\cE /\cF} = -\omega_E^2 \br_{QP} +\ba_{C/\cF} \end{equation*}
where point \(Q\) is the projection of \(P\) on the Earth’s axis of rotation. To find \(\ba_{C/\cF}\) we apply Newton’s second law for the motion of \(C\) relative to referential \(\cF\)
\begin{equation*} M \ba_{C/\cF} = M ( \bG_M (C) + \bG_S (C) ) \end{equation*}
under the gravitational attraction of the Moon and the Earth.

Definition 10.2.12. Tidal Acceleration.

The terrestrial gravitational acceleration can be expressed in the form
\begin{align*} \bg (P) = \bG_E (P)+ (\bG_M (P)-\bG_M (C)) \amp + (\bG_S (P)- \bG_S (C)) \\ \amp + \omega_E^2 \br_{QP} \end{align*}
where \(G_M\) and \(G_S\) are the gravitational fields generated by the Moon and the Sun, respectively. The differential terms \(\bG_M (P)-\bG_M (C)\) and \(\bG_S (P)- \bG_S (C)\) are called tidal accelerations.

Remark 10.2.13.

The contributions of the tidal accelerations to \(\bg(P)\) are responsible for the tidal phenomena of the Earth’s oceans, and for other celestial phenomena.
The following problem attempts to give an explanation of the ocean’s tides.

Example 10.2.14.

Consider the system Earth-Moon shown in Figure 10.2.15. The plane of the figure is assumed to be the equatorial plane of Earth containing the mass center \(M\) of the Moon. Both bodies are assumed spherical and of homogeneous mass. The lines \(CM\) intersects Earth’s surface at points \(P_1\) and \(P_2\text{.}\) Find the tidal accelerations \(\bG_M (P)-\bG_M (C)\) at these two diametrically opposed points to show that two effects of opposed directions and equal magnitude are obtained, thus explaining the bulging of the oceans at these locations (high tides). Assume that the distance \(D_M = |CM|\) is much larger than Earth radius \(R\text{.}\) How many daily high tides are obtained? Where are low tides obtained?
To assess the strength of the Moon’s tidal acceleration relative to that of the Sun, evaluate the ratio \(|\bG_M (P_1)-\bG_M (C)|/ |\bG_S (P_1)-\bG_S (C)|\text{.}\) Use \(D_M / D_S= 2.5 \times 10^{-3}\) and \(m_M /m_S = 3.65 \times 10^{-8}\text{.}\)
Figure 10.2.15.
Solution.
Denote by \(\bu\) the unit vector \(\br_{CM} /|\br_{CM}|\text{.}\) Assuming that the Moon is spherical and homogeneous, its gravitational field at points \(P_1\text{,}\) \(P_2\) and \(C\) are given by (denoting \(m_M\) the mass of the Moon)
\begin{equation*} \bG_M (P_1) = {G m_M \over (D_M+R)^2} \bu, \quad \bG_M (P_2) = {G m_M \over (D_M-R)^2} \bu, \quad \bG_M (C) = {G m_M \over D_M^2} \bu \end{equation*}
This leads to the expressions
\begin{align*} \bG_M (P_1)-\bG_M (C) \amp = \left({G m_M \over (D_M+R)^2}-{G m_M \over D_M^2} \right)\bu \\ \amp = {G m_M \over D_M^2} \left( (1 - 2{R\over D_M} + \cdots) -1 \right)\bu \\ \amp = -{2G m_M R\over D_M^3} \bu \end{align*}
\begin{align*} \bG_M (P_2)-\bG_M (C) \amp = \left({G m_M \over (D_M-R)^2}-{G m_M \over D_M^2} \right)\bu\\ \amp = {G m_M \over D_M^2} \left( (1 + 2{R\over D_M} + \cdots) -1 \right)\bu \\ \amp = {2G m_M R\over D_M^3} \bu \end{align*}
where we have used the fact that \(R/D_M \ll 1\text{.}\) We have thus obtained two opposed tidal accelerations of equal magnitude at \(P_1\) and \(P_2\text{.}\) This explains the rising of Earth’s ocean surface.
Tidal phenomena vary in time due to the Earth’s rotation about its axis. Assuming the Moon stationary relative to a geocentric referential, a given meridian will reach the position \(P_1 P_2\) twice a day, thus leading to two high tides per day. In reality, the Moon has a period of 27 days around the Earth, leading to high tides every 12h25min. Low tides are obtained at the locations of a meridian which is perpendicular to line \(P_1 P_2\text{.}\)
We can replace the role played by the Moon by that of the Sun:
\begin{equation*} \bG_S (P_1)-\bG_S (C) = -{2G m_S R\over D_S^3} \bu , \qquad \bG_S (P_2)-\bG_S (C) = {2G m_S R\over D_S^3} \bu \end{equation*}
We obtain the ratio of tidal accelerations
\begin{equation*} {|\bG_M (P_1)-\bG_M (C)| \over |\bG_S (P_1)-\bG_S (C)|} = {m_M \over m_S} {D_S^3 \over D_M^3} \approx 2.3 \end{equation*}
This shows that the Moon’s tidal force is stronger than that of the Sun by a factor of 2.3. These two effects are superimposed when the Sun and Moon become aligned with Earth: this happens about twice a month, at the new and full moons, leading to a maximum range between low and high tides.