|Program name||Package id||Status||Status date|
|Package ID||Orig. computer||Test computer|
|IAEA1382/04||PC Windows||Linux-based PC,PC Windows|
SRNA-2020 performs Monte Carlo transport simulation of proton in 3D source and 3D geometry of arbitrary materials. The proton transport based on condensed history model, and on model of compound nuclei decays that creates in nonelastic nuclear interaction by proton absorption.
The SRNA-2020 package is developed for time independent simulation of proton transport by Monte Carlo techniques for numerical experiments in complex geometry, using PENGEOM from PENELOPE with different material compositions, and arbitrary spectrum of proton generated from the 3D source. This package developed for 3D proton dose distribution in proton therapy and dosimetry, and it was based on the theory of multiple scattering. The compound nuclei decay was simulated by our and Russian MSDM models using ICRU 49 and ICRU 63 data. If protons trajectory is divided on great number of steps, protons passage can be simulated according to Berger's Condensed Random Walk model. The conditions of angular distribution and fluctuation of energy loss determinate step length.
Physical picture of these processes is described by stopping power, Moliere's angular distribution, Vavilov's distribution with Sulek's correction per all electron orbits, and Chadwick's cross sections for nonelastic nuclear interactions, obtained by his GNASH code. According to physical picture of protons passage and with probabilities of protons transition from previous to next stage, which is prepared by SRNADAT program, simulation of protons transport in all SRNA codes runs according to usual Monte Carlo scheme.
The proton from the spectrum prepared for random choice of energy, position and space angle is emitted from the source; proton is loosing average energy on the step; on that step, proton experience a great number of collisions, and it changes direction of movement randomly chosen from angular distribution; random fluctuation is added to average energy loss; protons step is corrected with data about protons position before and after scattering; there is final probability on step for nonelastic nuclear interaction to happen, and for proton to be absorbed.
Compound nucleus decays with emission of protons, neutrons, deuterons, tritons, alpha particles or photons. Particular decay particle is sampled from Poisson's distribution with appropriate average values of multiplication factor of each particle. Energy and angle of particle emission and factors of multiplication are determined from the cross section that obtained by the integration of differential cross section for nonelastic nuclear interaction. Energy and angle of secondary neutron are sampled from emission spectrum. Neutron and photon transport are not included in the current model. They are registered in data file and can be used by other code to simulate their transport. Emitted deuteron, triton and alpha particles are absorbed at the place their creation.
Proton kinetic energies have to be in the range from 100 keV to 250 MeV. No more than 128 geometry zones and less than 32 materials with no more than 15 elements in each material. In these conditions user can obtain geometry image by gview2d.exe from PENELOPE code.
The adopted parameters (energy cutoffs, geometry zones, etc.) have an influence on the computing time. As an example, a pencil beam depth-dose distribution of 250 MeV protons incident on a water phantom, obtained by simulating 100.000 histories, can be obtained with a running time of 1.4 minutes on a DEL 755, IntelCeleron 430 1.8GHz, 2GB DDR2.
SRNADAT auxiliary program to generate probabilities distributions of proton, gview11.exe from PENELOPE (http://www.oecd-nea.org/tools/abstract/detail/nea-1525/) for viewing a geometry image, the GRAF program for output data preparation and package VOXELVIEW (not free) for images of results presentations.
|Package ID||Status date||Status|
|IAEA1382/04||07-JUN-2020||Tested at NEADB|
- Radovan D. Ilić et al: The Monte Carlo SRNA-VOX Code for 3D Proton Dose Distribution in Voxelized Geometry Using CT Data, Phys. Med. Biol. 50 (2005) 1011-1017 (Feb. 2005)
No specific requirements.
Tested at the NEA Data Bank with:
- COMPUTER: AMD A6-9210 RADEON R4 2.40 GHz CPU Processor, RAM: 4.0 GB
- OPERATING SYSTEM: Microsoft Windows 10
- COMPILER: MinGW GNU Fortran and C/C++ Compiler v9.2.0
- COMPUTER: Intel® Core™ NUC7i7BNH i7-7567U processor 4.0 GHz Turbo Dual Core CPU, RAM: 16.0 GB
- OPERATING SYSTEM: Ubuntu 18.04 LTS
- COMPILER: gfortran v7.4, gcc v7.5
|Package ID||Computer language|
Keywords: Monte Carlo method, geometry, protons, three-dimensional, transport.