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.
In this study, the transient response of a reactor in the case of reactivity insertion with and without reactivity feedback was explored. The uncertainties resulting from using various numerical solvers to simulate the reactor response using the point reactor kinetics model were investigated.
Solution of reactor kinetics equation numerically unconditionally stable for all values of reactivity or integration step size
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.
Accurately predicting the criticality of an experiment before interacting with the experimental components is very important for criticality safety. Radiation transport software can be utilized to calculate the effective neutron multiplication factor of a nuclear system. Because of the integral nature of the effective neutron multiplication factor, the value calculated contains various sources of nuclear-data induced uncertainty. The sensitivity analysis and data assimilation technique presented in this paper exhibit one possible method of identifying and reducing the effective neutron multiplication factor nuclear-data induced uncertainty. The results presented in this work show that it is possible to use relative sensitivity coefficients of the prompt neutron decay constant and the effective delayed neutron fraction to 239 Pu nuclear data to reduce nuclear-data induced uncertainties in the effective neutron multiplication factor. This work has been utilized by members of the Los Alamos National Laboratory project EUCLID (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) for optimally designing a new experiment, which will be used to reduce compensating errors in 239 Pu nuclear data.
Sensitivity coefficients of prompt neutron decay constant (α) and effective delayed neutron fraction (β eff ) to Pu-239 nuclear data were calculated using a newly available tool (ACEtk). Uncertainty calculations showed the following trends: (1) The prompt neutron decay constant can be used to reduce nuclear data-induced uncertainty in Pu-239(n,f) and (2) the effective delayed neutron fraction can be used to effectively reduce nuclear data-induced uncertainty in Pu-239 total fission $\tilde{v}$. Sensitivity/uncertainty analysis can be used to determine optimal experiment/detector setup to identify compensating errors in nuclear data. Upcoming work includes investigating other nuclide-reaction pairs and response sensitivities, such as Δ$\rho$, neutron leakage spectra, and reaction rate measurements to constrain nuclear data of interest.
Abstract not provided.
Abstract not provided.
The use of electrons as main reagent for the recovery and recycling of critical metals from spent lithium-ion batteries (LIBs) is a process electrification strategy that can be used to close the life-cycle loop of LIBs through more sustainable methods. Electrochemical leaching, a process that uses a reductant that is constantly regenerated electrochemically for the leaching of lithium-ion battery black mass (LIBBM), has shown high extraction efficiencies and sustainable scores. However, slow kinetics, reactor design challenges and lack of deeper understanding of the underlying processes are barriers to the optimization, scale-up, and market adoption of this technology. In this paper, a kinetic study and mathematical model for dissolving LIBBM is presented to better understand the underlying mechanisms aiming to reduce the processing time and make predictions for future design and scale-up. The effect of acid and electrochemically mediated reductant concentrations, LIBBM loading, and cathode/reactor designs were explored. As a result, the leaching time was reduced from 7h to under 1h at a pulp density of 73 g/L, without external heating. A novel reactor with parallel baffle electrodes (PBE) was developed, which significantly reduced the leaching time by improving convection in a stirred slurry electrochemical reactor. Dimensionless numbers were deduced from an unsteady state model, which can be used in dimensional analysis for future process design and scale-up.
The systems analysis of anticipated operating occurrences and design basis accidents for pebble-bed reactor systems requires knowledge of neutron point kinetics equations (PKE) parameters. Typically, the generation of PKE parameters is performed in a global manner using standalone neutronics calculations, without the inclusion of thermal fluid distributions. We utilize Griffin and Pronghorn for generating global and local PKE parameters which includes the use of thermal fluid distributions to account for localized effects. This work establishes a methodology for calculating PKE parameters for a pebble bed reactor with a coupled neutronics/thermal fluids analysis. PKE parameters generated on a global and local basis are compared against a diffusion solve for a typical load-following transient. Locally-generated neutron kinetic parameters are able to reduce the maximum error in the transient power level from 5% to below 1.5%; along with this, bulk temperature errors were reduced from 12 K to 4 K.
This work applies Lie Group Theory to the nonlinear point-reactor neutron kinetics equations for several reactivity feedback and reactivity insertion models. The goal of this work is to connect known analytical solutions to their corresponding Lie groups as well as to obtain new analytical solutions and the circumstances from which they arise. The reactivity models we consider are the Nordheim-Fuchs Model, an arbitrary-order polynomial reactivity insertion model, and the Fuchs Ramp-Insertion Model. In all three models, we generalize the analysis when possible and obtain analytical solutions in every case. For the first two models, we succeed in finding the Lie groups corresponding to classical solutions and find additional solutions that were otherwise unknown. Finally, for the final model, we obtain a reduced-order nonlinear kinetics equation which is not analytically solvable, but by applying a second-order truncation, we find a novel closed-form solution and demonstrate its validity with a numerical example.
Ketonization of volatile fatty acids (VFAs) produced via arrested methanogenesis of wet waste represents the only unit operation of wet waste upgrading to sustainable (net-zero or negative life cycle CO2 emissions) aviation fuel not currently in industrial practice. Ketone product yields of close to 100% have been obtained during gas-phase reactions of VFAs over oxide catalysts at the laboratory scale, but design of an industrial ketonization reactor requires understanding of the impacts of reactant and product partial pressures, deactivation, and heat and mass transport phenomena on observed ketonization rates. Our work leverages rigorous kinetic analysis of ketonization processes to inform reactor scale-up efforts through packed-bed reactor modeling. We first present results of a kinetic study of ketonization of a model VFA, hexanoic acid, over an industrial ZrO2 catalyst performed in a packed-bed microreactor in conditions free of significant heat or mass transfer gradients. Major findings of the analysis include: (i) hexanoic acid saturates all catalyst active sites at relatively low partial pressure (~10 kPa) and (ii) ketonization products 6-undecanone, water, and CO2 inhibit reaction rates. Kinetic data are used to fit a rate expression quantifying the functional dependence of ketonization rate on partial pressures of VFA reactants and ketone, water, and CO2 products. A packed-bed reactor model describing vapor-phase hexanoic acid ketonization over ZrO2 extrudate catalysts is developed based on the kinetic model. The effects of (i) bed- and pellet-scale mass- and heat-transfer limitations and (ii) axial pressure drop guide our development of recommendations for optimal sizing, temperature, and influent flow composition of an industrial-scale ketonization reactor. The quantitative understanding of VFA ketonization developed in this study represents an advance toward derisking the VFA ketonization step of wet waste upgrading to sustainable aviation fuel.
Delayed neutron precursor (DNP) group data are important for modeling reactor dynamics. Although the data for individual DNPs have been developed over time, the DNP group data present in the Evaluated Nuclear Data Files (ENDF) have not been updated in the past 20 years. In this work, we use SCALE to recreate the Godiva experiment that was used to generate the original DNP group structure for fast fission of 235U. However, each DNP is modeled using up-to-date data, and the results are then converted into a newly updated group structure. This conversion uses an iterative linear least squares solver to minimize chi-squared. The approaches used in this work also enable energy spectrum generation and uncertainty tracking. The method used in this paper for fast 235U fission DNP group structure updating can be applied to different energies and fissile nuclides. Demonstration of the uncertainty tracking in reactor kinetics and dynamics simulations is shown using point reactor kinetics simulations. Results show that there are data discrepancies between the International Atomic Energy Agency database and data used in ORIGEN, which are currently being fixed. Results also show that the proposed method for group spectra generation performs well.
Molten Salt Reactors (MSRs) offer significant versatility supported by their design flexibility and safety advantages due to high temperatures, design customization capabilities, strong negative reactivity feedback effects, and salt chemistry characteristics. Different control schedules and different mass flow rates in the MSR primary systems impact neutronics and thermal hydraulics coupling. The result is two-fold. First, it leads to adaptability in various applications when the customization of a MSR unit is desired. Second, it calls for the quantification of dynamics and safety characteristics to assess the operational domain and make sure reactor stability and inherent safety characteristics are maintained. This paper provides an in-depth transient analysis of the operational implications of MSRs having heat-generating fuel and heat-transporting salt mixed forming liquid fuel salts contained within primary systems. A coupling model was developed to simulate fuel flow in molten salt reactors. A zero-dimensional reactor kinetics model was used with a one-dimensional heat transfer model in the reactor core. The temperature reactivity feedbacks of the fuel salt resulting from the thermal-hydraulics model were used in the reactor kinetics model to complete the coupled code. The model was applied to simulate the Molten Salt Breeder Reactor (MSBR) response at steady state and transients. The developed model was shown to be a suitable tool for the dynamic analysis of molten salt reactors.
Not Available
Department of Energy National Nuclear Security Administration Office of Material Management and Minimization’s mission includes the conversion of civilian research reactors from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel. Analyses have shown that the Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) will need to operate at 95 MW for the LEU silicide dispersion fuel designs to match key performance metrics obtained with HEU fuel at 85 MW. To prove safe operation of HFIR after installation of plant modifications to increase power, a high power HEU test cycle was proposed. Neutronics model updates and reactor physics analyses are performed to support the development of safety design reports for the high power (HP) HEU test cycle. Reactor physics metrics evaluated herein include fuel depletion, actinide production, cycle length, fission rate density distributions, axial power peaking factors, and reactor kinetic parameters. These reactor physics analyses support the development of future LEU safety design reports by providing key input for future HP HEU HFIR thermal hydraulics and reactor transient safety analyses.
The equivalent fundamental-mode source refers to a “source” which is identically distributed in space, energy, and angle as that of a fundamental-mode fission source distribution. It has previously been shown that the system fixed source multiplication ($M_{fs}$) relates to the effective multiplication factor ($\kappa_{eff}$) as $M_{fs} = \frac{g^*}{1-k_{eff}}$ The term $\mathcal{g}$* therefore relates the system source distribution to the system equivalent fundamental-mode source. As previously shown, the $\mathcal{g}$* term is approximately 1 for deeply subcritical systems and then diverges as the system multiplication increases. For a system with a point source located in the center of an assembly (such a sphere), $\mathcal{g}$* will be greater than 1; this happens because there is a higher probability of induced fission when all of the starter neutrons originate from the center of the sphere. For a uniform source, $\mathcal{g}$* is less than 1, but does not diverge from 1 as dramatically as seen by a point source. One reason that $\mathcal{g}$* is an important parameter is that it is used in measurement methods to infer reactor kinetics parameters ($β_{eff}$ being one example). This work will discuss multiple ways to simulate $\mathcal{g}$* using MCNP®6.2, which include use of traditional methods versus a single input file method that was recently published. A discussion related to methods to simulate $M_{fs}$ will be presented. Spherical systems of Highly Enriched Uranium (HEU) will be compared with previous works. In addition, results of plutonium including multiple reflector materials will be investigated. Simulated systems will include both parametric studies as well as configurations used in recent experiments performed at the National Criticality Experiments Research Center (NCERC). These studies will include subtleties associated with system geometry that have not been previously shown. Last, this work will investigate the relationship between $\mathcal{g}$* and reactor kinetics parameters (such as $β_{eff}$).
The National Criticality Experiments Research Center (NCERC) is a general purpose criticality experiments facility located inside the Device Assembly Facility (DAF) at the Nevada National Security Site (NNSS). Critical experiments containing any special nuclear material, any enrich ment/separation, most physical forms, and any configuration are possible within the constraints of the defined safety basis. NCERC draws upon physical assets and experimental knowledge to solve some of the most difficult problems with respect to criticality safety, reactor physics, and reactor kinetics. In terms of physical assets, NCERC houses hundreds of kilograms of special nuclear material with a majority consisting of highly enriched uranium (HEU) and weapons grade plutonium (WGPu). NCERC is home to four critical assembly machines: Comet, Planet, Flattop, and Godiva IV. To support various derivative diagnostics on fissioning systems, NCERC houses a count room to measure irradiated samples and dosimeters. This paper will step through each of these capabilities explaining recently completed work and upgrades.
The presentation discusses benchmarking in comparison to experimental truth, the integral experiments that were performed, and that benchmarks are evaluated experiments. It also discusses validation testing, and that the ultimate goal of the validation end product is to improve evaluated nuclear data for applications It does state that validation highlights errors in the Nuclear Data Pipeline that affect applications, for example: Missing Cd Capture Gammas in ENDF/B-VII.1. It acknowledges that additional types of experiments are needed to test data used in applications, and that validation experiments do not have to be complicated and expensive. The presentation also discusses different types of measurements including activation foil and fission chamber measurements, reactor kinetics measurements, and subcritical measurements. The presentation concludes with an experiment summary and states that all applications need validation.