Function Of Burnup Research Articles

Potential of mixed oxide fuel for higher fissile conversion ratio in light water reactors

The prospect of adopting various mixed oxide fuels to increase the conversion ratio (CR) of conventional LWRs has been explored in this study. The Westinghouse AP1000 fuel assembly (FA) was taken as the reference design. While the geometry of the FA was kept the same, five assembly models were constructed utilizing different MOX fuels in the fuel rods. Depleted UO2 was mixed with PuO2 for the first model, the next two models featured a combination of UO2 with ThO2, the fourth one utilized a mixture of PuO2 and ThO2, and the final model employed a combination of UO2, PuO2, and ThO2. Depletion analysis of these fuel assembly models was then carried out using the lattice physics code Dragon 5. The burnup-dependent infinite multiplication factor and reactivity were determined and a significant improvement in cycle length was observed at the end of the burnup period because of fertile to fissile conversion. A higher fast fission rate was obtained for plutonium-based fuel models compared to the reference, whereas the thorium-based fuel models demonstrated higher thermal fission rates. The linear reactivity model (LRM) was utilized to calculate fuel cycle performance for a three-batch refueling scheme. The three fuel models utilizing PuO2 achieved improved discharge burnup and cycle length, and required lower charge fuel mass to obtain the same thermal power. All of the five proposed models attained higher conversion ratios compared to the reference, while the highest increment in CR was obtained for the last model utilizing all three oxides. The CR of conventional UO2 fuel in AP1000 was 0.64, whereas it was calculated to be 0.73 for Model 5. Key safety parameters including the fuel and moderator temperature coefficients of reactivity and the effective delayed neutron fraction of the models were also determined. Finally, the isotopic evolution of minor actinides and fission product poisons of the models were analyzed as a function of burnup and compared with the reference.

Neutronic characteristic analysis of the GAMA microreactor

The GAMA Microreactor is a low-power nuclear power plant dedicated for remote area with a power output of 300 kWe. The primary characteristics of the GAMA Microreactor are its system compactness, design simplicity, and safe operation without moving parts. The fuel is uranium hydride powder, which works simultaneously as a moderator, contained within a stainless-steel vessel and surrounded by a graphite reflector. Liquid sodium-filled heat pipes are embedded in the fuel powder to extract fission-generated heat via passive coolant flow due to capillary action. As a part of the design process, neutronic parameters of the GAMA Microreactor have been determined using SCALE6.2 and OpenMC codes. Among the analyzed parameters are the excess reactivity as a function of burn-up, shutdown margin, control rod worth, temperature coefficient of reactivity (TCR), fuel density coefficient of reactivity (DCR), and reactivity coefficient due to hydrogen dissociation. From the calculation results, the GAMA Microreactor has a maximum excess reactivity of 0.049 at the operational temperature (700 °C) and 0.0788 at ambient temperature. The TCR value is negative (−2.267 pcm/°C) at the beginning of the cycle (BOL) and −4.9618 pcm/°C at the end of the cycle (EOC). Hydrogen dissociation imposes negative reactivity on the reactor both at the BOL (−0.2631 Δk/k-%dissociation) and EOC (−0.3524 Δk/k-%dissociation). Therefore, the GAMA Microreactor maintains its inherent safety characteristics over a considerably long operation time.

TRISO Burnup-Dependent Failure Analysis of a HTGR Design-Basis Accident Using BISON

This work assesses the failure behavior of the silicon carbide (SiC) layer in TRIstructural ISOtropic (TRISO) fuel particles with BISON during both steady-state and transient conditions for the MHTGR-350 design by General Atomics. A one-dimensional BISON model for a uranium oxycarbide–bearing TRISO is developed to simulate a power level of 50 mW/particles up to 12.4% fissions per initial metal atom (FIMA). Stress predictions for the materials of interest are presented as a function of burnup, along with the SiC failure probability computed using Weibull statistics under the assumption of realistic SiC quality. A design-basis accident (DBA) that involves control rod withdrawal is then simulated in BISON at various burnup levels. Boundary conditions during the DBA are obtained from RELAP5-3D. The predicted SiC failure probability is 3 × 10−14%. This value does not increase during the transient as a function of burnup since irradiation induces a compressive hoop stress state that reaches a maximum absolute value of around 300 MPa. This compressive stress compensates the tensile loading conditions introduced by the transient. The tensile components appear amplified at higher burnups since they increase from 100 to 140 MPa going from 6% FIMA to 12.4% FIMA. Nonetheless, burnup provides a negligible impact on the overall failure probability predictions for SiC, as transient conditions do not translate to a stress increase toward tensile values at any of the considered burnups. Additionally, a two-dimensional (2-D) exploratory BISON model is developed to determine the effects of inner pyrolytic carbon (IPyC) crack formation on SiC. IPyC failure is predicted using Weibull statistics at 105 days from the beginning of irradiation, which is set as the time for crack initiation using XFEM in the 2-D BISON model. Stress concentrations induced by the crack cause the SiC failure probability to increase with respect to the intact particle case by approximately 10 orders of magnitude. Maximum stress concentrations are found at the time of crack formation, after which stress relaxation is observed during remaining steady-state irradiations and subsequent transient simulations. Considerations are made on the results of the 2-D model being contingent on both mechanical and anisotropy properties for pyrolytic carbon. Recommendations are provided for future work involving higher dimensionality models, further investigations on uncertainties, material properties, and additional designs, such as the fluoride salt-cooled high-temperature reactor, that operate at higher burnups.

Verification of a novel fuel burnup algorithm in the RAPID code system based on Serpent-2 simulation of the TRIGA Mark II research reactor

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.

Nuclear data uncertainty quantification on PWR spent nuclear fuel as a function of burnup

Nuclear data uncertainty analysis on the spent nuclear fuel inventory was performed on the Takahama-3 NT3G23 assembly, where the sample SF95-4 was irradiated up to a burnup of approximately 36GWdtaccording to the SFCOMPO benchmark. The cross-section covariance matrices stored in the ENDF/B-VIII.0, JEFF-3.3 and JENDL-4.0u evaluated nuclear data libraries were propagated with the stochastic sampling algorithms implemented in the SANDY code. A comparison of the concentration uncertainty differences obtained using data from the three libraries is reported. Similarities were found with the fuel composition uncertainty results obtained for the Calvert Cliffs MKP109 sample P SFCOMPO benchmark. Such a similarity was also found when comparing concentration uncertainties along the sample irradiation. Therefore, the main contributors to the concentration uncertainty of a number of nuclides were identified at different burnup levels in the two samples. To complement the similarity analysis, a correlation study of the concentration distributions predicted by the two models was performed. The reported results hint a dominance of the common uncertainty propagation mechanisms over the model differences in the determination of concentration uncertainty.

Study on the benefit of HLLW partitioning on the high-level waste glasses from the viewpoint of waste management

The volume and characteristics of heat-generating waste glasses produced through spent fuel reprocessing PUREX process + HLLW partitioning TRPO process are evaluated as a function of burn-up and cooling period of spent fuel in this study. The waste loading of glass is assumed to be restricted by the heat generation rate, MoO3 content, and noble metal content. The results show that: 1) For spent fuel of 8–20 years cooling, the volume of waste glasses after HLLW partitioning is 39%–25% of that before partitioning. 2) After removing >99.967% of α-nuclides in HLLW through TRPO process, the TRPO raffinate stream can be immobilized into waste glasses. After 95 years of decay storage, intermediate-depth disposal rather than geological disposal of TRPO raffinate waste glasses could be engineered. 3) Separation of Sr and Cs from TRPO raffinate hardly reduces the volume of waste glasses. SrCs separation before vitrification is unnecessary. 4) The insoluble residues generated during the shearing and dissolution of spent fuel may not be suitable for mixing with the TRPO raffinate stream, nor for immobilizing into a borosilicate glass matrix. Alternative solutions rather than vitrification are encouraged.

Multiphysics simulation of VVER-1200 fuel performance during normal operating conditions

Nuclear fuel performance modeling and simulation are critical tasks for nuclear fuel design optimization and safety analysis under normal and transient conditions. Fuel performance is a complicated phenomenon that involves thermal, mechanical, and irradiation mechanisms and requires special multiphysics modules. In this study, a fuel performance model was developed using the COMSOL Multiphysics platform. The modeling was performed for a 2D axis-symmetric geometry of a UO2 fuel pellet in the E110 clad for VVER-1200 fuel. The modeling considers all relevant phenomena, including heat generation and conduction, gap heat transfer, elastic strain, mechanical contact, thermal expansion, grain growth, densification, fission gas generation and release, fission product swelling, gap/plenum pressure, and cladding thermal and irradiation creep. The model was validated using a code-to-code evaluation of the fuel pellet centerline and surface temperatures in the case of constant power, in addition to validation of fission gas release (FGR) predictions. This prediction proved that the model could perform according to previously published VVER nuclear fuel performance parameters. A sensitivity study was also conducted to assess the effects of uncertainty on some of the model parameters. The model was then used to predict the VVER-1200 fuel performance parameters as a function of burnup, including the temperature profiles, gap width, fission gas release, and plenum pressure. A compilation of related material and thermomechanical models was conducted and included in the modeling to allow the user to investigate different material/performance models. Although the model was developed for normal operating conditions, it can be modified to include off-normal operating conditions.

Model for determining rupture area in Zircaloy cladding under LOCA conditions

The dispersal of nuclear fuel through a cladding rupture opening during a transient event is the primary safety concern preventing the United States nuclear industry from extending burnup beyond a peak rod average burnup of 62 GWd/tU. The majority of research thus far has prioritized understanding, or at a minimum prioritized quantifying, fuel fragmentation and pulverization as a function of burnup and to a lesser extent onset temperature. Early high burnup loss-of-coolant accident tests performed at the Halden Boiling Water Reactor evaluated the impact of relocation on cladding temperature, whereas historical and recent cladding burst tests focused primarily on generating data related to cladding balloon and burst timing. However, very little effort has been put forth to develop a relationship for determining the rupture opening dimensions for Zircaloy cladding tubes bursting during a simulated light water reactor loss-of-coolant accident. This is significant because high burnup loss-of-coolant accident data suggest dispersal could be limited by the size of the rupture opening or rupture dimensions, thereby minimizing the safety concern. Therefore, the purpose of this research is to build a database composed of all simulated loss-of-coolant accident data to develop an empirical model that conservatively captures the rupture opening dimensions and translates those dimensions to rupture area. Analysis of the experimental data indicates the rupture dimensions are strongly correlated to the peak balloon strain. This relationship is leveraged to develop a rupture length empirical correlation as a function of peak balloon strain and a rupture width empirical correlation as a function of rupture length. Outliers in the experimental data are identified and discussed. Finally, the rupture opening empirical model is presented in such a way for general implementation into high fidelity fuel performance codes or fuel performance codes used for safety analyses.

Inter-Code Comparison of Computational VERA Depletion Benchmark Using OpenMC, OpenMC-ONIX and DRAGON

This research focuses on the comparative analysis of the PWR fuel assembly based on VERA depletion benchmark problems using community-developed open source Monte Carlo code OpenMC, python based burnup code system ONIX (a coupling interface for Monte Carlo code OpenMC), and deterministic DRAGON code. The depletion analysis was performed using OpenMC and ONIX with ENDF/B-VII.1 nuclear data library, and DRAGON with SHEM-361 based DRAGLIB format library (ENDF/B-VII.1). The code-to-code analysis on the evolution of k eff , atom number density, and power distribution as a function of burnup has been performed and the result shows a good agreement with the maximum difference within 200 pcm at EOC. However small discrepancy around 90 pcm has been observed in k eff calculated by DRAGON compared to OpenMC in the presence of integral fuel burnable absorbers (IFBA). The above-mentioned codes have been validated successfully for the first time against PWR fuel assembly based on VERA depletion benchmark problems. It can be concluded that initial implementation of these codes at the Department of Nuclear Science and Engineering under Military Institute of Science and Technology, Dhaka, was successful and that further research works are to be performed to utilize these codes for depletion/neutronics calculation of existing 3MW TRIGA Mark-II research reactor and VVER-type power reactor that is to be commissioned in Bangladesh.

ANALYSIS OF THORIUM PIN CELL BURN UP OF THE PWR USING WIMS CODE

A thorium-fueled benchmark comparison was made in this study between state-of-the-art codes, WIMSD-5B code to MOCUP (MCNP4B + ORIGEN2), and CASMO-4 for burnup calculations as part of efforts to examine the possible benefits of using thorium in PWR fuel. WIMSD-5B calculations employ the same model as a reference, MOCUP, and CASMO, however, there are some variances in methodology and cross-section libraries. On a PWR pin cell model, eigenvalue and isotope concentrations were examined up to high burnup. The eigenvalue comparison as a function of burnup is good, with a maximum difference of less than 5% and an average absolute difference of less than 1%. The isotope concentration comparisons outperform a set of ThO2-UO2 fuel benchmarks and are comparable to a set of uranium fuel benchmarks previously published in the literature. As a function, the eigenvalue comparison The actinide and fission product data sources for a typical thorium fuel are reported in the WIMSD-5B burnup calculations. The reasons for discrepancies in coding are examined and explored.Keywords: Thorium, PWR Fuel, Burn up, Pin Cell, WIMSD-5B

Neutronic Analysis of SiC/SiC Sandwich Cladding Design in APR-1400 under Normal Operation Conditions

Our aim is to study the neutronic behaviour of potential accident-tolerant fuel (ATF) claddings in a pressurised water reactor under normal operations. This work compares ATF silicon carbide composite (SiC/SiC) cladding to conventional ZIRLOTM cladding in APR-1400. Additionally, a “sandwich” cladding design developed by the CEA is used for SiC/SiC. The design structure includes a liner in between two layers of the composite to ensure leak tightness. The two proposed liners are Niobium (Nb) and Tantalum (Ta). Serpent 2, a Monte Carlo reactor physics lattice code, is employed to model both cladding materials in APR-1400 at three different levels: pin cell, fuel assembly, and core. The criticality, neutron spectrum, actinide inventory, and power distribution as a function of burnup are investigated. The simulations show that SiC/SiC with the Nb liner displays a far superior performance than the Ta liner across all examined characteristics. Ta leads to a harder neutron spectrum and increased Pu-239 content throughout the cycle, while Nb presents negligible effects. In fact, SiC/SiC with the Nb liner performs very similarly to ZIRLOTM at all model levels. The results indicate that, in terms of neutronics, the adoption of the SiC/SiC composite would entail little to no changes to current APR-1400 operations.

Neutronic Performance of a Compact 10 MWe Nuclear Reactor with Low Enrichment (ThxU1\u2212x)N Fuel

Recent interest in compact nuclear reactors for applications in space or in remote locations drives innovation in nuclear fuel design, especially non-oxide ceramic nuclear fuels. This work details neutronic modeling designed to support the development of a new nuclear fuel concept based on a mixture of thorium and uranium nitride. A Monte Carlo N-Particle Version 6.2 (MCNP-6) model of a compact 10 MWe reactor design which incorporates (ThxU1−x)N fuel is presented. In this context, a “compact” reactor is a completely assembled reactor which may be emptied of coolant and transported by specialized commercial vehicle, deployed by a C130J aircraft, or launched into space. Core geometry, reflector barrels, and the heat exchange zones are designed to support reduction of overall reactor volume of core components while maintaining criticality with a fixed total fuel mass of 4500 kg. Dense mixed nitrides of thorium nitride (ThN) additions in uranium nitride (UN) in 5 wt.% increments between 0.05 le x le 0.5 have been considered for calculation of k_{infty } and k_{{{text{effective}}}}. ThN additions in UN results in a slight increase in the magnitude of the temperature coefficient of reactivity, which is negative by design. The isotopic distribution of the principal actinide inventory as a function of burnup, time, and initial fuel composition is presented and discussed within the context of the proliferation risk of this core design.

Calculation of PWR Thorium Pin Cell Burnup and Isotope Prediction Using WIMSD-5B Code

The amount of thorium in the world is very abundant and the benefits of using thorium as fuel for nuclear reactors such as PWR are increasingly widespread, so it needs to be analyzed. In the initial phase, the PWR pin cell fuel analysis was carried out using the WIMSD-5B code. The PWR ThO2-UO2 fuel in the cladding with the fuel assembly model 17 x 17 is analyzed by its fuel burn up, neutron flux and isotope composition after maximum operation. The pin cell model consisting of fuel, cladding and coolant with a square pitch of 12.6 cm are calculated using different libraries. Eigenvalues, neutron fluxes and isotope concentrations were compared with the PWR pin cell model to high burnup. The eigenvalues and flux neutron as a function of burnup are good, the maximum difference is within 7.41% and the mean absolute difference is less than 3.12 %. The best comparison of fission product isotope concentrations is to use the ENDFB7.1 library with Thorium fuel and is comparable to the uranium fuel reported in the literature. Data sources for actinides and fission products used infuel depletion calculations for thorium fuel are all documented.

Burnup computations of multi-pass fuel loading scenarios in HTR-10 using a pre-generated fuel composition library

This study aims to precisely investigate the core characteristics of HTR-10 using an online-refueling fuel loading scenario through an innovative approach that incorporates a pre-generated fuel composition library. HTR-10 is an experimental high-temperaturegas- cooled pebble-bed reactor constructed in China. Applying multi-pass fuel loading scenario, the pebbles in HTR-10 pass through the core several times via onlinerefueling. In this study, we proposed a model that simulates the dynamic operation of HTR-10 by using continuously static calculations. The model was constructed by dividing the core into 3 flow channels with different amounts of horizontal zones according to the residence time of fuel pebbles. However, random distribution of the pebbles with various burnup values complicates the refueling process. In order to simulate online-refueling more conveniently, a comprehensive method based on multipass fuel loading scenarios was thus developed. Firstly, a pre-generated fuel composition library containing the content of both actinides and fission products of fuel pebbles as a function of burnup, was built by implementing full core depletion calculations. During the refueling processes, a specific ratio of discharged fuel pebbles to fresh fuel pebbles was reloaded into the top zones. Thus, the average burnup values in the top zones can be resulted by considering the supplement of fresh fuel pebbles. Eventually, the fuel composition in the top zones were obtained by interpolating the average burnup value from the pre-generated library. In this study, all the calculations were performed using the MCNP6 computer code together with the ENDF/B-VII continuous energy neutron data library. The results revealed that the core reached a state of equilibrium after seven refueling processes with an increased rate of volume fraction. Our findings correspond to those of other studies. Another model constructed by locating the actual distribution of fresh pebbles also confirms the verification of the refueling method. Eventually, the burnup characteristics as well as the power and flux distributions in the HTR-10 core were analyzed using the simulation method.

Simulation of a NuScale core design with the CASL VERA code

Three-dimensional (3D) full-core calculations are an integral part of fuel reload design for today’s light water reactors (LWRs). The current approaches are typically based on the nodal diffusion codes that calculate criticality state points and power distributions as a function of burnup. The CASL VERA (Consortium for Advanced Simulation of Light Water Reactors, Virtual Environment for Reactor Applications) code represents one of the latest advancements in 3D full-core calculation analysis based on transport theory methods employing the Method of Characteristics (MOC) and coupled multi-physics. In this article, a publicly released version of the NuScale reactor core is analysed with the VERA code (version 3.9 and 4.1) and contrasted against some static Serpent and Polaris based simulations. The analysis shows an excellent agreement for the lattice-level calculations as well as with some of the 3D full-core models. However, larger deviations were found in cases with heavy reflector models, whereby the reflector composition was found to impact the differences between the VERA and Serpent results. In this analysis, it was determined that greater than 90% stainless steel content in the heavy reflector leads to higher deviations between the VERA results and the Serpent results. The burnup calculations showed that the presence of the heavy reflector extends the cycle length and also leads to a flatter power distribution in the core, which can generally be interpreted as more efficient.

Effect of Axial Distribution of Gadolinium Burnable Poison in Advanced Pressurized Water Reactor Assembly

The radial and axial power distribution in power reactors are determined mainly by the patterns of the fuel assembly and the burnable absorber at the beginning of cycle. In Advanced Pressurized Water Reactor (APWR), gadolinium burnable absorber is used to decrease the relative power of fresh fuel assemblies. In this paper, the effect of the axial distribution of gadolinium (Gd) on the power of the APWR assembly is studied. Three models of APWR assemblies are simulated using MCNP6 code. In the first model, UO2 fuel is distributed uniformly in all the fuel rods. In the other two models some of the UO2 fuel rods are replaced by UO2-Gd2O3 rods in part length distribution. Two gadolinium concentrations 6% and 10% are used. The main neutronic parameters are estimated for the three models: the multiplication factor (K-infinity) as a function of burnup (GWd/MTU), the radial and axial power distributions. The results show that the distribution of the gadolinium absorber in the central region of fuel rod (part-length absorber) leads to flattening of axial power, which means additional axial power distribution control.

Whole core analysis of the single-fluid double-zone thorium molten salt reactor (SD-TMSR)

The SD-TMSR (2,250 MWth) is a Single-fluid Double-zone Thorium-based Molten Salt Reactor. The active core of the SD-TMSR is divided into the inner zone (486 fuel tubes) and the outer zone (522 fuel tubes) to improve the Th-U3 breeding performance. This work adopted the SERPENT-2 Monte Carlo code to analyze the whole core model of the SD-TMSR. Built-in SERPENT-2 capabilities simulated online reprocessing and refueling and calculated the multiplication factor and Breeding Ratio (BR). We found that the molten salt Temperature Coefficient of Reactivity (αT) was negative for initial and equilibrium states. This study investigated the variation of the multiplication factor, BR, and build-up of important nuclides in the core as a function of burnup. Under online reprocessing and refueling, we studied the variation of the reactivity during 60 years of reactor operation. Additionally, the neutron spectrum shift during the reactor operation was calculated. Finally, these simulations determined the appropriate 232Th and 233U feed rates for maintaining criticality and enabled analysis of the overall SD-TMSR fuel cycle performance.

Cross-Section Generation Using TXT2NTAB Code for Uncertainty Propagation with Burnup Dependence

One of the challenges of studying the neutronics of reactors is to generate reliable parameterized libraries that contain information to simulate the core in all possible operational and transient conditions. These libraries must include tables of cross sections and other neutronic and kinetic parameters and are obtained by simulating all the segments in a transport code. At the lattice level, one can use branch calculations to change “instantaneously” the feedback parameters as a function of burnup. When using random sampling for the lattice calculations, one can obtain statistical information about the output parameters and use it in a core simulation to characterize the accuracy of data estimating uncertainties when simulating a heterogeneous system at different scales of detail. This work presents the methodology to generate NEMTAB libraries from data obtained in the SCALE code system to be used in PARCS simulations. The code TXT2NTAB is used to reorder the cross-section tables in NEMTAB format and generate another NEMTAB of standard deviation. With these libraries, the authors perform a steady-state calculation for a light water reactor to propagate several uncertainties at the core level. The methodology allows obtaining statistical information of the most important output parameters: multiplication factor keff, axial power peak Pz, and axial peak node Nz.

Fuel burnup analysis and fuel management for TRIGA Mark III research reactor

The fundamental advantage in using Monte Carlo methods for burnup calculations is to formulate an effective optimal fuel management strategy for the TRR-1/M1 research reactor. The core management study has been performed by utilizing the essentially parameters including multiplication factor, power peaking, neutron flux and burnup calculation based on the Monte Carlo calculation. The fuel element burnup was calculated after reshuffling the reactor core. The fuel cycle length and core parameters such as core excess reactivity, neutron flux, axial and radial power factors and other parameters are determined. The core excess reactivity was calculated as a function of burnup. The maximum excess reactivity shall not exceed 6.3% Δk/k. The maximum fuel temperature shall not exceed 930 °C during steady-state operation. Typically, a core loading operated with the maximum burnup between 100 to 200 MWD depending on the utilization requirements. The thermal neutron flux in the irradiation positions is within the order of 1011 - 1013 n/cm2-sec. The study gives valuable results into the behaviour of the TRR-1/M1 research reactor and will ensure optimized utilization and operation of the reactor during its life time. It will establish the strategic planning for fuel management in the reshuffling and reloading schemes patterns and its safe implementation in the future.

Modeling and evaluation of fissile material utilization of the UMLRR using Monte Carlo MCNP6 code

Accurate knowledge of the neutron flux and temporal nuclide inventory is essential for reactor operations to meet the requirements of criticality safety, safeguards, and spent fuel storage. The University of Massachusetts Lowell Research Reactor (UMLRR) reactor has been modeled with the deterministic-based code, VENTURE and with an older version of MCNP (MCNP5). The Monte Carlo N-Particle (MCNP6) code with integrated buildup depletion code (CINDER90) provides a high-fidelity tool that can be used to perform 3D, full core simulation to evaluate fissile material utilization, and nuclide inventory calculations as a function of burnup. This paper provides a comprehensive analysis of the UMLRR core, including the development of a 3D model of core assemblies for criticality calculation, the analysis of the axial and radial neutron flux profiles and the depletion evaluation of the UMLRR fuel. The benchmarked analysis of the UMLRR core using MCNP6 will be used to track core performance, for control blade worth evaluation, the design of experiments, fuel lifetime and optimize fuel utilization.