Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fissile”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Discriminating Uranium Isotopes Based on Fission Signatures Induced by Delayed Neutrons

The use of active interrogation to induce delayed neutron emission is a well-established technique for the characterization of special nuclear materials. Delayed neutrons have isotope-characteristic spectral and temporal signatures, which provide the basis for isotope identification. However, in bulk materials that contain an appreciable fissile (e.g., 235 U or 233 U) fraction, such as highly enriched uranium (HEU), delayed neutrons have a high probability of inducing additional fissions. As a result, the overall delayed neutron signature consists of two distinct components: the “primary” delayed neutrons (emitted directly by fission fragments) and the “secondary prompt” fission neutrons produced in fission induced by primary delayed neutrons. These prompt products differ from “primary” delayed neutrons both in their energy spectra and in the presence of coincident radiation released by the parent fission event. The presence and relative quantity of prompt products from delayed fission depend on the cross section of the material in the energy range of delayed neutrons, which may differ significantly between isotopes, thus providing an exploitable means for isotope differentiation. Here, we demonstrate two experimental approaches for discriminating between 235 U and 238 U isotopes based on the measurement of delayed neutron-induced fission products. First, HEU and depleted uranium objects are differentiated through the detection of high-energy prompt neutrons from delayed fission using both recoil-based organic liquid scintillators and thermalization spectra from a custom-built capture-gated composite detector. Secondly, coincident radiation measurements are used as the basis for discrimination by comparing the overall rates and time evolution of fission events when delayed neutrons are present.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2023

Passive and active neutron signatures of 233 U for nondestructive assay

The thorium fuel cycle is emerging as an attractive alternative to conventional nuclear fuel cycles, as it does not require the enrichment of uranium for long-term sustainability. The operating principle of this fuel cycle is the irradiation of 232 Th to produce 233 U, which is fissile and sustains the fission chain reaction. 233 U poses unique challenges for nuclear safeguards, as it is associated with a uniquely extreme γ-ray environment from 232 U contamination, which limits the feasibility of the γ-ray-based assay, as well as more conservative accountability requirements than for 235 U set by the International Atomic Energy Agency. Consequently, instrumentation used for safeguarding 235 U in traditional fuel cycles may be inapplicable. It is essential that the nondestructive signatures of 233 U be characterized so that nuclear safeguards can be applied to thorium fuel-cycle facilities as they come online. In this work, a set of 233 U 3 O 8 plates, containing 984 g 233 U, was measured at the National Criticality Experiments Research Center. A high-pressure 4 He gaseous scintillation detector, which is insensitive to γ-rays, was used to perform a passive fast neutron spectral signature measurement of 233 U 3 O 8 , and was used in conjunction with a pulsed deuterium-tritium neutron generator to demonstrate the differential die-away signature of this material. Furthermore, an array of 3 He detectors was used in conjunction with the same neutron generator to measure the delayed neutron time profile of 233 U, which is unique to this nuclide. These measurements provide a benchmark for future nondestructive assay instrumentation development, and demonstrate a set of key neutron signatures to be leveraged for nuclear safeguards in the thorium fuel cycle.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Extension of the Hauser-Feshbach fission fragment decay model to multichance fission

The Hauser-Feshbach fission fragment decay model, HF 3 D, which calculates the statistical decay of fission fragments, has been expanded to include multichance fission, up to neutron incident energies of 20 MeV. The deterministic decay takes as input prescission quantities—fission probabilities and the average energy causing fission—and postscission quantities—yields in mass, charge, total kinetic energy, spin, and parity. From these fission fragment initial conditions, the full decay is followed through both prompt and delayed particle emissions, allowing for the calculation of prompt neutron and γ properties, such as multiplicity and energy distributions, both independent and cumulative fission yields, and delayed neutron observables. Here, we describe the implementation of multichance fission in the HF 3 D model, and show an example of prompt and delayed quantities beyond first-chance fission, using the example of neutron-induced fission on 235 U. This expansion represents significant progress in consistently modeling the emission of prompt and delayed particles from fissile systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Nuclear fission reaction simulations in compact stars

Type-Ia supernovas are powerful stellar explosions that provide important distance indicators in cosmology. Recently, we proposed a new Type-Ia supernova mechanism that involves a nuclear fission chain-reaction in an isolated white dwarf [Phys. Rev. Lett. 126, 131101 (2021)]. Here, in this work, we perform novel reaction network simulations of the actinide-rich first solids in a cooling white dwarf. The network includes neutron-capture and fission reactions on a range of U and Th isotopes with various possible values for 235 U enrichment. We find, for modest 235 U enrichments, neutron capture on 238 U and 232 Th can breed additional fissile nuclei so that a significant fraction of all U and Th nuclei may fission during the chain reaction. Finally, we compute the energy release from the fission chain reaction for various uranium enrichments; a novel result that is a necessary input for thermal diffusion simulations of carbon ignition.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Nuclear rotation at the fission limit in Rf 254

A ground-state rotational band in the fissile nucleus Rf254 was observed for the first time. Levels up to spin 14 and excitation energy of 1.56 MeV were observed. The Rf254 nuclei were produced using the Pb206(Ti50,2n) fusion-evaporation reaction. It is the weakest reaction channel ever studied using in-beam γ-ray spectroscopic methods. The reaction products were separated from the beam in the Argonne gas-filled analyzer (AGFA). The Rf254 nuclei were implanted into a double-sided Si strip detector at the AGFA focal plane and tagged with subsequent ground-state spontaneous fission decays using temporal and spatial correlations. Prompt γ rays in coincidence with the Rf254 recoils were detected in the Gammasphere array of Ge detectors. In order to identify the ground-state rotational band in Rf254, a method for identifying rotational bands in low statistics γ-ray spectra was developed. The deduced Rf254 kinematic moment of inertia is smaller compared to neighboring even-even nuclei. This is most likely associated with a slightly lower quadrupole deformation and stronger pairing correlations in Rf254. The behavior of the moment of inertia as a function of rotational frequency is similar to that of the lighter N=150 isotones Fm250 and No252.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Constraining calculations of Pu 239 ( n , γ ) cross sections with multistep cascade γ -ray spectra measured at the DANCE facility

Measurements of multistep cascade gamma-ray spectra from several low-lying 1 + resonances in the Pu 239 ( n , γ ) reaction were made using the DANCE detector. The measured spectra are compared to spectra simulated using several models for the photon strength function and the nuclear level density, representing the first analysis of this nature on a fissile nucleus. In general, models of the photon strength function that included a strong M 1 “scissors mode” contribution coupled with a Gilbert-Cameron nuclear level density provided a reasonable description of the spectra. Simulations using photon strength functions and nuclear level densities specifically for Pu 240 , extracted by the Oslo method, resulted in a very good representation of the measured spectra. Neutron-capture cross section calculations that included the scissors mode provided a good estimate of the previously measured values, without additional normalization. However, the scissors-mode energy that produced the best description of the multistep cascade spectra was different from the value that best reproduced the capture cross section. In addition, the 〈 Γ γ 〉 inferred from the multistep cascade calculations was somewhat less than the tabulated value, which was determined from resonance analysis. Published by the American Physical Society 2025

Ullmann, J. L. (ORCID:0000000207784405)↗
DOE OSTI2025

Comprehensive Measurement of the Reactor Antineutrino Spectrum and Flux at Daya Bay

This Letter reports the precise measurement of the reactor antineutrino spectrum and flux based on the full dataset of 4.7 ×10 6 inverse-beta-decay candidates collected at Daya Bay near detectors. Expressed in terms of the inverse-beta-decay yield per fission, the antineutrino spectra from all reactor fissile isotopes and the specific 235 U and 239 Pu isotopes are measured with 1.3%, 3%, and 8% uncertainties, respectively, near the 3 MeV spectrum peak in reconstructed energy, reaching the best precision in the world. The total antineutrino flux and isotopic 235 U and 239 Pu fluxes are precisely measured to be 5.84 ± 0.07, 6.16 ± 0.12, and 4.16 ± 0.21 in units of 10 −43 cm 2 /fission. These measurements are compared with the Huber-Mueller model, the reevaluated conversion model based on the Kurchatov Institute measurement, and the latest summation model (SM2023). Furthermore, the Daya Bay flux shows good consistency with the Kurchatov Institute and SM2023 models but disagrees with the Huber-Mueller model. The Daya Bay spectrum, however, disagrees with all model predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Colloquium : Neutrino detectors as tools for nuclear security

For over 40 years, physicists have discussed possible uses for neutrino detectors in nuclear nonproliferation, arms control, and fissile materials security. Neutrinos are an attractive fission signature because they readily pass through matter. The same property makes neutrinos challenging to detect in systems that would be practical for nuclear security applications. This Colloquium presents a broad overview of several potential neutrino applications, including the near-field monitoring of known reactors, far-field monitoring of known or discovery of undeclared reactors, detection of reactor waste streams, and detection of nuclear explosions. Recent detector advances have made near-field monitoring feasible, whereas farther-field reactor detection and waste stream detection monitoring may be possible in some cases with further research and development. Very long-range reactor monitoring and nuclear explosion detection do not seem feasible for the foreseeable future due to considerable physical and/or practical constraints.

42 ENGINEERING↗
DOE OSTI2023

Review—Fundamental Uranium Electrochemistry and Spectroscopy in Molten Salt Systems

Uranium is a key element used for nuclear energy production. Some advanced reactor designs, specifically molten salt reactors, will continue to use uranium as the fissile material for energy production. These new technologies require an intimate understanding of uranium chemistry during and after energy production. This review covers contemporary research on the coordination chemistry and behavior of uranium with the coolant and pyroprocessing salts as proposed for use in future reactor designs. Discussed topics include the nature of U redox reactions involving the reduction of U(III) to U metal and oxidation of U(III) to U(IV). These systems have been interrogated using cyclic voltammetry, chronopotentiometry, and optical and X-ray absorption spectroscopies. Insights obtained into the electrode potentials, the uranium species, and their diffusion coefficients in alkali halide melts from decades of research are summarized selectively. Further, perspectives are provided on the importance of unifying studies for comparison across multiple institutions. The application of synchrotron radiation research and multimodal approaches involving two (or more) probes, such as the widespread combination of UV–visible spectroscopy and electroanalysis known as spectroelectrochemistry, can provide new knowledge about the main process of uranium electrorefining—diffusion, as will be demonstrated in this review through the lack of comparable results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Agn-201 Digital Twin

This is the repository for all code related to the AGN-201 Nuclear Reactor Digital Twin at Idaho State University. The goal of this code repository is to consolidate all pieces required to run the AGN-201 Digital Twin in the [DeepLynx](https://github.com/idaholab/Deep-Lynx) ecosystem. This is the first successfully launched digital twin of a fissile nuclear reactor that we are aware of. While the code is not complex, the problems of networking, policy, and initial groundwork were significant to overcome.

Darrington, JohnW.↗
DOE OSTI2022

Space-dependent calculation of the multiplicity moments for shells

In earlier work, we extended the methodology of multiplicity counting in nuclear safeguards, by elaborating the one-speed stochastic transport theory of the calculation of the so-called multiplicity moments, i.e. the factorial moments of the number of neutrons emitted from a fissile item, following a source event from an internal neutron source (spontaneous fission and (α, n) reactions). Calculations were made for spheres and cylinders of various shapes. In all our work so far, the material of the items was homogeneous, and the distribution of the internal source was assumed to be uniformly distributed within the item, with the neutron emission assumed to be isotropic. In the present work the calculations are extended to the case of a point source inside either a solid sphere or in a spherical shell. This necessitates the extension of the theory to non-homogeneous items and non-uniform and non-isotropic sources. This work describes the extension of the theory and provides some quantitative results. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Adjoint sensitivity analysis and data assimilation for verification of dry storage cask contents

Dry cask storage is a method for interim storage of spent fuel assemblies which contain fissile isotopes of uranium and plutonium. These can present a proliferation concern and consequently there is a need for non-destructive testing methods to verify a dry cask's contents for proliferation protection. We present an application of adjoint sensitivity analysis and data assimilation to a multigroup diffusion model of dry cask storage. Adjoint sensitivity analysis allows the efficient calculation of sensitivities for use in data assimilation to calibrate imprecisely known parameter values and data consistency tests to detect diversion scenarios. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

VADER: A Tool for Criticality Safety Validation

The purpose of criticality safety is to prevent any inadvertent criticality from occurring during the handling or storage of fissile material. Calculations are frequently used to demonstrate that a sufficient subcritical margin exists. Validation is a key aspect of the evaluation process, establishing the suitability, accuracy, and associated uncertainty of the computational method and data to be used for the intended application. The validation process is performed by comparing the results of critical experiments with the calculated results from models of the experiments using the computational method to be validated. Laboratory critical experiments are controlled systems that achieve a k eff of approximately 1 in order to investigate the parameters at which such a critical condition is achieved. The validation parameters that are traditionally applied to safety analysis calculations are the bias and the bias uncertainty . The bias is the deviation of the average k eff of the validation suite from unity. The bias uncertainty accounts for the statistical uncertainty in the bias based on the standard deviation, sample size, and distribution of k eff values of the validation suite. The values of bias and bias uncertainty ensure that the systems predicted to be subcritical by the computational method will indeed be subcritical. The bias and bias uncertainty are often combined to determine an upper subcritical limit (USL) or computational margin that can then be applied to safety analysis calculations. Many methods have been developed by different organizations to calculate the bias and bias uncertainty for various types of criticality analyses. Each of these methods typically requires that the validity of various underpinning statistical assumptions be confirmed to demonstrate that the method is appropriate for the analysis of a given validation suite. An example of the validation decision making flow is shown in Fig.1. As shown in Fig. 1, the analyst performing the validation fits a trend line to the data and performs a test to determine if the trend was a statistically better representation of the data than if it were treated as an uncorrelated sample. If the trend line is a better representation of the data, then the analyst uses any one of a number of trending techniques to determine the bias and bias uncertainty. If a trend is not an appropriate representation of the data, then the analyst proceeds to perform a normality assessment for the data. If the normal assumption can be shown to be acceptable, then the analyst calculates the bias and bias uncertainty with the parametric technique. If the assumption of normality cannot be justified, then the nonparametric technique is used. Once the decision flow has been followed and the appropriate technique has been selected, the bias and bias uncertainty is typically combined with an administrative margin to determine a USL below which calculated values of k eff for safety analysis models can be considered subcritical. The calculations used in each decision are often performed with spreadsheets or with small programs available at various sites performing criticality analyses. Expertise in understanding and interpreting the results must be maintained to perform these calculations. This can often be an error-prone process. Oak Ridge National Laboratory (ORNL) is currently developing the Validation and Data Evaluation Resource (VADER) to simplify and automate the criticality safety validation process and to provide a software quality assurance pedigree to the calculational methods used. This paper discusses the use of the Fulcrum user interface with VADER, the anticipated initial capabilities of VADER to perform validation analyses, and the output from the code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Comparison of NMC estimates with trans-stilbene, EJ-309, and He-3 detection systems

Neutron multiplicity counting (NMC) is a technique for the assay of fissile material. In this work, three detection systems are utilized for active interrogation assay of shells of the Rocky Flats shells (highly enriched uranium, 93% 235U) stacked from 13.25-54.92 kg assemblies. The singles ($R_1$) and doubles ($R_2$) rates are calculated with each system to estimate two sample parameters: $M_L$- the leakage multiplication and F - the sample fission rate. The estimated mass, m, is found by dividing F by a constant activity per unit mass. Since we are interrogating HEU, the α ratio of (α; n) to fission neutrons is taken to be zero. The system of equations to calculate these quantities was originally derived. The Neutron Multiplicity 3 He Array Detector (NOMAD) consists of 15 3 He tubes inside polyethylene and is the traditional, capture-based detection system for NMC. The polyethylene moderates incoming neutrons for thermal capture in individual tubes. The low gamma background, discrete capture signals, and high efficiency of the NOMAD are beneficial for NMC. However, the time to slow down neutrons to thermal energies leads to a system die away time on the scale of microseconds. Comparatively, the Rossi-alpha Measurements – Rapid Organic (n, γ) Discrimination Detector (RAMRODD) and the Organic Scintillator Array (OSCAR) are scatter-based detection systems. RAMRODD consists of 8, 5.08 by 5.08 cm EJ-309 liquid scintillators in 4 pairs, evenly-spaced with 90 degrees separation about the center of each assembly. OSCAR consists of a single array of 12, 5.08 by 5.08 cm trans-stilbene crystals aligned with the center of the assembly. Scatter-based systems detect fast neutrons without any moderation, leading to system die away times on the scale of tens of nanoseconds. This allows much shorter timing gate widths compared to thermal systems, thus increasing counting statistics of true correlated fission events. However, scatter-based systems are susceptible to neutron cross-talk when an incident neutron scatters off one detector and interacts in an adjacent detector, causing two seemingly correlated detection signals. The equations developed adjust for cross-talk to conduct NMC with scatter-based systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Performance of the Initial Implementation of the Shift Monte Carlo Code in SCALE 6.3

The Shift Monte Carlo code will be introduced in SCALE 6.3 as an alternative to the KENO V.a, KENO-VI, and Monaco codes. Calculations were performed to establish the performance of Shift for criticality safety analyses within the criticality safety analyses sequence (CSAS) based on models in the Verified, Archived Library of Inputs and Data (VALID). This test suite contains over 600 critical experiment models covering a broad range of fissile materials and neutron energy spectra. The comparisons presented include calculated k eff values and runtime performance for serial calculations and a selection of parallel calculations. Comparisons are presented for multigroup (MG) and continuous-energy (CE) calculations for KENO V.a and KENO-VI models. Results generated with a beta version of SCALE 6.3 indicate excellent agreement in k eff values between KENO and Shift. The largest differences in the average k eff value calculated for the 15 categories of KENO V.a models are 0.00020 ± 0.00011 Δ k for MG calculations and 0.00011 ± 0.00005 Δ k for CE calculations. Similar comparisons in three categories using KENO-VI result in the largest differences for MG calculations: as 0.00004 ± 0.00003 Δ k , and -0.00002 ± 0.00003 Δ k for CE. The preliminary results also indicate that Shift is faster than KENO on a per particle basis, especially for fast spectrum systems. The uncertainty per history is also higher, however, so the Monte Carlo figure of merit is higher for KENO for thermal and intermediate spectrum systems. As expected, Shift generally has better speedup than KENO for parallel calculations, regardless of neutron energy spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Update of the Nuclear Criticality Slide Rule Calculations: Plutonium systems – Delayed Fission Gamma

IRSN (France), LLNL (USA) and ORNL (USA) began a long-term collaboration effort in 2015 to update the nuclear criticality Slide Rule for the emergency response to a nuclear criticality accident. The Slide Rule permits the estimation of neutron and gamma dose rates and integrated doses based upon estimated fission yields, as a function of distance from the fission source, and time after criticality accidents for different critical systems. This paper presents results from the fourth phase of the current update of the Slide Rule project, in which delayed fission-product gamma (DFG) dose rates of unreflected plutonium critical systems were compared by several modern 3D radiation transport codes (MCNP, COG, SCALE), using updated flux-to-dose conversion factors. Dose rates are calculated for fissile material at five moderation ratios (H/Pu), at 1 m above the ground as a function of distance (between 30 cm and 1,200 m) from the external surface of the source to the center of the detector, and for periods between 1 s and 1,000 min after the critical instantaneous event. Further efforts have been devoted to the determination of the delayed gamma source, by comparing the time-dependent energy spectra obtained from several methods. Overall, DFG dose rates calculated by each participant led to consistent results. Extra effort is under way to identify the cause of the remaining differences, by comparing precisely the gamma source, and particularly nuclides inventories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Impact of Increased Latent Generations on Sensitivity Calculations with SCALE

Analyses of cross section sensitivity data from systems with fissile material allow analysts to associate an importance for each material, nuclide, reaction, and neutron energy by simulating real world criticality scenarios. Although criticality safety validation efforts can be guided by the cross-section sensitivity and uncertainty data generated for a particular system, these calculations can often be computationally expensive and sometimes cumbersome without proper guidance. The TSUNAMI suite within the SCALE code package has several methods for generating sensitivity data, including multigroup and continuous energy (CE) capabilities. The release of SCALE 6.3 has three different CE methods for generating cross section sensitivity data: (1) the Iterated Fission Probability (IFP) method with the KENO Monte Carlo transport solver, (2) the IFP method with the Shift Monte Carlo transport solver, and (3) the Contributon-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method with the KENO Monte Carlo transport solver. Although the CLUTCH method has additional parameters for generating sensitivity data files relative to the IFP method, all three methods use latent generations, which are the generations between an event (i.e., fission) and the assessment of importance based on the asymptotic population of progeny neutrons. Increasing the number of latent generations in a calculation leads to increased discrimination of the sensitivity coefficients but at the cost of the increased uncertainty associated with those generated values. Analysts must balance the accuracy of the sensitivity calculations and its uncertainty with the associated computational cost involved in generating the values. This paper discusses the impact of adjusting the latent generation parameter for a range of sensitivity values and how these changes compare with the direct perturbation values obtained from a change of ±0.5% Δ k in both benchmark and safety application models. Two benchmarks from the International Handbook of Evaluated Criticality Safety Benchmark Experiments and the MPC-32 dual purpose canister for spent nuclear fuel are used for analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Current Progress of the Final Design of a Subcritical Assembly at the Oak Ridge National Laboratory

During the past decade, the US Department of Energy (DOE)/National Nuclear Security Administration (NNSA) Nuclear Criticality Safety Program (NCSP) has been coordinating nuclear criticality safety (NCS) courses to provide training and qualification to new staff, supervisors, managers, regulators, and other professionals involved with NCS. To address the growing need for more NCS training possibilities, a subcritical assembly was designed at Oak Ridge National Laboratory (ORNL). The Oak Ridge Subcritical Assembly (ORSA), operated at ORNL, could be used in a new course for fissile material handlers/operators or as a backup training facility for currently existing NCSP hands-on NCS training. This paper describes some of the design steps necessary for obtaining a final ORSA design, beginning with a successful feasibility study performed in 2020. Numerous parameters within each main element of the assembly were considered, such as the reflector used, the type of fuel coating or plating used, and the number of fuel plates used. The design choices related to each parametric study described in this paper are detailed and justified at the end of each subsection. An overview of the current final design characteristics of ORSA and the experiments’ possibilities is also provided. ORSA will be built with fuel plates from legacy Aerojet General Nucleonics (AGN)-201M research reactors currently stored at the Y-12 National Security Complex (Y-12). The cylindrical fuel plates comprise UO 2 particles embedded in polyethylene and have a nominal enrichment of 19.5 ± 0.5 wt % 235 U. The basis for the final design was the feasibility study from 2020. In the feasibility study, ORSA achieved k eff = 0.95, corresponding to a multiplication factor of 20, with approximately 620 g of 235 U (four 4 cm plates, three 2 cm plates, and two 1 cm plates) and an 8.5 cm thick graphite reflector radially and on top of the assembly. The bottom of the assembly has a graphite thickness of 5 cm, and the chosen graphite density was 2.3 g/cm 3 . All design calculations were performed using SCALE 6.2.4 KENO-VI and with the ENDF/B VII.1 Continuous Energy cross-section library.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗