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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Nuclear Thermal Rocket Emulator for a Hardware-in-the-Loop Test Bed

To support NASA’s mission to use nuclear thermal rockets for future Mars missions, an instrumentation and control test bed has been built at Oak Ridge National Laboratory. The system is designed as a hardware-in-the-loop test bed for testing control elements and autonomous control algorithms for nuclear thermal propulsion rockets. The mock reactor system consists of a modular and scalable framework, using inexpensive components and open-source software. The hardware system consists of a two-phase flow loop and a mock reactor with six control drums. A single-board computer (NVIDIA Jetson) handles reactor core emulation and hosts a message queuing telemetry transport broker that allows user-deployed control algorithms to interact with the system hardware. The reactor emulator receives sensor data from the hardware and provides the simulated performance of the reactor under steady-state, transient, and fault conditions. The emulator uses a reactivity lookup table and the point kinetics equations to solve for the reactor dynamics in real time. Emulated reactor dynamics and sensor input inform the autonomous control algorithm’s decision-making in a closed-loop manner. The current system is capable of operating at 10 Hz, but faster cycle rates are an area of ongoing research. This test bed will enable NASA and other space vendors to rigorously test their autonomous control systems for NTP rockets under transient (reactor startup and shutdown), steady-state, and fault conditions to reduce development time and risk for autonomous control systems in future missions.

autonomous control↗
DOE OSTI2023

A Predictive Transient Model of the TREAT-SIRIUS Experiments

To support experiments modeling and transient analysis of the NASA-sponsored SIRIUS experiments in the Transient Reactor Test Facility (TREAT), an innovative predictive transient model has been developed to simulate various transient experiments. The developed model is based on a new approach that utilizes steady-state Monte Carlo calculations along with a surrogate model based on polynomial regression to determine the core reactivity. The reactivity is supplied to an exact point kinetics model that determines the reactor power, and it is coupled to an adiabatic feedback model to determine the fuel average temperature. In this work, the newly developed model is introduced along with validation test results considering SIRIUS-1 experiments at different power levels. The initial test results show a very good agreement with experimental data at different power levels, which indicates that the model can be used to predict the power of other transient tests.

42 ENGINEERING↗
DOE OSTI2025

Progress in modeling hydrogen assisted ammonia oxidation with new experiments and a further reconciliation of the NH 3 + OH rate constant

Hydrogen-assisted oxidation of ammonia in a premixed, laminar flow tubular reactor under reducing conditions was investigated experimentally and through chemical kinetic modeling. Due to its impact on the competition for OH among ammonia and hydrogen, the rate constant for NH 3 + OH (R1) was determined through state-of-the-art theoretical kinetics calculations, employing composite energies that include the effects of higher order electronic excitations on the electronic energies along the variational reaction path and treating the limitations in the kinetics posed by the passage through a hydrogen-bonded complex. The resulting rate constant was in close agreement with the recent experimental value from Zaczek et al. (2025), settling a long-term dispute about the high-temperature value of k 1 and confirming within 20% the value previously used in modeling. The chemical kinetic model, with no other changes, captured well measured concentrations of NH 3 , H 2 , NO, and N 2 O from flow reactor oxidation of NH 3 /H 2 at slightly reducing conditions over a range of temperature (900-1350 K) and NH 3 /H 2 ratios (0.5-2.0). Comparison of the present results with reported data from a non-premixed setup indicates that for laminar flow tubular reactors, the reactor configuration may have implications for the observed H 2 consumption due to the possibility of preferential oxidation during mixing.

Ab initio theory↗
DOE OSTI2023

Non-equilibrium low-temperature plasma-assisted combustion of iso-octane: Perturbing pyrolysis and oxidation kinetics

Here, in this study, a plasma-coupled flow reactor facility is used to examine the effects of non-equilibrium low-temperature plasmas on perturbing the pyrolysis and oxidation kinetics of iso-octane. Experiments were performed in highly dilute reactive mixtures of nitrogen, at near isothermal conditions for temperatures ranging from 523 K to 1203 K. Experiments cumulatively demonstrated enhanced chemical reactivity with the plasma for temperatures below 900 K, and a lowering of the hot-ignition temperature. Detailed kinetic insight was derived from a 0D plasma-coupled kinetic model, utilizing a constructed mechanism that combined both plasma-specific chemistry and the neutral combustion chemistry. For pyrolysis conditions, the model displayed relatively good agreement with fuel consumption and the formation of most intermediates compared to the experimental data, demonstrating the model is able to accurately predict primary radical formation from the plasma directly interacting with the fuel. Enhanced reactivity was attributed to collisional quenching of excited-states of N 2 with fuel, which led to efficient fuel fragmentation and enhancement of the H-radical flux. For oxidation conditions, the model displayed satisfactory agreement with the experiments. Model predictions were able to accurately predict fuel consumption and most intermediate speciation data for T > 800 K, but most discrepancies were towards T < 800 K in particular with oxygenated intermediates. In the presence of oxygen, plasma effects were predominantly spent on efficient enhancement of O- and H-radical fluxes, leading to further fuel fragmentation and initiation of both the OH- and HO 2 -radical pools. Subsequent reactivity of iso-octane was then dictated by the response of the temperature-dependent neutral chemistry. At low-temperatures (T = 643 K), enhanced fuel radicals and O 2 -additon chemistry lead to the formation of oxygenated species, while at intermediate temperatures (T = 843 K) net decrease in OH-radical reactivity led to an increase in hydrocarbon speciation. Near the self-ignition threshold (T = 1163 K), radicals generated by high-temperature branching reactions dominate the oxidation process and effectively ignition. This study ultimately demonstrated that the enhancement of radicals afforded by the plasma causes a deviation in known understanding of iso-octane kinetics in some regards and warrants future studies to reconcile these discrepancies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

A kinetic line-driven radiation operator and its application to Gyrokinetics

A velocity dependent, kinetic model for line radiation is developed for continuum kinetic codes. It has been implemented in the full-f gyrokinetic code Gkeyll. The total radiation for a charge state is modeled as an advection in velocity space with a form of $\nabla_v \cdot(v\nu(v)f(v))$, guaranteeing particle conservation. The velocity dependence (in the form of an effective frequency $\nu(v)$) is found through fitting the energy loss of the operator, i.e. the second velocity moment, to the radiation data in the OpenADAS database. Therefore, each individual transition does not need to be evaluated every time step, significantly reducing the computational cost of including line radiation in a kinetic model. The dependence on velocity instead of the usual, temperature, allows the radiation to be computed from non-Maxwellian electron distribution functions: We benchmark the model against a collisional radiative model using isotropic non-Maxwellian distribution functions. A velocity dependent model of radiation can more accurately describe the radiation in the more kinetic regimes expected in reactor-scale devices. The velocity dependence qualitatively captures the quantum mechanical need for a minimum velocity before any radiation occurs.

kinetic↗
DOE OSTI2022

Stability analysis of a one-dimensional multiphysics model of a molten salt fast reactor

Reactor designs with a fast neutron spectrum have recently moved into the research focus for possible Generation IV reactor concepts - a representative of which is the molten salt fast reactor (MSFR). In order to assess the safety and reliability of this reactor concept, it is necessary to obtain knowledge of the system behaviour. An important part of this research concerns the dynamical stability of the MSFR. Research regarding the MSFR stability so far has either been using linearized models for stability criteria or only been examining certain transients at fixed parameter values. This work delivers a comprehensive stability analysis for a non-linearized MSFR model, and a wide range of values for all relevant parameters. For this purpose, a one-dimensional MSFR model was set up, taking neutron kinetics and thermal hydraulics in the form of a system of differential equations into account. The stability of this model was investigated by means of the numerical tool MATCONT, which was used to monitor the evolution of a so-called fixed-point solution, here referring to the steady state at operating conditions of the system, while system parameters got varied. The MATCONT analyses showed no loss of stability for any of the considered parameter variations and no solutions that might co-exist in parallel with the steady-state fixed-point solution. Therefore these results indicate a stable fixed point to which all solution transients converge, making the considered MSFR model insusceptible to deviations from the equilibrium. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Investigations into Plasma-Mediated Decomposition of Organoiodide Species as a Pretreatment for Mitigation of Radioiodine Emissions

A plasma-mediated pretreatment step is proposed to decompose methyl iodide (CH 3 I) for facilitating the capture of the iodine radioisotopes present in off-gas streams from used nuclear fuel reprocessing operations. Simulated CH 3 I gas streams were exposed to dielectric barrier discharge (DBD) plasma in a quartz glass continuous-flow reactor using copper electrodes. The kinetics of decomposition of CH 3 I in nitrogen was investigated by varying the applied voltage, inlet concentration, and residence time. The rate of decomposition was proportional to CH 3 I concentration and increased as the applied potential was increased. This effect was incorporated in the kinetics through the dependence of the rate constant on the applied electric field. Decomposition of CH 3 I was not affected by moisture, and the DBD pretreatment is likely to be similarly effective in an air environment as well; however, the results were confounded by a plasma-mediated reaction between nitrogen and oxygen.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Kinetic study of reaction C 2 H 5 + HO 2 in a photolysis reactor with time-resolved Faraday rotation spectroscopy

The rate constant and branching ratios of ethyl reaction with hydroperoxyl radical, C 2 H 5 + HO 2 (1), a key radical-radical reaction for intermediate temperature combustion chemistry, were measured in situ for the first time in a photolysis Herriott cell by using mid-IR Faraday rotation spectroscopy (FRS) and UV-IR direct absorption spectroscopy (DAS). The microsecond time-resolved diagnostic technique in this work enabled the direct rate measurements of the target reaction at 40 and 80 mbar and reduced the experimental uncertainty considerably. C 2 H 5 and HO 2 radicals were generated by the photolysis of (COCl) 2 /C 2 H 5 I/CH 3 OH/O 2 /He mixture at 266 nm. By direct measurements of the transient profiles of C 2 H 5 , HO 2 and OH concentrations, the overall rate constant for this reaction at 297 K was determined as k 1 (40 mbar) = (3.8 ± 0.8) × 10 –11 cm 3 molecule –1 s –1 and k 1 (80 mbar) = (4.1 ± 1.0) × 10 –11 cm 3 molecule –1 s –1 . As a result, the direct observation of hydroxyl radical (OH) indicated that OH formation channel was the major channel with a branching ratio of 0.8 ± 0.1.

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

Machine learning-based ethylene and carbon monoxide estimation, real-time optimization, and multivariable feedback control of an experimental electrochemical reactor

Electrochemical reduction of CO 2 gas is a novel CO 2 utilization technique that has the potential to mitigate the global climate crisis caused by anthropogenic CO 2 emissions, and enable the large-scale storage of energy generated from renewable sources in the form of carbon-based chemicals and fuels. However, due to the complexity of the electrochemical reactions, the explicit first-principles models for CO2 reduction are not available yet, and there has been a limited effort to develop process modeling, optimization and control of CO 2 electrochemical reactors. To this end, a rotating cylinder electrode (RCE) reactor has been constructed at UCLA to understand the mass transfer and reaction kinetics effects separately on the productivity. In the RCE reactor, the applied potential strongly influences the reaction energetics and the electrode rotation speed affects the hydrodynamic boundary layer and modifies the film mass transfer coefficient, which involves convective and diffusive transport. Further, the present work aims to develop a multi-input multi-output (MIMO) control scheme for the RCE reactor that integrates techniques from artificial and recurrent neural network modeling, nonlinear optimization, and process controller design. Specifically, production rates of two products from the experimental reactor, ethylene and carbon monoxide, are controlled by manipulating two inputs, applied potential and catalyst rotation speed. Process dynamics and controllability are analyzed, a feedback control strategy is designed and the controllers are tuned accordingly. The experimental electrochemical cell is employed to gather data for process modeling and implement the multivariable control system. Finally, the experimental results are presented which demonstrate excellent closed-loop performance by the control system and regulation of the outputs at three different set-points including an economically-optimal set-point.

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

Kinetics of Hydrogen Desorption from Yttrium Hydride: The Development of Material Characterization Methods

Yttrium hydride possesses attractive neutronic and thermal properties for moderator applications as a high-temperature moderator material in advanced thermal neutron spectrum reactors that require small core volumes. For safe operation of the nuclear reactors, it is critical to understand the kinetics of H desorption from yttrium hydride at elevated temperatures. In previous studies, the H desorption flux as a function of temperature was measured by thermal desorption spectroscopy. To fully understand the H desorption kinetics, the amount of H within the specimen must be locally evaluated, especially if this uniformity strongly affects the thermomechanical properties and irradiation resistance of the bulk material.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Discerning the effect of various irradiation modes on the corrosion of Zircaloy-4

Using proton irradiation, this study investigates the individual influence of several factors on the corrosion kinetics of Zircaloy-4 in a hydrogenated water environment simulating a Pressurized Water Reactor (PWR). Using both simultaneous irradiation-corrosion and autoclave corrosion, we separately examine: (i) the effect of pre-irradiation on modifying the structure of the material, (ii) the impact of irradiation on creating defects in the growing oxide layer during corrosion, and (iii) the influence of irradiation on increasing the corrosion potential through radiolysis during corrosion. To replicate neutron-irradiated microstructure, two proton pre-irradiation schedules were employed: Schedule 1 (isothermal irradiation at 350 °C to 5 dpa) to simulate high-temperature PWR conditions, and Schedule 2 (two-step process: irradiation to 2.5 dpa at -10 °C followed by 2.5 dpa at 350 °C) to simulate lower temperature PWR and Boiling Water Reactor (BWR) conditions. Long-term autoclave corrosion testing over 360 days at 320 °C revealed no significant difference between unirradiated samples and those pre-irradiated according to either schedule, with all samples exhibiting sub-cubic kinetics within the pre-transition regime. Irradiated samples underwent Simultaneous Irradiation Corrosion (SIC) tests, corroding in 320 °C water while being irradiated with protons. Corrosion was found to accelerate in all SIC-tested samples relative to autoclave conditions, with the greatest increase observed in non-pre-irradiated samples. Pre-irradiation with either schedule resulted in a slower corrosion rate compared to non-pre-irradiated regions under SIC conditions. The degree of radiolysis observed in the SIC tests surpassed typical PWR conditions, approaching levels found in BWRs. Radiolysis products were identified as the primary contributors to accelerated corrosion, corroborated by radiolysis bar tests. Furthermore, these findings underscore the intricate interactions between irradiation, corrosion, and water chemistry in determining Zircaloy-4 corrosion kinetics within nuclear reactor environments.

Pressurized Water Reactor↗
DOE OSTI2021

Effective kinetics driven by dynamic concentration gradients under coupled transport and reaction

Biogeochemical reaction kinetics are generally established from batch reactors where concentrations are uniform. In natural systems, many biogeochemical processes are characterized by spatially and temporally variable concentration gradients that often occur at scales which are not resolved by field measurements or biogeochemical and reactive transport models. Yet, it is not clear how these sub-scale chemical gradients affect reaction kinetics compared to batch kinetics. Here we investigate this question by studying the paradigmatic case of localized pulses of solute reacting with a solid or a dissolved species in excess. Additionally, we consider non-linear biogeochemical reactions, representative of mineral dissolution, adsorption and redox reactions, which we quantify using simplified power-law kinetics. The combined effect of diffusion and reaction leads to effective kinetics that differ quantitatively and qualitatively from the batch kinetics. Depending on the nonlinearity (reaction order) of the local kinetics, these effects lead to either enhancement or decrease of the overall reaction rate, and result in a rich variety of reaction dynamics. We derive analytical results for the effective kinetics, which are validated by comparison to direct numerical simulations for a broad range of Damköhler numbers and reaction order. Our findings provide new insights into the interpretation of imperfectly mixed lab experiments, the effective kinetics of field systems characterized by intermittent reactant release and the integration of sub-scale concentration gradients in reactive transport models.

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

Low Temperature CO 2 Hydrogenation on Unsupported Mo 2 C Catalysts

CO 2 hydrogenation to methanol, a key reaction for decarbonizing the fuel and chemical industries, requires catalyst formulations that hydrogenate CO 2 selectively to methanol at temperatures where methanol conversion is not significantly equilibrium limited (<423 K). Herein we report continuous CO 2 hydrogenation at low temperatures (348-408 K, H 2 /CO 2 = 0.1-50, 5-35 bar) with high selectivity to methanol (up to ca. 80%) over unsupported β-Mo 2 C catalysts. Active site density quantification via titration with trifluoroacetic acid at reaction temperatures enables an assessment of site-specific rates. Methanation and reverse water gas shift (RWGS) occur concurrently with methanol synthesis during CO 2 hydrogenation over Mo 2 C. Reaction pathway analysis, product cofeeds, and reversibility formalisms show that all products form through primary reaction pathways from CO 2 , but secondary reactions of CO contribute significantly to rates of methanation. Dependences of forward rates on reactant and product concentration determined by independently varying the CO 2 , H 2 , CO, H 2 O, CH 3 OH, and CH 4 pressure in conjunction with reversibility formalisms reveal that all products form through H-assisted CO 2 activation and involve partially hydrogenated CO 2 -derived intermediates. Here, these inferences were verified by quantitative agreement between measured site-time yields and site-time yields predicted by closed form kinetic rate expressions in an integral reactor model over widely varying conditions (85-2000 kPa H 2 , 80-1500 kPa CO 2 , 0-45 kPa H 2 O, 0-21 kPa CO, 0-25 kPa CH 3 OH, 0-75 kPa CH 4 , 5-87 mol Mo s s mol CO 2 -1 ). Coverages calculated based on the kinetic model reveal that the Mo 2 C surface is covered with bidentate CO- and CO 2 -derived intermediates of the stoichiometry H 2 CO 2 and H 2 CO, indicating that H 2 and CO x do not compete for surface occupancy but instead adsorb cooperatively to form partially hydrogenated intermediates. Hydrogenation of the CO-derived H 2 CO** intermediate favors methanation, while hydrogenation of CO 2 -derived H 2 CO 2 ** favors methanol synthesis. Together, these findings demonstrate the ability of unsupported Mo 2 C to catalyze the hydrogenation of CO 2 to methanol at low temperatures and provide insight into the reaction network and mechanisms involved in its formation.

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

Simulating Catalysis with Realistic Pellet Geometries Using Mesoflow: A Case Study of Catalytic Propane Dehydrogenation

We present a case study of catalytic propane dehydrogenation with our open-source multiphysics solver, Mesoflow. The solver was developed to simulate reactive flow coupled to heterogeneous catalytic reactions and deactivation in the context of complex, mesoscale geometry. The method leverages cartesian block-structured adaptive mesh refinement to capture realistic catalyst microstructural features acquired directly from X-ray computed tomography data. A kinetic model for propane dehydrogenation and catalyst deactivation was developed based on temporal analysis of products (TAP) reactor experiments. The TAP reactor experiments allow for precise characterization of intrinsic kinetic reaction steps which are implemented into Mesoflow simulations to model the spatial and temporal evolution of reactants, products, and catalyst active sites. The short-term and long-term deactivation behavior is studied by using XCT data collected from fresh and aged catalyst pellets, which exhibit different microstructural features. This study employs time-splitting algorithms to connect disparate reaction and flow timescales, enabling the simulations to achieve realistic deactivation timescales on the order of minutes while the flow time-scales for small particles (100 microns) are several milliseconds. We also introduce a flexible automated python script that writes the necessary files to construct a Mesoflow simulation from user-created chemical mechanisms. We will also introduce a few new features that are added to Mesoflow such as higher order schemes, implicit chemistry integrators and the ability to run on AMD and NVIDIA graphics-processing-units.

AMReX↗
DOE OSTI2021

Advancing Catalytic Fast Pyrolysis Through Integrated Experimentation and Multi-Scale Computational Modeling

This webinar will highlight recent results from a multi-disciplinary research effort in which integrated reaction testing was coupled with particle- and reactor-scale computational modeling to advance catalytic fast pyrolysis (CFP) for the production of renewable hydrocarbon fuels. Data will be presented from a series of ex situ CFP experiments in which a fixed bed of Pt/TiO2 was utilized with co-fed H2 to upgrade woody biomass pyrolysis vapors. Further discussion will include the application of these data towards the development of (1) a multiscale simulation framework to de-couple apparent kinetics from both intraparticle and reactor-scale transport phenomena and (2) a finite element computational model to understand and predict thermal excursions during catalyst regeneration. Throughout the presentation, the speakers will emphasize synergistic outcomes derived from the collaborative approach and highlight ongoing research efforts to accelerate technology maturation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Electrochemical Control for Corrosion in Molten Chlorides During CSP Plant Operation

The Liquid Pathway of the Concentrating Solar Power Generation 3 (CSP Gen3) program proposed low-cost molten chloride salt for energy storage. However, online corrosion control was identified as a remining major risk of the Liquid Pathway approach. This project addressed that risk. Electrochemical solutions for corrosion mitigation during CSP plant operation were investigated and their feasibility and scalability were evaluated. The leading cause of corrosion in molten chloride salt systems was identified as corrosive impurities that form within the salt upon exposure to trace amounts of air and moisture. Leveraging electrochemistry, reduction/oxidation reactions can be employed to remove these corrosive impurities. In Phase 1 of this project, a bench-scale batch electrochemical reactor was designed, fabricated, and used to assess the kinetics and thermodynamics of electrochemical salt purification. In Phase 2, a laboratory-scale flow reactor was designed, fabricated, and used to assess the efficacy of the electrochemical method under flowing conditions. Results show that under proposed operating conditions for the Liquid Pathway Gen3 Pilot Plant, the electrochemical method is significantly more effective at removing impurities than alternative chemical and thermal methods, and that the electrochemical method produces less harmful byproducts. A key advance made in the course of this project was the development of a 2-electrode method for electrochemical purification that is more scalable than previously developed 3 electrode methods. This novel method is based on Magnesium (Mg) electrowinning. A provisional patent based on this invention has been submitted (USPTO Application No. 63/480,355). Additional key advances made during this project include assessment of the effect of dissimilar alloys on corrosion, kinetic and thermodynamic evaluation of thermolysis reactions of impurities within the molten salt, characterization of byproducts of purification reactions, and generation of IP focused on isolating value-added products using molten salt-based electrochemistry that could be deployed to valorize the process (USPTO Application No. 63/478,806). Ultimately, this project represented a step toward feasibility of Liquid Pathway Gen3 CSP. The method developed under this project could significantly reduce capital expenses and operating costs and increase plant profitability by enabling use of less expensive alloys, decreasing maintenance, and increasing plant longevity. Key focus areas for follow-on work have been identified as 1) evaluation of the efficacy of the electrochemical method under turbulent conditions in a larger flow system, such as the FASTR loop, 2) development of methods for removal of purification byproducts, 3) modeling pilot and industrial scale performance of electrochemical salt purification during plant operation and 4) further assessment of the effect of impurities on salt vapor phase.

14 SOLAR ENERGY↗