gaussian_pulse
Linearized telegraph-equation reference; finite-amplitude nonlinear M1 can have a model error floor. Energy excess uses a logarithmic display with a labeled floor; signed errors and fluxes remain linear.
Convergence becomes available after two resolutions are complete.
| N | E L1 | E L2 | E L∞ | p(L1) | p(L2) | p(L∞) |
|---|---|---|---|---|---|---|
| 32 | 1.627470e-06 | 2.164338e-06 | 6.541292e-06 | — | — | — |
32³ cells
Central cell layer at z=dx/2; ray at y=z=dx/2. Paired values CSV · Profiles PDF
Flux maps
Radiation conservation
R = Q(t) − Q(t_ref) + B(t) − S(t). B is cumulative outward transport; S is the signed source change. Flux integrals divided by c² give radiation momentum. The reference is the first recorded time (zero unless labeled as a resumed segment). All plotted totals and residuals use its fixed reference magnitude |Q(t_ref)|. Zero/negligible reference flux uses the labeled fixed scale c|E(t_ref)|. Neither denominator changes during the run.
| Field | Scale | Final R/scale | Final |R|/scale | Maximum |R|/scale |
|---|---|---|---|---|
| er | abs(initial_energy) | 0.000000e+00 | 0.000000e+00 | 1.667907e-16 |
| fx | c * abs(initial_energy) | -9.744115e-23 | 9.744115e-23 | 1.182696e-22 |
| fy | c * abs(initial_energy) | 2.205480e-23 | 2.205480e-23 | 1.664544e-22 |
| fz | c * abs(initial_energy) | -3.138153e-22 | 3.138153e-22 | 3.427853e-22 |
Budget summary CSV · Raw history CSV · Fraction history CSV
