last modified: 18-MAY-2011 | catalog | categories | new | search |

IAEA1382 SRNA-2K5.

SRNA-2K5, Proton Transport Using 3-D by Monte Carlo Techniques

top ]
1. NAME OR DESIGNATION OF PROGRAM:  SRNA-2K5.
top ]
2. COMPUTERS
To submit a request, click below on the link of the version you wish to order. Only liaison officers are authorised to submit online requests. Rules for requesters are available here.
Program name Package id Status Status date
SRNA-2K5 IAEA1382/03 Tested 18-MAY-2011

Machines used:

Package ID Orig. computer Test computer
IAEA1382/03 Linux-based PC Linux-based PC,PC Windows
top ]
3. DESCRIPTION OF PROGRAM OR FUNCTION

SRNA-2K5 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.
top ]
4. METHODS

The SRNA-2K5 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. 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:
(i) proton from the spectrum prepared for random choice of energy, position and space angle is emitted from the source;
(ii) proton is loosing average energy on the step;
(iii) on that step, proton experience a great number of collisions, and it changes direction of movement randomly chosen from angular distribution;
(iv) random fluctuation is added to average energy loss; (v) protons step is corrected with data about protons position before and after scattering;
(vi) 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.
top ]
5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

Proton kinetic energies have to be in the range from 100 keV to 250 MeV. No more then 128 geometry zones and less then 32 materials with no more then 15 elements in each material. In these conditions user can obtain geometry image by gview2d.exe.
top ]
6. TYPICAL RUNNING TIME

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 4.8 minutes on a Pentium IV 1.67 GHz 512 MB RAM.
top ]
7. UNUSUAL FEATURES

Because of heaving no limitations in 3D proton simulation, SRNA-2K5 is suitable for numerical experiments in dosimetry, especially for protontherapy treatment planning calculations.
top ]
8. RELATED OR AUXILIARY PROGRAMS

SRNADAT auxiliary program to generate probabilities distributions of proton, gview2d.exe from PENELOPE for viewing a geometry image, and the GRAF program for output data preparation.
top ]
9. STATUS
Package ID Status date Status
IAEA1382/03 18-MAY-2011 Tested at NEADB
top ]
10. REFERENCES
IAEA1382/03, included references:
- Radovan D. Ilic:
SRNA-2K5 - Protons Transport Simulation by Monte Carlo Techniques
Version 2K5, March, 2005 (User's Guide)
- Radovan D. Ilic et al.:
SRNA - Monte Carlo Codes for Proton Transport Simulation in Combined and
Voxelized Geometries
NT&RP 2002 XVII 1-2 p. 27-36 http://ntrp.vin.bg.ac.yu/
- Radovan D. Ilic 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)
top ]
11. HARDWARE REQUIREMENTS

There are no specific requirements concerning the SRNA codes.
IAEA1382/03
This package was tested at the NEA Databank on:
- COMPUTER : apple iMac i7 2800MHz
             IBM T61p Intel Core2Duo T7700 @2400MHz
- OPERATING SYSTEM : MS Windows XP Professional
                     ubuntu Linux 10.4 Kernel 2.6
- COMPILER : Compac Visual Fortran 6.1
             GNU gfortran Version 4.4.3 (based on gcc 4.4.3)
top ]
12. PROGRAMMING LANGUAGE(S) USED
Package ID Computer language
IAEA1382/03 FORTRAN-77
top ]
13. SOFTWARE REQUIREMENTS

FORTRAN 77 with MinGW under Windows or gcc 3.3 or later under Linux.
top ]
14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS
top ]
15. NAME AND ESTABLISHMENT OF AUTHORS

Radovan D. Ilic
Institute of Nuclear Sciences VINCA
Physics Laboratory (010)
11001 Beograd, P.O.Box 522
Serbia and Montenegro
WEB:   http://www.vin.bg.ac.yu/~rasa/
top ]
16. MATERIAL AVAILABLE
IAEA1382/03
SRNADOC\ Documentation directory including manual and references in PDF
SRLINPEN\ SRNA-2K5 files on the working folder
GEOMPEN\ PENGEOM examples and viewer
SRNATEST\ Input and output SRNA-2K5 test case  files
top ]
17. CATEGORIES
  • J. Gamma Heating and Shield Design

Keywords: Monte Carlo method, geometry, protons, three-dimensional, transport.