c5133401e0
Multi-session work bundle on a draft branch. Splits into a clean
sequence of commits later; pushed here so it isn't lost on a reboot.
Reach work
- code/scripts/reach/reach_scram_pj.jl: shutdown_margin halfspace
X_exit (replaces "n <= 1e-4 AND T_f bound" framing); per-step
envelope extraction added.
- code/scripts/reach/reach_scram_pj_fat.jl: per-step envelope
extraction added; shutdown_margin discharge logic mirrored from the
tight scram script. 3 probes (10/30/60s) all discharge from the
fat union polytope.
- code/scripts/reach/reach_scram_full_fat.jl (NEW): full nonlinear
PKE scram reach with fat entry. Hits the stiffness wall at
~1.5 s plant time as expected; saves NaN-tolerant per-step
envelopes. Demonstrates concretely why PJ is the right tool for
the longer-horizon proof.
- code/scripts/reach/reach_heatup_pj.jl: T_REF_START_C constant
(entry-conditioned ramp) replaces T_STANDBY-init that was making
the FL controller command cooling at t=0. Per-step extraction
already in place.
- code/configs/heatup/tight.toml: bumped maxsteps; probe horizon
parameterized.
Hot-standby SOS barrier
- code/scripts/barrier/barrier_sos_2d_shutdown.jl (NEW): mirrors the
operation SOS machinery on the hot-standby thermal projection.
Includes the eps-slack pattern (so feasibility doesn't silently
collapse to B == 0).
- code/scripts/barrier/barrier_sos_2d.jl: refactored to use the same
helper.
- code/src/sos_barrier.jl (NEW): solve_sos_barrier_2d helper module
factoring out the SOS construction; eps-slack with eps_cap=1.0 to
avoid unbounded primal.
Library
- code/src/pke_states.jl (NEW): single source of truth for canonical
initial-condition vectors per DRC mode (op, shutdown, heatup) keyed
off plant + predicates.
- code/scripts/sim/{main_mode_sweep,validate_pj}.jl, code/CLAUDE.md:
migrated to pke_states.
Predicates + invariants
- reachability/predicates.json: new shutdown_margin predicate (1%
dk/k tech-spec floor, expressed as alpha_f*T_f + alpha_c*T_c
halfspace). Used as scram X_exit.
Plot script
- code/scripts/plot/plot_reach_tubes.jl: plot_tubes_scram_pj() with
variant=:fat|:tight knob; plot_tubes_scram_full() for full-PKE
3-panel (T_c, T_f, rho); plot_tubes_heatup_pj() reads results/
not reachability/.
Journal + memory
- journal/entries/2026-04-27-shutdown-sos-and-scram-X_exit.tex (NEW):
long-form entry on the SOS hot-standby barrier and the scram X_exit
refactor.
- journal/journal.tex: input chain updated.
- claude_memory/ — three new session notes:
* 2026-04-27-scram-X_exit-shutdown-margin.md
* 2026-04-28-DICE-2026-conference-intel.md (people, sessions,
strategic notes for the May 12 talk)
* 2026-04-28-path1-sos-pj-sketch.md (sketch of nonlinear-SOS via
polynomial multiply-through; saved for an overnight session)
Docs
- docs/model_cheatsheet.md (NEW): one-page reference of state vector,
dynamics, constants, modes, predicates, sanity numbers — the talk
prep cheatsheet Dane asked for.
- docs/figures/reach_*_tubes.png: regenerated with the new mat data.
- presentations/prelim-presentation/outline.md: revised arc per the
April-28 review pass (cuts: Lyapunov-fails standalone slide,
operation-tube standalone slide, SOS standalone; adds: scopes-of-
control framing, scram on the headline result slide).
- app/predicate_explorer.jl: minor.
Hacker-Split: end-of-session scratch bundle
39 lines
1.3 KiB
TOML
39 lines
1.3 KiB
TOML
# Tight X_entry heatup reach: T_c in [285, 291] (6 K vs 14 K baseline).
|
|
# Produces the clean "all 6 inv1_holds halfspaces discharged at T=300s"
|
|
# result from 2026-04-20 overnight session.
|
|
|
|
name = "tight"
|
|
description = "Narrow X_entry on T_c/T_f/T_cold; n in [1e-3, 2e-3]"
|
|
|
|
use_predicates_entry = false
|
|
|
|
[entry]
|
|
n_range = [1.0e-3, 2.0e-3]
|
|
T_f_range_C = [285.0, 291.0]
|
|
T_c_range_C = [285.0, 291.0]
|
|
T_cold_range_C = [278.0, 285.0]
|
|
|
|
[tmjets]
|
|
orderT = 4
|
|
orderQ = 2
|
|
abstol = 1e-9
|
|
maxsteps = 1000000
|
|
|
|
[probes]
|
|
# Single probe at the nominal heatup completion time.
|
|
# At T_REF_START_C = 285, RAMP_RATE = 28 C/hr, T_ref reaches T_c0 = 308.35
|
|
# at t = (308.35 - 285) / (28/3600) = 3001 s. Probing here checks whether
|
|
# the tube has entered X_exit (t_avg_in_range).
|
|
#
|
|
# Why not also probe at T_min = 7714 s as the formal obligation requires:
|
|
# the demo heatup controller has no clamp on T_ref. Ramping past T_c0
|
|
# would drive the tube past t_avg_high_trip before T_min. Discharging the
|
|
# full obligation needs a clamped controller (a smooth-min compatible with
|
|
# @taylorize) — flagged as the next thrust. Tonight we discharge the
|
|
# nominal-heatup-time entry into X_exit; the controller redesign follows.
|
|
horizons_seconds = [3000.0]
|
|
|
|
[output]
|
|
save_per_step = true
|
|
result_file = "reach_heatup_pj_tight.mat"
|