All problem comparisons · Campaign

Marshak: Euler and ADER2

Reference: Su–Olson diffusion reference. The supplied Su–Olson radiation/material reference solves the linear thermal-capacity diffusion problem in normalized units. The saved parameter values define its match to the source model; diffusion and finite-speed M1 transport can differ. Reference values are evaluated at cell centers, without an additional cell-volume average.

Su and Olson (1996), “Benchmark results for the non-equilibrium Marshak diffusion problem,” JQSRT 56, 337–351.

Coupled ADER2 accuracy remains under investigation. A completed run or a small difference does not establish uniform second-order accuracy.

Revision 6 validation limits. Streaming-wave temporal-order thresholds failed at limiter transitions even with resolved references. Optically thick radiation temporal probes remain unresolved under the unchanged 1e-12 measurement guard. Strict stiff coupled temporal tests have not been rerun for revision 6 and remain an open limitation.

Level 2

Euler (first order): complete results, parameters and reference discussion

Coupled ADER2: complete results, parameters and reference discussion

matched. No other refinement level or physical time is substituted for a missing partner.

Level 2 method and available reference profile overlays

Matched profile coordinates and values (CSV)

ADER2 minus Euler: slice discrepancy

FieldUnitsL₁L₂L∞
Radiation energy densitynormalized energy density8.318004e-023.046462e-011.447607e+00
Radiation flux xnormalized flux2.273613e-038.824560e-034.698007e-02
Radiation flux ynormalized flux0.000000e+000.000000e+000.000000e+00
Radiation flux znormalized flux0.000000e+000.000000e+000.000000e+00
Material energy densitynormalized energy density8.267613e-023.033556e-011.443975e+00

Euler (first order) versus Su–Olson diffusion reference

Area-weighted errors on this method’s original leaf slice.

FieldUnitsL₁L₂L∞
Radiation energy densitynormalized energy density5.825905e-021.963974e-018.907747e-01
Radiation flux xnormalized flux4.134736e-031.438703e-027.059959e-02
Radiation flux ynormalized flux0.000000e+000.000000e+000.000000e+00
Radiation flux znormalized flux0.000000e+000.000000e+000.000000e+00
Material energy densitynormalized energy density5.796755e-021.956904e-018.869638e-01
Solver-reported full-volume reference norms
FieldUnitsL₁L₂L∞
Radiation energy densitynormalized energy density5.825905e-021.963974e-018.907747e-01
Radiation flux xnormalized flux4.134736e-031.438703e-027.059959e-02
Radiation flux ynormalized flux0.000000e+000.000000e+000.000000e+00
Radiation flux znormalized flux0.000000e+000.000000e+000.000000e+00
Material energy densitynormalized energy density5.796755e-021.956904e-018.869638e-01

Source slice CSV · Run metadata

Coupled ADER2 versus Su–Olson diffusion reference

Area-weighted errors on this method’s original leaf slice.

FieldUnitsL₁L₂L∞
Radiation energy densitynormalized energy density2.939480e-021.522552e-018.594192e-01
Radiation flux xnormalized flux2.941168e-031.359233e-027.580205e-02
Radiation flux ynormalized flux0.000000e+000.000000e+000.000000e+00
Radiation flux znormalized flux0.000000e+000.000000e+000.000000e+00
Material energy densitynormalized energy density2.915356e-021.515333e-018.553697e-01
Solver-reported full-volume reference norms
FieldUnitsL₁L₂L∞
Radiation energy densitynormalized energy density2.939480e-021.522552e-018.594192e-01
Radiation flux xnormalized flux2.941168e-031.359233e-027.580205e-02
Radiation flux ynormalized flux0.000000e+000.000000e+000.000000e+00
Radiation flux znormalized flux0.000000e+000.000000e+000.000000e+00
Material energy densitynormalized energy density2.915356e-021.515333e-018.553697e-01

Source slice CSV · Run metadata

Quantitative comparison and input hashes (JSON)

Level 3

Euler (first order): complete results, parameters and reference discussion

Coupled ADER2: complete results, parameters and reference discussion

matched. No other refinement level or physical time is substituted for a missing partner.

Level 3 method and available reference profile overlays

Matched profile coordinates and values (CSV)

ADER2 minus Euler: slice discrepancy

FieldUnitsL₁L₂L∞
Radiation energy densitynormalized energy density4.822940e-022.034221e-011.104716e+00
Radiation flux xnormalized flux2.263366e-038.809519e-034.676294e-02
Radiation flux ynormalized flux0.000000e+000.000000e+000.000000e+00
Radiation flux znormalized flux0.000000e+000.000000e+000.000000e+00
Material energy densitynormalized energy density4.801267e-022.029451e-011.105102e+00

Euler (first order) versus Su–Olson diffusion reference

Area-weighted errors on this method’s original leaf slice.

FieldUnitsL₁L₂L∞
Radiation energy densitynormalized energy density2.876046e-021.117458e-016.246538e-01
Radiation flux xnormalized flux2.266185e-031.044080e-027.504904e-02
Radiation flux ynormalized flux0.000000e+000.000000e+000.000000e+00
Radiation flux znormalized flux0.000000e+000.000000e+000.000000e+00
Material energy densitynormalized energy density2.864960e-021.115211e-016.247201e-01
Solver-reported full-volume reference norms
FieldUnitsL₁L₂L∞
Radiation energy densitynormalized energy density2.876046e-021.117458e-016.246538e-01
Radiation flux xnormalized flux2.266185e-031.044080e-027.504904e-02
Radiation flux ynormalized flux0.000000e+000.000000e+000.000000e+00
Radiation flux znormalized flux0.000000e+000.000000e+000.000000e+00
Material energy densitynormalized energy density2.864960e-021.115211e-016.247201e-01

Source slice CSV · Run metadata

Coupled ADER2 versus Su–Olson diffusion reference

Area-weighted errors on this method’s original leaf slice.

FieldUnitsL₁L₂L∞
Radiation energy densitynormalized energy density2.065564e-021.108698e-017.965644e-01
Radiation flux xnormalized flux1.788943e-038.431122e-035.432214e-02
Radiation flux ynormalized flux0.000000e+000.000000e+000.000000e+00
Radiation flux znormalized flux0.000000e+000.000000e+000.000000e+00
Material energy densitynormalized energy density2.054846e-021.106992e-017.967297e-01
Solver-reported full-volume reference norms
FieldUnitsL₁L₂L∞
Radiation energy densitynormalized energy density2.065564e-021.108698e-017.965644e-01
Radiation flux xnormalized flux1.788943e-038.431122e-035.432214e-02
Radiation flux ynormalized flux0.000000e+000.000000e+000.000000e+00
Radiation flux znormalized flux0.000000e+000.000000e+000.000000e+00
Material energy densitynormalized energy density2.054846e-021.106992e-017.967297e-01

Source slice CSV · Run metadata

Quantitative comparison and input hashes (JSON)

Level 4

Euler (first order): pending.

Coupled ADER2: pending.

Quantitative comparison and input hashes (JSON)

How these comparisons are measured

All field values use the saved units (CGS for the seven original problems; normalized units for Marshak). The central profiles share x coordinates and the positive row adjacent to y = z = 0. Adaptive leaf-cell averages are prolonged as constants onto the finest common Cartesian x-y partition, with no extrapolation or time interpolation. Different leaf sizes may average different z extents adjacent to the plane. The plotted curves therefore represent the saved finite-volume slices.

Slice L₁ is mean absolute difference, L₂ is root-mean-square difference, and L∞ is maximum absolute difference, with area weighting. These are two-dimensional slice diagnostics. The separately labeled solver norms cover the full volume.

Marshak numerical values are leaf-cell averages, while the Su–Olson reference is evaluated at those cell centers without a cell-volume average. The reported reference errors include this sampling difference as well as transport/model and discretization effects.