Comprehensive Study Notes on Universal Gravitation

Newton’s Law of Gravitation

Newton’s law of gravitation was formulated by Isaac Newton in 1665. It states that the force of attraction between two masses, m1m_1 and m2m_2, separated by a distance rr, is directly proportional to the product of their masses and inversely proportional to the square of the distance between them. This is expressed by the formula F=Gm1m2r2F = \frac{Gm_1m_2}{r^2}. In this equation, GG is known as the universal gravitational constant and its value is approximately 6.67×1011Nm26.67 \times 10^{-11}\,Nm^{-2}.

Gravitational Field Intensity

Gravitational field intensity at a specific point is defined as the gravitational force per unit mass placed at that point. For the Earth, the gravitational field intensity is denoted as gg. The mathematical representation for field intensity is I=Fm\vec{I} = \frac{\vec{F}}{m}, where the test mass mm is considered to be exceptionally small. This ensures that the test mass itself does not significantly alter the gravitational field it is measuring.

Gravitational Potential (VgV_g)

Gravitational potential at a particular point is the amount of work done to bring a unit mass from infinity to that point under the influence of the gravitational field produced by a given mass MM. The potential is represented by the formula Vg=GMrV_g = -\frac{GM}{r}. The negative sign indicates that the work is being done by the gravitational field or that the potential represents an attractive system.

Potential and Field for Mass Distributions

For a system of discrete mass distributions, the total gravitational potential is the algebraic sum of individual potentials, given as V=V1+V2+V3+V = V_1 + V_2 + V_3 + \dots, which can be expressed in summation notation as V=i=1NViV = \sum_{i=1}^N V_i. Similarly, the total gravitational field intensity is the vector sum of individual intensities: I=I1+I2+I3++IN\vec{I} = \vec{I}_1 + \vec{I}_2 + \vec{I}_3 + \dots + \vec{I}_N, or I=i=1NIi\vec{I} = \sum_{i=1}^N \vec{I}_i.

In the case of continuous mass distributions, the total potential and field intensity are determined through integration. The total potential is given by V=dVV = \int dV, where (dV)(dV) is the potential due to an elementary mass dMdM. The total field intensity is given by I=dI\vec{I} = \int d\vec{I}, where dId\vec{I} is the field intensity generated by the elementary mass dMdM.

Gravitational Potential Energy (UgU_g)

The gravitational potential energy of a two-mass system is defined as the amount of work done to bring a mass mm from infinity to a point PP under the influence of the gravitational field of a primary mass MM. It is calculated as Ug=GMmrU_g = -\frac{GMm}{r}. It is important to note the relationship between gravitational potential energy and gravitational potential: Ug=Vg×mU_g = V_g \times m.

In general terms, the gravitational potential energy of a system is the work done against the gravitational force in order to assemble the system from its reference configuration. For mass systems, the standard reference configuration is one where there is infinite mutual separation between the masses.

Gravitational Field Intensity of a Ring

For a ring of radius RR and mass MM, the gravitational field intensity at any point PP on the axial line at a distance xx from the centre of the ring is given by Eg=GMx(R2+x2)3/2E_g = \frac{GMx}{(R^2 + x^2)^{3/2}}. This field is always directed towards the centre of the ring. At the centre of the ring, the field intensity is at its minimum, which is Eg=0E_g = 0. The value of EgE_g reaches its maximum when the distance is related to the radius by x=R2x = \frac{R}{\sqrt{2}}.

Relation Between Field and Potential

There is a direct mathematical relationship between the gravitational field intensity and the gravitational potential, expressed as the negative gradient of the potential. Specifically, I=dVdr\vec{I} = -\frac{dV}{dr}. In Cartesian coordinates, this field component is broken down as I=(Vxi^+Vyj^+Vzk^)\vec{I} = -\left(\frac{\partial V}{\partial x}\hat{i} + \frac{\partial V}{\partial y}\hat{j} + \frac{\partial V}{\partial z}\hat{k}\right). Alternatively, the potential difference can be written in integral form as dV=IdrdV = -\vec{I} \cdot d\vec{r}.

Work Done and Gravitational Force

The work done against gravitational force to change the configuration of a system is equal to the difference between the potential energy in the final configuration and the potential energy in the initial configuration. The equation is Wagainst gravitational force=U2U1W_{\text{against gravitational force}} = U_2 - U_1. Conversely, the work done by the gravitational force is the negative of this value: Wby gravitational force=(U2U1)W_{\text{by gravitational force}} = -(U_2 - U_1).

Variation of Gravity with Height (hh)

The acceleration due to gravity changes as one moves away from the Earth's surface. The primary formula is g=g[RR+h]2g' = g\left[\frac{R}{R + h}\right]^2, which is applicable when the height hh is significant (specifically if h>R10h > \frac{R}{10}). For smaller heights where h<R10h < \frac{R}{10}, the approximation g=g(12hR)g' = g\left(1 - \frac{2h}{R}\right) is used. It should be noted that the value of gg never actually becomes zero with height, but rather approaches zero (g0g \rightarrow 0) as height approaches infinity (hh \rightarrow \infty).

Variation of Gravity with Depth (dd)

As one moves below the surface of the Earth, the acceleration due to gravity changes according to the linear relationship g=g(1dR)g' = g\left(1 - \frac{d}{R}\right). Here, gg represents the acceleration due to gravity at the Earth's surface. As the depth increases, the effective mass attracting the body decreases, leading to a decrease in gg'.

Variation of Gravity with Rotation and Latitude

The rotation of the Earth affects the apparent acceleration due to gravity because of centrifugal forces. This variation is represented by the formula g=g(1Rω2cos2(λ)g)g' = g\left(1 - \frac{R\omega^2 \cos^2(\lambda)}{g}\right). This effect means that the acceleration due to gravity is at its maximum at the poles and at its minimum at the equator.

Escape Velocity (vev_e)

Escape velocity is the minimum velocity required for a mass to escape from the gravitational field of a planet or the Earth. If a mass is launched with a velocity greater than or equal to the escape velocity, it will never return to the planet. The formula for escape velocity is ve=2GMRv_e = \sqrt{\frac{2GM}{R}}.

Planetary Motion and Orbital Velocity

Orbital velocity (vov_o) is the speed at which a planet or satellite moves in its orbit. It is derived from the requirement that the gravitational force provides the necessary centripetal force: GMmr2=mvo2r\frac{GMm}{r^2} = \frac{mv_o^2}{r}, which simplifies to vo=GMrv_o = \sqrt{\frac{GM}{r}}.

The orbital time period (TT) is defined as T=2πrvoT = \frac{2\pi r}{v_o}. Substituting the orbital velocity, the square of the time period can be expressed as T2=4π2r3GMT^2 = \frac{4\pi^2 r^3}{GM}.

The energy components of an orbiting body are as follows:

  1. Kinetic Energy (KEKE): KE=12mvo2=GMm2rKE = \frac{1}{2}m v_o^2 = \frac{GMm}{2r}.
  2. Potential Energy (PEPE): PE=GMmrPE = -\frac{GMm}{r}.
  3. Net Energy (EE): E=KE+PE=GMm2rE = KE + PE = -\frac{GMm}{2r}.

Kepler’s Laws of Planetary Motion

Kepler's First Law (Law of Orbits) states that all planets revolve around the Sun in elliptical orbits, with the Sun located at one of the two foci of the ellipse.

Kepler’s Second Law (Law of Areas) states that a radial line connecting a planet to the Sun sweeps out equal areas in equal intervals of time. This law is a direct consequence of the conservation of angular momentum because the torque about the axis of rotation is zero. The areal velocity is constant: dAdt=L2m=constant\frac{dA}{dt} = \frac{L}{2m} = \text{constant}. This implies that I1ω1=I2ω2I_1\omega_1 = I_2\omega_2, or (mr12)v1r1=(mr22)v2r2(mr_1^2)\frac{v_1}{r_1} = (mr_2^2)\frac{v_2}{r_2}, simplifying to v1r1=v2r2v_1r_1 = v_2r_2. Consequently, vperihelionvaphelion=raphelionrperihelion\frac{v_{\text{perihelion}}}{v_{\text{aphelion}}} = \frac{r_{\text{aphelion}}}{r_{\text{perihelion}}}. A planet moves faster when it is closer to the sun.

Kepler’s Third Law (Law of Periods) states that the square of the time period of a planet is proportional to the cube of the semi-major axis of its orbit (T2a3T^2 \propto a^3) or the cube of its mean orbital radius (T2r3T^2 \propto r^3).

If the eccentricity of the orbit is ee, the relationship between the distance and the semi-major axis is given by:

  • raphelion=rmax=a(1+e)r_{\text{aphelion}} = r_{\text{max}} = a(1 + e)
  • rperihelion=rmin=a(1e)r_{\text{perihelion}} = r_{\text{min}} = a(1 - e)

Weightlessness in a Satellite

In a satellite, the net force towards the center is the resultant of gravitational force and the normal contact force (NN). The equation of motion is GMmr2N=mv2r\frac{GMm}{r^2} - N = m\frac{v^2}{r}. Since the orbital acceleration is provided solely by the gravitational field (v2r=GMr2\frac{v^2}{r} = \frac{GM}{r^2}), the equation becomes GMmr2N=m(GMr2)\frac{GMm}{r^2} - N = m\left(\frac{GM}{r^2}\right). This results in N=0N = 0, meaning the surface of the satellite exerts no force on the body, resulting in apparent weightlessness.

Gravitational Potential of Specific Shapes

For a ring of mass MM and radius RR, the gravitational potential at any point on its axis at distance xx is V=GMR2+x2V = -\frac{GM}{\sqrt{R^2 + x^2}}. For this ring, the potential at the center is V=GMRV = -\frac{GM}{R}.

For a spherical shell (hollow sphere):

  • Inside and on surface (xRx \le R): Vin=Vsur=GMRV_{\text{in}} = V_{\text{sur}} = -\frac{GM}{R}.
  • Outside (x>Rx > R): Vout=GMxV_{\text{out}} = -\frac{GM}{x}.

For a solid sphere of radius RR:

  • Inside (0xR0 \le x \le R): Vin=GM(3R2x2)2R3V_{\text{in}} = -\frac{GM(3R^2 - x^2)}{2R^3}. Specifically, at the center (x=0x=0), V=3GM2RV = -\frac{3GM}{2R}.
  • On the surface (x=Rx = R): Vsur=GMRV_{\text{sur}} = -\frac{GM}{R}.
  • Outside (x>Rx > R): Vout=GMxV_{\text{out}} = -\frac{GM}{x}.

Gravitational Field Intensity of Specific Shapes

For a disc of mass MM and radius RR, the gravitational field intensity at a point PP on the axis at distance xx is E=2GMR2[1xR2+x2]E = \frac{2GM}{R^2} \left[1 - \frac{x}{\sqrt{R^2 + x^2}}\right], which can also be written as E=2GMR2(1cos(θ))E = \frac{2GM}{R^2} (1 - \cos(\theta)).

For a solid sphere:

  • Inside (x<Rx < R): Ein=GMxR3E_{\text{in}} = \frac{GMx}{R^3}.
  • On surface (x=Rx = R): Esur=GMR2E_{\text{sur}} = \frac{GM}{R^2}.
  • Outside (x>Rx > R): Eout=GMx2E_{\text{out}} = \frac{GM}{x^2}.

For a hollow sphere:

  • Inside (x<Rx < R): Ein=0E_{\text{in}} = 0.
  • On surface (x=Rx = R): Esur=GMR2E_{\text{sur}} = \frac{GM}{R^2}.
  • Outside (x>Rx > R): Eout=GMx2E_{\text{out}} = \frac{GM}{x^2}.