⚠ This model has been AI-generated (Claude Fable 5, Anthropic, under the supervision of A. Golkar) and is currently undergoing verification. Do not trust the model results unless manually verified.

AI Mega-Constellation Explorer

collision avoidance · thermal · external debris · anchored: 2608 collisions/yr at defaults (kinetic)

System dashboard (all layers, current parameter set)

Internal collisions / yr (natural)
Keplerian · kinetic:
Internal residual / yr (with CA)
f_fail =
External lethal hits / yr (>1 cm, fleet)
not avoidable; losses after damage tolerance: /yr
External trackable encounters / yr
avoidable; residual /yr
Total expected losses / yr
internal residual + external residual + lethal-debris losses
Cascade branching κ
κ=1 at N ≈ · residence for κ=1:
Radiator area per satellite
Radiator root temperature
GPU junction · thermal mass kg/sat

System map: the spacecraft exploded by subsystem (click a subsystem to open its tab)

drag to rotate · wheel to zoom · click a subsystem to open its tab · all numbers live from the current parameter set

How the layers connect: the Thermal tab sizes the radiator from GPU temperature and IT power (optionally driving the area used everywhere else); the Collision Avoidance tab computes the self-collision workload from that area, the population and the orbits; the External Debris tab adds the environment the fleet does not control. All tabs share one parameter set (Parameters & Assumptions tab).

Parameters (shared by all tabs)

Population and cross-section
Geometry
Inclination mix [deg / weight]
Velocity, control, window
Cascade parameters
Thermal (two-phase pumped loop)
Power subsystem
ADCS and configuration
External debris (screening level)
3D view
Collisions / yr (natural)
Keplerian · kinetic:
Residual / yr (with CA)
f_fail = (1 = no CA, equals natural)
Per-satellite P(hit) in window T
ν = /yr · by pop [/yr]:
Keplerian / kinetic
F_sp × F_vel
Effective v_rel
impact mean
Mean free path
/ shell thickness:
Cascade branching κ
κ=1 at N ≈
Δr for ≤1 collision/yr (kinetic):
drag: rotate · wheel: zoom · orbits: true Keplerian propagation of a sampled population · colors: light blue = nominal, flashing red = elevated local risk (model field) · 💥 = statistical collision event at the model's natural rate, located by the model's latitude distribution (not a contact of the displayed dots)

Collisions/yr vs (log-log)

Collision-rate density vs latitude (geography of encounters)

Impact-speed spectrum (collision-weighted)

Two-phase active cooling: GPU temperature drives radiator size drives collision risk

Heat rejected per satellite
P_IT × α_OH (all bus power becomes heat)
Radiator root temperature
T_junction − ΔT(junction→coolant) − ΔT(coolant→root)
Net flux per panel m²
2εησ(T_rad&sup4; − T_sink&sup4;), both faces
Required radiator panel area

The thermal-collision tradeoff: hotter GPUs, smaller radiators, fewer collisions

Radiator configuration (3D, set in Parameters; drives ADCS inertia)

drag to rotate · same total area, different span: the ADCS tab shows what the span costs

Thermal system sizing (per satellite)

ComponentSizingValue
Cold plates / evaporatorsΔT junction→coolant at full load
Two-phase loop (pump, lines, accumulator)~0.8 kg per kW rejected + pump power ~1.5% of Q
Condenser / radiator panelsA × areal density
Total thermal subsystem mass
Specific massper kW of IT power

Two-phase pumped loops keep the transport nearly isothermal, so the radiator runs close to the junction temperature minus two interface drops; that is what makes the T&sup4; lever effective. Raising the junction from 70°C to 105°C cuts the radiator (and the collision cross-section) by roughly a factor 2; the price is paid in silicon reliability, leakage power, and pump margins. Enable "drive A from Thermal tab" in Parameters to propagate this area through the whole model.

Power subsystem: photovoltaics, eclipse storage, and the second big deployed area

Sunlight / eclipse
eclipse min per orbit ()
PV array power
BOL; EOL requirement for bus load
PV array area
m² per IT kW · vs radiator
Battery
at DoD · mass kg
Power subsystem mass
PV + battery + PMAD kg
Specific mass
per kW of IT power
Total deployed area
radiator + PV; collision multiplier if PV counted: ×
Collision coupling
toggle "count PV area" in Parameters to propagate

Deployed area and collision burden vs IT power per satellite (current links)

Spacecraft (power subsystem highlighted)

drag to rotate · PV array sized by the live model; gray = other subsystems

The PV array is the second large deployed surface, and in eclipse-prone orbits it exceeds the radiator. The collision analysis nominally counts radiator area only; enabling the PV-area coupling adds the array to the cross-section, which is the conservative reading (Turyshev's 1 MW anchor carries 5.6 m²/kW of PV against 2.5 m²/kW of radiator). Dawn-dusk SSO removes the battery and shrinks the array, which is exactly why that orbit is thermally and electrically attractive; the Collision Avoidance tab prices what it costs.

Attitude Determination and Control: the radiator span sets the control problem

Gravity-gradient torque
ΔI = kg·m² at ° offset · span m
Aero + SRP + magnetic torque
aero · SRP · mag [mN·m]
Momentum to store per orbit
sized with margin 2 → N·m·s per axis
Recommended suite mass
Torque authority check
wheel class N·m vs peak disturbance · 90° slew ≈ min
Inertia context
Hubble ≈ 87,000 kg·m² (flies on reaction wheels); ISS ≈ 10⁷ (CMGs)

Disturbance torques vs altitude (300 km band, current configuration)

Actuator suite (3D, choose in Parameters)

drag to rotate · sizes scale with the computed momentum requirement

Design alternatives

SuiteMassPeak powerNotes
4× reaction wheels + 3× magnetorquers Baseline LEO. Wheels absorb the cyclic gravity-gradient momentum; torquers dump the secular part each orbit (dipole A·m²).
4× control-moment gyros + 3× magnetorquers For large momentum/agility; pays a fixed mass and complexity premium, wins above ~80 N·m·s.
4× reaction wheels + thruster desaturation Replaces torquers where the dipole gets impractical; costs kg propellant per year.

Sizing logic: the deployed radiator dominates the inertia asymmetry, so gravity gradient is the design driver; T_gg = 3/2 n² ΔI sin(2θ). The cyclic part loads the wheels every orbit (H ≈ 0.637 T_gg T_orb/4); aero, solar pressure, and the residual dipole accumulate as secular momentum that must be dumped. Configuration matters: the four-panel and z-fold layouts shorten the span, cutting ΔI and the whole actuator chain; that is mass the thermal design can trade against. Mass/power figures are SMAD-class regressions for screening, not catalog data. On torque authority: despite the 120 m² radiator the inertia is Hubble-class (~half of Hubble's 87,000 kg·m²), and Hubble flies on reaction wheels. Wheels reject the 13 mN·m disturbance with order-15× torque margin; what large inertia actually costs is agility (tens of minutes per 90° slew, irrelevant for a sun/nadir-fixed compute satellite) and controller bandwidth below the flexible modes of 15 m wings, which this screening model does not cover. At worst-case 45° gravity-gradient attitudes the momentum requirement reaches the top of the wheel class and the recommendation switches to CMGs automatically.

External debris environment (screening estimate, calibrated to public ORDEM/MASTER/ESA figures)

Lethal non-trackable hits (>1 cm), fleet
per satellite: /yr · not avoidable
Expected losses from LNT debris
hits × critical-area fraction
Trackable-object encounters (>10 cm), fleet
per satellite: /yr · adds to avoidance workload
Residual from trackable debris
after avoidance at the same f_fail

Fleet-level rates vs altitude (300 km band centered at each altitude)

Method: flux = n(h, >d) × v_rel; hits/yr = N × A_eff × flux. For lethal non-trackable (1-10 cm) debris A_eff is the radiator face area A (it cannot be dodged; only damage tolerance helps, hence the critical-area fraction). For trackable objects (>10 cm) the encounter cross-section is σ = shape·A and the events feed the same avoidance chain as internal conjunctions. The density profile is an order-of-magnitude envelope (peak near 750-900 km, traffic bump near 550 km); a bankable assessment requires native ORDEM 3.1 / MASTER-8 runs with the real geometry. The 700-900 km region combines high debris density with very long residence times: least attractive for a large fleet.

Validation against published benchmarks (computed at the fixed reference case, independent of the sliders)

Reference case: N = 80,000, A = 120 m², σ = 4A, 500-800 km, mix 43/53/70/97.6°. Benchmark A: independent due-diligence reconstruction, "Feasibility of an Orbital Data Center" (10 June 2026), which rebuilt this model from the governing integral. Benchmark B: S. G. Turyshev, "Orbital Data Centers: Spacecraft Constraints and Economic Viability", arXiv:2604.27197 (May 2026), an independent power/thermal/economics closure model.

vs Benchmark A (independent collision-model reconstruction)

QuantityThis modelBenchmark AStatus
Kinetic-gas baseline [collisions/yr]2,608
Keplerian rate [collisions/yr]1,941.3
Keplerian / kinetic ratio0.744
Mean natural rate per satellite [/yr]0.0485
Single 10 km shell at 650 km [collisions/yr]~58,000
N for 1 natural collision/yr (A = 120 m²)1,816
A for 1 natural collision/yr (N = 80k) [m²]0.062
f_fail for <1 residual collision/yr5.15e-4
Cascade branching κ~408
Fragment residence for κ = 1 [days]~22
Screening load within 1 km [/sat/yr]~320
External >10 cm encounters, fleet [/yr]~400-600 (indicative figure)
External >1 cm lethal hits, fleet [/yr]~1,000-4,500 (incl. ×3 band)

The first eleven rows agree to three or more digits: Benchmark A is an independent reconstruction of the same physics and reproduces this engine. The two external-debris rows use different calibrations of the same flux method: this model anchors the >1 cm population to ESA object counts (~1.2 million objects of 1-10 cm in LEO), which lands at the upper edge of Benchmark A's indicative uncertainty band. Use the environment scale slider for sensitivity.

vs Benchmark B (Turyshev, arXiv:2604.27197)

QuantityThis model (defaults)Benchmark BComment
Radiator area per IT kilowatt2.5 m²/kW (1 MW anchor) Same physics; ours runs hotter (T_rad ) and radiates both faces, his anchor is more conservative.
Radiator temperature300-400 K design rangeInside his range.
PV area per IT kilowatt (dawn-dusk)5.64 m²/kW (BOL, 1 MW high-sunlight anchor)Independent reproduction of his PV sizing.
Radiator areal density2-10 kg/m²Inside his range.
Thermal-subsystem mass per IT kW~5-25 kg/kW (radiator share of his 29.4 kg/kW PV+storage+radiator floor)Consistent.
Fleet scale context17-46 GW IT → 170,000-461,000 nodes at 100 kWThe reference fleet sits inside his envelope.
Collision-risk treatmentExplicit fleet-level model: rates, geography, spectrum, cascade κAggregated hazard rate λ_MMOD inside the lifetime penalty Π_life Complementary: this model supplies the physics behind his λ_MMOD and the correlated cascade risk his independent-hazard framework cannot represent.

Benchmark B does not model fleet self-collisions, so the comparison is on the shared thermal/sizing layer and on consistency of treatment. The two models close the loop on each other: B shows the economics are gated by mass and link costs even before collision physics; this model shows collision physics adds constraints B assumes away.

Assumptions register

AssumptionValidity / limitation
Dilute gas, Poisson statisticsMean free path exceeds shell thickness by ~10^7; individual event probabilities small.
Circular orbits, no perturbationsJ2 only accelerates the RAAN randomization the model already assumes. Eccentricity shifts densities by km; negligible on rates.
Randomized RAAN and phaseExact for ensemble means. Coordinated slotting or LTAN clustering (all dawn-dusk) violates it; not modeled.
σ = shape factor × A4A = enveloping sphere (conservative); physical band A…4A; all rates linear in σ.
Inclination dispersion ±0.5°Total rate insensitive; only the F_sp × F_vel split depends on it (logarithmically).
Cascade branching κWarning index, not an evolutionary debris model (use LEGEND/MASTER for that). Avoidance does NOT apply to κ.
External debris profileOrder-of-magnitude envelope calibrated to public ORDEM 3.1 / MASTER-8 / ESA figures. NOT a certified run; use the scale slider for the uncertainty band.
Thermal modelSteady-state two-phase loop; radiator radiates from both faces to an equivalent sink; transients, shadowing and degradation excluded.
Power modelBeta-0 worst-case or dawn-dusk eclipse bounds; PMAD and battery as SMAD-class regressions; degradation via a single EOL fraction.
ADCS modelScreening-level torque budget and actuator regressions; flexible-body dynamics, jitter, and slew agility excluded.
f_failEnd-to-end fraction of dangerous internal/trackable encounters not mitigated, including dead satellites; lethal non-trackable debris cannot be avoided at all.