Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactivity Coefficients”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Dramatic Slowdown of Bimolecular Reactions of Criegee Intermediates in Organic Aerosols

Although Criegee intermediates (CIs) play a central role in driving the transformation of organic molecules, their condensed-phase reaction kinetics remains poorly understood. While CIs undergo unimolecular decomposition to produce OH radicals, stabilized CIs participate in a host of other bimolecular reactions forming larger molecular weight products, including secondary ozonides (kSOZ), alkoxyalkyl hydroperoxides (kAlAHP), and acyloxyalkyl hydroperoxides (kAcAHP). Here, we demonstrate, using prior literature and new measurements of heterogeneous organic aerosol reactions using O3 or OH radicals, that these condensed-phase bimolecular reaction rate coefficients kSOZ (10-20 cm3·molecule-1·s-1), kAlAHP, and kAcAHP (10-19 cm3·molecule-1·s-1) are 7-10 orders of magnitude slower than analogous gas-phase reactions (∼ 10-10-10-13 cm3·molecule-1·s-1). This contrasts with current evidence that the rate coefficients for unimolecular steps in the condensed phase are similar in magnitude to those in the gas phase. Kinetic models using gas-phase-like bimolecular coefficients consistently overestimate the importance of bimolecular relative to unimolecular pathways while underpredicting the aerosol reactivity by up to 30-fold. The slowing of bimolecular reactivity in aerosols is likely due to solvation effects and the increased steric hindrance of larger molecular weight CIs and their scavengers in the condensed phase. These results challenge the conventional view that the high reactant density of an organic aerosol should favor bimolecular over unimolecular reactions, underscoring the need for more accurate CI condensed-phase rate constants for reliable modeling of oxidative processes across a diverse set of atmospheric and environmental systems.

Zeng, Meirong↗
DOE OSTI2020

A Multidimensional Operando Study Showing the Importance of the Electrode Macrostructure in Lithium Sulfur Batteries

Lithium sulfur batteries are one of the most promising next-generation energy storage technologies, because of their impressive theoretical energy density, low materials cost, and relative safety. However, incomplete understanding of their underlying operation mechanisms has hindered their further development and commercialization. In this work, to gain a better understanding of the operation mechanisms in the lithium sulfur battery, three macroscopically different (woven and nonwoven) and microscopically similar (made from the same carbon fibers) free-standing carbon felts are used as the conductive matrices to show how the macrostructure influences the redox reactions and the electrochemical performance. Additionally, an operando radiography and simultaneous in situ electrochemical impedance spectroscopy study are performed to highlight the differences of the carbon hosts and how they compare electrochemically in a lithium sulfur battery. The electrochemical results show that the carbon host with a more open structure results in increased capacity as well as a higher diffusion coefficient. The operando radiography shows that the open structures are more conducive to efficient redox reactivity in forming solid sulfur species, mostly in the form of elongated needlelike structures identified as β-S8, near the end of charge, which is corroborated by the impedance data. The data also reveal a “breathing” mechanism, where the polysulfides are pushed to the edges during discharge and pulled back during charge. Furthermore, it is also revealed that the “breathing” process is a limiting factor in the high rate performance. Overall, this study shows that there is a close relationship between the macrostructure of the carbon matrix and the electrochemical performance of the sulfur cathode.

25 ENERGY STORAGE↗
DOE OSTI2021

Predicting CaO-(MgO)-Al2O3-SiO2 glass reactivity in alkaline environments from force field molecular dynamics simulations

In this investigation, force field-based molecular dynamics (MD) simulations have been employed to generate detailed structural representations for a range of amorphous quaternary CaO-MgO-Al{sub 2}O{sub 3}-SiO{sub 2} (CMAS) and ternary CaO-Al{sub 2}O{sub 3}-SiO{sub 2} (CAS) glasses. Comparison of the simulation results with select experimental X-ray and neutron total scattering and literature data reveals that the MD-generated structures have captured the key structural features of these CMAS and CAS glasses. Based on the MD-generated structural representations, we have developed two structural descriptors, specifically (i) average metal oxide dissociation energy (AMODE) and (ii) average self-diffusion coefficient (ASDC) of all the atoms at melting. Both structural descriptors are seen to more accurately predict the relative glass reactivity than the commonly used degree of depolymerization parameter, especially for the eight synthetic CAS glasses that span a wide compositional range. Hence these descriptors hold great promise for predicting CMAS and CAS glass reactivity in alkaline environments from compositional information.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Energy, angle, and temperature dependencies of the sticking of D atoms on Li surfaces

Detailed experimental and computational information on the response of lithium surfaces to irradiation by slow hydrogenic particles (ions, atoms, molecules) is sparse and mainly speculative. In this work, we present a computational study of the reflection and retention of deuterium (D) atoms at crystalline and amorphous lithium surfaces at 300 and 500 K, where the D atoms have an impact energy in the range of 0.025–5 eV and incident angles of 0° (perpendicular incidence) or 85° (near-grazing incidence). Classical molecular dynamics simulations are performed with the reactive bond-order force field (ReaxFF) potentials. This study provides quantitative information on the deuterium sticking probability and recycling coefficient for lithium surfaces. Further, our results support the ongoing work at the Lithium Tokamak eXperiment-β fusion experiment as well as relevant experiments in the laboratory setting.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗
DOE OSTI2021

Inelastic, exchange, and reactive processes in rovibrationally excited collisions of HD with H

ABSTRACT The HD molecule is an important coolant in early universe chemistry models and a tracer of H2 in star-forming regions. Rate coefficients for collisional excitation and de-excitation of HD rotational and vibrational levels form important ingredients in astrophysical models. While collisions with He, H2, and H are the most important, available data for H + HD collisions are largely limited to temperatures less than 1000 K for the vibrational ground state, low-lying rotational levels of the v = 1 HD vibrational level, or computed without reactive contributions. Here, through explicit quantum scattering calculations, we report extensive data for rovibrational transitions in HD induced by H atoms for a range of rotational levels in v = 1 and some v = 0 levels for temperatures up to 1000 K. The significance of the computed results for astrophysical modeling is discussed.

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

Diffusion of Atoms in Glassy Mixtures of Deuterium and Lithium

The diffusion coefficients of D and Li were calculated in amorphous, glassy Li+D mixtures for various concentrations of D in lithium over the temperature range 100–1000 K. The densities of the mixtures as a function of temperature were also determined. The diffusion coefficients were obtained by the analysis of the mean-squared displacement using molecular dynamics with ReaxFF, a reactive force field. In conclusion, the diffusion and density data obtained for Li+D mixtures were compared with the available experimental and calculated data reported in the literature.

Diffusion↗
DOE OSTI2023

High-fidelity multiphysics load following and accidental transient modeling of microreactors using NEAMS tools: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The feasibility of modeling microreactors using high-fidelity models with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tools is investigated in this report. Three overarching questions guided this research: can NEAMS tools readily be applied for high-fidelity multiphysics modeling of different types of transients in microreactor designs; how accurate are the results obtained; and are improvements needed in accuracy or user experience of NEAMS tools, especially considering newly developed capabilities? This work builds upon FY-2022 work, and two microreactor concepts considering heat pipe (HP-MR) and gas-cooled (GC-MR) technologies were further analyzed using high-fidelity multiphysics simulations. The NEAMS tools considered and coupled within the MultiApp environment are Griffin for neutronics, BISON for thermo-mechanics, Sockeye for heat pipe modeling (in HP-MR), SAM for 1D Fluid – 3D solid modeling of coolant channels and system modeling of balance of plant components (in GC-MR), and the SWIFT code for hydrogen redistribution in hydride moderator. The Heat Pipe MicroReactor (HP-MR) concept was further analyzed in FY-2023 to demonstrate the stochastic TRISO failure modeling capability in BISON to check operational limits of the TRISO fuel. A new full-core Gas-Cooled MicroReactor (GC-MR) model was developed based on the initial assembly-model used in Y-2022 and used for steady-state and accidental depressurization transient simulations. Accuracy of the simulations performed was assessed through 1) verification analyses completed on the different physics with code-to-code comparison, and 2) validation of the multiphysics simulations based on modeling of the Kilopower Reactor Using Stirling Technology (KRUSTY) experiment. In FY-2023, the mesh and model of KRUSTY was updated to closely match publicly available data, and the neutronic model was verified and validated against experimental control rod worth measurements. The multiphysics model of KRUSTY was developed and used for steady-state analysis and for modeling reactivity insertion transient. The calculated power increase and stabilization agrees well with experimental data following adjustment in fuel thermal expansion coefficient. As an important component of this project, the ANL team gathered experience with a wide range of NEAMS tools: the MOOSE Mesh System, Griffin, BISON, SWIFT, Sockeye, SAM, Workbench, and the MOOSE MultiApp System, and provided assessment of new capabilities. Noteworthy are the user assessment of the “vapor-only” flow model in Sockeye and development of a multiphysics startup transient in HP-MR unit cell for use as tutorial in Sockeye. The full-core GC-MR model was used for assessment of SAM for balance of plant modeling and for demonstrating the SWIFT code capability for hydrogen redistribution modeling in multiphysics transient analyses. In this process, several bugs/issues were identified and reported to developers. Finally, the assembly GC-MR model developed in FY-2022 coupling Griffin, BISON and SAM through flow blockage and rod ejection transients was published to the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB). The Heat Pipe MicroReactor (HP-MR) concept high-fidelity multiphysics coupling of Griffin/BISON/Sockeye in load-following and heat pipe failure transients was also published on the VTB. Those submissions are enabling thorough review of these models as well as wide distribution to industry, regulator, and university users. In this analysis, several new research questions were uncovered, and follow-up analyses are recommended to further improve some models, consider additional transients, and continue development of VTB models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Xenon compatibility in magmatic processes: Hadean to current contexts

Xenon (Xe) behaviour in petrological processes, albeit essential to constrain mantle ingassing and degassing models, is elusive due to its volatile nature, and lack of direct investigation at the pressures (P) and temperatures (T) relevant to magma formation and crystallisation at depth. Xenon stands out amongst noble gases due to its unique reactivity with silicates of the lower crust and upper mantle, which could at least partially explain that published mineral/melt partitioning coefficients span up to six orders of magnitude. We report partition coefficients of Xe using in situ X-ray fluorescence at high P and T, and mass spectrometry analyses. Xenon is found to be moderately incompatible in anorthite-clinopyroxene mix in equilibrium with basalt (partition coefficient value of 0.16 ± 0.06), and compatible in olivine in equilibrium with basalt (partition coefficients in the range 88 ± 22 to 302 ± 46). While Xe is, thus, concentrated in basaltic melts coexisting with crystallising pyroxenes and feldspars, it is strongly retained in olivine at depth. Consequently, Xe originally contained in solid Earth has been preferentially retained at depth throughout Earth’s history, from the magma ocean stages to present day partial mantle melting processes.

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

Quantum reactive scattering calculations for the cold and ultracold Li + LiNa → Li 2 + Na reaction

A first-principles based quantum dynamics study of the Li + LiNa(ν = 0, j = 0) → Li 2 (ν', j') + Na reaction is reported for collision energies spanning the ultracold (1 nK) to cold (1 K) regimes. A full-dimensional ab initio potential energy surface for the ground electronic state of Li 2 Na is utilized that includes an accurate treatment of the long-range interactions. The Li + LiNa reaction is barrierless and exoergic and exhibits a deep attractive potential well that supports complex formation. Thus, significant reactivity occurs even for collision temperatures approaching absolute zero. The reactive scattering calculations are based on a numerically exact time-independent quantum dynamics methodology in hyperspherical coordinates. Total and rotationally resolved rate coefficients are reported at 56 collision energies and include all contributing partial waves. Several shape resonances are observed in many of the rotationally resolved rate coefficients and a small resonance feature is also reported in the total rate coefficient near 50 mK. Furthermore, of particular interest, the angular distributions or differential cross sections are reported as a function of both the collision energy and scattering angle. Unique quantum fingerprints (bumps, channels, and ripples) are observed in the angular distributions for each product rotational state due to quantum interference and shape resonance contributions. The Li + LiNa reaction is under active experimental investigation so that these intriguing features could be verified experimentally when sufficient product state resolution becomes feasible for collision energies below 1 K.

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

Bimolecular Reaction of Methyl-Ethyl-Substituted Criegee Intermediate with SO 2

Methyl-ethyl-substituted Criegee intermediate (MECI) is a four-carbon carbonyl oxide that is formed in the ozonolysis of some asymmetric alkenes. MECI is structurally similar to the isoprene-derived methyl vinyl ketone oxide (MVK-oxide) but lacks resonance stabilization, making it a promising candidate to help us unravel the effects of size, structure, and resonance stabilization that influence the reactivity of atmospherically important, highly functionalized Criegee intermediates. We present experimental and theoretical results from the first bimolecular study of MECI in its reaction with SO 2 , a reaction that shows significant sensitivity to the Criegee intermediate structure. Using multiplexed photoionization mass spectrometry, we obtain a rate coefficient of (1.3 ± 0.3) × 10 –10 cm 3 s –1 (95% confidence limits, 298 K, 10 Torr) and demonstrate the formation of SO 3 under our experimental conditions. Through high-level theory, we explore the effect of Criegee intermediate structure on the minimum energy pathways for their reactions with SO 2 and obtain modified Arrhenius fits to our predictions for the reaction of both syn and anti conformers of MECI with SO 2 (k syn = 4.42 × 10 11 T –7.80 exp(-1401/T) cm 3 s –1 and k anti = 1.26 × 10 11 T –7.55 exp(-1397/T) cm 3 s –1 ). In conclusion, our experimental and theoretical rate coefficients (which are in reasonable agreement at 298 K) show that the reaction of MECI with SO 2 is significantly faster than MVK-oxide + SO 2 , demonstrating the substantial effect of resonance stabilization on Criegee intermediate reactivity.

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

Combining direct ink writing with reactive melt infiltration to create architectured thermoelectric legs

Here, we present a new additive-reactive synthesis method where inks – cast into molds or 3D-additively extruded into architectured shapes – are reacted into intermetallic thermoelectric compounds. The new method, as demonstrated for equiatomic TiNiSn, combines: (i) extrusion printing (or casting) of inks containing Ni and Ti powders, (ii) debinding and reactive sintering to form a porous NiTi network, (iii) network infiltration with liquid Sn and subsequent reaction to synthesize the TiNiSn phase. Thin plates, created through this method, show high phase purity and low residual porosity. A thermoelectric figure of merit zT = 0.47 ± 0.05 is achieved at 800 K, within the broad range of values (zT= 0.3–1.0) reported for TiNiSn created via traditional powder hot-pressing methods. Layered TiNiSn microlattices – after 3D ink-extrusion printing and infiltration – exhibit high relative densities and minimal undesirable secondary phase content, with a Seebeck coefficient on par with the ink-cast TiNiSn plates, demonstrating that thermoelectric legs with far-ranging architectural freedom can be created additively with this novel reactive manufacturing method. Numerous other thermoelectric (and other) compounds, currently limited to basic geometries due to brittleness, are amenable to this new method based on a sequence of reactive sintering/infiltration of cast or ink-extruded precursors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na + , Ca 2+ , and Mg 2+ cations and the CO 3 2– anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me–O–C (Me = metal) three-body valence angle parameters and Me–C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na + predominantly forms carbonates due to weaker angular strain, while Ca 2+ and Mg 2+ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. Here, the nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca 2+ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na + and Mg 2+ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

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

An age- and sex-specific biokinetic model for radon *

Publication 137 of the International Commission on Radiological Protection (ICRP) describes a biokinetic model for radon used to derive dose coefficients for occupational intake of radon isotopes. The model depicts transfer of inhaled or ingested radon to blood, exchange of radon between blood and tissues, and gradual loss of radon from the body based on physical laws governing transfer of a non-reactive and soluble gas between materials. Here, this paper describes an age- and sex-specific variation of that model developed for use in an upcoming ICRP series of reports on environmental intake of radionuclides by members of the public titled ‘Dose Coefficients for Intakes of Radionuclides by Members of the Public’. The proposed model modifies the model structure and transfer coefficients presented in Publication 137 to allow more realistic dosimetric treatment of bone marrow and breast and expands the model to address pre-adult ages, based on the physical principles used in the development of the model of Publication 137 together with anatomical and physiological changes occurring during human development.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

The neutron number probability distribution in coupled lumped assemblies

Here, the validity of the gamma distribution in describing the neutron number probability distribution function for both isolated and coupled multiplying assemblies when constrained to reproduce the true mean and variance is investigated in lumped geometry by numerical comparison with kinetic Monte Carlo simulations. The mean and variance are obtained from numerical solution of moment equations constructed from the relevant forward Master equation with assembly coupling coefficients obtained from a view factor method. Numerical results for a two-group, two coupled assemblies model, with static and dynamic reactivity insertion, show that except for subcritical assemblies, the gamma distribution well-approximates the number distribution. Differences in the fast and thermal neutron population shapes are explained in terms of effective source strengths due to downscatter and coupling.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

New capabilities of the MORET 6 Monte Carlo neutron transport code

The MORET code is a simulation tool that solves the transport equation for neutrons using the Monte Carlo method. It allows users to model complex three-dimensional geometrical configurations in a user-friendly way. New features have been introduced to extend the application field of MORET beyond the usual criticality calculations for which it has been initially designed. The most important change is the addition of an analog fixed source mode which allows studies of systems of any reactivity combined with very flexible outputs. Other useful improvements have been added concerning the geometric part, the fission matrix, the multigroup sensitivity coefficients and the outputs. This paper presents an overview of these new features. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Experimental investigation of power transient flow boiling

To advance the mechanistic understanding of power transients in flow boiling, experimental investigations are performed to capture the thermal response of four different kinds of cladding materials, including FeCrAls, to Fuchs reactivity initiated accident transients and linear ramp power transients. The transition boiling regime is phenomenologically reflected in power-transient flow boiling, which is different from steady-state boiling. In light of this, the critical heat flux that is usually featured with temperature overshooting does not indicate thermal safety margin of cladding materials as conservatively as the maximum heat flux, which is greater than critical heat flux in the power transient boiling curve. This is physically attributed to the thermal energy deposition on the cladding wall. In addition, the thermal energy released from the cladding wall results in different boiling heat transfer coefficient for the decreasing and increasing power stages respectively during the Fuchs power transient. The maximum heat flux difference gap, which appears due to changes in cladding materials and transient time scales, is appreciable under an intermediate heat convection regime. However, this difference gap is gradually weakened by the progressive increasing of mass flux and/or inlet subcooling due to the enhanced dominance of heat convection over heat conduction. Finally, the small difference gap of maximum heat flux is found to be insignificant under the weak heat convection because of power transient induced annular flow instability.

42 ENGINEERING↗
DOE OSTI2022

Sensing performance of sub-100-nm vanadium oxide films for room temperature thermal detection applications

Vanadium oxide films are widely employed as thermal detectors in uncooled infrared detection systems due to their high temperature coefficient of resistance near room temperature. One strategy toward maximizing detectivity and reducing the thermal time constant in these systems is to minimize the system platform dimensions. This approach necessitates thinner film thicknesses ([Formula: see text]100 nm), for which there is little information regarding thermal sensing performance. Herein, we report on the sensitivity of reactively sputtered vanadium oxide thin film resistive thermometers nominally ranging from 100 to 25 nm and assess the influence of thermal annealing. Here we demonstrate that films in this minimum limit of thickness maintain a high temperature coefficient while additionally providing an enhancement in characteristics of the noise equivalent power.

36 MATERIALS SCIENCE↗