Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactor Kinetics”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Putting multipoint kinetics to the test: preliminary analysis of the MUSE-4 experiment with APOLLO3

Modeling neutron kinetics is a challenging task and over the years, several approaches have been developed offering different tradeoffs between precision and cost. Probably the most popular one is point kinetics (PK), which, despite being very cheap, provides good results for many applications. One of the limitations of PK is that it is not able to model effects that happen before the prompt jump time scale. In practice, applications, where those time scales are interesting, are rare and those for which experimental data are available are even rarer. One of those cases is the MUSE-4 experiment for which the focus is on the very fast behavior and where we have no thermal feedback. As we have shown in a recent paper, a simple extension of PK, known as multipoint kinetics (MPK), can overcome this limitation. In this paper, which is intended as a preliminary study, we show how a simple MPK model can qualitatively reproduce the fast behavior that we cannot model with PK. Moreover, we will show how to exploit its low cost and compact form to model also slow effects as the precursors built up, behaviours that are not affordable to model with more expensive techniques, such as kinetic Monte Carlo simulations. For the periodic pulse transients with no feedbacks we will also provide an analytical solution in closed form. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Quantifying the impact of temporal analysis of products reactor initial state uncertainties on kinetic parameters

The temporal analysis of products (TAP) reactor, a transient kinetic tool, provides users with information as the catalyst state evolves. However, the state of the art for TAP uncertainty quantification only considers the experimental noise present in the outlet flow signal. Additional sources of uncertainty, including, but not limited to, initial surface coverages, the catalyst zone location, the inert void fraction, and the gas pulse intensity and pulse delay, are not included. For this reason, a framework for quantifying all uncertainty sources present in TAP experiments is presented and applied to a carbon monoxide oxidation case study. Herein, two methods for quantifying these sources of uncertainty are introduced. The first utilizes initial state sensitivities to approximate the parameter variances, as well as to provide insights into the structural certainty of the model. The second generates parameter confidence distributions through an ensemble-based sampling algorithm. This initial state covariance matrix can ultimately be merged with the experimental noise covariance matrix, providing a unified description of the parameter uncertainties for a single TAP experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Inferring small-sample-reactivity worth from oscillation experiment with auto-corrected local flux

To estimate reactivity changes as a function of time from the signal of detectors count rates, one usually uses point kinetic equations. However, when detectors are located inside the core, it is often necessary to compute a correction factor to account for the flux-shape deformation (i.e. the change of detector efficiency that may have been induced by neutron absorbers movements). Various simulation methods have been developed depending on transients observed. They range from the resolution of static equations up to full-time-dependent neutron transport in 3D. However, those methods require high computational-cost calculations of local reaction rates that rely on nuclear data evaluations. Therefore, a 'measured' reactivity is never determined through pure measurements. Taking advantage of the known oscillation movement of a small sample at the center of the MINERVE reactor, the innovative algorithm presented in this paper derives the correction factor needed to get no time drift on reactivity steps. Derived correction factors are in good agreement with those obtained independently with TRIPOLI-4 through sensitivity analysis. Low sensitivity to the sample composition and moderate sensitivity to the detectors position reinforce the reliability of those results. This application paves the way to almost pure nuclear data reactivity measurements. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Machine Learning-Assisted Recovery of Delicate Kinetic Information from Transient Reactor Experiments

Identifying active sites and their roles in chemical reaction steps remains a vital challenge in heterogeneous catalysis. Transient experiments offer a unique way to probe active sites and distinguish subtle kinetic features. Although physics-based analysis methods may be well-developed, they can be highly susceptible to experimental noise, and smoothing methods may erase or even distort important features; a smooth curve is not always the best curve. We demonstrate a new workflow for the direct interpretation of intrinsic kinetic information from exit flux curves measured in transient reactor experiments. This workflow contains three artificial neural networks (ANNs), including a noise reducer, a concentration predictor, and a rate predictor to analyze experimental data, followed by the virtual TAP (VTAP) physics-based reactor model and density functional theory (DFT) calculations of adsorption energies on specific sites. We use this workflow to analyze the data from experiments titrating Pt/Al 2 O 3 and Pt/SiO 2 catalysts with carbon monoxide (CO) in the temporal analysis of products (TAP) reactor. Our workflow separates the time-evolving chemical reaction and mass transfer information contained in the TAP pulse response. The existence of strong- and weak-binding sites on the Pt/Al 2 O 3 catalyst is observed in the catalyst titration experiment in the transient reactor. The structures of the strong- and weak-binding sites are then identified by using DFT calculations. We find that the Pt/SiO 2 catalyst has only strong-binding sites, which aligns with the inactive support effect of SiO 2 . We demonstrate how machine learning methods provide unique insights with high-resolution data analysis that cannot be achieved by using state-of-the-art physics-based methods.

Adsorption↗
DOE OSTI2026

Mechanistic Insights for Plasma-Catalytic CO 2 Reduction over TiO 2 in a Dielectric Barrier Discharge Reactor

Reaction kinetics experiments coupled with phenomenological kinetic modeling and parameter estimation are used to elicit insights into the mechanism and active sites for the plasma-catalytic dissociation of CO 2 on TiO 2 . Experimental and model insights showed that gas-phase reactions contribute at least two-thirds of the overall product formation at explored conditions; weak temperature dependence, strong sensitivity to specific energy input (SEI), apparent first order in CO 2 , and positive influence of cofed argon (Ar) and oxygen (O 2 ) for the gas-phase contributions all suggest that expected plasma reaction steps such as electron-impact and high-energy collisions are the dominant modes for CO 2 dissociation. The Arrhenius-like expression for gas contributions resulted in a preexponential of 4.40 × 10 –3 s –1 , an E SEI,g of 7.90 × 10 –4 mol/kJ, and an E a,g of 1.00 × 10 –3 J/mol. For surface contributions, the small apparent barrier of 16.3 kJ/mol, relatively weaker dependence on SEI, first-order dependence on CO 2 , and insensitivity to cofed Ar and O 2 all point to CO 2 dissociation on TiO 2 surface facets without vacancies and aided by plasma (leading to vibrationally excited CO 2 and/or a reactive surface with significant surface charge accumulation). The Arrhenius-like expression resulted in a preexponential of 7.81 × 10 –2 s –1 , an E SEI,s of 1.90 × 10 –3 mol/kJ, and an E a,s of 1.63 × 10 4 J/mol. The derived kinetic model further enabled a systematic evaluation of the effect of inputs (plasma power, flow rate, CO 2 inlet concentration, and temperature) to identify process trends and optimal operating conditions.

catalyst↗
DOE OSTI2021

Versatile semi-batch apparatus for manometric measurement of gas–solid reaction rates

In this study, we describe the design, construction, operation, and performance of a simple and versatile semi-batch reactor that is especially useful for measurement of gas/solid reaction rates at pressures in the range of 1 mTorr to 1500 Torr. The reactor operates by repeatedly imposing small AC modulations of reactant gas pressure on top of a much larger DC pressure background. Based on the rate of pressure relaxation following each AC pulse, the reaction rate is determined. Our design is characterized by modular construction from off-the-shelf, ultra-high-vacuum-compatible components, which facilitate easy retrofitting and adaptation to a range of experimental conditions. Automated experiment control and data acquisition is accomplished via a custom National Instruments © LabView virtual instrument. Data analysis is automated using a custom series of Mathworks © Matlab scripts. We demonstrate reactor performance through measurements of hydrogenation kinetics for a composite H 2 getter material consisting of 1,4-bis(phenylethynyl)benzene mixed with a palladium/carbon catalyst.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Unstructured mesh-based multi-physics reactor transient simulation with iMC code

In this study, we present a time-dependent Monte Carlo simulation result about a complex fuel-loaded light water reactor system with thermal-hydraulic feedback using the KAIST iMC code. The intra-pin temperature distribution is calculated based on Monte Carlo-tallied detailed power distribution and the finite element heat transfer method. The temperature effect of fuel, absorber, and coolant are adequately reflected into the dynamic response by the on-the-fly Doppler broadened cross-section correction. For the time-dependent Monte Carlo transport simulation, the predictor-corrector quasi-static Monte Carlo method is used for this study. We tested the coupled analysis framework with a 17-by-17 fuel assembly loaded with the centrally-shielded burnable absorber (CSBA) fuel element. The temperature feedback effect on the material cross-section and the moderator density incurred negative reactivity in response to the external positive reactivity insertion, showing the inherent stability of the reactor system. The pin-wise power and temperature change during the system transient are presented and discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Stability margin analysis of the Holos-quad microreactor design

A simplified point model of the Holos-Quad microreactor is introduced. The model is based on the point kinetics equations coupled to three heat balance equations, representing the mean temperatures of the fuel particles, the graphite moderator and the helium coolant. The differential equations are converted to the frequency domain, enabling the construction of the closed-loop reactor transfer function. Using this function the stability margins of the core design is analyzed for various power levels. It is shown that the gain margins approaches infinity, demonstrating the stability of the reactor for all power levels. The phase margin at nominal power is about 60 degrees, however it shows a non-monotonous dependence on power, with minimal value obtained for about 10% of nominal power. As power level is further increased, the phase margin also increases, demonstrating the reactor becomes more stable. This behavior may be of high importance in load-follow scenarios, where the power level of the reactor changes with time. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Review and Summary of Oxide Thickness Data for Aluminum-Clad Spent Nuclear Fuel

The U.S. Department of Energy (DOE) owns a large inventory of aluminum-clad spent nuclear fuel (ASNF) in interim storage pending ultimate disposition, with more being generated by currently operating research reactors. Dry storage in sealed DOE Standard Canisters is being investigated as an approach for long-term interim storage and/or disposition in a repository for ASNF. A primary challenge for ASNF storage is the presence of aluminum (oxy)hydroxide layers formed on the cladding surfaces during water exposure in the reactor and in wet storage, which forms a reservoir of chemisorbed water not readily removed at low (<100°C) drying temperatures. Free, physisorbed, and chemisorbed water are all susceptible to radiolytic breakdown under irradiation and can release hydrogen gas. Identifying the likely range of (oxy)hydroxide loadings on ASNF that may be placed in dry storage will help to ensure that the impact of the (oxy)hydroxide is adequately accounted for while avoiding over conservatism and enable mitigation strategies to be implemented where needed. The current report summarizes information on (oxy)hydroxide thicknesses and characteristics from the literature and from recent measurements taken in the present research campaign. In general, corrosion studies have indicated that the corrosion and oxide buildup on aluminum are affected by numerous conditions, including pH, temperature and heat flux, coolant flowrates, irradiation (in-reactor vs. unirradiated tests), duration of water exposure, and the amount of oxide already on the surface. Some of these factors are interrelated (e.g., the local temperature, heat flux, and coolant flowrate). Observations of the impact of irradiation (in-reactor versus ex-reactor measurements) on corrosion and oxide thickness are mixed, but multiple studies indicate that there is a significant difference between in-reactor and ex-reactor corrosion kinetics, even when other operating conditions such as pH, flowrate, heat flux, etc., simulate those in a reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Physics-Informed Neural Network Solution of Point Kinetics Equations for a Nuclear Reactor Digital Twin

A digital twin (DT) for nuclear reactor monitoring can be implemented using either a differential equations-based physics model or a data-driven machine learning model. The challenge of a physics-model-based DT consists of achieving sufficient model fidelity to represent a complex experimental system, whereas the challenge of a data-driven DT consists of extensive training requirements and a potential lack of predictive ability. We investigate the performance of a hybrid approach, which is based on physics-informed neural networks (PINNs) that encode fundamental physical laws into the loss function of the neural network. We develop a PINN model to solve the point kinetic equations (PKEs), which are time-dependent, stiff, nonlinear, ordinary differential equations that constitute a nuclear reactor reduced-order model under the approximation of ignoring spatial dependence of the neutron flux. The PINN model solution of PKEs is developed to monitor the start-up transient of Purdue University Reactor Number One (PUR-1) using experimental parameters for the reactivity feedback schedule and the neutron source. The results demonstrate strong agreement between the PINN solution and finite difference numerical solution of PKEs. We investigate PINNs performance in both data interpolation and extrapolation. For the test cases considered, the extrapolation errors are comparable to those of interpolation predictions. Extrapolation accuracy decreases with increasing time interval.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Development and benchmarking of transient nodal code SIMULATE5-K neutron kinetics solver

SIMULATE5-K is Studsvik's next generation best estimate transient code. The time dependent diffusion equation is solved with a nodal method consistent with that implemented in the licensed core design code SIMULATE5. Arbitrary number of neutron and delayed neutron precursor groups can be used. For the solution of the spatial problem, the coupling coefficients used to relate the node leakages are found by first converting the time dependent diffusion equation to a static diffusion equation with the use of flux and delayed neutron precursor dynamic frequencies. Once the static-like equations are obtained, the multi-group analytical nodal model is used to obtain the coupling coefficients, expressing the node leakage in terms of adjacent node average fluxes. The coupling coefficients are then inserted into the time dependent nodal balance equation. For the time integration, the time dependent neutron balance equation is solved with the frequency transformation method. The treatment of the temporal dependence yields a fixed source problem which can be solved utilizing the existing fixed-source methodology. The primary purpose of this paper is to describe the neutron kinetics methodology implemented in SIMULATE5-K. The accuracy of the method is demonstrated for a series of well-known, neutronic-only benchmark problems. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Utilizing Physics-Informed Synthetic Data to Train a Digital Twin for Predicting Reactor Operations

Understanding techniques to strengthen the nuclear nonproliferation regime is crucial in reducing the creation of nuclear weaponry on the basis of advancements in the nuclear energy industry. Prior to construction of a nuclear power plant, it is necessary to understand the proliferation potential of the plant’s reactor. Digital twins serve as a unique solution to recognizing reactor behavior indicative of nuclear proliferation. The following research conducted serves as a validation to training a digital twin on synthetic data fabricated via means of reactor physics simulations based on parameters of Idaho State University’s AGN-201 reactor. The synthetic data is utilized to train a long short-term memory (LSTM) recurrent neural network model. The accuracy of the predicted data is measured against real operational data to verify the reliability of the synthetic data creation methods and if these methods should be used in the future to inform inspectors of a reactor’s proliferation capabilities.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2021

Fuel Properties and Kinetics

The Co-Optimization of Fuels and Engines (CoOptima) program is focused on a two-pronged approach to develop new, lower net-carbon fuels and advanced engine technologies in tandem. This effort brings nine national labs and several university and industry partners together to simultaneously identify new fuel candidates and to further engine technology development in the light-, medium-, and heavy-duty vehicle sectors with a focus on fuel economy and emission reduction. The fuel properties and kinetics project is aimed at achieving a better understanding of fuel property kinetics and how they impact engine performance such that more efficient and lower emitting engines can be made a reality. Several promising biofuel candidates were identified through tiered screening approaches and were evaluated utilizing unique bench-scale measurement techniques. Kinetic models were developed to describe autoignition and soot pre-cursor formation mechanisms to identify the most promising candidates that advanced various combustion strategies.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Moving from Fundamental Knowledge of Kinetics and Mechanisms on Surfaces to Prediction of Catalyst Performance in Reactors

Catalysis is defined by kinetics, but after over a century of catalyst development our kinetic and mechanistic knowledge governing real-world catalysis is still severely lacking. This Perspective considers how we can use the precise knowledge available from both theoretical and fundamental experimental techniques to understand and even predict catalytic performance under operational conditions. Here, we describe advances that link fundamental and "real-world" measurements, crucial elements required to successfully recreate a reaction network, and avenues which appear to offer viable routes to success.

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

High-order accurate solutions of the point kinetics equations with the spectral deferred correction method

Solving initial value problems with higher-order methods can improve the accuracy of the simulation results or the efficiency of the calculation. In this paper, we apply the spectral deferred correction (SDC) method to solve the initial value problem of the point kinetics equations (PKE). SDC is a stable, robust, and efficient high-order time-integration method capable of an arbitrary order of accuracy. For our implementation we show that it is A-stable for orders up to 8 and the order of accuracy is verified for PKE problems with a range of different reactivities. A 5.-order SDC method was then implemented to solve the exact PKE (EPKE) in the Transient Multilevel (TML) method of MPACT. The error from solutions of the EPKE is shown to be negligible. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

The impact of the third O 2 addition reaction network on ignition delay times of neo-pentane

In this work, we studied the oxidation of neo-pentane by combining experiments, theoretical calculations, and mechanistic developments to elucidate the impact of the 3rd O 2 addition reaction network on ignition delay time predictions. The experiments are based on photoionization mass spectrometry in jet-stirred and time-resolved flow reactors allowing for sensitive detection of the keto-hydroperoxide (KHP) and keto-dihydroperoxide (KDHP) intermediates. With neo-pentane exhibiting a unique symmetric molecular structure, which consequently results only in single KHP and KDHP isomers, theoretical calculations of ionization and fragment appearance energies and of absolute photoionization cross sections enabled the unambiguous identification and quantification of the KHP intermediate. Its temperature and time-resolved profiles together with calculated and experimentally observed KHP-to-KDHP signal ratios were compared to simulation results based on a newly developed mechanism that describes the 3rd O 2 addition reaction network. A satisfactory agreement has been observed between the experimental data points and the simulation results, thus adding confidence to the model's overall performance. Finally, this mechanism was used to predict ignition delay times reported previously in shock tube and rapid compression machine experiments (J. Bugler et al., Combust. Flame 163 (2016) 138–156). While the model accurately reproduces the experimental data, simulations with and without the 3rd O 2 addition reaction network included reveal only a negligible effect on the predicted ignition delay times at 10 and 20 atm. According to model calculations, low temperatures and high pressures promote the importance of the 3rd O2 addition reactions.

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

Rapid depletion analysis of flowing-pebble reactor systems at equilibrium using SCALE

Several high-temperature gas-cooled reactor concepts (and more recently, salt-cooled designs such as the fluoride salt-cooled high-temperature reactor) feature core designs employing continuously circulating fuel pebbles. These reactor designs permit both continuous online refueling of fuel elements as well as higher overall achievable discharge burnups. However, rapid calculation of time-dependent fuel isotopic inventories proves challenging for this class of dynamic systems with current analysis tools. While iterative approaches employing coupled neutron transport have been developed to solve this issue, rapid depletion analysis techniques are needed to calculate time-dependent inventories for individual pebbles and batches (and thus the construction of full- core inventory at equilibrium). We propose a depletion analysis strategy for this type of system for cores at equilibrium. Drawing upon previous neutronic analysis of the PBMR-400 equilibrium core, we demonstrate the viability of developing collapsed one-group cross section libraries suitable for performing rapid depletion analyses with SCALE. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗