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 361 records · Page 20

Thermal diffusivity recovery and defect annealing kinetics of self-ion implanted tungsten probed by insitu transient grating spectroscopy

Tungsten is a promising candidate material for plasma-facing armour components in future fusion reactors. A key concern is irradiation-induced degradation of its normally excellent thermal transport properties. In this comprehensive study, thermal diffusivity degradation in ion-implanted tungsten and its evolution from room temperature (RT) to 1073 K is considered. Five samples were exposed to 20 MeV self-ions at RT to achieve damage levels ranging from 3.2 x 10 -4 to 3.2 displacements per atom (dpa). Transient grating spectroscopy with in-situ heating was then used to study thermal diffusivity evolution as a function of temperature. Using a kinetic theory model, an equivalent point defect density is estimated from the measured thermal diffusivity. The results showed a prominent recovery of thermal diffusivity between 450 K and 650 K, which coincides with the onset of mono-vacancy mobility. After 1073 K annealing samples with initial damage of 3.2 x 10 -3 dpa or less recover close to the pristine value of thermal diffusivity. For doses of 3.2 x 10 -2 dpa or higher, on the other hand, a residual reduction in thermal diffusivity remains even after 1073 K annealing. Transmission electron microscopy reveals that this is associated with extended, irradiation-induced dislocation structures that are retained after annealing. A sensitivity analysis shows that thermal diffusivity provides an efficient tool for the assessing total defect content in tungsten up to 1000 K.

36 MATERIALS SCIENCE↗

Low-Temperature Plasma Assisted Kinetics Study of Ethanol

Non-equilibrium plasma assisted kinetics of ethanol at 0.5 atm over a temperature range of 523-1203 K were studied at stoichiometric conditions by using a plasma flow reactor (PFR). The effects of plasma chemistry on fuel and its influence on neutral chemistry is decoupled with both experimental and modelling efforts. Enhanced reactivity was observed experimentally in both plasma-assisted oxidation and pyrolysis cases with fuel consumption beginning at temperatures as low as 523 K. Absolute fuel consumption was observed at temperatures nearly 200 K earlier than pure thermal reactions highlighting the effectiveness of non-equilibrium plasma to promote reactivity. Furthermore, modelling results illustrated the effect of electron impact reactions towards accelerating fuel consumption and promoting low temperature chemistry by generating excited state N2 and reactive O and OH radicals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Griffin Software Development Plan

Griffin is a MOOSE-based reactor physics application for advanced reactor multiphysics modeling and simulation. The application is developed in a consistent multiphysics environment with strong software quality assurance. Griffin inherited most of the capabilities of MAMMOTH/Rattlesnake and is adopting the capabilities from PROTEUS that are needed in the code. The toolset includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, ad-joint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. The code contains the cross-section preparation capabilities applicable to fast and thermal reactors, including TRISO-fueled reactors. Core management capabilities include core performance, fuel depletion and shuffling, equilibrium core calculation, pebble-bed reactor run-in and equilibrium core, molten-salt reactor delayed neutron precursor drift, and control rod and drum movement with cusping correction. This software development plan presents the current and future capabilities and features in Griffin for the design and analysis of non-light-water reactor systems in steady-state and transient conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.

97 MATHEMATICS AND COMPUTING↗

Mobility assessment of the BCC and carbide phases in the C-Nb, C-U and Nb-U systems

Uranium carbides with refractory metal additions are considered for Gen IV nuclear reactors and nuclear thermal propulsion as fuels for their high-temperature and corrosion resistant properties. Understanding kinetic effects that dictate microstructural evolution during fabrication and operating conditions is essential to advance technological development of these fuels. This work presents the development of an atomic mobility database for C-Nb-U systems based off available experimental data supported with ab-initio methods. The mobility assessments and uncertainty quantification (using Markov chain Monte Carlo) were conducted in the Kawin software. Carbon diffusion is considered dominant, as metal diffusion is much slower, with niobium diffusion being even slower and rate limiting than uranium metal. We provide a comprehensive and self-consistent thermo-kinetic database that is validated by diffusion couple simulations through Kawin. In conclusion, this enables prediction of microstructural and phase evolution critical for the development and lifetime assessment of next generation nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Selective Methane Oxidation to Methanol on ZnO/Cu 2 O/Cu(111) Catalysts: Multiple Site-Dependent Behaviors

Due to the abundance of natural gas in our planet, a major goal is to achieve a direct methane to methanol conversion at medium to low temperatures using mixtures of methane and oxygen. Here, we report an efficient catalyst, ZnO/Cu 2 O/Cu(111), for this process investigated using a combination of reactor testing, scanning tunneling microscopy, ambient-pressure x-ray photoemission spectroscopy, density functional calculations, and kinetic Monte Carlo simulations. The catalyst is capable of methane activation at room temperature and transforms mixtures of methane and oxygen to methanol at 450 K with a selectivity of ~30%. This performance is not seen for other heterogenous catalysts which usually require the addition of water to enable a significant conversion of methane to methanol. The unique coarse structure of the ZnO islands supported on a Cu 2 O/Cu(111) substrate provides a collection of multiple centers that display different catalytic activity during the reaction. ZnO-Cu 2 O step sites are active centers for methanol synthesis when exposed to CH 4 and O 2 due to an effective O-O bond dissociation, which enables a methane-to-methanol conversion with a reasonable selectivity. Upon addition of water, the defected O-rich ZnO sites, introduced by Zn vacancies, show superior behavior towards methane conversion and enhance the overall methanol selectivity to over 80%. Thus, in this case, the surface sites involved in a direct CH 4 → to CH 3 OH conversion are different from those engaged in methanol formation without water. Furthermore, the identification of the site-dependent behavior of ZnO/Cu 2 O/Cu(111) opens a design strategy for guiding efficient methane reforming with high methanol selectivity.

36 MATERIALS SCIENCE↗

A kinetics study of the homogeneous and heterogeneous components of the HCl + ClONO2 reaction

The kinetics of the reaction HCl + ClONO2 to Cl2 + HNO3 were investigated at 298 K using a flow reactor with FTIR analysis to assess the importance of this reaction for stratospheric chemistry. The observed reaction was characteristic of a heterogeneous process; an upper limit of 5 x 10 to the -18th cu cm/molecule per s was obtained for the homogeneous gas phase rate constant. From calculations of the first order wall rate constant, estimates were made of the reaction rate on stratospheric aerosols. Because both HCl and ClONO2 need to be adsorbed on the particle surface, the reaction will be of negligible importance under most stratospheric conditions.

Friedl, Randall R.↗

Rare-isotope and kinetic studies of Pt/SnO2 catalysts

Closed-cycle pulsed CO2 laser operation requires the use of an efficient CO-O2 recombination catalyst for these dissociation products which otherwise would degrade the laser operation. The catalyst must not only operate at low temperatures but also must operate efficiently for long periods. In the case of the Laser Atmospheric Wind Sounder (LAWS) laser, an operational lifetime of 3 years is required. Additionally, in order to minimize atmospheric absorption and enhance aerosol scatter of laser radiation, the LAWS system will operate at 9.1 micrometers with an oxygen-18 isotope CO2 lasing medium. Consequently, the catalyst must not only operate at low temperatures but must also preserve the isotopic integrity of the rare-isotope composition in the recombination mode. Several years ago an investigation of commercially available and newly synthesized recombination catalysts for use in closed-cycle pulsed common and rare-isotope CO2 lasers was implemented at the NASA Langley Research Center. Since that time, mechanistic efforts utilizing both common and rare oxygen isotopes have been implemented and continue. Rare-isotope studies utilizing commercially available platinum-tin oxide catalyst have demonstrated that the catalyst contributes oxygen-16 to the product carbon dioxide thus rendering it unusable for rare-isotope applications. A technique has been developed for modification of the surface of the common-isotope catalyst to render it usable. Results of kinetic and isotope label studies using plug flow, recycle plug flow, and closed internal recycle plug flow reactor configuration modes are discussed.

Upchurch, Billy T.↗

Oxidation kinetics of model compounds of metabolic waste in supercritical water

In this NASA-funded study, the oxidation kinetics of methanol and ammonia in supercritical water have been experimentally determined in an isothermal plug flow reactor. Theoretical studies have also been carried out to characterize key reaction pathways. Methanol oxidation rates were found to be proportional to the first power of methanol concentration and independent of oxygen concentration and were highly activated with an activation energy of approximately 98 kcal/mole over the temperature range 480 to 540 C at 246 bar. The oxidation of ammonia was found to be catalytic with an activation energy of 38 kcal/mole over temperatures ranging from 640 to 700 C. An elementary reaction model for methanol oxidation was applied after correction for the effect of high pressure on the rate constants. The conversion of methanol predicted by the model was in good agreement with experimental data.

Webley, Paul A.↗

Chemical Kinetics Underlying the Sooting Tendency and Auto-Ignition Characteristics of Linear, Branched, and Cyclic Ether Compounds

Biofuels present opportunities for improving the performance and reducing emissions from internal combustion engines by incorporating oxygenated functional groups to the fuels. Among various oxygenates, ethers have been recognized as promising candidates for an alternative to conventional diesel fuel owing to their higher reactivity and lower sooting tendency. The detailed guidelines for designing ethers, however, have not been fully discussed, even though their combustion characteristics are sensitive to the molecular structure. This study was devoted to exploring the structure-property relationships, particularly focusing on the cetane number and yield sooting index, using five linear, branched, and cyclic ethers: di-amyl-ether, 4-butoxy-heptane, 3,3-dimethyl-oxetane, 2-ethyl-4-methyl-1,3-dioxolane, and 2-isopropyl-4-methyl-1,3-dioxolane. First, we examined the chemical kinetics underlying the sooting tendency of the test fuels. The combustion product distribution was measured from flow reactor experiments at 750-1100 K, F=3, at atmospheric pressure. As a result, it was revealed that the sooting tendency is closely related to the size of hydrocarbon intermediates in the high-temperature regime (>1000 K); that is, larger hydrocarbons lead to more soot precursor formation. The underlying chemistry determining the size of the hydrocarbons from the tested fuels was analyzed using reaction pathway analysis and quantum mechanics calculations, which showed that the branched and cyclic ether structures form abundant C3-C4 compounds. Moreover, the auto-ignition characteristics of the test fuels were studied using the flow reactor at low-temperature (400-700 K) and F=1. We found a clear difference in the combustion-product distribution from high and low reactivity fuels, which was then correlated to the systematic analysis of the key reaction energy barriers with the varying molecular structure.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗

Preservation of kinetics parameters generated by Monte Carlo calculations in two-step deterministic calculations

The generation of accurate kinetic parameters such as mean generation time Λ and effective delayed neutron fraction β eff via Monte Carlo codes is established. Employing these in downstream deterministic codes warrants another step to ensure no additional error is introduced by the low-order transport operator when computing forward and adjoint fluxes for bilinear weighting of these parameters. Another complexity stems from applying superhomogenization (SPH) equivalence in non-fundamental mode approximations, where reference and low-order calculations rely on a 3D full core model. In these cases, SPH factors can optionally be computed for only part of the geometry while preserving reaction rates and K-effective, but the impact of such approximations on kinetics parameters has not been thoroughly studied. This paper aims at studying the preservation of bilinearly-weighted quantities in the Serpent–Griffin calculation procedure. Diffusion and transport evaluations of IPEN/MB-01, Godiva, and Flattop were carried out with the Griffin reactor physics code, testing available modeling options using Serpent-generated multigroup cross sections and equivalence data. Verifying Griffin against Serpent indicates sensitivities to multigroup energy grid selection and regional application of SPH equivalence, introducing significant errors; these were demonstrated to be reduced through the use of a transport method together with a finer energy grid.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A wide range experimental study and further development of a kinetic model describing propane oxidation

This work presents an experimental and kinetic modeling study of propane oxidation. Ignition delay times of propane were measured in a high-pressure shock tube and in rapid compression machines in the temperature range 689 - 1700 K at equivalence ratios of 0.5, 1.0 and 2.0 in 'air', for a wide range of pressures from 20 to 90 bar. CO and H 2 O mole fraction profiles for propane oxidation were measured in a shock tube behind reflected shock waves in the temperature range 1370-1840 K at equivalence ratios of 0.5, 1.0 and 2.0 and at a pressure of approximately 1.3 atm. Moreover, propane oxidation was studied using a jet-stirred reactor coupled to a synchrotron vacuum ultraviolet photoionization mass spectrometer at low temperatures in the range 565 - 690 K and at a pressure of 1 atm. Further, this wide range of experimental datasets for propane oxidation was used to reoptimize and update our previous kinetic mechanisms, AramcoMech3.0 and NUIGMech1.1. In the current mechanism, NUIGMech1.3, the thermo-chemical parameters of all species relevant to low-temperature propane oxidation chemistry, including propyl-peroxyl, hydroperoxyl-propyl, hydroperoxyl-propyl-peroxyl, and carbonyl-hydroperoxide radicals, are updated based on newly calculated values at the CCSD(T)-F12/TZ-F12//B2PLYPD3/TZ///B2PLYP-D3/TZ level of theory. The improvements made in the thermochemical values and in the kinetic parameters for the low-temperature propane oxidation reactions in NUIGMech1.3 result in better model agreement with the new IDTs and speciation data, including carbon monoxide, formaldehyde, propene, acetaldehyde and various minor products such as ethylene, acetic acid, acrolein as well as various hydroperoxide and cyclic ether species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of Planar-Cavity Receiver Reactor for Solar Thermochemical Dry-Reforming

Concentrating solar thermal (CST) systems can be leveraged to provide not only heat for power generation, but also for energy storage and thermochemical fuel production. Such solar thermochemical processes have been studied conceptually, from solar thermochemical hydrogen production (STCH) and thermochemical energy storage (TCES), to gasification, reforming, and fuel upgrading by various means. The solar receiver and reactor are critical components in the conversion of solar energy into chemical energy in the form of 'solar fuels'. For effective conversion of solar energy within a coupled solar receiver-reactor, extremely high temperatures are required, thereby demanding a high solar concentration ratio (CR) for efficient operation. This creates a design challenge for the receiver-reactor, as many thermochemical processes involve gas or gas-solid systems that are limited by low heat transfer coefficients. A unique receiver design is proposed that has the potential to incorporate various high-temperature thermochemical processes such as TCES-assisted power generation, methane reforming, or STCH processes. Modeling this receiver and its potential applications requires a full three-dimensional model that accurately captures the interconnected effects of receiver geometry, spatial solar irradiance, complex radiation, reaction kinetics, fluid dynamics, and heat transfer. In this paper we analyze a CST system integrated with this unique planar-cavity receiver-reactor design using the developed model. The model presented in this paper showed where improved thermal management was needed to achieve suitable receiver performance when a dry-methane reforming process is implemented.

14 SOLAR ENERGY↗

Isomeric effects on the reactivity of branched alkenes: An experimental and kinetic modeling study of methylbutenes

Here, a detailed experimental study of the low-to-intermediate temperature combustion of methylbutene isomers, i.e., branched C 5 alkenes, has been undertaken with multiple experimental facilities. Ignition delay times were measured at equivalence ratios 0.5–2.0, 685–1020 K and up to 45 bar condition from two rapid compression machines and showed slight deviation from an Arrhenius behavior for all three isomers, while their reactivity order differs as temperature changes. Sampled intermediates formed during the oxidation process of mixtures at 900–1150 K and 0.82 bar from a flow reactor and at 730 K and 20 bar from a rapid compression machine were analyzed using gas chromatography techniques. Trends in the formation and consumption of sampled intermediates were modeled using a kinetic model developed in this work for all three isomers. Rate of production and sensitivity analyses emphasize the role of double bond-specific reactions governing the global reactivity of these fuels. Additional studies of the addition reactions of HO 2 radicals to the double bond and to allylic radicals may improve the model performance.

2-Methyl-1-butene↗

Ion irradiation effects on Cr-coated zircaloy-4 surface wettability and pool boiling critical heat flux.

The concept of coating the currently used nuclear fuel cladding (zirconium-based alloy, typically Zircaloy-4 or Zirc-4) with an oxidation preventive layer is a progressing Accident tolerant Fuel (ATF) candidate alloys. The coated Zirc-4-based alloys could be a solution to suppress undesirable fast reaction kinetics with high-temperature steam. Zirc-4 has been the most preferred cladding material in pressurized water reactors (PWRs). Chromium (Cr) based alloys as a coating material provides excellent corrosion protection and good strength and wear resistance. In this paper, we present the surface wettability measurements and pool boiling Critical Heat Flux (CHF) for Cr-coated Zirc-4 claddings pre- and post-exposure to an ion irradiation environment. The wettability measurements, including static contact angle (contact angle, θ) and average surface roughness (surface roughness, Ra), are introduced for samples of different coating thicknesses (5–30 μm thick). The coatings fabricated by the cold spray of Cr-Al particles to 10 mm × 10 mm × 1.95 mm Zirc-4 substrates. Post fabrication, a Pilgering (cold rolling) process, was applied to finalize the coating thickness and resulted in a significant reduction in surface roughness of initially fabricated rough surfaces. The process produced three distinguished samples 5-μm unpolished (as machined), 5-μm, and 30-μm polished (cold rolled). The measurements are presented for the three surfaces and bare Zirc-4 as a baseline surface. The contact angle analyses were implemented in theoretical models from the literature to predict pool boiling CHF. Pool boiling experiments were conducted to measure the pool boiling CHF values and compare them to the predicted values. Scanning Electron Microscope (SEM) images and Energy Dispersive X-ray Spectroscopy (EDS) analysis was performed to characterize the surfaces for better understanding and interpreting the results. The SEM images showed localized surface damage due to ion irradiation. No recognized change in the measured surface roughness due to ion irradiation. The contact angles of irradiated Cr-coated surfaces are consistently higher (10°) than pre-irradiated surfaces. Decreasing the Cr-coating layer thickness resulted in lower contact angle pre- and post- ion irradiation. The predicted pool boiling CHF using the Kandlikar model is in good agreement with the experimentally measured CHF values within ±12% for all samples.

36 MATERIALS SCIENCE↗

Irradiation-induced damage in concrete-forming aggregates: revisiting literature data through micromechanics

The radiation-induced volumetric expansion (RIVE) of aggregate-forming minerals causes damage in concrete exposed to high levels of fast neutrons fluence ( \(>\,\sim \,10^{19}\, \hbox {n}\,\hbox {cm}^{-2}\) at kinetics energy above 0.1 MeV). Historical post-irradiation RIVE and Young’s modulus data obtained in test reactors were revisited using a polycrystalline homogenization model (self-consistent scheme) accounting for the aggregates’ minerals content and the formation of voids/cracks during irradiation. It was found that the formation of extra voids/cracks can contribute to the aggregate expansion more substantially than the cumulated expansions of aggregate-forming minerals. Furthermore, the rate of void creation appears to be inversely correlated to the silica content, although aggregates of higher silicate contents exhibit higher RIVEs. The loss of Young’s modulus decreases exponentially with the aggregates’ expansion. At a comparable RIVE level, the relative loss of modulus is more significant in aggregates of higher silica content.

36 MATERIALS SCIENCE↗

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗

Improved Breathing Cartridge (IBC)

The Fire Cartridge (FC) is designed to work with the Emergency Mask (EM) to protect the crew by cleaning the inhaled air in a post fire environment. The Improved Breathing Cartridge (IBC) is designed as a next generation FC with the lessons learned from years of building and testing the FC design, as well as novel catalyst materials to improve the performance and robustness of the cartridge. The FC was designed with minimal margin to save on mass, while the IBC was sized to have over 100% excess catalyst to achieve the required conversion, and still reduce the production cost. Kinetic data taken from another project was used to size the catalyst bed from the standard plug flow reactor conversion equation. The cartridge has been 3D printed and will soon be filled with the required levels of media. Testing has not yet been performed to determine final performance.

Fire Response↗