Laboratory · Historical model comparison · Explanatory preview

Does the wave equation keep its form?

Apply two coordinate substitutions to the same forward-travelling plane wave and compare the analytic residuals.

Change one setting and compare the consequences

Static worked example
Model inputs · SI units
Predict first

Your next accepted calculation shows the before-and-after values beside this note. No score or answer gate is used. The note stays in this tab and is not uploaded or saved.

Supported range: -0.95 to 0.95. Scientific notation is accepted: type 10⁻⁹ as 1e-9.
Experiment settings angular wavenumber k
Supported range: 0.001 to 1000000 rad/m. Scientific notation is accepted: type 10⁻⁹ as 1e-9.

Changes here apply with Apply settings.

Residual normalized by k² · one common linear scale
Galilean substitution0.19
Lorentz substitution0

Bar range: −0.19 to +0.19. The center line is zero. The scale is recalculated after Apply.

Accepted calculation 0. Every model row shares these inputs.

Inputs used for these results
Signed frame speed divided by c
0.1
Angular wavenumber k
1 rad/m
Computed values · shared inputs · not experimental observations
ObservableGalilean substitutionLorentz substitution
Residual normalized by k²0.190
Maximum absolute residual (1/m²)0.190
Mixed-derivative coefficient (s/m)6.671 × 10⁻¹⁰0

This checks the same forward-travelling plane wave after each coordinate substitution, not a measured wave or a grid simulation. At zero relative speed both residuals vanish. The Lorentz zero is an analytic identity under the stated map; it does not independently prove a field transformation or select an interpretation of space and time. Reversing the frame velocity while holding the forward wave fixed is a different relative configuration.

A light wave obeys a single wave equation. Describe the same wave from a moving frame with the everyday change of coordinates and the equation no longer holds in its old form; use the change of coordinates Einstein derived, and it does.

Take a wave travelling at c along x, φ = cos(k(x − ct)), which satisfies ∂2φ/∂x2 − (1/c2)∂2φ/∂t2 = 0. Describe it from a frame moving at v. With the Galilean substitution x′ = x − vt, t′ = t, the wave no longer satisfies an equation of the old form in the new coordinates: it leaves a residual of relative size |β(2 − β)|, 0.19 at β = 0.1, and the correct operator there is a different one, (1 − β2)∂2/∂x′2 + (2v/c2)∂2/∂x′∂t′ − (1/c2)∂2/∂t′2. With the transformation of §3 of the relativity paper the equation keeps its form and the residual is 0 at every speed. Section 6 shows the same for the full Maxwell–Hertz equations of empty space, with the fields transforming as well; this instrument checks only the scalar wave operator.

Modern qualifications

Voigt noticed in 1887 that a transformation of this kind leaves the wave equation unchanged, and Lorentz and Poincaré worked with the full set for electrodynamics by 1904 and 1905. Einstein's §6 applies his transformation to the Maxwell–Hertz equations for empty space and finds them unchanged in form, and §7 draws the Doppler and aberration formulas from that. Before 1905 the invariance was a property of the equations; in the paper it follows from two principles about measurement.

These are host reference calculations with modern SI calibration and illustrative settings. They are not historical measurements, a historical dataset, or publication of the strict 1904 shelf. No FrankenSim WASM calculation is claimed.

What is being calculated?

The diagnostic is the scalar forward-travelling wave cos(k(x − ct)) and the vacuum operator ∂²/∂x² − c⁻²∂²/∂t². The six-component Maxwell–Hertz field equations are not modeled: this residual belongs to the scalar wave alone. The largest absolute residual is the analytic amplitude over phase, not a sampled numerical maximum.

x′=x−vt,t′=t⟹(1−β2)∂x′2+2vc2∂x′∂t′−1c2∂t′2x'=x-vt,\quad t'=t\quad\Longrightarrow\quad (1-\beta^2)\partial_{x'}^2+\frac{2v}{c^2}\partial_{x'}\partial_{t'}-\frac{1}{c^2}\partial_{t'}^2
max⁡∣□′ϕ∣k2=∣β(2−β)∣(Galilean),0(Lorentz)\frac{\max|\square'\phi|}{k^2}=|\beta(2-\beta)|\quad\text{(Galilean)},\qquad 0\quad\text{(Lorentz)}

Hold the frame speed fixed and double k. The absolute residual changes while the normalized one does not. Then set the frame speed to zero. A preserved wave operator alone does not fix the physical meaning or normalization of all transformed fields.

Continue to the Maxwell–Hertz field-equation laboratory · Construct the coordinate map from its constraints

Prediction is not a measurement

No observed points, digitized fringe shifts or experimental confidence bounds have been added to these plots. A source-specific dataset, with its geometry, wavelength, protocol and uncertainty, is still needed before making a numerical comparison with a historical experiment.

The reference implementation is src/physics/reference/shelfOptics.ts. Its functions own the arm times, fringe shifts, drag speeds and wave residuals; the interface only projects their results. Numerical precision in a table does not imply measurement accuracy.

Inspect the reference calculations

This comparison belongs to the special-relativity discovery route.

Embed this laboratory on another page. It starts from its worked defaults, not your current settings.