Search NASA⌕ Search

SEARCH · Search NASA

Results for “core reactivity”

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 91 records · Page 5

Core-flood Effluent and Shale Surface Chemistries in Predicting Interaction between Shale, Brine, and Reactive Fluid

We report field and laboratory observations to date indicate that the efficiency of hydraulic fracturing, as it relates to hydrocarbon recovery, depends significantly on geochemical alterations to rock surfaces that diminish accessibility by partial or total plugging of the pore and fracture networks. This is caused by mineral scale deposition such as coating of fracture surfaces with precipitates, particle migration, and/or crack closure due to dissolution under stress. In reactive flow-through experiments, mineral reactions in response to acidic fluid injection significantly reduced system porosity and core permeability. The present study focuses on changes to fluid chemistry and shale surfaces (inlet and fracture walls) resulting from shale-fluid interactions and integrating these findings for an improved estimate of transport-related consequences. The reacted shale surfaces were examined by spatially-resolved scanning electron microscopy - energy dispersive spectroscopy (SEM-EDS) analysis. Importantly, inductively coupled plasma - mass spectrometry/optical emission spectroscopy (ICP-MS/OES) was utilized to probe the chemical evolution of the core-flood effluents. The three study cores selected from the Marcellus formation represent different mineralogies and structural features. In flow-through experiments, lab-generated brine and HCl-based fracture fluid (pH=2) were injected sequentially under effective stress (up to 500 psi) at reservoir temperature (80°C). SEM-EDS results confirmed by the ICP concentration trends showed significant Fe hydroxide precipitates in clay- and pyrite-rich outcrop samples due to partial oxidation of Fe-bearing phases in the case of intrusion of low salinity water-based fluids. Porosity reduction in the MSEEL (Marcellus Shale Energy and Environmental Laboratory) carbonate-rich sample is related to compaction of cores under stress due to matrix softening with substantial dissolution, and pore-filling by hydroxides, as well as barite and salts. Despite the same fluid compositions and experimental conditions used for both MSEEL samples, barite precipitation was much more intense in the MSEEL clay-rich sample due to its greater sorption capacity and additional sulfate source as well as fissile nature with multiple lengthwise cracks. ICP tests revealed time-resolved concentration trends in produced brine and reactive fluids that in turn complemented the pre-/post-reaction SEM-EDS observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conceptual core design of the molten metal fueled microreactor with self-regulating capability

The microreactor concept has received significant attention in the United States for its lower capital investment, siting flexibility and high mobility. A microreactor aims to provide reliable electricity in remote communities, mining areas, or isolated islands, which will save the high fuel transportation cost. Thus, minimizing the staff level for operating and maintaining a microreactor is desirable. The objective of the study is to design a self-regulated Molten Metal fueled microReactor core (MMR) and to demonstrate the self-regulating capability of the MMR core. The control of the reactor is solely relying on reactivity feedbacks from fuel temperature perturbation. Molten UMn (liquid form at the temperature over 720°C) was utilized as the fuel material in MMR due to the larger thermal expansion coefficient compared to the conventional solid fuel. The MMR core has a fast neutron spectrum, which corresponds to a small burnup reactivity swing. The core design parameters were selected based on DAKOTA, REBUS, and ANLHTP calculations. The MMR is designed to run for 10 years at a power level of 15 MWth. Furthermore, no refueling is needed during the reactor lifetime. The system dynamic analyses were performed, and it was found that the reactor can be self-regulating within a temperature range of 800 ± 80°C for a loss of heat sink transient.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Reactions of U + with H 2 , D 2 , and HD Studied by Guided Ion Beam Tandem Mass Spectrometry and Theory

The kinetic energy-dependent reactions of the atomic actinide uranium cation (U + ) with H 2 , D 2 , and HD were examined by guided ion beam tandem mass spectrometry. An average 0 K bond dissociation energy of D 0 (U + – H) = 2.48 ± 0.06 eV is obtained by analysis of the endothermic product ion cross sections. Quantum chemistry calculations were performed for comparison with experimental thermochemistry, including high-level CASSCF–CASPT2–RASSI calculations of the spin–orbit corrections. CCSD(T) and the CASSCF levels show excellent agreement with experiment, whereas B3LYP and PBE0 slightly overestimate and the M06 approach badly underestimates the bond energy for UH + . Theory was also used to investigate the electronic structures of the reaction intermediates and potential energy surfaces. The experimental product branching ratio for the reaction of U + with HD indicates that these reactions occur primarily via a direct reaction mechanism, despite the presence of a deep-well for UH 2 + formation according to theory. The reactivity and hydride bond energy for U + are compared with those for transition metal, lanthanide, and actinide cations, and periodic trends are discussed. Furthermore, these comparisons suggest that the 5f electrons on uranium are largely core and uninvolved in the reactive chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiphysics Simulations of MSRE with NEAMS Thermal Hydraulics Tools

This report documents the benchmarks being developed and simulations performed using tools and codes developed under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, utilizing MSRE experimental data. In FY23, three main work scopes were investigated under the NEAMS MSR work package at ANL. The first scope investigated the Griffin-SAM coupling model for simulating the pump startup transient experiment of MSRE. The analyses start with a simple model (single-channel, single-lattice), gradually adding more details (multi-channel, full-core) into the model. The results show that the reactivity loss curve is very sensitive to the axial boundary conditions and the radial core discretization. The simple model can predict a similar reactivity trend as that of the more sophisticated model, which is likely due to error cancellation. Accurately modeling the axial boundary condition may further improve the reactivity trend but would require significant efforts to generate the mesh of the MSRE inlet and upper plenum. The core channel radial discretization for the Griffin-SAM coupled model also depends on the flow distribution. Given the complex geometry in the inlet plenum, the flow distribution needed to be calculated from CFD analysis, which was performed using the NekRS code. This analysis employed a MSRE CAD model developed by Copenhagen Atomics. The CAD model was disassembled to keep the inlet plenum region only, which was subsequently cleaned and modified so that the mesh generated is under the memory limit. The results are merged to a few radial regions to show that the flow rate is highest in the central region. This would be useful for future improvement of the Griffin-SAM coupling model of the MSRE core. The last task investigated is tritium transport modeling using the standalone SAM code. This task aimed to initiate the effort to demonstrate and validate the tritium transport model implemented in SAM. The preliminary investigation employed an MSRE model consisting of the primary loop. Three tritium transport pathways were examined including the retention in the graphite, the permeation through the HX tube wall, and the removal from the off-gas system. The results compare well with the MSRE data, but improvements are still needed on the initial conditions (i.e., the present state may not have reached equilibrium), the boundary conditions, the off-gas system modeling, and a better numerical strategy to reach the equilibrium state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Versatile Test Reactor Conceptual Core Design

The VTR is a 300-MW(thermal) sodium-cooled fast reactor (SFR) designed for the specific purpose of delivering unique testing capabilities to enable the advancement of all reactor technologies. With its flux level, irradiation volume, and operational flexibility, the VTR will enable accelerated testing of materials, fuels, and various components needing irradiation testing. Proven SFR technologies and design approaches have been leveraged in designing the VTR core, ensuring the highest possible readiness level. This resulted in the VTR using ternary metallic fuel and delivering fast flux levels in excess of 4 x 10 15 n/cm 2 ∙ s over large useful volumes, corresponding to about 60 dpa/year in steel. As part of the design efforts, the VTR core performance has been determined for a representative configuration, ensuring that the reactivity control systems offer sufficient shutdown margins, that the core can be safely cooled in all situations, and that reactivity feedback coefficients are conducive to a favorable safety behavior. Furthermore, the incorporation of features such as fuel assembly storage in the shield region supports the flexible and reliable operation of the VTR. Additional design work has been ongoing as well. This includes thorough shielding performance evaluations to ensure safe operation of the VTR, verification and validation of the design tools used to achieve compliance with Nuclear Quality Assurance (NQA-1) requirements, early assessment of the impact of irradiation experiments on the core performance envelope and associated margins, and in-depth uncertainty quantification efforts to quantify the anticipated range of performance characteristics. An experimental program supporting the VTR core design has been set up, with the current focus being on thermal-hydraulic experiments. The purpose of this experimental program is to obtain confirmatory measurements to serve directly as part of the core design basis or as part of the validation cases supporting the simulation tools used.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Aerosol Optical Tweezers Elucidate the Chemistry, Acidity, Phase Separations, and Morphology of Atmospheric Microdroplets

Aerosol particles represent unique chemical environments because of their high surface area-to-volume ratio that promotes the effects of interfacial chemistry in confined environments. Properties such as viscosity, diffusivity, water content, pH, and morphology—following liquid–liquid phase separation—can strongly alter how a particle interacts with condensable vapors and reactive trace gases, thus modifying its continual evolution and environmental effects. Our understanding of this chemical evolution of atmospheric particulate matter and its environmental impacts is largely limited by our ability to directly observe how these critical particle properties respond to the addition or reactive uptake of new chemical components. Aerosol optical tweezers (AOT) stably trap particles in focused laser beams, providing positional control and the retrieval of many of these critical properties required to understand and predict the chemistry of aerosolized microdroplets. The analytical power of the AOT stems from the retrieval of the cavity-enhanced Raman spectrum induced by the trapping laser. Analysis of the whispering gallery modes (WGMs) that resonate as a standing wave around the droplet’s interface, provide high accuracy measurements of the droplet’s size, refractive index (and thus a measurement of composition), and can distinguish between core–shell, partially engulfed, and homogeneous morphologies. In this study, we have advanced the ability to determine the properties of the core and shell phases in biphasic droplets, including obtaining high-accuracy pH measurements. These capabilities were applied to perform AOT physical chemistry experiments on authentic secondary organic aerosol (SOA) produced directly in the AOT chamber by ozonolysis of terpene vapors. The propensity of the SOA to phase separate as a shell from a wide range of nonpolar to polar core phases was observed, along with the discovery of a stable emulsified state of SOA particles in an aqueous salt droplet. Micron-thick SOA shells did not impede the gain or loss of water or squalane from the core to the surrounding air, indicating no significant diffusional limitations to condensational growth or partitioning even under dry conditions. These experiments formed the foundation of a new framework that predicts how the phase-separated morphology of complex aerosols containing organic carbon evolves during continual atmospheric oxidation processes. Increases in oxidation state will quickly drive conversion from a partially engulfed to core–shell morphology that has dramatically different chemical reactivity since the core phase is completely concealed by the shell. The recent advances in the experimental capabilities of the AOT technique such as presented here enable novel experimental methodologies that provide insights into the chemistry and multidimensional properties of aerosol microdroplets, and how these coevolve and respond to continual chemical reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A New Proposal Generalized Predictive Control Algorithm With Polynomial Reference Tracking Applied for Sodium Fast Reactors

This paper proposes a generalized predictive control (GPC) with constraints and orthonormal Laguerre functions using the simplified model of the primary system (reactor core and intermediate heat exchanger (IHX)) of a prototypical sodium fast reactor (SFR). This paper develops a multiple-input multiple-output (MIMO) GPC with input constraints able to track polynomial references of any degree applied in coolant temperature difference across the core and fractional power. The manipulated variables of the GPC-SFR are the reactivity and the sodium flow rate of the primary and secondary pipes. Moreover, orthonormal Laguerre functions and step down condition number techniques were also applied to avoid the numerical ill-conditioning issue in quadratic programming of large systems. Thus, a GPC type-2 was designed to control fractional power, coolant temperature difference across the core and sodium tank temperature of the SFR primary system when temperature references change according to a linear ramp after reaching their steady-state operation, sustaining 100% power operation on the reactor. In order to analyze the load tracking capability of the GPC-SFR type-2, the load following from 100% fractional power (FP) to 60% FP at 0.8% FP/min rate is simulated. Constraints on the rate of coolant temperature difference across the core and reactivity were applied for the design safety. For comparison criteria, this paper compares the GPC-SFR type-2 with the GPC-SFR type-1, i.e., standard model predictive control (MPC), to verify the viability and superior performance of the proposal regarding: (a) ramp-tracking capability of temperature and load; (b) the rejections of a reactivity disturbance of -1 cent and a secondary sodium inlet temperature disturbance of +10°F; and (c) a simulation with uncertainty in reactor design. The simulations show that the GPC-SFR type-2 overcome the GPC-SFR type-1 robustness and performance.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Modeling a generic TRISO-fueled heat pipe microreactor using SCALE: Depletion, transportation criticality, and shielding

This paper demonstrates the applicability of the SCALE code system to tristructural-isotropic (TRISO)-fueled heat pipe microreactors through depletion, transportation criticality, and shielding analyses of a generic reference design. The study conducted supports US Nuclear Regulatory Commission code readiness efforts for advanced non–light-water reactor technologies and is intended as a code capability demonstration rather than as an optimization of a specific microreactor design. The modeled reactor employs high-assay low-enriched uranium (HALEU) uranium oxycarbide (UCO) TRISO fuel and beryllium oxide (BeO) reflectors and operates at 7.5 MWth with a nominal lifetime of about 3 effective full power years. Representative cases for fresh and irradiated cores were selected to exercise SCALE methods relevant to reactor operation and post-irradiation transport. The discharged-core decay heat is approximately 6% of operating power immediately after shutdown. Transportation criticality calculations show that internal water ingress is the dominant reactivity effect, with fully flooded fresh core and irradiated core configurations remain above the subcriticality criterion, even with the available control mechanisms. Shielding calculations for a simplified transportation package indicate that normal-condition dose rates are governed mainly by shielding thickness and cooling time, whereas the breached hypothetical accident case is governed primarily by cooling time. Overall, the study shows that SCALE supports depletion, transportation criticality, and shielding evaluations efficiently for TRISO-fueled heat pipe microreactors within a single code system.

Criticality↗

Pnictogen Bonding at the Core of a Carbene-Stiborane-Gold Complex: Impact on Structure and Reactivity

Our interest in the design of ambiphilic ligands and their coordination to gold has led us to synthesize an indazol-3-ylidene gold chloride complex functionalized at the 4-position of the indazole backbone by a stibine functionality. The antimony center of this new complex cleanly reacts with o-chloranil to afford the corresponding stiborane derivative. Structural analysis indicates that the stiborane coordination environment is best described as a distorted square pyramid whose open face is oriented toward the gold center, allowing for the formation of a long donor–acceptor, or pnictogen, Au → Sb bonding interaction. The presence of this interaction, which has been probed computationally, is also manifested in the enhanced catalytic activity of this complex in the cyclization of N-propargyl-4-fluorobenzamide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scalable modular dynamic molten salt reactor system model with decay heat

Nodal dynamic models allow rapid reduced-order simulations of complex systems. This paper presents a publicly available system dynamic model of a molten salt reactor system with a topology of the Molten Salt Reactor Experiment (MSRE), in MATLAB/Simulink, building upon a model previously published and validated with MSRE data. It adds scaling of nominal power, dynamic representation of decay heat generation, a generic decay heat removal system, improved documentation, modularity, and test cases. One hour of system transient runs in less than 10 min on a laptop. This model is intended for first order engineering calculations, developing insights into system transient behavior, equipment sizing, parameter sensitivity studies. The model contains an easy to use toolkit which can be customized and expanded to describe any molten salt reactor system. Finally, scenarios involving off-normal transients are presented. Relation between shutdown reactivity insertions, decay heat removal, core recriticality, and safe shutdown conditions are discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Data and Insights for the KRUSTY Nuclear-Powered Tests

This document, and the associated files, describe all of the data from the KRUSTY warm criticals and the KRUSTY nuclear system test. The warm criticals were a group of tests where fission power heated the core and the feedback response was studied. There were 3 experiments, the 15-cent run (which was a true “free run”), followed by 30-cent and 60-cent runs, where reactivity was added to heat the core to higher temperatures. The nuclear system test (often referred to as the full-run or the final-run), engaged the power conversion system and evaluated system performance. These experiments are described in a pair of NUCLEAR TECHNOLOGY papers, which includes plots of the data that this document-set includes. Note that a lot of information in this document is based on my interpretation and could be wrong, but I’ve included it to potentially save some people a lot of time and headaches by trying to figure out which instrumentation to believe, but beware treating anything I say as fact. I’ve also included several post-irradiation photos, to show some thermocouple locations and indicate possible changes in emissivity and/or geometry (although no changes in geometry were indicated).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

SAS4A/SASSYS-1 Validation with FFTF Tests

The FFTF liquid sodium cooled test reactor was built to assist the development and testing of advanced fuels and material for fast reactors. The goals of the FFTF Passive Safety Testing program included confirming the safety margins of FFTF as a liquid metal reactor, providing data for computer code validation, and demonstrating the inherent and passive safety benefits of its specific design features. The program included thirteen unprotected loss of flow without scram tests. In preparation for the Passive Safety Testing program, a series of individual reactivity feedback tests were performed which were designed to simulate and validate specific features and reactivity feedbacks of the FFTF core.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Validation of the SCALE/Polaris–PARCS Code Procedure With the ENDF/B-VII.1 AMPX 56-Group Library: Boiling Water Reactor

The SCALE/Polaris–PARCS code procedure has been used in the confirmatory analysis for boiling water reactors by the US Nuclear Regulatory Commission. In this study, the SCALE/Polaris v6.3.0–PARCS v3.4.2 code procedure with the Evaluated Nuclear Data File (ENDF)/B-VII.1 AMPX 56-group library was validated by comparing the simulated results with the measured data for operating boiling water reactors, including Peach Bottom Unit 2 cycles 1–3, Hatch Unit 1 cycles 1–3, and Quad Cities Unit 1 cycles 1–3. The uncertainties and biases of the SCALE/Polaris–PARCS code package for boiling water reactor physics analysis were evaluated in the validation for key nuclear parameters such as reactivity and traversing in-core probe data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Review of Tools for Modeling Core Radial Expansion in Liquid Metal-Cooled Fast Reactors

Core radial expansion in liquid-metal cooled fast reactor systems is a well-known phenomenon that produces strong reactivity feedback effects. It is crucial from a safety standpoint to design a reactor with the proper core restraint system to guide expansion of the fuel assemblies into a formation that produces negative reactivity feedback in accident conditions. However, the modeling of core radial expansion and its associated reactivity feedback is extremely challenging. Liquid metal-cooled fast reactor designers (both industry and commercial) have pointed to improved modeling of core radial expansion as a key modeling challenge for these types of reactors. The Department of Energy Nuclear Energy Advanced Modeling and Simulation (DOE-NEAMS) program has an important role to play in helping to develop a robust, predictive tool for this phenomenon. A literature review has been performed to identify current simulation capabilities applicable to the prediction of core radial expansion in liquid metal-cooled fast reactors and the resultant reactivity feedback effects. An overview as well as a detailed description of the multi-physics phenomena underlying core radial expansion is given in the first two chapters. The objective of this report is to identify existing software that has been or could be applied to the core radial expansion problem, and then to identify possible paths to develop a robust, coupled Multiphysics tool for advanced analysis, in which all three physics (radiation transport, structural mechanics, thermal fluids) are taken into account considering feedback effects from each physics to the other. Several computer codes and workflows have been developed in the United States and abroad to perform individual aspects of radial expansion calculations or to perform loosely coupled analysis. No code system currently exists which tightly and robustly couples the neutronics, thermal hydraulics, and thermal mechanical physics with enough detail to fully resolve the complex core radial expansion reactivity feedback effects. The future availability of such a code system is of vital importance to fully understanding the reactivity feedback effects that occur due to radial expansion, and consequently to optimizing the design of the core restraint system. A potential code development path forward using MOOSE-based tools is suggested in order to leverage the natural tight coupling and robustness of MOOSE-based applications. Additionally, potential gaps for modeling and common assumptions are listed which require further investigation to ensure accuracy for prediction of deformations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Impact of U-10Mo HALEU fuel element tolerances on the Massachusetts Institute of Technology reactor safety and operational performance – Thermal hydraulics

The U.S. is coordinating efforts for the conversion of six U.S. High Performance Research Reactors (USHPRR) including one critical facility from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts of Institute of Technology Reactor (MITR), and achieve similar performance, high-assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being evaluated. The impact of the fabrication specification and tolerances was assessed following the preliminary design of the MITR LEU fuel elements using the U-10Mo monolithic alloy. This research focuses on the analysis of fabrication specification impact on thermal hydraulics (TH) characteristic of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., coolant channel gap thickness, fuel plate thickness, etc.). The analyses are performed based on an all-fresh LEU fuel conversion plan identified in a preliminary safety analysis report submitted to the Nuclear Regulatory Commission. The reactor power margin to the onset of nucleate boiling (ONB) is assessed under the limiting safety system settings (LSSS), where a scram occurs, to ensure there is sufficient margin to the reactor safety limit, which is defined by the onset of flow instability that occurs after the ONB. The best estimate plus uncertainty approach is employed to analyze this TH characteristic, which yields realistic results while maintaining adequate conservatism, utilizing a statistical uncertainty propagation method with the STAT7 code. The TH characteristic is analyzed as a function of the variability of the specification parameters resulting from the fabrication process. The main findings of this study show that the MITR core can meet the TH safety and operational requirements at the all-LEU initial core startup (cycle 1), selected transition cycles (most reactive cycle and most limiting cycle: cycle 3 and 5, respectively) and equilibrium (cycle 14) cores under all limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core power margin to ONB on those specification parameters that have the most direct impact on TH performance is non-linear but monotonically decreasing within the specification tolerances. The third order polynomial fit curves are reported in detail for selected limiting cases and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion↗

Experiments to Measure the Effect of Tantalum on Critical Systems

Sandia National Laboratories (SNL) and Oak Ridge National Laboratory (ORNL) collaborated to develop a capability to test the epithermal/intermediate cross sections of materials at the SNL critical experiment facility using the Seven Percent Critical Experiment (7uPCX) fuel. As a result, a new set of critical experiments has been designed to target the epithermal cross sections of tantalum (Ta) and is scheduled to be performed at SNL in 2023. These critical experiments will be evaluated for inclusion in the 2024 edition of the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook. The focus of these critical experiments is to measure the effects of Ta on the critical array size. The critical array size will be determined by an approach-to-critical experiment with the number of fuel rods in the array as a free parameter. The core configurations are designed to optimize the reactivity worth of Ta and the overall percentage of Ta absorption rates in the epithermal/intermediate energy range (0.625 eV – 100 keV). The baseline core configuration includes 7uPCX fuel rods set at a triangular pitch of 1.016 cm and a central dry test region that utilizes a cadmium liner for filtering out thermal neutrons. The central test region has locations for 85 Ta rods set at a triangular pitch of 0.813 cm. The Ta reactivity worth for the case with 85 Ta rods is approximately 2.55%, with the percentage of Ta absorption rates within the intermediate energy range at nearly 90%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Assessment of Core Physics Characteristics of Extended Enrichment and Higher Burnup LWR Fuels using the Polaris/PARCS Two-Step Approach. Vol. I: PWR Fuel

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.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗