Experimental measurements of isothermal reactivity coefficient and temperature-dependent reactivity changes with associated uncertainty evaluations
Not provided.
SEARCH · Engineering Papers
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.
Not provided.
Nuclear fuel with extended enrichment (235U enrichment within 5-8 wt%) is one of the evolutionary changes that have been pursued in recent years by commercial light water reactor operators and fuel vendors to improve the fuel cycle economy and operation performance of a nuclear plant. This work assesses the performance of the Polaris/PARCS two-step approach in core physics modeling of the pressurized water reactor cores with extended enrichment fuel, referred to as “LEU+” in this report. A representative LEU+ core with a 24-month fuel cycle developed by Southern Nuclear Company (SNC) was modeled using this two-step approach. A representative LEU core with an 18-month fuel cycle was also modeled to provide a reference for the LEU+ core. As expected, significantly more burnable poison absorbers were used in the LEU+ core to accomodate its higher fuel enrichment. Nine different fuel assembly types were modeled using Polaris for each core to generate the assembly cross sections, which were then processed by GenPMAX to prepare the cross-section data for PARCS. The average specific powers of each fuel batch in each core were derived from VERA results and higher specific powers in fresh assemblies were found in the LEU core due to its less total uranium loading included in the VERA LEU model, given that the total core power was assumed to be the same for both cores. PARCS models were developed to simulate the steady-state operations of both cores. PARCS results on the LEU+ core were first compared with the VERA results for verification purpose; good agreements were seen in soluble boron and burnup distribution results, indicating that the Polaris/PARCS modeling and simulation were correctly implemented. Core physics parameters calculated by PARCS, at zero power physics tests, beginning of cycle (BOC), and end of cycle conditions (EOC), were compared between the LEU+ core and the LEU core, including soluble boron concentration, burnup distributions, assembly and pin power peaking factors, fuel temperature reactivity coefficients, moderator temperature and density reactivity coefficients, control rod worth, and shut down margin. The main differences in PARCS results between the LEU+ and the LEU cores are summarized below: 1)The critical boron concentrations were found to be much higher in the LEU+ core than in the LEU core (1582 vs. 1335 ppm for peak values). 2)Higher assembly radial power peaking factors (1.4 vs. 1.3 for peak values), 2D pin peaking factors (1.53 vs. 1.42 for peak values), and 3D pin peaking factors (1.89 vs. 1.81 for peak values) were found in the LEU+ core than in the LEU core. 3)Significantly higher reactivity coefficients of moderator temperature (and density) were found in LEU+ than in LEU.4)Significantly lower control rod worth at EOC were found in LEU+ than LEU for all but one control banks.5)Significantly lower shut down margins were found in the LEU+ core than in the LEU core. The Polaris/GenPMAX/PARCS code suite was found to be capable of modeling the LEU+ PWR core for steady-state operations and no unexpected results in core physics parameters were observed, in spite of that a) several bugs in PARCS were identified and workarounds were used; b) several features were found lacking in the current version of PARCS that would be useful for core modeling. A list of requests for bug fixes and feature upgrades for PARCS originated from this work were transmitted to the code developers. The assessments on the performance of the Polaris/PARCS two-step approach in core modeling for boiling water reactor with LEU+ fuel is ongoing.
Deimos was the first critical experiment using high-assay low-enriched uranium (HALEU) TRistructural ISOtropic (TRISO) fuel in over 40 years. HALEU TRISO is the desired fuel form for many of the advanced reactor designs in development; however, very little experimental data are available for this fuel type. Deimos was designed to utilize existing HALEU TRISO fuel in a large graphite moderator to obtain nuclear and reactor physics data to fill the gaps surrounding this fuel type and enrichment. In addition to cold critical data, three separate heated experiments were conducted to measure the temperature reactivity coefficient for this type of system. These measured coefficients were then compared to simulated coefficients to a first level order of fidelity. This comparison showed very good agreement for the experiment where only the inner core was heated and good agreement for the other two configurations, which included heating portions of the outer core. Less agreement when the outer core was heated is attributed to potential heating in the beryllium reflector, which has a positive temperature reactivity coefficient and was unaccounted for in the first-order models. Future heated experiments with Deimos will include temperature monitoring of the beryllium reflector to account for beryllium heating in the simulations.
Currently, studies on small modular reactors (SMRs) present an important development due to the potential they represent in terms of safety, operational flexibility, economy, and non-energy applications. Furthermore, there is increasing interest in studying the use of thorium as fuel, as an effective way to solve problems such as the shortage of uranium reserves, reduction of nuclear waste and nuclear proliferation. Also, the use of thorium combined with highly enriched uranium TRISO particles has been studied, showing proper performance. In this work, the concept of ThO{sub 2}-UO{sub 2} duplex fuel is used, for the core configuration of a PWR type SMR that uses TRISO fuel, designed to achieve extended fuel cycles. Three distribution cases of ThO{sub 2} and UO{sub 2} in TRISO particles inside the fuel rods are compared. First, particles composed by a homogeneous mixture of ThO{sub 2} and UO{sub 2} are distributed inside of fuel elements. Second, the fuel zone of the fuel elements is divided into two radially, an internal one where the TRISO particles composed of ThO{sub 2} are distributed and the external one where the TRISO particles composed of UO{sub 2} are distributed. The third case is like the previous one, except that the particles that contain UO{sub 2} are distributed in the inner zone and those that contain ThO{sub 2} in the outer zone of the fuel elements. The comparison of the cases is carried out in terms of cycle main isotopes' mass transmutation, power distribution and temperature reactivity coefficients. (authors)
As part of a US Nuclear Regulatory Commission–sponsored project to assess the modeling and simulation capabilities for accident progression, source term, and consequence analysis for advanced reactor technologies with SCALE and MELCOR, SCALE was used for the modeling and simulation of a molten salt-fueled reactor (MSR). SCALE capabilities for the modeling of MSR physics were demonstrated based on the Molten Salt Reactor Experiment (MSRE). Of primary interest were the determination of the system’s nuclide inventory, as well as the inventories in the various regions of the loop, considering that the fuel is continuously pumped through the system. This report contains discussions on the following: 1. Determination of the system-average fuel salt inventory considering fission gas removal in the off-gas system and noble metal removal through plating out at the heat exchanger using recent enhancements in SCALE’s depletion sequence TRITON, 2. Assessment of the nuclide spatial distribution throughout the loop using SCALE’s depletion solver ORIGEN, 3. Calculation of the core’s power profile, flux profile, temperature reactivity coefficients, and xenon reactivity using full-core calculations with SCALE’s Monte Carlo code KENO-VI. The results obtained with SCALE were post-processed to provide the MELCOR team with the core inventory and decay heat of the system, as well as the inventory and decay heat of individual regions in the loop, a zone-wise power profile, temperature feedback coefficients, and the xenon worth.
As part of a US Nuclear Regulatory Commission–sponsored project to assess the modeling and simulation capabilities for accident progression, source term, and consequence analysis for advanced reactor technologies with SCALE and MELCOR, SCALE was used for the modeling and simulation of a fluoride salt-cooled high-temperature reactor (FHR). Based on the preconceptual design for a small modular 236 MWth FHR developed by the University of California, Berkeley (PB-FHR-Mk1), a SCALE model of the PB-FHR-Mk1 reactor core was developed. The reactor was modeled at equilibrium state with different fuel compositions in different regions of the reactor. An iterative approach was used to interpolate and mix the burnup-dependent fuel compositions obtained through the depletion calculation of a core slice model. After demonstrating the applicability of SCALE’s multigroup (MG) approach for the simulation of the PB-FHR-Mk1, the resulting equilibrium core was studied in terms of the power profile, the flux profile, temperature reactivity coefficients, and the xenon reactivity. Furthermore, the tritium production rate in the salt coolant was determined, and the dependence of the one-group cross sections on the burnup and location in the reactor core was studied. The results obtained with SCALE were post-processed to provide the MELCOR team with the core inventory and decay heat of the equilibrium core, a zone-wise power profile, temperature feedback coefficients, the tritium production rate, and the xenon worth.
The system analysis of anticipated operating occurrences and design basis accident for pebblebed reactor systems requires knowledge of the neutron kinetic parameters. These parameters take into account various types of feedback from the pebble-bed core and are typically integrated into system analysis tools such as the System Analysis Module (SAM) in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tool package. These parameters are typically generated by a higher fidelity full-core coupled neutronics/thermal fluids analysis, such as using Griffin and Pronghorn, in 2D or 3D. The delayed neutron fractions (beta) and the neutron lifetimes (lambda) are typically generated by applying the adjoint solution of the neutron flux. Temperature reactivity coefficients are the most important feedback for all anticipated operating occurrences and design basis accident, then depending on the transient other effects might be important as rector control system worth, s-curves, and Xenon generation. This work package will establish a methodology and process for calculating the reactor kinetics parameters for high-temperature gas-cooled reactor and provide them in a suitable form for system level analyses of typical anticipated operating occurrences and design basis accident.
This paper compiles and analyzes refined results of the coordinated research project (CRP) on Neutronics Benchmark of China Experimental Fast Reactor (CEFR) Start-Up Tests conducted by the International Atomic Energy Agency (IAEA) since 2018. Twenty-eight research organizations participated with various code systems. The China Institute of Atomic Energy (CIAE) provided the benchmark specifications and participants conducted the benchmark through blind and refined phases. This paper presents the benchmark results on six experimental measurements, criticality, control rod worth, temperature reactivity coefficients, sodium void worth, assembly swap reactivities, and foil activations. Except for a few outlier results, the simulation results show good agreement with the measured data within 1-σ experiment uncertainties. (authors)
When designing a novel reactor, nuclear experiments are essential to validate the predictive capability of modeling/simulation tools and to justify the investment of a full-scale prototype. Being able to accurately predict temperature reactivity coefficients is of high importance to reactor designers as it impacts the safety and performance of the system. Deimos is a proposed graphite-moderated, beryllium-reflected, high assay low-enriched uranium (HALEU) tri-structural isotropic (TRISO) fueled experiment for NCERC. Deimos is a valuable experiment for validating the predictive performance of modeling and simulation tools since it comprises a relatively thermal energy spectrum with 84% thermal neutrons causing fission, 15% epithermal neutrons causing fission, and 1% fast neutrons causing fission. Additionally, by electrically heating the experiment, reactivity can be validated at various temperatures. Deimos will serve as a testbed for future advanced reactor concepts.
The work presented in this article is preliminary to downstream analysis of a generic fluoride salt-cooled high-temperature reactor (gFHR) core performed by the University of Tennessee in collaboration with Kairos Power (KP). A Monte Carlo transport model of the publicly available gFHR equilibrium core is developed in SCALE/KENO with multigroup energy treatment. Several output quantities of interest are used to verify the simulations against a benchmark model developed by KP using the continuous energy Monte Carlo code Serpent 2. Good agreement is seen in flux and fission rate profiles with a maximum relative difference of 1.4% and 2.8% respectively. Furthermore, an effective multiplication factor bias of 44 pcm was observed between the two simulations. The fuel temperature reactivity coefficient calculated with SCALE is within uncertainty to the reference model. This verification acts as a publicly reproducible benchmark for the gFHR in SCALE/KENO. A simplified depletion model is also presented where a single fuel pebble is depleted to discharge burnup through the equilibrium core while the equilibrium core is assumed to be invariant. This method produces results that intercept the equilibrium core concentrations in every case, however, an interesting artifact of this particular depletion model is uncovered. The phenomenon is shown to be a fundamental feature of the differential rate equations and inspires questions about how this system behaves when the time evolution of the equilibrium core is considered. (authors)
This work investigates nuclear reactor performance and safety characteristics of UO 2 with high thermal conductivity Mo insert structures by using multiphysics modeling techniques. Additionally, the purpose of this study is to use scoping analyses to quantify the impact of using Mo inserts from neutronic and heat transfer standpoints. Attention is given to reactor performance parameters, such as cycle length, maximum fuel temperature, temperature gradients in the fuel, and stored energy in the fuel. The finite-element code BISON and the Monte Carlo particle transport code Serpent were used to perform sensitivity analyses on the Mo insert geometry to optimize the insert design and inform larger scale modeling that required the homogenization of the UO 2 and Mo. Although BISON is often used as a fuel performance analysis tool, it is used in this context for heat transfer analysis only. Fuel performance optimization is outside the scope of the current study, but would be important for future work focused on this concept. The results showed that the insert had little impact on neutronic performance and that homogenizing the UO 2 and Mo was acceptable for reactor physics calculations. Reactivity temperature coefficients calculated using homogeneous UO 2 -Mo were shown to be relatively similar to UO 2 , but higher Mo content and 235 U enrichment can reduce the worth of soluble boron and control rods. The effect of insert geometry on heat transfer was much greater, and an approximately 15–20% difference in maximum fuel temperature was predicted between the best and worst performing heat transfer geometries. Furthermore, thermal conductivity calibration based on the finite element analysis results was performed to improve the accuracy of temperature predictions in reactor analysis models that homogenized the UO 2 -Mo fuel. Compared with UO 2 in a pressurized water reactor (PWR), the optimized UO 2 -Mo design increased the margin to fuel melt by 13–32% across the fuel cycle, but it requires the 235U enrichment to exceed 5% to match the cycle length of conventional UO 2 .
A preliminary reference core design was previously developed for the Versatile Test Reactor. The main features are that it is a 300 MW th sodium-cooled fast reactor using ternary metallic fuel U-20Pu-10Zr and able to achieve peak fast fluxes (E n > 0.1 MeV) in excess of 4.0x10 15 n/cm 2 -s. The plutonium in the fuel was assumed to be “reactor grade”, with a fissile quality of about 72% (i.e., 239 Pu and 241 Pu comprise 72% of the plutonium isotopes), and the uranium was assumed to be low-enriched uranium with 5% 235 U. This preliminary core design has been used through the CD-0 phase of the VTR project. Progression of the work after CD-0 led to revising, updating and refining the reference VTR core design, based on considerations from the various VTR teams: fuel, experiment, safety, and plant design teams. The objective of this report is to summarize the various changes made and details added to the model, to provide a description of the updated reference VTR core design and of its performance characteristics. This includes the traditional reactor physics characteristics, fuel cycle details, reactivity coefficients, temperature distributions, control rod worths and shutdown requirements. All of these define the updated “reference VTR core” intended to be used through CD-1.
Advanced reactor initiatives are growing significantly through programs nationwide. This research area includes small modular reactors, microreactors, and space reactors. Many of the reactors being designed are untested concepts. They include unique moderators, varying fuel types, high temperatures, and compact configurations. The shift in fuel type, from highly enriched uranium (HEU) to high assay, low enriched uranium (HALEU), is particularly important as it has driven many of the other changes. For example, lower enrichment requires advanced moderators, which in turn require different reflectors to make the systems compact. The change in materials including the transition from HEU to HALEU affects the temperature feedback of the systems. Additionally, these advanced reactor concepts generally have a thermal neutron spectrum in contrast to earlier fast spectrum advanced reactor. With the extensive changes from previous reactor designs, validation experiments are needed. The National Criticality Experiments Research Center (NCERC) is uniquely equipped to perform such experiments. The Deimos experiment, designed for execution at NCERC, will serve as a testbed for advanced reactor concepts. It will use HALEU fuel in a graphite matrix, provide the ability to use advanced moderators, and allow measurements of temperature reactivity coefficients (TRCs).
High temperature moderators for nuclear reactors and processes for their production are disclosed. The moderators include at least one hydrided metal and/or hydride metal allow, such as yttrium hydride, thorium hydride, yttrium-cerium hydride, yttrium-gadolinium hydride, yttrium calcium hydride, cerium hydride, etc. Such metal hydrides and/or hydride alloys may have high thermal stability, a relatively low thermal neutron absorption cross section, the ability to retain hydrogen over a large temperature range, and have good mechanical properties. Such moderators may induce spectral shift in reactors which, in turn, magnifies the Doppler reactivity temperature coefficient. Such moderators to thermalize neutrons may also enhance fuel utilization and cost-effectiveness of the reactor while keeping the core portable.
The capability to utilize neural networks to predict macroscopic and microscopic cross section parametric spaces has been developed for the Griffin reactor physics application. The LibTorch interface enables Griffin's MOOSE-based materials to interact with LibTorch-trained models, allowing for the evaluation of complex macroscopic or microscopic cross section spaces, which are then used to evaluate the neutronic properties of the Griffin finite element model. This study benchmarks traditional ISOXML-formatted tabulation libraries against neural network-based models for 279 nuclides on 20,160 grid points for zero-dimensional and two-dimensional reactor models. Benchmark metrics include the fundamental mode eigenvalue, fission and absorption rates, and various temperature coefficients of reactivity (isothermal, fuel, and moderator). From the perspective of storage space, the complete set of LibTorch models uses 11 MB on disk, compared to the 10 GB for the ISOXML multigroup library that covers the same grid space. For the two-dimensional performance case considered in Griffin, the Torch model uses 97% less RAM than the reference ISOXML dataset while runtime increases by a factor of 3 when using the LibTorch model compared to the ISOXML dataset with multi-linear interpolation. The LibTorch model consistently yields errors within 0.01% for most analyzed quantities except for the temperature coefficients of reactivity where the maximum discrepancies are up to 0.3 $\frac{pcm}{K}$. Due to the neural network attempting to best predict quantities with no regard for a positive or negative bias for any given quantity, predictions may experience random fluctuations, resulting in both positive and negative errors. Future work will entail both depletion and coupled transient analysis to determine the predictive capabilities of Griffin with neural network-based cross sections.
Godiva IV is a fast burst critical assembly located at the Nation Criticality Experiment Research Center (NCERC) in the Nevada National Security Site (NNSS). It is constructed of approximately 65 kg of highly enriched uranium (HEU) fuel alloyed with 1.5% molybdenum for strength. The assembly can be operated at delayed critical or can be used to perform super-prompt critical bursts with temperature rises of up to 250 °C. Several projects with the objective of characterizing Go diva IV are ongoing including characterizations of its radiation emission and thermomechanics. One of the characterization projects is a study of the variability in the relative source emission (or fluence in absolute terms) between bursts that are nominally the same size. There is some inherent variation in the number of fissions, or "size" of the burst, even for bursts where the operators perform identical processes. The source of the variation can be attributed to several factors. The largest contributor is the precision with which the operators are able to determine delayed critical (DC). Operators find DC for every burst operation. The impact of being slightly above or below DC can change the super-prompt critical reactivity and therefore the source term. The temperature of the fuel is another factor. Establishing delayed critical prior to performing a burst act as a control of the bulk temperature coefficient of reactivity.
Many advanced reactor designs require fuel enriched between 5% and 20% 235 U. To assist in producing fuel at these enrichment levels, government-owned inventories of highly enriched uranium can be downblended. However, fuel produced from these inventories contain uranium impurities that are not often found when enriching natural uranium or accounted for when modeling reactor cores. To address this concern, this work models reactor designs like the X-energy Xe-100 and the Ultra Safe Nuclear Company’s Micro Modular Reactor, and compares their performance with fuel from enriching natural uranium to fuel from downblended highly enriched uranium. This paper evaluates the models based on the effective neutron multiplication factor, k eff , effective delayed neutron fraction, β eff , and energy- and spatially dependent neutron flux, ϕ, as well as the fuel, coolant, moderator, and total reactivity temperature feedback coefficients, α F , α C , α M , and α T . The results show that the fuel from downblended highly enriched uranium inventories leads to differences in each of the metrics, especially in the keff values. In the Xe-100–like and Micro Modular Reactor–like models, k eff changes by about 1400 pcm and up to 1200 pcm, respectively. Total reactivity feedback coefficients α T are negative with the impure fuels and the keff values remain above 1 for each core configuration and fuel composition. These results show that the impure fuel compositions do not necessarily prevent achieving key design parameters, such as cycle length, or from operating in a safe condition.
Fluoride high-temperature reactors (FHRs) are high-temperature, low-pressure reactor concepts that use tri-structural isotropic (TRISO) fuel and molten fluoride salt coolant. These reactors have the potential to provide both electrical power and high-temperature process heat. We used generic FHR parameters for a pebble-bed FHR to develop an initial model with fresh fuel for a generic FHR (gFHR) in MELCOR and SCALE (NEWT and KENO). In this paper, we present the development of our gFHR models, which will serve as the baseline for a sensitivity and uncertainty analysis to quantify the range of possible source terms for FHRs in severe accidents. We present MELCOR results for fuel and coolant temperatures through the core, a nodalization study for the steady-state thermal hydraulic model, and development of reactor physics models in SCALE. As this work progresses, these models will be used to calculate source terms for a loss-of-forced-flow accident and to conduct a sensitivity study on this accident to establish a range of possible source terms. SCALE will provide reactor physics parameters like isotopic inventory, decay heat generation, and temperature coefficients of reactivity. Using the uncertainty quantification tools within SCALE, we will generate distributions for those parameters and will use the uncertainty quantification code RAVEN or DAKOTA to sample those distributions in MELCOR to quantify the impact of reactor physics and thermal hydraulic uncertainties on FHR source terms. (authors)