This 2D version replaces the flat "one fixed emissivity times ΔT⁴" bookkeeping with a real spectral calculation: at every frame it numerically integrates Planck's law over wavelength to find exactly what fraction of the surface's thermal glow actually falls inside the 8-13 µm atmospheric window — the left chart draws that integral directly, shading the escaping (orange) versus reabsorbed (brown) area under the curve.
M(λ,T) = 2πhc² / { λ⁵ [exp(hc/λkT) − 1] } [Planck spectral exitance]
F(T) = ∫₈¹³ M(λ,T) dλ / σT⁴ [window-band fraction, Simpson's rule]
P_solar = (1 − R_solar) · S [absorbed sunlight, S ≈ 900 W/m², day only]
P_window = ε_win · [ F(T_s)·σT_s⁴ − F(T_sky)·σT_sky⁴ ]
P_other = ε_other · [ (1−F(T_s))·σT_s⁴ − (1−F(T_amb))·σT_amb⁴ ]
P_conv = h · (T_amb − T_s)
dT_s/dt = ( P_solar + P_conv − P_window − P_other ) / C_eff
Cross-check against the 3D twin: re-deriving this ODE from the stated steady-state balance and testing it numerically (Planck-integral verified to 1.2×10⁻⁵ against standard blackbody-fraction tables; see verification script) turned up a sign error in the flat-emissivity 3D model — it subtracts the absorbed solar term from the rate of change instead of adding it, so in that version a more reflective coating paradoxically ends up warmer, and reflectance sweeps run backwards. This 2D model uses the corrected sign shown above, so higher Rsolar and higher εwin both correctly push the steady state colder, exactly as the passive-cooling literature (and Test 7 in the verification script) predicts.
- F(T) — because it comes from the real Planck curve rather than a constant, it shifts with temperature (Wien's law): a cooler surface's peak sits closer to the 8-13 µm band, so F rises slightly as Ts falls — a genuine physical detail the flat-emissivity model cannot show.
- Rsolar — fraction of sunlight the coating bounces away; every percent left over is absorbed heat (Psolar).
- εwin — how strongly the coating emits specifically inside the window; higher means more of F(Ts)·σTs⁴ actually escapes to space.
- Tsky ≈ Tamb − 30 K (the window's effective view of space); the non-window band exchanges instead with the atmosphere near Tamb (εother fixed at 0.15).