Characterization and optimization of irradiation facilities in a research reactor is important for optimal performance. Nowadays this is commonly done with advanced Monte Carlo neutron transport ...computer codes such as MCNP. However, the computational model in such calculations should be verified and validated with experiments. In the paper we describe the irradiation facilities at the JSI TRIGA reactor and demonstrate their computational characterization to support experimental campaigns by providing information on the characteristics of the irradiation facilities.
► TRIGA reactor at JSI suitable for irradiation under well defined conditions. ► It features irradiation channels of different fluxes, spectra, and dimensions. ► Computational model has been developed and experimentally verified. ► The model used for optimization of experiments and evaluation of uncertainties.
Abstract
Collection and interpolation of radiation observations is of vital importance to support routine operations in the nuclear sector globally, as well as for completing surveys during crisis ...response. To reduce exposure to ionizing radiation that human workers can be subjected to during such surveys, there is a strong desire to utilise robotic systems. Previous approaches to interpolate measurements taken from nuclear facilities to reconstruct radiological maps of an environment cannot be applied accurately to data collected from a robotic survey as they are unable to cope well with irregularly spaced, noisy, low count data. In this work, a novel approach to interpolating radiation measurements collected from a robot is proposed that overcomes the problems associated with sparse and noisy measurements. The proposed method integrates an appropriate kernel, benchmarked against the radiation transport code MCNP6, into the Gaussian Process Regression technique. The suitability of the proposed technique is demonstrated through its application to data collected from a bespoke robotic system used to conduct a survey of the Joz̆ef Stefan Institute TRIGA Mark II nuclear reactor during steady state operation, where it is shown to successfully reconstruct gamma dosimetry estimates in the reactor hall and aid in identifying sources of ionizing radiation.
The operation of a nuclear power plant relies on precalculated nuclear design predictions based on core calculations of various reactor states. The fuel temperature is a crucial factor in determining ...the reactor fuel behavior, but assessing the temperature variation in a fuel pellet taking into account neutron transport is challenging. Detailed simulation of the temperature behavior within the fuel pellet can be obtained by coupling of Monte Carlo neutron transport codes with thermal-hydraulics solvers. However, this approach is not practical for standard nuclear design calculations, and computationally cheaper and faster methods must be used. In nuclear core simulators, a concept of a single "effective temperature" that yields the same neutron response as in the case of the actual temperature shape is mainly applied. This paper evaluates various fuel temperature models used in nuclear core simulation calculations, ultimately recommending a new effective temperature model that considers the burnup correction.
Reactivity worth of fuel rods at the JSI TRIGA research reactor was measured. Differently burned fuel rods were chosen to validate fuel burnup calculations. Two methods of measuring reactivity worth ...of fuel rods are used, traditional method is compared to newly introduced method using fuel rods swapping. Connection between both methods is described theoretically and the theory is validated experimentally. Fuel rod worth calculated using the newly introduced fuel rod swap method was within 1σ of worth measured using the traditional method. In addition to the recently performed experiments, weekly measurements of reactor core reactivity throughout the operational history are used for validation. The measured data were used to validate the fuel burnup and core criticality calculations. Fuel burnup calculations are performed using three different computer codes: the deterministic TRIGLAV, the Monte Carlo Serpent-2, and the hybrid RAPID. Great agreement was observed for Serpent-2 and RAPID by simulating fuel rod worth and its burnup, indicating that the fuel burnup and criticality calculations are accurate and that reactivity changes due to small burnup differences on the order of 10 pcm can be accurately simulated. In addition it was shown using ex-core detectors and large fission chamber that detector response changes due to fuel swapping are evident for fuel rod burnup differences of 20 MWd/kg. Fuel burnup calculations were further validated on excess reactivity measurements for three mixed TRIGA cores. The calculated burnup reactivity coefficient ΔρBU using Serpent-2 and RAPID was within 1σ of the measurements, showing both codes are capable of calculating burnup for different TRIGA fuel types.
Potential changes in patient anatomy during proton radiotherapy may lead to a deviation of the delivered dose. A dose estimate can be computed through a deformable image registration (DIR) driven ...dose accumulation. The present study evaluates the accumulated dose uncertainties in a patient subject to an inadvertent breathing associated motion.
A virtual lung tumour was inserted into a pair of single participant landmark annotated computed tomography images depicting opposite breathing phases, with the deep inspiration breath-hold the planning reference and the exhale the off-reference geometry. A novel Monte Carlo N-Particle, Version 6 (MCNP6) dose engine was developed, validated and used in treatment plan optimization. Three DIR methods were compared and used to transfer the exhale simulated dose to the reference geometry. Dose conformity and homogeneity measures from International Committee on Radioactivity Units and Measurements (ICRU) reports 78 and 83 were evaluated on simulated dose distributions registered with different DIR algorithms.
The MCNP6 dose engine handled patient-like geometries in reasonable dose calculation times. All registration methods were able to align image associated landmarks to distances, comparable to voxel sizes. A moderate deterioration of ICRU measures was encountered in comparing doses in on and off-reference anatomy. There were statistically significant DIR driven differences in ICRU measures, particularly a 10% difference in the relative D
for planning tumour volume and in the 3 mm/3% gamma passing rate.
T he dose accumulation over two anatomies resulted in a DIR driven uncertainty, important in reporting the associated ICRU measures for quality assurance.
The Real-time Analysis for Particle-transport and In-situ Detection (RAPID) Code System, developed based on the Multi-stage Response-function Transport (MRT) methodology, enables real-time simulation ...of nuclear systems such as reactor cores, spent nuclear fuel pools and casks, and sub-critical facilities. This paper presents the application of a novel fission matrix-based burnup methodology to the well-characterized JSI TRIGA Mark II research reactor. This methodology allows for calculation of nuclear fuel depletion by combination and interpolation of RAPID's burnup dependent fission matrix (FM) coefficients to take into account core changes due to burnup. The methodology is compared to experimentally validated Serpent-2 Monte Carlo depletion calculations. The results show that the burnup methodology for RAPID (bRAPID) implemented into RAPID is capable of accurately calculating the keff burnup changes of the reactor core as the average discrepancies throughout the whole burnup interval are 37 pcm. Furthermore, capability of accurately describing 3D fission source distribution changes with burnup is demonstrated by having less than 1% relative discrepancies compared to Serpent-2. Good agreement is observed for axially and pin-wise dependent fuel burnup and nuclear fuel nuclide composition as a function of burnup. It is demonstrated that bRAPID accurately describes burnup in areas with high gradients of neutron flux (e.g. vicinity of control rods). Observed discrepancies for some isotopes are explained by analyzing the neutron spectrum. This paper presents a powerful depletion calculation tool that is capable of characterization of spent nuclear fuel on the fly while the reactor is in operation.
During an earthquake on December 29th 2020, the Krško NPP automatically shutdown due to the trigger of the negative neutron flux rate signal on the power range nuclear instrumentation. From the time ...course of the detector signal, it can be concluded that the fluctuation in the detector signal may have been caused by the mechanical movement of the ex-core neutron detectors or the pressure vessel components rather than the actual change in reactor power. The objective of the analysis was to evaluate the sensitivity of the neutron flux at the ex-core detector position, if the detector is moved in the radial or axial direction. In addition, the effect of the core barrel movement and core inside the baffle movement in the radial direction were analysed. The analysis is complemented by the calculation of the thermal and total neutron flux gradient in radial, axial and azimuthal directions. The Monte Carlo particle transport code MCNP was used to study the changes in the response of the ex-core detector for the above-mentioned scenarios. Power and intermediate-range detectors were analysed separately, because they are designed differently, positioned at different locations, and have different response characteristics. It was found that the movement of the power range ex-core detector has a negligible effect on the value of the thermal neutron flux in the active part of the detector. However, the radial movement of the intermediate-range detector by 5 cm results in 7%–8% change in the thermal neutron flux in the active part of the intermediate-range detector. The analysis continued with an evaluation of the effects of moving the entire core barrel on the ex-core detector response. It was estimated that the 2 mm core barrel radial oscillation results in ∼4% deviation in the power and intermediate-range detector signal. The movement of the reactor core inside baffle can contribute ∼6% deviation in the ex-core neutron detector signal. The analysis showed that the mechanical movement of ex-core neutron detectors cannot explain the fluctuations in the ex-core detector signal. However, combined core barrel and reactor core inside baffle oscillations could be a probable reason for the observed fluctuations in the ex-core detector signal during an earthquake.
An evaluation of validity of theoretical pulse models (Fuchs-Hansen and Nordheim-Fuchs model) in reactor pulse modelling was performed by developing so called Improved Pulse Model. The effect of each ...of the five assumptions on the most important pulse physical parameters, maximal power, total released energy and full width at half maximum was studied. In the Improved Pulse Model the assumptions are disposed out with the improvements, where to account the delayed neutrons the six point kinetic equations are solved, the temperature dependences of the temperature reactivity coefficient of fuel and specific heat are taken into account, also the final ejected time of transient control rod from reactor core, whose value of reactivity varies in height and the heat dissipation from the fuel are considered. It is found that the theoretical models predict a higher maximum power, lower total released energy and full width at half maximum than the Improved Pulse Model.
A closed-water activation loop is being built at the Jožef Stefan Institute TRIGA reactor in Slovenia to serve as a well-defined and stable source of high-energy gamma rays and neutrons. The radial ...piercing port, which penetrates the graphite reflector and touching the reactor core was chosen for the installation of the closed-water loop due to the high neutron flux and favourable shielding conditions of the surrounding concrete bioshield. The main objective of this work is to assess the neutron and gamma dose field outside the port to obtain important details for the final design of the inner part of the irradiation facility and to assess the background noise for the detectors. The main part of the work consists of the design of the shielding plugs and the construction of a detailed MCNP model of the whole irradiation facility. The dose field calculations were performed with a two-step hybrid transport approach using ADVANTG for variance reduction and MCNP for particle transport. Such deep penetration and shielding calculations are challenging and computationally intensive. The results showed that the dose rate using shielding plugs is more than 7 orders of magnitude lower compared to an empty open port. To reduce the computational uncertainty, further optimisation of the ADVANTG input is essential. The final design of the shielding plugs is described. Additional lead shielding blocks will be added outside the port afterwards if needed.
We report on influence of neutron radiation on majority and minority carrier traps in n-type 4H-SiC. Together with the increase of the well-known carbon vacancy (VC) majority carrier related trap, ...neutron irradiation has introduced two deep traps, labeled as EH1 and EH3 with the activation energies for electron emission estimated as 0.4 and 0.7 eV bellow the conduction band, respectively. Based on Laplace deep level transient spectroscopy (DLTS) results, we have assigned EH1 trap to silicon vacancy (VSi). Two minority carrier traps labelled as B and D-center were detected by minority transient spectroscopy (MCTS) and assigned to substitutional boron BSi and BC, respectively. Activation energies for hole emission for B and D-center are estimated as 0.27 and 0.60 eV above the valence band, respectively. We have identified two emission lines for D-center by Laplace-MCTS measurements and assigned them to BC sitting at hexagonal (−h) and cubic (−k) lattice sites.