TY - JOUR

T1 - Numerical evaluation of the three-dimensional searchlight problem in a half-space

AU - Kornreich, D. E.

AU - Ganapol, B. D.

PY - 1997/11

Y1 - 1997/11

N2 - The linear Boltzmann equation for the transport of neutral particles is investigated with the objective of generating a benchmark-quality calculation for the three-dimensional searchlight problem in a semi-infinite medium. The derivation assumes stationarity, one energy group, and Isotropic scattering. The scalar flux (both surface and interior) and the current at the surface are the quantities of interest. The source considered is a pencil-beam incident at a point on the surface of a semi-infinite medium. The scalar flux will have two-dimensional variation only if the beam is normal; otherwise, it is three-dimensional. The solutions are obtained by using Fourier and Laplace transform methods. The transformed transport equation is formulated so that it can be related to a one-dimensional pseudo problem, thus providing some analytical leverage for the inversions. The numerical inversions use standard numerical techniques such as Gauss-Legendre quadrature, summation of infinite series, H-function iteration and evaluation, and Euler-Knopp acceleration. The numerical evaluations of the scalar flux and current at the surface are relatively simple, and the interior scalar flux is relatively difficult to calculate because of the embedded two-dimensional Fourier transform inversion, Laplace transform inversion, and H-function evaluation. Comparisons of these numerical solutions to results from the MCNP probabilistic code and the THREE-DANT discrete ordinates code are provided and help confirm proper operation of the analytical code.

AB - The linear Boltzmann equation for the transport of neutral particles is investigated with the objective of generating a benchmark-quality calculation for the three-dimensional searchlight problem in a semi-infinite medium. The derivation assumes stationarity, one energy group, and Isotropic scattering. The scalar flux (both surface and interior) and the current at the surface are the quantities of interest. The source considered is a pencil-beam incident at a point on the surface of a semi-infinite medium. The scalar flux will have two-dimensional variation only if the beam is normal; otherwise, it is three-dimensional. The solutions are obtained by using Fourier and Laplace transform methods. The transformed transport equation is formulated so that it can be related to a one-dimensional pseudo problem, thus providing some analytical leverage for the inversions. The numerical inversions use standard numerical techniques such as Gauss-Legendre quadrature, summation of infinite series, H-function iteration and evaluation, and Euler-Knopp acceleration. The numerical evaluations of the scalar flux and current at the surface are relatively simple, and the interior scalar flux is relatively difficult to calculate because of the embedded two-dimensional Fourier transform inversion, Laplace transform inversion, and H-function evaluation. Comparisons of these numerical solutions to results from the MCNP probabilistic code and the THREE-DANT discrete ordinates code are provided and help confirm proper operation of the analytical code.

UR - http://www.scopus.com/inward/record.url?scp=0001466608&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0001466608&partnerID=8YFLogxK

U2 - 10.13182/NSE97-A1938

DO - 10.13182/NSE97-A1938

M3 - Review article

AN - SCOPUS:0001466608

VL - 127

SP - 317

EP - 337

JO - Nuclear Science and Engineering

JF - Nuclear Science and Engineering

SN - 0029-5639

IS - 3

ER -