RSICC CODE PACKAGE CCC-734
1. NAME AND TITLE
TRIPOLI-4 version 4.3.3: Code System for Coupled Neutron, Photon, Electron, Positron, 3-D, Time Dependent, Monte-Carlo, Transport Calculations.
DATA LIBRARIES
Related libraries in a TRIPOLI-4 specific format based on JEF2, ENDFB6R4, and ENDL are included. These libraries were produced using NJOY.
2. CONTRIBUTORS
Commissariat à l'énergie atomique, CEA/SACLAY, Cedex, France, through the OECD Nuclear Energy Agency Data Bank, Issy-les-Moulineaux, France.
3. CODING LANGUAGE AND COMPUTER
Fortran 77 and C; SUN, IBM, HP, Digital, SGI, Cray, Compaq and Linux based PCs
(RSICC ID: C00734MNYCP00). (NEADB identifier is NEA-1716/02.)
4. NATURE OF PROBLEM SOLVED
TRIPOLI 4.3 solves the transport equation for neutral particles in general three-dimensional geometrical configurations. The following problems can be treated:
The TRIPOLI-4 code is used essentially for four major classes of applications: shielding studies, criticality studies, core physics studies and instrumentation studies. It computes particle fluxes and currents and several related physical quantities such as, reaction rates, dose rates, heating, energy deposition, effective multiplication factor, perturbation effects due to density, concentration or partial cross-section variations.
Types of particles Neutron, gamma, electron, positron.
Nuclear data ENDF-6 format. Cross sections processed
by NJOY code system. Pointwise cross
section representation. Probability
tables (PT) representation possible in
unresolved energy range (PT produced by
CALENDF code).
Energy ranges For neutrons and gamma the energy range
is the same as range specified into the
evaluation in ENDF-6 format.
Neutrons : 0. - 150 MeV
Gamma : 0 - 100 MeV
For Electrons/positrons :
1 MeV - several GeV (with bremsstrahlung
gamma rays).
Geometry 3-D surface and combinatorial, network
and network of networks.
Sources General factorised description: space,
energy, angle, time
Calculated physical See table below.
quantities, estimators
Biasing Exponential Transform,
splitting/roulette. Automatisation of
biasing.
Time dependent transport for neutrons.
Perturbation Density, concentration of isotopes,
partial cross-sections.
Coupled particle (neutron, gamma), (electron, positron,
transport gamma)
Qualification Shielding (SINBAD benchmarks),
criticality (ICSBEP benchmarks)
------------------------------------------------------------------
Estimated physical quantities Type of estimator
volume surface point
Total flux Track length yes yes
and collision
Integrated flux Track length yes yes
per energy group and collision
and time interval
Integrated flux Track length yes yes
per energy group and and collision
angular mesh
Integrated flux per Track length yes yes
energy group and collision
Total current yes
Integrated current yes
per energy group and
time interval
Integrated current yes
per energy group and
angular mesh
Integrated flux per yes
energy group
Reaction rate on Track length yes
whole energy range and collision
Reaction rate per Track length yes
energy group and collision
time interval
Reaction rate per Track length yes
energy group collision
Deposited energy collision
Equivalent dose rate Track length yes yes
collision
Gamma production yes
Multiplication factor yes (kstep, kcoll, ktrack, kij)
Eigen value of a yes
criticality
configuration
Cross section yes
(partial or total)
perturbations
Density or yes
concentration
perturbation
Data libraries distributed with the TRIPOLI-4.3:
- ENDFB6R4: neutron, gamma data library
- ENDL: gamma data library
- JEF2: neutron, gamma data library
- Mott-Rutherford: electron, positron cross-section library
- Qfission: energy release during fission library
5. METHOD OF SOLUTION
The geometry package allows the user to describe a three dimensional configuration by means of surfaces (as in the MCNP code) and also through predefined shapes combine with operators (union, intersection, subtraction...). It is also possible to repeat a pattern to build a network of networks.
TRIPOLI-4 can use four different types of cross-sections representation:
Neutron and photon transport is governed by the Boltzmann equation. The code resolves this equation by the Monte Carlo method. This requires a random generator: TRIPOLI4 uses drand48 in non-parallel calculations and gfsr for parallel calculations. The calculation is carried out by division into several batches of particles, each batch containing an equal number of source particles.
Within a batch, the code monitors particles in the following manner:
1) The source particles are first generated by random selection according to the distribution supplied by the user (in space, energy, angle and time).
2) Each particle is monitored individually:
a. Between two impacts the trajectory is a straight line (since they are neutral particles). The distance between two impacts depends on the total cross-section of the medium.
b. When an impact occurs, sampling concerns the impacted isotope (as a function of the concentration), then the type of interaction (as a function of the isotope) and finally the energy and direction after impact (as a function of the interaction and initial energy).
The mean results for each quantity required (based on the number of source particles per batch) are calculated over all the batches, each of which can be considered to be a sample. The final scores are obtained by calculating the average for the number of batches making up the total results. Standard deviation can then be calculated for the final values. If the number of batches is low and the central limit theorem does not apply, the code can use suitable bootstrap techniques.
If there is high attenuation, biasing techniques must be used to avoid excessive calculation times. In this case, the code performs an exponential transformation. (An importance factor is assigned to each point of phase space so that biased cross-sections can be used).
Several biasing techniques are implemented in TRIPOLI-4: exponential biasing scheme, quota sampling and collision biasing. The biasing is automatised from a Dijkstra algorithm.
The communication library has been written to achieve message passing between the different processes involved in the simulation. This library only cares of data that are exchange between processes while the parallelism library deals with actions that must be achieved by the processes.
6. RESTRICTIONS OR LIMITATIONS
None noted.
7. TYPICAL RUNNING TIME
Run times vary. Test cases completed in approximately 20 minutes.
8. COMPUTER HARDWARE REQUIREMENTS
TRIPOLI-4 runs on Unix workstation (SUN, IBM, HP, Digital, SGI); on Linux based PCs; and on massively parallel machines (Cray, Compaq).
9. COMPUTER SOFTWARE REQUIREMENTS
TRIPOLI runs under Linux or Unix workstations. X11 is required. Fortran and C compilers are required if one cannot use the developers’ executables which are included for the systems listed below. A Fortran to C converter (F2C) may be required for some Linux systems. At RSICC, TRIPOLI was tested by running the developers’ executables on an AMD Opteron under Red Hat Enterprise Linux WS release 4 version 2.6.9-34 and on IBM RS/6000 under AIX Version 5.1. At NEADB, TRIPOLI was tested on a Pentium III under RedHat Linux 7.2 with f77 vers. 0.5.26 and gcc vers. 2.96.
The developers compiled and tested on the following architectures. Executables for each of these systems is included in the package.
* Solaris-sparc (5.7)
+ GNU make
+ cc, CC, f77 of SUNWspro6.2
* linux-intel (2.2.xx or 2.4.xx)
+ GNU make
+ gcc, g++ version 2.95 or 2.96
+ f2c-19991025
* aix (5.1)
+ GNU make
+ xlC, xlc, xlf version 5
* hpux (B10.20)
+ GNU make
+ aCC, c89, f77
* osf1 (V4.0)
+ GNU make
+ cxx, cc, f77
10. REFERENCES
10.a included in package in electronic files:
- J.-P. Both, A. Mazzolo, Y. Peneliau, O. Petit and B. Roesslinger:
User Manual for Version 4.3 of the TRIPOLI-4 Monte-Carlo Method Particle
Transport Computer Code
ISSN 0429 - 3460 Rapport CEA-R-6044 (November 2003)
- J.-P. Both, A. Mazzolo, Y. Peneliau, O. Petit and B. Roesslinger:
Notice d'Utilisation du Code TRIPOLI-4 Version 4.3: Code de Transport de
Particules par la Methode de Monte Carlo
ISSN 0429 - 3460 Rapport CEA-R-6043 (Novembre 2003)
10.b references concerning theory and functionalities - not included in pkg
J. P. Both, A. Mazzolo, O. Petit, Y. Peneliau, B. Roesslinger, User Manual for
version 4.3 of the TRIPOLI-4 Monte Carlo method particle transport computer
code, CEA-Report : CEA-R-6044, DTI, CEA/Saclay, France, 2003
J.P. Both, H. Derriennic, B. Morillon, J.C. Nimal, " A Survey of TRIPOLI-4 ",
Proceedings of the 8th International Conference on Radiation Shielding,
Arlington, Texas, USA, 24-28 avril 1994, pp. 373-380.
B. Morillon, " Methode de Monte Carlo non analogue - Application à la
simulation des neutrons ", Note CEA-N-2805, CEA/Saclay, janvier 1996.
J.P. Both, Y.K. Lee, " Computations of Homogenised Multigroup Cross Sections
with the Monte Carlo Code TRIPOLI 4. ", Proceedings SARATOGA 1997, JIC MM
&SNA, Saratoga Springs, New York, USA, october, 1997, p. 439.
Y. Peneliau, J.P. Both, " Parallelization of the Monte Carlo Code TRIPOLI-4 ",
Mathematical and Computation, Reactor Physics and Environmental Analysis in
Nuclear Applications, Madrid, Spain, Sept. 1999, p. 412.
J. P. Both, Treatment of Cross Section Uncertainties in the Transport Monte
Carlo Code TRIPOLI-4, ICRS'9, Tsukuba, Ibaraki, Japan, october 17-22, 1999,
Journal of Nuclear Science and Technology, Supplement 1, pp. 420-42, March
2000, pp. 420-426.
Y. Peneliau, " Electron Photon Shower Simulation TRIPOLI-4 in Monte Carlo Code
", Advanced Monte Carlo for Radiation Physics, Particle Transport Simulation
and Applications, Lisbon, Portugal, Oct. 2000.
J. P. Both, Y.K. Lee, A. Mazzolo, O. Petit, Y. Peneliau, B. Roesslinger, M.
Soldevila, TRIPOLI-4 - A Three Dimensional Polykinetic Particle Transport
Monte Carlo Code, SNA'2003, Paris Sept. 2003.
10.c references concerning qualification - not included in pkg
Y.K. Lee, G. Neron, J.P. Both, Y. Peneliau, C. Diop, " Validation of Monte
Carlo Code TRIPOLI-4 with PWR Critical Lattices by using JEF 2.2 and ENDF/B-VI
Evaluations ". JIC MM &SNA, Saratoga Springs, New York, USA, October, 1997, p.
253.
Y.K. Lee, S.H. Zheng, G. Neron, J.P. Both, Y. Peneliau, C. Diop, " ICNC'9,
Cristal, Criticality Safety Package Validation : Tripoli-4 Monte Carlo Code,
JEF2.2 Library and ICSBEP Experiments ", Sixth International Conference on
Nuclear Criticality Safety. Versailles, France, 20-24 September 1999.,
Y. K. Lee, « Analysis of the NRC PCA Pressure Vessel Dosimetry Benchmark Using
TRIPOLI-4.3 Monte Carlo Code and ENDF/B-VI, JEF-2.2 and IRDF-90 Libraries,
Nuclear Mathematical and Computational Science », Gatlinburg, Tennessee, USA,
April 6-11 2003.
E. Gagnier, Y. K. Lee, L. Aguiar, N. Védrenne, Validation of the 3D Transport
Monte Carlo Code TRIPOLI-4.3 for moderated and unmoderated metallic fissile
media configurations with Jef-2.2 and ENDF/B6.r4 cross section evaluations ;
Y. K. Lee, Analysis of the LEU-COMP-THERM-049 Maracas Critical Configurations
Using TRIPOLI-4.3 3D Lattices Geometry and JEFF-3.0 Library, ICNC'2003,
Tokai-Mura, October 2003.
Y. K. Lee, Analysis of the LEU-COMP-THERM-049 Maracas Critical Configurations
Using TRIPOLI-4.3 3D Lattices Geometry and JEFF-3.0 Library.
11. CONTENTS OF CODE PACKAGE
The package is transmitted on a CD which contains the reports listed in section 10.a above, source code, scripts, executable files, data files and test cases in a Unix tar file.
12. DATE OF ABSTRACT
March 2007.
KEYWORDS:
COMPLEX GEOMETRY; COUPLED NEUTRON-GAMMA-RAY; ELECTRON CROSS SECTIONS; CRITICALITY CALCULATIONS; MONTE CARLO; NEUTRON; POSITRON; TIME-DEPENDENT