Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic methods”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

A kinetic-theory approach to turbulent chemically reacting flows

The paper examines the mathematical and physical foundations for the kinetic theory of reactive turbulent flows, discussing the differences and relation between the kinetic and averaged equations, and comparing some solutions of the kinetic equations obtained by the Green's function method with those obtained by the approximate bimodal method. The kinetic method described consists essentially in constructing the probability density functions of the chemical species on the basis of solutions of the Langevin stochastic equation for the influence of eddies on the behavior of fluid elements. When the kinetic equations are solved for the structure of the diffusion flame established in a shear layer by the bimodal method, discontinuities in gradients of the mean concentrations at the two flame edges appear. This is a consequence of the bimodal approximation of all distribution functions by two dissimilar half-Maxwellian functions, which is a very crude approximation. These discontinuities do not appear when the solutions are constructed by the Green's function method described here.

Chung, P. M.↗
DOE OSTI2025

A kinetic-based regularization method for data science applications

We propose a physics-based regularization technique for function learning, inspired by statistical mechanics. By drawing an analogy between optimizing the parameters of an interpolator and minimizing the energy of a system, we introduce corrections that impose constraints on the lower-order moments of the data distribution. This minimizes the discrepancy between the discrete and continuum representations of the data, in turn allowing to access more favorable energy landscapes, thus improving the accuracy of the interpolator. Our approach improves performance in both interpolation and regression tasks, even in high-dimensional spaces. Unlike traditional methods, it does not require empirical parameter tuning, making it particularly effective for handling noisy data. We also show that thanks to its local nature, the method offers computational and memory efficiency advantages over Radial Basis Function interpolators, especially for large datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Scaling methods for accelerating kinetic Monte Carlo simulations of chemical reaction networks

Various kinetic Monte Carlo algorithms become inefficient when some of the population sizes in a system are large, which gives rise to a large number of reaction events per unit time. In this work, we present a new acceleration algorithm based on adaptive and heterogeneous scaling of reaction rates and stoichiometric coefficients. The algorithm is conceptually related to the commonly used idea of accelerating a stochastic simulation by considering a subvolume λΩ (0 < λ < 1) within a system of interest, which reduces the number of reaction events per unit time occurring in a simulation by a factor 1/λ at the cost of greater error in unbiased estimates of first moments and biased overestimates of second moments. Our new approach offers two unique benefits. First, scaling is adaptive and heterogeneous, which eliminates the pitfall of overaggressive scaling. Second, there is no need for an a priori classification of populations as discrete or continuous (as in a hybrid method), which is problematic when discreteness of a chemical species changes during a simulation. The method requires specification of only a single algorithmic parameter, N c , a global critical population size above which populations are effectively scaled down to increase simulation efficiency. The method, which we term partial scaling, is implemented in the open-source BioNetGen software package. We demonstrate that partial scaling can significantly accelerate simulations without significant loss of accuracy for several published models of biological systems. Finally, these models characterize activation of the mitogen-activated protein kinase ERK, prion protein aggregation, and T-cell receptor signaling.

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

Mechanistic Studies of Continuous Partial Methane Oxidation on Cu−Zeolites Using Kinetic and Spectroscopic Methods

Abstract Over the past few decades, a significant amount of research effort has focused on investigating the active site requirements and reaction mechanisms for partial methane oxidation (PMO) to methanol over copper–exchanged zeolites during stoichiometric and stepwise chemical looping routes. More recently, research efforts have expanded to include investigating the PMO reaction in a continuous catalytic regime, primarily focusing on determining the influence of catalyst composition on Cu speciation and structure and, in turn, on PMO rate and selectivity. The structures of candidate Cu active sites are commonly studied using a combination of ex situ and in situ spectroscopic approaches. In this perspective, we critically examine the prior literature on catalytic PMO over Cu–zeolites to identify key knowledge gaps that remain in our understanding as motivation for future research efforts. We identify opportunities for future research to address these gaps by adapting analogous interrogation techniques that have been successfully used to elucidate the active site requirements and mechanistic details of another catalytic redox reaction cycle on Cu–zeolites, the selective catalytic reduction (SCR) of nitrogen oxides (NO x ).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1997

Using Non-Linear Statistical Methods with Laboratory Kinetic Data

This paper will demonstrate the usefulness of standard non-linear statistical analysis on ICR and SIFT kinetic data. The specific systems used in the demonstration are the isotopic and change transfer reactions in the system of H2O+/D2O, H30+/D2O, and other permutations.

non-linear analysis kinetic isotopic permutations ↗
DOE OSTI2022

A Cartesian-diffusion Langevin method for hybrid kinetic-fluid Coulomb scattering in particle-in-cell plasma simulations

A novel, drag-diffusion Langevin method of hybrid, kinetic-fluid Coulomb scattering in plasmas is presented. Unlike previous methods, the frictional drag is always applied in the simulation frame of reference. The velocity-space diffusion is performed in the stationary-fluid frame of reference when anisotropic, and in the laboratory frame when isotropic. While the general method is mass-ratio independent, we focus on interactions of kinetic-ions and fluid-electrons to show first-order modifications to the electron velocity distribution function that are an important correction for the accurate calculation of electric resistivity. Inclusion of sub-cycling and a limit to the maximum collision frequency is shown to allow for arbitrarily large timesteps without numerical failure. Here the Langevin method is compared with a grid-based binary method and found to require a much less restrictive timestep in cases of ion–electron slowing and temperature equilibration; this finding differs from previous work and is dependent on the mass ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Optimization of Quenching Parameters for the Reduction of Titaniferous Magnetite Ore by Lean Grade Coal Using the Taguchi Method and Its Isothermal Kinetic Study

In the present study, a unique method is adopted to achieve higher reducibility of titaniferous magnetite lump ore (TMO). In this method, TMO is initially heated followed by water quenching. The quenching process generates cracks due to thermal shock in the dense TMO lumps, which, in turn, increases the extent of reduction (EOR) using the lean grade coal as a reductant. The optimum combination of parameters found by using Taguchi’s L{sub 27} orthogonal array (OA) (five factors, three levels) is − 8 + 4 mm of particle size (PS{sub 1}), 1423 K of quenching temperature (Qtemp{sub 2}), 15 minutes of quenching time (Qtime{sub 3}), 3 times the number of quenching {(No. of Q)_3}, and 120 minutes of reduction time (Rtime{sub 3}) at fixed reduction temperature of 1473 K. At optimized levels of the parameters, 92.39 pct reduction is achieved. Isothermal reduction kinetics of the quenched TMO lumps at the optimized condition reveals mixed controlled mechanisms [initially contracting geometry (CG3) followed by diffusion (D3)]. Activation energies calculated are 69.895 KJ/mole for CG3 and 39.084 KJ/mole for D3.

36 MATERIALS SCIENCE↗
NASA NTRS1977

Modelling the effects of turbulence on acoustic wave propagation

A review of recent theoretical modeling of the effects of turbulence on acoustic wave propagation is presented. The Born approximation limit of the smoothing method of Keller and of the kinetic theory method of Howe are compared with the classical perturbation results of Chernov. Secular growth is observed in the higher order terms of the kinetic theory method of Howe when multiple scaling is ignored. The forward propagation - beam broadening model of Brown and Clifford is discussed as well as the single scatter model of DeLoach and of Rudd. Each model is examined with respect to its prediction of conservation of energy.

Harris, W. L.↗
DOE OSTI2019

Iron Ore Reduction by Hydrogen Using a Laboratory Scale Fluidized Bed Reactor: Kinetic Investigation—Experimental Setup and Method for Determination

The reduction kinetics of hematite iron ore fines to metallic iron by hydrogen using a laboratory fluidized bed reactor were investigated in a temperature range between 873 K and 1073 K, by measuring the weight change of the sample portion during reduction. The fluidization conditions were checked regarding plausibility within the Grace diagram and the measured pressure drop across the material during experiments. The apparent activation energy of the reduction was determined against the degree of reduction and varied along an estimated two-peak curve between 11 and 55 kJ mol{sup −1}. Conventional kinetic analysis for the reduction of FeO to metallic iron, using typical models to describe gas–solid reactions, does not show results with high accuracy. Multistep kinetic analysis, using the Johnson–Mehl–Avrami model, shows that the initial stage of reduction from Fe{sub 2}O{sub 3} to Fe{sub 3}O{sub 4}, and partly to FeO, is controlled by diffusion and chemical reaction, depending on the temperature. Further reduction can be described by a combination of nucleation and chemical reaction, whereby the influence of nucleation increases with an increasing reduction temperature. The results of the kinetical analysis were linked to the shape of the curve from apparent activation energy against the degree of reduction.

36 MATERIALS SCIENCE↗
NASA NTRS1982

A numerical method for integrating the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes

An extension of the method of moments is developed for the numerical integration of the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes. The number density function n sub k (x,t) in each separate droplet packet between droplet mass grid points (x sub k, x sub k+1) is represented by an expansion in orthogonal polynomials with a given weighting function. In this way droplet number concentrations, liquid water contents and other moments in each droplet packet are conserved and the problem of solving the kinetic equations is replaced by one of solving a set of coupled differential equations for the number density function moments. The method is tested against analytic solutions of the corresponding kinetic equations. Numerical results are obtained for different coalescence/breakup and condensation/evaporation kernels and for different initial droplet spectra. Also droplet mass grid intervals, weighting functions, and time steps are varied.

Emukashvily, I. M.↗
NASA NTRS1987

Decoupled direct method for sensitivity analysis in combustion kinetics

An efficient, decoupled direct method for calculating the first order sensitivity coefficients of homogeneous, batch combustion kinetic rate equations is presented. In this method the ordinary differential equations for the sensitivity coefficients are solved separately from , but sequentially with, those describing the combustion chemistry. The ordinary differential equations for the thermochemical variables are solved using an efficient, implicit method (LSODE) that automatically selects the steplength and order for each solution step. The solution procedure for the sensitivity coefficients maintains accuracy and stability by using exactly the same steplengths and numerical approximations. The method computes sensitivity coefficients with respect to any combination of the initial values of the thermochemical variables and the three rate constant parameters for the chemical reactions. The method is illustrated by application to several simple problems and, where possible, comparisons are made with exact solutions and those obtained by other techniques.

Radhakrishnan, Krishnan↗
DOE OSTI2019

Rate/Concentration Kinetic Petals: A Transient Method to Examine the Interplay of Surface Reaction Processes

Transient pulse response experiments are used to construct rate/concentration (RC) kinetic dependencies—RC petals—and provide a new method to distinguish the timing and interplay of adsorption, surface reaction, and product formation on complex (industrial) materials. A petal shape arises as the dynamic “reaction-diffusion” experiment forces the concentration and reaction rate to return to zero. In contrast to the typical steady-state “Langmuir-type” RC dependence, RC petals have two branches, which arise as a result of decoupled gas and surface concentrations in the non-steady-state regime. To demonstrate this approach, the characteristics of petal shapes using ammonia decomposition as a probe reaction are presented. Ammonia, hydrogen, and nitrogen transformation rates are compared on three simple materials: Fe, Co, and a bimetallic CoFe preparation when ammonia is pulsed at 550°C in a low-pressure diffusion reactor. All materials demonstrate a two-branch kinetic RC dependence for ammonia adsorption, and rate constants are quantified in the low-coverage regime. We found that H 2 and N 2 product formation was dependent on the concentration of surface intermediates for all materials with one exception: for Co, an additional fast hydrogen generation process was observed; the rate of which coincided with ammonia adsorption. Nitrogen generation was only significant for CoFe and Co and on the CoFe catalyst, a self-inhibition property was observed. A method for estimating the number of active sites based on the RC petals is presented and was applied to the Fe and CoFe samples. The surface coverage and rate of formation/conversion of surface intermediates are interpreted from the examination of shape characteristics of the RC petals for each material.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Gas-Kinetic Theory Based Flux Splitting Method for Ideal Magnetohydrodynamics

A gas-kinetic solver is developed for the ideal magnetohydrodynamics (MHD) equations. The new scheme is based on the direct splitting of the flux function of the MHD equations with the inclusion of "particle" collisions in the transport process. Consequently, the artificial dissipation in the new scheme is much reduced in comparison with the MHD Flux Vector Splitting Scheme. At the same time, the new scheme is compared with the well-developed Roe-type MHD solver. It is concluded that the kinetic MHD scheme is more robust and efficient than the Roe- type method, and the accuracy is competitive. In this paper the general principle of splitting the macroscopic flux function based on the gas-kinetic theory is presented. The flux construction strategy may shed some light on the possible modification of AUSM- and CUSP-type schemes for the compressible Euler equations, as well as to the development of new schemes for a non-strictly hyperbolic system.

Xu, Kun↗
NASA NTRS1984

A comparison of the efficiency of numerical methods for integrating chemical kinetic rate equations

The efficiency of several algorithms used for numerical integration of stiff ordinary differential equations was compared. The methods examined included two general purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes were applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code available for the integration of combustion kinetic rate equations. It is shown that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient then evaluating the temperature by integrating its time-derivative.

Radhakrishnan, K.↗
DOE OSTI2020

A Digital Imaging Method for Evaluating the Kinetics of Vapochromic Response

This work describes the use of a cell phone camera and the L*a*b method (color space specified by the International Commission on Illumination) to characterize the color change in different vapochromic systems. In this study we have developed a semi-automatic color change analysis software that digitally analyzes images (e.g., video frames) collected while a vapochromic material is absorbing vapor. The advantages of using this method, as compared to reflectance spectroscopy or transmission spectroscopy through a thin film, include low cost, convenience, portability, ease of sample preparation, the absence of need for specialized equipment, and the ease of simultaneously collecting data on different samples under identical conditions. In addition, this method arguably provides direct insight into what a human would observe when monitoring these color changes by eye. Limitations of the method also are discussed in this paper.

cell phone camera, colorimetric sensing, vapor det↗