Van Allen Beltinator 3D Fokker-Planck Radiation Belt Explorer
Rendering T+0.0h
Physics Parameters
L* (Drift Shell) 4.5
Energy (MeV) 1.0
Equatorial Pitch Angle α₀ (°) 45
DL*L* (Radial Diffusion) 1.0e-6 ⟲
Dpp (Momentum Diffusion) 1.0e-7
Dα₀α₀ (Pitch-Angle Diffusion) 1.0e-5
Dα₀p (Mixed Diffusion) 0 (off)
Loss Timescale τ (days) 10.0
Geomagnetic Conditions
Kp Baseline 3
Dst Index (nT) -30
Solar Wind Vsw (km/s) 400
Quiet Conditions
Kp(t) — 3-Hour Cadence Profile
0h12h24h
Visualization
Sim Speed (×) 60
Flux Intensity (cm⁻²s⁻¹sr⁻¹MeV⁻¹)
10²10⁴10⁶10⁸
Presets
Satellites (Click in 3D to zoom)
Van Allen Probes (RBSP A/B)
THEMIS A / D / E
GOES (GEO)
LEO Satellites
Lomonosov

Computed Values

T(α₀) Bounce Factor 1.0802
Phase Space Density f 2.34e-10
Differential Flux j 4.56e+04
DL*L* (Kp-driven) 1.0e-6
Kp(t) Current 3.0
Magnetopause Lmp 10.2 RE
Last Closed Drift Shell 7.8 RE
Bounce Loss Cone αBLC 3.2°
Drift Loss Cone αDLC 5.1°
Particle Status Trapped

Governing Equations

Eq.1 — Generalized Diffusion Equation
∂f/∂t = ∂/∂Jᵢ (Dᵢⱼ ∂f/∂Jⱼ)
i,j = 1,2,3

Phase space density f evolves via diffusion in the three adiabatic invariants (J₁, J₂, J₃), corresponding to gyration, bounce, and azimuthal drift motion.

Eq.2 — Fokker-Planck (no mixed terms)
∂f/∂t = L*² ∂/∂L*|μ,J (DL*L* L*⁻² ∂f/∂L*|μ,J)
+ 1/p² ∂/∂p|y,L (p² ⟨Dpp(y,p)⟩ba ∂f/∂p|y,L)
+ 1/(T(y)·y) ∂/∂y|p,L (T(y)·y·⟨Dyy(y,p)⟩ba ∂f/∂y|p,L)
− f/τ

Separates radial diffusion (L*), momentum diffusion (p), and pitch-angle diffusion (y = sin α₀) with bounce-averaged coefficients. DL*L* is now time-dependent via Kp(t).

Eq.3 — Fokker-Planck with Mixed Terms
∂f/∂t = L*² ∂/∂L*|μ,J (DL*L* L*⁻² ∂f/∂L*|μ,J)
+ 1/p² ∂/∂p|α₀,L p²(⟨Dpp⟩ ∂f/∂p + ⟨Dpα₀⟩ ∂f/∂α₀)
+ 1/(T(α₀)sin2α₀) ∂/∂α₀|p,L T(α₀)sin2α₀
  ·(⟨Dα₀α₀⟩ ∂f/∂α₀ + ⟨Dα₀p⟩ ∂f/∂p) − f/τ

Dα₀p is the mixed energy–pitch angle diffusion coefficient. Cross-diffusion terms couple momentum and pitch-angle evolution.

Bounce Period Factor T(α₀)
T(α₀) = 1.3802 − 0.3198·(sin(α₀) + sin²(α₀))

Approximation for the normalized bounce period (Schulz & Lanzerotti 1974). Particles near the loss cone have T ≈ 1.38; those near 90° have T ≈ 0.74.

DL*L*(Kp) — Brautigam & Albert 2000
DL*L* = D₀ · 10^(0.506·Kp − 9.325) · L*¹⁰

Radial diffusion coefficient parameterized by Kp index and L-shell. Higher Kp dramatically increases radial transport, driving inward diffusion and energization during storms. Kp is updated every 3 hours in this simulation.

Radiation Belt Physics

Relativistic electrons (>500 keV) undergo three periodic motions: gyration (μ), bounce (J), and azimuthal drift (Φ = L*). ULF waves drive radial diffusion, chorus waves produce local acceleration, and magnetopause shadowing causes loss on the dayside.

What this means in layman's terms: High-energy particles get trapped by Earth's magnetic field and form two donut-shaped belts around the planet. These particles constantly spiral around magnetic field lines, bounce back and forth between the north and south poles, and slowly drift around the Earth. During magnetic storms, powerful waves can push particles deeper into the belts (making them more energetic and dangerous to satellites), while other processes can knock particles out of the belts entirely. It's a constant tug-of-war between forces that fill the belts and forces that drain them.

Magnetic Field Notes

The local B-field determines energy from μ. Because electron spectra are steep, small ΔB errors produce large flux errors — especially at high L. At LEO, populations include Trapped, BLC, and DLC electrons requiring accurate magnetic mapping to classify.

What this means in layman's terms: Earth's magnetic field acts like an invisible container that holds these dangerous particles in place. The strength of the magnetic field at any location determines how energetic the particles are. Even small errors in measuring the magnetic field can lead to very wrong estimates of how many particles are present — this matters enormously for predicting radiation damage to satellites and astronauts. At low altitudes where spacecraft like the ISS orbit, some particles are stably trapped, while others are on their way to hitting the atmosphere.

L* 4.5 RE
E 1.0 MeV
α₀ 45°
Kp(t) 3.0
DLL 1.0e-6
Model Eq.2
Sim T+0.0h