Case Study 1 — Threading a Satellite Through the Sky

DataField.Dev

Case Study 1 — Threading a Satellite Through the Sky

Field: Aerospace engineering / astrodynamics Calculus used: Vector-valued functions as trajectories (§28.4), integrating acceleration to recover motion (§28.4), arc length (§28.5), and the planarity of orbits from the cross-product rule (§28.11).

A satellite is a vector-valued function

Mei Tanaka is a flight dynamics engineer for an Earth-imaging company. Her satellite — call it Heron-3 — circles the planet roughly fifteen times a day at an altitude near 600 km, photographing farmland for crop-health analytics. To Mei, Heron-3 is not a piece of hardware floating in the abstract; it is a vector-valued function. Its position at every instant is

$$\mathbf{r}(t) = \langle x(t),\, y(t),\, z(t)\rangle,$$

a single point tracing a single curve through three-dimensional space — exactly the object Chapter 28 opens with. Everything her team does, from pointing the camera to dodging space debris, is reading derivatives and integrals of this one function. This is the chapter's threshold idea made operational: one input (time), one vector output, and all of single-variable calculus carried over component-wise.

The governing law is Newtonian gravity, written as a second-order vector ODE (§28.11):

$$\mathbf{r}''(t) = -\frac{GM}{|\mathbf{r}(t)|^3}\,\mathbf{r}(t).$$

The acceleration points back toward Earth's center (the $-\mathbf{r}$ direction) with magnitude $GM/|\mathbf{r}|^2$, the inverse-square law. With Earth's gravitational parameter $\mu = GM \approx 398{,}600\ \text{km}^3/\text{s}^2$ and an orbital radius $r = 6378 + 600 \approx 6978\ \text{km}$, Mei first checks the circular-orbit speed. Setting the gravitational pull equal to the centripetal requirement $v^2/r = \mu/r^2$ gives

$$v = \sqrt{\frac{\mu}{r}} = \sqrt{\frac{398{,}600}{6978}} \approx \sqrt{57.12} \approx 7.56\ \text{km/s}.$$

That is the number every textbook quotes as "about 7.5 km/s for low Earth orbit," and here it falls straight out of equating two magnitudes. The orbital period follows from the circumference over the speed:

$$T = \frac{2\pi r}{v} = \frac{2\pi(6978)}{7.56} \approx \frac{43{,}840}{7.56} \approx 5800\ \text{s} \approx 97\ \text{minutes}.$$

Fifteen orbits a day, as advertised. Mei has not solved the ODE yet — she has only read off the magnitudes of velocity and acceleration. That is the power of the vector viewpoint: the geometry yields physical numbers before any integration.

Why the orbit is flat — a one-line proof

Heron-3's customers care about coverage: which strips of ground the camera can see. That depends on a structural fact about orbits that the chapter proves with nothing more than the cross-product rule (§28.3, §28.11). Define the angular momentum per unit mass

$$\mathbf{L}(t) = \mathbf{r}(t) \times \mathbf{r}'(t).$$

Differentiate using the cross-product rule, keeping the factor order intact (the §28.3 pitfall):

$$\mathbf{L}'(t) = \mathbf{r}' \times \mathbf{r}' + \mathbf{r} \times \mathbf{r}''.$$

The first term is $\mathbf{0}$ because any vector crossed with itself vanishes. In the second, substitute the gravity law $\mathbf{r}'' = -\frac{\mu}{|\mathbf{r}|^3}\mathbf{r}$:

$$\mathbf{r} \times \mathbf{r}'' = \mathbf{r} \times \Big(\!-\frac{\mu}{|\mathbf{r}|^3}\mathbf{r}\Big) = -\frac{\mu}{|\mathbf{r}|^3}\,(\mathbf{r}\times\mathbf{r}) = \mathbf{0}.$$

So $\mathbf{L}'(t) = \mathbf{0}$: angular momentum is a constant vector. Because $\mathbf{r}$ and $\mathbf{r}'$ are always perpendicular to the fixed vector $\mathbf{L}$, the whole trajectory lies in a single fixed plane through Earth's center. Orbits are planar — and Mei's team just proved it from the product rule, not from a chart. The constancy of $|\mathbf{L}|$ is also exactly Kepler's Second Law (equal areas in equal times), since the rate at which the position sweeps out area is $\tfrac12|\mathbf{r}\times\mathbf{r}'| = \tfrac12|\mathbf{L}|$. And the shape of the path inside that plane — Heron-3's nearly circular orbit, or the long ellipse of a transfer trajectory — is one of the conic sections from Chapter 27, emerging when you actually solve the ODE. Part VI quietly hands the baton back to Part V.

Integrating a thruster burn

Heron-3's orbit slowly decays: even at 600 km, the thin upper atmosphere exerts drag, lowering the altitude a little each week. Left alone, the satellite would re-enter within a few years. So once a month Mei commands a short prograde burn — a brief thrust along the direction of motion — to nudge the velocity back up. Designing that burn is a direct application of "integrate acceleration twice to recover the trajectory" (§28.4).

During the burn, the thruster adds a constant acceleration $\mathbf{a}_{\text{thrust}}$ on top of gravity. Over the few seconds of firing, gravity barely changes the velocity, so Mei models the burn locally as constant acceleration. Suppose, in a frame where $x$ runs along the current velocity, the thruster delivers

$$\mathbf{a}(t) = \langle 0.40,\; 0,\; 0\rangle\ \text{m/s}^2,$$

with the satellite entering the burn at $\mathbf{v}(0) = \langle 7560,\, 0,\, 0\rangle$ m/s and (in this local frame) $\mathbf{r}(0) = \langle 0,0,0\rangle$. Integrate once for velocity (§28.4):

$$\mathbf{v}(t) = \int \mathbf{a}(t)\,dt = \langle 0.40\,t,\; 0,\; 0\rangle + \mathbf{C}_1,$$

and the initial condition fixes $\mathbf{C}_1 = \langle 7560, 0, 0\rangle$, so $\mathbf{v}(t) = \langle 7560 + 0.40\,t,\; 0,\; 0\rangle$. Integrate again for position:

$$\mathbf{r}(t) = \langle 7560\,t + 0.20\,t^2,\; 0,\; 0\rangle + \mathbf{C}_2, \qquad \mathbf{C}_2 = \mathbf{0}.$$

For a 30-second burn, the velocity gain is the delta-v

$$\Delta v = 0.40 \times 30 = 12\ \text{m/s},$$

and the distance covered during the burn is $7560(30) + 0.20(30)^2 = 226{,}800 + 180 = 226{,}980$ m, about 227 km of arc — virtually all of it the coasting term $7560t$, with the thrust contributing only the tiny $0.20t^2$ correction. This is the constant of integration doing real work: $\mathbf{C}_1$ encodes "where the velocity started," and the new velocity is the old velocity plus the accumulated thrust. An inertial navigation unit on board performs this same double integration in reverse — reading accelerometers and integrating to know where Heron-3 is without any GPS fix, the Net Change Theorem of Chapter 14 running in real time on a vector.

How far has it traveled?

For mission accounting, Mei sometimes needs the distance flown, not just the displacement — the arc length of the trajectory (§28.5). Over one orbit of the (very nearly circular) path, the speed is essentially the constant $v \approx 7.56$ km/s, so

$$L = \int_0^{T} |\mathbf{r}'(t)|\,dt \approx \int_0^{5800} 7.56\,dt \approx 43{,}800\ \text{km},$$

which matches the circumference $2\pi r = 2\pi(6978) \approx 43{,}840$ km. The small discrepancy is exactly the rounding in the speed. The lesson the chapter stresses is that arc length is accumulated speed, and that it depends on the curve, not on how fast you drive it — re-clocking the same orbit at twice the parameter rate would halve the time interval and double the integrand, leaving $L$ untouched.

The harder reality, and the numerical answer

The clean two-body picture is only the first draft. Earth is not a perfect sphere; its equatorial bulge perturbs $\mathbf{r}''$, slowly rotating the orbital plane — the so-called $J_2$ precession that Mei actually exploits to keep Heron-3's lighting consistent for imaging. Atmospheric drag, the Moon, the Sun, and even sunlight pressure each add their own term to the acceleration vector. None of these has a closed-form solution. So in practice the trajectory is computed exactly as §28.10 describes: feed the full perturbed acceleration into a numerical integrator and step the vector ODE forward thousands of times per orbit. The pseudocode mirrors the chapter's np.gradient workflow, integration replacing differentiation:

# Sketch: integrate the two-body vector ODE r'' = -mu r / |r|^3 forward in time.
# (Hand-checked outputs shown as comments; not executed here.)
mu = 398600.0          # km^3/s^2
r  = [6978.0, 0.0, 0.0]      # initial position (km)
v  = [0.0, 7.56, 0.0]        # initial velocity (km/s), perpendicular -> near-circular
# After one period (~5800 s) the integrator returns r back near its start:
# r ~ [6978, 0, 0], having traced an arc of length ~43,840 km.

Every crop-health image Heron-3 returns rests on this chain of reasoning: a satellite is a vector-valued function; its acceleration is a vector law; integrating that law gives the path; the cross-product rule guarantees the path is planar; and arc length measures the journey. The mathematics of Chapter 28 is not a metaphor for the mission — it is the mission.

Discussion Questions

Mei found the orbital speed by setting $v^2/r = \mu/r^2$ without solving any differential equation. What did she trade away by working with magnitudes only, and what could she not learn this way (think: shape and orientation of the orbit)?
The proof that orbits are planar used only $\mathbf{r}\times\mathbf{r} = \mathbf{0}$ and the cross-product rule. Where exactly would the argument fail if gravity did not point along $-\mathbf{r}$ (a so-called non-central force)?
During the 30-second burn, the thrust contributed only 180 m out of nearly 227 km of arc. Why is the quadratic term $0.20t^2$ so small compared with the linear term, and what does that say about the relative sizes of velocity and acceleration here?
Explain why doubling the parameter rate (driving the same orbit "twice as fast" in $t$) leaves the arc length unchanged. Which §28.5 pitfall does this guard against?
The real orbit needs numerical integration because of perturbations. Connect this to the §28.6 computational note: why can so few trajectories be reparametrized by arc length in closed form?

A Short Annotated Reading

Curtis, H. (2020). Orbital Mechanics for Engineering Students (4th ed.), Ch. 2–3. Derives the two-body equation $\mathbf{r}'' = -\mu\mathbf{r}/r^3$ and the conservation of $\mathbf{L}$ in exactly the vector language of this chapter; the best next step if the planarity proof intrigued you.
Bate, Mueller & White (1971). Fundamentals of Astrodynamics, Dover. A famously inexpensive classic; Chapter 1 builds orbits from vector calculus alone, with the angular-momentum argument front and center.
Vallado, D. (2013). Fundamentals of Astrodynamics and Applications (4th ed.). The professional reference; see its treatment of $J_2$ perturbation and numerical propagation to see where the clean two-body model gives way to the integrators of §28.10.