Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetically constrained models”

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 109 records · Page 6

Transition from Gaseous Compounds to Aerosols in Titan's Atmosphere

We investigate the chemical transition of simple molecules like C2H2 and HCN into aerosol particles in the context of Titan's atmosphere. Experiments that synthesize analogs (tholins) for these aerosols can help understand and constrain these polymerization mechanisms. Using information available from these experiments, we suggest chemical pathways that can link simple molecules to macromolecules, that will be the precursors to aerosol particles: polymers of acetylene and cyanoacetylene, polycyclic aromatics (PAHs), polymers of HCN and other nitriles, and polynes. Although our goal here is not to build a detailed kinetic model for this transition, we propose parameterizations to estimate the production rates of these macromolecules, their C/N and C/H ratios, and the loss of parent molecules (C2H2, HCN, HC3N and other nitriles, C6H6) from the gas phase to the haze. We use a 1-dimensional photochemical model of Titan's atmosphere to estimate the formation rate of precursors macromolecules. We find a production zone slightly lower than 200 km altitude with a total production rate of 4 x 10(exp -14) g/ sq cm s and a C/N approx. = 4. These results are compared with experimental data, and to microphysical models requirements. The Cassini/Huygens mission will bring a detailed picture of the haze distribution and properties, that will be a great challenge for our understanding of those chemical processes.

Lebonnois, Sebastien↗

Cluster/Peace Electrons Velocity Distribution Function: Modeling the Strahl in the Solar Wind

We present a study of kinetic properties of the strahl electron velocity distribution functions (VDF's) in the solar wind. These are used to investigate the pitch-angle scattering and stability of the population to interactions with electromagnetic (whistler) fluctuations. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer. Our study focuses on the mechanisms that control and regulate the pitch-angle and stability of strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the electron heat-flux and temperature anisotropy. The goal is to check whether the strahl electrons are constrained by some instability (e.g., the whistler instability), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by fitting the VDF's to a spectral spherical harmonic model from which the moments are derived directly from the model coefficients.

Figueroa-Vinas, Adolfo↗

Constraining the Low-Temperature Oxidation Mechanism of n -Hexanol through the Detection and Identification of C 6 Elusive Intermediates

Alcohol-based fuels are currently considered to be viable energy carriers for the transportation sector. Consequently, a comprehensive mechanistic understanding of the low-temperature oxidation of alcohols is essential for application in advanced low-temperature compression engines. Here, in this work, a multidimensional approach involving experimental investigations, kinetic modeling, and theoretical calculations was used to provide new insights into the low-temperature oxidation mechanism of a C 6 alcohol, n -hexanol (CH 3 (CH 2 ) 5 OH), through the detection and identification of elusive C 6 intermediates. The oxidation of n-hexanol was investigated in a jet-stirred reactor under stoichiometric conditions (ϕ = 1.0), an initial fuel concentration of 2%, a residence time of 2 s, a temperature range between 500 and 660 K, and a pressure of 700 Torr. The reactants, intermediates, and final products were detected and identified by means of molecular-beam mass spectrometry coupled with single-photon ionization employing tunable synchrotron-generated vacuum ultraviolet radiation. Chemical kinetic simulations were performed using a previously published kinetic model (Togbé et al., Energy Fuels 2010, 11, 5859−5875) to predict the reactivity of n-hexanol and elucidate the predominant formation pathways of the observed low-temperature species. Experimental photoionization efficiency curves in conjunction with ab initio calculations, enabled the identification of important low-temperature species, such as C 6 unsaturated alcohols, C 6 olefinic hydroperoxides, C 6 cyclic ethers, C 6 diones, and C 6 ketohydroperoxides. The results of this study provide valuable insight into the mechanism of the low-temperature oxidation chemistry of n -hexanol, contributing to the development of kinetic models for the low-temperature oxidation of n-hexanol and other long-chain linear alcohols.

alcohols↗

Modeling the Strahl in the Solar Wind Using CLUSTER/PEACE Electron Distribution Functions

We present a study of kinetic properties of the strahl electron velocity distribution functions (VDF's) in the solar wind. These are used to investigate the pitch-angle scattering and stability of the population to interactions with electromagnetic (whistler) fluctuations. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer. Our study focuses on the mechanisms that control and regulate the pitch-angle and stability of strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the strahl-electron density, temperature anisotropy, and electron heat-flux. The goal is to check whether the strahl electrons are constrained by some instability (e.g., the whistler instability), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by modeling the 3D-VDFs to a spectral spherical harmonic model from which the moments are obtained directly from the spectral coefficients.

Vinas, Adolfo F.↗

Reactive Transport and Mass Balance Modeling of the Stimson Sedimentary Formation and Altered Fracture Zones Constrain Diagenetic Conditions at Gale Crater, Mars

On a planet as cold and dry as present-day Mars, evidence of multiple aqueous episodes offers an intriguing view into very different past environments. Fluvial, lacustrine, and eolian depositional environments are being investigated by the Mars Science Laboratory Curiosityin Gale crater, Mars. Geochemical and mineralogical observations of these sedimentary rocks suggest diagenetic processes affected the sediments. Here, we analyze diagenesis of the Stimson formation eolian parent material, which caused loss of olivine and formation of magnetite. Additional, later alteration in fracture zones resulted in preferential dissolution of pyroxene and precipitation of secondary amorphous silica and Ca sulfate. The ability to compare the unaltered parent material with the reacted material allows constraints to be placed on the characteristics of the altering solutions. In this work we use a combination of a mass balance approach calculating the fraction of a mobile element lost or gained, τ, with fundamental geochemical kinetics and thermodynamics in the reactive transport code CrunchFlow to examine the characteristics of multiple stages of aqueous alteration at Gale crater, Mars. Our model results indicate that early diagenesis of the Stimson sedimentary formation is consistent with leaching of an eolian deposit by a near-neutral solution, and that formation of the altered fracture zones is consistent with a very acidic, high sulfate solution containing Ca, P and Si. These results indicate a range of past aqueous conditions occurring at Gale crater, Mars, with important implications for past Martian climate and environments.

Gale crater↗

Reactive Transport and Mass Balance Modeling of the Stimson Sedimentary Formation and Altered Fracture Zones Constrain Diagenetic Conditions at Gale Crater, Mars

On a planet as cold and dry as present-day Mars, evidence of multiple aqueous episodes offers an intriguing view into very different past environments. Fluvial, lacustrine, and eolian depositional environments are being investigated by the Mars Science Laboratory Curiosity in Gale crater, Mars. Geochemical and mineralogical observations of these sedimentary rocks suggest diagenetic processes affected the sediments. Here, we analyze diagenesis of the Stimson formation eolian parent material, which caused loss of olivine and formation of magnetite. Additional, later alteration in fracture zones resulted in preferential dissolution of pyroxene and precipitation of secondary amorphous silica and Ca sulfate. The ability to compare the unaltered parent material with the reacted material allows constraints to be placed on the characteristics of the altering solutions. In this work we use a combination of a mass balance approach calculating the fraction of a mobile element lost or gained, tau, with fundamental geochemical kinetics and thermodynamics in the reactive transport code CrunchFlow to examine the characteristics of multiple stages of aqueous alteration at Gale crater, Mars. Our model results indicate that early diagenesis of the Stimson sedimentary formation is consistent with leaching of an eolian deposit by a near-neutral solution, and that formation of the altered fracture zones is consistent with a very acidic, high sulfate solution containing Ca, P and Si. These results indicate a range of past aqueous conditions occurring at Gale crater, Mars, with important implications for past martian climate and environments.

Gale↗

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 ↗

Quasilinear theory: the lost ponderomotive effects and why they matter

Quasilinear theory (QLT) has been used for modeling wave–plasma interactions for decades but remains largely heuristic. Plasma inhomogeneity, ponderomotive effects, microscopic fluctuations, and collisions are not easily accommodated from first principles in QLT, and typically are ignored entirely, due to the limitations of the standard Fourier–Laplace global-mode approach. This results in inconsistencies, for example, violation of the action conservation for nonresonant waves. However, these issues can be avoided, and the theory can be substantially generalized and corrected, if QLT is formulated using more suitable analytical tools, particularly, the Weyl symbol calculus. Here, an attempt is made to deliver an accessible review of this modern formulation, provide intuitive calculations for special cases, and elaborate on the connection with the ‘oscillation-center QLT’ originally proposed by Dewar (Phys Fluids 16:1102, 1973). A Fokker–Planck equation for a ‘dressed’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, ponderomotive forces, and interactions with background fields for a generic Hamiltonian, so many known formulations of QLT for specific plasma models become corollaries of a single unifying theory. Also, waves are allowed to be off-shell (not constrained by a dispersion relation), which allows them to accommodate microscopic fluctuations. This leads to a collision integral of the Balescu–Lenard type that has all the usual properties but is not restricted to any specific plasma model. For on-shell waves, a generalized version of the classic oscillation-center QLT is obtained. Finally, combined with the wave-kinetic equation, this formulation not only conserves particles, momentum, and energy, like the classic QLT but also reinstates the action conservation for nonresonant waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Weak lensing combined with the kinetic Sunyaev–Zel’dovich effect: a study of baryonic feedback

ABSTRACT Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmology Telescope DR5 kinematic Sunyaev–Zel’dovich (kSZ) to jointly constrain cosmological and astrophysical baryonic feedback parameters using a flexible analytical model, ‘baryonification’. First, using WL only, we compare the $S_8$ constraints using baryonification to a simulation-calibrated halo model, a simulation-based emulator model, and the approach of discarding WL measurements on small angular scales. We find that model flexibility can shift the value of $S_8$ and degrade the uncertainty. The kSZ provides additional constraints on the astrophysical parameters, with the joint WL + kSZ analysis constraining $S_8=0.823^{+0.019}_{-0.020}$. We measure the suppression of the non-linear matter power spectrum using WL + kSZ and constrain a mean feedback scenario that is more extreme than the predictions from most hydrodynamical simulations. We constrain the baryon fractions and the gas mass fractions and find them to be generally lower than inferred from X-ray observations and simulation predictions. We conclude that the WL + kSZ measurements provide a new and complementary benchmark for building a coherent picture of the impact of gas around galaxies across observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Asteroidal Space Weathering: The Major Role of FeS

Space weathering (SW) effects on the lunar surface are reasonably well-understood from sample analyses, remote-sensing data, and experiments, yet our knowledge of asteroidal SW effects are far less constrained. While the same SW processes are operating on asteroids and the Moon, namely solar wind irradiation, impact vaporization and condensation, and impact melting, their relative rates and efficiencies are poorly known, as are their effects on such vastly different parent materials. Asteroidal SW models based on remote-sensing data and experiments are in wide disagreement over the dominant mechanisms involved and their kinetics. Lunar space weathering effects observed in UVVIS-NIR spectra result from surface- and volume-correlated nanophase Fe metal (npFe(sup 0)) particles. In the lunar case, it is the tiny vapor-deposited npFe(sup 0) that provides much of the spectral reddening, while the coarser (largely melt-derived) npFe(sup 0) produce lowered albedos. Nanophase FeS (npFeS) particles are expected to modify reflectance spectra in much the same way as npFe(sup 0) particles. Here we report the results of experiments designed to explore the efficiency of npFeS production via the main space weathering processes operating in the asteroid belt.

Keller, L. P.↗

Data-Driven Kinetic Reaction Networks for Separation Chemistry

Understanding complex, multistep chemical reactions at the molecular level is a major challenge whose solution would greatly benefit the design and optimization of numerous chemical processes. The separation of rare-earth (4f) and actinide (5f) elements is an example where improving our chemical understanding is important for designing and optimizing new chemistries, even with a limited number of observations. Here, in this work, we leverage data-driven artificial intelligence and machine-learning approaches to develop kinetic reaction networks that describe the liquid–liquid extraction mechanism of uranium using N,N-di-2-ethylhexyl-isobutyramide (DEHiBA). Specifically, we compare and contrast the properties of two classes of models: (1) purely data-driven models that are regularized using chemistry-agnostic, L1 regression and (2) chemistry-informed models that are regularized using relative reaction energies provided by quantum mechanical calculations. We observe that purely data-driven models are unbiased, simple, and accurate in their predictions of experimental measurements when provided with sufficient data but are difficult to fully constrain and interpret. In contrast, chemistry-informed models exhibit significantly improved chemical interpretability and consistency, providing a detailed description of the separation process while achieving high accuracy through ensemble averaging. Overall, the dominant species predicted to be extracted into the organic phase is UO 2 (NO 3 ) 2 (DEHiBA) 2 , agreeing with experimental slope analysis, thermodynamic modeling, EXAFS, and crystal structures. This work demonstrates that leveraging the fundamental structure of the problem can lead to efficient learning schemes that provide both accurate predictions and chemical insights at a low computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Crystallization by Amorphous Particle Attachment: On the Evolution of Texture

Abstract Crystallization by particle attachment (CPA) is a gradual process where each step has its own thermodynamic and kinetic constrains defining a unique pathway of crystal growth. An important example is biomineralization of calcium carbonate through amorphous precursors that are morphed into shapes and textural patterns that cannot be envisioned by the classical monomer‐by‐monomer approach. Here, a mechanistic link between the collective kinetics of mineral deposition and the emergence of crystallographic texture is established. Using the prismatic ultrastructure in bivalve shells as a model, a fundamental leap is made in the ability to analytically describe the evolution of form and texture of biological mineralized tissues and to design the structure and crystallographic properties of synthetic materials formed by CPA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insights into Constraining Rate Coefficients in Fuel Oxidation Mechanisms Using Genetic Algorithm Optimization

Accurate fuel oxidation mechanisms can enable predictive capabilities that aid in advancing combustion technologies. High-level computational kinetics can yield reasonable rate coefficients with uncertainties, in some cases, below a factor of 2. Computed rate coefficients can be constrained further by optimizing against experimental data. Here, we explore the application of genetic algorithm (GA) optimization to constrain computed rate coefficients in complex fuel oxidation mechanisms in conjunction with temperature-dependent species mole fractions from jet-stirred reactor (JSR) measurements. Cyclohexane is a model candidate for understanding the reactivity of cyclic fuels. In this work, we optimize the rate coefficients of the most recent literature cyclohexane mechanism, which incorporates theoretically computed rate coefficients for the reaction networks stemming from the first and second O 2 addition pathways, against the experimental results of two separate literature JSR studies. Optimization consistency is evaluated by carrying out three GA optimizations: fitting to the temperature-dependent species mole fractions in each JSR experiment separately and simultaneously fitting the species mole fractions in both experiments. Local sensitivity analyses are used to identify five influential low-temperature oxidation reactions for optimization. Although the three optimizations do not yield identical rate coefficients, the direction of change in all five rate coefficients is consistent among the three optimizations. Performance of the models from the three optimizations is assessed against literature ignition delay times with differences in the level of agreement observed among the different optimizations. In conclusion, comparisons are made with our recent optimization work of a cyclopentane oxidation master-equation model against time-resolved species concentrations, and insights and improvements of the strategy for constraining rate coefficients using GA optimization are discussed.

33 ADVANCED PROPULSION SYSTEMS↗

Triple Oxygen Isotope Kinetics of Sulfate From Pyrite Oxidation

Sulfates from low temperature terrestrial and marine environments are out of equilibrium with their oxygen sources. Instead, oxygen isotope compositions of sulfate imply they are dominated by kinetic signatures of oxidative and reductive processes (Killingsworth et al., 2022) that are under-constrained compared to equilibrium effects (Hayles and Killingsworth, 2022). Consequently, the specific role of pyrite oxidation in setting the oxygen isotope composition of sulfate is also under-constrained. This scenario inhibits a more nuanced understanding of small variations in the oxygen isotopes of sulfate and their implications for tracing oxygen sources and sulfur cycle transformations, particularly in low-temperature relatively sulfate-limited freshwater systems such as rivers, lakes, aquifers, and mining-affected environments. To alleviate these knowledge gaps, we conduct pyrite oxidation experiments at time scales spanning minutes to years, analyze the triple oxygen isotope compositions (∆ 17 O and δ 18 O) of sulfate products, and interpret results using molecular modeling. Our preliminary results show isotopic relationships of 1) pyrite-derived sulfate ∆ 17 O reaching above the meteoric water line as its δ 18 O becomes lighter than water, and 2) sulfate ∆ 17 O and δ 18 O extending between water and O 2 . From these results we interpret a kinetic isotope effect between sulfate and water and mixing between sulfate end members containing most of their oxygen atoms from water and from O 2 .

Oxygen isotopes↗

A Limited-Ranged Two-Phase Iron Equation of State Model

This report describes the construction of a two-phase EOS for Fe that includes both (hcp) and liquid phases, and which is meant to be valid over the restricted range: V ϵ [4 : 7] °A 3 /atom, and T ϵ [1000 : 20000] K. It was recently used as the basis for a theoretical study of liquid → solid phase transformation kinetics in Fe which appeared in the Supplemental Material for a paper describing a series of NIF shock + ramp compression experiments aimed at constraining the high-P melt curve of Fe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Markov chain Monte Carlo (MCMC) Bayesian inference approach to analyze apparent activation barriers and reaction orders from microreactor data

Statistical analysis of steady-state catalytic kinetic data is often limited by data sparsity due to the slow pace at which the data is collected. Data sparsity and limitations in statistical analysis make it difficult to differentiate between mechanistic models and catalytic sites. A Bayesian inference tool is reported for catalysis researchers to estimate error in the determination of reaction orders from steady state microreactor data. The benefits of a Bayesian inference approach are discussed, as an alternative to the more common frequentist approach. The approach incorporates prior knowledge of the system and the data collected to form an error estimate on reaction orders. We investigated the effects of three distinct data treatments—individual fitting of trials, pooled analysis, and constrained regression methods—on the precision and uncertainty of reaction order determinations. To assess the robustness of our findings, we conducted sensitivity analyses to evaluate the influence of Bayesian parameters on uncertainty estimation. Additionally, we utilized synthetic data to illustrate how data quality impacts the precision of uncertainty assessments. We show Bayesian analysis can obtain a more precise estimation of error with a sparse data set than a frequentist analysis. Finally, this work provides strong evidence that the adoption of Bayesian analysis of kinetic data may help researchers make more precise arguments as to the strength of their evidence for a particular mechanistic hypothesis, or in comparing across different catalysts.

42 ENGINEERING↗

Reionization effective likelihood from Planck 2018 data

We release relike (reionization effective likelihood), a fast and accurate effective likelihood code based on the latest Planck 2018 data that allows one to constrain any model for reionization between 6 < z < 30 using five constraints from the CMB reionization principal components (PC). We tested the code on two example models which showed excellent agreement with sampling the exact Planck likelihoods using either a simple Gaussian PC likelihood or its full kernel density estimate. This code enables a fast and consistent means for combining Planck constraints with other reionization data sets, such as kinetic Sunyaev-Zeldovich effects, line-intensity mapping, luminosity function, star formation history, quasar spectra, etc., where the redshift dependence of the ionization history is important. Since the PC technique tests any reionization history in the given range, we also derive model-independent constraints for the total Thomson optical depth τ PC = $0.0619$$^{+0.0056}_{–0.0068}$ and its 15 ≤ z ≤ 30 high redshift component τ PC (15,30) < 0.020 (95% C.L.). Furthermore, the upper limits on the high-redshift optical depth is a factor of ~3 larger than those reported in the Planck 2018 cosmological parameter paper using the FlexKnot method and we validate our results with a direct analysis of a two-step model which permits this small high-z component.

79 ASTRONOMY AND ASTROPHYSICS↗

Opaque Minerals in the Matrix of the Bishunpur (LL3.1) Chondrite: Constraints on the Chondrule Formation Environment

The chemistry and mineralogy of a group of opaque mineral assemblages in the matrix of the Bishunpur LL3.1 ordinary chondrite provide insight into the nebular environment in which they formed. The assemblages consist of a kamacite (Fe,Ni) core that is rimmed by troilite (FeS) and fayalite (Fe2Si04). Accessory phases in the rims include silica (Si02), chromite (FeCr204), whitlockite (Ca3(P04)2), maricite (FeNaP04), magnetite (Fe304), and tetrataenite (FeNi). We suggest that the metal melted in and equilibrated with an igneous chondrule under high- temperature, reducing conditions. In this environment the molten alloys incorporated varied amounts of Si, Ni, P, Cr, and Co, depending on the oxygen fugacity and temperature of the melt. Some of the metal was subsequently expelled from the chondrule interiors into the surrounding nebular gas. As the temperature dropped, the alloy solidified and volatile elements corroded the metal. The main reaction products were troilite and fayalite. Thermodynamic equilibrium calculations are used to constrain the conditions under which these two phases can form simultaneously in the solar nebula. Kinetic factors are used to place a lower limit on the formation temperature. We determine that the metal corroded between 1173 and 1261 K at a total pressure in the range of 1025.0 to 1024.1 bars and a dust/gas ratio of 302 to 355 x relative to solar composition. These conditions are consistent with our model that the metal corroded in a dust- rich region of the solar nebula that was cooling after a chondrule formation event.

Lauretta, D. S.↗