Marshak wave · Octo-TIGER

The experiment

Physical setup

The normalized cube spans −1 ≤ x,y,z ≤ 1. The material density and momentum are fixed; radiation exchanges energy with the material through a linear equilibrium law E_eq = e. The propagation speed is c = 1, absorption and extinction are both χ = 100, and the default final time is t = 1. Equal material and radiation energy capacities correspond to ε = 1. Gravity and hydrodynamic transport are disabled.

This is a heating test in normalized units, rather than an ideal-gas temperature benchmark. Plots show E/4 and e/4. The saved metadata and source record the actual boundary treatment and reference quadrature.

This case uses normalized units with c = 1, opacity 100, and equal material/radiation energy capacities. Hydrodynamics and gravity are disabled. The default cube spans −1 to +1, and the default final time is 1. Adaptive runs retain their actual leaf spacing in exported data.

The comparison

Comparison and interpretation

The Su–Olson reference is the solution of a linear, gray, nonequilibrium diffusion problem in a semi-infinite stationary medium:

∂_tE − (1/3χ)∂_x²E = χ(e − E)
∂_te = χ(E − e)

With ξ = √3χ(x + 1), τ = χt, u = E/4, and v = e/4, the reference has u − (2/√3)∂_ξu = 1 at the illuminated boundary and zero initial u and v. The existing integral quadrature supplies radiation energy, material energy, and diffusion flux. Reference values are evaluated at cell centers; the exported comparison does not apply an additional cell-volume average.

What this test measures

Inspect the penetration depth and delay between radiation and material heating. The full-volume L1, L2, and maximum norms include material energy as well as all four radiation moments. Nonlinear M1 has finite propagation speed and a closure different from diffusion near the illuminated boundary and early front. Consequently these norms include model and boundary differences; they are not an exact-M1 order test. A zero transverse-flux error has no logarithmic convergence order. The numerical boundary combines a prescribed incoming half-isotropic bath with outgoing moments from the interior M1 intensity. Its finite-grid ghost-cell/HLL realization and early causal transient can differ from the diffusion boundary solution.

Measured results

All resolutions

Level	Grid	Batch status
2	—³	Ready
3	—³	Ready

These runs have different comparison settings or duplicate resolutions; convergence orders are not computed.

Full-volume norms and orders

Grid	Field	Norm	Error	Order	Error units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	er	L1	5.825905e-02	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	er	L2	1.963974e-01	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	er	Linf	8.907747e-01	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	fx	L1	4.134736e-03	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	fx	L2	1.438703e-02	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	fx	Linf	7.059959e-02	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	fy	L1	0.000000e+00	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	fy	L2	0.000000e+00	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	fy	Linf	0.000000e+00	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	fz	L1	0.000000e+00	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	fz	L2	0.000000e+00	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	fz	Linf	0.000000e+00	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	material_egas	L1	5.796755e-02	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	material_egas	L2	1.956904e-01	—	code units
AMR levels 1–2; 18432 active cells; finest equivalent 32³	material_egas	Linf	8.869638e-01	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	er	L1	2.876046e-02	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	er	L2	1.117458e-01	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	er	Linf	6.246538e-01	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	fx	L1	2.266185e-03	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	fx	L2	1.044080e-02	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	fx	Linf	7.504904e-02	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	fy	L1	0.000000e+00	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	fy	L2	0.000000e+00	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	fy	Linf	0.000000e+00	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	fz	L1	0.000000e+00	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	fz	L2	0.000000e+00	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	fz	Linf	0.000000e+00	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	material_egas	L1	2.864960e-02	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	material_egas	L2	1.115211e-01	—	code units
AMR levels 2–3; 90112 active cells; finest equivalent 64³	material_egas	Linf	6.247201e-01	—	code units

This problem: norms/orders CSV JSON

AMR levels 1–2; 18432 active cells; finest equivalent 32³ / level 2

Recorded setup and run provenance

Origin	Octo-TIGER application
Domain	Cube of side 2.000000e+00 (code length); centered at the origin
Resolution	AMR levels 1–2; 18432 active cells; finest equivalent 32³; finest Δx = 6.250000e-02 (code length)
Slice geometry	AMR: positive-side cells touching z=0, projected as constant leaf averages. Centers have z=leaf_dx/2; the ray similarly touches y=0. This is not a common cell-center plane. Raw CSV preserves native cells; display CSV repeats coarse values on finest pixels, retaining source dx/z. Full-volume norms use actual leaf volumes; uniform-grid convergence orders are not reported.
Final comparison time	1.000000e+00 (code time)
Radiation integration	Legacy first-order radiation update
Radiation backend	cpu
Physical light speed	1.000000e+00 (code speed)
Reduced-speed factor α	1.000000e+00
Radiation transport speed	1.000000e+00 (code speed)
Periodic boundaries	off
Hydrodynamics	off
Gravity	off
Implicit radiation source	on
Uniform grid	off
CFL	0.4
Timestep cap	0.01 (code time)
Output interval	0.050000000000000003 (code time)

{
  "adaptive": true,
  "background": 0.0,
  "c": 1.0,
  "case": "marshak",
  "cells": 32,
  "command": [
    "/home/dmarce1/workspace/octotiger/radiation_results/results/local-euler-ader2-20261011-restart-01/.campaign/launch-solver",
    "/home/dmarce1/workspace/octotiger/radiation_results/results/local-euler-ader2-20261011-restart-01/.campaign/runtime/octotiger",
    "--config_file=/home/dmarce1/workspace/octotiger/radiation_results/results/local-euler-ader2-20261011-restart-01/runs/l2/marshak/euler/marshak/l2/run.ini",
    "--hpx:threads=12"
  ],
  "comparison_signature": {
    "atomic_mass": "1",
    "atomic_number": "1",
    "cfl": "0.4",
    "code_to_cm": "1",
    "code_to_g": "1",
    "code_to_s": "1",
    "correct_am_hydro": "off",
    "disable_analytic": "off",
    "disable_diagnostics": "on",
    "eos": "IDEAL",
    "gravity": "off",
    "hard_dt": "0.01",
    "hydro": "off",
    "hydro_device_kernel_type": "OFF",
    "hydro_host_kernel_type": "LEGACY",
    "monopole_device_kernel_type": "OFF",
    "multipole_device_kernel_type": "OFF",
    "n_species": "1",
    "omega": "0",
    "periodic": "off",
    "problem": "MARSHAK",
    "rad_implicit": "on",
    "radiation": "on",
    "radiation_integrator": "euler",
    "radiation_reduced_c_factor": "1",
    "stop_step": "1000000",
    "stop_time": "1",
    "study_config_sha256": "a7ffc569ea95d062e7f4d4b96786783a2c678d14b32981f7e48f1a25899973e6",
    "unigrid": "off",
    "xscale": "1"
  },
  "config": {
    "atomic_mass": "1",
    "atomic_number": "1",
    "cfl": "0.4",
    "code_to_cm": "1",
    "code_to_g": "1",
    "code_to_s": "1",
    "correct_am_hydro": "off",
    "datadir": "/home/dmarce1/workspace/octotiger/radiation_results/results/local-euler-ader2-20261011-restart-01/runs/l2/marshak/euler/marshak/l2/",
    "disable_analytic": "off",
    "disable_diagnostics": "on",
    "disable_output": "off",
    "eos": "IDEAL",
    "gradient_refine_buffer": "0",
    "gradient_refine_coarsen_factor": "0.80000000000000004",
    "gradient_refine_field": "radiation_energy",
    "gradient_refine_floor": "9.9999999999999995e-21",
    "gradient_refine_mode": "absolute",
    "gradient_refine_threshold": "4",
    "gravity": "off",
    "hard_dt": "0.01",
    "hydro": "off",
    "hydro_device_kernel_type": "OFF",
    "hydro_host_kernel_type": "LEGACY",
    "max_level": "2",
    "min_level": "1",
    "monopole_device_kernel_type": "OFF",
    "multipole_device_kernel_type": "OFF",
    "n_species": "1",
    "odt": "0.050000000000000003",
    "omega": "0",
    "periodic": "off",
    "problem": "MARSHAK",
    "rad_implicit": "on",
    "radiation": "on",
    "radiation_backend": "cpu",
    "radiation_integrator": "euler",
    "radiation_reduced_c_factor": "1",
    "stop_step": "1000000",
    "stop_time": "1",
    "unigrid": "off",
    "xscale": "1"
  },
  "conservation": [
    {
      "boundary": -2.1342175106414416,
      "boundary_fraction": -0.06669429720754505,
      "field": "er",
      "final": 1.0701993626107305,
      "final_fraction": 0.03344373008158533,
      "initial": 0.0,
      "initial_fraction": 0.0,
      "integral_units": "code units",
      "max_abs_residual": 8.881784197001252e-16,
      "max_normalized_error": 2.7755575615628914e-17,
      "normalization_basis": "bath_energy * volume",
      "normalization_scale": 32.0,
      "normalized_error": 2.7755575615628914e-17,
      "normalized_residual": -2.7755575615628914e-17,
      "raw_change": 1.0701993626107305,
      "residual": -8.881784197001252e-16,
      "source": -1.0640181480307103,
      "source_fraction": -0.033250567125959696,
      "zero_initial_flux_relative_tolerance": 1.4210854715202004e-14
    },
    {
      "boundary": -4.242203394003569,
      "boundary_fraction": -0.13256885606261154,
      "field": "fx",
      "final": 0.04599314944779051,
      "final_fraction": 0.0014372859202434534,
      "initial": 0.0,
      "initial_fraction": 0.0,
      "integral_units": "code units",
      "max_abs_residual": 4.440892098500626e-16,
      "max_normalized_error": 1.3877787807814457e-17,
      "normalization_basis": "c * bath_energy * volume",
      "normalization_scale": 32.0,
      "normalized_error": 0.0,
      "normalized_residual": 0.0,
      "raw_change": 0.04599314944779051,
      "residual": 0.0,
      "source": -4.196210244555779,
      "source_fraction": -0.1311315701423681,
      "zero_initial_flux_relative_tolerance": 1.4210854715202004e-14
    },
    {
      "boundary": -2.056145933965384e-22,
      "boundary_fraction": -6.425456043641825e-24,
      "field": "fy",
      "final": 0.0,
      "final_fraction": 0.0,
      "initial": 0.0,
      "initial_fraction": 0.0,
      "integral_units": "code units",
      "max_abs_residual": 2.7257675730861304e-22,
      "max_normalized_error": 8.518023665894157e-24,
      "normalization_basis": "c * bath_energy * volume",
      "normalization_scale": 32.0,
      "normalized_error": 6.425456043641825e-24,
      "normalized_residual": -6.425456043641825e-24,
      "raw_change": 0.0,
      "residual": -2.056145933965384e-22,
      "source": 0.0,
      "source_fraction": 0.0,
      "zero_initial_flux_relative_tolerance": 1.4210854715202004e-14
    },
    {
      "boundary": -2.320843729982353e-22,
      "boundary_fraction": -7.252636656194852e-24,
      "field": "fz",
      "final": 0.0,
      "final_fraction": 0.0,
      "initial": 0.0,
      "initial_fraction": 0.0,
      "integral_units": "code units",
      "max_abs_residual": 4.1964200808574826e-22,
      "max_normalized_error": 1.3113812752679633e-23,
      "normalization_basis": "c * bath_energy * volume",
      "normalization_scale": 32.0,
      "normalized_error": 7.252636656194852e-24,
      "normalized_residual": -7.252636656194852e-24,
      "raw_change": 0.0,
      "residual": -2.320843729982353e-22,
      "source": 0.0,
      "source_fraction": 0.0,
      "zero_initial_flux_relative_tolerance": 1.4210854715202004e-14
    }
  ],
  "dx": 0.0625,
  "executable_sha256": "fc0e982bc397d82e25c7105cff6a97a514f09c57439e579587f4ddf6f6a8fa14",
  "inx": 8,
  "length": 2.0,
  "level": 2,
  "marshak_capacity_ratio": 1.0,
  "marshak_opacity": 100.0,
  "max_level": 2,
  "mesh": {
    "full_mesh": {
      "last_regrid_step": 195.0,
      "last_regrid_time": 1.0,
      "leaf_blocks": 36,
      "leaf_cells": 18432,
      "levels": {
        "1": {
          "leaf_blocks": 4,
          "leaf_cells": 2048
        },
        "2": {
          "leaf_blocks": 32,
          "leaf_cells": 16384
        }
      }
    },
    "geometry": "Positive-side leaf cells touching geometric z=0; cell centers have z=leaf_dx/2",
    "native_slice_cells": 640,
    "projection": "Piecewise-constant leaf cell averages onto finest display pixels; raw CSV retains native cells",
    "slice_cells_by_level": {
      "1": 128,
      "2": 512
    }
  },
  "mesh_description": "Adaptive Cartesian leaves; cells/dx denote configured finest equivalent resolution",
  "min_level": 1,
  "movie_capture": {
    "hard_dt": 0.01,
    "odt": 0.05,
    "requested_snapshots": 21,
    "steps_per_output_check": 5
  },
  "norms": {
    "er": {
      "L1": 0.05825905,
      "L2": 0.1963974,
      "Linf": 0.8907747
    },
    "fx": {
      "L1": 0.004134736,
      "L2": 0.01438703,
      "Linf": 0.07059959
    },
    "fy": {
      "L1": 0.0,
      "L2": 0.0,
      "Linf": 0.0
    },
    "fz": {
      "L1": 0.0,
      "L2": 0.0,
      "Linf": 0.0
    },
    "material_egas": {
      "L1": 0.05796755,
      "L2": 0.1956904,
      "Linf": 0.8869638
    }
  },
  "origin": "Octo-TIGER application",
  "physical_c": 1.0,
  "radiation_ader_version": 0,
  "radiation_backend": "cpu",
  "radiation_equations": "physical-c",
  "radiation_integrator": "euler",
  "radiation_reduced_c_factor": 1.0,
  "reference_kind": "Su-Olson nonequilibrium diffusion reference",
  "status": "complete",
  "study_config_sha256": "a7ffc569ea95d062e7f4d4b96786783a2c678d14b32981f7e48f1a25899973e6",
  "time": 1.0,
  "transport_c": 1.0,
  "units": {
    "c": 1.0,
    "energy_density": "normalized energy density",
    "flux": "normalized flux",
    "length": "normalized length",
    "system": "dimensionless",
    "time": "normalized time"
  },
  "validation_scope": "Normalized stationary material, equal radiation/material energy capacities. Analytic diffusion comparison is not an exact solution of nonlinear M1 transport."
}

Final fields and profiles

AMR: positive-side cells touching z=0, projected as constant leaf averages. Centers have z=leaf_dx/2; the ray similarly touches y=0. This is not a common cell-center plane. Raw CSV preserves native cells; display CSV repeats coarse values on finest pixels, retaining source dx/z. Full-volume norms use actual leaf volumes; uniform-grid convergence orders are not reported.

Numerical, reference, and signed-error energy slices. AMR: positive-side cells touching z=0, projected as constant leaf averages. Centers have z=leaf_dx/2; the ray similarly touches y=0. This is not a common cell-center plane. Raw CSV preserves native cells; display CSV repeats coarse values on finest pixels, retaining source dx/z. Full-volume norms use actual leaf volumes; uniform-grid convergence orders are not reported. · PDF

Numerical/reference line profiles and signed errors for all four radiation fields · PDF

Normalized radiation and material energies with their semi-analytic diffusion references · PDF

All three flux maps

fx: numerical, reference, and signed-error slice · PDF

fy: numerical, reference, and signed-error slice · PDF

fz: numerical, reference, and signed-error slice · PDF

Movies

No completed movie is available for this resolution.

Conservation

All values below are dimensionless fractions of the fixed normalization scale at t_ref = 0.

Field	Scale	Q(t_ref)/scale	Q(t)/scale	Boundary B/scale	Source S/scale	Residual R/scale	Maximum \|R\|/scale
er	E_bath V	0.000000e+00	3.344373e-02	-6.669430e-02	-3.325057e-02	-2.775558e-17	2.775558e-17
fx	c E_bath V	0.000000e+00	1.437286e-03	-1.325689e-01	-1.311316e-01	0.000000e+00	1.387779e-17
fy	c E_bath V	0.000000e+00	0.000000e+00	-6.425456e-24	0.000000e+00	-6.425456e-24	8.518024e-24
fz	c E_bath V	0.000000e+00	0.000000e+00	-7.252637e-24	0.000000e+00	-7.252637e-24	1.311381e-23

Integrated radiation budgets and their transport/source-corrected residuals · PDF

Data and reproducibility

Effective configuration Run metadata Solver log Paired slice values CSV Slice units Raw conservation history Budget summary CSV Conservation fractions CSV Fraction units Conservation units L1 norms L2 norms L∞ norms Final numerical Silo Reference Silo

AMR levels 2–3; 90112 active cells; finest equivalent 64³ / level 3

Recorded setup and run provenance

Origin	Octo-TIGER application
Domain	Cube of side 2.000000e+00 (code length); centered at the origin
Resolution	AMR levels 2–3; 90112 active cells; finest equivalent 64³; finest Δx = 3.125000e-02 (code length)
Slice geometry	AMR: positive-side cells touching z=0, projected as constant leaf averages. Centers have z=leaf_dx/2; the ray similarly touches y=0. This is not a common cell-center plane. Raw CSV preserves native cells; display CSV repeats coarse values on finest pixels, retaining source dx/z. Full-volume norms use actual leaf volumes; uniform-grid convergence orders are not reported.
Final comparison time	1.000000e+00 (code time)
Radiation integration	Legacy first-order radiation update
Radiation backend	cpu
Physical light speed	1.000000e+00 (code speed)
Reduced-speed factor α	1.000000e+00
Radiation transport speed	1.000000e+00 (code speed)
Periodic boundaries	off
Hydrodynamics	off
Gravity	off
Implicit radiation source	on
Uniform grid	off
CFL	0.4
Timestep cap	0.01 (code time)
Output interval	0.050000000000000003 (code time)

{
  "adaptive": true,
  "background": 0.0,
  "c": 1.0,
  "case": "marshak",
  "cells": 64,
  "command": [
    "/home/dmarce1/workspace/octotiger/radiation_results/results/local-euler-ader2-20261011-restart-01/.campaign/launch-solver",
    "/home/dmarce1/workspace/octotiger/radiation_results/results/local-euler-ader2-20261011-restart-01/.campaign/runtime/octotiger",
    "--config_file=/home/dmarce1/workspace/octotiger/radiation_results/results/local-euler-ader2-20261011-restart-01/runs/l3/marshak/euler/marshak/l3/run.ini",
    "--hpx:threads=12"
  ],
  "comparison_signature": {
    "atomic_mass": "1",
    "atomic_number": "1",
    "cfl": "0.4",
    "code_to_cm": "1",
    "code_to_g": "1",
    "code_to_s": "1",
    "correct_am_hydro": "off",
    "disable_analytic": "off",
    "disable_diagnostics": "on",
    "eos": "IDEAL",
    "gravity": "off",
    "hard_dt": "0.01",
    "hydro": "off",
    "hydro_device_kernel_type": "OFF",
    "hydro_host_kernel_type": "LEGACY",
    "monopole_device_kernel_type": "OFF",
    "multipole_device_kernel_type": "OFF",
    "n_species": "1",
    "omega": "0",
    "periodic": "off",
    "problem": "MARSHAK",
    "rad_implicit": "on",
    "radiation": "on",
    "radiation_integrator": "euler",
    "radiation_reduced_c_factor": "1",
    "stop_step": "1000000",
    "stop_time": "1",
    "study_config_sha256": "a7ffc569ea95d062e7f4d4b96786783a2c678d14b32981f7e48f1a25899973e6",
    "unigrid": "off",
    "xscale": "1"
  },
  "config": {
    "atomic_mass": "1",
    "atomic_number": "1",
    "cfl": "0.4",
    "code_to_cm": "1",
    "code_to_g": "1",
    "code_to_s": "1",
    "correct_am_hydro": "off",
    "datadir": "/home/dmarce1/workspace/octotiger/radiation_results/results/local-euler-ader2-20261011-restart-01/runs/l3/marshak/euler/marshak/l3/",
    "disable_analytic": "off",
    "disable_diagnostics": "on",
    "disable_output": "off",
    "eos": "IDEAL",
    "gradient_refine_buffer": "0",
    "gradient_refine_coarsen_factor": "0.80000000000000004",
    "gradient_refine_field": "radiation_energy",
    "gradient_refine_floor": "9.9999999999999995e-21",
    "gradient_refine_mode": "absolute",
    "gradient_refine_threshold": "4",
    "gravity": "off",
    "hard_dt": "0.01",
    "hydro": "off",
    "hydro_device_kernel_type": "OFF",
    "hydro_host_kernel_type": "LEGACY",
    "max_level": "3",
    "min_level": "2",
    "monopole_device_kernel_type": "OFF",
    "multipole_device_kernel_type": "OFF",
    "n_species": "1",
    "odt": "0.050000000000000003",
    "omega": "0",
    "periodic": "off",
    "problem": "MARSHAK",
    "rad_implicit": "on",
    "radiation": "on",
    "radiation_backend": "cpu",
    "radiation_integrator": "euler",
    "radiation_reduced_c_factor": "1",
    "stop_step": "1000000",
    "stop_time": "1",
    "unigrid": "off",
    "xscale": "1"
  },
  "conservation": [
    {
      "boundary": -1.7196938735909348,
      "boundary_fraction": -0.053740433549716714,
      "field": "er",
      "final": 0.8621385721670598,
      "final_fraction": 0.02694183038022062,
      "initial": 0.0,
      "initial_fraction": 0.0,
      "integral_units": "code units",
      "max_abs_residual": 5.551115123125783e-16,
      "max_normalized_error": 1.734723475976807e-17,
      "normalization_basis": "bath_energy * volume",
      "normalization_scale": 32.0,
      "normalized_error": 1.0408340855860843e-17,
      "normalized_residual": 1.0408340855860843e-17,
      "raw_change": 0.8621385721670598,
      "residual": 3.3306690738754696e-16,
      "source": -0.8575553014238754,
      "source_fraction": -0.026798603169496105,
      "zero_initial_flux_relative_tolerance": 1.4210854715202004e-14
    },
    {
      "boundary": -4.419372396959257,
      "boundary_fraction": -0.13810538740497677,
      "field": "fx",
      "final": 0.04778446277834607,
      "final_fraction": 0.0014932644618233146,
      "initial": 0.0,
      "initial_fraction": 0.0,
      "integral_units": "code units",
      "max_abs_residual": 8.881784197001252e-16,
      "max_normalized_error": 2.7755575615628914e-17,
      "normalization_basis": "c * bath_energy * volume",
      "normalization_scale": 32.0,
      "normalized_error": 0.0,
      "normalized_residual": 0.0,
      "raw_change": 0.04778446277834607,
      "residual": 0.0,
      "source": -4.37158793418091,
      "source_fraction": -0.13661212294315345,
      "zero_initial_flux_relative_tolerance": 1.4210854715202004e-14
    },
    {
      "boundary": 4.435108929788294e-22,
      "boundary_fraction": 1.385971540558842e-23,
      "field": "fy",
      "final": 0.0,
      "final_fraction": 0.0,
      "initial": 0.0,
      "initial_fraction": 0.0,
      "integral_units": "code units",
      "max_abs_residual": 5.688763182772694e-22,
      "max_normalized_error": 1.777738494616467e-23,
      "normalization_basis": "c * bath_energy * volume",
      "normalization_scale": 32.0,
      "normalized_error": 1.385971540558842e-23,
      "normalized_residual": 1.385971540558842e-23,
      "raw_change": 0.0,
      "residual": 4.435108929788294e-22,
      "source": 0.0,
      "source_fraction": 0.0,
      "zero_initial_flux_relative_tolerance": 1.4210854715202004e-14
    },
    {
      "boundary": 1.0539701850429638e-21,
      "boundary_fraction": 3.293656828259262e-23,
      "field": "fz",
      "final": 0.0,
      "final_fraction": 0.0,
      "initial": 0.0,
      "initial_fraction": 0.0,
      "integral_units": "code units",
      "max_abs_residual": 1.2574311848598425e-21,
      "max_normalized_error": 3.929472452687008e-23,
      "normalization_basis": "c * bath_energy * volume",
      "normalization_scale": 32.0,
      "normalized_error": 3.293656828259262e-23,
      "normalized_residual": 3.293656828259262e-23,
      "raw_change": 0.0,
      "residual": 1.0539701850429638e-21,
      "source": 0.0,
      "source_fraction": 0.0,
      "zero_initial_flux_relative_tolerance": 1.4210854715202004e-14
    }
  ],
  "dx": 0.03125,
  "executable_sha256": "fc0e982bc397d82e25c7105cff6a97a514f09c57439e579587f4ddf6f6a8fa14",
  "inx": 8,
  "length": 2.0,
  "level": 3,
  "marshak_capacity_ratio": 1.0,
  "marshak_opacity": 100.0,
  "max_level": 3,
  "mesh": {
    "full_mesh": {
      "last_regrid_step": 385.0,
      "last_regrid_time": 1.0,
      "leaf_blocks": 176,
      "leaf_cells": 90112,
      "levels": {
        "2": {
          "leaf_blocks": 48,
          "leaf_cells": 24576
        },
        "3": {
          "leaf_blocks": 128,
          "leaf_cells": 65536
        }
      }
    },
    "geometry": "Positive-side leaf cells touching geometric z=0; cell centers have z=leaf_dx/2",
    "native_slice_cells": 1792,
    "projection": "Piecewise-constant leaf cell averages onto finest display pixels; raw CSV retains native cells",
    "slice_cells_by_level": {
      "2": 768,
      "3": 1024
    }
  },
  "mesh_description": "Adaptive Cartesian leaves; cells/dx denote configured finest equivalent resolution",
  "min_level": 2,
  "movie_capture": {
    "hard_dt": 0.01,
    "odt": 0.05,
    "requested_snapshots": 21,
    "steps_per_output_check": 5
  },
  "norms": {
    "er": {
      "L1": 0.02876046,
      "L2": 0.1117458,
      "Linf": 0.6246538
    },
    "fx": {
      "L1": 0.002266185,
      "L2": 0.0104408,
      "Linf": 0.07504904
    },
    "fy": {
      "L1": 0.0,
      "L2": 0.0,
      "Linf": 0.0
    },
    "fz": {
      "L1": 0.0,
      "L2": 0.0,
      "Linf": 0.0
    },
    "material_egas": {
      "L1": 0.0286496,
      "L2": 0.1115211,
      "Linf": 0.6247201
    }
  },
  "origin": "Octo-TIGER application",
  "physical_c": 1.0,
  "radiation_ader_version": 0,
  "radiation_backend": "cpu",
  "radiation_equations": "physical-c",
  "radiation_integrator": "euler",
  "radiation_reduced_c_factor": 1.0,
  "reference_kind": "Su-Olson nonequilibrium diffusion reference",
  "status": "complete",
  "study_config_sha256": "a7ffc569ea95d062e7f4d4b96786783a2c678d14b32981f7e48f1a25899973e6",
  "time": 1.0,
  "transport_c": 1.0,
  "units": {
    "c": 1.0,
    "energy_density": "normalized energy density",
    "flux": "normalized flux",
    "length": "normalized length",
    "system": "dimensionless",
    "time": "normalized time"
  },
  "validation_scope": "Normalized stationary material, equal radiation/material energy capacities. Analytic diffusion comparison is not an exact solution of nonlinear M1 transport."
}

Final fields and profiles

All three flux maps

Movies

No completed movie is available for this resolution.

Conservation

All values below are dimensionless fractions of the fixed normalization scale at t_ref = 0.

Field	Scale	Q(t_ref)/scale	Q(t)/scale	Boundary B/scale	Source S/scale	Residual R/scale	Maximum \|R\|/scale
er	E_bath V	0.000000e+00	2.694183e-02	-5.374043e-02	-2.679860e-02	1.040834e-17	1.734723e-17
fx	c E_bath V	0.000000e+00	1.493264e-03	-1.381054e-01	-1.366121e-01	0.000000e+00	2.775558e-17
fy	c E_bath V	0.000000e+00	0.000000e+00	1.385972e-23	0.000000e+00	1.385972e-23	1.777738e-23
fz	c E_bath V	0.000000e+00	0.000000e+00	3.293657e-23	0.000000e+00	3.293657e-23	3.929472e-23

Data and reproducibility

This problem: all conservation summaries CSV JSON

Reading the results

Errors, budgets, and movies

Every norm uses the full three-dimensional domain, with cell-volume weighting. Here e = numerical − reference and V is the domain volume:

L1 = Σ |e| ΔV / V, L2 = √(Σ e² ΔV / V), L∞ = max |e|
p = ln(error_coarse / error_fine) / ln(Δx_coarse / Δx_fine)

Orders require two completed runs with matching physical time, configuration, units, executable, and origin. An exact zero error has no logarithmic order. The central slice and line profiles are visual diagnostics: the report slice is the cell layer at z = Δx/2, and its line profile is at y = z = Δx/2.

Conservation is assessed using the measured domain integrals and the signed cumulative boundary and source budgets:

A labeled checkpoint continuation uses its checkpoint as t_ref, with counters and fixed normalization starting there. Its reported residuals cover that resumed segment only; original histories remain separate.

R(t) = Q(t) − Q(0) + B(t) − S(t)

B is outward transport and S is the source change. Conservation plots and this page's budget tables show fractions of the initial totals: Q(t)/D, [Q(0) − B(t) + S(t)]/D, and the signed residual R(t)/D. The denominator D = |Q(0)| is fixed for the entire history. It never grows with the current total or accumulated transport/source budgets. A negative initial flux therefore starts at −1, and the plotted fraction retains its sign.

A component whose initial net flux is zero or negligible cannot be divided by its own initial total. It instead uses the fixed scale D = c |E₀|, where E₀ is the initial integrated radiation energy. This fallback is labeled on its panels and in the table. “Negligible” means |Q(0)|/(c |E₀|) ≤ 64ε ≈ 1.421085 × 10⁻¹⁴, with ε the double-precision machine epsilon. This prevents roundoff left by cancellation in a symmetric field from becoming the denominator. For the initially cold Marshak case, the prescribed bath sets D = E_bath V for energy and D = c E_bath V for flux, with E_bath = 4/c. These scales remain fixed even though radiation enters the initially empty domain.

Raw histories and dimensional summary columns remain in physical units: energy integrals in erg, physical-flux integrals in erg cm/s for CGS runs. Dividing a flux integral by c² gives radiation momentum. Fraction-history CSV files are dimensionless except for their time column; summary files record each denominator and its basis. Maximum fractional residuals use the same fixed denominator as the final residual.

For runs labeled RSLA, these budgets check the modified radiation equations. They do not establish conservation of physical gas-plus-radiation energy. With reduced-speed factor α, the closed-system coupled invariants are Egas + Er/α and gas momentum + F/(αc²), with physical c in the closure. Compare multiple α values before interpreting an RSLA result as an approximation to physical-c evolution.

Movies show the numerical Silo states with a fixed color range for each movie. Recorded states are held for their actual simulation-time gaps; intermediate solutions are not interpolated. The 3D option shows the domain's exterior surface. Movies and plots retain their recorded physical units.

Sources & provenance

Where this problem comes from

Su & Olson (1996), “Benchmark results for the non-equilibrium Marshak diffusion problem,” JQSRT 56, 337–351, gives the nonequilibrium diffusion benchmark. The local reference is implemented in src/test_problems/marshak/marshak.cpp; the run records source hashes for reproducibility.

Radiation method context

Skinner & Ostriker (2013), “A Two-moment Radiation Hydrodynamics Module in Athena Using a Time-explicit Godunov Method,” ApJS 206, 21, provides context for the two-moment M1 system and source integration. For the radiation-only Gaussian and sphere cases, Euler uses backward-Euler damping: Fⁿ⁺¹ = F^*/(1 + cχΔt), where F^* is the flux after transport. This uses the algebra of their equations (42b)–(43) with θ = 1 and prescribed χ. Their paper neglects scattering; the suite reuses that damping formula for its prescribed scattering medium. ADER2 revision 6 uses the same two-stage Radau IIA local integration for all eight problems. Transport is evaluated on the old state, a predicted time average and a corrected time average, with reconstruction recomputed each time. Prescribed-medium constraints are applied within the implicit stages: the Gaussian and sphere tests hold the material fixed; Marshak holds density and velocity fixed while thermal energy evolves; streaming disables interaction. The coupled problems evolve gas and radiation together on one shared timestep. This implementation does not yet include gravity, and uniform stiff coupled second-order accuracy remains unresolved. Historical revision 5 used exponential affine-forcing solutions for prescribed scattering, SDIRK2 for stationary Marshak heating, and Radau stages for coupled hydro–radiation. Those saved results retain their original method revision and do not validate revision 6. Here the formulas use physical light speed; reduced-speed runs retain their recorded transport speed.

The solver's M1 implementation also cites Hanawa & Audit (2014), “Reformulation of the M1 model of radiative transfer,” JQSRT 145, 9–16. These references describe the numerical/physical framework; the exact test definitions are the source functions identified above.

Reproducing this batch

Descriptions document the reviewed source implementation. Older or edited solver versions can use different profile constants or averaging rules; these are not all encoded in run.ini. Use the run metadata, executable hash, and saved batch source hashes to identify the implementation behind a result. A source path or descriptive formula alone is not proof of the executable's contents.

Batch settings and source hashes Batch status