Critical discretization issues in 3-D SN simulations relevant to dosimetry and medical physics

Ahmad Khaled Al-Basheer, Glenn E. Sjoden, Monica Ghtta

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

Dosimetry problems inherently involve dose determinations among widely varying materials and densities, and may require complex, detailed investigations of the angular, spatial, and energy behavior of the applied radiation transporting throughout the simulation geometry. Traditionally, Monte Carlo codes have been implemented in solving these types of problems using voxelized geometries and phantoms. The motivation of this work is to investigate the discretization requirements for deterministic radiation transport simulations for these problems via direct solutions of the linear Boltzmann transport equation, focusing on the discrete ordinales (SN) method. The SN method can yield accurate global solutions, provided the inherent discretizations among the angular, spatial, and energy domains properly represent problem physics. In this paper, the Ss approach is implemented using a three-dimensional (3-D) 60Co Photon transport simulation to highlight the critical issues encountered in performing deterministic photon simulations in dosimetry problems. Calculations were performed using the PENTRAN parallel SN code to obtain a 3-D distribution of flux and dose computed using a collisional kerma approximation. For an acceptable result, we determined that a minimum angular Legendre-Chebychev quadrature of S 32 with P3 anisotropy is required, with block-adaptive meshes on the order of 1 cm, even in air regions, implemented with an adaptive differencing scheme (implemented in the PENTRAN code) to yield optimal solution convergence. Also, photon cross-section libraries should be carefully evaluated for the problem studied; for our test problem, the BUGLE-96 photon library yielded the closest results to Monte Carlo (MCNP5) among those tested. Overall, this work details the levels of discretization involved in performing deterministic computations in dosimetry problems and will be useful in enabling future efforts to perform rapid deterministic computations of phantom doses.

Original languageEnglish (US)
Pages (from-to)252-262
Number of pages11
JournalNuclear Technology
Volume169
Issue number3
DOIs
StatePublished - Jan 1 2010

Fingerprint

Dosimetry
dosimeters
Photons
Physics
physics
photons
dosage
simulation
Radiation
Boltzmann transport equation
radiation transport
Geometry
geometry
quadratures
mesh
Anisotropy
Fluxes
requirements
anisotropy
energy

Keywords

  • Accuracy
  • Discrete ordinates
  • Photons

ASJC Scopus subject areas

  • Nuclear and High Energy Physics
  • Nuclear Energy and Engineering
  • Condensed Matter Physics

Cite this

Critical discretization issues in 3-D SN simulations relevant to dosimetry and medical physics. / Al-Basheer, Ahmad Khaled; Sjoden, Glenn E.; Ghtta, Monica.

In: Nuclear Technology, Vol. 169, No. 3, 01.01.2010, p. 252-262.

Research output: Contribution to journalArticle

@article{0a9723d5b3f14783a3037931cd2e39c9,
title = "Critical discretization issues in 3-D SN simulations relevant to dosimetry and medical physics",
abstract = "Dosimetry problems inherently involve dose determinations among widely varying materials and densities, and may require complex, detailed investigations of the angular, spatial, and energy behavior of the applied radiation transporting throughout the simulation geometry. Traditionally, Monte Carlo codes have been implemented in solving these types of problems using voxelized geometries and phantoms. The motivation of this work is to investigate the discretization requirements for deterministic radiation transport simulations for these problems via direct solutions of the linear Boltzmann transport equation, focusing on the discrete ordinales (SN) method. The SN method can yield accurate global solutions, provided the inherent discretizations among the angular, spatial, and energy domains properly represent problem physics. In this paper, the Ss approach is implemented using a three-dimensional (3-D) 60Co Photon transport simulation to highlight the critical issues encountered in performing deterministic photon simulations in dosimetry problems. Calculations were performed using the PENTRAN parallel SN code to obtain a 3-D distribution of flux and dose computed using a collisional kerma approximation. For an acceptable result, we determined that a minimum angular Legendre-Chebychev quadrature of S 32 with P3 anisotropy is required, with block-adaptive meshes on the order of 1 cm, even in air regions, implemented with an adaptive differencing scheme (implemented in the PENTRAN code) to yield optimal solution convergence. Also, photon cross-section libraries should be carefully evaluated for the problem studied; for our test problem, the BUGLE-96 photon library yielded the closest results to Monte Carlo (MCNP5) among those tested. Overall, this work details the levels of discretization involved in performing deterministic computations in dosimetry problems and will be useful in enabling future efforts to perform rapid deterministic computations of phantom doses.",
keywords = "Accuracy, Discrete ordinates, Photons",
author = "Al-Basheer, {Ahmad Khaled} and Sjoden, {Glenn E.} and Monica Ghtta",
year = "2010",
month = "1",
day = "1",
doi = "10.13182/NT10-A9377",
language = "English (US)",
volume = "169",
pages = "252--262",
journal = "Nuclear Technology",
issn = "0029-5450",
publisher = "American Nuclear Society",
number = "3",

}

TY - JOUR

T1 - Critical discretization issues in 3-D SN simulations relevant to dosimetry and medical physics

AU - Al-Basheer, Ahmad Khaled

AU - Sjoden, Glenn E.

AU - Ghtta, Monica

PY - 2010/1/1

Y1 - 2010/1/1

N2 - Dosimetry problems inherently involve dose determinations among widely varying materials and densities, and may require complex, detailed investigations of the angular, spatial, and energy behavior of the applied radiation transporting throughout the simulation geometry. Traditionally, Monte Carlo codes have been implemented in solving these types of problems using voxelized geometries and phantoms. The motivation of this work is to investigate the discretization requirements for deterministic radiation transport simulations for these problems via direct solutions of the linear Boltzmann transport equation, focusing on the discrete ordinales (SN) method. The SN method can yield accurate global solutions, provided the inherent discretizations among the angular, spatial, and energy domains properly represent problem physics. In this paper, the Ss approach is implemented using a three-dimensional (3-D) 60Co Photon transport simulation to highlight the critical issues encountered in performing deterministic photon simulations in dosimetry problems. Calculations were performed using the PENTRAN parallel SN code to obtain a 3-D distribution of flux and dose computed using a collisional kerma approximation. For an acceptable result, we determined that a minimum angular Legendre-Chebychev quadrature of S 32 with P3 anisotropy is required, with block-adaptive meshes on the order of 1 cm, even in air regions, implemented with an adaptive differencing scheme (implemented in the PENTRAN code) to yield optimal solution convergence. Also, photon cross-section libraries should be carefully evaluated for the problem studied; for our test problem, the BUGLE-96 photon library yielded the closest results to Monte Carlo (MCNP5) among those tested. Overall, this work details the levels of discretization involved in performing deterministic computations in dosimetry problems and will be useful in enabling future efforts to perform rapid deterministic computations of phantom doses.

AB - Dosimetry problems inherently involve dose determinations among widely varying materials and densities, and may require complex, detailed investigations of the angular, spatial, and energy behavior of the applied radiation transporting throughout the simulation geometry. Traditionally, Monte Carlo codes have been implemented in solving these types of problems using voxelized geometries and phantoms. The motivation of this work is to investigate the discretization requirements for deterministic radiation transport simulations for these problems via direct solutions of the linear Boltzmann transport equation, focusing on the discrete ordinales (SN) method. The SN method can yield accurate global solutions, provided the inherent discretizations among the angular, spatial, and energy domains properly represent problem physics. In this paper, the Ss approach is implemented using a three-dimensional (3-D) 60Co Photon transport simulation to highlight the critical issues encountered in performing deterministic photon simulations in dosimetry problems. Calculations were performed using the PENTRAN parallel SN code to obtain a 3-D distribution of flux and dose computed using a collisional kerma approximation. For an acceptable result, we determined that a minimum angular Legendre-Chebychev quadrature of S 32 with P3 anisotropy is required, with block-adaptive meshes on the order of 1 cm, even in air regions, implemented with an adaptive differencing scheme (implemented in the PENTRAN code) to yield optimal solution convergence. Also, photon cross-section libraries should be carefully evaluated for the problem studied; for our test problem, the BUGLE-96 photon library yielded the closest results to Monte Carlo (MCNP5) among those tested. Overall, this work details the levels of discretization involved in performing deterministic computations in dosimetry problems and will be useful in enabling future efforts to perform rapid deterministic computations of phantom doses.

KW - Accuracy

KW - Discrete ordinates

KW - Photons

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

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

U2 - 10.13182/NT10-A9377

DO - 10.13182/NT10-A9377

M3 - Article

VL - 169

SP - 252

EP - 262

JO - Nuclear Technology

JF - Nuclear Technology

SN - 0029-5450

IS - 3

ER -