Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reaction rates”

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 775 records · Page 43

Formation of phenanthrenyl radicals via the reaction of acenaphthyl with acetylene

Ab initio G3(MP2,CC)/B3LYP/6–31G** calculations of the potential energy surface for the reaction of 3-acenaphthyl radical with acetylene were combined with Rice-Ramsperger-Kassel-Marcus Master Equation calculations of the temperature- and pressure-dependent reaction rate constants and product branching ratios with the goal to examine the transformation of an edge five-membered ring into a six-membered one. The reaction mechanism has been shown to contain two stages: (1) the formation of 3-ethynylacenaphthylene + H in the 3-acenaphthyl + C 2 H 2 reaction via acetylene addition to the radical site and an H loss, (2) the formation of the 1-phenanthrenyl radical in the 3-ethynylacenaphthylene + H reaction via H addition to the five-membered ring followed by the ring flip. The results of kinetic modeling show this pathway being faster than the pathway proposed earlier featuring the CCH chemisorption on the six-membered ring without involvement of the five-membered ring. A kinetic test demonstrated the dominance of the former pathway for the formation of pyrene from naphthalene. Here, it is expected that the new mechanism plays a major role in the transformation of five-membered rings on edges of PAH structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Antarctic stratospheric chemistry of chlorine nitrate, hydrogen chloride, and ice - Release of active chlorine

The reaction rate between atmospheric hydrogen chloride (HCl) and chlorine nitrate (ClONO2) is greatly enhanced in the presence of ice particles; HCl dissolves readily into ice, and the collisional reaction probability for ClONO2 on the surface of ice with HCl in the mole fraction range from about 0.003 to 0.010 is in the range from about 0.05 to 0.1 for temperatures near 200 K. Chlorine is released into the gas phase on a time scale of at most a few milliseconds, whereas nitric acid (HNO3), the other product, remains in the condensed phase. This reaction could play an important role in explaining the observed depletion of ozone over Antarctica; it releases photolytically active chlorine from its most abundant reservoir species, and it promotes the formation of HNO3 and thus removes nitrogen dioxide from the gas phase. Hence it establishes the necessary conditions for the efficient catalytic destruction of ozone by halogenated free radicals.

Molina, Mario J.↗

Elucidating the Role of Water on Limonene Oxidation with H 2 O 2 over γ -Al 2 O 3

Limonene oxide, which is produced from limonene epoxidation, is a valuable molecule that can be applied in flavor, fragrance, and renewable polymer applications. A catalytic reaction system using H 2 O 2 with γ-Al 2 O 3 and ethyl acetate (EtOAc) as the solvent has been explored as an effective system for this reaction. In these previous studies, a number of postulates have been proposed as to how water affects the reaction; therefore, the focus of this work is to elucidate the role of water in limonene epoxidation. While not impacting the selectivity to limonene oxide, the amount of water in the reaction system is shown to significantly impact the limonene reactivity. Furthermore, through both addition of excess water and removal of water with a Dean–Stark apparatus, the control of the H 2 O 2 /H 2 O ratio is demonstrated to be the primary factor controlling reactivity. In contrast, changes in limonene concentrations for a specific H 2 O 2 /H 2 O ratio are shown to have little impact on the reaction rate. This study shows that the competitive adsorption of H 2 O 2 and water on the catalyst surface is key in explaining the water impact on the reaction performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transverse liquid fuel jet breakup, burning, and ignition

An analytical study of the breakup, burning, and ignition of liquid fuels injected transversely into a hot air stream is conducted. The non-reacting liquid jet breakup location is determined by the local sonic point criterion. Two models, one employing analysis of an elliptical jet cross-section and the other employing a two-dimensional blunt body to represent the transverse jet, were used for sonic point calculations. An auxiliary criterion based on surface tension stability is used as a separate means of determining the breakup location. For the reacting liquid jet problem, a diffusion flame supported by a one-step chemical reaction within the gaseous boundary layer is solved along the ellipse surface in subsonic cross flow. Typical flame structures and concentration profiles were calculated for various locations along the jet cross-section as a function of upstream Mach numbers. The integration reaction rate along the jet cross-section is used to predict ignition position, which is found to be situated near the stagnation point. While a multi-step reaction is needed to represent the ignition process more accurately, the present calculation does yield reasonable predictions concerning ignition along a curved surface.

Li, Hsi-Shang↗

A Green's Function Approach to Simulate DNA Damage by the Indirect Effect

The DNA damage is of fundamental importance in the understanding of the effects of ionizing radiation. DNA is damaged by the direct effect of radiation (e.g. direct ionization) and by indirect effect (e.g. damage by.OH radicals created by the radiolysis of water). Despite years of research, many questions on the DNA damage by ionizing radiation remains. In the recent years, the Green's functions of the diffusion equation (GFDE) have been used extensively in biochemistry [1], notably to simulate biochemical networks in time and space [2]. In our future work on DNA damage, we wish to use an approach based on the GFDE to refine existing models on the indirect effect of ionizing radiation on DNA. To do so, we will use the code RITRACKS [3] developed at the NASA Johnson Space Center to simulate the radiation track structure and calculate the position of radiolytic species after irradiation. We have also recently developed an efficient Monte‐Carlo sampling algorithm for the GFDE of reversible reactions with an intermediate state [4], which can be modified and adapted to simulate DNA damage by free radicals. To do so, we will use the known reaction rate constants between radicals (OH, eaq, H,...) and the DNA bases, sugars and phosphates and use the sampling algorithms to simulate the diffusion of free radicals and chemical reactions with DNA. These techniques should help the understanding of the contribution of the indirect effect in the formation of DNA damage and double‐strand breaks.

Plante, Ianik↗

Inorganic Sulfur Species Formed upon Heterogeneous OH Oxidation of Organosulfates: A Case Study of Methyl Sulfate

In this report recent studies reveal that organosulfates at the particle surface can be oxidized by gas-phase OH radicals with significant rates. Inorganic sulfur species, such as the bisulfate ion (HSO 4 - ) and sulfate ion (SO 4 2- ), can be formed upon these heterogeneous oxidation processes through the formation and subsequent reactions of sulfate radical anions (SO 4 •- ) in the particle phase. However, the amount of inorganic sulfur species produced in these heterogeneous oxidation reactions is not known. We investigate the heterogeneous OH oxidation of sodium methyl sulfate (CH 3 SO 4 Na), the smallest organosulfate detected in atmospheric particles, using an oxidation flow reactor at a relative humidity of 75%. We quantify the kinetics by measuring the decay of CH 3 SO 4 Na and the amount of HSO 4 - and SO 4 2- formed upon oxidation using ion chromatography. Kinetic measurements determine the heterogeneous OH reaction rate to be (5.72 ± 0.14) x 10 -13 cm 3 molecule -1 s -1 , with an effective OH uptake coefficient, γ eff , of 0.31 ± 0.06. The molar yield of inorganic sulfur species, defined as the total number of moles of HSO 4 - and SO 4 2- formed per mole of CH 3 SO 4 Na consumed upon oxidation, is found to be significant and has an average value of 0.62 ± 0.18 upon oxidation. A kinetic model is developed to describe the kinetics and inorganic sulfur species formation upon oxidation. Model simulations suggest that CH 3 SO 4 Na tends to decompose rapidly into formaldehyde and SO 4 •- , and the reaction of SO 4 •- with CH 3 SO 4 Na plays a significant role in both governing the kinetics and the formation of inorganic sulfur species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting Thermal Response in a Li-ion Cell on a UAV Fight Profile

As the energy storage devices continue to "pack" more energy in a small space, any damage, battery component failure, manufacturing defect, or electrically abusing the battery can lead to catastrophic thermal runaway events. A catastrophic thermal event in a cell leads to high temperature, in some instances spewing of battery materials due to gas development from side reactions initiated due to high internal temperatures. Also, a thermal runaway event can propagate from a single "failed" cell to the pack in a battery pack, leading to a more significant event. Mitigating a thermal runaway event is important in the commercial and automotive sectors. However, preventing such events in an electric aircraft (or air taxis) is paramount due to the lack of alternatives in the event of a failure. Battery prognostics algorithms allow the prediction of state-of-charge (SOC) and end-of-life (EOL) of a Li-ion battery in a UAV (unmanned air vehicle) [1]. For this presentation, we will extend this two-level battery predictive algorithm to predict SOC, EOL, and estimated maximum temperature during a simulated flight. The model is extended by integrating a lumped physics-driven thermal model for high current densities [2]. The parameters used to control SOC and EOL are maximum storable charge, time constant for Li-ion diffusivity in the carbon particles, and internal cell resistance. Cycling leads to an increase in the heat generated by an aged Li-ion cell with a LiyCoO2 (LCO) cathode and a LixC6 (MCMB) anode. The aging of a cell leads to increase in SEI layer thickness, the diffusion time for the lithium ions inside the electrodes, and the local reaction rates, in addition to the thermodynamic abuse caused by fixed cycling voltages controlled by a Battery Management System. As the battery ages, the cell resistance increases, while the onset temperature of the thermal runaway decreases (depends on the cell chemistry and cell abuse history). Any large deviation of the cell temperature from the estimated (expected) value can identify a faulty cell. Since SEI decomposition has the lowest onset temperature in the series of reactions leading to thermal runaway, the model considers the self-heating rate of the SEI decomposition as onset temperature (similar to Ref. [3]). The parameters in the Arrhenius equation for the SEI heating rate depend on the number of cycles, the cell's operating temperature, and the cell's abuse history [4,5]. Coupling the electrochemical, thermal, and aging model allow the prognostic algorithm to estimate a typical cell voltage and temperature as a function of age (cycling and calendar), whose departure from measured values from the BMS is used to identify a safety event. In addition, we will present the results from two simulated flight scenarios for a UAV: typical and extreme, since the power requirements vary significantly during take-off, landing, and changing altitudes, while the power requirements remain low during the cruise. For this presentation, the power requirement for a battery pack in a UAV is scaled to a single cell. This cell is cycled through a simulated profile, and the data is collected and used to predict a safety event.

Li-ion↗

Solvent effects on the heterogeneous growth of TiO 2 nanostructure arrays by solvothermal synthesis

One-dimensional titanium dioxide (TiO 2 ) nanostructure arrays (nanoarrays) are important metal oxide nanomaterials that can be synthesized via facile solvothermal methods. The properties of the organic solvents can impose significant effects on the microstructures and properties of the final products. However, the discussions were limited to the homogeneously nucleated TiO 2 nanomaterials in free-standing powder form, while the solvent effects are less understood during the heterogeneous growth of TiO 2 nanoarrays on a substrate surface. In this work, six organic compounds, namely 2-butanone, n-decane, n-hexane, toluene, ethylene glycol and ethanol, were selected as the solvents for the solvothermal synthesis of TiO 2 nanoarrays on the cordierite monolithic substrates. Special attentions are paid to the morphology, crystallinity, specific surface area, and porosity of the samples. The heterogeneous growth of TiO 2 nanoarrays on substrate surfaces is found to favor the solvents with moderate dielectric constants, which can be partially dissolved in the aqueous solution and modulate the reaction rate during the solvothermal synthesis. Organic solvents with low dielectric constants may result in a complete separation between the precursors and aqueous solution, and therefore slow down the overall reaction, causing the insufficient growth of the nanoarrays. The TiO 2 nanoarrays are obtained with optimum morphology from the combination of 2-butanone and titanium (IV) butoxide as solvent and precusor, respectively, with a high specific surface area up to 56 m 2 /g including cordierite substrate given a micron thickness. When loaded with Pt catalyst, the TiO 2 nanoarray-based monolithic catalysts show excellent low-temperature catalytic activity and hydrothermal stability for the CO and hydrocarbon oxidation under the simulated exhausted conditions. This article shall shed light on a better understanding of the growth mechanism and rational design of TiO 2 nanoarrays for high-performance catalytic converters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Oxidation kinetics of biochar from woody and herbaceous biomass

The goal of this study was to determine oxidation kinetics for biochar produced from fast pyrolysis of various biomass feedstocks. In particular, the role of inherent ash content on the oxidation rate was evaluated. Thermogravimetric analysis (TGA) and fluidized bed combustion experiments were used to explore oxidation kinetics of six fast pyrolysis produced biochars with diverse ash content. Here, the reaction rates varied by a factor of three under chemical kinetic-limited conditions, demonstrating inorganic content impacts oxidation rate. Chemical kinetic rate coefficients were proposed as a function of compositional parameters to determine overall fit and impact. Potassium content was found to have a positive correlation, best describing the differences in the oxidation kinetic rate coefficients. Additionally, feedstocks were subjected to a 1 M HCl acid wash mitigating the catalytic activity of the metals. Acid washed biochars had lower oxidation kinetic rates compared to their unwashed counterparts, indicating the removal of catalytically active metals reduced oxidation rate. Gas composition (CO/CO 2 ) was measured during fluidized bed experiments for both acid-washed and unwashed biochars, which varied between the six biochars. Formation of CO 2 was greatly affected by catalytic metals, finding potassium content to correlate well with a higher percentage of CO 2 formation as compared to CO. Comparison of oxidation rates were made between the two experimental apparatuses to measure the effect of attrition on biochar oxidation.

09 BIOMASS FUELS↗

Initial Assessment of TeF 6 Adsorption

Research and development that supports the management of off-gases from nuclear fuel reprocessing has historically been focused on the off-gas streams that arise from aqueous reprocessing technology. However, as Gen-IV reactor development pathways move toward deployment, alternative spent nuclear fuel (SNF) processing and disposition pathways have been considered more actively. This work is focused upon aspects of fluoride volatility (FV) processing. The versatility of this method for deployment against multiple types of spent fuel encourages the continued advancement of both the primary separations processes and the secondary processes, including waste treatment, material control and accountability, and engineering designs. Recent work noted that TeF 6 , the most highly volatile fluorinated compound produced during FV processing, did not have a clear abatement technology recommended in the literature. Thus, this work performed scoping tests on activated alumina and copper shot to assess whether they could be used to remove TeF 6 from gas streams that bear F 2 . Previous work on this topic was not well described in the literature and was not performed with excess F 2 in the stream as would be typical of spent fuel processing via FV. To support an understanding of the concentration of TeF 6 contacting the adsorbent beds, a series of preliminary testing identified the TeF 6 production rate and the equilibrium concentration of TeF 6 in the gas stream contacting the adsorbent. Excess F 2 was determined to not affect the ability of activated alumina to remove TeF 6 from the gas stream quickly and completely. A determination of whether excess F 2 affected the distribution depth of Te in the sorbent bed is still pending analysis of the used sorbent. Literature suggests that when activated alumina is near saturation, TeF 6 could migrate from the sorbent bed. Future testing should investigate this possibility. Copper metal did not adsorb TeF 6 in the presence of F 2 across the sorbent temperature range of 50 to 335°C. F 2 was fully removed by the copper bed. Although preliminary thermodynamics would indicate adsorption to be energetically favorable, other factors that impact adsorption (e.g., slow kinetics, excess fluorine on the copper surface, an unfavorable transition state) are likely preventing the adsorption of TeF 6 at an easily measurable rate. The work also investigated TeF 6 production from multiple forms of Te in the temperature range of 100 to 250°C. No previous study had assessed the initial reaction rates for this process. A carefully designed study allowed determination of the activation energy for the production of TeF 6 from Te metal. The substantial amount of data collected during this study merits analysis beyond what is described here. A more in-depth kinetic analysis will be pursued. Additional analytical results will provide the ability to benchmark adsorption coefficients for TeF 6 , understand the distribution of TeF 6 within the alumina bed, and better understand the effect of F 2 partial pressure on Te fluorination. The data from this report, as supplemented by these additional analyses, will be submitted to a peer-reviewed journal.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Photochemical Phosphorus-Hydrogen-Oxygen Network for Hydrogen-dominated Exoplanet Atmospheres

Due to the detection of phosphine (PH 3 ) in the solar system gas giants Jupiter and Saturn, PH 3 has long been suggested to be detectable in exosolar substellar atmospheres too. However, to date, direct detection of phosphine has proven to be elusive in exoplanet atmosphere surveys. We construct an updated phosphorus-hydrogen-oxygen (PHO) photochemical network suitable for the simulation of gas giant hydrogen-dominated atmospheres. Using this network, we examine PHO photochemistry in hot Jupiter and warm Neptune exoplanet atmospheres at solar and enriched metallicities. Our results show for HD 189733b-like hot Jupiters that HOPO, PO, and P 2 are typically the dominant P carriers at pressures important for transit and emission spectra, rather than PH 3 . For GJ1214b-like warm Neptune atmospheres our results suggest that at solar metallicity PH 3 is dominant in the absence of photochemistry, but is generally not in high abundance for all other chemical environments. At 10 and 100 times solar, small oxygenated phosphorus molecules such as HOPO and PO dominate for both thermochemical and photochemical simulations. The network is able to reproduce well the observed PH 3 abundances on Jupiter and Saturn. Despite progress in improving the accuracy of the PHO network, large portions of the reaction rate data remain with approximate, uncertain, or missing values, which could change the conclusions of the current study significantly. Improving understanding of the kinetics of phosphorus-bearing chemical reactions will be a key undertaking for astronomers aiming to detect phosphine and other phosphorus species in both rocky and gaseous exoplanetary atmospheres in the near future.

atmospheric composition↗

Simplified Modeling of Oxidation of Hydrocarbons

A method of simplified computational modeling of oxidation of hydrocarbons is undergoing development. This is one of several developments needed to enable accurate computational simulation of turbulent, chemically reacting flows. At present, accurate computational simulation of such flows is difficult or impossible in most cases because (1) the numbers of grid points needed for adequate spatial resolution of turbulent flows in realistically complex geometries are beyond the capabilities of typical supercomputers now in use and (2) the combustion of typical hydrocarbons proceeds through decomposition into hundreds of molecular species interacting through thousands of reactions. Hence, the combination of detailed reaction- rate models with the fundamental flow equations yields flow models that are computationally prohibitive. Hence, further, a reduction of at least an order of magnitude in the dimension of reaction kinetics is one of the prerequisites for feasibility of computational simulation of turbulent, chemically reacting flows. In the present method of simplified modeling, all molecular species involved in the oxidation of hydrocarbons are classified as either light or heavy; heavy molecules are those having 3 or more carbon atoms. The light molecules are not subject to meaningful decomposition, and the heavy molecules are considered to decompose into only 13 specified constituent radicals, a few of which are listed in the table. One constructs a reduced-order model, suitable for use in estimating the release of heat and the evolution of temperature in combustion, from a base comprising the 13 constituent radicals plus a total of 26 other species that include the light molecules and related light free radicals. Then rather than following all possible species through their reaction coordinates, one follows only the reduced set of reaction coordinates of the base. The behavior of the base was examined in test computational simulations of the combustion of heptane in a stirred reactor at various initial pressures ranging from 0.1 to 6 MPa. Most of the simulations were performed for stoichiometric mixtures; some were performed for fuel/oxygen mole ratios of 1/2 and 2.

Bellan, Josette↗

A-Priori Tuning of Modified Magnussen Combustion Model

In the application of CFD to turbulent reacting flows, one of the main limitations to predictive accuracy is the chemistry model. Using a full or skeletal kinetics model may provide good predictive ability, however, at considerable computational cost. Adding the ability to account for the interaction between turbulence and chemistry improves the overall fidelity of a simulation but adds to this cost. An alternative is the use of simple models, such as the Magnussen model, which has negligible computational overhead, but lacks general predictive ability except for cases that can be tuned to the flow being solved. In this paper, a technique will be described that allows the tuning of the Magnussen model for an arbitrary fuel and flow geometry without the need to have experimental data for that particular case. The tuning is based on comparing the results of the Magnussen model and full finite-rate chemistry when applied to perfectly and partially stirred reactor simulations. In addition, a modification to the Magnussen model is proposed that allows the upper kinetic limit for the reaction rate to be set, giving better physical agreement with full kinetic mechanisms. This procedure allows a simple reacting model to be used in a predictive manner, and affords significant savings in computational costs for simulations.

Norris, A. T.↗

REACTER 2.0: Quantum-Informed Reaction Constraints and Automated Interaction Typing

REACTER is a heuristic method for modeling chemical reactions in classical molecular dynamics simulations, implemented in LAMMPS as fix bond/react. The authors recently extended LAMMPS to support alphanumeric labels for atom types, bond types etc., which enables the pre- and post-reaction templates required by the REACTER protocol to be portable between different simulations and greatly simplifies the task of creating simulation-ready reaction templates. To further increase the generality of reaction templates, support for wildcard characters within atom types has been added, along with the automatic assignment of interaction types for new bonds, angles, etc. based on the involved atom types. In some cases, this feature can express a class of reactions with one pair of reaction templates, where previously dozens may have been required. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, as well as a custom constraint that leverages LAMMPS’ powerful built-in variable framework. Other new features include variable support for various inputs (e.g., to allow reaction rates or cutoffs to be dependent on overall conversion), on-the-fly update of molecule IDs, and the ability to create new atoms positioned with respect to the reaction site. The new features are applied to modeling polymeric, thermosetting and composite materials, and advanced applications of the new reaction constraints are demonstrated. For example, REACTER is shown to accurately reproduce mechanically-induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a constraint on the total potential energy of the involved atoms.

polymer simulations↗

AGC-4 As-Run Uncertainty Analysis

The concept of measuring neutron flux (or fluence rate) with activation foils has existed since the late 1940s and is a well-used methodology. Activation foils are selected for their unique isotopic cross sections that create reaction products with desirable properties: relatively simple radioactive decay scheme, measurable quantities of gamma-rays or betas, and half-lives appropriate for the irradiation and decay durations of the experiment. The activation materials, deemed flux wires, fluence wires, dosimeters, or flux foils, are fabricated to have a mass appropriate for the irradiation as well so the measurement using a radiation detector following the irradiation is sufficient but not overwhelming in terms of total activity. Activation measurement results are expressed in terms of reaction rate, which can then be used to “adjust” the neutron spectrum. The Advanced Graphite Creep (AGC)-4 experiment contained 26 flux wire packages, each containing an iron and a niobium wire inside a vanadium capsule. For all stackup positions, the value of fast to thermal neutrons is higher than expected, but the total neutron fluence is near the predicted value of ~2 x 1022 n/cm2. As compared to the predicted neutron spectrum, there are four takeaways in the results as compared to the predicted ATR irradiation: 1) Additional fast (fission) neutrons present, 2) Spectrum hardening due to temperature, 3) Spectrum hardening due to container attenuation, 4) Total neutron fluence is close to predicted values. Since the neutron energy spectrum is a key parameter when attributing neutron damage in materials, the adjusted energy-dependent spectrum and the total neutron fluence results from the dosimeters can be incorporated for future calculations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Absolute rate of the reaction of bromine atoms with ozone from 200-360 K

The rate constant for the reaction Br + O3 yields BrO + O2 was measured from 200 to 360 K by the technique of flash photolysis coupled to time resolved detection of bromine atoms by resonance fluorescence (FP-RF). Br atoms were produced by the flash photolysis of CH3Br at lambda 165nm.O3 was monitored continuously under reaction conditions by absorption at 253.7 nm. At each of five temperatures the results were independent of substantial variations in O3, total pressure and limited variations in flash intensity. The measured rate constants obeyed the Arrhenius expression, where the error quoted is two standard deviations. Results are compared with previous determinations which employed the discharge flow-mass spectrometric technique.

Michael, J. V.↗

Analysis of thermochemical energy storage in metal carbonates: characterizing cycling-induced degradation

A solid–gas reaction rate modeling framework is applied to characterize degradation in thermochemical materials induced by thermal cycling. Time constants and a kinetic conversion ratio are established for a representative carbonation-calcination reaction, which provides insight into degradation mechanisms and reveals a mitigation strategy by tuning reaction duration.

25 ENERGY STORAGE↗