M1 RK1 and FLD

Movies

19 movies available from 19 completed runs. Open the movie gallery.

Saved states are shown with zero-order hold and no solution interpolation.

Streaming Wave · L2 · M1 RK1

Streaming Wave · L2 · FLD backward Euler

Streaming Front · L2 · M1 RK1

Streaming Front · L2 · FLD backward Euler

Gaussian Pulse · L2 · M1 RK1

Gaussian Pulse · L2 · FLD backward Euler

Equilibrium Sphere · L2 · M1 RK1

Equilibrium Sphere · L2 · FLD backward Euler

Subcritical Shock · L2 · M1 RK1

Subcritical Shock · L2 · FLD backward Euler

Supercritical Shock · L2 · M1 RK1

Supercritical Shock · L2 · FLD backward Euler

Radiating Blast · L2 · M1 RK1

Radiating Blast · L2 · FLD backward Euler

Marshak · L2 · M1 RK1

Marshak · L2 · FLD backward Euler

Streaming Wave · L3 · M1 RK1

Streaming Wave · L3 · FLD backward Euler

Streaming Front · L3 · M1 RK1

Status: ready. 48 uniform-grid runs: eight problems, both models, levels 2/3/4 (32/64/128 cells per side). 121 requested physical snapshots and 60-second movies. Positive streaming optical depth 1e-4 for both models. Full shock/blast durations, longer Gaussian/Marshak evolution. Preserve validated samples and initial/final 3D states; retire intermediate 3D files only after verified extraction. First-order moving FLD uses a serial root-gathered composite solve and remains subject to memory limits. Full-resolution runtime and timestep convergence are unverified. Start only on explicit command. Model discrepancies are not FLD error estimates. FLD uses a serial root-gathered multigrid solve. No moving-gas or scaling claim.

M1 evolves radiation energy and flux using RK1 with radiation substeps. FLD solves energy and material heating together implicitly; its flux is a derived diagnostic. Cross-model differences are not truncation-error estimates.

Model comparison plots

Problem	Model / integrator	Status	Final time	Maximum timestep
streaming_wave · L2	M1 RK1	completed	0.2	0.001
streaming_wave · L2	FLD backward Euler	completed	0.2	0.001
streaming_front · L2	M1 RK1	completed	0.2	0.001
streaming_front · L2	FLD backward Euler	completed	0.2	0.001
gaussian_pulse · L2	M1 RK1	completed	1.5	0.00625
gaussian_pulse · L2	FLD backward Euler	completed	1.5	0.00625
equilibrium_sphere · L2	M1 RK1	completed	0.2	0.001
equilibrium_sphere · L2	FLD backward Euler	completed	0.2	0.001
subcritical_shock · L2	M1 RK1	completed	38300	191.5
subcritical_shock · L2	FLD backward Euler	completed	38300	191.5
supercritical_shock · L2	M1 RK1	completed	7500	37.5
supercritical_shock · L2	FLD backward Euler	completed	7500	37.5
radiating_blast · L2	M1 RK1	completed	22000	110
radiating_blast · L2	FLD backward Euler	completed	22000	110
marshak · L2	M1 RK1	completed	60	0.05
marshak · L2	FLD backward Euler	completed	60	0.05
streaming_wave · L3	M1 RK1	completed	0.2	0.001
streaming_wave · L3	FLD backward Euler	completed	0.2	0.001
streaming_front · L3	M1 RK1	completed	0.2	0.001
streaming_front · L3	FLD backward Euler	pending	0.2	0.001
gaussian_pulse · L3	M1 RK1	pending	1.5	0.00625
gaussian_pulse · L3	FLD backward Euler	pending	1.5	0.00625
equilibrium_sphere · L3	M1 RK1	pending	0.2	0.001
equilibrium_sphere · L3	FLD backward Euler	pending	0.2	0.001
subcritical_shock · L3	M1 RK1	pending	38300	191.5
subcritical_shock · L3	FLD backward Euler	pending	38300	191.5
supercritical_shock · L3	M1 RK1	pending	7500	37.5
supercritical_shock · L3	FLD backward Euler	pending	7500	37.5
radiating_blast · L3	M1 RK1	pending	22000	110
radiating_blast · L3	FLD backward Euler	pending	22000	110
marshak · L3	M1 RK1	pending	60	0.05
marshak · L3	FLD backward Euler	pending	60	0.05
streaming_wave · L4	M1 RK1	pending	0.2	0.001
streaming_wave · L4	FLD backward Euler	pending	0.2	0.001
streaming_front · L4	M1 RK1	pending	0.2	0.001
streaming_front · L4	FLD backward Euler	pending	0.2	0.001
gaussian_pulse · L4	M1 RK1	pending	1.5	0.00625
gaussian_pulse · L4	FLD backward Euler	pending	1.5	0.00625
equilibrium_sphere · L4	M1 RK1	pending	0.2	0.001
equilibrium_sphere · L4	FLD backward Euler	pending	0.2	0.001
subcritical_shock · L4	M1 RK1	pending	38300	191.5
subcritical_shock · L4	FLD backward Euler	pending	38300	191.5
supercritical_shock · L4	M1 RK1	pending	7500	37.5
supercritical_shock · L4	FLD backward Euler	pending	7500	37.5
radiating_blast · L4	M1 RK1	pending	22000	110
radiating_blast · L4	FLD backward Euler	pending	22000	110
marshak · L4	M1 RK1	pending	60	0.05
marshak · L4	FLD backward Euler	pending	60	0.05

Time is in seconds except for normalized Marshak time. Gaussian, equilibrium sphere and Marshak have analytic diffusion-limit references for FLD energy and flux; Marshak also has material energy. These are not exact nonlinear LP FLD solutions. Thin positive-opacity streaming cases compare closures and conservation without claiming an analytic FLD profile. Reference discrepancies remain separate from cross-model differences. Prepared runs are not validation evidence.

Supported scope

All eight problems are planned for both models at physical light speed on CPU. FLD uses a mixed-frame moving-gas update for shocks/blast. Streaming opacity is explicitly positive, and zero-opacity advection curves are limit references. Any legacy exclusions are listed here:

Plan and provenance · Prepared commands

Execution requires a separate explicit start. This workflow does not publish to a remote website.