Benchmarks for Infinite Medium, Time Dependent Transport Problems with Isotropic Scattering

Abstract

The widely used AZURV1 transport benchmarks package provides a suite of solutions to isotropic scattering transport problems with a variety of initial conditions. Most of these solutions have an initial condition that is a Dirac delta function in space; as a result these benchmarks are challenging problems to use for verification tests in computer codes. Nevertheless, approximating a delta function in simulation often leads to low orders of convergence and the inability to test the convergence of high-order numerical methods. While there are examples in the literature of integration of these solutions as Green’s functions for the transport operator to produce results for more easily simulated sources, they are limited in scope and briefly explained. For a sampling of initial conditions and sources, we present solutions for the uncollided and collided scalar flux to facilitate accurate testing of source treatment in numerical solvers. The solution for the uncollided scalar flux is found in analytic form for some sources. Since integrating the Green’s functions is often nontrivial, discussion of integration difficulty and workarounds to find convergent integrals is included. Additionally, our uncollided solutions can be used as source terms in verification studies, in a similar way to the method of manufactured solutions.

Keywords:

• Benchmarks
• transport theory
• neutron transport
• photon transport
• transport equation

Acknowledgment

The authors would like to thank Minwoo Shin for providing independent numerical solutions to the cylindrical Gaussian pulse (Sec. 5).

Notes

1 There is contention whether one should use “Green ’ s function ” or “ Green function .” We follow the convention recommended by Wright (2006 ) in retaining the possessive.

2 github.com/wbennett39/transport_benchmarks.

Additional information

Funding

This work was supported by a Los Alamos National Laboratory under contract #599008, “Verification for Radiation Hydrodynamics Simulations.”