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 data validation is often performed today using criticality measurements. The gold standard for criticality measurements is the International Criticality Safety Benchmark Experiment Project (ICSBEP). The validation specifically focuses on the effective multiplication factor ($k_{eff}$). $K_{eff}$ is a relatively easy parameter to infer and has reduced uncertainty due to being at or above critical. However, while $k_{eff}$ is the most documented parameter and its uncertainties and sensitivities have been evaluated in great detail, it cannot be used as a standalone metric to determine inaccuracies in nuclear data (e.g., cross section data, PFNS, nu), which is based on theory, physics, and differential measurements. The Experiments Underpinned by Computational Learning for Improvements in nuclear Data (EUCLID) project aims to identify compensating errors in specific isotope nuclear data by optimally designing experiments that are, or are not sensitive to a suite of measurement parameters beyond $k_{eff}$. By identifying parameters that are sensitive to each other, oppositely sensitive, or have substantial magnitude differences in sensitivity, experiments can be designed to constrain questionable nuclear data. One sensitivity that is of particular interest to this project includes the sensitivity of reactivity coefficients. Reactivity coefficients compare reactivity, which is related to $k_{eff}$ at two different states therefore being sensitive to small changes in the system. The most common type of reactivity coefficient measurements is comparison to void for a small sample within the assembly. It is key that the sample sizes are small enough to not affect the flux of the full system. Reactivity coefficients were evaluated for many early experiments to better understand transport corrected cross sections. In fact, ICSBEP includes reactivity coefficient results as “Supplemental Measurements” in appendices for a handful of older benchmarks. One of those benchmarks is Jezebel, the bare Pu critical assembly. This paper compares new simulations of reactivity coefficients for Jezebel, and explores the sensitivity of reactivity coefficients to small changes in nuclear data.
Core design studies on a 20-MW thermal power low enriched uranium (LEU) fueled research reactor are underway at the NIST Center for Neutron Research (NCNR) with the primary purpose of providing quality neutron beams for scientific research. The 'tank-in-pool' design pattern is adopted for the new reactor: A compact core composed of two horizontally split halves are surrounded by a cylindrical D{sub 2}O reflector tank, and the tank is itself surrounded by a large H{sub 2}O pool, which functions as thermal and biological shielding. The split core consists of 18 MTR-type fuel elements which are evenly distributed into two regions. The LEU fuel used in the study is the U{sub 3}Si{sub 2}-Al dispersion fuel with U-235 enrichment of 19.75%, the only NRC certified LEU fuel in the U.S. The core is cooled and moderated by H{sub 2}O, and separated from D{sub 2}O by zircaloy core boxes. Hafnium control blades are utilized as both criticality and safety control elements for the reactor. They are controlled by a mechanical driver located at the bottom of the core but with the fully withdrawn positions at the top of the core. In a standard reactor calculation, it is necessary to evaluate the reactivity coefficient with respect to some physical parameters such as fuel or moderator temperature. For nearly all reactor designs, negative reactivity coefficients are required for safety. As a light-water cooled thermal reactor, the new reactor is designed to work in an under-moderated mode with the most achievable reactivity. The core is required to have negative reactivity coefficients under any normal or abnormal circumstances such as over power transient, loss of coolant or reflector, etc. Therefore the reactivity coefficient addresses reactor safety concerns as well as provides valuable quantitative parameters for reactor safety analyses. In this summary, three typical reactivity coefficients for the low-power LEU core are calculated using the Monte-Carlo code MCNP-6: the moderator temperature coefficient (MTC), the fuel temperature coefficient (FTC), and the void coefficient (VC). A direct perturbation methodology is applied to calculate all these coefficients. The perturbed case is produced by manually perturbing a single physical property with respect to the reference case (the one with normal operational conditions), as the multiplication factor (i.e. k{sub eff}) of the perturbed case is calculated by MCNP. The corresponding reactivity coefficient is thereby determined by dividing the reactivity change over the associated parameter changes. Core design studies are underway at NIST to develop a low-power LEU fueled beam reactor to advance the neutron source capability at NCNR for the next century. The reactivity effect caused by moderator temperature, fuel temperature and void is investigated though a direct perturbation approach using the MCNP code. The MTC, FTC, and VC for a representative startup core is obtained and the negative effect characteristics of these coefficients are observed. The reactivity effects under several postulated abnormal scenarios are also investigated. The overall negative feedback to the reactivity at each situation is confirmed. (authors)
The activity coefficients of a reactor facility with VVER in normal operating regimes associated with a change of the external load are studied. The model used to evaluate the reactivity coefficients in terms of the moderator density and fuel temperature aided by the connected KORSAR/GP code is described. A three-dimensional model of the neutron kinetics and a realistic model of the main regulators are used. Calculations of the reactivity coefficients for a regime with unloading of the turbogenerator and the regime of primary regulation of the power system frequency are studied as test problems. A physical analysis of the results is performed.
Compensating errors between several nuclear data observables in a nuclear data library can adversely impact application simulations. The primary goal of the EUCLID project (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) is to reduce compensating errors between fast (0.1–5 MeV) 239Pu nuclear data for prompt fission neutron spectra (PFNS), average prompt fission neutron multiplicities, and neutron induced fission, capture, elastic, and inelastic cross sections. This work will focus on the design and execution of void reactivity coefficient measurements in the EUCLID experiment, performed on the Planet vertical lift critical assembly machine at the National Criticality Experiments Research Center (NCERC). Two different base configurations were designed and measured, one with high neutron leakage, and one with low neutron leakage. Both were primarily made up of plutonium metal (Zero Power Physics Reactor plates) without interstitial moderators and reflected by half-inch aluminum. Design optimization showed that void reactivity coefficient measurements in three locations per configuration was most impactful to reduce nuclear data uncertainties due to the varying impacts from elastic and inelastic scattering, as well as fission and capture. The locations for measurements were chosen based on preliminary studies which balanced measurement uncertainty and measurement practicality. The measurements were also selected to have sensitivities maximally complementary to previous arrangements. Comparisons across nuclear data libraries highlight the potential impact.
We present an efficient method to estimate the fuel temperature reactivity coefficient (FTC) by the Monte Carlo adjoint-weighted correlated sampling method. In this method, a fuel temperature change is regarded as variations of the microscopic cross sections and the temperature in the free gas model which is adopted to correct the asymptotic double differential scattering kernel. The effectiveness of the new method is examined through the continuous energy MC neutronics calculations for PWR pin cell problems. The isotope-wise and reaction-type-wise contributions to the FTCs are investigated for two free gas models - the constant scattering cross section model and the exact model. It is shown that the proposed method can efficiently predict the reactivity change due to the fuel temperature variation. (authors)
A new set of critical experiments exploring the temperature-dependence of the reactivity in a critical assembly is described. In the experiments, the temperature of the critical assembly will be varied to determine the temperature that produces the highest reactivity in the assembly. This temperature is the inversion point of the isothermal reactivity coefficient of the assembly. An analysis of relevant configurations is presented. Existing measurements are described and an analysis of these experiments presented. The overall experimental approach is described as are the modifications to the critical assembly needed to perform the experiments.
Two experiment series are planned to measure temperature effects in the Sandia Critical Experiments. The first series (IER-304) will measure the critical size of a fuel rod configuration at several temperatures. The second series (IER-452) will measure the inversion temperature of the isothermal reactivity coefficient in the Sandia critical assembly. Each experiment in the second series will be preceded by one or more experiments in the first series. This presentation addresses the conceptual design of the inversion temperature experiments.
Nuclear data validation is often performed today using criticality measurements. The gold standard for criticality measurements is the International Criticality Safety Benchmark Experiment Project (ICSBEP). The validation specifically focuses on the effective multiplication factor (k eff ). K eff is a relatively easy parameter to infer and has reduced uncertainty due to being at or above critical. However, while k eff is a well-documented parameter with detailed sensitivity and uncertainty analysis, it cannot be used as a standalone metric to determine inaccuracies in nuclear data (e.g., cross section data, PFNS, nu), which is based on theory, physics, and differential measurements. The Experiments Underpinned by Computational Learning for Improvements in nuclear Data (EUCLID) project aims to identify compensating errors in specific nuclear data by optimally designing experiments that are sensitive to a suite of measurement parameters beyond k eff . By identifying how each parameter's nuclear data sensitivity differs from others, experiments can be designed to constrain questionable nuclear data. One sensitivity that is of particular interest to this project includes reactivity coefficient sensitivities.
An assessment of the coolant reactivity feedback influence on a small Lead-cooled Fast Reactor (LFR) dynamics has been made aimed at providing both qualitative and quantitative insights into the system transient behavior depending on the sign of the above mentioned coefficient. The need of such an investigation has been recognized since fast reactors cooled by heavy liquid metals show to be characterized by a strong coupling between primary and secondary systems. In particular, the coolant density and radial expansion coefficients have been attested to play a major role in determining the core response to any perturbed condition on the Steam Generator (SG) side. The European Lead-cooled System (ELSY)-based demonstrator (DEMO) has been assumed as the reference LFR case study. As a first step, a zero-dimensional dynamics model has been developed and implemented in MATLAB/SIMULINK{sup R} environment; then typical transient scenarios have been simulated by incorporating the actual negative lead density reactivity coefficient and its opposite. In all the examined cases results have shown that the reactor behaves in a completely different way when considering a positive coolant feedback instead of the reference one, the system free dynamics resulting moreover considerably slower due to the core and SG mutually conflicting reactions. The outcomes of the present analysis may represent a useful feedback for both the core and the control system designers. (authors)
One of the research activities to support the commercial radioisotope production program is a safety research target irradiation FPM (Fission Product Molybdenum). FPM targets form a tube made of stainless steel in which the nuclear degrees of superimposed high-enriched uranium. FPM irradiation tube is intended to obtain fission. The fission material widely used in the form of kits in the world of nuclear medicine. Irradiation FPM tube reactor core would interfere with performance. One of the disorders comes from changes in flux or reactivity. It is necessary to study a method for calculating safety terrace ongoing configuration changes during the life of the reactor, making the code faster became an absolute necessity. Neutron safety margin for the research reactor can be reused without modification to the calculation of the reactivity of the reactor, so that is an advantage of using perturbation method. The criticality and flux in multigroup diffusion model was calculate at various irradiation positions in some uranium content. This model has a complex computation. Several parallel algorithms with iterative method have been developed for the sparse and big matrix solution. The Black-Red Gauss Seidel Iteration and the power iteration parallel method can be used to solve multigroup diffusion equation system and calculated the criticality and reactivity coeficient. This research was developed code for reactivity calculation which used one of safety analysis with parallel processing. It can be done more quickly and efficiently by utilizing the parallel processing in the multicore computer. This code was applied for the safety limits calculation of irradiated targets FPM with increment Uranium.
Abstract not provided.
A plan is being developed for the conversion of the NIST research reactor (NBSR) from high-enriched uranium (HEU) fuel to low-enriched uranium (LEU) fuel. The LEU fuel may be a monolithic foil (LEUm) of U10Mo (10% molybdenum by weight in an alloy with uranium) or a dispersion of U7Mo in aluminum (LEUd). A previous report provided neutronic calculations for the LEUm fuel and this report presents the neutronics parameters for the LEUd fuel. The neutronics parameters for the LEUd fuel are compared to those previously obtained for the present HEU fuel and the proposed LEUm fuel. The results show no significant differences between the LEUm and the LEUd other than the LEUd fuel requires slightly less uranium than the LEUm fuel due to less molybdenum being present. The calculations include kinetics parameters, reactivity coefficients, reactivity worths of control elements and abnormal configurations, and power distributions under normal operation and with misloaded fuel elements.
For a nuclear reactor, the temperature coefficient of reactivity (TCR) is of significance for both power control and safety analysis, but it is not always a constant. For example, the TCR may be dependent upon temperature due to Doppler broadening, and may be also dependent upon time scale due to heat transfer and thermal expansion. A fast burst reactor is a favorable nuclear facility to simulate the time scale effects on the TCR. Fast burst reactors are generally unmoderated research reactors of enriched uranium metal or uranium alloys. The neutrons are fast and the Doppler broadening effects can be neglected, and the TCR is mainly contributed by fuel thermal expansion. Fast burst reactors can operate on both steady state power mode and burst power mode. In steady state power mode, the systems are assumed to be in thermal equilibrium. In burst mode, the full width at half maximum (FWHM) of the power is from tens of microseconds to a few milliseconds, and the heat transfer is generally neglected. If FWHM is short enough, thermal expansion will have a time delay due to inertial effect. This paper is aimed to study the exact time scale effects on the TCR due to thermal expansion. The TCR contributed by the Doppler broadening effects was neglected, because it is relatively small (The k{sub eff} as a function of temperature was also calculated using the Monte Carlo particle transport code JMCT by our group, and no definite difference due to the temperature change was obtained near room temperature), and there is no method to evaluate the time scale dependence of the Doppler broadening effects yet. In the time scale from 1 x 10{sup -6} s to 1 x 10{sup 5} s, the TCR of a Godiva-like spherical fast burst reactor was calculated with perturbation theory. A Godiva-like spherical fast burst reactor was studied, the calculation of the time scale effects on TCR involved thermal conduction, thermoelasticity and reactivity variation analysis. Both the thermal and the thermoelastic problems were solved using ANSYS software. ANSYS is software implementation of finite element simulations, which is widely used in engineering field. The reactivity variation due to thermal expansion was calculated with perturbation theory.
The present paper discusses applicability of the measured data of Joyo cores from a view point of integral validation for the 4S nuclear design methodology. Through the evaluation of isothermal reactivity coefficients and reactivity losses due to burnup, the results confirm that those MK-I and MK-II database are effective in order to increase the dataset for uncertainty estimation for the prediction. Discussions on the 4S design method validation are also done through the analyses of criticality, power distributions and reactivity loss due to burn-up. The C/E values for criticality and reaction rate distributions are confirmed to be consistent with those obtained from the physics benchmark experiments. Through an analysis of burnup coefficient of the MK-I core by the detailed Monte Carlo calculations, the C/E value is 1.1, which is close to 1.06 obtained by the deterministic transport analysis. (authors)
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.