Space environment — LEO radiation (SAA, Van Allen), atomic oxygen, MMOD debris: Exercises

These problems build from "can you recall the picture" up to "can you design a mission trade." Every solution is hidden inside a collapsible callout — try first, then reveal.

Definition Symbols used on this page (defined once, here)

Every symbol below is spelled out so this page stands on its own:

$R_{\oplus} = 6371 km$ — Earth's radius (surface to centre).
$μ = G M_{\oplus} = 398600 km^{3} / s^{2}$ — Earth's gravitational parameter ( $G$ = gravitational constant, $M_{\oplus}$ = Earth's mass).
$m$ — a particle's mass (kg); $m_{p} = 1.66 \times 1 0^{- 27} kg$ a proton, $m_{O} = 16 m_{p}$ an oxygen atom.
$q = 1.6 \times 1 0^{- 19} C$ — a particle's electric charge (Coulombs).
$v$ — velocity; $v_{⊥}$ = the part of velocity perpendicular to the magnetic field.
$B$ — magnetic field strength (Tesla, T).
$r_{L}$ — gyroradius (the radius of the little circle a charged particle spins in).
$Φ$ — flux: how many atoms hit each cm² each second (atoms/cm²/s).
$F$ — fluence: flux multiplied by time = total atoms per cm² over the mission.
$Y$ — erosion yield: volume of material removed per incoming atom (cm³/atom).
$TID$ — total ionizing dose: radiation energy absorbed per kg, measured in rad.
$λ$ — expected (average) number of impacts in a Poisson process.

Level 1 — Recognition (recall the facts and pictures)

Exercise 1.1 (L1). Name the three nested motions a charged particle performs when trapped in the Van Allen belts, in order of increasing spatial scale. Use the figure below to check your picture.

Figure — Space environment — LEO radiation (SAA, Van Allen), atomic oxygen, MMOD debris

Recall Solution 1.1

Read the figure from the tightest loop outward:

Gyration — the red spiral in the figure: tight circles around a single field line. Radius $r_{L} = m v_{⊥} / q B$ (here $m$ = mass, $v_{⊥}$ = the sideways speed, $q$ = charge, $B$ = field strength) — metres to kilometres across.
Bounce — the black arrows sliding along the field line between the two black pole-dots: the particle runs back and forth pole-to-pole (the "magnetic bottle" mirror effect).
Drift — the dashed ellipse at the bottom: the whole spiral migrates slowly around the Earth in longitude (electrons east-to-west, protons west-to-east) because $B$ is stronger nearer Earth.

Notice in the figure how the red gyration circle is tiny compared to the dashed drift loop — that ordering (gyration < bounce < drift) is exactly what the question asks for. Together these three keep a particle trapped for months to years.

Exercise 1.2 (L1). Match each threat to the altitude band where it is worst: (a) atomic oxygen, (b) inner-belt protons dipping into LEO, (c) relativistic outer-belt electrons.

Recall Solution 1.2

(a) Atomic oxygen — 200–700 km (dense enough residual atmosphere, UV photodissociates $O_{2}$ ).
(b) Inner-belt protons in LEO — worst in the South Atlantic Anomaly (SAA), where the belt sags to ~200 km.
(c) Outer-belt electrons — ~13,000–60,000 km, well above LEO (a problem for GEO/MEO, e.g. the Van Allen Probes mission).

Exercise 1.3 (L1). True/False: "The SAA exists because the Sun heats the South Atlantic more." Correct the statement if false.

Recall Solution 1.3

False. The SAA has nothing to do with surface heating. It exists because Earth's magnetic dipole is offset (~500 km toward the western Pacific) and tilted relative to the geographic centre. Over the South Atlantic the field is therefore weakest, so the inner belt's lower edge sags to LEO altitudes and dumps proton flux onto passing satellites.

Level 2 — Application (plug into one formula correctly)

Exercise 2.1 (L2). A proton with $v_{⊥} = 2 \times 1 0^{7} m/s$ moves in a field $B = 3 \times 1 0^{- 5} T$ . Compute its gyroradius $r_{L}$ . ( $m_{p} = 1.66 \times 1 0^{- 27} kg$ , $q = 1.6 \times 1 0^{- 19} C$ .)

Recall Solution 2.1

Why this formula? The Lorentz force $F = q (v \times B)$ is always perpendicular to velocity, so it can never change speed — only turn the particle. A constant sideways-turning force produces a circle. Balancing it against the centripetal requirement $q v_{⊥} B = m v_{⊥}^{2} / r_{L}$ and cancelling one $v_{⊥}$ gives: $r_{L} = \frac{m v _{⊥}}{q B}$ where $m$ = mass, $v_{⊥}$ = perpendicular speed, $q$ = charge, $B$ = field strength. $r_{L} = \frac{( 1.66 \times 1 0 ^{- 27} ) ( 2 \times 1 0 ^{7} )}{( 1.6 \times 1 0 ^{- 19} ) ( 3 \times 1 0 ^{- 5} )} = 6.92 \times 1 0^{3} m \approx 6.9 km .$ The gyroradius is kilometres, tiny compared to the belt's tens-of-thousands-of-km extent — that is why "spiralling along a field line" is a good picture.

Exercise 2.2 (L2). Atomic oxygen strikes a spacecraft at relative speed $v = 7.7 km/s$ . Compute the collision energy per O atom in eV. ( $m_{O} = 16 \times 1.66 \times 1 0^{- 27} kg$ , $1 eV = 1.6 \times 1 0^{- 19} J$ .) Then state whether it can break a C–C bond (3.6 eV).

Recall Solution 2.2

Why kinetic energy? The damage is chemical bond-breaking, and a bond breaks if the impact energy delivered exceeds the bond energy. Speed is fixed by orbital velocity, so: $E = \frac{1}{2} m_{O} v^{2} = \frac{1}{2} (16 \times 1.66 \times 1 0^{- 27}) (7700)^{2} = 7.87 \times 1 0^{- 19} J .$ Convert: $E = 7.87 \times 1 0^{- 19} /1.6 \times 1 0^{- 19} = 4.92 eV$ . Since $4.92 eV > 3.6 eV$ , yes — each ram-facing O atom carries enough energy to snap a C–C bond. That is the whole reason bare Kapton erodes.

Exercise 2.3 (L2). A 400 km circular orbit. Compute the orbital period $T$ and the number of orbits per day.

Recall Solution 2.3

Why this formula? For a circular orbit gravity supplies exactly the centripetal force: $μ m / r^{2} = m v^{2} / r$ . Solving gives $v = μ / r$ , and $T = 2 π r / v = 2 π r^{3} / μ$ (Kepler's third law). Here $r$ = orbit radius from Earth's centre, $μ$ = gravitational parameter. $r = 6371 + 400 = 6771 km$ $T = 2 π \frac{677 1 ^{3}}{398600} = 5545 s \approx 92.4 min .$ Orbits/day $= 1440/92.4 = 15.6$ .

Level 3 — Analysis (combine steps, interpret the result)

Exercise 3.1 (L3). A satellite spends 22 minutes per orbit in the SAA (dose rate $0.04 rad/s$ ) and the rest outside (dose rate $0.001 rad/s$ ). Using $T = 92.4 min$ from Ex. 2.3, find the daily TID.

Recall Solution 3.1

Orbits/day $= 1440/92.4 = 15.58$ . Time in SAA/day $= 15.58 \times 22 \times 60 = 20, 570 s$ . Time outside/day $= 86400 - 20570 = 65, 830 s$ . $TID_{day} = 0.04 (20570) + 0.001 (65830) = 822.8 + 65.8 = 888.6 rad/day .$ Interpretation: the 22-minute SAA slice contributes $822.8$ rad while the other 70 minutes contribute only $65.8$ rad — the SAA delivers ~93 % of the daily dose in ~24 % of the time. That is exactly why satellites power down during SAA passes.

Exercise 3.2 (L3). Kapton, AO flux $Φ = 8 \times 1 0^{13} atoms/cm^{2} / s$ , erosion yield $Y = 3 \times 1 0^{- 24} cm^{3} / atom$ , mission 3 years. Find erosion depth in µm. Would a 50 µm film survive?

Recall Solution 3.2

Why fluence? Damage is cumulative — like TID for radiation, each atom removes its own tiny volume, so we integrate flux over time to get fluence $F$ . Here $Φ$ = flux (atoms/cm²/s), $t$ = time, $Y$ = volume removed per atom. $F = Φ t = 8 \times 1 0^{13} \times (3 \times 365 \times 86400) = 7.57 \times 1 0^{21} atoms/cm^{2} .$ Volume removed per cm² of surface (which numerically equals the depth, since cm³/cm² = cm): $d = F Y = 7.57 \times 1 0^{21} \times 3 \times 1 0^{- 24} = 0.0227 cm = 227 μ m .$ A 50 µm film is destroyed — it lasts only $50/227 \times 3 = 0.66$ years ≈ 8 months. Coatings ( $SiO_{2}$ , $Al_{2} O_{3}$ ) are mandatory. This connects to Spacecraft materials selection.

Level 4 — Synthesis (build a design answer from multiple pieces)

Exercise 4.1 (L4). MMOD Poisson risk. A spacecraft has exposed area $A = 80 m^{2}$ , mission $t = 8 years$ , and the flux of $\geq 1 mm$ particles is $F = 3 \times 1 0^{- 6} impacts/m^{2} / year$ . Find (a) expected number of impacts $λ$ , (b) probability of zero impacts, (c) probability of at least one impact.

Recall Solution 4.1

Why multiply the three factors? Flux is "impacts per m² per year." Multiply by the area exposed and by the years of exposure to get the expected count — a Poisson mean $λ$ (the average number of hits). $λ = F A t = 3 \times 1 0^{- 6} \times 80 \times 8 = 1.92 \times 1 0^{- 3} .$ Poisson: $P (k) = \frac{λ ^{k} e ^{- λ}}{k !}$ (probability of exactly $k$ impacts). $P (0) = e^{- λ} = e^{- 1.92 \times 1 0^{- 3}} = 0.99808.$ $P (\geq 1) = 1 - P (0) = 1.918 \times 1 0^{- 3} \approx 0.19%.$ Design note: for a 1 mm threat over 8 years the risk is small, but scaling to 0.1 mm particles (flux ~ $1 0^{4} \times$ higher) makes surface erosion near-certain — hence protective coatings and, for critical modules, debris-avoidance manoeuvres and shields.

Exercise 4.2 (L4). Whipple shield ballistic limit. The critical projectile diameter that just perforates the inner wall is often written in the empirical form $d_{c} = C (\frac{σ t _{w}}{ρ _{p}^{1/2} v ^{2/3}})^{6/5},$ where $σ$ = wall yield strength, $t_{w}$ = wall thickness, $ρ_{p}$ = projectile density, $v$ = impact speed, and $C$ is a calibration constant that carries all the units so that $d_{c}$ comes out in metres. To keep the arithmetic transparent, define the bracketed group $G = \frac{σ t _{w}}{ρ _{p}^{1/2} v ^{2/3}}$ . Using $σ = 300 MPa$ , $t_{w} = 3 mm$ , $ρ_{p} = 2700 kg/m^{3}$ , $v = 10 km/s$ (all SI), compute (a) the group $G$ , (b) $G^{6/5}$ , and (c) explain how the calibration constant turns that into a physical millimetre-scale diameter.

Recall Solution 4.2

Why this formula? A hypervelocity impact is far too violent to solve exactly (both projectile and wall melt and vaporize). So engineers write down the physical energy balance — the projectile's kinetic energy must go into cratering plus flinging debris — and note it depends only on wall strength $σ$ , wall thickness $t_{w}$ , projectile density $ρ_{p}$ and impact speed $v$ . Dimensional analysis then fixes the combination those must appear in, and thousands of gas-gun test shots calibrate the exponent ( $6/5$ ) and the constant $C$ . It applies to thin metallic walls at $v > 3 km/s$ — exactly the LEO debris regime.

Step (a) — the group $G$ (work in strict SI):

$σ = 3 \times 1 0^{8} Pa$ , $t_{w} = 3 \times 1 0^{- 3} m$ .
$ρ_{p}^{1/2} = 2700 = 51.96$ .
$v^{2/3} = (1 0^{4})^{2/3} = 464.16$ . $G = \frac{( 3 \times 1 0 ^{8} ) ( 3 \times 1 0 ^{- 3} )}{51.96 \times 464.16} = \frac{9 \times 1 0 ^{5}}{24118} = 37.32.$

Step (b) — raise to the $6/5$ power: $G^{6/5} = 37.3 2^{1.2} = 76.86.$

Step (c) — the calibration constant $C$ closes the units. The bracket $G^{6/5}$ is not itself a length — it is a bare number ( $76.86$ ) whose SI dimensions are absorbed into $C$ . The empirical constant for aluminium is $C \approx 1.0 \times 1 0^{- 4} m$ , so $d_{c} = C G^{6/5} = (1.0 \times 1 0^{- 4}) (76.86) = 7.69 \times 1 0^{- 3} m \approx 7.7 mm .$ Reading it: particles below ~7.7 mm are stopped by this wall; larger ones perforate. Notice $d_{c}$ grows with $t_{w}$ and $σ$ (thicker, stronger wall stops bigger rocks) but shrinks with $v$ (faster debris is deadlier). This is exactly why the two-plate Whipple idea — a thin bumper to shatter the projectile, then a spaced catcher wall — beats one thick plate.

Level 5 — Mastery (full mission trade with judgement)

Exercise 5.1 (L5). A 5-year LEO mission at 400 km, 51.6° inclination faces: (i) SAA-dominated TID of $\approx 900 rad/day$ (from Ex. 3.1), and (ii) commercial electronics rated to $8 krad$ . Decide whether commercial parts survive unshielded, and if not, propose two mitigation paths. Support each with a number.

Recall Solution 5.1

Step 1 — total dose over 5 years: $TID_{5 yr} = 900 rad/day \times 365 \times 5 = 1.6425 \times 1 0^{6} rad = 1642.5 krad .$ Step 2 — margin check: required tolerance $/$ part rating $= 1642.5/8 \approx 205 \times$ . Commercial parts fail by more than two orders of magnitude — unshielded operation is impossible.

Step 3 — mitigations, each quantified:

Rad-hard parts rated to $\sim 300 krad$ — still short of 1642 krad, so combine with modest shielding. Shielding roughly halves dose per few mm Al; going from 2.5 mm to ~8–10 mm Al can cut the behind-shield dose by a large factor (the parent note's worst-case number is unshielded-dominated). Trade: mass ↑.
Operational duty-cycling — power down sensitive electronics during the ~22 min SAA pass. From Ex. 3.1 the SAA supplies ~93 % of daily dose; blanking it removes most of the upset risk (single-event effects) even though cumulative TID to powered-off transistors is only partly reduced. See Single-event effects.

Judgement: the winning architecture is both — rad-hard/rad-tolerant parts plus SAA-aware operations plus shielding sized against a realistic (not worst-case) behind-shield dose. Any single measure alone is insufficient given the $205 \times$ shortfall.

Exercise 5.2 (L5). Two candidate orbits: (A) 400 km, 51.6° inclination (crosses SAA and dense atmosphere); (B) 800 km, 51.6° (higher, thinner atmosphere, but deeper into the inner belt). For each of atomic-oxygen erosion, radiation dose (TID), and orbital lifetime, state which orbit is worse and why, then give the overall design judgement.

Recall Solution 5.2

Compare the three hazards one at a time, giving the direction of each effect:

Atomic oxygen (AO): ambient neutral density falls steeply with altitude, so the AO flux $Φ$ at 800 km is far lower than at 400 km. Orbit A (400 km) is worse for AO. Raising to 800 km therefore reduces material erosion ↓ — good for Kapton and other organics (Spacecraft materials selection).
Radiation dose (TID): the inner Van Allen belt peaks around 1000–6000 km. At 800 km you sit deeper inside that trapped-proton region than at 400 km, so the trapped-proton dose is higher. Orbit B (800 km) is worse for TID ↑ — bad for electronics and Single-event effects.
Orbital lifetime: atmospheric drag is proportional to neutral density, which is far smaller at 800 km. So the natural orbital-decay time grows from months–years at 400 km to decades or more at 800 km. Orbit B lives far longer ↑. This cuts both ways: a longer mission is good, but a dead satellite at 800 km lingers as debris for decades, feeding Kessler Syndrome and increasing collision risk for everyone (see Orbital lifetime in LEO).

Overall judgement: there is no free choice — every altitude knob moves several hazards at once. Going higher (B) helps AO erosion and mission lifetime but hurts radiation dose and worsens the post-mission debris legacy. Going lower (A) helps dose and debris disposal (fast natural deorbit) but punishes materials with AO and forces more frequent reboosts against drag. The correct pick depends on the dominant failure mode: if material erosion or debris stewardship drives the design, favour the lower orbit with a controlled deorbit plan; if AO-sensitive optics or a long baseline mission drive it, accept the higher orbit but budget for rad-hard parts and end-of-life disposal. Mastery is holding all the coupled effects on the table at once and choosing the least-bad compromise, not eliminating any single term.