Deterministic radiation transport simulations for diagnostic imaging applications

Monica Ghita, Glenn E. Sjoden, Ahmad Khaled Al-Basheer, Manuel M. Arreola, Wesley Bolch, Choonsik Lee

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Abstract

Monte Carlo methods have been established along the time as the gold standard for computer simulations in the medical physics community. Depending on the problem and user's needs, deterministic radiation transport simulations may provide a more detailed and faster solution. In this work we investigate the possibility of using deterministic radiation transport simulations as a viable and more convenient tool for real clinical applications. Therefore, the discrete ordinates PENTRAN code is used to calculate average organ doses in voxelized human phantoms and the results are compared with state-of-the-art MCNP5 Monte Carlo simulations in the diagnostic energy range (50-140 keV). Generally, good agreement for the average organ scalar fluxes, less than 6% difference, is obtained provided adequate quadrature order, mesh size and energy group structure is used in the deterministic calculations. The energy group structure, particularly for the diagnostic energy range, has a major impact on the deterministic solution for the average organ doses since the interaction and mass energy absorption coefficients are highly energy dependent in the diagnostic range. Though an optimization of the group structure is possible, it is problem (namely x-ray source spectrum) and organ dependent, which impose serious limitations of the deterministic solution for practical application in diagnostic medical physics.

Original languageEnglish (US)
Title of host publicationAmerican Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009
Pages1822-1834
Number of pages13
StatePublished - Dec 1 2009
EventInternational Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009 - Saratoga Springs, NY, United States
Duration: May 3 2009May 7 2009

Publication series

NameAmerican Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009
Volume3

Other

OtherInternational Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009
CountryUnited States
CitySaratoga Springs, NY
Period5/3/095/7/09

Fingerprint

radiation transport
Diagnostics
organs
Radiation
Imaging
Imaging techniques
Physics
Energy
Simulation
simulation
Energy absorption
Dose
energy
Range of data
Discrete Ordinates
Monte Carlo methods
Energy Absorption
dosage
physics
Dependent

Keywords

  • Average organ dose
  • Deterministic discrete ordinates
  • Monte Carlo
  • Radiation transport
  • Voxelized phantom

ASJC Scopus subject areas

  • Nuclear Energy and Engineering
  • Computational Mathematics
  • Nuclear and High Energy Physics

Cite this

Ghita, M., Sjoden, G. E., Al-Basheer, A. K., Arreola, M. M., Bolch, W., & Lee, C. (2009). Deterministic radiation transport simulations for diagnostic imaging applications. In American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009 (pp. 1822-1834). (American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009; Vol. 3).

Deterministic radiation transport simulations for diagnostic imaging applications. / Ghita, Monica; Sjoden, Glenn E.; Al-Basheer, Ahmad Khaled; Arreola, Manuel M.; Bolch, Wesley; Lee, Choonsik.

American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009. 2009. p. 1822-1834 (American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009; Vol. 3).

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Ghita, M, Sjoden, GE, Al-Basheer, AK, Arreola, MM, Bolch, W & Lee, C 2009, Deterministic radiation transport simulations for diagnostic imaging applications. in American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009. American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009, vol. 3, pp. 1822-1834, International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009, Saratoga Springs, NY, United States, 5/3/09.
Ghita M, Sjoden GE, Al-Basheer AK, Arreola MM, Bolch W, Lee C. Deterministic radiation transport simulations for diagnostic imaging applications. In American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009. 2009. p. 1822-1834. (American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009).
Ghita, Monica ; Sjoden, Glenn E. ; Al-Basheer, Ahmad Khaled ; Arreola, Manuel M. ; Bolch, Wesley ; Lee, Choonsik. / Deterministic radiation transport simulations for diagnostic imaging applications. American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009. 2009. pp. 1822-1834 (American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009).
@inproceedings{72041591e37e4d0688311a7dacb7603e,
title = "Deterministic radiation transport simulations for diagnostic imaging applications",
abstract = "Monte Carlo methods have been established along the time as the gold standard for computer simulations in the medical physics community. Depending on the problem and user's needs, deterministic radiation transport simulations may provide a more detailed and faster solution. In this work we investigate the possibility of using deterministic radiation transport simulations as a viable and more convenient tool for real clinical applications. Therefore, the discrete ordinates PENTRAN code is used to calculate average organ doses in voxelized human phantoms and the results are compared with state-of-the-art MCNP5 Monte Carlo simulations in the diagnostic energy range (50-140 keV). Generally, good agreement for the average organ scalar fluxes, less than 6{\%} difference, is obtained provided adequate quadrature order, mesh size and energy group structure is used in the deterministic calculations. The energy group structure, particularly for the diagnostic energy range, has a major impact on the deterministic solution for the average organ doses since the interaction and mass energy absorption coefficients are highly energy dependent in the diagnostic range. Though an optimization of the group structure is possible, it is problem (namely x-ray source spectrum) and organ dependent, which impose serious limitations of the deterministic solution for practical application in diagnostic medical physics.",
keywords = "Average organ dose, Deterministic discrete ordinates, Monte Carlo, Radiation transport, Voxelized phantom",
author = "Monica Ghita and Sjoden, {Glenn E.} and Al-Basheer, {Ahmad Khaled} and Arreola, {Manuel M.} and Wesley Bolch and Choonsik Lee",
year = "2009",
month = "12",
day = "1",
language = "English (US)",
isbn = "9781615673490",
series = "American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009",
pages = "1822--1834",
booktitle = "American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009",

}

TY - GEN

T1 - Deterministic radiation transport simulations for diagnostic imaging applications

AU - Ghita, Monica

AU - Sjoden, Glenn E.

AU - Al-Basheer, Ahmad Khaled

AU - Arreola, Manuel M.

AU - Bolch, Wesley

AU - Lee, Choonsik

PY - 2009/12/1

Y1 - 2009/12/1

N2 - Monte Carlo methods have been established along the time as the gold standard for computer simulations in the medical physics community. Depending on the problem and user's needs, deterministic radiation transport simulations may provide a more detailed and faster solution. In this work we investigate the possibility of using deterministic radiation transport simulations as a viable and more convenient tool for real clinical applications. Therefore, the discrete ordinates PENTRAN code is used to calculate average organ doses in voxelized human phantoms and the results are compared with state-of-the-art MCNP5 Monte Carlo simulations in the diagnostic energy range (50-140 keV). Generally, good agreement for the average organ scalar fluxes, less than 6% difference, is obtained provided adequate quadrature order, mesh size and energy group structure is used in the deterministic calculations. The energy group structure, particularly for the diagnostic energy range, has a major impact on the deterministic solution for the average organ doses since the interaction and mass energy absorption coefficients are highly energy dependent in the diagnostic range. Though an optimization of the group structure is possible, it is problem (namely x-ray source spectrum) and organ dependent, which impose serious limitations of the deterministic solution for practical application in diagnostic medical physics.

AB - Monte Carlo methods have been established along the time as the gold standard for computer simulations in the medical physics community. Depending on the problem and user's needs, deterministic radiation transport simulations may provide a more detailed and faster solution. In this work we investigate the possibility of using deterministic radiation transport simulations as a viable and more convenient tool for real clinical applications. Therefore, the discrete ordinates PENTRAN code is used to calculate average organ doses in voxelized human phantoms and the results are compared with state-of-the-art MCNP5 Monte Carlo simulations in the diagnostic energy range (50-140 keV). Generally, good agreement for the average organ scalar fluxes, less than 6% difference, is obtained provided adequate quadrature order, mesh size and energy group structure is used in the deterministic calculations. The energy group structure, particularly for the diagnostic energy range, has a major impact on the deterministic solution for the average organ doses since the interaction and mass energy absorption coefficients are highly energy dependent in the diagnostic range. Though an optimization of the group structure is possible, it is problem (namely x-ray source spectrum) and organ dependent, which impose serious limitations of the deterministic solution for practical application in diagnostic medical physics.

KW - Average organ dose

KW - Deterministic discrete ordinates

KW - Monte Carlo

KW - Radiation transport

KW - Voxelized phantom

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

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

M3 - Conference contribution

SN - 9781615673490

T3 - American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009

SP - 1822

EP - 1834

BT - American Nuclear Society - International Conference on Mathematics, Computational Methods and Reactor Physics 2009, M and C 2009

ER -