Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Kinetic-Scale Magnetic Turbulence and Finite Larmor Radius Effects at Mercury

We use a nonstationary generalization of the higher-order structure function technique to investigate statistical properties of the magnetic field fluctuations recorded by MESSENGER spacecraft during its first flyby (01/14/2008) through the near-Mercury space environment, with the emphasis on key boundary regions participating in the solar wind - magnetosphere interaction. Our analysis shows, for the first time, that kinetic-scale fluctuations play a significant role in the Mercury's magnetosphere up to the largest resolvable timescale (approx.20 s) imposed by the signal nonstationariry, suggesting that turbulence at this plane I is largely controlled by finite Larmor radius effects. In particular, we report the presence of a highly turbulent and extended foreshock system filled with packets of ULF oscillations, broad-band intermittent fluctuations in the magnetosheath, ion-kinetic turbulence in the central plasma sheet of Mercury's magnetotail, and kinetic-scale fluctuations in the inner current sheet encountered at the outbound (dawn-side) magnetopause. Overall, our measurements indicate that the Hermean magnetosphere, as well as the surrounding region, are strongly affected by non-MHD effects introduced by finite sizes of cyclotron orbits of the constituting ion species. Physical mechanisms of these effects and their potentially critical impact on the structure and dynamics of Mercury's magnetic field remain to be understood.

Uritsky, V. M.↗
NASA NTRS2012

Chemical Kinetics of the TPS and Base Bleeding During Flight Test

The present research deals with thermal degradation of polyurethane foam (PUF) during flight test. Model of thermal decomposition was developed that accounts for polyurethane kinetics parameters extracted from thermogravimetric analyses and radial heat losses to the surrounding environment. The model predicts mass loss of foam, the temperature and kinetic of release of the exhaust gases and char as function of heat and radiation loads. When PUF is heated, urethane bond break into polyol and isocyanate. In the first stage, isocyanate pyrolyses and oxidizes. As a result, the thermo-char and oil droplets (yellow smoke) are released. In the second decomposition stage, pyrolysis and oxidization of liquid polyol occur. Next, the kinetics of chemical compound release and the information about the reactions occurring in the base area are coupled to the CFD simulations of the base flow in a single first stage motor vertically stacked vehicle configuration. The CFD simulations are performed to estimate the contribution of the hot out-gassing, chemical reactions, and char oxidation to the temperature rise of the base flow. The results of simulations are compared with the flight test data.

Osipov, Viatcheslav↗
NASA NTRS2013

Microwave Kinetic Inductance Detector with Selective Polarization Coupling

A conventional low-noise detector requires a technique to both absorb incident power and convert it to an electrical signal at cryogenic temperatures. This innovation combines low-noise detector and readout functionality into one device while maintaining high absorption, controlled polarization sensitivity, and broadband detection capability. The resulting far-infrared detectors can be read out with a simple approach, which is compact and minimizes thermal loading. The proposed microwave kinetic inductance detector (MKID) consists of three basic elements. The first is the absorptive section in which the incident power is coupled to a superconducting resonator at far-infrared frequency above its superconducting critical frequency (where superconductor becomes normal conductor). This absorber's shape effectively absorbs signals in the desired polarization state and is resonant at the radio frequency (RF) used for readout of the device. Control over the metal film used in the absorber allows realization of structures with either a 50% broadband or 100% resonance absorptance over a 30% fractional bandwidth. The second element is a microwave resonator - which is realized from the thin metal films used to make the absorber as transmission lines - whose resonance frequency changes due to a variation in its kinetic inductance. The resonator's kinetic inductance is a function of the power absorbed by the device. A low-loss dielectric (mono-crystalline silicon) is used in a parallel-plate transmission line structure to realize the desired superconducting resonators. There is negligible coupling among the adjacent elements used to define the polarization sensitivity of each detector. The final component of the device is a microwave transmission line, which is coupled to the resonator, and allows detection of changes in resonance frequency for each detector in the focal plane array. The spiral shape of the detector's absorber allows incident power with two polarizations to couple to the detector equally. A stepped impedance resonator was used that allows the incident power absorbed in the detecting membrane area to be uniformly distributed in the detector's transmission line at the RF readout frequency. This maximizes the sensitivity of the detector. The signal is read out via a frequency multiplexing technique that requires a minimum number of interface transmission lines for readout. This reduces the packaging complexity and coupling to the device's thermal environment.

Wollack, Edward↗
NASA NTRS2012

Test of Shi et al. Method to Infer the Magnetic Reconnection Geometry from Spacecraft Data: MHD Simulation with Guide Field and Antiparallel Kinetic Simulation

When analyzing data from an array of spacecraft (such as Cluster or MMS) crossing a site of magnetic reconnection, it is desirable to be able to accurately determine the orientation of the reconnection site. If the reconnection is quasi-two dimensional, there are three key directions, the direction of maximum inhomogeneity (the direction across the reconnection site), the direction of the reconnecting component of the magnetic field, and the direction of rough invariance (the "out of plane" direction). Using simulated spacecraft observations of magnetic reconnection in the geomagnetic tail, we extend our previous tests of the direction-finding method developed by Shi et al. (2005) and the method to determine the structure velocity relative to the spacecraft Vstr. These methods require data from four proximate spacecraft. We add artificial noise and calibration errors to the simulation fields, and then use the perturbed gradient of the magnetic field B and perturbed time derivative dB/dt, as described by Denton et al. (2010). Three new simulations are examined: a weakly three-dimensional, i.e., quasi-two-dimensional, MHD simulation without a guide field, a quasi-two-dimensional MHD simulation with a guide field, and a two-dimensional full dynamics kinetic simulation with inherent noise so that the apparent minimum gradient was not exactly zero, even without added artificial errors. We also examined variations of the spacecraft trajectory for the kinetic simulation. The accuracy of the directions found varied depending on the simulation and spacecraft trajectory, but all the directions could be found within about 10 for all cases. Various aspects of the method were examined, including how to choose averaging intervals and the best intervals for determining the directions and velocity. For the kinetic simulation, we also investigated in detail how the errors in the inferred gradient directions from the unmodified Shi et al. method (using the unperturbed gradient) depended on the amplitude of the calibration errors. For an accuracy of 3 for the maximum gradient direction, the calibration errors could be as large as 3% of reconnection magnetic field, while for the same accuracy for the minimum gradient direction, the calibration errors could only be as large as 0.03% of the reconnection magnetic field. These results suggest that the maximum gradient direction can normally be determined by the unmodified Shi et al. method, while the modified method or some other method must be used to accurately determine the minimum gradient direction. The structure velocity was found with magnitude accurate to 2% and direction accurate to within 5%.

Denton, R.↗
NASA NTRS2013

Kinetic and Potential Sputtering of Lunar Regolith: Contribution of Solar-Wind Heavy Ions

Sputtering of lunar regolith by protons as well as solar-wind heavy ions is considered. From preliminary measurements of H+, Ar+1, Ar+6 and Ar+9 ion sputtering of JSC-1A AGGL lunar regolith simulant at solar wind velocities, and TRIM simulations of kinetic sputtering yields, the relative contributions of kinetic and potential sputtering contributions are estimated. An 80-fold enhancement of oxygen sputtering by Ar+ over same-velocity H+, and an additional x2 increase for Ar+9 over same-velocity Ar+ was measured. This enhancement persisted to the maximum fluences investigated is approximately 1016/cm (exp2). Modeling studies including the enhanced oxygen ejection by potential sputtering due to the minority heavy ion multicharged ion solar wind component, and the kinetic sputtering contribution of all solar wind constituents, as determined from TRIM sputtering simulations, indicate an overall 35% reduction of near-surface oxygen abundance. XPS analyses of simulant samples exposed to singly and multicharged Ar ions show the characteristic signature of reduced (metallic) Fe, consistent with the preferential ejection of oxygen atoms that can occur in potential sputtering of some metal oxides.

Meyer, F. W.↗
NASA NTRS2014

Bidirectional Energy Cascades and the Origin of Kinetic Alfvenic and Whistler Turbulence in the Solar Wind

The observed steep kinetic scale turbulence spectrum in the solar wind raises the question of how that turbulence originates. Observations of keV energetic electrons during solar quiet time suggest them as a possible source of free energy to drive kinetic turbulence. Using particle-in-cell simulations, we explore how the free energy released by an electron two-stream instability drives Weibel-like electromagnetic waves that excite wave-wave interactions. Consequently, both kinetic Alfvénic and whistler turbulence are excited that evolve through inverse and forward magnetic energy cascades.

plasmas↗
NASA NTRSDate not supplied

Zircon (U-TH)/He Impact Crater Thermochronometry and the Effects of Shock Microstructures on Helium Diffusion Kinetics

Accurate age determinations of hyper-velocity impact and cratering events remains difficult and often controversial, while less than half of all known impact craters are regarded as accurately and precisely dated. Besides 40Ar/39Ar and U-Pb methods, zircon (U-Th)/He dating of impactites is a burgeoning technique to date large- to medium-sized impact structures. Zircon (U-Th)/He ages can be fully reset in minutes at 1000°C, T commonly reached in and directly adjacent to impact melt domains, whereas complete resetting of zircon (U-Th)/He at <300°C, which might be encountered near the crater margins or persist in post-impact hydrothermal systems, may take >103-4 years. However, there is a critical need to test the reliability of (U-Th)/He impact dating in shock deformed zircon, and to quantify helium diffusion kinetics in well-characterized grains with a broad spectrum of shock-induced defect substructures. For this purpose, we investigated samples from two impact structures, the 66 Ma Chicxulub multi-ring basin and the 15 Ma Ries complex crater, to compare zircon diffusion kinetics from impact structures with varying parameters, including size, age, and hydrothermal system longevity. Shock microstructures were characterized by backscattered-electron (BSE) imaging prior to determination of diffusion step-heating fractional release experiments using light-bulb furnace with prograde and retrograde incrementally 10°C steps from 300°C to 600°C. We find that zircon with low-level shock microstructures exhibit no significant deviation from helium diffusion kinetics of undamaged zircon. In contrast, zircon grains with planar deformation features and granular textures classified by SEM are characterized by a dramatic decrease in helium retentivity, similar to radiation damage, due to the reduction in the effective domain size and the introduction of fast diffusion pathways. This likely renders shocked grains more susceptible to impact-induced hydrothermal resetting. Hence, characterization of impact microstructure is critical for determining accurate impact ages, but also offers the opportunity to determine the magnitude and duration of post-impact hydrothermal circulation.

Catherine Ross↗
NASA NTRS2017

Factors Limiting Li+ Charge Transfer Kinetics in Li-ion Batteries

Understanding the factors limiting Li+ charge transfer kinetics in Li-ion batteries is essential in improving the rate performance, especially at lower temperatures. The Li+ charge transfer process involved in the lithium intercalation of graphite anode includes the step of de-solvation of the solvated Li+ in the liquid electrolyte and the step of transport of Li+ in the preformed solid electrolyte interphase (SEI) on electrodes until the Li+ accepts an electron at the electrode and becomes a Li in the electrode. Whether the de-solvation process or the Li+ transport through the SEI is a limiting step depends on the nature of the interphases at the electrode and electrolyte interfaces. Several examples involving the electrode materials such as graphite, lithium titanate (LTO), lithium iron phosphate (LFP), lithium nickel cobalt aluminum oxide (NCA) and solid Li+ conductor such as lithium lanthanum titanate or Li-Al-Ti-phosphate are reviewed and discussed to clarify the conditions at which either the de-solvation or the transport of Li+ in SEI is dominating and how the electrolyte components affect the activation energy of Li+ charge transfer kinetics. How the electrolyte additives impact the Li+ charge transfer kinetics at both the anode and the cathode has been examined at the same time in 3-electrode full cells. The resulting impact on Li+ charge transfer resistance, Rct, and activation energy, Ea, at both electrodes are reported and discussed.

Delp, Samuel A.↗
NASA NTRSDate not supplied

Zircon (U-Th)/He Impact Crater Thermochronometry and the Effects of Shock Microstructures on Helium Diffusion Kinetics

Absolute age determination of impact cratering events remains difficult and often controversial; a challenge that has resulted in < 50% of known impact craters regarded as accurately and precisely dated. Besides conventional 40Ar/39Ar and U-Pb methods, zircon (U-Th)/He (ZHe) dating of impactites has been applied to large- to medium-sized impact structures. ZHe dates can be fully reset in minutes at 1000°C, which is commonly reached in central sections of the melt sheet, whereas resetting of ZHe at <300°C, which might be encountered near the crater margins or persist in post-impact hydrothermal systems, may take >103-6 years. There is a critical need to quantify the effects of shock-induced microstructures and impact metamorphism on helium diffusion kinetics in well-characterized, variably shocked zircon to further establish the reliability of (U-Th)/He for dating impacts. For this purpose, we investigated suevite and impact melt samples from two impact structures, the Chicxulub multi-ring basin and the Ries complex crater, which enables us to compare zircon helium diffusion kinetics from impact structures with differing sizes, ages, and hydrothermal system longevities. Shock microstructures were characterized by backscattered-electron (BSE) imaging prior to diffusion step-heating fractional release experiments using light-bulb furnace with prograde and retrograde incremental 10°C steps from 250°C to 600°C. Afterward, we characterize the diffusion domain sizes and their interconnectivity within the shocked zircon grains using electron backscatter diffraction (EBSD). We find that zircon with few shock microstructures exhibit no significant deviation from helium diffusion kinetics of undamaged zircon. In contrast, zircon grains with planar deformation features (PDFs) and granular textures classified by BSE and EBSD are characterized by a dramatic decrease in helium retentivity, similar to radiation damaged grains, due to a reduction in the effective domain size and the introduction of interconnected fast diffusion pathways created by shock microstructures. A subset of grains were dated by ZHe after the external morphology of the grains was determined by BSE imaging. The euhedral grains yielded a weighted mean age within the uncertainty of the accepted impact ages, whereas the grains with PDFs or granular textures exhibited younger ages. Thus, these diffusion experiments and ZHe dates suggest that the dramatic decrease in domain size likely renders shocked grains more susceptible to impact-induced hydrothermal resetting and subsequent overprinting. Hence, characterization of shock microstructures is critical for determining accurate impact ages using ZHe methods especially when applied to previously unconstrained craters. The thermochronometer also offers the opportunity to determine the magnitude and duration of post-impact hydrothermal circulation.

zircon↗
NASA NTRS2022

After DART: Using the first full-scale test of a kinetic impactor to inform a future planetary defense mission

NASA’s Double Asteroid Redirection Test (DART) is the first full-scale test of an asteroid deflection technology. Results from the hypervelocity kinetic impact and Earth-based observations, coupled with LICIACube and the later Hera mission, will result in measurement of the momentum transfer efficiency accurate to ~10% and characterization of the Didymos binary system. But DART is a single experiment; how could these results be used in a future planetary defense necessity involving a different asteroid? We examine what aspects of Dimorphos’s response to kinetic impact will be constrained by DART results; how these constraints will help refine knowledge of the physical properties of asteroidal materials and predictive power of impact simulations; what information about a potential Earth impactor could be acquired before a deflection effort; and how design of a deflection mission should be informed by this understanding. We generalize the momentum enhancement factor β, showing that a particular direction-specific β will be directly determined by the DART results, and that a related direction-specific β is a figure of merit for a kinetic impact mission. The DART β determination constrains the ejecta momentum vector, which, with hydrodynamic simulations, constrains the physical properties of Dimorphos’s near-surface. In a hypothetical planetary defense exigency, extrapolating these constraints to a newly discovered asteroid will require Earth-based observations and benefit from in-situ reconnaissance. We show representative predictions for momentum transfer based on different levels of reconnaissance and discuss strategic targeting to optimize the deflection and reduce the risk of a counterproductive deflection in the wrong direction.

Asteroids↗
NASA NTRS2023

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗
DOE OSTI2013

Ignition Quality Tester (IQT): An Alternative for Characterizing the Combustion Kinetics of Low Volatility Fuels

The Ignition Quality Tester (IQT) is a constant volume spray combustion system that can be heated and pressurized to conditions that are similar to a diesel engine at top dead center. With no moving parts and the ability to handle low volatility fuels, the IQT can be a bridge between engines and traditional methods for studying chemical kinetics. By comparing experimental data with model predictions, the IQT has been used to validate skeletal kinetic models of ignition. CFD modeling of the IQT using KIVA-3V was used to predict ignition of n-heptane accurately. Operating the IQT in a regime where chemical kinetics dominates (long ignition delays) allowed NTC behavior to be observed for some isomers of heptane. Experimental results with the low volatility fuel heptamethylnonane also show NTC behavior. At long ignition delays, experimental results can be compared with 0-D detailed chemical mechanisms.

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

High-Pressure Turbulent Flame Speeds and Chemical Kinetics of Syngas Blends with and without Impurities

This Topical Report documents the first year of the project, from October 1, 2013 through September 30, 2014. Efforts for this project included experiments to characterize the atmospheric-pressure turbulent flame speed vessel over a range of operating conditions (fan speeds and turbulent length scales). To this end, a new LDV system was acquired and set up for the detailed characterization of the turbulence field. Much progress was made in the area of impurity kinetics, which included a numerical study of the effect of impurities such as NO2, NO, H2S, and NH3 on ignition delay times and laminar flame speeds of syngas blends at engine conditions. Experiments included a series of laminar flame speed measurements for syngas (CO/H2) blends with various levels of CH4 and C2H6 addition, and the results were compared to the chemical kinetics model of NUI Galway. Also, a final NOx kinetics mechanism including ammonia was assembled, and a journal paper was written and is now in press. Overall, three journal papers and six conference papers related to this project were published this year. Finally, much progress was made on the design of the new high-pressure turbulent flame speed facility. An overall design that includes a venting system was decided upon, and the detailed design is in progress.

Peterson, Eric↗
DOE OSTI2011

Turbulent Flame Speeds and NOx Kinetics of HHC Fuels with Contaminants and High Dilution Levels

This progress report documents the first year of the project, from October 1, 2010 through September 30, 2011. Laminar flame speeds and ignition delay times have been measured for hydrogen and various compositions of H 2 /CO (syngas) at elevated pressures and elevated temperatures. Two constant-volume cylindrical vessels were used to visualize the spherical growth of the flame through the use of a schlieren optical setup to measure the laminar flame speed of the mixture. Hydrogen experiments were performed at initial pressures up to 10 atm and initial temperatures up to 443 K. A syngas composition of 50/50 was chosen to demonstrate the effect of carbon monoxide on H 2 -O 2 chemical kinetics at standard temperature and pressures up to 10 atm. All atmospheric mixtures were diluted with standard air, while all elevated-pressure experiments were diluted with a He:O 2 of 7:1 to minimize hydrodynamic instabilities. The laminar flame speed measurements of hydrogen and syngas are compared to available literature data over a wide range of equivalence ratios where good agreement can be seen with several data sets. Additionally, an improved chemical kinetics model is shown for all conditions within the current study. The model and the data presented herein agree well, which demonstrates the continual, improved accuracy of the chemical kinetics model. A high-pressure shock tube was used to measure ignition delay times for several baseline compositions of syngas at three pressures across a wide range of temperatures. The compositions of syngas (H 2 /CO) presented in this study include 80/20, 50/50, 40/60, 20/80, and 10/90, all of which are compared to previously published ignition delay times from a hydrogen-oxygen mixture to demonstrate the effect of carbon monoxide addition. Generally, an increase in carbon monoxide increases the ignition delay time, but there does seem to be a pressure dependency. At low temperatures and pressures higher than about 12 atm, the ignition delay times appear to be indistinguishable with an increase in carbon monoxide. However, at high temperatures the composition of H 2 and CO has a strong influence on ignition delay times. Model agreement is good across the range of the study, particularly at the elevated pressures. Also an increase in carbon monoxide causes the activation energy of the mixture to decrease.

Petersen, Eric↗
DOE OSTI2012

Turbulent flame speeds and NOx kinetics of HHC fuels with contaminants and high dilution levels

This progress report documents the second year of the project, from October 1, 2011 through September 30, 2012. Characterization of the new turbulent flame speed vessel design was completed. Turbulence statistics of three impellers with different geometric features were measured using particle image velocimetry inside a Plexiglas model (~1:1 scale) of a cylindrical flame speed vessel (30.5 cm ID × 35.6 cm L). With four impellers arranged in a central-symmetric configuration, turbulence intensities between 1.2 and 1.7 m/s with negligible mean flow (0.1u´) were attained at the lowest fan speeds. Acceptable ranges for homogeneity and isotropy ratios of the velocity fields were set within a narrow bandwidth near unity (0.9-1.1). Homogeneity ratios were unaffected by changes to the impeller geometry, and the prototype with the higher number of blades caused the flow to become anisotropic. The integral length scale of the flow fields varied between 27 and 20 mm, which correlates well with those typically observed inside a gas turbine combustor. The mechanism to independently vary the intensity level and the integral length scale was established, where turbulence intensity level was dependent on the rotational speed of the fan, and the integral length scale decreased with increasing blade pitch angle. Ignition delay times of H₂/O₂ mixtures highly diluted with Ar and doped with various amounts of N₂O (100, 400, 1600, 3200 ppm) were measured in a shock tube behind reflected shock waves over a wide range of temperatures (940-1675 K). The pressure range investigated during this work (around 1.6, 13, and 30 atm) allows studying the effect of N₂O on hydrogen ignition at pressure conditions that have never been heretofore investigated. Ignition delay times were decreased when N₂O was added to the mixture only for the higher nitrous oxide concentrations, and some changes in the activation energy were also observed at 1.5 and 30 atm. When it occurred, the decrease in the ignition delay time was proportional to the amount of N₂O added and depended on pressure and temperature conditions. A detailed chemical kinetics model was developed using kinetic mechanisms from the literature. This model predicts well the experimental data obtained during this study and from the literature. The chemical analysis using this model showed that the decrease in the ignition delay time was mainly due to the reaction N₂O +M ⇌ N₂ + O +M which provides O atoms to strengthen the channel O + H₂ ⇌ OH + H. Ignition delay times have been measured behind reflected shock waves at 1.5, 12 and 30 atm for a mixture representative of a syngas produced from biomass (0.29659% CO / 0.29659% H₂ / 0.15748% CO₂ / 0.08924% CH₄ / 0.20997% H₂O / 0.95013% O₂ in 98% Ar (mol.%)) and for the same biomass-derived syngas mixture doped with 200 ppm of NH₃. The importance of the various constituents on the ignition delay time was investigated by comparing the results with data from various baseline mixtures (H₂/O₂/Ar, H₂/CO/O₂/Ar and H₂/CO/O₂/Ar with one of the other constituent of the syngas (i.e. CO₂, H₂O, CH₄ or NH₃)). The equivalence ratio was set to 0.5 during this study. Several recent detailed kinetics mechanisms from the literature were computed against these data, with fair agreement. Results showed that the mixture composition can have an important effect on the ignition delay time, with most of the effect being due to CH₄ addition through the reaction CH₄+OH⇌CH₃+H₂O. The ammonia impurity had very little effect on the ignition delay time over the range of conditions studied.

Petersen, Eric↗
DOE OSTI2017

An Evaluation of Systematic Biases of the Point Kinetics Model for Active Interrogation of Uranium

The point kinetics model relies upon a series of assumptions that simplify the process of neutron fission chain propagation. Tagged neutron interrogation of uranium using 14.1 MeV neutrons can potentially benefit from the use of the point kinetics model, but the appropriateness of the model's assumptions have yet to be validated for this application. The two assumptions considered by this work are the assumption that a single distribution can be used to describe the neutron multiplicity of a chain-starting event and the assumption that parasitic neutron capture within uranium is negligible. Monte Carlo simulations are used to evaluate the appropriateness of each assumption, as well as the consequences on the emitted doubles rate. This work demonstrates that the use of the point kinetics model for this application is possible, but corrections to the source multiplicity assumption are necessary that would likely require prior knowledge of sample isotopics.

Tweardy, Matthew↗
DOE OSTI2022

Development Of Thermodynamic and Kinetic Simulation Tools and Testing Procedures for Enhanced Durability of Concrete Containing Industrial By-Products

This project developed screening tools that enable evaluation of alternative cementitious binders that create concretes to significantly reduce energy and emissions while remaining cost competitive on both initial and long-term costs. The team began with the viewpoint that acceptance of new cementitious binder products has a substantially greater chance of successful implementation when capital investment is not excessive and the end product has customers that have experience and resources to use this. As such, a binder system that is based on portland cement with blended using industrial by-products (alternative cementitious materials) has the potential for dramatic and meaningful impact. The team has focused on developing implementable solutions in specifications and current practice. This however requires three main factors: 1) ability to screen byproducts and alternative materials for success, 2) ability to ‘treat’ materials chemically to enhance kinetics, and 3) ability to provide predictions of performance of both binders and concrete from first principles. The project developed/refined a state of the art and scientifically based screening test for SCM called the pozzolanic reactivity test. The team developed kinetic models to simulate these materials as well as experimental approaches to alter selected reactions. Simulation tools were developed that enable the performance of concrete to be predicted based on the chemistry and reactivity of the cement and alternative SCM. Specifically, this project: • enhanced the kinetic reactivity models for use in multi-scale computational programs that use thermodynamics to predict reaction products. 2 • developed scaling models to extend thermodynamic modeling to link these models with pore structure. This enables strength, transport property, and coupled transport prediction. • developed tools to predict performance of cementitious materials using the pozzolanic reactivity test and chemical composition. The predicted properties are consistent with AASHTO R101 and the CEB-FIP model code and can be measured using associated test procedures. • demonstrated mechanical and fracture based modeling tools that thermodynamic predictions and inputs to predict concrete service life. These results have been used to demonstrate the value of enabling specifications to include ASTM C 595 cement as well ASTM C150 cement. In addition, these products are being used to expedite the evaluation of alternative SCM to aide in determining which materials have potential value and what ‘compositions’ of blended cements merit further investment.

42 ENGINEERING↗
DOE OSTI2023

Learning heterogeneous reaction kinetics from X-ray movies pixel-by-pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries 1 and electrocatalysts 2. Experimental characterizations of such materials by operando microscopy produce rich image datasets 3–6, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation 7. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

Chueh, William↗