Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Kinetic-energy matrix elements for atomic Hylleraas-CI wave functions

Hylleraas-CI is a superposition-of-configurations method in which each configuration is constructed from a Slater-type orbital (STO) product to which is appended (linearly) at most one interelectron distance r{sub ij}. Computations of the kinetic energy for atoms by this method have been difficult due to the lack of formulas expressing these matrix elements for general angular momentum in terms of overlap and potential-energy integrals. It is shown here that a strategic application of angular-momentum theory, including the use of vector spherical harmonics, enables the reduction of all atomic kinetic-energy integrals to overlap and potential-energy matrix elements. The new formulas are validated by showing that they yield correct results for a large number of integrals published by other investigators.

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

Effective computation of stochastic protein kinetic equation by reducing stiffness via variable transformation

The stochastic protein kinetic equations can be stiff for certain parameters, which makes their numerical simulation rely on very small time step sizes, resulting in large computational cost and accumulated round-off errors. For such situation, we provide a method of reducing stiffness of the stochastic protein kinetic equation by means of a kind of variable transformation. Theoretical and numerical analysis show effectiveness of this method. Its generalization to a more general class of stochastic differential equation models is also discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

First-principles thermodynamics and defect kinetics guidelines for engineering a tailored RRAM device

Resistive Random Access Memories are among the most promising candidates for the next generation of non-volatile memory. Transition metal oxides such as HfOx and TaOx attracted a lot of attention due to their CMOS compatibility. Furthermore, these materials do not require the inclusion of extrinsic conducting defects since their operation is based on intrinsic ones (oxygen vacancies). Using Density Functional Theory, we evaluated the thermodynamics of the defects formation and the kinetics of diffusion of the conducting species active in transition metal oxide RRAM materials. The gained insights based on the thermodynamics in the Top Electrode, Insulating Matrix and Bottom Electrode and at the interfaces are used to design a proper defect reservoir, which is needed for a low-energy reliable switching device. The defect reservoir has also a direct impact on the retention of the Low Resistance State due to the resulting thermodynamic driving forces. The kinetics of the diffusing conducting defects in the Insulating Matrix determine the switching dynamics and resistance retention. The interface at the Bottom Electrode has a significant impact on the low-current operation and long endurance of the memory cell. Our first-principles findings are confirmed by experimental measurements on fabricated RRAM devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Stepwise crystallization and the layered distribution in crystallization kinetics of ultra-thin poly(ethylene terephthalate) film

Crystallization is an important property of polymeric materials. In conventional viewpoint, the transformation of disordered chains into crystals is usually a spatially homogeneous process (i.e., it occurs simultaneously throughout the sample), that is, the crystallization rate at each local position within the sample is almost the same. Here, we show that crystallization of ultra-thin poly(ethylene terephthalate) (PET) films can occur in the heterogeneous way, exhibiting a stepwise crystallization process. We found that the layered distribution of glass transition dynamics of thin film modifies the corresponding crystallization behavior, giving rise to the layered distribution of the crystallization kinetics of PET films, with an 11-nm-thick surface layer having faster crystallization rate and the underlying layer showing bulk-like behavior. The layered distribution in crystallization kinetics results in a particular stepwise crystallization behavior during heating the sample, with the two cold-crystallization temperatures separated by up to 20 K. Meanwhile, interfacial interaction is crucial for the occurrence of the heterogeneous crystallization, as the thin film crystallizes simultaneously if the interfacial interaction is relatively strong. We anticipate that this mechanism of stepwise crystallization of thin polymeric films will allow new insight into the chain organization in confined environments and permit independent manipulation of localized properties of nanomaterials.

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

Comparison of the growth kinetics of In{sub 2}O{sub 3} and Ga{sub 2}O{sub 3} and their suboxide desorption during plasma-assisted molecular beam epitaxy

We present a comprehensive study of the In{sub 2}O{sub 3} growth kinetics during plasma-assisted molecular beam epitaxy and compare it to that of the related oxide Ga{sub 2}O{sub 3} [P. Vogt and O. Bierwagen, Appl. Phys. Lett. 108, 072101 (2016)]. The growth rate and desorbing fluxes were measured during growth in-situ by a laser reflectometry set-up and line-of-sight quadrupole mass spectrometer, respectively. We extracted the In incorporation as a function of the provided In flux, different growth temperatures T{sub G}, and In-to-O flux ratios r. The data are discussed in terms of the competing formation of In{sub 2}O{sub 3} and desorption of the suboxide In{sub 2}O and O. The same three growth regimes as in the case of Ga{sub 2}O{sub 3} can be distinguished: (i) In-transport limited, O-rich (ii) In{sub 2}O-desorption limited, O-rich, and (iii) O-transport limited, In-rich. In regime (iii), In droplets are formed on the growth surface at low T{sub G}. The growth kinetics follows qualitatively that of Ga{sub 2}O{sub 3} in agreement with their common oxide and suboxide stoichiometry. The quantitative differences are mainly rationalized by the difference in In{sub 2}O and Ga{sub 2}O desorption rates and vapor pressures. For the In{sub 2}O, Ga{sub 2}O, and O desorption, we extracted the activation energies and frequency factors by means of Arrhenius-plots.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A practical approach to the sensitivity analysis for kinetic Monte Carlo simulation of heterogeneous catalysis

Lattice kinetic Monte Carlo simulations have become a vital tool for predictive quality atomistic understanding of complex surface chemical reaction kinetics over a wide range of reaction conditions. In order to expand their practical value in terms of giving guidelines for atomic level design of catalytic systems, it is very desirable to readily evaluate a sensitivity analysis for a given model. The result of such a sensitivity analysis quantitatively expresses the dependency of the turnover frequency, being the main output variable, on the rate constants entering the model. In the past the application of sensitivity analysis, such as Degree of Rate Control, has been hampered by its exuberant computational effort required to accurately sample numerical derivatives of a property that is obtained from a stochastic simulation method. Here in this study we present an efficient and robust three stage approach that is capable of reliably evaluating the sensitivity measures for stiff microkinetic models as we demonstrate using CO oxidation on RuO 2 (110) as a prototypical reaction. In a first step, we utilize the Fisher Information Matrix for filtering out elementary processes which only yield negligible sensitivity. Then we employ an estimator based on linear response theory for calculating the sensitivity measure for non-critical conditions which covers the majority of cases. Finally we adopt a method for sampling coupled finite differences for evaluating the sensitivity measure of lattice based models. This allows efficient evaluation even in critical regions near a second order phase transition that are hitherto difficult to control. The combined approach leads to significant computational savings over straightforward numerical derivatives and should aid in accelerating the nano scale design of heterogeneous catalysts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Collisionless kinetic theory of oblique tearing instabilities

The linear dispersion relation for collisionless kinetic tearing instabilities is calculated for the Harris equilibrium. In contrast to the conventional 2D geometry, which considers only modes at the center of the current sheet, modes can span the current sheet in 3D. Modes at each resonant surface have a unique angle with respect to the guide field direction. Both kinetic simulations and numerical eigenmode solutions of the linearized Vlasov-Maxwell equations have recently revealed that standard analytic theories vastly overestimate the growth rate of oblique modes. In this paper, we find that this stabilization is associated with the density-gradient-driven diamagnetic drift. The analytic theories miss this drift stabilization because the inner tearing layer broadens at oblique angles sufficiently far that the assumption of scale separation between the inner and outer regions of boundary-layer theory breaks down. The dispersion relation obtained by numerically solving a single second order differential equation is found to approximately capture the drift stabilization predicted by solutions of the full integro-differential eigenvalue problem. Finally, a simple analytic estimate for the stability criterion is provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Wave-kinetic approach to zonal-flow dynamics: Recent advances

The basic physics of drift-wave turbulence and zonal flows has long been studied within the framework of the wave-kinetic theory. Recently, this framework has been reexamined from first principles, which has led to more accurate yet still tractable “improved” wave-kinetic equations. In particular, these equations reveal an important effect of the zonal-flow “curvature” (the second radial derivative of the flow velocity) on the dynamics and stability of drift waves and zonal flows. Here, we overview these recent findings and present a consolidated high-level picture of (mostly quasilinear) zonal-flow physics within reduced models of drift-wave turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Probabilistic effects of porosity and chemical kinetics on the shock initiation of an octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) based explosive

In this work, we investigate the effects of porosity and chemical reactivity on the shock initiation of an HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine)-based polymer bonded explosive due to thin flyer impact using two-dimensional mesoscale simulations. Microstructural samples are computationally generated, and the calculated shock response is compared to experiments. The sample dimensions are 0.1×1.05 mm 2 , and each sample accounts for grain scale heterogeneities, including randomized grain morphology and pore spacing. Simulations are performed using the Lawrence Livermore National Laboratory (LLNL) arbitrary Lagrangian–Eulerian multi-physics code ALE3D. Chemistry and equation of states are handled with LLNL's thermochemical code Cheetah. We employ a modified Arrhenius-based chemical kinetics model that scales the temperature exponent fitting parameter with flyer shock conditions. This new model is calibrated to the mean response of experimental shock initiation thresholds for three separate Kapton flyer thicknesses. A simplified shock initiation criterion is proposed, which classifies the macroscale initiation response based on cross-sectional temperature profiles of the shocked sample. This criterion is used to further investigate the probabilistic response of the sample to variations in the mean porosity and chemical kinetics rate, and we find an inverse correlation between the shock pressure and the confidence interval of the experimental threshold velocity. Furthermore, the methodology presented in this work for quantifying the probabilistic nature of shock initiation is material agnostic. These mesoscale studies may be used to inform continuum scale models and expedite the development of future, tailored energetic formulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations

Traditional Kinetic Monte Carlo (KMC) approaches, rooted in Gillespie’s stochastic simulation algorithm, become computationally demanding in systems with a large range of timescales. The goal of this work is to propose and study new adaptive lattice-KMC time integration strategies for spatially non-uniform systems. To that end, two novel adaptive tau-leaping methods and their corresponding time integration strategies are developed based on the idea of the “n-fold” direct KMC method. These strategies allow for the simultaneous execution of multiple reactions, advancing time by adaptively selected coarse increments. We present numerical experiments comparing the proposed methods with existing approaches in a catalytic surface kinetics application involving ammonia decomposition.

Bimolecular reactions↗
DOE OSTI2024

Kinetic corrections to heat-flow and Nernst advection for laser heated plasmas

Reduced models for approximating the impact of kinetic electron behavior on the transport of thermal energy and magnetic field are investigated. The thermal flux limiter has improved agreement with Vlasov–Fokker–Planck data when a harmonic form is used that adjusts the electron mean free path to account for electron–electron collisions; these results apply to both unmagnetized and magnetized plasmas. Once a magnetic field is incorporated, the mean free path should also be modified using the electron gyroradius. A flux limiter on Nernst advection of magnetic fields is also required; a form that limits Nernst by the same fraction as the thermal heat-flow best reproduces kinetic simulations. A flux limiter form for the cross terms (Righi–Leduc and cross-gradient-Nernst) is also suggested. Hohlraum simulations relevant to fusion experiments on the National Ignition Facility are found to be sensitive to all of these details.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Plasma hydrodynamics from mean force kinetic theory

Mean force kinetic theory is used to evaluate the electrical conductivity, thermal conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity of a two-component plasma. The results are compared with molecular dynamics simulations. These simulations are made possible by assuming a repulsive Coulomb force for all interactions. Good agreement is found for all coefficients up to a Coulomb coupling strength of Γ ≈ 20. This is over 100-times larger than the coupling strength at which traditional theories break down. It is concluded that mean force kinetic theory is a promising framework for extending hydrodynamics to dense plasmas.

Electrical conductivity↗
DOE OSTI2020

Kinetics of the luminescence decay of Fe{sup 2+} impurity centres in polycrystalline ZnSe upon excitation by an electron beam

The kinetics of the decay of the luminescence of Fe{sup 2+} ions is measured at nitrogen temperature in polycrystalline ZnSe excited by a short pulse of accelerated electrons. The time dependence of the luminescence intensity differs from the exponential one, observed upon excitation of luminescence by a short light pulse. The obtained nonexponential dependence is theoretically described. The explanation is based on the quenching effect of the excited state of the Fe{sup 2+} ion by free electrons of the volume charge of the current of accelerated electrons in the sample (Auger effect). It was shown that the relaxation of the volume charge after the electron-accelerating voltage is removed makes a significant contribution to the decay kinetics of the impurity luminescence. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Quantitative imaging of anion exchange kinetics in halide perovskites

Significance Fundamental understanding of chemical transformation mechanism in solid-state semiconductors can offer guidelines for engineering these materials with various functionalities. Microscopically resolving such transformation process on the single-particle level can provide a direct picture of the dynamics, but has always been challenging due to the deficiencies of characterization techniques and material dimensions. This work presents a unique design on realizing the quantitative imaging of anion exchange reaction kinetics in halide perovskites by taking advantage of their high and tunable photoluminescence and characteristic lattice dynamics, coupled with controllable chemical treatments and optical microscopic techniques. Our study on kinetically visualizing ion exchange behavior in solid-state materials could also inspire mechanistic insights into other phenomena in a broad context of physical science.

2D heterostructures↗
DOE OSTI2021

Kinetic description of site ensembles on catalytic surfaces

Significance The kinetic assessment of catalytic reactions is essential to the design and application of thermochemical processes useful in petrochemical synthesis, energy conversion, and environmental remediation. We demonstrate that the ubiquitously used Langmuir–Hinshelwood formalism is incomplete and fails to correctly describe even simple catalytic reactions (e.g., A + A → A 2 ). An accurate kinetic description requires explicit rate equations for multisite configurations (e.g., A*–A* pairs) to account for the nonrandom clustering/isolation of surface species. We rigorously describe clustering/isolation phenomena by derivation of coverage-dependent stoichiometric coefficients, which capture the distinct influence each elementary step exerts on each catalytic site ensemble—thereby overcoming limitations of the mean-field (i.e., Hinshelwood) assumption in describing reactions of mobile Langmuirian adsorbates on catalytic surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electron heating in kinetic-Alfvén-wave turbulence

We report analytical and numerical investigations of subion-scale turbulence in low-beta plasmas using a rigorous reduced kinetic model. We show that efficient electron heating occurs and is primarily due to Landau damping of kinetic Alfvén waves, as opposed to Ohmic dissipation. This collisionless damping is facilitated by the local weakening of advective nonlinearities and the ensuing unimpeded phase mixing near intermittent current sheets, where free energy concentrates. The linearly damped energy of electromagnetic fluctuations at each scale explains the steepening of their energy spectrum with respect to a fluid model where such damping is excluded (i.e., a model that imposes an isothermal electron closure). In conclusion, the use of a Hermite polynomial representation to express the velocity-space dependence of the electron distribution function enables us to obtain an analytical, lowest-order solution for the Hermite moments of the distribution, which is borne out by numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Evaporation and condensation kinetics of Corundum: The origin of the 13 μm feature of oxygen-rich AGB stars

Corundum is predicted to be the most abundant refractory circumstellar dust and has been suggested as a candidate emitting the 13 μm feature in infrared spectra of oxygen-rich asymptotic giant branch stars. Crystalline dust has morphologies reflecting the anisotropic crystal structure and formation conditions or processes and may emit characteristic infrared features. In this study, evaporation and condensation experiments of corundum along different crystallographic orientations were performed in vacuum to determine key kinetic parameters controlling morphological anisotropy: evaporation and condensation coefficients (α{sub e} and α{sub c}), the degrees of kinetic hindrance on evaporation and condensation. Plausible shapes of corundum condensates were estimated from the experimental results to evaluate the infrared features of corundum condensates. The evaporation coefficients of corundum are 0.02–0.2 at 1600–1785°C, which increase with temperature and show notable anisotropy. The order of α{sub e} along the crystallographic c, a, and m axes is α {sub e}{sup m} ≫ α {sub e}{sup a} > α {sub e}{sup c} irrespective of temperature. The obtained condensation coefficients along the c, a, and m axes at 1575°C and a supersaturation ratio of ∼4 are 0.04–0.06, 0.06–0.08, and 0.1–0.2, respectively. The morphology of condensed circumstellar corundum expected from our experiments is oblate and slightly flattened to the c axis and is consistent with the fact that no presolar corundum with eccentric shapes has been found. The mass absorption coefficient of oblate corundum slightly flattened to the c axis shows a peak at 13 μm without any accompanying strong peaks. These results strongly indicate that corundum condensed anisotropically in circumstellar environments is a carrier of the unidentified 13 μm feature around oxygen-rich evolved stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Localized surface plasmon sensing based investigation of nanoscale metal oxidation kinetics

The localized surface plasmon resonance (LSPR) of nanoparticles can be a powerful and sensitive probe of chemical changes in nanoscale volumes. We used the LSPR of silver (Ag) to study the oxidation kinetics of nanoscopic volumes of cobalt (Co) metal. Bimetal nanoparticles of the immiscible Co-Ag system prepared by pulsed laser dewetting were aged in ambient air and the resulting changes to the LSPR signal and bandwidth were used to probe the oxidation kinetics. We found that Co preferentially oxidizes first. This resulted in a significant enhancement by a factor of 8 or more in the lifetime of stable Ag plasmons over that of pure Ag. Furthermore, theoretical modeling based on optical mean field approximation was able to predict the oxidation lifetimes and could help design stable Ag-based plasmonic nanoparticles for sensing applications.

2-COMPONENT COMPOSITE↗