Add __main__ demos to power.py and thermal.py

Now runnable: python -m orbital_compute.power (battery + aging demo)
and python -m orbital_compute.thermal (orbit thermal cycling demo).
7 more core modules still need __main__ blocks.

Co-Authored-By: Claude Opus 4.6 (1M context) <noreply@anthropic.com>

Author: noreply

orbital_compute/power.py

@@ -195,3 +195,35 @@ def predict_eol(self, cycles_per_day: float, avg_dod: float = 30.0,
 
         cycles_to_eol = remaining_fade / fade_per_cycle
         return cycles_to_eol / cycles_per_day
+
+
+if __name__ == "__main__":
+    print("=" * 60)
+    print("  POWER SUBSYSTEM DEMO")
+    print("=" * 60)
+
+    # Simulate 3 orbits
+    pm = PowerModel(PowerConfig(solar_panel_watts=2000, battery_capacity_wh=5000))
+    print(f"\n  Solar: {pm.config.solar_panel_watts}W, Battery: {pm.config.battery_capacity_wh}Wh")
+    print(f"\n  Simulating 3 LEO orbits (95 min each)...")
+    for orbit in range(3):
+        # Sunlit phase (60 min)
+        for _ in range(60):
+            s = pm.step(60, in_eclipse=False, compute_load_w=500)
+        sun_batt = s.battery_pct
+        # Eclipse phase (35 min)
+        for _ in range(35):
+            s = pm.step(60, in_eclipse=True, compute_load_w=500)
+        ecl_batt = s.battery_pct
+        print(f"  Orbit {orbit+1}: sunlit→{sun_batt:.0%} eclipse→{ecl_batt:.0%} "
+              f"solar={s.solar_output_w:.0f}W load={s.load_w:.0f}W")
+
+    # Battery aging
+    print(f"\n  Battery Aging (5 years, 15 cycles/day, 30% DoD):")
+    ba = BatteryAgingModel(initial_capacity_wh=5000)
+    for year in range(1, 6):
+        for _ in range(365 * 15):
+            ba.cycle(dod_pct=30, c_rate=0.5, temp_c=28)
+        ba.calendar_age(8760)
+        print(f"  Year {year}: SoH={ba.soh_pct:.1f}% capacity={ba.capacity_wh:.0f}Wh")
+    print(f"\n{'=' * 60}")

orbital_compute/thermal.py

@@ -104,3 +104,22 @@ def max_sustainable_heat_w(self, target_temp_c: float = 70.0) -> float:
         q_solar = c.solar_absorptivity * 1361.0 * c.radiator_area_m2 * 0.25
         q_earth = c.earth_ir_w_m2 * c.radiator_area_m2 * 0.5
         return max(0.0, q_rad - q_solar - q_earth)
+
+
+if __name__ == "__main__":
+    print("=" * 60)
+    print("  THERMAL SUBSYSTEM DEMO")
+    print("=" * 60)
+    tm = ThermalModel(ThermalConfig(radiator_area_m2=4.0), initial_temp_c=20.0)
+    print(f"\n  Radiator: {tm.config.radiator_area_m2} m², Start: 20°C")
+    print(f"  Max sustainable heat: {tm.max_sustainable_heat_w(70):.0f}W at 70°C")
+    print(f"\n  Simulating sunlit orbit with 500W GPU heat...")
+    for minute in range(95):
+        eclipse = minute >= 60
+        s = tm.step(60, 500 if not eclipse else 100, eclipse)
+        if minute % 10 == 0:
+            phase = "ECLIPSE" if eclipse else "SUNLIT"
+            print(f"  t={minute:3d}min {phase:7s} temp={s.temp_c:5.1f}°C "
+                  f"heat={s.heat_generated_w:.0f}W rad={s.heat_radiated_w:.0f}W "
+                  f"throttle={s.throttle_pct:.0%}")
+    print(f"\n{'=' * 60}")