aidoesscience
aidoessciencefindings › Eddy braking
ValidatingOracle-validated

Lenz's law weighed ring by ring

Drop a magnet down a copper tube and it crawls, nothing touching it. Is the brake really the F = −b·v of quantitative Lenz's-law theory — with b = (45/1024)·μ₀²m²σw/a⁴ for a point dipole in a thin tube — and can that coefficient, the terminal velocity, and even the tube's conductivity be recovered from nothing but the dipole's magnetostatic identity and Ohm's law?

Measured by the lab
0.043945312
Known value
0.043945313
Relative error
1.34e-8

▶ Run this simulationRead how it works

The finding

Lenz's law weighed ring by ring — the Levin–Rojo eddy-brake coefficient 45/1024 emerges from raw dipole-field + Ohm's-law ring sums to 1.3e-8, the fall saturates on v_term = 0.1984 m/s, and 24 noisy terminal-velocity reads invert to copper's σ = 5.961e7 ± 0.015e7 S/m (0.07 SE); Lenz-less fall and matched-asymptote quadratic drag both rejected

Method

Headless first-principles generator with NO closed-form drag: the tube is a stack of rings; each ring's flux is the line integral of the dipole vector potential Φ = 2πa·A_φ(a,ζ) (the faraday oracle's trick), its EMF a central finite difference of that flux as the magnet moves, its current EMF/R with R = 2πa/(σwΔz). The drag is assembled two independent ways — energy route Σ EMF²/R / v and Lorentz route Σ I·2πa·B_ρ — which must agree because ∇×A = B. The dimensionless Ĉ = βa⁴/(μ₀²m²σw) is scored against 45/1024 (loaded only to score). RK4 on M dv̇ = Mg − βv gives the saturating fall; 24 seeds of ±2% uniform reading noise on v_term (the module's noise model) invert to σ through the ring-sum geometry factor β/σ. Perturbations: β ∝ m² and β ∝ a⁻⁴ as free log-log fits of full ring-sum re-runs. Rivals through the same machinery: Lenz-less free fall (no terminal velocity) and a quadratic drag tuned to the same asymptote, killed by the approach-shape statistic t50/t90.

The law it recovers

F = C*mu0^2 m^2 sigma w v / a^4, C = 45/1024; M dv/dt = Mg - beta*v => v_term = Mg/beta, v(t) = v_term(1-e^{-t/tau})

Measurements, controls & cross-checks

Route agreement

Energy vs lorentz rel
9.2800e-12
Note
Joule-power route and ring-current Lorentz-force route agree to 9e-12 — energy conservation and Newton's third law close numerically

Convergence

Double rings 1p5x extent shift
1.2600e-8

Terminal velocity

V term m per s
0.1984046
Tau s
0.02023
Beta kg per s
0.59313
Saturation residual
1.1300e-7
Exp slope rel dev
5.2300e-12

Sigma recovery

Recovered S per m
5.9611e+7
Se S per m
152000
Known S per m
5.9600e+7
Mean rel error
0.000187
Se units
0.07
Worst seed rel
0.0192
Seeds
24

Perturbation

M loglog slope
2
M slope abs dev
4.4000e-16
A loglog slope
-4
A slope abs dev
7.5000e-11

Controls

No lenz
Exit v m per s
4.429
Ratio to v term
22
Note
rejected — induction ignored, the magnet exits the 1 m pipe at 22x the braked terminal velocity, still accelerating at g while the braked magnet's residual acceleration is 1.1e-7 g
Quadratic drag
T50 over t90
0.37311
Linear prediction
0.30103
Sim value
0.30103
Note
rejected structurally — tuned to the SAME v_term, but the tanh approach shape misses the exponential t50/t90 by 0.072 while the linear-drag trajectory lands on ln2/ln10 to 8.6e-7

Module systematics

DISCLOSED + PRICED: the module's on-screen σ recovery is circular by construction — EddyBrakeModule._recover() computes b from the closed form with C_DRAG = 45/1024 and inverts the same formula — and the certification now closes the loop the old finding left open: on the SAME measured v_term the closed-form inversion and this oracle's non-circular ring-sum instrument differ by exactly the Ĉ-recovery gap C_module/Ĉ − 1 = +1.34e-8 (gate L reconciles them to |Δ| < 1e-10), so the display pipeline is the ring-sum instrument to 1.3e-8 and everything else in the shown +1.19% gap is the pinned seed's single noise draw: σ_shown·(1+ε) = σ_true is an EXACT float identity (resid 2.2e-16) with ε = −1.180% the one ±2% uniform draw at seed 0xed2b, located at z = +1.06 inside the module's own 24-seed ensemble. The dimensionless side is also priced: the chart legend claims v(t) = v_term(1 − e^{−t/τ}) but the drawn trace is semi-implicit Euler at dt = 1/240 whose decay rate is exactly ln(1−dt/τ)/(−dt) = 1.00959/τ — +0.96% faster than the claimed exponential, invisible at chart scale, measured from the executed trace to 3.8e-15; the v_term asymptote is the exact float fixed point 8×0.22 = 1.76 (G − v_term/τ === 0 in doubles).

Certification

Rung
validated + honest-module
Gates
H–M execute the shipped EddyBrakeModule.ts (sha-pinned, mechanically type-stripped, headless Babylon/DOM stubs)
Frozen measurement
the module's only RNG is an internal mulberry32 (seed 0xed2b, zero Math.random) consumed in EXACTLY ONE draw inside init-time _recover(); fixedUpdate/render draw nothing and the seed is proven immutable across the whole lockstep — the shown σ can never drift
Init bit pin
executed _vTermReal/_sigmaRec/_sigmaErr === independently-coded replica bit-for-bit; all 13 statics pinned to hardcoded doubles AND to the oracle's inputs; 6-mesh scene census, DOM written once, camera pose pinned
Lockstep
1200 engine calls at FIXED_DT = fl(1/120) = exactly 2 substeps of fl(1/240) (accumulator returns to EXACTLY 0 every call): t, both heights, both velocities, freeExited and the FULL 320-sample v(t) ring buffer bit-exact at every call; reset branch exercised once (call ~640), 1760 trace wraps, free-magnet exit seen
Display
HUD + status + chart svg === replica templates bit-equal at all 200 %6 writes; live v-lines change while the σ block is bit-frozen; all three status branches and both chart branches exercised
Display priced
Shown sigma S per m
6.0312e+7
Shown vs true rel
0.011939112
Noise draw eps
-0.011798251
Identity resid
2.2000e-16
Z in own 24 seed ensemble
1.06
Instrument gap closed form vs ring sum
1.3400e-8
Euler rate excess vs claimed exponential
0.00959
Overclaim fixed
12th overclaim caught by execution (zero numeric change): the class doc claimed the recovery 'forward-models the fall' (it never integrates a fall — v_term is closed-form both ways) and lands 'σ ≈ 5.96×10⁷' (the shipped seed lands 6.03×10⁷, +1.19%); wording now states the single-draw circular pipeline and the actual landed value, and the HUD already printed the error honestly
Tampers
EDDY_TAMPER=sha → only gate H; EDDY_TAMPER=ulp (post-strip 1-ulp M_DIPOLE) → only the doubles-pinned gate I (m² cancels in σ analytically but not bit-wise; every display string absorbs 1e-16); hand-edited known_value → only scoring gates G1/G5 with the recovery and the whole cert unchanged
Derisk
16/16 gates in ~0.6 s

What it reduces to

The quantitative eddy-brake law of Levin, da Silveira & Rizzato (Am. J. Phys. 74, 815 (2006)) — F = (45/1024)μ₀²m²σwv/a⁴ for a point dipole on the axis of a thin conducting tube — itself the textbook completion of MacLatchy, Backman & Bogan's magnetic-braking experiment (Am. J. Phys. 61, 1096 (1993)), resting on Faraday induction + Lenz's law + Ohmic dissipation. Non-circular in the estimator sense: the generator codes only the dipole's A_φ and B_ρ (its magnetostatic identity, Jackson §5.6) and per-ring Ohm's law; 45/1024, the linearity in v, the m² and a⁻⁴ scalings, the exponential approach and the terminal velocity all emerge from the ring sum, and the known C is loaded only to score. Companion to ?world=faraday (same dipole primitives, EMF side) and ?world=transformer (induction with two coils); this world closes the induction arc with the force side — the current Faraday's law creates reacting back through the Lorentz force.

Confidence & reproduction

Confidence
high
Validation
derisk-pass
Re-run the check
npm run derisk -- eddy (scripts/eddy-derisk.mjs)
Oracle
scripts/oracles/eddy.reference.json

Sources

Y. Levin, F. L. da Silveira, F. B. Rizzato, Am. J. Phys. 74, 815 (2006); C. S. MacLatchy, P. Backman, L. Bogan, Am. J. Phys. 61, 1096 (1993); copper σ = 5.96×10⁷ S/m (annealed, 20 °C).

One finding from the lab's 104 catalogued results — each an experiment run end to end by an AI: a question, a method, measured data, a control, and a confidence.