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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2021

A boundary value “reservoir problem” and boundary conditions for multi-moment multifluid simulations of sheaths

Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma–wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remain an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2004

Magnetic-Nozzle Studies for Fusion Propulsion Applications: Gigawatt Plasma Source Operation and Magnetic Nozzle Analysis

This project has been a multiyear effort to assess the feasibility of a key process inherent to virtually all fusion propulsion concepts: the expansion of a fusion-grade plasma through a diverging magnetic field. Current fusion energy research touches on this process only indirectly through studies of plasma divertors designed to remove the fusion products from a reactor. This project was aimed at directly addressing propulsion system issues, without the expense of constructing a fusion reactor. Instead, the program designed, constructed, and operated a facility suitable for simulating fusion reactor grade edge plasmas, and to examine their expansion in an expanding magnetic nozzle. The approach was to create and accelerate a dense (up to l0(exp 20)/m) plasma, stagnate it in a converging magnetic field to convert kinetic energy to thermal energy, and examine the subsequent expansion of the hot (100's eV) plasma in a subsequent magnetic nozzle. Throughout the project, there has been a parallel effort between theoretical and numerical design and modelling of the experiment and the experiment itself. In particular, the MACH2 code was used to design and predict the performance of the magnetoplasmadynamic (MPD) plasma accelerator, and to design and predict the design and expected behavior for the magnetic field coils that could be added later. Progress to date includes the theoretical accelerator design and construction, development of the power and vacuum systems to accommodate the powers and mass flow rates of interest to out research, operation of the accelerator and comparison to theoretical predictions, and computational analysis of future magnetic field coils and the expected performance of an integrated source-nozzle experiment.

Gilland, James H.↗

DOE OSTI2020

Pulse Radiolysis Study of Radiation Effect on Molten Salt

The molten salt reactor concept, where nuclear fuel is dissolved in a molten salt that also serves as the heat transfer fluid, is a leading candidate for next generation nuclear reactors. This necessitates a thorough understanding of radiation effects on molten salt media to support the design, development, and deployment of such reactors. Early, pioneering pulse radiolysis molten salt experiments observed key primary radiolysis products such as the solvated electron (esolv?) and dichloride radical anion (Cl2??), but left many fundamental mechanistic and reactivity questions unanswered. Here we report on the reaction kinetics of esolv? and Cl2?? in molten LiCl-KCl eutectic salt doped with Zn2+ ion, using electron pulse radiolysis to observe the transient behavior from nanosecond to microsecond time scales. Experiments were performed at the BNL Laser-Electron Accelerator Facility using a recently-developed high-temperature sample holder. Prompt formation of esolv? and Cl2?? are observed; esolv? decays within hundreds of nanoseconds while Cl2?? decays more slowly by second-order kinetics, more likely by cross-recombination rather than disproportionation. This work was supported as part of the Molten Salts in Extreme Environments Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science.

38 - RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCL↗

DOE OSTI2021

Effects of non-thermal termolecular chemistry on detonation development in hydrogen (H2) / methane (CH4) - air mixtures

In combustion simulations, it is usually assumed that highly-energized collision complexes (either radical or stable species) formed in exothermic reactions are always under thermal equilibrium as they evolve through competing networks of reactions. However, in practical flames, non-trivial amounts of reactive radicals such as H, O and OH are present apart from O2. As a result, collisions of these reactive species with the energized collision complexes have been shown in recent studies to induce non-thermal reactivity. These studies have also demonstrated that such non-thermal reactions can be suitably represented in macroscopic kinetics models as chemically termolecular reactions. The present work was focused on identifying and quantifying the effects of including such reactions on the evolution of an initial deflagration front to a developing detonation in H2/CH4-air mixtures under boosted internal combustion (IC) engine conditions. Specifically, fully resolved simulations, with and without non-thermal reactivity, were performed for a constant volume reactor containing stoichiometric H2/CH4-air mixture. It was found that inclusion of termolecular chemistry resulted in a delayed onset of end-gas auto-ignition. Concurrently, the developing detonation intensity was observed to be significantly higher. Chemical explosive mode analysis (CEMA) was performed to identify the dominant species/reactions responsible for the observed variation in the results.

Detonation↗

DOE OSTI2021

Non-Boltzmann Effects in Chain Branching and Pathway Branching for Diethyl Ether Oxidation

Low-temperature (LT) engine applications have several potential benefits, including reduced emissions and increased efficiency. Attaining these benefits requires accurate kinetic modeling of LT chain branching, which depends heavily on ketohydroperoxide (KHP) decomposition. For diethyl ether (DEE), a promising biofuel, current estimates of the KHP decomposition rate constant are largely based on empirical fits to data. In this study, we investigate the most important KHP isomer in DEE LT oxidation by applying variable reaction coordinate transition state theory to the main pathway for KHP decomposition: OO bond fission to produce •OH and a keto-alkoxy radical, •OQ'O. We also use ab initio kinetics methods to investigate the decomposition of •OQ'O, where we find dominant branching to acetic acid, with the remaining flux going to CH 3 C(O)OCHO. Additionally, new time-resolved measurements of DEE and acetic acid concentrations during LT (450–600 K) DEE oxidation are obtained in a laser photolysis flow reactor coupled with multiplexed photoionization mass spectrometry. These new experimental data, along with jet-stirred reactor data in the literature, are compared with the predictions of a recent DEE mechanism (Tran et al. Proc. Comb. Inst. 2019, 37, 511-519) that was modified with the newly calculated ab initio rate constants for KHP and •OQ'O decomposition. The predictions of the modified mechanism are quite poor when compared to the experimental data; this is primarily due to the new KHP ⇌ •OQ'O + •OH rate constant, which is 1–2 orders of magnitude slower than empirical values employed in recent mechanisms. To reconcile the new KHP rate constant and the experimental data, we explore and quantify the possible role of non-Boltzmann (nB) reaction sequences. The nB reactions have a substantial effect on both the overall mechanism reactivity and the •OQ'O branching to acetic acid. We also provide guidance on the proper implementation of nB reactions in kinetic mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Microwave-Assisted Plastic Upcycling: Dynamic Data Reconciliation, Parameter Estimation, and Kinetic Modeling

Microwave (MW)-assisted catalytic pyrolysis offers a promising pathway for efficient plastic upcycling. This work develops an integrated modeling framework combining dynamic data reconciliation, a temperature-dependent rate model, and a yield model to represent the time-varying production rate of components in MW-assisted LDPE pyrolysis conducted in a batch reactor. An Arrhenius-type rate model with a temperature-dependent reaction order is developed. A biexponential correlation is proposed for the yield of gaseous products that enables to capture the evolving product formation behavior during conversion. In the yield correlation, one term is used to represent the initial increase in yield, reflecting the rapid formation of intermediate or primary products at the early stages of the reaction when a larger fraction of the reactant remains available. As conversion progresses, the influence of this term gradually diminishes. The other term accounts for the subsequent decrease in the predicted yield, representing secondary reactions such as further cracking or coke formation that reduce the concentration of certain products at higher conversion. The model is found to accurately represent reconciled experimental flow rate profiles from an in-house MW-assisted catalytic batch reactor for major products, including ethylene, ethane, 1-butene, and benzene, across 250−350 °C. Ethylene remains the dominant product but decreases from about 41.95% at 250 °C to 30.14% at 350 °C, while heavier products increase significantly, with 1-butene rising to nearly 8.37% and benzene reaching 2.17% at intermediate temperatures. The model shows that the ethylene production rate can be maximized at around 270 °C. The models developed in this work can be utilized for process optimization, reactor design and scale-up of microwave-assisted plastic conversion technologies, and economic analysis.

Damahe, Harish [West Virginia Univ., Morgantown, W↗

NASA NTRS1991

Mechanisms of ilmenite reduction and their impact on the design of effective reactor systems

One of the first activities at a lunar base could be oxygen recovery from ilmenite (FeTiO3). Oxygen produced from lunar soils could be used to fuel transportation vehicles operating in near-earth space. The first step in developing a suitable reactor system for lunar operation is to determine the mechanisms and rates of oxygen removal from ilmenite. In-situ gravimetric measurements and microscopic examinations were used to determine the hydrogen reduction mechanisms of synthetic ilmenite discs between 823 to 1353 K. A shrinking core reaction model, modified to account for the growth of an iron film on the surface of discs, was capable of predicting experimentally observed conversion-time relationships. The observed reduction mechanism, kinetic rates, and associated activation energy established a base line from which comparisons could be made to improve oxygen yield and removal rates. One proposed technique to improve the rate and extent of oxygen removal was to preoxidize ilmenite. Preoxidation is commonly used during the reduction of ilmenite ores in the iron industry and has been employed for many years to lower operating temperatures and increase reduction rates. This technology could prove beneficial for oxygen production facilities on the Moon as less massive reactor vessels and/or less energy could be associated with a process including preoxidation. X ray diffraction and energy dispersive spectroscopy were utilized to follow the progression of ilmenite oxidation at 1123 and 1140 K and the reduction of pseudobrookite (Fe2TiO5) at 873 and 973 K. Structures formed during the progress of oxidation were related to the system's phase diagrams. Results indicated that after initially producing ilmenite-hematite solutions and rutile (TiO2), pseudobrookite was the end product of oxidation at all temperatures examined (1049 to 1273 K). Initial results from the reduction of pseudobrookite indicate a series of phases are produced including ferropseudobrookite (FeTi2O5), ulvospinel (Fe2TiO4), and ilmenite. Rates of pseudobrookite reduction were typically 50 to 200 times that of ilmenite.

Briggs, R. A.↗

DOE OSTI2025

Validation of Numerical Tools for Calculating Reactivity Feedback in Sodium Fast Reactors Using SEFOR Experimental Data

The Southwest Experimental Fast Oxide Reactor (SEFOR) was an experimental sodium-cooled fast breeder reactor operated from 1969 to 1972 with experiments designed to measure Doppler reactivity feedback in a wide temperature range from around 350 °F to temperatures approaching the melting point of mixed oxide fuel of around 5000 °F, providing valuable data for code validations. Co-supported by the Department of Energy (DOE) Fast Reactor Program (FRP) and the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the SEFOR benchmark project focused on using the experimental data to validate numerical tools that are used in industry and academia to design and license sodium-cooled fast reactors (SFRs). By the end of FY-25, substantial progress was achieved in the SEFOR benchmark study. A variety of numerical tools commonly used for modeling SFRs were applied to develop models for SEFOR core configurations I-D, I-E, I-I, and I-J. These included Monte Carlo codes such as MCNP, Serpent, and Shift; deterministic codes such as the legacy Argonne Reactor Computation (ARC) suite and the high-fidelity NEAMS code Griffin; and the system analysis code SAS4A/SASSYS-1 (SAS). Using these models, both SEFOR zero-power experiments and power-ascending tests were successfully simulated. Comparisons were performed against experimental measurements of core criticalities, reflector worth, kinetics parameters (Λ/βeff), isothermal reactivity feedback (from 350 °F to 760 °F at zero power), and power-ascending reactivity feedback (as power increased from 0.4 MW to 17 MW). In general, these comparisons demonstrated very good agreement between numerical results and experimental data. In Fiscal Year 26 (FY-26), the SEFOR benchmark project will continue to address the modeling issues identified in FY-25. Effort will focus on the simulation of reactivity insertion transients in SEFOR core II using the ARC/SAS model. Future work will also focus on incorporating BISON into the SEFOR core modeling process to enable the first Multiphysics simulations of the isothermal tests based on the MOOSE framework.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Investigating biomass composition and size effects on fast pyrolysis using global sensitivity analysis and CFD simulations

It is notoriously difficult to build an accurate universal model for biomass pyrolysis due to its sensitivity to a wide number of critical material attributes such as chemical species and physical sizes. In this work, a biomass pyrolysis kinetics with 32 heterogeneous reactions and 59 species was implemented in an open-source multiphase computational fluid dynamics (CFD) software MFiX and validated against two different experimental pyrolysis data sets that provided detailed data describing chemical component yields. The reaction scheme was then used to build a surrogate model and assess the sensitivity of pyrolysis yields to feedstock compositions. The sensitivity analysis determined that the yield of bio-char showed a strong positive sensitivity to the carbon-rich lignin and tannin pseudo-species in the reaction scheme while the bio-oil and bio-gas were correlated to oxygen-rich lignin pseudo-species. The reaction scheme was then integrated into a coarse-grained discrete element model to simulate fast pyrolysis in a bubbling fluidized bed over a range of feedstock particle sizes. The reactor simulations showed further sensitivity to particle size and hydrodynamics. Notably, particles under 0.5 mm have small heat transfer limitations but left the reactor before completely converting and thus reduced the bio-oil yield. Results from this study can be used to guide future development of highly accurate models for fast pyrolysis reactors with a variety of feedstock properties and operating conditions.

09 BIOMASS FUELS↗

DOE OSTI2022

Formation of PAHs, phenol, benzofuran, and dibenzofuran in a flow reactor from the oxidation of ethylene, toluene, and n-decane

The necessity to reveal the formation mechanism of not only polycyclic aromatic hydrocarbons (PAHs) but also oxygenated PAHs (OPAHs) during combustion is increasing. Although many studies on PAHs have been conducted, fundamental studies investigating OPAH formation are still limited. Phenol, benzofuran, and dibenzofuran were selected as OPAHs in this study. We experimentally investigated fuel-rich oxidation in a flow reactor at atmospheric pressure, mean gas temperatures from 1050 to 1350 K, residence times from 0.2 to 1.5 s, and equivalence ratios from 3.0 to 12.0. Ethylene, toluene, and n-decane were used as fuels. Three kinds of OPAHs as well as 23 kinds of PAHs including monocyclic structures were quantitatively measured through a direct sampling using gas chromatography mass spectrometry. The results showed that concentrations of PAHs and OPAHs were strongly affected by the fuel type, with increasing concentrations in the following order: ethylene < n-decane < toluene. A chemical kinetic model for OPAH formation was developed based on a recent PAH growth model. The model showed that the predicted concentrations of OPAHs and PAHs were in reasonable agreement with the measured data in this study. Some modifications were made to the previous model based on recent literature studies and on the comparison of the simulated and measured results. The effect of the fuel type on the formation of PAHs and OPAHs was investigated through kinetic analysis using the model to discuss their reaction pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Development and preliminary validation of a mechanistic multiscale model for fuel-cladding chemical interaction in metallic nuclear fuels

Despite decades of fuel rod material and design improvements, fuel-cladding chemical interaction (FCCI) remains the single-most lifetime-limiting behavior for modern metallic fuel rods. Constraining fuel lifetime increases operating costs, limiting the economic viability of commercializing metallic nuclear fuel technology. A mechanistic multiscale model utilizing the finite element method-based MARMOT and BISON codes was developed to more confidently predict cladding-side FCCI and its impact on fuel performance. The new BISON model incorporates mesoscale models for the effects of fuel microstructure evolution on the transport of wastage-inducing lanthanides through the fuel and for the kinetics of cladding wastage layer growth. The mesoscale models, in turn, build on lanthanide transport property data obtained from the atomistic scale. Preliminary validation studies using wastage thickness and cladding profilometry data from four fuel rods irradiated in Experimental Breeder Reactor II experiment X447 and one fuel rod from Fast Flux Test Facility experiment IFR1 show that the new model predicts cladding wastage and its effects on cladding deformation as well as existing empirical FCCI correlations. The new model is expected to aid in the design of new metallic fuel concepts, including fuel additives, cladding liners, and sodium-free annular fuel geometries. In conclusion, future work will focus on broader validation and refinement of the model’s treatment of different fuel alloys and cladding materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Evaluation of lithium as a tritium storage medium for betavoltaics

Lithium foils were demonstrated to absorb surrogate protium for tritium-powered betavoltaics. 20μm thick lithium foils were hole-punched from a ribbon of electrodeposited lithium on copper foil. The lithium foils were loaded with hydrogen in a custom Sievert apparatus where the pressure drop showed full hydriding at a hydrogen pressure of 2 bar and at all loading temperatures above the lithium melting point at 190, 200, 225, 250, and 300. Lithium hydride formation was confirmed with Raman spectroscopy after hydrogen loading. The kinetics of experimental hydride formation was compared to the diffusion-limited Mintz–Bloch model. While the Mintz–Bloch model showed good fit with the experimental loadings, the model overpredicted the loading kinetics starting at 250 °C and at higher temperatures. In conclusion, The overprediction was either caused by lithium hydride outgassing due to some reduction with some residual lithium hydroxide created from brief air exposure when sealing the lithium in the reactor or a transition from diffusion-limited hydride growth to surface or metal–hydride interface-limited hydride growth.

25 ENERGY STORAGE↗

NASA NTRS1989

Production of oxygen from lunar ilmenite

The overall objective of this project was to develop a novel carbothermal reduction process for production of oxygen from lunar ilmenite. The specific objective was to use a reaction sequence in which a wide variety of carbonaceous compounds (including carbonaceous wastes) can be used as reducing agents. During the first phase, two reactor systems were designed, constructed, and operated to study the reaction fundamental important in this process. One system is a small fluidized bed, and the other is a thermo-gravimetric reactor system. Preliminary experiments on synthetic ilmenite are conducted to study the effect of carbon type, carbon loading, temperature, and gas flow rate. Results indicate that a reaction path based on carbon gasification can be used to promote the overall kinetics. A unique temperature and concentration-programmed reaction procedure was being developed for rapid parametric study of the process.

Shadman, F.↗

DOE OSTI2021

NRC Multiphysics Analysis Capability Deployment (FY2021 - Part 2)

This report details the progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” Task 4b was completed for this report. INL developed a sample problem showing how to use Serpent 2 to calculate macroscopic cross sections for use in Griffin for a typical sodium fast reactor Unprotected Loss of Flow (ULOF) transient. The complicating factor is the relatively large axial gradient of the coolant (and hence fuel) temperature and its evolution during the transient. A 3D Griffin model of the Advanced Burner Test Reactor (ABTR) is coupled to the System Analysis Module (SAM) to perform the ULOF transient. The model includes various explicit feedback mechanisms for fast reactors, including Doppler, radial expansion from the displacement of the support plate, and axial expansion from the displacement of the fuel pins. The results for the ULOF transient are consistent with published values. We observed small effects in the transient results that arise from the superhomogenization equivalence correction of uniform and nonuniform temperature data sets, but they diminish as more dominant reactivity mechanisms are added to the model. Potential improvements to the neutronics model include: adding the differential control rod and sodium feedback mechanisms, implementing cross section re-homogenization based on the partial volume of the various materials, and enhancing the kinetics parameters. Furthermore, the fluids model only includes four subchannels to represent the entire active core due to performance limitations in SAM when attempting to model one subchannel per assembly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Development of the TREAT M2 experiment as a transient benchmark

A transient benchmark based on the Transient Reactor Test Facility (TREAT) M2 Calibration experiment (M2-CAL) is under development. TREAT, at Idaho National Laboratory, is a graphite moderated air-cooled research reactor which has been used extensively for fuel material testing under extreme and accident conditions. Accurate benchmark models are a beneficial component in the operations, experimental planning, and development of TREAT. In this work, we present a transient benchmark model for the M2-CAL experiment core loading. The benchmark model incorporates the coupling between the Monte Carlo Code SERPENT and the computational fluid dynamics code OpenFOAM to capture the temperature feedback mechanism. The M2-CAL transient 2580 was simulated in this work. A pre-transient analysis was performed to determine the optimum core composition and conditions before the beginning of the transient. The analysis was validated against historic TREAT kinetic measurements and the worth of the transient rod T-2. The axial power distribution in the flux wire for the M2-CAL experiment was determined and contrasted with the experiment. The model was then used to simulate the M2-CAL transient 2580 experiment based on the reported pre-transient and transient conditions. Several transient observables were calculated and compared to the experiment. The model shows good agreement with the experimental power traces, period, and average increase of power at the power ramp (time > 7.8 s). Inverse point kinetics analysis was performed during the period of the power ramp for the model and the experiment based on that, the temperature feedback component was isolated and contrasted with the experimental feedback. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Experimental and computational study of particle formation kinetics in UF6 hydrolysis

The formation and growth of UO 2 F 2 particles by gas-phase UF 6 hydrolysis remains of interest to actinide chemistry researchers. The total number concentration of the UO 2 F 2 aerosol particles that can be produced in the reaction is regulated primarily by the availability of water molecules under our reactor conditions. An increase in water molecule concentration corresponds with a higher amount and larger size of UO 2 F 2 aerosol particles produced. The growth rates of aerosol particles appear to approach a single number in the range of [0.05 ± 0.03–0.08 ± 0.04] (nm s -1 ), as the molar ratio of water to UF 6 decreases below 1. The size of primary particles produced from the UF 6 hydrolysis under water-deprived conditions was estimated to be 3.6 ± 0.4 nm. As the molar ratio became greater than 1.7, the size of primary particles increased with increased availability of water molecules. The primary particle model developed in this work predicted a size range for the UO 2 F 2 primary particles similar to that estimated based on the data from gas-phase UF 6 hydrolysis experiments. This result suggests that the volume-driven coalescence process assumption used in the derivation of the primary particle model was reasonable. The ability to precisely control the availability of water molecules and reaction time could lead to the production of nearly monodispersed aerosol particles. This finding has significant implications in the engineering and manufacturing of fuel powder materials and possibly the future development and deployment of environmental sampling apparatus.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1972

Fusion power for space propulsion.

Principles of operation, interplanetary orbit-to-orbit mission capabilities, technical problems, and environmental safeguards are examined for thermonuclear fusion propulsion systems. Two systems examined include (1) a fusion-electric concept in which kinetic energy of charged particles from the plasma is converted into electric power (for accelerating the propellant in an electrostatic thrustor) by the van de Graaf generator principle and (2) the direct fusion rocket in which energetic plasma lost from the reactor has a suitable amount of added propellant to obtain the optimum exhaust velocity. The deuterium-tritium and the deuterium/helium-3 reactions are considered as suitable candidates, and attention is given to problems of cryogenic refrigeration systems, magnet shielding, and high-energy particle extraction and guidance.

Roth, R.↗

NASA NTRS1990

A unifying picture of gas-phase formation and growth of PAH (Polycyclic Aromatic Hydrocarbons), soot, diamond and graphite

A variety of seemingly different carbon formation processes -- polycyclic aromatic hydrocarbons and diamond in the interstellar medium, soot in hydrocarbon flames, graphite and diamond in plasma-assisted-chemical vapor deposition reactors -- may all have closely related underlying chemical reaction mechanisms. Two distinct mechanisms for gas-phase carbon growth are discussed. At high temperatures it proceeds via the formation of carbon clusters. At lower temperatures it follows a polymerization-type kinetic sequence of chemical reactions of acetylene addition to a radical, and reactivation of the resultant species through H-abstraction by a hydrogen atom.

Frenklach, Michael↗