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

D. E. Kornreich, Barry D Ganapol

Research output: Contribution to journalArticle

9 Citations (Scopus)

Abstract

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.

Original languageEnglish (US)
Pages (from-to)317-337
Number of pages21
JournalNuclear Science and Engineering
Volume127
Issue number3
StatePublished - Nov 1997

Fingerprint

Searchlights
Fluxes
Laplace transforms
Fourier transforms
Boltzmann equation
Function evaluation
Scattering

ASJC Scopus subject areas

  • Nuclear Energy and Engineering

Cite this

Numerical evaluation of the three-dimensional searchlight problem in a half-space. / Kornreich, D. E.; Ganapol, Barry D.

In: Nuclear Science and Engineering, Vol. 127, No. 3, 11.1997, p. 317-337.

Research output: Contribution to journalArticle

@article{86b69be542b445239b2bfef7f3146d50,
title = "Numerical evaluation of the three-dimensional searchlight problem in a half-space",
abstract = "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.",
author = "Kornreich, {D. E.} and Ganapol, {Barry D}",
year = "1997",
month = "11",
language = "English (US)",
volume = "127",
pages = "317--337",
journal = "Nuclear Science and Engineering",
issn = "0029-5639",
publisher = "American Nuclear Society",
number = "3",

}

TY - JOUR

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

AU - Kornreich, D. E.

AU - Ganapol, Barry 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

M3 - 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 -