Radiation physics, environmental contamination, and computational science. Researcher working for the Japan Atomic Energy Agency, in Kashiwa, Japan

Spotted an error on this site?
Corrections appreciated, please email alex (at) alexmalins.com or leave a comment.

Continuous energy adjoint transport for photons in PHITS

Calculation of a gamma ray spectrum using conventional forward and adjoint Monte Carlo method


Adjoint Monte Carlo can be an effcient algorithm for solving photon transport problems where the size of the tally is relatively small compared to the source. Such problems are typical in environmental radioactivity calculations, where natural or fallout radionuclides spread over a large area contribute to the air dose rate at a particular location. Moreover photon transport with continuous energy representation is vital for accurately calculating radiation protection quantities. Here we describe the incorporation of an adjoint Monte Carlo capability for continuous energy photon transport into the Particle and Heavy Ion Transport code System (PHITS). An adjoint cross section library for photon interactions was developed based on the JENDL- 4.0 library, by adding cross sections for adjoint incoherent scattering and pair production. PHITS reads in the library and implements the adjoint transport algorithm by Hoogenboom. Adjoint pseudo-photons are spawned within the forward tally volume and transported through space. Currently pseudo-photons can undergo coherent and incoherent scattering within the PHITS adjoint function. Photoelectric absorption is treated implicitly. The calculation result is recovered from the pseudo-photon flux calculated over the true source volume. A new adjoint tally function facilitates this conversion. This paper gives an overview of the new function and discusses potential future developments.

A. Malins, M. Machida & K. Niita
EPJ Web of Conferences 153, 06001 (2017)

DOI: 10.1051/epjconf/201715306001
PDF: download

Page last updated: 27th September 2017