Dane Sabo c5133401e0 Session work scratch: scram X_exit refactor, hot-standby SOS, fat scram tubes, model cheatsheet, journal entry
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
2026-05-02 23:02:50 -04:00

pwr-hybrid-3-demo

Preliminary example for the HAHACS thesis — a verified hybrid controller for a small modular PWR startup. Composes three layers into one demonstrable pipeline:

  • Discrete layer (fret-pipeline/): FRET natural-language requirements → LTL → synthesized AIGER controller → state-machine diagram.
  • Continuous layer (code/): 10-state point kinetic equation + thermal-hydraulics PWR model with bounded steam-generator heat removal as the disturbance input. Controllers, linearization, LQR, reach-tube propagator, Lyapunov barrier — all Julia.
  • Verification artifacts (reachability/): predicate concretizations (single source of truth in predicates.json) and the standalone reach analysis writeup (WALKTHROUGH.md).
  • Research context (thesis/): the HAHACS PhD proposal.
  • Lab journal (journal/): chronological invention log in LaTeX.
  • Predicate explorer app (app/): Pluto.jl notebook bridging FRET predicates and continuous-state halfspaces.

Layout

pwr-hybrid-3-demo/
  CLAUDE.md                AI-facing context and architecture map
  docs/
    architecture.md        How the layers compose
    figures/               Shared figures for thesis + talks
  fret-pipeline/           FRET → ltlsynt → AIGER → state machine
  code/                    Plant model, controllers, reach (all Julia)
  reachability/            predicates.json + WALKTHROUGH.md
  app/                     Pluto.jl predicate explorer
  journal/                 LaTeX lab notebook
  hardware/                Ovation HIL artifacts (TBD)
  claude_memory/           Short AI-context notes
  thesis/                  [submodule] PhD proposal
  presentations/
    2026DICE/              [submodule] DICE 2026 abstract

Quickstart

Clone with submodules:

git clone --recurse-submodules <url>
cd pwr-hybrid-3-demo

Run the controller synthesis pipeline:

cd fret-pipeline
python3 scripts/fret_to_synth.py pwr_hybrid_3.json specs/synthesis_config_v3.json
bash scripts/synthesize.sh specs/synthesis_config_v3.json circuits
python3 scripts/trace_aiger.py circuits/PWR_HYBRID_3_DRC.aag diagrams
dot -Tpng diagrams/PWR_HYBRID_3_DRC_states.dot -o diagrams/PWR_HYBRID_3_DRC_states.png

Run the plant model and reach analysis:

cd code
julia --project=. -e 'using Pkg; Pkg.instantiate()'    # first time only

julia --project=. scripts/main_mode_sweep.jl           # all 5 DRC modes
julia --project=. scripts/reach_operation.jl           # operation-mode linear reach
julia --project=. scripts/barrier_lyapunov.jl          # Lyapunov barrier
julia --project=. scripts/barrier_compare_OL_CL.jl     # OL vs CL barrier
julia --project=. scripts/reach_heatup_nonlinear.jl    # nonlinear heatup (10s cap)

Open the predicate explorer:

cd app
julia --project=. -e 'using Pluto; Pluto.run()'
# Browser opens; navigate to predicate_explorer.jl

Soundness note: the current reach tubes are over-approximations of the LINEAR model, not sound over-approximations of the nonlinear plant. See reachability/README.md and reachability/WALKTHROUGH.md.

Prerequisites

  • Julia 1.10+ (via juliaup).
  • Python 3.10+ (FRET pipeline only).
  • Spot for ltlsynt (brew install spot).
  • Graphviz for dot (brew install graphviz).
  • LaTeX (via latexmk) for the thesis + journal builds.

Further reading

  • CLAUDE.md — orientation for AI agents working in this repo
  • docs/architecture.md — how the layers compose
  • code/CLAUDE.md — code architecture, conventions, validity range
  • code/README.md — usage and dependencies
  • reachability/README.md — reach scope, soundness status
  • reachability/WALKTHROUGH.md — standalone analysis writeup
  • journal/README.md — journal format conventions
  • journal/journal.tex — the journal itself, dated entries
  • thesis/CLAUDE.md — the thesis project structure
  • fret-pipeline/README.md — FRET naming conventions and pipeline details
Description
No description provided
Readme 6.7 MiB
Languages
Julia 46.5%
TeX 34.5%
Python 17.5%
Shell 1.5%