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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Convergence in simulating global soil organic carbon by structurally different models after data assimilation

Abstract Current biogeochemical models produce carbon–climate feedback projections with large uncertainties, often attributed to their structural differences when simulating soil organic carbon (SOC) dynamics worldwide. However, choices of model parameter values that quantify the strength and represent properties of different soil carbon cycle processes could also contribute to model simulation uncertainties. Here, we demonstrate the critical role of using common observational data in reducing model uncertainty in estimates of global SOC storage. Two structurally different models featuring distinctive carbon pools, decomposition kinetics, and carbon transfer pathways simulate opposite global SOC distributions with their customary parameter values yet converge to similar results after being informed by the same global SOC database using a data assimilation approach. The converged spatial SOC simulations result from similar simulations in key model components such as carbon transfer efficiency, baseline decomposition rate, and environmental effects on carbon fluxes by these two models after data assimilation. Moreover, data assimilation results suggest equally effective simulations of SOC using models following either first‐order or Michaelis–Menten kinetics at the global scale. Nevertheless, a wider range of data with high‐quality control and assurance are needed to further constrain SOC dynamics simulations and reduce unconstrained parameters. New sets of data, such as microbial genomics‐function relationships, may also suggest novel structures to account for in future model development. Overall, our results highlight the importance of observational data in informing model development and constraining model predictions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Combined density functional theory/kinetic Monte Carlo investigation of surface morphology during cycling of Li-Cu electrodes

In this work, we employed a combined density functional theory and kinetic Monte Carlo approach in order to investigate the nanoscale behavior of lithium stripping and plating on lithium-copper nanoalloys. Density functional theory calculations were used to develop an atom-by-atom understanding of the nucleation behavior of lithium on copper surfaces with different surface morphologies. We found that nanometric pits resulted in significantly more favorable lithium deposition when compared to pristine or lightly-defective Cu surfaces, which is in agreement with previous experimental studies. We then used density functional theory and constrained ab initio molecular dynamics to develop interaction and reaction parameters for a larger-scale kinetic Monte Carlo model in order to provide insight into the stripping and plating behavior of nanostructured Li-Cu nanoalloy anodes. Substantial losses in capacity and strong indications of morphological change were observed for a wide range of Cu:Li ratios over 100 cycles. We observed a complex relationship between the nanoparticle composition and the properties of the formed solid-electrolyte interphase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Theoretical and kinetic modeling study of H 2 S pyrolysis

Hydrogen sulfide pyrolysis was investigated theoretically and through chemical kinetic modeling. Reactions on the SHH potential energy surface, primarily S + H 2 (+Ar) ⇌ H 2 S (+ Ar) (R1) and S + H 2 ⇌ SH + H (R6b) were characterized by ab initio calculations. Results for k 1 were in good agreement with experiment, but deviated strongly below 2000 K from values previously used in modeling. Collider efficiencies for H 2 S, S 2 , and N 2 compared to Ar were calculated for R1. Hydrogen sulfide decomposition experiments reported in literature were re-examined in terms of an updated detailed chemical kinetic model. Concentration profiles for the atomic S at high temperature in shock tubes supported the present value of k 1 and served to constrain the rate constants for reaction of S with SH and H 2 S. To explain results from batch and flow reactors, conducted at high H 2 S concentrations in the 900–1400 K range, a very fast rate constant was required for HSS + H ⇌ SH + SH. Under dilute conditions, the gas-phase chemistry was too slow to compete and the decomposition of H 2 S was controlled by loss on the reactor surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2011

Lunar Solar Origins Exploration (LunaSOX)

The Moon offers a unique vantage point from which to investigate the Sun and its interaction via the solar wind magnetic fields, plasma, and energetic particles with the geospace system including the Moon itself. The lunar surface and exosphere provide in part a record of solar coronal plasma material input and resultant space weathering over billions of years. The structure and dynamics of solar wind interactions with the Moon provide an accessible near-Earth laboratory environment for study of general solar wind interactions with the vast multitude of airless asteroidal bodies of the inner solar system. Spacecraft in lunar orbit have the often simultaneous opportunity, except when in the Earth's magnetosphere, to make in-situ compositional measurements of the solar wind plasma and to carry out remote observations from the Moon of the solar corona, potentially enabled by lunar limb occultation of the solar disk. The LunaSOX project at NASA Goddard Space Flight Center is addressing these heliophysical science objectives from and of the Moon with support from NASA's Lunar Advanced Science and Exploration Research (LASER) program: (1) specify history of solar wind parameters at and sunward of the Moon through enhanced access (http://lunasox.gsfc.nasa.gov/) to legacy and operational mission data products from the Apollo era to the present, (2) model field and plasma interactions with the lunar surface, exosphere, and wake, as constrained by the available data, through hybrid kinetic code simulations, and (3) advance mission concepts for heliophysics from and of the Moon.

Cooper, John F.↗

DOE OSTI2020

Toward Higher Voltage Solid-State Batteries by Metastability and Kinetic Stability Design

The energy density of battery systems is limited largely by the electrochemical window of the electrolyte. Here, the combined thermodynamic and kinetic effects of mechanically induced metastability are shown to greatly widen the operational voltage window of solid-state batteries based on ceramic-sulfide electrolytes. Solid electrolyte voltage stability up to 10 V is achieved with minimal degradation, far beyond the capability of organic liquid electrolytes. Furthermore, combined experiment, ab initio computation, and theoretical modeling identify the nature of mechanically constrained Li 10 GeP 2 S 12 decomposition both within the bulk and at interfaces with cathode materials at very high voltages. Previously unclear kinetic processes are identified that, when properly implemented, can potentially allow solid-state full cells with remarkably high operational voltages.

25 ENERGY STORAGE↗

NASA NTRS2004

A practical six-degree of freedom solar sail dynamics model for optimizing solar sail trajectories with torque constraints

Controlled flight of a solar sail-propelled spacecraft ('sailcraft') is a six-degree-of-freedom dynamics problem. Current state-of-the-art tools that simulate and optimize the trajectories flown by sailcraft do not treat the full kinetic (i.e. force and torque-constrained) motion, instead treating a discrete history of commanded sail attitudes, and either neglecting the sail attitude motion over an integration timestep, or treating the attitude evolution kinematically with a spline or similar treatment. The present paper discusses an aspect of developing a next generation sailcraf trajectory designing optimization tool JPL, for NASA's Solar Sail Spaceflight Simulation Software (SS). The aspect discussed in an experimental approach to modeling full six-degree-of-freedom kinetic motion of a solar sail in a trajectory propagator. Early results from implementing this approach in a new trajectory propagation tool are given.

six-degree-of-freedom↗

DOE OSTI2022

Structural magnetic glassiness in the spin ice Dy 2 Ti 2 O 7

The origin and nature of glassy dynamics presents one of the central enigmas of condensed-matter physics across a broad range of systems ranging from window glass to spin glasses. The spin-ice compound Dy 2 Ti 2 O 7 , which is perhaps best known as hosting a three-dimensional Coulomb spin liquid with magnetically charged monopole excitations, also falls out of equilibrium at low temperature. How and why it does so remains an open question. Based on an analysis of low-temperature diffuse neutron-scattering experiments employing different cooling protocols alongside recent magnetic noise studies, combined with extensive numerical modeling, we argue that upon cooling, the spins freeze into what may be termed a “structural magnetic glass,” without an a priori need for chemical or structural disorder. Specifically, our model indicates the presence of frustration on two levels, first producing a near-degenerate constrained manifold inside which phase ordering kinetics is in turn frustrated. A remarkable feature is that monopoles act as sole annealers of the spin network and their pathways and history encode the development of glass dynamics, allowing the glass formation to be visualized. Our results suggest that spin ice Dy 2 Ti 2 O 7 provides one prototype of magnetic glass formation specifically and a setting for the study of kinetically constrained systems more generally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2022

Spectral Power-law Formation by Sequential Particle Acceleration in Multiple Flare Magnetic Islands

We present a first-principles model of pitch-angle and energy distribution function evolution as particles are sequentially accelerated by multiple flare magnetic islands. Data from magnetohydrodynamic (MHD) simulations of an eruptive flare/coronal mass ejection provide ambient conditions for the evolving particle distributions. Magnetic islands, which are created by sporadic reconnection at the self-consistently formed flare current sheet, contract and accelerate the particles. The particle distributions are evolved using rules derived in our previous work. In this investigation, we assume that a prescribed fraction of particles sequentially "hops" to another accelerator and receives an additional boost in energy and anisotropy. This sequential process generates particle number spectra that obey an approximate power law at mid-range energies and presents low- and high-energy breaks. We analyze these spectral regions as functions of the model parameters. We also present a fully analytic method for forming and interpreting such spectra, independent of the sequential acceleration model. The method requires only a few constrained physical parameters, such as the percentage of particles transferred between accelerators, the energy gain in each accelerator, and the number of accelerators visited. Our investigation seeks to bridge the gap between MHD and kinetic regimes by combining global simulations and analytic kinetic theory. The model reproduces and explains key characteristics of observed flare hard X-ray spectra as well as the underlying properties of the accelerated particles. Our analytic model provides tools to interpret high-energy observations for missions and telescopes, such as RHESSI, FOXSI, NuSTAR, Solar Orbiter, EOVSA, and future high-energy missions.

magnetic reconnection↗

DOE OSTI2024

System-size dependence of the hadronic rescattering effect at energies available at the CERN Large Hadron Collider

The first measurements of K * ( 892 ) 0 resonance production as a function of charged-particle multiplicity in Xe-Xe collisions at s N N = 5.44 TeV and p p collisions at s = 5.02 TeV using the ALICE detector are presented. The resonance is reconstructed at midrapidity ( | y | < 0.5) using the hadronic decay channel K * 0 → K ± π ∓ . Measurements of transverse-momentum integrated yield, mean transverse-momentum, nuclear modification factor of K * 0 , and yield ratios of resonance to stable hadron ( K * 0 / K ) are compared across different collision systems ( p p , p -Pb, Xe-Xe, and Pb-Pb) at similar collision energies to investigate how the production of K * 0 resonances depends on the size of the system formed in these collisions. The hadronic rescattering effect is found to be independent of the size of colliding systems and mainly driven by the produced charged-particle multiplicity, which is a proxy of the volume of produced matter at the chemical freeze-out. In addition, the production yields of K * 0 in Xe-Xe collisions are utilized to constrain the dependence of the kinetic freeze-out temperature on the system size using the hadron resonance gas–partial chemical equilibrium model. ©2024 CERN, for the ALICE Collaboration 2024 CERN

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

A hybrid-kinetic simulation tool for non-thermal warm x-ray z-pinch sources, with gas-puff and wire array exemplars

Increasing the fluence of z-pinch x-ray radiation sources above ∼ 10 keV has been a long-standing goal for scientists at Sandia National Laboratories’ Z Machine. Optimizing sources for non-thermal “cold Kα” emission in higher atomic-number materials appears to be a promising path to increase warm x-ray yield. However, this emission is generated by supra-thermal electrons, which are not treated in the magnetohydrodynamic (MHD) codes that are typically used in z-pinch source development. MHD codes do not allow for charge separation or space-charge-generated electric fields, and constrain particle kinematics to Maxwellian distributions. The kinetic codes which do accommodate discrete, non-thermal energy distributions are computationally prohibitive when modeling plasmas near solid density and when modeling/tracking higher ionization states. Thus, modeling non-thermal z-pinch sources requires a new simulation tool. In this report, we present a new hybrid modeling capability that uses the fast features of MHD-type particles to the greatest extent possible, then transitions to the slower but more complete kinetic particle treatment to correctly capture the particle energy spectra that generate non-thermal emission. This capability is founded on the fully-relativistic particle-in-cell code Chicago, which already includes fluid particle treatments. The governing equations and hybrid methodology presented here are applied in simulations of an argon gas-puff and a molybdenum wire-array to provide preliminary code validation. The argon simulation is compared to measured implosion times and yields from Jones et al., Phys. Plasmas 22, 020706 (2015). The simulated x-ray yield is within 25% of measurements and the implosion times agree within a few percent. The molybdenum wire array simulation captures the implosion timing reported in Hansen et al., Phys. Plasmas 21, 031202 (2014), but work is needed to verify the available EOS table. These exemplar simulations represents the type of non-thermal sources that will be developed using the hybrid code capability going forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2014

Towards Co-evolution of Membrane Transport and Metabolism

Protocellular boundaries were inextricably connected to the metabolism they encapsulated: to be inheritable, early metabolism must have led to an increased rate of growth and division of vesicles and, similarly, transport through vesicle boundaries must have supported the evolution of metabolism. Even though explaining how this coupling emerged and evolved in the absence of the complex machinery of modern cells is one of the key issues in studies on the origin of life, little is known about the biochemical and biophysical processes that might have been involved. This gap in our knowledge is a major impediment in efforts to construct scenarios for the origin of life and laboratory models of protocells. A combination of experimental and computational studies carried out by us and our collaborators is aimed at helping to close this gap. Properties of membranes might have contributed to the selection of RNA as an early biopolymer. A kinetic mechanism was proposed (Sacerdote & Szostak, 2005) in which ribose was supplied more quickly than other aldopentoses to primordial cells for preferential incorporation of ribonucleotides into nucleic acids. This proposal is based on a finding that ribose permeates membranes an order of magnitude faster than its diastereomers, arabinose and xylose. Our computer simulations, which yield permeation rates in excellent agreement with experiment, and kinetic modeling explain this phenomenon in terms of inter- and intramolecular interactions involving exocyclic hydroxyl groups attached to carbon atoms of the pyranose ring (Wei and Pohorille, 2009). They also constrain scenarios for the formation of the earliest nucleic acids (Wei and Pohorille, 2013). In one scenario, sugars permeate protocellular walls and subsequently are used to synthesize nucleic acids inside protocells. As long as this process proceeds at the rate faster than 6x10(exp -3)/s, ribose derivatives will be available for synthesis easier than their diastereomers. If nucleosides or their activated derivatives are synthesized outside protocells and subsequently transported across protocellular membranes the kinetic mechanism does not apply because all diastereomers, which have their sugars in the furanose rather than pyranose form, permeate the membrane at approximately the same rate. Properties of membranes might have been also coupled to metabolism involving peptides. Recently, Adamala and Szostak (2013) have shown that a dipeptide inside fatty-acid vesicles catalyzes the formation of another dipeptide that binds to vesicle walls and, by doing so, promotes their growth at the expense of other vesicles. This coupling of metabolism, permeability of vesicles and their growth is the first demonstration of evolutionary advantage imparted by small, membrane-bound peptides. Building on this work we have calculated the rate at which different blocked amino acids are delivered to a protocell for synthesis of dipeptides. We have further shown that the dipeptides are located at the water-membrane interface rather than in the center of the bilayer. On these basis it is anticipated that other dipeptides containing aromatic, but not necessarily hydrophobic amino acids (e.g. tyrosine) could have the same catalytic effects. Insight from these studies allows for estimating the rate of vesicle growth and the rates of dipeptide synthesis required to keep the system in balance. These results, in combination with our earlier studies, lead to a general scenario for evolution from membrane-bound dipeptides to ion channels in the origin of life.

Metabolism↗

DOE OSTI2022

Kinetics-Informed Neural Networks

Chemical kinetics and reaction engineering consists of the phenomenological framework for the disentanglement of reaction mechanisms, optimization of reaction performance and the rational design of chemical processes. Here, we utilize feed-forward artificial neural networks as basis functions to solve ordinary differential equations (ODEs) constrained by differential algebraic equations (DAEs) that describe microkinetic models (MKMs). We present an algebraic framework for the mathematical description and classification of reaction networks, types of elementary reaction, and chemical species. Under this framework, we demonstrate that the simultaneous training of neural nets and kinetic model parameters in a regularized multi-objective optimization setting leads to the solution of the inverse problem through the estimation of kinetic parameters from synthetic experimental data. We analyze a set of scenarios to establish the extent to which kinetic parameters can be retrieved from transient kinetic data, and assess the robustness of the methodology with respect to statistical noise. Furthermore, this approach to inverse kinetic ODEs can assist in the elucidation of reaction mechanisms based on transient data.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Nitrification is a minor source of nitrous oxide (N 2 O) in an agricultural landscape and declines with increasing management intensity

Abstract The long‐term contribution of nitrification to nitrous oxide (N 2 O) emissions from terrestrial ecosystems is poorly known and thus poorly constrained in biogeochemical models. Here, using Bayesian inference to couple 25 years of in situ N 2 O flux measurements with site‐specific Michaelis–Menten kinetics of nitrification‐derived N 2 O, we test the relative importance of nitrification‐derived N 2 O across six cropped and unmanaged ecosystems along a management intensity gradient in the U.S. Midwest. We found that the maximum potential contribution from nitrification to in situ N 2 O fluxes was 13%–17% in a conventionally fertilized annual cropping system, 27%–42% in a low‐input cover‐cropped annual cropping system, and 52%–63% in perennial systems including a late successional deciduous forest. Actual values are likely to be <10% of these values because of low N 2 O yields in cultured nitrifiers (typically 0.04%–8% of NH 3 oxidized) and competing sinks for available in situ. Most nitrification‐derived N 2 O was produced by ammonia‐oxidizing bacteria rather than archaea, who appeared responsible for no more than 30% of nitrification‐derived N 2 O production in all but one ecosystem. Although the proportion of nitrification‐derived N 2 O production was lowest in annual cropping systems, these ecosystems nevertheless produced more nitrification‐derived N 2 O (higher V max ) than perennial and successional ecosystems. We conclude that nitrification is minor relative to other sources of N 2 O in all ecosystems examined.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Data from: Nitrification is a minor source of nitrous oxide (N2O) in an agricultural landscape and declines with increasing management intensity

The long-term contribution of nitrification to nitrous oxide (N2O) emissions from terrestrial ecosystems is poorly known and thus poorly constrained in biogeochemical models. Here, using Bayesian inference to couple 25 years of in situ N2O flux measurements with site-specific Michaelis-Menten kinetics of nitrification-derived N2O, we test the relative importance of nitrification-derived N2O across six cropped and unmanaged ecosystems along a management intensity gradient in the U.S. Midwest. We found that the maximum potential contribution from nitrification to in situ N2O fluxes was 13-17% in a conventionally fertilized annual cropping system, 27-42% in a low-input cover-cropped annual cropping system, and 52-63% in perennial systems including a late successional deciduous forest. Actual values are likely to be less than 10% of these values because of low N2O yields in cultured nitrifiers (typically 0.04 to 8% of NH3 oxidized) and competing sinks for available NH4+ in situ. Most nitrification-derived N2O was produced by ammonia oxidizing bacteria (AOB) rather than archaea (AOA), who appeared responsible for no more than 30% of nitrification-derived N2O production in all but one ecosystem. Although the proportion of nitrification-derived N2O production was lowest in annual cropping systems, these ecosystems nevertheless produced more nitrification-derived N2O (higher Vmax) than perennial and successional ecosystems. We conclude that nitrification is minor relative to other sources of N2O in all ecosystems examined.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Dynamic Evolution of Mass and Physical Properties of Atmospheric Organic Aerosol under Solar Irradiance

Organic aerosol (OA) particles constitute a substantial fraction of sub-micron particulate mass in the atmosphere and play a critical role in climate system. OA undergoes dynamic aging processes in the atmosphere, with photolytic aging induced by ultraviolet solar irradiance being an important yet poorly characterized mechanism. Knowledge gaps persist regarding the role of volatility transformations during photolytic aging on the OA mass decay kinetics and the evolution of climate-relevant properties, such as hygroscopicity, hindering the model evaluation of OA spatiotemporal distributions and atmospheric budgets. In this study, we conduct isothermal photolytic aging experiments on both laboratory-generated secondary organic aerosols and ambient-collected particles from urban Atlanta, utilizing a high-sensitivity Quartz Crystal Microbalance. Our results reveal that photolytic aging reduces 40–66% of the low-volatility OA mass with lifetimes ranging from 8 to 200 hours under solar irradiance, and 44–92% of the photolytic mass loss is through slow evaporation of semi- or intermediate-volatile products, kinetically limited by their volatility. We observe up to ±50% changes in OA hygroscopicity with the transformation of fresh OA to photo-recalcitrant low-volatility products, associated with changes in oxygen-to-carbon ratio and molecular weight. A kinetic model incorporating photolytic volatility transformation provides a cohesive explanation for the observed photolysis-induced changes in mass, volatility, and hygroscopicity. Our results can help constrain model representation of the dynamic evolutions of mass and climate-relevant properties during photolytic aging processes of the ambient OA, improving our understanding of OA atmospheric behavior and climate impact.

Bai, Bin↗

NASA NTRS1992

Evolution of a terrestrial magma ocean: Thermodynamics, kinetics, rheology, convection, differentiation

The evolution of an initially totally molten magma ocean is constrained on the basis of analysis of various physical problems in the magma ocean. First of all an equilibrium thermodynamics of the magma ocean is developed in the melting temperature range. The equilibrium thermodynamical parameters are found as functions only of temperature and pressure and are used in the subsequent models of kinetics and convection. Kinematic processes determine the crystal size and also determine a non-equilibrium thermodynamics of the system. Rheology controls all dynamical regimes of the magma ocean. The thermal convection models for different rheological laws are developed for both the laminar convection and for turbulent convection in the case of equilibrium thermodynamics of the multiphase system. The evolution is estimated on the basis of all the above analysis.

Solomatov, V. S.↗

DOE OSTI2024

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

36 MATERIALS SCIENCE↗