Octo-TIGER radiation results

Octo-TIGER application

Equilibrium sphere · Gaussian pulse · Streaming front · Streaming wave · All norms/orders CSV · All conservation budgets CSV · Conservation JSON · Movie gallery

Norms use the entire 3D domain: L1 = Σ|e|ΔV/V, L2 = √(Σe²ΔV/V), L∞ = max|e|. Orders are p = log(error coarse / error fine) / log(Δx coarse / Δx fine), measured at the same final time. Zero norms have no reported order. CGS: positions and cell sizes in cm, time in s, E and its errors in erg/cm³, Fx,Fy,Fz and their errors in erg/(cm² s). CSV numbers use these same units. The central cell layer moves toward z=0 with refinement; each numerical value is compared to its own matched cell average.

The first-order time update can limit smooth-case convergence. A moving discontinuity can have different orders in different norms. These figures show measurements; they do not apply the CTest pass/fail thresholds.

Equilibrium sphere

Linearized / diffusion-limit reference; finite-amplitude nonlinear M1 has a model error floor.

Convergence PDF

Four-field convergence plot

Energy errors (erg/cm³) and observed orders

NL1L2L∞p(L1)p(L2)p(L∞)
322.018819e-055.148698e-056.422827e-04
644.416456e-061.156222e-053.466218e-042.192551e+002.154789e+008.898460e-01
32³ cells — level 2

Slice at z = 9.375000e+08 cm. Profiles at y = z = Δx/2. Paired values CSV · Profiles PDF · CSV units

Energy numerical, reference, and signed error Radiation profiles and signed errors
Flux slice maps

Fx PDF

Fx slice

Fy PDF

Fy slice

Fz PDF

Fz slice

Radiation conservation

Full-domain integrals Q = ∫(E, Fx, Fy, Fz) dV. The budget residual is R(t) = Q(t) − Q(0) + B(t) − S(t), where B is cumulative net outward boundary transport and S is the cumulative signed change actually applied by radiation source updates. Physical sources and boundary transport can change the radiation totals. Radiation momentum is ∫F dV / c²; division by c² gives the same relative budget errors.

Normalized errors are |R(t)| / scale(t). For E, scale = max(|E₀|, |E(t)|, |B_E(t)|, |S_E(t)|). For each flux component, scale also includes c times this energy scale and that component's |Q₀|, |Q(t)|, |B(t)| and |S(t)|. This remains defined for zero net flux. A zero scale with zero residual reports zero; no fixed physical-unit floor is used. Maximum errors cover all recorded timesteps.

IntegralInitialFinal Raw changeBoundary BSource SResidual R Final |R|/scaleMax |R|/scale
E (erg)2.190278e+322.190273e+32-4.391900e+261.079692e+301.079253e+30-2.322169e+161.060216e-161.124470e-16
Fx (erg cm/s)-9.407839e+212.307171e+242.316579e+24-4.603658e+25-4.364027e+25-7.973144e+221.214254e-201.231883e-20
Fy (erg cm/s)-5.829171e+213.042569e+243.048398e+24-7.144645e+25-6.823781e+25-1.602492e+232.440483e-204.352114e-20
Fz (erg cm/s)-3.689349e+19-3.314880e+23-3.314511e+238.044846e+247.530694e+241.827012e+232.782412e-202.848641e-20

Conservation history CSV · Budget summary CSV · Budget summary JSON · CSV units · History PNG · History PDF

Energy and three flux totals with source and boundary corrected conservation residuals
64³ cells — level 3

Slice at z = 4.687500e+08 cm. Profiles at y = z = Δx/2. Paired values CSV · Profiles PDF · CSV units

Energy numerical, reference, and signed error Radiation profiles and signed errors
Flux slice maps

Fx PDF

Fx slice

Fy PDF

Fy slice

Fz PDF

Fz slice

Radiation conservation

Full-domain integrals Q = ∫(E, Fx, Fy, Fz) dV. The budget residual is R(t) = Q(t) − Q(0) + B(t) − S(t), where B is cumulative net outward boundary transport and S is the cumulative signed change actually applied by radiation source updates. Physical sources and boundary transport can change the radiation totals. Radiation momentum is ∫F dV / c²; division by c² gives the same relative budget errors.

Normalized errors are |R(t)| / scale(t). For E, scale = max(|E₀|, |E(t)|, |B_E(t)|, |S_E(t)|). For each flux component, scale also includes c times this energy scale and that component's |Q₀|, |Q(t)|, |B(t)| and |S(t)|. This remains defined for zero net flux. A zero scale with zero residual reports zero; no fixed physical-unit floor is used. Maximum errors cover all recorded timesteps.

IntegralInitialFinal Raw changeBoundary BSource SResidual R Final |R|/scaleMax |R|/scale
E (erg)2.190278e+322.190277e+32-9.207343e+251.079345e+301.079253e+30-1.956251e+168.931520e-171.188727e-16
Fx (erg cm/s)-3.320414e+218.274103e+238.307307e+23-1.397909e+25-1.316113e+251.276745e+221.944394e-213.622956e-21
Fy (erg cm/s)-1.660207e+21-6.972869e+22-6.806849e+223.723291e+243.672198e+24-1.697562e+222.585268e-211.347029e-20
Fz (erg cm/s)-2.213609e+201.797082e+231.799295e+233.234526e+243.371555e+244.290021e+226.533404e-211.543611e-20

Conservation history CSV · Budget summary CSV · Budget summary JSON · CSV units · History PNG · History PDF

Energy and three flux totals with source and boundary corrected conservation residuals

Gaussian pulse

Linearized / diffusion-limit reference; finite-amplitude nonlinear M1 has a model error floor.

Convergence PDF

Four-field convergence plot

Energy errors (erg/cm³) and observed orders

NL1L2L∞p(L1)p(L2)p(L∞)
322.983866e-073.652407e-079.418877e-07
648.208483e-081.002988e-072.674563e-071.861995e+001.864543e+001.816252e+00
32³ cells — level 2

Slice at z = 9.375000e+08 cm. Profiles at y = z = Δx/2. Paired values CSV · Profiles PDF · CSV units

Energy numerical, reference, and signed error Radiation profiles and signed errors
Flux slice maps

Fx PDF

Fx slice

Fy PDF

Fy slice

Fz PDF

Fz slice

Radiation conservation

Full-domain integrals Q = ∫(E, Fx, Fy, Fz) dV. The budget residual is R(t) = Q(t) − Q(0) + B(t) − S(t), where B is cumulative net outward boundary transport and S is the cumulative signed change actually applied by radiation source updates. Physical sources and boundary transport can change the radiation totals. Radiation momentum is ∫F dV / c²; division by c² gives the same relative budget errors.

Normalized errors are |R(t)| / scale(t). For E, scale = max(|E₀|, |E(t)|, |B_E(t)|, |S_E(t)|). For each flux component, scale also includes c times this energy scale and that component's |Q₀|, |Q(t)|, |B(t)| and |S(t)|. This remains defined for zero net flux. A zero scale with zero residual reports zero; no fixed physical-unit floor is used. Maximum errors cover all recorded timesteps.

IntegralInitialFinal Raw changeBoundary BSource SResidual R Final |R|/scaleMax |R|/scale
E (erg)2.160120e+322.160120e+320.000000e+000.000000e+000.000000e+000.000000e+000.000000e+001.667907e-16
Fx (erg cm/s)0.000000e+00-1.527621e+19-1.527621e+190.000000e+00-6.088281e+205.935519e+209.165583e-231.865357e-22
Fy (erg cm/s)0.000000e+00-2.594073e+18-2.594073e+180.000000e+003.919032e+19-4.178440e+196.452314e-242.472506e-22
Fz (erg cm/s)0.000000e+002.154522e+192.154522e+190.000000e+00-3.461016e+195.615538e+198.671471e-241.594747e-22

Conservation history CSV · Budget summary CSV · Budget summary JSON · CSV units · History PNG · History PDF

Energy and three flux totals with source and boundary corrected conservation residuals
64³ cells — level 3

Slice at z = 4.687500e+08 cm. Profiles at y = z = Δx/2. Paired values CSV · Profiles PDF · CSV units

Energy numerical, reference, and signed error Radiation profiles and signed errors
Flux slice maps

Fx PDF

Fx slice

Fy PDF

Fy slice

Fz PDF

Fz slice

Radiation conservation

Full-domain integrals Q = ∫(E, Fx, Fy, Fz) dV. The budget residual is R(t) = Q(t) − Q(0) + B(t) − S(t), where B is cumulative net outward boundary transport and S is the cumulative signed change actually applied by radiation source updates. Physical sources and boundary transport can change the radiation totals. Radiation momentum is ∫F dV / c²; division by c² gives the same relative budget errors.

Normalized errors are |R(t)| / scale(t). For E, scale = max(|E₀|, |E(t)|, |B_E(t)|, |S_E(t)|). For each flux component, scale also includes c times this energy scale and that component's |Q₀|, |Q(t)|, |B(t)| and |S(t)|. This remains defined for zero net flux. A zero scale with zero residual reports zero; no fixed physical-unit floor is used. Maximum errors cover all recorded timesteps.

IntegralInitialFinal Raw changeBoundary BSource SResidual R Final |R|/scaleMax |R|/scale
E (erg)2.160120e+322.160120e+323.602880e+160.000000e+000.000000e+003.602880e+161.667907e-161.667907e-16
Fx (erg cm/s)0.000000e+001.297037e+181.297037e+180.000000e+00-1.005814e+211.007111e+211.555173e-221.611538e-22
Fy (erg cm/s)0.000000e+00-2.954361e+18-2.954361e+180.000000e+00-5.016852e+204.987309e+207.701363e-239.167078e-23
Fz (erg cm/s)0.000000e+00-7.205759e+17-7.205759e+170.000000e+00-1.674168e+201.666962e+202.574110e-234.639644e-23

Conservation history CSV · Budget summary CSV · Budget summary JSON · CSV units · History PNG · History PDF

Energy and three flux totals with source and boundary corrected conservation residuals

Streaming front

Exact nonlinear M1 streaming reference.

Convergence PDF

Four-field convergence plot

Energy errors (erg/cm³) and observed orders

NL1L2L∞p(L1)p(L2)p(L∞)
328.321504e-021.450355e-014.179089e-01
644.231950e-021.015654e-014.205968e-019.755217e-015.139971e-01-9.249390e-03
32³ cells — level 2

Slice at z = 9.375000e+08 cm. Profiles at y = z = Δx/2. Paired values CSV · Profiles PDF · CSV units

Energy numerical, reference, and signed error Radiation profiles and signed errors
Flux slice maps

Fx PDF

Fx slice

Fy PDF

Fy slice

Fz PDF

Fz slice

Radiation conservation

Full-domain integrals Q = ∫(E, Fx, Fy, Fz) dV. The budget residual is R(t) = Q(t) − Q(0) + B(t) − S(t), where B is cumulative net outward boundary transport and S is the cumulative signed change actually applied by radiation source updates. Physical sources and boundary transport can change the radiation totals. Radiation momentum is ∫F dV / c²; division by c² gives the same relative budget errors.

Normalized errors are |R(t)| / scale(t). For E, scale = max(|E₀|, |E(t)|, |B_E(t)|, |S_E(t)|). For each flux component, scale also includes c times this energy scale and that component's |Q₀|, |Q(t)|, |B(t)| and |S(t)|. This remains defined for zero net flux. A zero scale with zero residual reports zero; no fixed physical-unit floor is used. Maximum errors cover all recorded timesteps.

IntegralInitialFinal Raw changeBoundary BSource SResidual R Final |R|/scaleMax |R|/scale
E (erg)1.080000e+321.080000e+32-1.801440e+160.000000e+000.000000e+00-1.801440e+161.668000e-161.668000e-16
Fx (erg cm/s)3.237759e+423.237759e+42-6.189700e+260.000000e+000.000000e+00-6.189700e+261.911724e-163.823448e-16
Fy (erg cm/s)0.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+00
Fz (erg cm/s)0.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+00

Conservation history CSV · Budget summary CSV · Budget summary JSON · CSV units · History PNG · History PDF

Energy and three flux totals with source and boundary corrected conservation residuals
64³ cells — level 3

Slice at z = 4.687500e+08 cm. Profiles at y = z = Δx/2. Paired values CSV · Profiles PDF · CSV units

Energy numerical, reference, and signed error Radiation profiles and signed errors
Flux slice maps

Fx PDF

Fx slice

Fy PDF

Fy slice

Fz PDF

Fz slice

Radiation conservation

Full-domain integrals Q = ∫(E, Fx, Fy, Fz) dV. The budget residual is R(t) = Q(t) − Q(0) + B(t) − S(t), where B is cumulative net outward boundary transport and S is the cumulative signed change actually applied by radiation source updates. Physical sources and boundary transport can change the radiation totals. Radiation momentum is ∫F dV / c²; division by c² gives the same relative budget errors.

Normalized errors are |R(t)| / scale(t). For E, scale = max(|E₀|, |E(t)|, |B_E(t)|, |S_E(t)|). For each flux component, scale also includes c times this energy scale and that component's |Q₀|, |Q(t)|, |B(t)| and |S(t)|. This remains defined for zero net flux. A zero scale with zero residual reports zero; no fixed physical-unit floor is used. Maximum errors cover all recorded timesteps.

IntegralInitialFinal Raw changeBoundary BSource SResidual R Final |R|/scaleMax |R|/scale
E (erg)1.080000e+321.080000e+32-1.801440e+160.000000e+000.000000e+00-1.801440e+161.668000e-163.336000e-16
Fx (erg cm/s)3.237759e+423.237759e+42-6.189700e+260.000000e+000.000000e+00-6.189700e+261.911724e-163.823448e-16
Fy (erg cm/s)0.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+00
Fz (erg cm/s)0.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+00

Conservation history CSV · Budget summary CSV · Budget summary JSON · CSV units · History PNG · History PDF

Energy and three flux totals with source and boundary corrected conservation residuals

Streaming wave

Exact nonlinear M1 streaming reference.

Convergence PDF

Four-field convergence plot

Energy errors (erg/cm³) and observed orders

NL1L2L∞p(L1)p(L2)p(L∞)
322.385068e-012.684704e-013.902711e-01
646.439468e-027.202155e-021.155435e-011.889017e+001.898262e+001.756041e+00
32³ cells — level 2

Slice at z = 9.375000e+08 cm. Profiles at y = z = Δx/2. Paired values CSV · Profiles PDF · CSV units

Energy numerical, reference, and signed error Radiation profiles and signed errors
Flux slice maps

Fx PDF

Fx slice

Fy PDF

Fy slice

Fz PDF

Fz slice

Radiation conservation

Full-domain integrals Q = ∫(E, Fx, Fy, Fz) dV. The budget residual is R(t) = Q(t) − Q(0) + B(t) − S(t), where B is cumulative net outward boundary transport and S is the cumulative signed change actually applied by radiation source updates. Physical sources and boundary transport can change the radiation totals. Radiation momentum is ∫F dV / c²; division by c² gives the same relative budget errors.

Normalized errors are |R(t)| / scale(t). For E, scale = max(|E₀|, |E(t)|, |B_E(t)|, |S_E(t)|). For each flux component, scale also includes c times this energy scale and that component's |Q₀|, |Q(t)|, |B(t)| and |S(t)|. This remains defined for zero net flux. A zero scale with zero residual reports zero; no fixed physical-unit floor is used. Maximum errors cover all recorded timesteps.

IntegralInitialFinal Raw changeBoundary BSource SResidual R Final |R|/scaleMax |R|/scale
E (erg)2.160000e+322.160000e+323.602880e+160.000000e+000.000000e+003.602880e+161.668000e-161.668000e-16
Fx (erg cm/s)3.461309e+423.461309e+420.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+00
Fy (erg cm/s)-5.191964e+42-5.191964e+420.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+00
Fz (erg cm/s)1.730655e+421.730655e+420.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+00

Conservation history CSV · Budget summary CSV · Budget summary JSON · CSV units · History PNG · History PDF

Energy and three flux totals with source and boundary corrected conservation residuals
64³ cells — level 3

Slice at z = 4.687500e+08 cm. Profiles at y = z = Δx/2. Paired values CSV · Profiles PDF · CSV units

Energy numerical, reference, and signed error Radiation profiles and signed errors
Flux slice maps

Fx PDF

Fx slice

Fy PDF

Fy slice

Fz PDF

Fz slice

Radiation conservation

Full-domain integrals Q = ∫(E, Fx, Fy, Fz) dV. The budget residual is R(t) = Q(t) − Q(0) + B(t) − S(t), where B is cumulative net outward boundary transport and S is the cumulative signed change actually applied by radiation source updates. Physical sources and boundary transport can change the radiation totals. Radiation momentum is ∫F dV / c²; division by c² gives the same relative budget errors.

Normalized errors are |R(t)| / scale(t). For E, scale = max(|E₀|, |E(t)|, |B_E(t)|, |S_E(t)|). For each flux component, scale also includes c times this energy scale and that component's |Q₀|, |Q(t)|, |B(t)| and |S(t)|. This remains defined for zero net flux. A zero scale with zero residual reports zero; no fixed physical-unit floor is used. Maximum errors cover all recorded timesteps.

IntegralInitialFinal Raw changeBoundary BSource SResidual R Final |R|/scaleMax |R|/scale
E (erg)2.160000e+322.160000e+320.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+00
Fx (erg cm/s)3.461309e+423.461309e+420.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+00
Fy (erg cm/s)-5.191964e+42-5.191964e+420.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+00
Fz (erg cm/s)1.730655e+421.730655e+420.000000e+000.000000e+000.000000e+000.000000e+000.000000e+000.000000e+00

Conservation history CSV · Budget summary CSV · Budget summary JSON · CSV units · History PNG · History PDF

Energy and three flux totals with source and boundary corrected conservation residuals