Overview
A python toolkit for efficiently sampling fusion products (neutrons and alpha particles) from parameterised or numerically calculated tokamak equilibria with parameterised kinetic profiles.
Fusion product sampling tools details
This python toolkit allows efficient Monte-Carlo sampling of fusion products (neutrons and alpha particles) using either parameterised (Miller) equilibria or numerical equilibria read from a GEQDSK file. The code uses parameterised temperature and density profiles and can use either a Sadler-Van-Belle or Bosch-Hale reactivity profile. The products are sampled isotropically in velocity space and subject to Doppler broadening in energy.
The code can run directly to generate a sample of fusion products which are then output in a pickle file. It can also generate and write cumulative distribution functions (CDFs). These can then be natively sampled from within the OpenMC neutronics code. Full details of the numerical methods employed are in section 3.1 of the upcoming paper “Full Tokamak Device Model for Non-linear Effects of Fusion Neutron Irradiation” (L.Reali et al).
Who it’s for
The Fusion product (neutron/alpha) sampling tools package is for Fusion scientists studying neutronics or fast particle physics in tokamaks.
Licence
MIT Licence
Release date
7 June 2026
Programming languages
- Python (Main scripts)
- C++ (OpenMC coupling class + requisites)
Dependencies
- Numpy
- math
- matplotlib
- scipy
- fortranformat
- mpl_toolkits
- json
- pickle
Compatible platforms / systems
All