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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A kinetic model for the transport of electrons in a graphene layer

In this article, we propose a new numerical scheme for the computation of the transport of electrons in a graphene device. The underlying quantum model for graphene is a massless Dirac equation, whose eigenvalues display a conical singularity responsible for non-adiabatic transitions between the two modes. We first derive a kinetic model which takes the form of two Boltzmann equations coupled by a collision operator modeling the non-adiabatic transitions. This collision term includes a Landau–Zener transfer term and a jump operator whose presence is essential in order to ensure a good energy conservation during the transitions. We propose an algorithmic realization of the semi-group solving the kinetic model, by a particle method. We give analytic justification of the model and propose a series of numerical experiments studying the influences of the various sources of errors between the quantum and the kinetic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Development of a point-kinetic verification scheme for nuclear reactor applications

In this paper, a new method that can be used for checking the proper implementation of time- or frequency-dependent neutron transport models and for verifying their ability to recover some basic reactor physics properties is proposed. This method makes use of the application of a stationary perturbation to the system at a given frequency and extraction of the point-kinetic component of the system response. Even for strongly heterogeneous systems for which an analytical solution does not exist, the point-kinetic component follows, as a function of frequency, a simple analytical form. The comparison between the extracted point-kinetic component and its expected analytical form provides an opportunity to verify and validate neutron transport solvers. The proposed method is tested on two diffusion-based codes, one working in the time domain and the other working in the frequency domain. As long as the applied perturbation has a non-zero reactivity effect, it is demonstrated that the method can be successfully applied to verify and validate time- or frequency-dependent neutron transport solvers. Although the method is demonstrated in the present paper in a diffusion theory framework, higher order neutron transport methods could be verified based on the same principles.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Numerical schemes for 3-wave kinetic equations: A complete treatment of the collision operator

In our previous work Walton and Tran (2023), numerical schemes for a simplified version of 3-wave kinetic equations, in which only the simple forward-cascade terms of the collision operators are kept, have been successfully designed, especially to capture the long time dynamics of the equation given the multiple blow-up time phenomenon. In this second work in the series, we propose numerical treatments for the complete 3-wave kinetic equations, in which the complete, much more complicated collision operators are fully considered based on a novel conservative form of the equation. Here we then derive an implicit finite volume scheme to solve the equation. The new discretization uses an adaptive time-stepping method which allows for the simulations to be carried to very long times. Our computed solutions are compared with previously derived long-time asymptotic estimates for the decay rate of total energy of time-dependent solutions of 3-wave kinetic equations and found to be in excellent agreement.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A conservative discontinuous Galerkin algorithm for particle kinetics on smooth manifolds

A novel, conservative discontinuous Galerkin algorithm is presented for particle kinetics on manifolds. The motion of particles on the manifold is represented using both canonical and non-canonical Hamiltonian formulations. Our schemes apply to both formulations, but the canonical formulation results in a particularly efficient scheme that also conserves particle density and energy exactly. The collisionless update is coupled to a Bhatnagar-Gross-Krook (BGK) collision operator that provides a simplified model for relaxation to local thermodynamic equilibrium. An iterative scheme is constructed to ensure collisional invariants (density, momentum and energy) are preserved numerically. Rotation of the manifold is incorporated by modifying the Hamiltonian while ensuring a canonical formulation. Several test problems, including a kinetic version of the classical Sod shock problem, Kelvin-Helmholtz instability on the surfaces of a sphere and a hyperboloid, with and without rotations, are presented. A prospectus for further development of this approach to simulation of kinetic theory in general relativity is presented.

Discontinuous Galerkin↗
DOE OSTI2011

Variation in Mutation Rates Caused by RB69pol Fidelity Mutants Can Be Rationalized on the Basis of Their Kinetic Behavior and Crystal Structures

We have previously observed that stepwise replacement of amino acid residues in the nascent base-pair binding pocket of RB69 DNA polymerase (RB69pol) with Ala or Gly expanded the space in this pocket, resulting in a progressive increase in misincorporation. However, in vivo results with similar RB69pol nascent base-pair binding pocket mutants showed that mutation rates, as determined by the T4 phage rI forward assay and rII reversion assay, were significantly lower for the RB69pol S565G/Y567A double mutant than for the Y567A single mutant, the opposite of what we would have predicted. To investigate the reasons for this unexpected result, we have determined the pre-steady-state kinetic parameters and crystal structures of relevant ternary complexes. We found that the S565G/Y567A mutant generally had greater base selectivity than the Y567A mutant and that the kinetic parameters for dNMP insertion, excision of the 3'-terminal nucleotide residue, and primer extension beyond a mispair differed not only between these two mutants but also between the two highly mutable sequences in the T4 rI complementary strand. Comparison of the crystal structures of these two mutants with correct and incorrect incoming dNTPs provides insight into the unexpected increase in the fidelity of the S565G/Y567A double mutant. Taken together, the kinetic and structural results provide a basis for integrating and interpreting in vivo and in vitro observations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Experimental and numerical investigation of hydrate-guest molecule exchange kinetics

The 2012 Ignik Sikumi #1 Field Trial in Alaska was designed to test at field scale a methane hydrate production methodology that involves injecting carbon dioxide (CO 2 ) in situ to exchange with methane (CH 4 ) within a hydrate structure, thus releasing the methane for production. Since the completion of the field trial, experimental and numerical investigations have sought to better understand observations in terms of CH 4 production and changes to the hydrate bearing formation. Collectively, these insights may lead to more effective strategies for producing CH 4 from naturally occurring hydrates. This study presents results from a laboratory experiment performed at conditions experienced during the field trial (8.27 MPa and 275.15 K) to collect data relating to the CH 4 molecular exchange process. The experiment intentionally isolated the kinetic guest molecule exchange process from additional hydrate formation and mechanical changes to the hydrate bearing sand. Data were used to inform development of numerical models that were then used to fit key hydrate parameters (bound water saturation, kinetic formation constant, kinetic exchange constant, and preferential exchange weighting factors) and analyze experimental results. Results were in agreement with observations of preferential CO 2 /CH 4 guest molecule exchange at Ignik Sikumi #1 and from previous experimental studies while providing quantitative estimates of changing hydrate compositions. Additionally, simulations confirmed hydrate behavior and composition profiles otherwise indirectly evidenced by experimental results.

03 NATURAL GAS↗
DOE OSTI2012

Computer Simulations of Interstitial Loop Growth Kinetics in irradiated bcc Fe

The growth kinetics of (001) [001] interstitial loops in bcc Fe is studied by phase-field modeling. The effect of defect (vacancy/interstitial) concentration, generation, recombination, sink strength, and elastic interaction on the growth kinetics of interstitial loops is systematically simulated. Results show that the elastic interaction between the defects and interstitial loops speeds up the growth kinetics and affects the morphology of the interstitial loops. Linear growth rate, i.e., the loop average radius is linear to time, under both aging and irradiation are predicted, which is in agreement with experimental observation. The results also show that the interstitial loop growth rate, which is directly related to the sink strength of the interstitial loop for interstitials, increases linearly with the initial interstitial concentration during aging while changing logarithmically with the interstitial generation rate under irradiation.

growth kinetics↗
DOE OSTI2011

Phosphate-intercalated Ca-Fe-layered double hydroxides: Crystal structure, bonding character, and release kinetics of phosphate

The nitrate-form of Ca-Fe-layered double hydroxide (Ca-Fe-LDH) was synthesized via co-precipitation method, and its phosphate-intercalates were prepared by ion-exchange reaction. According to X-ray diffraction analysis, the Ca-Fe-LDH-NO{sub 3}{sup -} compound and its H{sub 2}PO{sub 4}{sup -}-intercalate showed hexagonal layered structures, whereas the ion-exchange reaction with HPO{sub 4}{sup 2-} caused a frustration of the layer ordering of LDH. Fe K-edge X-ray absorption spectroscopy clearly demonstrated that the Ca-Fe-LDH lattice with trivalent iron ions was well-maintained after the ion-exchange with HPO{sub 4}{sup 2-} and H{sub 2}PO{sub 4}{sup -}. Under acidic conditions, phosphate ions were slowly released from the Ca-Fe-LDH lattice and the simultaneous release of hydroxide caused the neutralization of acidic media. Fitting analysis based on kinetic models indicated a heterogeneous diffusion process of phosphates and a distinct dependence of release rate on the charge of phosphates. This study strongly suggested that Ca-Fe-LDH is applicable as bifunctional vector for slow release of phosphate fertilizer and for the neutralization of acid soil. -- Graphical abstract: We synthesized phosphate-intercalated Ca-Fe-LDH materials that can act as bifunctional inorganic vectors for the slow release of phosphate fertilizer and also the neutralization of acid soil. Fitting analysis based on kinetic models indicated a heterogeneous diffusion process of phosphates and a distinct dependence of release rate on the charge of phosphates. Display Omitted Research Highlights: {yields} The phosphate forms of Ca-Fe-layered double hydroxide (Ca-Fe-LDH) were synthesized via co-precipitation method. The crystal structure, bonding character, and release kinetics of phosphate of the phosphate-intercalates were investigated. These Ca-Fe-LDH materials are applicable as bifunctional vector for slow release of phosphate fertilizer and for the neutralization of acid soil.

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

Isothermal nitridation kinetics of TiSi{sub 2} powders

The aim of the present work is to determine the kinetics of reaction between TiSi{sub 2} powder and gaseous nitrogen. Isothermal nitridation of TiSi{sub 2} powders with fine (1.4 µm) and medium (4.5 µm) particle size has been studied in pure nitrogen atmosphere from 1000 to 1200 °C for duration up to 50 h. The isothermal nitridation kinetics of TiSi{sub 2} powders were investigated by thermogravimetry. The nitridation rate strongly depends on the particle size and temperature. Smaller size particle exhibits higher nitridation rate due to its larger surface area. The conversion process is complex with nucleation and growth of TiN at the surface of the grain and Si{sub 3}N{sub 4} inside the grain promoted by the Kirkendall effect with an influence of the volume increase. - Graphical abstract: Backscattered electrons image of a transverse TiSi{sub 2} grain nitrurated at 1100 °C for 50 h. - Highlights: • Influence of grain size on TiSi{sub 2} powder nitridation. • Influence of temperature on TiSi{sub 2} powder nitridation. • Experimental measurements of the nitridation kinetics. • An explanation of the nitridation mechanism.

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

Oxygen exchange reaction kinetics for cerium(IV) oxide at 1000 °C

Bulk oxygen exchange rate kinetics on CeO{sub 2} at 1000 °C were observed to have a first order dependence on the fraction of reaction remaining and to be independent of oxygen partial pressure, total pressure, particle size, and specific surface area. This suggests that the exchange reaction is dominated by an internal chemical reaction that is occurring throughout the bulk of the material, and not at the material surface. Oxygen exchange rates were limited by this internal chemical reaction for all CeO{sub 2} powders studied (15 nm to −325 mesh), and had a rate constant of 1.19×10{sup −2} s{sup −1} with a time to completion of 617 s. These results are similar to the exchange rates observed previously on PuO{sub 2}, suggesting that oxygen exchange on PuO{sub 2} may also be dominated by an internal chemical reaction under similar conditions. This work will help guide future experiments on {sup 238}PuO{sub 2} oxygen exchange reactions. - Graphical abstract: Oxygen exchange kinetics on CeO{sub 2} at 1000 °C are independent of a wide range of experimental conditions and exhibit first-order chemical reaction kinetics. - Highlights: • Stable oxygen exchange rates obtained on a variety of CeO{sub 2} powders at 1000 °C. • Exchange rates are independent of atmospheric composition and specific surface area. • Exchange rates are limited by an internal chemical reaction, not a surface reaction. • CeO{sub 2} exchange rates appear similar to the rates observed on PuO{sub 2} at 1000 °C.

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

Crystallization kinetics of the Cu{sub 50}Zr{sub 50} metallic glass under isothermal conditions

Amorphous structure of the melt-spun Cu{sub 50}Zr{sub 50} amorphous alloy ribbons were confirmed by X-ray diffraction (XRD) and high-resolution transmission electron microscopy (HR-TEM). Isothermal crystallization kinetics of these alloy ribbons were investigated using differential scanning calorimetry (DSC). Besides, Arrhenius and Johnson-Mehl-Avrami (JMA) equations were utilized to obtain the isothermal crystallization kinetic parameters. As shown in the results, the local activation energy E{sub α} decreases by a large margin at the crystallized volume fraction α<0.1, which proves that crystallization process is increasingly easy. In addition, the local activation energy E{sub α} is basically constant at 0.1<α<0.9. Therefore, it turns out that the unchanged barrier is overcome in the crystallization process. Finally, E{sub α} rapidly decreases at 0.9<α<1, implying that crystallization becomes easier and easier to proceed. Nucleation activation energy E{sub nucleation} is greater than growth activation energy E{sub growth}, so nucleation is harder than growth in isothermal process. In terms of the local Avrami exponent n(α), it ranges 1.1–7.4, revealing that isothermal crystallization mechanism is interface-controlled one- two- or three-dimensional growth with different nucleation rates. - Graphical abstract: The local Avrami exponent n(α), it ranges 1.1–7.4, revealing that isothermal crystallization mechanism is interface-controlled one- two- or three-dimensional growth with different nucleation rates. - Highlights: • Isothermal crystallization kinetics of Cu{sub 50}Zr{sub 50} metallic glass was investigated. • The relationship between the local activation energy E{sub α} and the crystallized volume fraction α were determined. • The nucleation activation energy E{sub nucleation} and grain growth activation energy E{sub growth} were obtained. • The local Avrami exponent n(α) was calculated in isothermal model.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mechanism and kinetic study of rare earth extraction from allanite by direct acid leaching

Direct acid leaching was conducted to recover rare earth elements (REEs) from an allanite ore assaying 9,861.4 ppm of REEs. The effects of operational variables such as acid type, sulfuric acid (H 2 SO 4 ) concentration, temperature, solid/liquid (S/L) ratio, and particle size on the total REE (TREE) recovery were systematically examined. The REE mineralogy and the involved leaching mechanisms were investigated by microscopic characterization and residue re-leaching tests. It was found that approximately 80% of REEs was extracted using 1 M H 2 SO 4 at 75 °C for 2 h. Increasing temperature significantly improved the REE extraction, while changing acid type, H 2 SO 4 concentration, and particle size did not. The TREE recovery also gradually decreased at higher S/L ratios. The low fraction of Si leached (<4.0%) suggested that the extraction did not significantly alter the crystal structure of the dominant mineral. The leaching kinetics of REEs were fast within the first 10 min, followed by a much slower stage (i.e., 10-120 min). The kinetic modeling results showed that each stage may be controlled by both chemical reaction and diffusion (i.e., mix control). The activation energies for REE leaching using 1 M H 2 SO 4 in the 0-10 min and 10-120 min were determined to be 20.3 and 10.8 kJ/mol, respectively. Finally, the residue re-leaching results suggested that the fast leaching kinetics at the beginning may be explained by the preferential decomposition of metamict allanite rather than the formation of silica gel, while the following slower leaching stage was due to the dissolution of difficult-to-leach REE-bearing minerals such as well-crystallized allanite.

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

Selective dopant segregation modulates mesoscale reaction kinetics in layered transition metal oxide

Incorporation of foreign elements into the cathode material is broadly adopted by both academia and industry to improve the battery performance. The lack of an in-depth understanding for the underlying mechanism, however, makes it a largely try-and-error process with unsatisfactory efficiency and effectiveness. This is particularly true for the electrochemical reaction kinetics that is heterogeneous over a broad range of length scales and is determined collectively by the cathode’s electronic structure, lattice configuration, and micro-morphology. In this work, we unveiled a facet-dependent dopant segregation effect in Zr-modified single-crystal LiNi 0.6 Co 0.2 Mn 0.2 O 2 cathode. By forming kinetically favored corners on the cathode particles, the presence of a trace amount of Zr critically modulates the mesoscale reaction kinetics. Our findings suggest that a delicately controlled dopant distribution is a viable strategy for designing the next-generation battery cathode with superior structural and chemical robustness.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Simulation of edge localized mode heat pulse using drift-kinetic ions and Boltzmann electrons

This article reports on 1D+2V heat pulse propagation studies using the COGENT guiding center kinetic code. The model uses magnetized kinetic ions and a simple Boltzmann electron model. Results agree with previous kinetic and fluid modeling benchmark studies that correspond to the parameters of edge localized modes (ELMs) observed on the JET tokamak. The plasma parameters for the edge pedestal and ensuing ELM dynamics are in the low collisionality regime. Hence, the dominant balance between the assumed Maxwellian ELM source and collisionless parallel advection causes the ion PDF to develop a significantly anisotropic velocity distribution. Adding nonlinear Coulomb ion-ion collisions to the model acts to smooth the sharp features of the ion distribution function, but the anisotropy remains robust due to the low collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Combustion chemistry in the twenty-first century: Developing theory-informed chemical kinetics models

Over the past 20 to 25 years theoretical chemistry (particularly theoretical chemical kinetics) has played an increasingly important role in developing chemical kinetics models for combustion. Theoretical methods of obtaining rate parameters are now competitive in accuracy with experiment, particularly for small molecules. Moreover, theoretical methods can deal with conditions that experiments frequently cannot. In addition to increased accuracy, theory has rejuvenated methods and discovered phenomena that were completely unappreciated, or at least underappreciated, in the 20th century. Our primary interest here is in molecular-level issues, i.e. in calculating rate and transport parameters. However, dealing with kinetics models that involve thousands of reactions and hundreds of species is important for practical applications and is relatively new to the 21st century. Theory, in a general sense, and theoretical methods development have a role to play here too. We discuss in this review all these topics in some detail with an emphasis on issues and methods that have emerged in the last 20 years or so. Even so, our review is selective, rather than comprehensive, out of necessity.

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

Kinetic modeling of ignition in miniature shock tube

The ignition delay times of toluene, n-heptane/iso-octane primary reference fuel blends, and toluene/n-heptane/iso-octane toluene primary reference fuel blends were studied in the miniature high repetition rate shock tube (HRRST). In the HRRST, the temperature and pressure change during the test time, complicating interpretation of ignition delay times. Two methods of chemical kinetic modeling were employed in order to account for the changing thermodynamic states in the HRRST. One was based on direct modeling using the pressure profile and the kinetic mechanism, and the other was based on a Livengood-Wu type approach using modeled constant volume ignition delay times. Two different kinetic mechanisms were used in modeling the data. Both methods showed agreement with experimental data over a range of temperatures, pressures, and blends. While overall agreement was fair, agreement was poorer for one of the mechanisms for the primary reference fuel blends. An evaluation of the ignition delay iso-contours from these mechanisms for representative HRRST experiment thermodynamic trajectories revealed that the HRRST measured IDTs are controlled by the high temperature chemistry and are more sensitive to temperature. Furthermore, the different performance of mechanisms is therefore attributed to their different high temperature reactivity.

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

An improved semi-global intrinsic kinetics model for high temperature carbon oxidation

Measurements of the oxidation rates of various forms of carbon (soot, graphite, coal char) have often shown an unexplained attenuation with increasing temperatures in the vicinity of 2000 K, even when accounting for diffusional transport limitations and gas-phase chemical effects (e.g. CO 2 dissociation). With the development of oxy-fuel combustion approaches for pulverized coal utilization with carbon capture, high particle temperatures are readily achieved in sufficiently oxygen-enriched environments. Here, in this work, a new semi-global intrinsic kinetics model for high temperature carbon oxidation is created by starting with a previously developed 5-step mechanism that was shown to reproduce all major known trends in carbon oxidation, except for its high temperature kinetic falloff, and incorporating a recently discovered surface oxide decomposition step. The predictions of this new model are benchmarked by deploying the kinetic model in a steady-state reacting particle code (SKIPPY) and comparing the simulated results against a carefully measured set of pulverized coal char combustion temperature measurements over a wide range of oxygen concentrations in N 2 and CO 2 environments. The results show that the inclusion of the spontaneous surface oxide decomposition reaction step significantly improves predictions at high particle temperatures. Furthermore, the simulations reveal that O atoms released from the oxide decomposition step enhance the radical pool in the near-surface region and within the particle interior itself. Incorporation of literature rates for O and OH reactions with the carbon surface results in a reduction in the predicted radical pool concentrations and a very minor enhancement of the overall carbon oxidation rate.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

On the prediction of pressure effects for the combination kinetics of two alkyl radicals with the geometric mean rule

Here, alkyl combination reactions, along with their reverse alkane unimolecular decomposition reactions, are widely present in pyrolysis and combustion systems. These reactions occur mainly at temperatures at which their rate coefficients depend on pressure, and only a few studies have explored this dependence for large size alkyls. The accurate calculation of these kinetic data is limited by the size of the alkyl fragments involved, since the theoretical methods used are computationally expensive. Therefore, little reliable data is available for the modeling of real fuels alkane decomposition kinetics. In this paper, we propose an alternative approach that allows an accurate estimation of pressure-dependent kinetic rates for the combination of long-chain alkyl radicals, based on tabulated reaction models. The rates of combination of linear alkyls radicals forming ethane to pentane were calculated using the CASPT2-F12/cc-pVDZ-F12 level within VRC-TST theory combined with master-equation simulations. The computed high-pressure limit rate coefficients showed that n-propyl could be a minimum-size representative structure for larger n-alkyl chain, but this approach fails to capture the pressure dependence of the reaction that remains size-dependent for larger alkyls. The Geometric Mean Rule (GMR), a simple rule that links the self and cross-reactions of three alkyl radicals, is used to address this issue. Our theoretical calculations show for the first time that the GMR applies not only at the high-pressure limit, but also for every pressure explored. The GMR is shown to be accurate, with a maximum 22 % deviation in the falloff region for the estimation of rates for small alkyl combinations. Using this method, we propose a tabulated reaction rate rule for the combination of alkyl radicals up to C5+C5. The use of computed rate coefficients with the GMR approximation opens ways to develop pressure-dependent reaction rate rules for large n-alkyl combinations, based on accurate and computationally expensive theoretical calculations on smaller alkyl fragments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗