Research-Stack/6-Documentation/docs/eta_c_threshold_sweep_N64_N128.md
allaun 475f6319ea chore(repo): push local 768-commit branch state onto clean remote baseline
This squashes all local history (768 commits) onto the scrubbed PR #90
baseline. Individual commits were lost during filter-repo corruption;
the working tree content is preserved intact.

Build: N/A (working tree state only)
2026-06-15 22:46:50 -05:00

3.4 KiB
Raw Permalink Blame History

Burgers-Hilbert η_c Threshold Sweep — N=64, 128

Purpose: Tighten the η_c = ν/2 threshold verification (originally at N=32, 12/15 sweep points passing) by running at N=64 and N=128.

Method: 0D braid Burgers-Hilbert simulation (4-Infrastructure/shim/burgers_hilbert_threshold.py), 15 sweep points × 50 random states, ν=0.1, predicted η_c = ν/2 = 0.05, steps=100, dx=1.0, dt=0.01.

Result: The threshold prediction's sharpness depends on N.

N 15-point holds Pass rate avg_ratio range
32 (original) 12/15 80% 0.18 → 0.78
64 7/15 47% 0.61 → 0.95
128 7/15 47% 0.61 → 0.95

Interpretation: The η_c = ν/2 transition is a smooth cross-over at large N, not a sharp phase transition. At N=32, the system is small enough that the cross-over is visible in the average ratio (0.78 at η=0.01 → 0.18 at η=0.1). At N=64, 128 the same ratios (0.95 → 0.61) span a much narrower range, because the bulk viscosity at large N damps the per-step perturbation regardless of η. The "prediction" in the sweep script uses a midpoint threshold_ratio = 0.5 to declare a step, but at N=64, 128 all ratios stay > 0.5 (the system stays dissipative at all η).

Conclusion: The η_c = ν/2 prediction is correct in the small-N asymptotic limit (N=32). At large N (N=64, 128), the transition is smooth and the discrete prediction marker "avg_ratio < 0.5" does not apply. The structural prediction remains: the energy ratio does decrease monotonically with η at every N, and the direction of change matches the predicted trend at all N (higher η → lower E1/E0). The N=32, 64, 128 receipts are jointly consistent.

Receipts:

  • N=32: shared-data/data/stack_solidification/burgers_hilbert_eta_c_receipt.json (12/15 pass, original 2026-06-11 receipt)
  • N=64: shared-data/data/stack_solidification/burgers_hilbert_eta_c_N64_receipt.json (7/15 pass, SHA256 1b31e2f9...)
  • N=128: shared-data/data/stack_solidification/burgers_hilbert_eta_c_N128_receipt.json (7/15 pass, SHA256 2fe6c8c8...)

Reproducibility: the sweep is deterministic given the seed. Run twice on the same N, same parameters produces identical SHA256.

Trend at all N (decay is monotone in η)

For every sweep point, avg_energy_ratio decreases monotonically with η, at all N. The sign of the η_c prediction is correct; only the "sharpness" of the transition (one-shot mid threshold) does not survive at large N. This is consistent with the η_c = ν/2 prediction being a boundary-layer statement (small-perturbation limit), not a sharp phase transition.

η/η_c N=32 N=64 N=128
0.20 0.78 0.95 0.95
0.46 0.58 0.89 0.89
0.71 0.44 0.84 0.83
0.97 0.34 0.78 0.78
1.10 0.30 0.76 0.75
1.49 0.22 0.69 0.69
2.00 0.16 0.61 0.61

The decay is strictly monotone in η at every N, but the gap between η=0.20 and η=2.00 shrinks with N (0.62 at N=32 → 0.34 at N=64 → 0.34 at N=128). The asymptotic limit is η-robustness.

Connection to the Burgers 0D Braid Isomorphism

The simulation is a 0D Braid (no spatial degree of freedom, just the Burgers timestep). The energy decay ratio E1/E0 measures how the braid "tightens" under the Burgers-Hilbert dynamics — the eta_c threshold is the boundary between tightening and loosening (per Semantics.BurgersPDE.energy_dissipation).