What JEE Adds to Gravitation

Sections 1 to 10 built this chapter properly, and that build is complete for the Board syllabus. But JEE — Main, and far more so Advanced — asks a family of gravitation questions the core never quite sets up. A rod pulls on a bead lying on its own axis. A stone is dropped down a hole bored clean through the Earth. A satellite is nudged out of one circle and has to arrive on another. Two stars circle a point in empty space between them.

None of it is harder physics. It is the same inverse-square law, asked from a place where you have to build the answer instead of quoting it.

Most of what follows sits outside the rationalised syllabus body text, but JEE Main and JEE Advanced ask it every year, so it is developed here from first principles.

The eight things this section teaches

# Skill Why it earns marks
1 Field and potential of an extended body by integration — rod, ring, disc The examiner picks a body with no formula. You have to build one
2 The shell and the solid sphere, inside and out, with BOTH graphs The single most drawn pair of graphs in the chapter
3 The tunnel through the Earth and the SHM it produces One set-up, one period, and a beautiful coincidence
4 Elliptical orbits with energy AND angular momentum together Two equations, two unknowns, every time
5 The vis-viva relation v2=GM(2r1a)v^2 = GM\left(\dfrac{2}{r}-\dfrac{1}{a}\right) One line that replaces both of the above
6 Orbit transfers and the Hohmann manoeuvre Two burns, and the arithmetic is always the same
7 Binary stars and the two-body reduction Both bodies move; the centre of mass never does
8 Gravitational self-energy, 3GM25R-\dfrac{3GM^2}{5R} Assembling a body, not just moving one

Running underneath all eight is one habit: decide whether you are computing a field or a potential, an energy or an energy per kilogram, and write the sign down before the number. More marks are lost in this chapter to a dropped minus sign than to any integral.

Conventions, fixed now

Constants. Throughout this section G=6.67×1011G = 6.67 \times 10^{-11} N m2^2/kg2^2, ME=5.97×1024M_E = 5.97 \times 10^{24} kg, RE=6.37×106R_E = 6.37 \times 10^{6} m, so GME=3.982×1014GM_E = 3.982 \times 10^{14} m3^3/s2^2. Where a surface value of gg is needed it is 9.8 m/s2^2. Every problem states which numbers it uses, and no problem mixes 9.8 with 10.

Notation. VV is the gravitational potential, measured in J/kg; UU is the gravitational potential energy of a body, measured in J, and U=mVU = mV. They are never the same object and the paper punishes anyone who treats them as one. g\vec{g} is the gravitational field (equivalently the acceleration a free particle would have), in N/kg or m/s2^2. aa is a semi-major axis, ee an eccentricity, rpr_p and rar_a the perigee and apogee distances.

The zero of potential is at infinity, always, in every line below. That is why VV and UU come out negative everywhere inside the universe of this chapter.

Key Point: Two relations connect everything in this section, and neither is new. V(r)=rgdrandg=dVdrr^V(r) = -\int_\infty^r \vec{g}\cdot d\vec{r} \qquad\text{and}\qquad \vec{g} = -\frac{dV}{dr}\,\hat{r} Potential is the running total of the field; the field is the slope of the potential. If you can compute either one, you already have the other — and the potential is almost always easier, because potentials add as plain numbers while fields add as vectors.

[Exam Tip] Before writing a single line, ask three questions. Am I inside the body or outside it? Is the quantity asked for a scalar (potential, energy) or a vector (field, force)? Which radius belongs in the formula — the body's RR, or the point's rr? Answer those three and most of these problems are already half solved.

Field and Potential of an Extended Body, by Integration

A point mass has g=GMr2r^\vec{g} = -\dfrac{GM}{r^2}\hat{r} and V=GMrV = -\dfrac{GM}{r}, and that is the end of it. A rod, a ring or a disc is a continuous smear of point masses, and you have to add up their contributions.

Rod, ring and disc with their mass elements, fields and potentials

The method, in four lines

  1. Choose an element dmdm every point of which is at the same distance from the field point PP. Get this right and the integral is easy; get it wrong and it is impossible.
  2. Express dmdm through the density: dm=λdxdm = \lambda\,dx for a wire, σdA\sigma\,dA for a lamina, ρdV\rho\,dV for a solid.
  3. Integrate the POTENTIAL first, because dV=GdmrdV = -\dfrac{G\,dm}{r} is a plain number with no direction to keep track of: V=GdmrV = -\int \frac{G\,dm}{r}
  4. Differentiate once to get the field: g=dVdxg = -\dfrac{dV}{dx} along the symmetry axis.

Integrating the field directly is legal, but then you must resolve every dgd\vec{g} into components and argue that the transverse ones cancel. The potential route skips that entirely.

The rod, along its own axis

A thin rod of mass MM and length LL lies on a line. PP is on that same line, a distance dd from the near end. An element of length dxdx at distance xx from PP has dm=MLdxdm = \dfrac{M}{L}dx:

V=dd+LGMLdxx=GMLln ⁣d+LdV = -\int_d^{d+L} \frac{G M}{L}\frac{dx}{x} = -\frac{GM}{L}\ln\!\frac{d+L}{d}

g=dVdd g=GMd(d+L) g = -\frac{dV}{dd} \quad\Longrightarrow\quad \boxed{\ g = \frac{GM}{d\,(d+L)}\ }

directed towards the rod. Notice that the answer is not GMd2\dfrac{GM}{d^2} and not GM(d+L/2)2\dfrac{GM}{(d + L/2)^2}: the near end pulls harder than the far end, so the effective distance is the geometric mean d(d+L)\sqrt{d(d+L)}, not the arithmetic one. Push dd out to 500 times LL and GMd(d+L)\dfrac{GM}{d(d+L)} sits 0.2% below GMd2\dfrac{GM}{d^2}, which is the point at which a distant rod finally counts as a point mass.

The ring, on its axis

Every point of a ring of radius aa is the same distance a2+x2\sqrt{a^2+x^2} from a point PP on the axis. So the potential needs no integration at all:

 V=GMa2+x2  g=GMx(a2+x2)3/2 \boxed{\ V = -\frac{GM}{\sqrt{a^2+x^2}}\ } \qquad\Longrightarrow\qquad \boxed{\ g = \frac{GMx}{(a^2+x^2)^{3/2}}\ }

Read off three things at once.

  • At the centre, x=0x = 0: the field is zero by symmetry, but the potential is GMa-\dfrac{GM}{a}, its most negative value. A ring is the cleanest example in the whole chapter of a place where g=0g = 0 and V0V \neq 0.
  • Far away, xax \gg a: gGMx2g \to \dfrac{GM}{x^2} and VGMxV \to -\dfrac{GM}{x}, as any body must.
  • In between, gg has a maximum. Setting dgdx=0\dfrac{dg}{dx} = 0 gives x=a2and theregmax=2GM33a2=0.385GMa2x = \frac{a}{\sqrt{2}} \qquad\text{and there}\qquad g_{max} = \frac{2GM}{3\sqrt{3}\,a^2} = 0.385\,\frac{GM}{a^2}

That maximum is asked directly, most years, in one form or another.

The disc, on its axis

A disc of mass MM and radius aa is a stack of rings. Take a ring of radius ss and width dsds; its share of the mass is dm=Mπa22πsds=2Msdsa2dm = \dfrac{M}{\pi a^2}\,2\pi s\,ds = \dfrac{2M s\,ds}{a^2}, and every point of it is s2+x2\sqrt{s^2+x^2} from PP:

V=0aGs2+x22Msdsa2=2GMa2(a2+x2x)V = -\int_0^a \frac{G}{\sqrt{s^2+x^2}}\cdot\frac{2Ms\,ds}{a^2} = -\frac{2GM}{a^2}\left(\sqrt{a^2+x^2}-x\right)

 g=2GMa2(1xa2+x2) \boxed{\ g = \frac{2GM}{a^2}\left(1 - \frac{x}{\sqrt{a^2+x^2}}\right)\ }

Two limits worth having on tap. Right at the centre of the face, x0x \to 0, the field is 2GMa2\dfrac{2GM}{a^2} — finite, not infinite. And writing σ=Mπa2\sigma = \dfrac{M}{\pi a^2} and letting aa \to \infty at fixed σ\sigma,

g2πGσg \to 2\pi G\sigma

independent of xx. An infinite sheet pulls with the same strength however far away you stand, exactly as a uniformly charged infinite sheet does in electrostatics. That parallel is worth carrying: every result in this block has an electrostatic twin with 14πϵ0\dfrac{1}{4\pi\epsilon_0} in place of G-G.

The table, so you never have to rebuild these under time pressure

Body Field point gg VV
Point mass MM distance rr GMr2\dfrac{GM}{r^2} GMr-\dfrac{GM}{r}
Rod, mass MM, length LL on its axis, dd from the near end GMd(d+L)\dfrac{GM}{d(d+L)} GMLlnd+Ld-\dfrac{GM}{L}\ln\dfrac{d+L}{d}
Ring, mass MM, radius aa on the axis, distance xx GMx(a2+x2)3/2\dfrac{GMx}{(a^2+x^2)^{3/2}} GMa2+x2-\dfrac{GM}{\sqrt{a^2+x^2}}
Ring at its centre 00 GMa-\dfrac{GM}{a}
Disc, mass MM, radius aa on the axis, distance xx 2GMa2(1xa2+x2)\dfrac{2GM}{a^2}\left(1-\dfrac{x}{\sqrt{a^2+x^2}}\right) 2GMa2(a2+x2x)-\dfrac{2GM}{a^2}\left(\sqrt{a^2+x^2}-x\right)
Infinite sheet, density σ\sigma anywhere 2πGσ2\pi G\sigma not defined (no zero at infinity)

Key Point: In every one of these, VV came first and gg came from dVdx-\dfrac{dV}{dx}. Check any result you derive by differentiating your own potential — if the two disagree, you have made an algebra slip, and you will find it in thirty seconds instead of losing the question.

[Advanced] The rod result generalises to a broken rod or a rod with a gap by simple subtraction, because both gg and VV are linear in the mass distribution. A rod from dd to d+Ld+L with a piece from d1d_1 to d2d_2 removed has g=GMfulld(d+L)g = \dfrac{GM_{full}}{d(d+L)} minus the missing piece's contribution, computed with its own limits. The same "add and subtract" trick is what makes the cavity problem in the next block work.

The Shell and the Solid Sphere: Both Graphs, Side by Side

Section 2 stated the two shell theorems and Section 4 drew gg against rr. What the paper actually asks for is both curves — the field and the potential — for both bodies, and the differences between them are exactly where the marks sit.

Field and potential against radius for a shell and a solid sphere

The four results for a uniform shell, mass MM, radius RR

rR:g=GMr2,V=GMrr<R:g=0,V=GMRr \ge R:\quad g = \frac{GM}{r^2},\quad V = -\frac{GM}{r} \qquad\qquad r < R:\quad g = 0,\quad V = -\frac{GM}{R}

Inside a shell the field vanishes but the potential does not. It sits at the constant value GMR-\dfrac{GM}{R}, the same as on the surface. That is not a contradiction: a constant potential has zero slope, and zero slope is zero field. Flat means constant, not zero.

The four results for a uniform solid sphere, mass MM, radius RR

rR:g=GMr2,V=GMrr \ge R:\quad g = \frac{GM}{r^2},\quad V = -\frac{GM}{r} r<R:g=GMrR3,V=GM(3R2r2)2R3r < R:\quad g = \frac{GMr}{R^3},\quad V = -\frac{GM\left(3R^2 - r^2\right)}{2R^3}

Inside, only the mass within radius rr pulls, and that mass is Mr3R3M\dfrac{r^3}{R^3}, so the field grows linearly from zero at the centre to GMR2\dfrac{GM}{R^2} at the surface. The potential is a downward parabola, deepest at the centre:

V(0)=3GM2R=1.5×V(surface)V(0) = -\frac{3GM}{2R} = 1.5 \times V(\text{surface})

Check the inside potential by differentiating it: dVdr=GMrR3-\dfrac{dV}{dr} = -\dfrac{GMr}{R^3}, whose magnitude is exactly the inside field. The two halves fit.

The whole thing in one table

Region Shell: gg Shell: VV Solid sphere: gg Solid sphere: VV
r>Rr > R GMr2\dfrac{GM}{r^2} GMr-\dfrac{GM}{r} GMr2\dfrac{GM}{r^2} GMr-\dfrac{GM}{r}
r=Rr = R GMR2\dfrac{GM}{R^2} GMR-\dfrac{GM}{R} GMR2\dfrac{GM}{R^2} GMR-\dfrac{GM}{R}
r<Rr < R 00 GMR-\dfrac{GM}{R} GMrR3\dfrac{GMr}{R^3} GM(3R2r2)2R3-\dfrac{GM(3R^2-r^2)}{2R^3}
r=0r = 0 00 GMR-\dfrac{GM}{R} 00 3GM2R-\dfrac{3GM}{2R}

Key Point — read the graphs, do not memorise them.

  • gg is discontinuous for a shell (it leaps from 00 to GMR2\dfrac{GM}{R^2} as you cross the surface) and continuous for a solid sphere.
  • VV is continuous for both. A jump in potential would mean an infinite field, and there is none.
  • gg is maximum at the surface of a solid sphere, never inside it and never outside it.
  • VV is most negative at the centre of a solid sphere and reaches zero only at infinity.
  • Both bodies look identical from outside. Nothing measured at r>Rr > R can tell a shell from a solid sphere of the same mass.

Superposition: concentric bodies add, cavities subtract

Because both g\vec{g} and VV are linear in mass, a compound body is handled by adding the pieces — being careful that each piece is evaluated at the same field point, with its own inside-or-outside rule.

Worked, in two lines. A solid sphere of mass MM and radius RR sits at the centre of a concentric thin shell of mass 2M2M and radius 3R3R.

At r=2Rr = 2R you are outside the sphere and inside the shell, so g=GM(2R)2+0=GM4R2,V=GM2R2GM3R=7GM6Rg = \frac{GM}{(2R)^2} + 0 = \frac{GM}{4R^2}, \qquad V = -\frac{GM}{2R} - \frac{2GM}{3R} = -\frac{7GM}{6R}

At the centre both fields vanish, but neither potential does: g=0,V=3GM2R2GM3R=13GM6Rg = 0, \qquad V = -\frac{3GM}{2R} - \frac{2GM}{3R} = -\frac{13GM}{6R}

That second line is the one students get wrong. The shell contributes nothing to gg and a great deal to VV.

The cavity. Scoop a sphere of radius R2\dfrac{R}{2} out of a uniform sphere of radius RR, with the hole's centre a distance dd from the big centre. Treat the hollow body as full sphere plus negative sphere. At any point PP inside the cavity, with r\vec{r} from the big centre and r=rd\vec{r}^{\,\prime} = \vec{r} - \vec{d} from the hole's centre,

g=43πGρr+43πGρr=43πGρd\vec{g} = -\frac{4}{3}\pi G\rho\,\vec{r} + \frac{4}{3}\pi G\rho\,\vec{r}^{\,\prime} = -\frac{4}{3}\pi G\rho\,\vec{d}

r\vec{r} has cancelled. The field everywhere inside the cavity is uniform, of magnitude 43πGρd=GMdR3\dfrac{4}{3}\pi G\rho d = \dfrac{GMd}{R^3} (with MM the mass of the original, unhollowed sphere), pointing from the cavity's centre towards the sphere's centre. A genuinely surprising result, and a JEE favourite.

[Exam Tip] The most expensive error in this block is using the wrong radius. In GMr2\dfrac{GM}{r^2} outside, rr is the distance of the field point. In GMR-\dfrac{GM}{R} inside a shell, RR is the shell's own radius and the field point's distance never appears. Write down which body you are inside before you write the formula, every single time.

The Tunnel Through the Earth

Bore a straight hole from one side of the Earth to the other, through the centre, evacuate it, and drop a stone in. What happens next is one of the most quoted results in the chapter — and, remarkably, one of the easiest to derive.

Tunnel through the Earth, the linear restoring force, and the resulting SHM

The force inside is proportional to the displacement

Treat the Earth as a uniform sphere of mass MEM_E and radius RER_E. At a distance rr from the centre, only the mass inside radius rr pulls, and the field is g(r)=GMErRE3g(r) = \dfrac{GM_E r}{R_E^3}, directed inwards. So for a stone of mass mm at position rr measured from the centre,

F=GMEmRE3rF = -\frac{GM_E m}{R_E^3}\,r

That is F=krF = -kr with k=GMEmRE3k = \dfrac{GM_E m}{R_E^3}: a linear restoring force, and therefore simple harmonic motion about the centre of the Earth.

ω2=km=GMERE3=gRE T=2πREg \omega^2 = \frac{k}{m} = \frac{GM_E}{R_E^3} = \frac{g}{R_E} \qquad\Longrightarrow\qquad \boxed{\ T = 2\pi\sqrt{\frac{R_E}{g}}\ }

The numbers

With RE=6.37×106R_E = 6.37 \times 10^6 m and g=9.8g = 9.8 m/s2^2:

T=2π6.37×1069.8=2π(806.2)=5.07×103 s=84.4 minutesT = 2\pi\sqrt{\frac{6.37\times 10^{6}}{9.8}} = 2\pi(806.2) = 5.07 \times 10^{3}\ \text{s} = 84.4\ \text{minutes}

Event Time
dropped at the surface, from rest 00
passing the centre, at maximum speed T/4=21.1T/4 = 21.1 min
momentarily at rest at the far surface T/2=42.2T/2 = 42.2 min
back where it started T=84.4T = 84.4 min

The maximum speed, at the centre, is the amplitude times ω\omega:

vmax=ωRE=gRE=(9.8)(6.37×106)=7.90×103 m/s=7.90 km/sv_{max} = \omega R_E = \sqrt{gR_E} = \sqrt{(9.8)(6.37\times 10^{6})} = 7.90 \times 10^{3}\ \text{m/s} = 7.90\ \text{km/s}

The coincidence that is not a coincidence

Look at that number. 7.90 km/s is exactly the orbital speed of a satellite skimming the Earth's surface, and 2πRE/g2\pi\sqrt{R_E/g} is exactly its orbital period. The stone in the tunnel and the surface-skimming satellite keep the same time.

That is not luck. A circular orbit of radius RER_E projected onto any straight line through the centre is simple harmonic motion of amplitude RER_E and angular frequency g/RE\sqrt{g/R_E}. The satellite's shadow on the tunnel and the stone in the tunnel obey the same equation, so they arrive together.

Key Point: T=2πREg84T = 2\pi\sqrt{\dfrac{R_E}{g}} \approx 84 minutes is the shortest period any object can have around the Earth, whether it goes through the planet or skims over it. Every higher orbit is slower. Remember the number: 84 minutes, or 5.07 kiloseconds.

The chord tunnel: the same 84 minutes

Now dig a tunnel that does not pass through the centre — a straight chord whose closest approach to the centre is bb. Let ss be the displacement along the tunnel from its midpoint. At that point r2=b2+s2r^2 = b^2 + s^2, the field has magnitude GMERE3r\dfrac{GM_E}{R_E^3}r pointing at the centre, and the component along the tunnel is that magnitude times sr\dfrac{s}{r}:

F=GMEmRE3rsr=GMEmRE3sF_{\parallel} = -\frac{GM_E m}{R_E^3}\,r \cdot \frac{s}{r} = -\frac{GM_E m}{R_E^3}\,s

The rr cancels again. Same spring constant, same ω\omega, same period of 84.4 minutes — for every chord, however short. Only the amplitude changes: it is the half-length of the chord, RE2b2\sqrt{R_E^2 - b^2}, so the maximum speed is ωRE2b2\omega\sqrt{R_E^2-b^2}, smaller for a shorter tunnel.

[Advanced] This is the physics behind the old "gravity train" proposal: a frictionless tunnel between any two points on Earth is a 42-minute one-way journey, whether it links Delhi and Mumbai or Delhi and Buenos Aires. The trip time is the same because the period is the same. If a question tells you the tunnel misses the centre and expects a different period, the question is wrong — but if it asks for the maximum speed, that genuinely does depend on the chord.

[Exam Tip] Two traps. First, TT depends only on RER_E and ggnot on the mass of the stone, not on the length of the tunnel. Second, the SHM picture uses the uniform-density Earth. The real Earth has a dense core, so the real field does not rise linearly from the centre and the real period would differ; every exam question assumes uniform density, and you should say so in one line if you are writing a full solution.

Elliptical Orbits: Two Conservation Laws, Used Together

A circular orbit has one unknown speed and one equation. An elliptical orbit has two speeds you care about — at perigee and at apogee — and needs two equations. Those two equations are always the same pair.

Elliptical orbit with perigee and apogee, and a Hohmann transfer between circles

The pair

Angular momentum. Gravity is a central force, so its torque about the focus is zero and L\vec{L} is constant. At perigee and apogee the velocity is perpendicular to the radius, so

 mvprp=mvaravpva=rarp \boxed{\ m v_p r_p = m v_a r_a \quad\Longrightarrow\quad \frac{v_p}{v_a} = \frac{r_a}{r_p}\ }

Energy. Gravity is conservative, so

 12vp2GMrp=12va2GMra \boxed{\ \frac{1}{2}v_p^2 - \frac{GM}{r_p} = \frac{1}{2}v_a^2 - \frac{GM}{r_a}\ }

Two equations, two unknowns. Solving them together, with rp=a(1e)r_p = a(1-e) and ra=a(1+e)r_a = a(1+e):

vp=GMa1+e1e,va=GMa1e1+e,vpva=GMav_p = \sqrt{\frac{GM}{a}\cdot\frac{1+e}{1-e}}, \qquad v_a = \sqrt{\frac{GM}{a}\cdot\frac{1-e}{1+e}}, \qquad v_p v_a = \frac{GM}{a}

That last relation is a lovely, and rarely quoted, check: the product of the perigee and apogee speeds equals the circular speed for radius aa, squared.

The geometry you must be able to write down cold

Quantity In terms of aa and ee In terms of rpr_p and rar_a
perigee distance rp=a(1e)r_p = a(1-e)
apogee distance ra=a(1+e)r_a = a(1+e)
semi-major axis aa a=rp+ra2a = \dfrac{r_p+r_a}{2}
eccentricity ee e=rarpra+rpe = \dfrac{r_a-r_p}{r_a+r_p}
semi-minor axis b=a1e2b = a\sqrt{1-e^2} b=rprab = \sqrt{r_p r_a}
total energy E=GMm2aE = -\dfrac{GMm}{2a} same, with aa from the sum
period T=2πa3GMT = 2\pi\sqrt{\dfrac{a^3}{GM}} same

Key Point: EE and TT depend on aa alone — not on ee. Two satellites with the same semi-major axis have the same total energy and the same period, however differently shaped their orbits are. A circle is just the special case e=0e = 0, a=ra = r, and every formula above collapses to the circular one when you put e=0e=0.

Vis-viva: the one line that does both jobs

Combine the two conservation laws and eliminate everything but rr:

 v2=GM(2r1a) \boxed{\ v^2 = GM\left(\frac{2}{r} - \frac{1}{a}\right)\ }

This is the vis-viva equation, and it is the single most useful formula in orbital mechanics. It gives the speed at any point of any orbit — circle, ellipse, parabola or hyperbola — from just two numbers.

Feed it the special cases and watch them fall out:

Put in You get Meaning
a=ra = r v=GMrv = \sqrt{\dfrac{GM}{r}} the circular orbit speed
aa \to \infty v=2GMrv = \sqrt{\dfrac{2GM}{r}} the escape speed — a parabola is an "ellipse" with aa infinite
r=rp=a(1e)r = r_p = a(1-e) vp=GMa1+e1ev_p = \sqrt{\dfrac{GM}{a}\dfrac{1+e}{1-e}} fastest point
r=ra=a(1+e)r = r_a = a(1+e) va=GMa1e1+ev_a = \sqrt{\dfrac{GM}{a}\dfrac{1-e}{1+e}} slowest point
r=ar = a (ends of the minor axis) v=GMav = \sqrt{\dfrac{GM}{a}} the same speed a circular orbit of radius aa would have

That last row is worth a moment. At the two points where the satellite is exactly aa from the focus, it is moving at precisely the circular speed for that distance — but not in a circle, because its velocity is not perpendicular to the radius there.

Areal velocity, in the form that gets asked

Kepler's second law is angular-momentum conservation wearing a hat:

dAdt=12r2dθdt=L2m=constant\frac{dA}{dt} = \frac{1}{2}r^2\frac{d\theta}{dt} = \frac{L}{2m} = \text{constant}

so the time to sweep any region equals its area divided by L2m\dfrac{L}{2m}. Sweeping the whole ellipse takes TT, and the whole area is πab\pi a b, which gives T=2πabmLT = \dfrac{2\pi a b m}{L} — another route to Kepler's third law if you ever need it.

[Exam Tip] When a problem gives you one speed and one distance on an ellipse and asks for another, do not set up the two conservation laws from scratch. Get aa from vis-viva at the point you know, then use vis-viva again at the point you want. Two substitutions, no simultaneous equations.

Orbit Transfers, the Hohmann Manoeuvre, and Catching a Satellite

A spacecraft in a circular orbit of radius r1r_1 has to end up in a circular orbit of radius r2r_2. It cannot simply "move outwards" — every path it takes is itself an orbit. The cheapest route, and the only one you are asked about, is the Hohmann transfer: half an ellipse that touches the inner circle at its perigee and the outer circle at its apogee.

The transfer ellipse

rp=r1,ra=r2, at=r1+r22 r_p = r_1, \qquad r_a = r_2, \qquad \boxed{\ a_t = \frac{r_1+r_2}{2}\ }

Burn 1, at r1r_1: speed up from the circular speed to the transfer ellipse's perigee speed. v1=GMr1,vp=GM(2r11at)=2GMr2r1(r1+r2)v_1 = \sqrt{\frac{GM}{r_1}}, \qquad v_p = \sqrt{GM\left(\frac{2}{r_1}-\frac{1}{a_t}\right)} = \sqrt{\frac{2GM r_2}{r_1(r_1+r_2)}}  Δv1=vpv1=GMr1(2r2r1+r21) \boxed{\ \Delta v_1 = v_p - v_1 = \sqrt{\frac{GM}{r_1}}\left(\sqrt{\frac{2r_2}{r_1+r_2}} - 1\right)\ }

Coast half an ellipse. That takes half the transfer orbit's period:  ttransfer=πat3GM \boxed{\ t_{transfer} = \pi\sqrt{\frac{a_t^{3}}{GM}}\ }

Burn 2, at r2r_2: the craft arrives at apogee moving too slowly to hold a circular orbit there, so it must speed up again. va=2GMr1r2(r1+r2),v2=GMr2v_a = \sqrt{\frac{2GM r_1}{r_2(r_1+r_2)}}, \qquad v_2 = \sqrt{\frac{GM}{r_2}}  Δv2=v2va=GMr2(12r1r1+r2) \boxed{\ \Delta v_2 = v_2 - v_a = \sqrt{\frac{GM}{r_2}}\left(1 - \sqrt{\frac{2r_1}{r_1+r_2}}\right)\ }

Key Point: Both burns are prograde — both speed the craft up — yet the craft ends up SLOWER than it started. A low orbit is a fast orbit. The engine adds kinetic energy twice, and the orbit's total energy rises from GMm2r1-\dfrac{GMm}{2r_1} to GMm2r2-\dfrac{GMm}{2r_2}, which is less negative; but the speed GM/r\sqrt{GM/r} has fallen. Energy up, speed down. Section 8 called this the orbit-raising paradox; here is the manoeuvre that actually does it.

The ledger, in one worked line

Going from a 7000 km orbit to the geostationary radius of 42 200 km, with GME=3.982×1014GM_E = 3.982\times10^{14}:

Step Value
circular speed at r1r_1 7.54 km/s
perigee speed on the transfer ellipse 9.88 km/s
Δv1\Delta v_1 2.34 km/s
apogee arrival speed 1.64 km/s
circular speed at r2r_2 3.07 km/s
Δv2\Delta v_2 1.43 km/s
total Δv\Delta v 3.77 km/s
coast time 19 200 s = 5.34 h

Those are real launch numbers: getting a communications satellite from a parking orbit to geostationary really does cost a little under 4 km/s of velocity change and a little over five hours of coasting.

Two satellites in the same orbit: the docking paradox

Here is the set-up that catches everybody. Two satellites are in the same circular orbit, and one is a few degrees behind the other. How does the chaser catch up?

Not by firing its engine forward. Speeding up raises the orbit, which lengthens the period, which makes the chaser fall further behind.

fire forward    higher orbit    longer period    you lose ground\text{fire forward} \;\Rightarrow\; \text{higher orbit} \;\Rightarrow\; \text{longer period} \;\Rightarrow\; \text{you lose ground}

The chaser must fire backwards. That drops it into a slightly smaller ellipse whose period is shorter, so it gains a little angle on every lap. After enough laps the gap is closed, and a second forward burn puts it back on the original circle, now alongside the target. This is called a phasing orbit, and it is exactly how a spacecraft docks with a space station.

Tphasing=T(1Δθ2πN),aphasing=(GMTphasing24π2)1/3T_{phasing} = T\left(1 - \frac{\Delta\theta}{2\pi N}\right), \qquad a_{phasing} = \left(\frac{GM\,T_{phasing}^2}{4\pi^2}\right)^{1/3}

where Δθ\Delta\theta is the angle to be made up and NN the number of laps you are prepared to spend. Fewer laps means a bigger drop and a bigger Δv\Delta v; the manoeuvre is cheap only if you are patient.

Two satellites in different orbits: when do they line up?

If two satellites have periods T1T_1 and T2T_2, their angular speeds differ by 2πT12πT2\left\lvert \frac{2\pi}{T_1} - \frac{2\pi}{T_2}\right\rvert, so they return to the same relative position after

 1Tsyn=1T11T2 \boxed{\ \frac{1}{T_{syn}} = \left\lvert\frac{1}{T_1} - \frac{1}{T_2}\right\rvert\ }

the synodic period. If they orbit in opposite senses, the angular speeds add instead: 1Tmeet=1T1+1T2\frac{1}{T_{meet}} = \frac{1}{T_1} + \frac{1}{T_2}

[Exam Tip] Read the sense of rotation before you pick a formula. "Same direction" gives a difference, "opposite directions" gives a sum, and swapping them is the single most common slip in this family of questions.

Binary Stars, Self-Energy, and the Traps That Cost Marks

Two bodies, one centre of mass

Everything so far assumed a tiny satellite around a huge fixed planet. Drop that assumption and both bodies move — around their common centre of mass, which stays put because no external force acts.

Two stars of masses m1m_1 and m2m_2, separated by a distance dd, in circular orbits. The centre of mass divides the separation in the inverse ratio of the masses:

r1=m2dm1+m2,r2=m1dm1+m2,r1+r2=dr_1 = \frac{m_2\,d}{m_1+m_2}, \qquad r_2 = \frac{m_1\,d}{m_1+m_2}, \qquad r_1 + r_2 = d

Both stars sweep the same angle in the same time, so they share one ω\omega and one period. For star 1, the gravitational pull supplies its centripetal need:

Gm1m2d2=m1ω2r1=m1ω2m2dm1+m2 ω2=G(m1+m2)d3 \frac{Gm_1m_2}{d^2} = m_1\omega^2 r_1 = m_1\omega^2\frac{m_2 d}{m_1+m_2} \quad\Longrightarrow\quad \boxed{\ \omega^2 = \frac{G(m_1+m_2)}{d^3}\ }

 T=2πd3G(m1+m2) \boxed{\ T = 2\pi\sqrt{\frac{d^3}{G(m_1+m_2)}}\ }

which is Kepler's third law with the total mass in the denominator. Everything else follows:

Quantity Result The sentence that goes with it
speeds v1=ωr1v_1 = \omega r_1, v2=ωr2v_2 = \omega r_2 v1v2=m2m1\dfrac{v_1}{v_2} = \dfrac{m_2}{m_1}: the heavier star moves slower
momenta m1v1=m2v2m_1 v_1 = m_2 v_2 equal and opposite; total momentum is zero
total KE 12μvrel2\dfrac{1}{2}\mu v_{rel}^2 with vrel=v1+v2v_{rel} = v_1+v_2 one particle again
total energy E=Gm1m22aE = -\dfrac{Gm_1m_2}{2a} aa is the semi-major axis of the relative orbit; a=da = d if circular

The reduced mass

The relative position r=r2r1\vec{r} = \vec{r}_2 - \vec{r}_1 obeys

μd2rdt2=Gm1m2r2r^,where μ=m1m2m1+m2 \mu\,\frac{d^2\vec{r}}{dt^2} = -\frac{Gm_1m_2}{r^2}\hat{r}, \qquad\text{where}\qquad \boxed{\ \mu = \frac{m_1m_2}{m_1+m_2}\ }

So the two-body problem is exactly a one-body problem: a single particle of mass μ\mu orbiting a fixed centre that exerts the same force. Every orbital result you know — vis-viva, T2a3T^2 \propto a^3, E=Gm1m22aE = -\dfrac{Gm_1m_2}{2a} — applies to the relative orbit unchanged.

Two sanity checks on μ\mu. If m2m1m_2 \gg m_1 then μm1\mu \to m_1, and the light body orbits a fixed heavy one, which is where we started. If m1=m2=mm_1 = m_2 = m then μ=m2\mu = \dfrac{m}{2}, and the two stars sit diametrically opposite at d2\dfrac{d}{2} each, moving at v=Gm2dv = \sqrt{\dfrac{Gm}{2d}}.

Gravitational self-energy

Everything so far has been the energy of a body in a field. Self-energy is different: it is the work done by gravity in assembling the body from dust scattered at infinity, or equivalently minus the work you must do to pull it apart and scatter it again.

Build a uniform sphere shell by shell. When the sphere already has radius rr, its mass is m=Mr3R3m = M\dfrac{r^3}{R^3}, and the next shell of thickness drdr carries dm=3Mr2R3drdm = \dfrac{3Mr^2}{R^3}dr. Bringing that shell from infinity onto a sphere of mass mm releases

dU=Gmdmr=3GM2r4R6drdU = -\frac{G\,m\,dm}{r} = -\frac{3GM^2 r^4}{R^6}\,dr

Integrate from 00 to RR:

 Uself=35GM2R \boxed{\ U_{self} = -\frac{3}{5}\frac{GM^2}{R}\ }

For a thin shell, where every scrap of mass ends up at the same radius RR, the same argument gives Uself=GM22RU_{self} = -\dfrac{GM^2}{2R}.

The Earth's number. With ME=5.97×1024M_E = 5.97\times10^{24} kg and RE=6.37×106R_E = 6.37\times10^{6} m,

Uself=3(6.67×1011)(5.97×1024)25(6.37×106)=2.24×1032 JU_{self} = -\frac{3(6.67\times10^{-11})(5.97\times10^{24})^2}{5(6.37\times10^{6})} = -2.24\times10^{32}\ \text{J}

so blowing the Earth apart and dispersing it to infinity would cost +2.24×1032+2.24\times10^{32} J — about seven days of the Sun's entire energy output, in every direction at once.

Key Point: Do not confuse the three energies. U=GMEmrU = -\dfrac{GM_Em}{r} is the energy of a body mm in the Earth's field. E=GMEm2rE = -\dfrac{GM_Em}{2r} is the total energy of that body in orbit. Uself=3GME25REU_{self} = -\dfrac{3GM_E^2}{5R_E} is the energy stored in the Earth itself, and it involves ME2M_E^2, not MEmM_E m. Look at which masses appear before you pick a formula.

The seven traps, collected

  1. Potential is not potential energy. VV is J/kg, U=mVU = mV is J. A question asking for "the potential at the centre of the sphere" wants 3GM2R-\dfrac{3GM}{2R}, not 3GMm2R-\dfrac{3GMm}{2R}.
  2. Zero field is not zero potential. Inside a shell, at the centre of a ring, at the centre of a two-mass system: g=0g = 0 and V0V \neq 0 at all three.
  3. Superposition of fields is a VECTOR sum. Two equal masses on opposite sides of a point give zero field and a doubled (more negative) potential. Adding their magnitudes is the classic sign error.
  4. The wrong radius. GMR-\dfrac{GM}{R} inside a shell uses the shell's radius; GMr2\dfrac{GM}{r^2} outside uses the field point's distance; GME(RE+h)2\dfrac{GM_E}{(R_E+h)^2} uses the sum. Three different radii, three different formulas.
  5. TT and EE of an ellipse depend on aa, not on rr. Do not put the current distance where the semi-major axis belongs.
  6. Both Hohmann burns speed you up, and you end up slower. If your working has a retro-burn in a raising manoeuvre, you have the sign backwards.
  7. Self-energy carries M2M^2. If your answer to a self-energy question is linear in MM, it is wrong before you check the number.

[Advanced] One more relation worth carrying. For any bound orbit under an inverse-square force, the time-averaged kinetic and potential energies satisfy U=2K\langle U\rangle = -2\langle K\rangle, so E=K=12UE = -\langle K\rangle = \dfrac{1}{2}\langle U\rangle. That is the virial theorem, and for a circle it is exactly the K:U:E=1:2:1K : U : E = 1 : -2 : -1 ratio Section 8 established. For an ellipse the instantaneous values wander, but the averages over one full period obey it exactly.

Solved Examples

Twelve problems at Advanced level. Try each on paper before reading the solution — the value here is in the set-up, not the arithmetic. Where numbers are used, G=6.67×1011G = 6.67\times10^{-11}, ME=5.97×1024M_E = 5.97\times10^{24} kg, RE=6.37×106R_E = 6.37\times10^{6} m and g=9.8g = 9.8 m/s2^2.

Example 1: A ring, on its own axis

A uniform ring of mass MM and radius aa lies in the yzyz-plane with its centre at the origin. For a point PP on the xx-axis a distance xx from the centre, find (a) the potential VV, (b) the field gg, (c) the value of xx at which gg is largest and the value of gg there, and (d) the speed with which a particle released from rest at x=3ax = \sqrt{3}\,a passes through the centre of the ring.

Solution:

  1. (a) The potential needs no integral. Every element of the ring is the same distance a2+x2\sqrt{a^2+x^2} from PP, so V=Ga2+x2dm=GMa2+x2V = -\frac{G}{\sqrt{a^2+x^2}}\int dm = -\frac{GM}{\sqrt{a^2+x^2}}

  2. (b) Differentiate once. g=dVdx=ddx[GM(a2+x2)1/2]=GMx(a2+x2)3/2g = -\frac{dV}{dx} = -\frac{d}{dx}\left[-GM\left(a^2+x^2\right)^{-1/2}\right] = -\frac{GMx}{(a^2+x^2)^{3/2}} The minus sign says the field points back towards the ring. Its magnitude is GMx(a2+x2)3/2\dfrac{GMx}{(a^2+x^2)^{3/2}}, which is zero at the centre and zero at infinity.

  3. (c) The maximum. Set the derivative to zero: ddx[x(a2+x2)3/2]=(a2+x2)3/2x3x(a2+x2)1/2(a2+x2)3=0  a2+x2=3x2\frac{d}{dx}\left[\frac{x}{(a^2+x^2)^{3/2}}\right]= \frac{(a^2+x^2)^{3/2} - x\cdot 3x(a^2+x^2)^{1/2}}{(a^2+x^2)^{3}} = 0 \ \Longrightarrow\ a^2 + x^2 = 3x^2 x=a2,gmax=GM(a/2)(32a2)3/2=2GM33a2=0.385GMa2x = \frac{a}{\sqrt{2}}, \qquad g_{max} = \frac{GM\,(a/\sqrt{2})}{\left(\tfrac{3}{2}a^2\right)^{3/2}} = \frac{2GM}{3\sqrt{3}\,a^2} = 0.385\,\frac{GM}{a^2}

  4. (d) Energy conservation, with potentials. At x=3ax = \sqrt{3}a, a2+3a2=2a\sqrt{a^2+3a^2} = 2a, so Vi=GM2a,Vf=GMaV_i = -\frac{GM}{2a}, \qquad V_f = -\frac{GM}{a} 12v2=ViVf=GM2a+GMa=GM2av=GMa\frac{1}{2}v^2 = V_i - V_f = -\frac{GM}{2a}+\frac{GM}{a} = \frac{GM}{2a} \quad\Longrightarrow\quad v = \sqrt{\frac{GM}{a}}

Final Answer: (a) GMa2+x2-\dfrac{GM}{\sqrt{a^2+x^2}}; (b) GMx(a2+x2)3/2\dfrac{GMx}{(a^2+x^2)^{3/2}} towards the ring; (c) x=a2x = \dfrac{a}{\sqrt{2}}, gmax=2GM33a2g_{max} = \dfrac{2GM}{3\sqrt{3}a^2}; (d) v=GMav = \sqrt{\dfrac{GM}{a}}.

Takeaway: Integrate the potential, then differentiate. The ring is the one body where the potential integral is free, and everything else follows from one derivative. And note part (d): the particle speeds up all the way to the centre even though the field is zero there — the field being zero says nothing about the potential.

Example 2: A rod, and the distance that is not the average

A thin uniform rod of mass MM and length LL lies along the xx-axis. A point PP lies on the same axis, a distance dd from the nearer end. Find (a) the potential at PP, (b) the field at PP, (c) both at d=L2d = \dfrac{L}{2}, and (d) the percentage by which the exact field differs from the point-mass estimate GMd2\dfrac{GM}{d^2} when d=500Ld = 500L.

Solution:

  1. (a) One integral. Take an element dxdx at distance xx from PP, with dm=MLdxdm = \dfrac{M}{L}dx. Every point of it is at xx, so V=dd+LGxMLdx=GML[lnx]dd+L=GMLln ⁣d+LdV = -\int_d^{d+L}\frac{G}{x}\cdot\frac{M}{L}dx = -\frac{GM}{L}\left[\ln x\right]_d^{d+L} = -\frac{GM}{L}\ln\!\frac{d+L}{d}

  2. (b) Differentiate with respect to dd. dVdd=GML(1d+L1d)=GMLLd(d+L)=GMd(d+L)\frac{dV}{dd} = -\frac{GM}{L}\left(\frac{1}{d+L}-\frac{1}{d}\right) = \frac{GM}{L}\cdot\frac{L}{d(d+L)} = \frac{GM}{d(d+L)} so the field has magnitude GMd(d+L)\dfrac{GM}{d(d+L)}, directed towards the rod. (Getting the same answer by integrating Gdmx2\dfrac{G\,dm}{x^2} directly is a useful check: dd+LGMLdxx2=GML(1d1d+L)\displaystyle\int_d^{d+L}\frac{GM}{L}\frac{dx}{x^2} = \frac{GM}{L}\left(\frac{1}{d}-\frac{1}{d+L}\right), the same thing.)

  3. (c) Put d=L/2d = L/2. g=GM(L2)(3L2)=4GM3L2=1.333GML2g = \frac{GM}{\left(\frac{L}{2}\right)\left(\frac{3L}{2}\right)} = \frac{4GM}{3L^2} = 1.333\,\frac{GM}{L^2} V=GMLln3=1.099GMLV = -\frac{GM}{L}\ln 3 = -1.099\,\frac{GM}{L}

  4. (d) The far-field check. With d=500Ld = 500L, gexact=GM500L501L=GM250500L2,gpoint=GM250000L2g_{exact} = \frac{GM}{500L\cdot 501L} = \frac{GM}{250\,500\,L^2}, \qquad g_{point} = \frac{GM}{250\,000\,L^2} gpointgexact=250500250000=1.002\frac{g_{point}}{g_{exact}} = \frac{250\,500}{250\,000} = 1.002 The point-mass estimate is 0.2% too big, even at five hundred rod-lengths away.

Final Answer: (a) GMLlnd+Ld-\dfrac{GM}{L}\ln\dfrac{d+L}{d}; (b) GMd(d+L)\dfrac{GM}{d(d+L)}; (c) 4GM3L2\dfrac{4GM}{3L^2} and GMLln3-\dfrac{GM}{L}\ln 3; (d) 0.2% high.

Takeaway: The effective distance to a rod on its own axis is the geometric mean d(d+L)\sqrt{d(d+L)}, not the arithmetic mean d+L2d + \dfrac{L}{2}. Every extended body approaches a point mass eventually — but "eventually" is further away than students expect.

Example 3: A disc, and the sheet it becomes

A uniform disc of mass MM and radius aa has a point PP on its axis at distance xx from the centre. (a) Derive VV and gg at PP. (b) Evaluate both at x=0.75ax = 0.75a. (c) Show that as aa \to \infty at fixed surface density σ\sigma, the field becomes independent of xx, and state its value.

Solution:

  1. (a) Slice into rings. A ring of radius ss and width dsds has area 2πsds2\pi s\,ds and therefore mass dm=Mπa22πsds=2Msdsa2dm = \frac{M}{\pi a^2}\,2\pi s\,ds = \frac{2Ms\,ds}{a^2} Every point of it is s2+x2\sqrt{s^2+x^2} from PP, so it contributes dV=Gdms2+x2dV = -\dfrac{G\,dm}{\sqrt{s^2+x^2}} and V=2GMa20asdss2+x2=2GMa2[s2+x2]0a=2GMa2(a2+x2x)V = -\frac{2GM}{a^2}\int_0^a \frac{s\,ds}{\sqrt{s^2+x^2}} = -\frac{2GM}{a^2}\left[\sqrt{s^2+x^2}\right]_0^a = -\frac{2GM}{a^2}\left(\sqrt{a^2+x^2}-x\right) g=dVdx=2GMa2(xa2+x21)g=2GMa2(1xa2+x2)g = -\frac{dV}{dx} = -\frac{2GM}{a^2}\left(\frac{x}{\sqrt{a^2+x^2}} - 1\right) \quad\Longrightarrow\quad \lvert g\rvert = \frac{2GM}{a^2}\left(1-\frac{x}{\sqrt{a^2+x^2}}\right)

  2. (b) At x=0.75ax = 0.75a. Then a2+0.5625a2=1.25a\sqrt{a^2+0.5625a^2} = 1.25a, so g=2GMa2(10.751.25)=2GMa2(0.40)=0.800GMa2g = \frac{2GM}{a^2}\left(1-\frac{0.75}{1.25}\right) = \frac{2GM}{a^2}(0.40) = 0.800\,\frac{GM}{a^2} V=2GMa2(1.25a0.75a)=2GMa2(0.5a)=GMaV = -\frac{2GM}{a^2}\left(1.25a - 0.75a\right) = -\frac{2GM}{a^2}(0.5a) = -\frac{GM}{a}

  3. (c) The infinite sheet. Write M=σπa2M = \sigma\pi a^2, so 2GMa2=2πGσ\dfrac{2GM}{a^2} = 2\pi G\sigma: g=2πGσ(1xa2+x2) a 2πGσ(10)=2πGσg = 2\pi G\sigma\left(1 - \frac{x}{\sqrt{a^2+x^2}}\right) \xrightarrow[\ a\to\infty\ ]{} 2\pi G\sigma\left(1-0\right) = 2\pi G\sigma The xx has vanished. An infinite sheet pulls equally hard at every distance.

Final Answer: (a) as boxed above; (b) g=0.800GMa2g = 0.800\dfrac{GM}{a^2}, V=GMaV = -\dfrac{GM}{a}; (c) g2πGσg \to 2\pi G\sigma, independent of xx.

Takeaway: A disc is a stack of rings, and each ring is already solved. Note also what happens at x=0x=0: the field on the face of a disc is a finite 2GMa2\dfrac{2GM}{a^2}, not infinite — the mass right under your feet is spread out, and its contributions cancel sideways.

Solved Examples (continued)

Example 4: The Earth, inside and out, in field and in potential

Model the Earth as a uniform sphere of mass 5.97×10245.97\times10^{24} kg and radius 6.37×1066.37\times10^{6} m. Find gg and VV at (a) the centre, (b) r=RE2r = \dfrac{R_E}{2}, (c) the surface, and (d) r=2REr = 2R_E. (e) State where gg is largest and where VV is smallest.

Solution:

  1. Set up the two constants once. GME=(6.67×1011)(5.97×1024)=3.982×1014 m3/s2GM_E = (6.67\times10^{-11})(5.97\times10^{24}) = 3.982\times10^{14}\ \text{m}^3\text{/s}^2 GMERE=3.982×10146.37×106=6.251×107 J/kg,GMERE2=9.81 m/s2\frac{GM_E}{R_E} = \frac{3.982\times10^{14}}{6.37\times10^{6}} = 6.251\times10^{7}\ \text{J/kg}, \qquad \frac{GM_E}{R_E^2} = 9.81\ \text{m/s}^2

  2. (a) At the centre, r=0r = 0, use the inside formulas: g=GME(0)RE3=0,V=3GME2RE=1.5(6.251×107)=9.38×107 J/kgg = \frac{GM_E(0)}{R_E^3} = 0, \qquad V = -\frac{3GM_E}{2R_E} = -1.5(6.251\times10^{7}) = -9.38\times10^{7}\ \text{J/kg}

  3. (b) At r=RE/2r = R_E/2, still inside: g=GMERE2rRE=9.81×0.5=4.91 m/s2g = \frac{GM_E}{R_E^2}\cdot\frac{r}{R_E} = 9.81\times 0.5 = 4.91\ \text{m/s}^2 V=GME(3RE2RE24)2RE3=118GMERE=1.375(6.251×107)=8.60×107 J/kgV = -\frac{GM_E\left(3R_E^2 - \frac{R_E^2}{4}\right)}{2R_E^3} = -\frac{11}{8}\frac{GM_E}{R_E} = -1.375(6.251\times10^{7}) = -8.60\times10^{7}\ \text{J/kg}

  4. (c) At the surface, both formulas agree: g=9.81 m/s2,V=GMERE=6.25×107 J/kgg = 9.81\ \text{m/s}^2, \qquad V = -\frac{GM_E}{R_E} = -6.25\times10^{7}\ \text{J/kg}

  5. (d) At r=2REr = 2R_E, outside: g=GME(2RE)2=9.814=2.45 m/s2,V=GME2RE=3.13×107 J/kgg = \frac{GM_E}{(2R_E)^2} = \frac{9.81}{4} = 2.45\ \text{m/s}^2, \qquad V = -\frac{GM_E}{2R_E} = -3.13\times10^{7}\ \text{J/kg}

  6. (e) The extremes. gg rises linearly to a maximum at the surface and falls as 1r2\dfrac{1}{r^2} beyond it. VV is most negative at the centre and climbs monotonically to zero at infinity. The two extremes are at different places, and that catches people.

Final Answer:

rr gg (m/s2^2) VV (J/kg)
00 00 9.38×107-9.38\times10^{7}
RE/2R_E/2 4.914.91 8.60×107-8.60\times10^{7}
RER_E 9.819.81 6.25×107-6.25\times10^{7}
2RE2R_E 2.452.45 3.13×107-3.13\times10^{7}

Takeaway: gg peaks at the surface; VV bottoms out at the centre. Going down a mine weakens the field but deepens the potential well, which is why you must still do work to climb back out. Check every inside-potential answer against V(0)=1.5V(R)V(0) = 1.5\,V(R).

Example 5: The cavity that holds a uniform field

A uniform solid sphere of mass M=6.0×1020M = 6.0\times10^{20} kg and radius R=4.0×105R = 4.0\times10^{5} m has a spherical cavity of radius R2\dfrac{R}{2} carved out of it, the cavity's centre lying a distance R2\dfrac{R}{2} from the sphere's centre. Find (a) the mass of the remaining body, (b) the gravitational field everywhere inside the cavity, and (c) its numerical value.

Solution:

  1. (a) The missing mass. The cavity has half the radius, hence (12)3=18\left(\tfrac{1}{2}\right)^3 = \tfrac{1}{8} of the volume, so it removes M8\dfrac{M}{8}: Mremaining=MM8=7M8=7(6.0×1020)8=5.25×1020 kgM_{remaining} = M - \frac{M}{8} = \frac{7M}{8} = \frac{7(6.0\times10^{20})}{8} = 5.25\times10^{20}\ \text{kg}

  2. (b) Superposition with a negative sphere. Model the hollowed body as a full sphere of density ρ\rho plus a sphere of density ρ-\rho filling the cavity. At a point PP inside the cavity, let r\vec{r} run from the big centre to PP and d\vec{d} from the big centre to the cavity's centre, so that r=rd\vec{r}^{\,\prime} = \vec{r}-\vec{d} runs from the cavity's centre to PP. Using ginside=43πGρr\vec{g}_{inside} = -\dfrac{4}{3}\pi G\rho\,\vec{r} for each piece: g=43πGρr  +  43πGρr=43πGρ(rr)=43πGρd\vec{g} = -\frac{4}{3}\pi G\rho\,\vec{r} \;+\; \frac{4}{3}\pi G\rho\,\vec{r}^{\,\prime} = -\frac{4}{3}\pi G\rho\left(\vec{r}-\vec{r}^{\,\prime}\right) = -\frac{4}{3}\pi G\rho\,\vec{d} r\vec{r} has cancelled. The field is the same vector at every point of the cavity: uniform, and directed from the cavity's centre towards the sphere's centre.

  3. (c) The number. Writing 43πGρ=GMR3\dfrac{4}{3}\pi G\rho = \dfrac{GM}{R^3} (with MM the mass of the original, unhollowed sphere) and d=R2d = \dfrac{R}{2}: g=GMR3R2=GM2R2=(6.67×1011)(6.0×1020)2(4.0×105)2=4.00×10103.2×1011=0.125 m/s2g = \frac{GM}{R^3}\cdot\frac{R}{2} = \frac{GM}{2R^2} = \frac{(6.67\times10^{-11})(6.0\times10^{20})}{2(4.0\times10^{5})^2} = \frac{4.00\times10^{10}}{3.2\times10^{11}} = 0.125\ \text{m/s}^2

Final Answer: (a) 5.25×10205.25\times10^{20} kg; (b) uniform, magnitude 43πGρd=GMdR3\dfrac{4}{3}\pi G\rho d = \dfrac{GMd}{R^3}, pointing from the cavity's centre towards the sphere's centre; (c) 0.125 m/s2^2.

Takeaway: Negative mass is a legal tool because both g\vec{g} and VV are linear in the mass distribution. And note the trap in the arithmetic: the MM in GMdR3\dfrac{GMd}{R^3} is the mass of the sphere before the hole was cut. Using 7M8\dfrac{7M}{8} there is the standard wrong answer.

Example 6: The tunnel, straight through and off-centre

A tunnel is bored through the Earth, treated as a uniform sphere with RE=6.37×106R_E = 6.37\times10^{6} m and surface g=9.8g = 9.8 m/s2^2. A stone is released at the surface with zero speed. Find (a) the period of the resulting motion, (b) the time to reach the centre, (c) the maximum speed, (d) how all three change if the tunnel is a chord whose closest approach to the centre is RE2\dfrac{R_E}{2}, and (e) why the answer to (a) equals the period of a surface-skimming satellite.

Solution:

  1. (a) Show it is SHM, then read off TT. Inside, only the mass within radius rr acts, giving g(r)=GMERE3r=gRErg(r) = \dfrac{GM_E}{R_E^3}r = \dfrac{g}{R_E}r and hence a=gRErω=gRE=9.86.37×106=1.240×103 rad/sa = -\frac{g}{R_E}\,r \quad\Longrightarrow\quad \omega = \sqrt{\frac{g}{R_E}} = \sqrt{\frac{9.8}{6.37\times10^{6}}} = 1.240\times10^{-3}\ \text{rad/s} T=2πω=5.07×103 s=84.4 minutesT = \frac{2\pi}{\omega} = 5.07\times10^{3}\ \text{s} = 84.4\ \text{minutes}

  2. (b) Quarter of a period. Released at an extreme, the stone reaches the centre — the equilibrium point — after T4=1.27×103 s=21.1 minutes\frac{T}{4} = 1.27\times10^{3}\ \text{s} = 21.1\ \text{minutes}

  3. (c) Maximum speed, at the centre. Amplitude times ω\omega: vmax=ωRE=gRE=(9.8)(6.37×106)=7.90×103 m/sv_{max} = \omega R_E = \sqrt{gR_E} = \sqrt{(9.8)(6.37\times10^{6})} = 7.90\times10^{3}\ \text{m/s}

  4. (d) The chord. Along the tunnel the restoring force is GMEmRE3s-\dfrac{GM_Em}{R_E^3}s, with ss measured from the tunnel's midpoint — the geometric factor cancels. So ω\omega, and therefore T=84.4T = 84.4 min and the 21.1 min to the midpoint, are unchanged. Only the amplitude falls, to the half-chord A=RE2(RE2)2=32REA = \sqrt{R_E^2 - \left(\frac{R_E}{2}\right)^2} = \frac{\sqrt{3}}{2}R_E vmax=ωA=32(7.90×103)=6.84×103 m/sv_{max} = \omega A = \frac{\sqrt{3}}{2}(7.90\times10^{3}) = 6.84\times10^{3}\ \text{m/s}

  5. (e) The coincidence. A satellite skimming the surface needs GMEmRE2=mv2RE\dfrac{GM_Em}{R_E^2} = \dfrac{mv^2}{R_E}, so v=gRE=7.90v = \sqrt{gR_E} = 7.90 km/s and T=2πREv=2πREgT = \dfrac{2\pi R_E}{v} = 2\pi\sqrt{\dfrac{R_E}{g}}the same expression. The reason is that uniform circular motion projected onto a diameter is simple harmonic motion with the same ω\omega; the satellite's shadow on the tunnel and the stone in the tunnel move identically.

Final Answer: (a) 84.4 min; (b) 21.1 min; (c) 7.90 km/s; (d) same period and same time to the midpoint, maximum speed down to 6.84 km/s; (e) circular motion projected on a diameter is SHM with the same ω\omega.

Takeaway: Every straight tunnel through a uniform Earth has the same 84-minute period, through the centre or not, long or short. What the geometry changes is the amplitude, and therefore the speed — never the clock.

Solved Examples (continued)

Example 7: An elliptical orbit, from the two conservation laws

A satellite orbits the Earth in an ellipse with perigee distance 7.00×1067.00\times10^{6} m and apogee distance 1.050×1071.050\times10^{7} m, both measured from the Earth's centre. Take GME=3.982×1014GM_E = 3.982\times10^{14} m3^3/s2^2. Find (a) the semi-major axis and the eccentricity, (b) the speed at perigee and at apogee, (c) the period, and (d) the total energy per kilogram.

Solution:

  1. (a) Geometry first. a=rp+ra2=7.00×106+1.050×1072=8.75×106 ma = \frac{r_p+r_a}{2} = \frac{7.00\times10^{6}+1.050\times10^{7}}{2} = 8.75\times10^{6}\ \text{m} e=rarpra+rp=3.50×1061.750×107=0.200e = \frac{r_a-r_p}{r_a+r_p} = \frac{3.50\times10^{6}}{1.750\times10^{7}} = 0.200

  2. (b) Both laws together. Angular momentum gives vprp=varav_p r_p = v_a r_a, so va=vprprav_a = v_p\dfrac{r_p}{r_a}. Substituting into energy conservation and solving, vp=2GMErarp(rp+ra)=2(3.982×1014)(1.050×107)(7.00×106)(1.750×107)=6.826×107=8262 m/sv_p = \sqrt{\frac{2GM_E\,r_a}{r_p\left(r_p+r_a\right)}} = \sqrt{\frac{2(3.982\times10^{14})(1.050\times10^{7})}{(7.00\times10^{6})(1.750\times10^{7})}} = \sqrt{6.826\times10^{7}} = 8262\ \text{m/s} va=vprpra=8262×7.0010.50=5508 m/sv_a = v_p\frac{r_p}{r_a} = 8262\times\frac{7.00}{10.50} = 5508\ \text{m/s} Cross-check with vis-viva at perigee: v2=GME(27.00×10618.75×106)=6.826×107v^2 = GM_E\left(\dfrac{2}{7.00\times10^{6}} - \dfrac{1}{8.75\times10^{6}}\right) = 6.826\times10^{7}. Agreed. A second check: vpva=(8262)(5508)=4.551×107v_pv_a = (8262)(5508) = 4.551\times10^{7}, and GMEa=4.551×107\dfrac{GM_E}{a} = 4.551\times10^{7}. Agreed again.

  3. (c) Period, from aa alone. T=2πa3GME=2π(8.75×106)33.982×1014=2π(1297)=8.15×103 s=136 minT = 2\pi\sqrt{\frac{a^3}{GM_E}} = 2\pi\sqrt{\frac{(8.75\times10^{6})^3}{3.982\times10^{14}}} = 2\pi(1297) = 8.15\times10^{3}\ \text{s} = 136\ \text{min}

  4. (d) Total energy per kilogram. Em=GME2a=3.982×10141.750×107=2.275×107 J/kg\frac{E}{m} = -\frac{GM_E}{2a} = -\frac{3.982\times10^{14}}{1.750\times10^{7}} = -2.275\times10^{7}\ \text{J/kg} Check it directly at perigee: 12(8262)23.982×10147.00×106=3.413×1075.689×107=2.276×107\tfrac{1}{2}(8262)^2 - \dfrac{3.982\times10^{14}}{7.00\times10^{6}} = 3.413\times10^{7} - 5.689\times10^{7} = -2.276\times10^{7}. Agreed.

Final Answer: (a) a=8.75×106a = 8.75\times10^{6} m, e=0.200e = 0.200; (b) 8262 m/s and 5508 m/s; (c) 8.15×1038.15\times10^{3} s, about 136 min; (d) 2.275×107-2.275\times10^{7} J/kg.

Takeaway: Two laws, two unknowns, and then three independent checks — vis-viva, vpva=GMav_pv_a = \dfrac{GM}{a}, and the direct energy sum. In an ellipse question you should never submit an answer you have not cross-checked; the checks cost fifteen seconds each.

Example 8: Vis-viva, used four times

A spacecraft is in an elliptical orbit about the Earth with a=8.75×106a = 8.75\times10^{6} m and e=0.200e = 0.200, so rp=7.00×106r_p = 7.00\times10^{6} m and ra=1.050×107r_a = 1.050\times10^{7} m. Using v2=GM(2r1a)v^2 = GM\left(\dfrac{2}{r}-\dfrac{1}{a}\right) with GME=3.982×1014GM_E = 3.982\times10^{14}, find its speed (a) at r=8.00×106r = 8.00\times10^{6} m, (b) at the two ends of the minor axis, where r=ar = a, (c) what its speed would be in a circular orbit of radius aa, and (d) the extra speed it would need at perigee to escape the Earth entirely.

Solution:

  1. (a) One substitution. v2=3.982×1014(28.00×10618.75×106)=3.982×1014(2.500×1071.143×107)v^2 = 3.982\times10^{14}\left(\frac{2}{8.00\times10^{6}} - \frac{1}{8.75\times10^{6}}\right) = 3.982\times10^{14}\left(2.500\times10^{-7} - 1.143\times10^{-7}\right) v2=5.404×107v=7351 m/sv^2 = 5.404\times10^{7} \quad\Longrightarrow\quad v = 7351\ \text{m/s}

  2. (b) At r=ar = a. Then 2r1a=2a1a=1a\dfrac{2}{r}-\dfrac{1}{a} = \dfrac{2}{a}-\dfrac{1}{a} = \dfrac{1}{a}: v=GMEa=3.982×10148.75×106=6746 m/sv = \sqrt{\frac{GM_E}{a}} = \sqrt{\frac{3.982\times10^{14}}{8.75\times10^{6}}} = 6746\ \text{m/s}

  3. (c) The circular comparison. A circular orbit of radius aa has v=GMEa=6746v = \sqrt{\dfrac{GM_E}{a}} = 6746 m/s — exactly the same number. But the spacecraft is not moving in a circle at that instant: on the ellipse its velocity is not perpendicular to the radius, so the same speed is being spent partly on moving outwards.

  4. (d) Escape from perigee. Escape needs E=0E = 0, that is vesc2=2GMErpv_{esc}^2 = \dfrac{2GM_E}{r_p}: vesc=2(3.982×1014)7.00×106=1.138×108=10666 m/sv_{esc} = \sqrt{\frac{2(3.982\times10^{14})}{7.00\times10^{6}}} = \sqrt{1.138\times10^{8}} = 10\,666\ \text{m/s} The satellite already has 8262 m/s there, so it needs Δv=106668262=2404 m/s\Delta v = 10\,666 - 8262 = 2404\ \text{m/s}

Final Answer: (a) 7351 m/s; (b) 6746 m/s; (c) the same 6746 m/s; (d) an extra 2404 m/s at perigee.

Takeaway: Vis-viva replaces both conservation laws with one substitution. Learn the four special cases — a=ra=r gives the circle, aa\to\infty gives escape, r=rpr=r_p and r=rar=r_a give the extremes — and most elliptical-orbit questions become a single line of arithmetic. And note (c): equal speed does not mean equal orbit.

Example 9: A Hohmann transfer to geostationary orbit

A satellite is in a circular orbit of radius r1=7.00×106r_1 = 7.00\times10^{6} m and must reach a circular orbit of radius r2=4.22×107r_2 = 4.22\times10^{7} m, the geostationary radius. Using a Hohmann transfer with GME=3.982×1014GM_E = 3.982\times10^{14}, find (a) the semi-major axis of the transfer ellipse, (b) both burns, (c) the total velocity change, and (d) the coast time. (e) Comment on the fact that both burns speed the satellite up.

Solution:

  1. (a) The transfer ellipse touches both circles. at=r1+r22=7.00×106+4.22×1072=2.46×107 ma_t = \frac{r_1+r_2}{2} = \frac{7.00\times10^{6}+4.22\times10^{7}}{2} = 2.46\times10^{7}\ \text{m}

  2. (b) Burn 1, at perigee. v1=GMEr1=3.982×10147.00×106=7542 m/sv_1 = \sqrt{\frac{GM_E}{r_1}} = \sqrt{\frac{3.982\times10^{14}}{7.00\times10^{6}}} = 7542\ \text{m/s} vp=GME(2r11at)=3.982×1014(2.857×1074.065×108)=9878 m/sv_p = \sqrt{GM_E\left(\frac{2}{r_1}-\frac{1}{a_t}\right)} = \sqrt{3.982\times10^{14}\left(2.857\times10^{-7}-4.065\times10^{-8}\right)} = 9878\ \text{m/s} Δv1=98787542=2336 m/s\Delta v_1 = 9878 - 7542 = 2336\ \text{m/s} Burn 2, at apogee. va=GME(2r21at)=3.982×1014(4.739×1084.065×108)=1639 m/sv_a = \sqrt{GM_E\left(\frac{2}{r_2}-\frac{1}{a_t}\right)} = \sqrt{3.982\times10^{14}\left(4.739\times10^{-8}-4.065\times10^{-8}\right)} = 1639\ \text{m/s} v2=GMEr2=3072 m/sΔv2=30721639=1433 m/sv_2 = \sqrt{\frac{GM_E}{r_2}} = 3072\ \text{m/s} \quad\Longrightarrow\quad \Delta v_2 = 3072 - 1639 = 1433\ \text{m/s}

  3. (c) Total. Δvtotal=2336+1433=3769 m/s3.77 km/s\Delta v_{total} = 2336 + 1433 = 3769\ \text{m/s} \approx 3.77\ \text{km/s}

  4. (d) Half an ellipse. t=πat3GME=π(2.46×107)33.982×1014=π(6115)=1.92×104 s=5.34 ht = \pi\sqrt{\frac{a_t^3}{GM_E}} = \pi\sqrt{\frac{(2.46\times10^{7})^3}{3.982\times10^{14}}} = \pi(6115) = 1.92\times10^{4}\ \text{s} = 5.34\ \text{h}

  5. (e) Why two prograde burns leave you slower. The satellite starts at 7542 m/s and finishes at 3072 m/s, having been accelerated twice. There is no contradiction: the engine raised the total energy from GMEm2r1-\dfrac{GM_Em}{2r_1} to GMEm2r2-\dfrac{GM_Em}{2r_2}, which is less negative, and in a higher orbit the required speed is smaller. Energy went up; speed went down. The difference went into potential energy, and then some.

Final Answer: (a) 2.46×1072.46\times10^{7} m; (b) Δv1=2336\Delta v_1 = 2336 m/s, Δv2=1433\Delta v_2 = 1433 m/s; (c) 3769 m/s; (d) 1.92×1041.92\times10^{4} s, about 5.34 hours; (e) both burns add energy, and a higher orbit is a slower orbit.

Takeaway: The Hohmann recipe never changes: at=r1+r22a_t = \dfrac{r_1+r_2}{2}, one vis-viva at each end, subtract, and the coast time is half the transfer period. Memorise that sequence rather than the closed-form Δv\Delta v expressions, which are easy to mis-copy under pressure.

Solved Examples (continued)

Example 10: A binary star system

Two stars of masses 3.0×10303.0\times10^{30} kg and 6.0×10306.0\times10^{30} kg move in circular orbits about their common centre of mass, separated by 2.0×10112.0\times10^{11} m. Take G=6.67×1011G = 6.67\times10^{-11}. Find (a) the distance of each star from the centre of mass, (b) the angular speed and the period, (c) each star's speed, (d) the reduced mass, and (e) the total energy of the system.

Solution:

  1. (a) The centre of mass splits the separation in the inverse mass ratio. r1=m2dm1+m2=(6.0×1030)(2.0×1011)9.0×1030=1.33×1011 mr_1 = \frac{m_2 d}{m_1+m_2} = \frac{(6.0\times10^{30})(2.0\times10^{11})}{9.0\times10^{30}} = 1.33\times10^{11}\ \text{m} r2=m1dm1+m2=(3.0×1030)(2.0×1011)9.0×1030=6.67×1010 mr_2 = \frac{m_1 d}{m_1+m_2} = \frac{(3.0\times10^{30})(2.0\times10^{11})}{9.0\times10^{30}} = 6.67\times10^{10}\ \text{m} The lighter star travels on the bigger circle. Their sum is 2.0×10112.0\times10^{11} m, as it must be.

  2. (b) One shared ω\omega. ω2=G(m1+m2)d3=(6.67×1011)(9.0×1030)(2.0×1011)3=6.003×10208.0×1033=7.504×1014\omega^2 = \frac{G(m_1+m_2)}{d^3} = \frac{(6.67\times10^{-11})(9.0\times10^{30})}{(2.0\times10^{11})^3} = \frac{6.003\times10^{20}}{8.0\times10^{33}} = 7.504\times10^{-14} ω=2.74×107 rad/s,T=2πω=2.29×107 s=265 days\omega = 2.74\times10^{-7}\ \text{rad/s}, \qquad T = \frac{2\pi}{\omega} = 2.29\times10^{7}\ \text{s} = 265\ \text{days}

  3. (c) Speeds. v1=ωr1=(2.74×107)(1.33×1011)=3.65×104 m/s=36.5 km/sv_1 = \omega r_1 = (2.74\times10^{-7})(1.33\times10^{11}) = 3.65\times10^{4}\ \text{m/s} = 36.5\ \text{km/s} v2=ωr2=1.83×104 m/s=18.3 km/sv_2 = \omega r_2 = 1.83\times10^{4}\ \text{m/s} = 18.3\ \text{km/s} Note v1v2=2=m2m1\dfrac{v_1}{v_2} = 2 = \dfrac{m_2}{m_1}: the heavier star crawls, the lighter one races. Their momenta are equal and opposite, so the centre of mass never moves.

  4. (d) Reduced mass. μ=m1m2m1+m2=(3.0×1030)(6.0×1030)9.0×1030=2.0×1030 kg\mu = \frac{m_1m_2}{m_1+m_2} = \frac{(3.0\times10^{30})(6.0\times10^{30})}{9.0\times10^{30}} = 2.0\times10^{30}\ \text{kg}

  5. (e) Total energy. K=12m1v12+12m2v22=2.00×1039+1.00×1039=3.00×1039 JK = \tfrac{1}{2}m_1v_1^2 + \tfrac{1}{2}m_2v_2^2 = 2.00\times10^{39} + 1.00\times10^{39} = 3.00\times10^{39}\ \text{J} U=Gm1m2d=(6.67×1011)(1.8×1061)2.0×1011=6.00×1039 JU = -\frac{Gm_1m_2}{d} = -\frac{(6.67\times10^{-11})(1.8\times10^{61})}{2.0\times10^{11}} = -6.00\times10^{39}\ \text{J} E=K+U=3.00×1039 JE = K + U = -3.00\times10^{39}\ \text{J} which matches the general result E=Gm1m22dE = -\dfrac{Gm_1m_2}{2d} for a circular relative orbit, and shows the familiar K:U:E=1:2:1K : U : E = 1 : -2 : -1 pattern surviving intact.

Final Answer: (a) 1.33×10111.33\times10^{11} m and 6.67×10106.67\times10^{10} m; (b) 2.74×1072.74\times10^{-7} rad/s, 2.29×1072.29\times10^{7} s (265 days); (c) 36.5 km/s and 18.3 km/s; (d) 2.0×10302.0\times10^{30} kg; (e) 3.00×1039-3.00\times10^{39} J.

Takeaway: In a binary, one ω\omega, one TT, two radii, two speeds — and the heavier star is always the slower one on the smaller circle. Kepler's third law still holds, with the total mass in the denominator, which is exactly how astronomers weigh stars they have never visited.

Example 11: Taking the Earth apart

(a) Derive the gravitational self-energy of a uniform sphere of mass MM and radius RR. (b) Evaluate it for the Earth, with ME=5.97×1024M_E = 5.97\times10^{24} kg and RE=6.37×106R_E = 6.37\times10^{6} m. (c) How long would the Sun, radiating 3.8×10263.8\times10^{26} W, take to supply that much energy? (d) Compare with the energy needed to lift a single 1 kg stone from the surface to infinity.

Solution:

  1. (a) Assemble it shell by shell. Suppose the sphere has already been built out to radius rr. Its mass so far is m=Mr3R3m = M\frac{r^3}{R^3} and the next shell, of thickness drdr, carries dm=ρ4πr2dr=3Mr2R3drdm = \rho\,4\pi r^2\,dr = \frac{3M r^2}{R^3}\,dr Bringing that shell in from infinity onto a body of mass mm changes the potential energy by dU=Gmdmr=GrMr3R33Mr2R3dr=3GM2R6r4drdU = -\frac{G\,m\,dm}{r} = -\frac{G}{r}\cdot M\frac{r^3}{R^3}\cdot\frac{3Mr^2}{R^3}dr = -\frac{3GM^2}{R^6}\,r^4\,dr Uself=3GM2R60Rr4dr=3GM2R6R55=35GM2RU_{self} = -\frac{3GM^2}{R^6}\int_0^R r^4\,dr = -\frac{3GM^2}{R^6}\cdot\frac{R^5}{5} = -\frac{3}{5}\frac{GM^2}{R}

  2. (b) The Earth's figure. Uself=3(6.67×1011)(5.97×1024)25(6.37×106)=7.132×10393.185×107=2.24×1032 JU_{self} = -\frac{3(6.67\times10^{-11})(5.97\times10^{24})^2}{5(6.37\times10^{6})} = -\frac{7.132\times10^{39}}{3.185\times10^{7}} = -2.24\times10^{32}\ \text{J} Dispersing the Earth to infinity therefore costs +2.24×1032+2.24\times10^{32} J.

  3. (c) In units of sunshine. t=2.24×10323.8×1026=5.9×105 s=6.8 dayst = \frac{2.24\times10^{32}}{3.8\times10^{26}} = 5.9\times10^{5}\ \text{s} = 6.8\ \text{days}

  4. (d) Against one stone. Lifting 1 kg from the surface to infinity costs GME(1)RE=6.25×107 J\frac{GM_E(1)}{R_E} = 6.25\times10^{7}\ \text{J} The ratio is 2.24×10326.25×107=3.6×1024\dfrac{2.24\times10^{32}}{6.25\times10^{7}} = 3.6\times10^{24} — comparable to, and in fact 35\tfrac{3}{5} of, the number of kilograms in the Earth. That is exactly what the factor 35\tfrac{3}{5} is telling you: on average, each kilogram of the Earth is bound with 35\tfrac{3}{5} of the energy that binds a kilogram sitting on the surface.

Final Answer: (a) Uself=3GM25RU_{self} = -\dfrac{3GM^2}{5R}; (b) 2.24×1032-2.24\times10^{32} J; (c) about 5.9×1055.9\times10^{5} s, roughly 6.8 days of the Sun's total output; (d) 3.6×10243.6\times10^{24} times the cost of freeing one surface kilogram.

Takeaway: Self-energy is about building a body, not moving one, and it always carries M2M^2. Contrast the three: GMEmr-\dfrac{GM_Em}{r} for a body in a field, GMEm2r-\dfrac{GM_Em}{2r} for the total energy of an orbit, 3GME25RE-\dfrac{3GM_E^2}{5R_E} for the Earth's own binding. A thin shell gives GM22R-\dfrac{GM^2}{2R} by the same argument, and 12<35\tfrac{1}{2} < \tfrac{3}{5} because a shell's mass sits further from itself.

Example 12: Catching a satellite, and one composite body

Part A. Two satellites are in the same circular orbit of radius 7.50×1067.50\times10^{6} m around the Earth, with the chaser 15°15° behind the target. The chaser will close the gap using a phasing orbit over 3 revolutions. Take GME=3.982×1014GM_E = 3.982\times10^{14}. Find (a) the original period, (b) the period and semi-major axis of the phasing orbit, (c) its perigee distance, and (d) the total velocity change required.

Part B. A uniform solid sphere of mass MM and radius RR sits at the centre of a concentric thin shell of mass 2M2M and radius 3R3R. Find gg and VV (e) at r=2Rr = 2R and (f) at the centre.

Solution:

  1. (a) The circular period. T=2πr3GME=2π(7.50×106)33.982×1014=2π(1029)=6467 s=108 minT = 2\pi\sqrt{\frac{r^3}{GM_E}} = 2\pi\sqrt{\frac{(7.50\times10^{6})^3}{3.982\times10^{14}}} = 2\pi(1029) = 6467\ \text{s} = 108\ \text{min}

  2. (b) Which way to go — and by how much. To catch a target ahead of it, the chaser needs a shorter period, so it must drop into a lower orbit and therefore fire backwards. Over 3 laps it must gain 15°15°, so each lap must be short by 5°: ΔT=5360(6467)=89.8 sT=646790=6377 s\Delta T = \frac{5}{360}(6467) = 89.8\ \text{s} \quad\Longrightarrow\quad T^{\prime} = 6467 - 90 = 6377\ \text{s} a=(GMET24π2)1/3=((3.982×1014)(4.067×107)39.48)1/3=7.430×106 ma^{\prime} = \left(\frac{GM_E\,T^{\prime 2}}{4\pi^2}\right)^{1/3} = \left(\frac{(3.982\times10^{14})(4.067\times10^{7})}{39.48}\right)^{1/3} = 7.430\times10^{6}\ \text{m}

  3. (c) The perigee. The burn happens at the original radius, which becomes the apogee of the phasing ellipse, so rp=2ar=2(7.430×106)7.50×106=7.361×106 mr_p = 2a^{\prime} - r = 2(7.430\times10^{6}) - 7.50\times10^{6} = 7.361\times10^{6}\ \text{m} That is 991 km above the surface — comfortably clear of the Earth, which matters: too big a phase angle in too few laps drives the perigee underground and the manoeuvre becomes impossible.

  4. (d) Two burns, equal and opposite. vcirc=GMEr=7287 m/s,vapo=GME(2r1a)=7252 m/sv_{circ} = \sqrt{\frac{GM_E}{r}} = 7287\ \text{m/s}, \qquad v_{apo} = \sqrt{GM_E\left(\frac{2}{r}-\frac{1}{a^{\prime}}\right)} = 7252\ \text{m/s} Δv=72877252=34.2 m/s per burn,total=68.4 m/s\lvert\Delta v\rvert = 7287-7252 = 34.2\ \text{m/s per burn}, \qquad \text{total} = 68.4\ \text{m/s} The first burn slows the chaser onto the ellipse; three laps later, alongside the target, an identical forward burn puts it back on the circle.

  5. (e) At r=2Rr = 2R: outside the sphere, inside the shell. g=GM(2R)2+0=GM4R2,V=GM2RG(2M)3R=7GM6Rg = \frac{GM}{(2R)^2} + 0 = \frac{GM}{4R^2}, \qquad V = -\frac{GM}{2R} - \frac{G(2M)}{3R} = -\frac{7GM}{6R}

  6. (f) At the centre: both fields vanish, neither potential does. g=0,V=3GM2R2GM3R=13GM6Rg = 0, \qquad V = -\frac{3GM}{2R} - \frac{2GM}{3R} = -\frac{13GM}{6R}

Final Answer: (a) 6467 s, about 108 min; (b) T=6377T^{\prime} = 6377 s, a=7.430×106a^{\prime} = 7.430\times10^{6} m; (c) 7.361×1067.361\times10^{6} m, i.e. 991 km altitude; (d) 34.2 m/s per burn, 68.4 m/s in total; (e) GM4R2\dfrac{GM}{4R^2} and 7GM6R-\dfrac{7GM}{6R}; (f) 00 and 13GM6R-\dfrac{13GM}{6R}.

Takeaway: To catch up, slow down. It is the most counter-intuitive sentence in orbital mechanics and it is exactly right. And in Part B, notice that the shell contributes nothing to gg inside it and a great deal to VV — zero field never means zero potential.