last modified: 21-JAN-2010 | catalog | categories | new | search | indexing |

NEA-1833 FLUKA2008.3C.

FLUKA2008.3c, Monte Carlo general purpose tool for calculations of particle transport and interactions with matter

top ]
1. NAME OR DESIGNATION OF PROGRAM:  FLUKA2008.3c
top ]
2. COMPUTERS
To submit a request for this package, please contact the liaison officer in your institution. Rules for requesters are available here.
Program name Package id Status Status date
FLUKA2008.3C NEA-1833/04 Arrived 21-JAN-2010

Machines used:

Package ID Orig. computer Test computer
NEA-1833/04 Linux-based PC
top ]
3. DESCRIPTION OF PROGRAM OR FUNCTION

This package contains:
FLUKA2008.3c
FLAIR-0.7-9.
  
FLUKA2008.3c
============
  
FLUKA is a Fortran 77 (f77) Monte Carlo computer code.
It is a general purpose tool for calculations of particle transport and interactions with matter, covering an extended range of applications spanning from proton and electron accelerator shielding to target design, calorimetry, activation, dosimetry, detector design, Accelerator Driven Systems, cosmic rays, neutrino physics, radiotherapy etc.
  
For details, see:
G. Battistoni, S. Muraro, P.R. Sala, F. Cerutti, A. Ferrari, S. Roesler, A. Fasso`, J. Ranft:
"FLUKA: a multi-particle transport code"
Proceedings of the Hadronic Shower Simulation Workshop 2006, Fermilab 6--8 September 2006, M. Albrow, R. Raja eds., AIP Conference Proceeding 896, 31-49, (2007)
and
"The FLUKA code: Description and benchmarking" A. Ferrari, P.R Sala, A. Fasso`, and J. Ranft, CERN-2005-10 (2005), INFN/TC_05/11, SLAC-R-773
  
website: http://www.fluka.org/

FLUKA is a general purpose Monte Carlo radiation transport code that tracks nearly all particles over an extended energy range. Various generations of FLUKA can be distinguished. Initially (1962-1987) the code, originated at Leipzig University and later modified in collaboration with CERN and Helsinki University of Technology, was transporting exclusively hadrons of energies larger than 50 MeV.  The present multi-particle, extended energy code started in 1989 as an effort mainly of INFN (Italy) and has been in continuous development since then, eventually in the frame of a collaboration between INFN and CERN which was formalized in 2003. A complete history is reported on the FLUKA manual, distributed with the code and available on line on the FLUKA web site.
  
The physical models of FLUKA include hadron-hadron and hadron-nucleus interactions and transport up to 10 PeV, nucleus-nucleus interactions and transport between 100 MeV/n and
10 PeV/n, electron, positron and photon interactions and transport between 1 keV and 10 PeV, neutrino interactions, neutron multigroup transport and interactions up to 20 MeV, charged particle transport including all relevant processes and transport in magnetic fields.
  
The program can be run in analog mode or with several variance reduction options.  Geometry description can be done with an advanced version of Combinatorial Geometry allowing lattice capabilities and voxel description.
The code has been extensively benchmarked and has found a large number of applications, such as cosmic ray physics, neutrino physics, accelerator design, calorimetry and particle detector simulation, shielding design, dosimetry and radiation protection, space radiation, hadron therapy, neutronics, ADS systems, waste transmutation, etc.
  
Information about FLUKA development can be found on the web site http://www.fluka.org/
  
The FLUKA package for LINUX platforms is distributed as two separated files. The release of the FLUKA source code is available under the licence established by the FLUKA Coordination Committee.
  
The distribution consists in a package containing a compiled library, user routines in source form, INCLUDE files, various unformatted and formatted data files and a number of scripts for compiling, linking and running the program on a given platform. A list of the contents is provided in a README file, and information on the current version, possibly overriding parts of the current manual, may be contained in a file RELEASE-NOTES.
A second package contains the source files.
No external library routines are required. The timing and other necessary service routines are already included.
NEA-1833/04
Changes compared to previous version:
------------------------------------
FLUKA2008.3c contains few bug fixes, in particular:
  
a) Tau+/- decays (whether "prompt" or after transport) are now controlled  by the PHYSICS card with SDUM=CHARMDEC as per charmed particle decays
  
and a significant functionality change/improvement which affects users of FLUGG and VMC (and which is the reason for bumping the patch level).
  
Functionality change with respect to the use of FLUGG and VMC (Virtual Monte Carlo):
  
- FLUGG: starting from this patch version, FLUGG is better supported by Flair. Most Flair functionalities, including geometry plotting, can be used with FLUGG, having care to point to the correct executable.
- VMC: due to several examples of (inadvertent) use of the FLUKA-VMC interface by people/groups who never asked for the required derogation to the license, the hooks for VMC have been removed.
top ]
4. METHODS

FLUKA treats an arbitrary three-dimensional configuration of materials in geometric regions bounded by first- and second-degree surfaces.  Derived from the Combinatorial Geometry package, it has been entirely rewritten.  A completely new, fast tracking strategy has been developed, with special attention to charged particle transport, especially in magnetic fields. New bodies have been introduced, resulting in increased rounding accuracy and speed. Input preparation has been made much easier by the possibility to use names instead of numbers, free format and nested parentheses. The distance to nearest boundary is taken into account for improved performance.
Repetitive structures (lattices) and voxel geometries can also be handled.
  
The FLUKA hadron-nucleon interaction models are based on resonance production and decay below a few GeV, and on the Dual Parton model above. Two models are used also in hadron-nucleus interactions. At momenta below 3-5 GeV/c the PEANUT package includes a very detailed Generalised Intra-Nuclear Cascade (GINC) and a preequilibrium stage, while at high energies the Gribov-Glauber multiple collision mechanism is included in a less refined GINC. Both modules are followed by equilibrium processes: evaporation, fission, Fermi break-up, gamma deexcitation. FLUKA can also simulate photonuclear interactions (described by Vector Meson Dominance, Delta Resonance, Quasi-Deuteron and Giant Dipole Resonance). The PEANUT model is set to become the default at all energies: please read the release notes for further details about this possibility.
  
Hadron elastic scattering is described by means of parameterised nucleon-nucleon cross sections, tabulated nucleon-nucleus cross sections and tabulated phase shift data for pion-proton and phase-shift analysis for kaon-proton scattering.  Detailed kinematics of elastic scattering is performed on hydrogen nuclei and transport of proton recoils.
  
Nuclear interactions ge
nerated by ions are treated through interfaces to external event generators: DPMJET-2.5 or DPMJET-3, with a special initialisation procedure, above 5 GeV per nucleon, modified Rqmd-2.4 between 0.1 and 5 GeV per nucleon, and BME (Boltzmann Master Equation) below 0.1 GeV per nucleon. Pre-compiled libraries for these event generators are included in the distributed packages: the source code is not yet included, pending finalization of proper licensing.
  
Particle transport includes time-dependence.
Transport of charged particles is based on an original treatment of multiple Coulomb scattering and of ionisation fluctuations which allows the code to handle accurately some challenging problems such as electron backscattering and energy deposition in thin layers even in the few keV energy range.
Energy loss of charged particles is based on the Bethe-Bloch theory, with optional delta-ray production and transport with account for spin effects and ionisation fluctuations. Shell and other low-energy corrections are derived from Ziegler, the density effect is according to Sternheimer.
For all charged particles (hadrons and muons as well as electrons and positrons) a special transport algorithm, based on Moliere's theory of multiple Coulomb scattering improved by Bethe, accounts for correlations between lateral and longitudinal displacement and the deflection angle, between projected angles, and between projected step length and total deflection.
The algorithm includes an accurate treatment of boundaries and curved trajectories in magnetic fields, an automatic control of the step, a path length correction, spin-relativistic effects at the level of the second Born approximation, nuclear size effects (scattering suppression) on option, and a correction for cross section variation with energy over the step.
  
Bremsstrahlung and electron pair production at high energy by heavy charged particles are treated as a continuous energy loss and deposition or
as discrete processes depending on user choice. Muon photonuclear interactions are simulated with or without transport of the produced secondaries.
  
Differences between positrons and electrons are taken into account concerning both stopping power and bremsstrahlung. The bremsstrahlung differential cross sections of Seltzer and Berger have been extended to include the finite value at "tip" energy, and the angular distribution of bremsstrahlung photons is sampled accurately. The Landau-Pomeranchuk-Migdal suppression effect and the Ter-Mikaelyan polarisation effect in the soft part of the bremsstrahlung spectrum are also implemented.

Positron annihilation is simulated both in flight and at rest. Delta-ray production by positrons and electrons is described via Bhabha and Moller scattering.  The lowest transport limit for electrons is 1 keV. Although in high-Z materials the Moliere multiple scattering model becomes unreliable below 20-30 keV, a single-scattering option is available which allows to obtain satisfactory results in any material also in this low energy range.
  
Photon interactions include pair production with actual angular distribution of electrons and positrons, Compton effect with account for atomic bonds through use of inelastic Hartree-Fock form factors, photoelectric effect with actual photoelectron angular distribution, detailed interaction on six K and L single sub-shells, optional emission of fluorescence photons and an approximate treatment of Auger electrons, and Rayleigh effect. Photon polarisation can be taken into account for Compton, Rayleigh and photoelectric effects. Photohadron production is modelled according to the Vector Meson Dominance Model, modified and improved using PEANUT below 770 MeV, Quasideuteron interactions and Giant Dipole Resonance. Photomuon production is described according to Tsai.
  
For neutrons with energy lower than 20 MeV, FLUKA uses its own neutron cross section libraries (P5 Legendre
angular expansion, 260 or 72 neutron energy groups), containing more than 200 different materials, selected for their interest in physics, dosimetry and accelerator engineering and derived from the most recently evaluated data. Gamma-ray generation and different temperatures are available. Doppler broadening is applied for temperatures above 0 K.
The neutron transport is based on standard multigroup transport with photon and fission neutron generation, detailed kinematics of elastic scattering on hydrogen nuclei, transport of proton recoils and protons from 14-N(n,p)14-C reaction.  Capture photons are generated according to the multigroup treatment, but transported with the more accurate electromagnetic package of FLUKA which performs continuous transport in energy and allows for secondary electron generation. The 2.226 MeV gamma line from capture in hydrogen is generated as an actual precise energy, and the whole gamma cascade is available for Xenon and Cadmium isotopes. For nuclei other than hydrogen, kerma factors are used to calculate energy deposition (including from low-energy fission).  Pointwise cross section transport is available for a few nuclei and reactions.
  
Electron, muon, and tau (anti)neutrinos are produced and tracked on option, without interactions, but neutrino interactions are implemented, independently from tracking.
  
Generation and transport are available (on user's request) of Cherenkov and scintillation radiation. Transport of light of given wavelength in materials can be simulated with user-defined optical properties.
  
FLUKA has extended scoring capabilities, requiring in most cases no user-written code. Quantities which can be scored include star (hadron inelastic interaction) density by producing particle and region, energy density by region, total or from electrons/photons only, and energy and momentum transfer density in a geometry-independent binning structure (Cartesian or cylindrical), averaged ove
r the run or event by event. The step size is independent of bin size.  Energy deposition can be weighted by a quenching factor (Birks law).  Scoring can be done in a time window. It is possible to simulate coincidences and anti-coincidences.
Fluence and current can be scored as a function of energy and angle, via boundary-crossing, collision and track-length estimators coincident with regions or region boundaries.
Track-length fluence can be scored in a binning structure (Cartesian or cylindrical) independent of geometry.
Particle yield from a target is available, or differential cross section with respect to several different kinematic variables.
Other scoring possibilities include residual nuclei, fission density, momentum transfer density, neutron balance, unweighted energy deposition.
All quantities from radioactive decay of residual nuclei can be scored according to user-defined irradiation and cooling time profiles (decay radiation transport is provided on request).
  
FLUKA can be used in analog mode or with a variety of variance reduction options.  These include: Leading particle biasing for electrons and photons: region dependent, below user-defined energy threshold and for selected physical effects; Russian Roulette and splitting at boundary crossing based on region relative importance; region-dependent multiplicity tuning in high energy nuclear interactions; region-dependent biased downscattering and non-analog absorption of low-energy neutrons; biased decay length for increased daughter production, biased inelastic nuclear interaction length; biased interaction lengths for electron and photon electromagnetic interactions; biased angular distribution of decay secondary particles; region-dependent weight window in three energy ranges (and energy group dependent for low energy neutrons).
top ]
5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

Neutrons in the FLUKA low energy libraries are available for about 200 materials or isotopes, temperature, and self-shielding combinations. All other particle interactions and transport are based on models and are not restricted by any material tabulation.  The upper energy limit for hadron-hadron and hadron-nucleus interactions and transport is 10 PeV, for nucleus-nucleus interactions and transport 10 PeV/n, when the interface with DPMJET-2.5 or DPMJET-3 is activated, 100 TeV otherwise. Electron, positron and photon interactions and transport are possible between 1 keV and 10 PeV.
top ]
6. TYPICAL RUNNING TIME

Runtimes vary greatly depending on computer speed, problem layout, and beam energy. On a 2.4 GHz Centrino CPU the compilation of FLUKA takes about 15 minutes.
top ]
7. UNUSUAL FEATURES
top ]
8. RELATED OR AUXILIARY PROGRAMS

RELATED DATA LIBRARIES

DATA LIBRARIES included in the distribution
===========================================
DPMJET-2.5 and DPMJET-3 libraries. rQMD-2.4 library.
Data files: Bremsstrahlung cross sections, Coherent atomic form factors, Fluorescence emission data, Photon cross sections, Low-energy neutron cross sections (72 and 260 groups), Nuclide masses, abundances and other data, Hadron elastic cross sections, Pion cross sections, Fission nuclide yields and neutron multiplicities, Silicon Damage tabulations.
  
AUXILIARY PROGRAMS included in the distribution
===============================================
Post-processing programs to analyse the user output.
The FLAIR graphical interface.
top ]
9. STATUS
Package ID Status date Status
NEA-1833/04 21-JAN-2010 Masterfiled restricted
top ]
10. REFERENCES

- G. Battistoni, S. Muraro, P.R. Sala, F. Cerutti, A. Ferrari, S. Roesler, A. Fasso', J. Ranft:
"FLUKA: a multi-particle transport code"
Proceedings of the Hadronic Shower Simulation Workshop 2006, Fermilab 6--8 September 2006, M. Albrow, R. Raja eds., AIP Conference Proceeding 896, 31-49, (2007)
NEA-1833/04, included references:
- Alfredo Ferrari, Paola R. Sala, Alberto Fasso', Johannes Ranft:
Fluka: a multi-particle transport code (Program version 2008)
CERN 2005-10 (2005), INFN/TC_05/11, SLAC-R-773
top ]
11. HARDWARE REQUIREMENTS

FLUKA runs under Linux operating system. It has been tested by the authors on various 32-bit and 64-bit workstations and personal computers. The compiled version of the code requires 370 Mbytes of memory.
top ]
12. PROGRAMMING LANGUAGE(S) USED
Package ID Computer language
NEA-1833/04 C-LANGUAGE, FORTRAN-77
top ]
13. SOFTWARE REQUIREMENTS

The gcc and g77 compilers are required to compile and link FLUKA. The GNU make utility is required to build the system on Linux platforms.
top ]
14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS

Users are bound to run the code only on the platforms and with the compiler, options included, approved by the authors: at present Linux operating systems with g77 (see the documentation for more details).
  
An self-contained image for running FLUKA in a virtual environment on Windows machines is available on the web site http://www.fluka.org/
top ]
15. NAME AND ESTABLISHMENT OF AUTHORS

Contributed by:
CERN, the European Organization for Nuclear Research,
CH-1211, Geneva 23
Switzerland
  
and
  
INFN, the National Institute of Nuclear Physics,
P.za dei Caprettari 70
00186 Rome
Italy
  
Developed by:
A. Fasso', SLAC National Accelerator Laboratory, USA
A. Ferrari, CERN
P.R. Sala, INFN
J. Ranft, Siegen University, Siegen, Germany
top ]
16. MATERIAL AVAILABLE
NEA-1833/04
readme and release notes
fluka2008.3c-linux_source.tar.gz contains the source files
fluka2008.3c-linuxAA.tar.gz contains:
libflukahp.a     compiled library with all Fluka routines
example inputs in various formats
random.dat       random number seed
datafiles
flutil/ contains:
lfluka  script to link Fluka (contains help) or auxiliary programs
scripts to link Fluka
programs to sum data files
flukapro/ contains all Fluka include files
usermvax/ contains fortran user routines
FM.pdf         manual, pdf version
CERN-2005-10.pdf  official reference
FLAIRflair-0.7-9
top ]
17. CATEGORIES
  • J. Gamma Heating and Shield Design

Keywords: ADS, Monte Carlo method, activation, antiparticles, biasing, calorimetry, charged particles, coincidences, complex geometry, cosmic rays, dosimetry, electrons, gamma ray, hadrotherapy, high energy, kaon, magnetic fields, neutrino, neutron, optical photon, pion, protons, radiation transport, shielding, spallation, time dependence.