HomeDocs-Technical WhitePaper36-EFT.WP.EDX.Current v1.0

Chapter 14 — Simulation Stack & Benchmark Cases (Methods.SimStack)


I. Chapter Objectives & Structure

  1. Objective: Build a multi-scale Simulation Stack (SimStack) that, under a unified dialect, executes the full forward → inversion → validation → release workflow from netlist/layout/binding to Z_eft(omega), T_arr, w_p, and kernels K_{·}; provide reproducible sim–meas alignment and a benchmark case library.
  2. Structure: Stack layers → Shared dialect & records → Numerical methods & stability → Benchmark cases → Sim–meas alignment → Falsifiability → Data structure & export → Compliance templates → Correspondence & degeneracy → Cross-chapter pointers & summary.
  3. Shared arrival dialect (two equivalent forms; path/measure explicit; record delta_form):
    • Constant-factored: T_arr = ( 1 / c_ref ) * ( ∫ n_eff d ell )
    • General: T_arr = ( ∫ ( n_eff / c_ref ) d ell )

II. Simulation Stack Layers (bottom-up)

  1. Physics Core (S20/S40/S50):
    • S20-*: continuity/power conservation & minimal effective conduction;
    • S40-*: causal kernel convolutions with path weighting;
    • S50-*: impedance mapping Z_eft = Z_ref + ΔZ_T (+ ΔZ_rad).
  2. Path Engine (I30 / Chapter 8):
    layout ↔ gamma(ell) binding; discretization {Δell_i, n_eff(i)}; multi-path weights w_p(omega).
  3. Boundary & Radiation:
    Positive-real radiation correction ΔZ_rad(omega) with Re{ΔZ_rad} ≥ 0 and K–K consistency.
  4. Solvers:
    Frequency-domain solvers (sparse linear / FFT convolution) and optional time-domain convolution; robust ∂Z/∂θ and sensitivity.
  5. UQ & Inversion (Mx-*):
    Prior/likelihood/evidence, HMC/NUTS or SMC, posterior predictive checks (PPC).
  6. I/O & Cards:
    Data/model/pipeline cards with mandatory arrival and qa_gates; versioning and reproducible seeds.

III. Shared Dialect & Records (mandatory)


IV. Numerical Methods & Stability

  1. Frequency-domain kernel convolution: J_T = ( K_s ⊛ ∇T_fil ) + ( K_t ⊛ ∂T_fil/∂t ) → multiplication in frequency with windowing/band-limits; apply smooth/band-limited priors to K_{·}(ω); forbid right-half-plane poles.
  2. Discrete arrival time:
    • n_over_c: T_arr ≈ (1/c_ref) * Σ_i n_eff(i) · Δell_i
    • one_over_c_times_n: T_arr ≈ Σ_i ( n_eff(i) / c_ref ) · Δell_i
  3. Grid & bandwidth: choose N_ω and B_ω by the coherence window and target phase linearity error E_phase; recommended numeric gate: small | d^2 (arg Z) / dω^2 | < κ_phase.
  4. Stability tactics: whitened residual spectral flatness checks; regularize or re-parameterize ill-conditioned sensitivities (e.g., log σ_eff).
  5. Parallelism & caching: parallelize over paths/frequencies; cache kernels/weights/geometry to avoid repeated convolutions.

V. Benchmark Case Library (minimal, reproducible)

Each case ships netlist/layout/binding, arrival, qa_gates, target gates, and reference outputs for regression.


VI. Simulation–Measurement Alignment Workflow (I30-3 / M10)


VII. Falsifiability Criteria (simulation side)


VIII. Data Structure & Export (minimal template)

simstack:

case_id: "SIM-01"

model_id: "EDX-Current-eft-ms"

freq_grid_Hz: [...]

layout_ref: "LAY-2025-001"

binding_ref: "LAY2PATH-xxxx"

deemb: {method:"TRL", version:"1.2"}

sync: {dt_sync_s: 2.0e-12}

arrival:

form: "n_over_c" # or "one_over_c_times_n"

gamma: "explicit"

measure: "d_ell"

c_ref: 299792458.0

Tarr_s: 1.234e-09

u_Tarr_s: 6.0e-12

outputs:

Z_eft: {real: [...], imag: [...]}

argZ: [...]

weights: {w_main: [...], w_side: [...]}

qa_gates: {check_dim:"pass", passivity:"pass", KK:"pass"}

seed: 20250915


IX. Compliance Templates (copy-ready)

api:

- id: "SimStack.build"

proto: "build(netlist, layout, binding, options) -> sim_handle"

- id: "SimStack.forward"

proto: "forward(sim_handle, theta, freq_grid) -> {Z_eft, T_arr_p, w_p}"

- id: "SimStack.invert"

proto: "invert(sim_handle, data, priors) -> posterior, logZ"

- id: "SimStack.ppc"

proto: "ppc(sim_handle, posterior) -> residual_spectrum, gates"

- id: "SimStack.export"

proto: "export(sim_handle, format:'cards|json') -> artifact"

phi = argZ[ω1:ω2] # coherence window

T_group = grad(phi, omega)

E_phase = max_abs(phi - (omega*Tarr + phi0_opt))

GDR = max_abs(T_group - median(T_group))

assert E_phase <= E_phase_gate and GDR <= GDR_gate


X. Correspondence & Degeneracy to the Classical Framework
With K_s = K_t = 0, w_p = 1, and ΔZ_rad = 0, the stack reduces to classical RLC/telegrapher solutions; the added value here is explicit paths/arrival/weights/positive-real radiation plus a unified metrology → inversion → evidence → release chain.


XI. Cross-Chapter Pointers & Summary


Copyright & License (CC BY 4.0)

Copyright: Unless otherwise noted, the copyright of “Energy Filament Theory” (text, charts, illustrations, symbols, and formulas) belongs to the author “Guanglin Tu”.
License: This work is licensed under the Creative Commons Attribution 4.0 International (CC BY 4.0). You may copy, redistribute, excerpt, adapt, and share for commercial or non‑commercial purposes with proper attribution.
Suggested attribution: Author: “Guanglin Tu”; Work: “Energy Filament Theory”; Source: energyfilament.org; License: CC BY 4.0.

First published: 2025-11-11|Current version:v5.1
License link:https://creativecommons.org/licenses/by/4.0/