SHBT Cold Fusion Reactor: Precision Simulation & Engineering Workbench

shbt-cf is the multi-domain physics simulation engine, hardware-in-the-loop (HIL) telemetry architecture, and numerical workbench for the Static Holographic Boundary Theory (SHBT) reactor specification. Primary physics derivations, structural bounds, and electrodynamic models are detailed in cf.pdf.

The repository is structured as a Rust Cargo workspace containing 9 crates paired with a native Python orchestration wrapper:

Cargo.toml                   # Workspace root (9 crates with strict DAG topology)
.cargo/config.toml           # Target CPU optimization, opt-level 3, strict IEEE 754 rules
crates/
  shbt-core-math             # Q64.64 fixed-point, MPFR, Lie algebra SO(3)/SE(3), Yoshida-6, DAZ/FTZ
  shbt-dielectric-floquet    # Floquet-Adler-Wiser dielectric matrix inversions & elastodynamics
  shbt-rcwa-optics           # 2D/3D vector RCWA, S-matrix recursion
  shbt-rcwa                  # High-level RCWA solver execution & Floquet LU verification
  shbt-fea-structural        # Belleville disc springs, Stoney stress, fatigue, viscoplastic kinetics
  shbt-contact               # EHD contact mechanics & projected damped active set (PDAS) solvers
  shbt-metrology-gum         # GUM Supplement 1 Monte Carlo uncertainty engine & dual numbers
  shbt-fabrication-hil       # Shared-memory HIL ring, ROM residual monitoring, CFT kinetics, transport
  shbt-py-bindings           # PyO3 FFI layer exposing Rust solvers to Python
python/shbt_cf               # Python orchestration suite, control, metrology, and workbench tools
scripts/
  check_dep_graph.py         # CI DAG topology validator
  verify_zerocopy.py         # Pointer-identity verification for PyO3 shared memory
sim_outputs/
  simulation_verification.json # Dynamic simulation outputs and closed-loop verification metrics
tests/
  test_full_multiphysics_pipeline.py  # Full multiphysics end-to-end simulation runner
  test_physics.py                     # Multi-domain physics validation suite
  test_py_bindings.py                 # FFI boundary integration tests
  test_update5_modules.py             # Module integration regression tests
  validate_benchmarks.py              # Automated performance benchmark validator
  verification_tests.rs               # Multi-crate integration verification suite

shbt-cf functions as the canonical cold fusion and thermal-hydraulic authority across the 10-repository Static Holographic Boundary Theory (SHBT) digital twin ecosystem. It defines the baseline energy generation models, non-equilibrium dynamic screening dynamics, dual-stage thermoelectric generation (TEG) recovery routines, and closed-loop two-phase thermal rejection circuits governing both stationary fusion power grids and deep-space instrumentation platforms.

                              +------------------------------------------+
                              |             [shbt-precision]             |
                              |      Computational Math & Cosmology      |
                              |     (512-bit MPFR / WZW Characters)      |
                              +--------------------+---------------------+
                                                   |
                 +---------------------------------+---------------------------------+
                 v                                 v                                 v
   +---------------------------+     +---------------------------+     +---------------------------+
   |       [shbt-power]        |     |         [shbt-cf]         |     |         [shbt-qc]         |
   |  Commercial Fusion Grid   |     |  1,800-Module LANR Array  |     | Bare-Metal Microkernel &  |
   |   (8,750 MW p-11B Twin)   |     |    & Thermal-Hydraulics   |     |   Photonic Quantum Bus    |
   +-------------+-------------+     +-------------+-------------+     +-------------+-------------+
                 |                                 |                                 |
                 +---------------------------------+---------------------------------+
                                                   |
                 +---------------------------------+---------------------------------+
                 v                                 v                                 v
   +---------------------------+     +---------------------------+     +---------------------------+
   |       [shbt-ghost]        |     |        [shbt-sglt]        |     |        [shbt-exotic]      |
   | Ghost Seed Propulsion &   |     | Synthetic Gravitational   |     | Multi-Protocol Spacetime  |
   |  Metric Interlocks (CCZ4) |     |  Lensing Telescope Stack  |     |  Co-Simulation & QI Audit |
   +-------------+-------------+     +-------------+-------------+     +-------------+-------------+
                 |                                 |                                 |
                 +---------------------------------+---------------------------------+
                                                   |
                           +-----------------------+-----------------------+
                           v                                               v
           +-------------------------------+               +-------------------------------+
           |          [shbt-recon]         |               |          [shbt-warp]          |
           |  Macro State Translocation &  |<=============>|   Holographic Warp Drive &    |
           |  Bipartite Fuel Foundry Twin  |  Dual-Gateway |    3+1D ADM Metric Engine     |
           |  (V_unified^macro / MAR-312C) |    Corridor   |  (Alcubierre / SHBT-MAR-312C) |
           +---------------+---------------+               +---------------+---------------+
                           ^                                               ^
                           |       +-------------------------------+       |
                           |       |          [shbt-swarm]         |       |
                           +=======|  Autonomous Von Neumann Swarm |=======+
                                   |  Universal Matter Foundry &   |
                                   |   5th-Gen SHBT-MAR-312X Twin  |
                                   +-------------------------------+
                                                   │
                 ┌─────────────────────────────────┼─────────────────────────────────┐
                 ▼                                 ▼                                 ▼
   ╭───────────────────────────╮     ╭───────────────────────────╮     ╭───────────────────────────╮
   │       [shbt-power]        │     │         [shbt-cf]         │     │         [shbt-qc]         │
   │  Commercial Fusion Grid   │     │  1,800-Module LANR Array  │     │ Bare-Metal Microkernel &  │
   │   (8,750 MW p-11B Twin)   │     │    & Thermal-Hydraulics   │     │   Photonic Quantum Bus    │
   ╰─────────────┬─────────────╯     ╰─────────────┬─────────────╯     ╰─────────────┬─────────────╯
                 │                                 │                                 │
                 └────────────────────────┬────────┴─────────────────────────────────┘
                                          ▼
   ╭──────────────────────────────────────────────────────────────────────────────────────────────╮
   │                                SPECIALIZED VEHICLE TWINS                                     │
   │                                                                                              │
   │  • shbt-ghost : Reactionless Propulsion & Local Gravity Wells (3+1 CCZ4 / PCSS Crowbars)     │
   │  • shbt-recon : Macroscopic State Translocation Gateway (Stinespring V_macro / 504 Gbps)     │
   │  • shbt-sglt  : Synthetic Gravitational Lensing Telescope (SE-L2 Swarm / TMSV Metrology)     │
   │  • shbt-warp  : Holographic Warp Metric & 3+1D Flight Twin (ADM α=1.0 / SHBT-MAR-312X Engine)│
   ╰──────────────────────────────────────────┬───────────────────────────────────────────────────╯
                                              │
                                              ▼
   ╭───────────────────────────────────────────────────────────────────────────────────────────╮
   │                                       shbt-exotic                                         │
   │                MULTI-PROTOCOL SPACETIME ENGINEERING CO-SIMULATION BENCH                   │
   │                                                                                           │
   │  • Cross-Protocol Field Coupling (Warp + Stasis + Translocation + Wells + Comms)          │
   │  • Global Energy Condition & Ford-Roman Quantum Inequality (QI) Dark-Ledger Auditing      │
   │  • Dynamic 5-Stage Multi-Technology Flight Director & Relativistic PDE Mesh Solvers       │
   ╰───────────────────────────────────────────────────────────────────────────────────────────╯

shbt-cf provides the solid-state low-energy nuclear reaction (LANR) baseload power, dynamic screening physics, and thermal-hydraulic models that support the initial deployment phases of sys1own/shbt-swarm:

1. 1,800-Module LANR Baseload Microgrid: shbt-cf's compact, reliable nickel-hydrogen cold fusion generator (999.054 kW net DC array @ 400 V) provides autonomous, unrefueled baseload power for shbt-swarm Phase 2 forward exploratory probes, communication relays, and initial planetary seed outposts.
2. Floquet–Adler–Wiser Dielectric Screening (U_eff = 352.48 eV): The non-equilibrium dynamic screening mechanisms and Bremsstrahlung suppression models developed in shbt-cf established the physical principles for non-thermal nuclear reaction control in matter-folding catalysts.
3. Dual-Stage Thermoelectric Generation (TEG) & Helium Thermal-Hydraulics: Dual-stage Seebeck enthalpy recovery manifolds and closed-loop supercritical helium cooling models provide reliable secondary heat reclamation across distributed swarm platforms.

In return, shbt-swarm closes the loop back into shbt-cf through Phase 3 universal 230-space-group matter folding (R_univ^SWARM), fabricating optimized superlattice geometries that enhance thermoelectric figure-of-merit and improve phonon extraction efficiency.

• Upstream Model & Ledger Integration: shbt-cf exports the 1,800-module Lattice-Assisted Nuclear Reaction (LANR) core specification (555.03 W net DC per cell, 999.054 kW net array) as continuous housekeeping power to sys1own/shbt-power for power-plant auxiliary cold-start bootstrapping, to sys1own/shbt-warp to service its boundary-emitter Landauer debt, and to sys1own/shbt-recon as the continuous baseload rail behind its pulsed ¹⁷⁸ᵐ²Hf graser bursts. It additionally exports the LANR ledger to sys1own/shbt-sglt to balance the 906.000 kW non-sheddable entropy debt of the deep-space synthetic gravitational lensing telescope array.
• Warp Drive Energy & Cryogenic Transfer: The 1,800-module LANR starter core (555.03 W net DC/cell, 999.054 kW net array) is additionally exported to sys1own/shbt-warp to balance the continuous 906.000 kW Landauer entropy debt of the boundary emitter array, maintaining a +93.054 kW raw electrical margin and +33.104 kW net continuous system surplus. shbt-cf further provides two-phase helium flow boiling and Kapitza thermal resistance models (ΔT_K = 3.546 K) ensuring ΔT ≥ 11.790 K cryogenic headroom below superconducting quench limits.
• Kinetics & Fast Recovery Transfers: The McNabb–Foster hydrogen transport and isotope trapping solvers (McnabbFosterSolver) developed in shbt-cf are imported by sys1own/shbt-ghost (ghost-lanr-interface) for long-term fuel retention, while shbt-cf integrates the sub-2.5 ns Photoconductive Semiconductor Switch (PCSS) optical trigger logic and 94.20% SiC inductive recovery shunts engineered in sys1own/shbt-ghost to harvest stored magnetic coil energy.
• Mathematical Foundations & Microkernel Contracts: Foundational MPFR arbitrary-precision arithmetic, condition-number tracking, and symplectic integrators (Yoshida6) are anchored in the canonical WZW affine branch models from sys1own/shbt-precision. Low-level register definitions and SECDED ECC standards derive from the bare-metal C11 shbt-os runtime in sys1own/shbt-qc, while boundary Conformal Field Theory state-vector formulations and the invariant rational dark ledger capacity partitioning (η_D = 23/33, η_A = 10/33) originate from sys1own/shbt-exotic. Finally, zero-copy lock-free POSIX SPSC circular ring buffers and 128-byte dual-cacheline C-ABI mapping conventions are shared with sys1own/shbt-recon.

The current baseline simulation state (sim_outputs/simulation_verification.json) validates complete physical closure and numerical convergence across all core physics modules, and audits the full 50-gate system verification matrix (GATE-01 – GATE-50, all PASS) defined by the cf3 engineering specification.

P_net = P_TEG − P_drive,net − P_aux = 1045.58 − 368.45 − 122.10 = +555.03 W > +550.00 W

• 3D Eulerian–Eulerian Thermal-Hydraulics (src/physics/thermal_hydraulics.rs): two-phase liquid/vapor conservation with Ishii–Zuber drag, Tomiyama lift, Antal–Frank wall lubrication, Burns turbulent dispersion, and RPI wall heat-flux partitioning with Hibiki–Ishii site density across 64 parallel OFHC-Cu micro-channels (D_h = 250 μm):

q"_tot = q"_1φ + q"_q + q"_e

• Dynamic Non-Equilibrium Deuterium Screening (src/physics/screening.rs): Fick–Soret transport with D₀ = 2.85×10^-7 m²/s, E_a = 0.224 eV, Q^* = 0.048 eV, driving 100 Hz recalculation of the 15,625 × 15,625 Floquet–Adler–Wiser dielectric matrix:

∂_t x = ∇ · [ D_D(T, x) (∇x + (x(1 − x) Q* / (k_B T²)) ∇T) ] + Ṡ_phase

ε_{G,G'} = δ_{G,G'} − v(q + G) χ⁰_{G,G'}(x)

• High-T_c Superconducting Coils + SiC Crowbar (src/physics/power.rs): thin-film MgB₂ micro-coils (T_c = 39.0 K) at 68.5 kHz with Bean critical-state and flux-flow losses (P_coil = 12.35 W); the SiC crowbar harvests the 16.62 mJ/cycle inductive energy at 94.20% efficiency, cutting parasitic drive from 422.22 W to 368.45 W.
• GST Optical Self-Healing (src/physics/optics_healing.rs): chalcogenide Ge₂Sb₂Te₅ buffer in the Pd_0.9132Ir_0.0868 grating stack; 50 ns electro-thermal pulses (27.9 mJ/cm²) trigger melt-quench recrystallization restoring R_a < 0.8 nm, A ≥ 98.40%, R_grating > 99.9%.
• C-ABI MMIO Map (crates/shbt-fabrication-hil/include/shbt_floquet_mmio.h): zero-copy 64-byte aligned shbt_mmio_control_t register map with 100 Hz HIL frame kernel (shbt_floquet_mmio_execute_frame).
• Active-Alloy McNabb–Foster Kinetics (src/physics/transport.rs): two-family diffusion and trapping model for Pd_0.9132Ir_0.0868D_x at x = 0.9132 with Soret thermophoresis and partial-molar-volume stress drift:

J_D = −D_D(T, x) ∇C_L + (V_H* / (R T)) D_D C_L ∇σ_h + (Q* / (R T²)) D_D C_L ∇T

with D₀ = 2.85×10^-7 m²/s, E_a = 0.224 eV, Q = 0.048 eV, V_H = 1.70×10^-6 m³/mol, dislocation traps N₁ = 1.50×10²⁴ m^-3 (E_t,1 = 0.23 eV) and grain-boundary traps N₂ = 5.00×10²³ m^-3 (E_t,2 = 0.15 eV).

• 3D Chaboche Thermoviscoplasticity & Joint FEA (src/physics/mechanics.rs, crates/shbt-fea-structural): two-term nonlinear kinematic hardening (C₁ = 45.2 GPa, γ₁ = 410, C₂ = 8.5 GPa, γ₂ = 62 at 298.15 K, linearly interpolated to the 623.15 K set) plus Voce isotropic hardening (R_∞ = 85 MPa, b = 12.46); VCCT energy release rates on the 3.5 μm Ni–Cu–Sn TLP bondline (G_IC = 25, G_IIC = 65, G_IIIC = 60 J/m², delamination factor f = 0.484); Morrow strain-life fatigue N_f ≈ 65,474 ≥ 52,400 cycles.
• Plant Stability Maps (src/physics/thermal_hydraulics.rs): system pump head H_pump(Q) = (65.0 − 1.25Q − 0.62Q²) kPa coupled to the two-phase core; Ledinegg excursive margin and Ishii–Zuber DWO ratio N_pch,crit = 1.45·N_sub + 2.5 evaluated across startup, nominal (ΔP = 42.8 kPa @ 4.85 L/min), surge, and low-flow regimes with Q < 3.95 L/min interlock.
• 3D Radiation Transport (src/physics/radiation.rs): 2.45 MeV DD neutrons plus secondary capture gammas (478 keV from ¹⁰B(n,α)⁷Li^*, 2.223 MeV from H(n,γ)D) through the H₂O / 316L / 5 wt% B-PE / Pb shielding stack; accessible surface dose (0.38 ± 0.012) μSv/h < 0.50 μSv/h at Q_N ≤ 10⁶ n/s.
• GUM Covariance Metrology (src/physics/metrology.rs, crates/shbt-metrology-gum): ISO/IEC 98-3 multi-variable input covariance Σ_X (8×8, correlations r_Tin,Tout = 0.85, r_Qt,QT = 0.92, r_Q,UA = 0.30) with analytic/dual-number Jacobian J, Σ_Y = JΣ_XJ^T, N = 10⁶ MC cross-check; propagated bounds u(P_thermal) ≈ 11.8 W, u(P_⁴He) ≈ 2.5×10^-9 mbar, u(FOM) ≈ 0.018 securing P_net = +555.03 W at 3σ.
• 128-byte C-ABI MMIO Record (crates/shbt-fabrication-hil/include/shbt_floquet_mmio.h): shbt_mmio_control_t widened to a dual-cacheline packed record — plastic_strain_max @ 0x0028, dose_surface_usv @ 0x0030, gpu_frame_count @ 0x0038, matrix_dim = 15625 @ 0x0040, d_floquet_ptr @ 0x0048, pad[16] @ 0x0050 — 64-byte aligned.

The core math engine links against system GMP and MPFR libraries. On Debian/Ubuntu distributions:

sudo apt-get install -y m4 libgmp-dev libmpfr-dev

Build, lint, and test all 9 crates across the workspace:

cargo check --workspace
cargo clippy --workspace --all-targets -- -D warnings
cargo test --workspace

Install the Python package in editable mode and run verification tools:

# Install Python orchestration layer
pip install -e python

# Run CI topology and self-tests
python3 scripts/check_dep_graph.py
python3 -m unittest discover -s python

# Inspect topology and run crate integration tests via shbt_cf CLI
PYTHONPATH=python python3 -m shbt_cf topology
PYTHONPATH=python python3 -m shbt_cf verify
PYTHONPATH=python python3 -m shbt_cf test shbt-core-math

To build and link the high-performance PyO3 native extension:

maturin build -m crates/shbt-py-bindings/Cargo.toml --release
maturin build -m crates/shbt-fabrication-hil/Cargo.toml --features python --release
pip install --user --no-deps target/wheels/shbt_py_bindings-*.whl
pip install --user --no-deps target/wheels/shbt_fabrication_hil-*.whl
python3 scripts/verify_zerocopy.py

To prevent non-deterministic floating-point drifting across microarchitectures and compiler toolchains:

1. Strict IEEE 754 Compliance: Fast-math flag sets are explicitly disabled across all workspace profiles (-C llvm-args=-enable-no-nans-fp-math=false -C llvm-args=-enable-no-signed-zeros-fp-math=false).
2. Fixed-Point Accumulation: Global energy and state accumulators utilize shbt_core_math::fixed::Q64x64 128-bit signed fixed-point integers (2^-64 fractional resolution).
3. Adaptive Precision Inversion: Ill-conditioned linear algebra operations dynamically select MPFR mantissa bit-widths via PrecisionMode::arbitrary_for_condition_number based on matrix condition number κ₁(A) and promotion thresholds (κ(A) · ε₆₄ ≥ θ_thresh).
4. Denormal Handling: Hardware denormals-are-zero (DAZ) and flush-to-zero (FTZ) modes are explicitly enforced via MXCSR (x86_64) and FPCR (aarch64).

Verification Gate: A double-well Duffing oscillator integrated for 10⁷ Yoshida6 steps maintains energy drift |Δ H / H₀ |< 10^-12 (measured ≈ 6 × 10^-14) while a 1-ulp shadow trajectory diverges to O(1).

Benchmark Gate: The 10 kHz telemetry channel pipeline processes over 61,000 frames across 3 synthetic streams with zero frame drops, zero mutex allocation, and sub-microsecond mean processing latency.

The native bindings layer (shbt-py-bindings and shbt-fabrication-hil) exports C-ABI functions and PyO3 native extensions:

• Zero-Copy Memory Access: Telemetry ring memory is exposed directly to Python using zero-copy array views. Pointer verification guarantees zero memory copy overhead (view.ctypes.data == ring.data_ptr()).
• Python Package (python/shbt_cf):
• workbench.py & workbench/: Parameter sweep engines, HUD models, and Coffin-Manson Monte Carlo wrappers.
• optimize.py: Pure-NumPy NSGA-III multi-objective optimization (Das-Dennis reference points, non-dominated sorting, polynomial mutation).
• hud_streamlit.py: Interactive Streamlit engineering interface.
• control/: System order reduction and state-space control routines.
• metrology/: Python-side telemetry analysis and calibration curve parsing.
• schemas/: JSON schemas for anisotropic permittivity, Pt100 calibration, and surface roughness topologies.

1. Physical Closure Audit (crates/shbt-fabrication-hil/tests/closure_audit.rs):

• Volume Mapping: Verifies active domain mapping ratios (V_domains / V_metal = 8000) according to cf.pdf Equations 29, 35, and 123.
• Calorimetric Uncertainty Budget: Validates GUM calorimetry Monte Carlo propagation (P = ṁ c_p Δ T + P_env), confirming standard uncertainty bounds (u_c ≈ 5.696 W < 5.70 W).
• Fatigue Life Ceilings: Evaluates Coffin-Manson strain-life predictions (ε'_f = 0.18, c = -0.62), verifying fatigue endurance bounds across cyclic plastic strain thresholds.

2. Dual-ISA Bitwise Regression (crates/shbt-fabrication-hil/tests/bitwise_regression.rs):

• Asserts exact bitwise kernel digest agreement (0x03de00d3acded967) across x86-64 (AVX-512) and ARM64 (NEON) architectures to prevent microarchitectural float discrepancies.

3. Multiphysics Integration Pipeline (tests/test_full_multiphysics_pipeline.py & tests/verification_tests.rs):

• Validates full end-to-end integration across Floquet dielectric matrix assembly, RCWA optical enhancement, structural FEA deformation, closed-loop thermal power ledger, and HIL telemetry streaming.