Search NASA⌕ Search

SEARCH · Search NASA

Results for “differential 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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Stoichiometrically-informed symbolic regression for extracting chemical reaction mechanisms from data

A data-driven computational method is introduced to extract chemical reaction mechanisms from time series chemical concentration data. It is realized through the use of dynamic symbolic regression in which a sparse analytical form for a dynamical system is discoverable from the underlying data. We specifically develop the stoichiometrically-informed symbolic regression (SISR) method to address a standing challenge in complex chemical reaction networks: given a time-series dataset of concentrations of several components, what is the mechanism and the associated rate constants? SISR finds the optimal mechanism, kinetic equations and rate constants by combining differential optimization with a genetic optimization approach that searches a symbolic space of possible reaction mechanisms. Use of SISR in several paradigmatic examples spanning linear and nonlinear reaction schemes results in excellent agreement between true and predicted mechanisms, including when the method is applied to noisy data. The advantages of a stoichiometrically-informed approach such as SISR to address reaction discovery is illustrated through comparison with the use of generic state-of-the-art data-driven approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Catalyst: Fast and flexible modeling of reaction networks

We introduce Catalyst.jl, a flexible and feature-filled Julia library for modeling and high-performance simulation of chemical reaction networks (CRNs). Catalyst supports simulating stochastic chemical kinetics (jump process), chemical Langevin equation (stochastic differential equation), and reaction rate equation (ordinary differential equation) representations for CRNs. Through comprehensive benchmarks, we demonstrate that Catalyst simulation runtimes are often one to two orders of magnitude faster than other popular tools. More broadly, Catalyst acts as both a domain-specific language and an intermediate representation for symbolically encoding CRN models as Julia-native objects. This enables a pipeline of symbolically specifying, analyzing, and modifying CRNs; converting Catalyst models to symbolic representations of concrete mathematical models; and generating compiled code for numerical solvers. Leveraging ModelingToolkit.jl and Symbolics.jl, Catalyst models can be analyzed, simplified, and compiled into optimized representations for use in numerical solvers. Finally, we demonstrate Catalyst’s broad extensibility and composability by highlighting how it can compose with a variety of Julia libraries, and how existing open-source biological modeling projects have extended its intermediate representation.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2000

The Effects of Thermal History on Nucleation of Tetragonal Lysozyme Crystals, or Hot Protein and Cold Nucleation

Chicken egg white lysozyme has a well characterized thermally driven phase transition. Between pH 4.2 and 5.2, the transition temperature, as defined by the point where the tetragonal and orthorhombic solubilities are equal, is a function of the pH, salt (precipitant) type and concentration, and most likely of the buffer concentration as well. This phase transition can be carried out with protein solution alone, prior to addition of precipitant solution. Warming a lysozyme solution above the phase transition point, then cooling it back below this point, has been shown to affect the subsequent nucleation rate, as determined by the numbers and size of crystals formed, but not the growth rate for the tetragonal crystal form . We have now measured the kinetics of this process and investigated its reversibility. The transition effects are progressive with temperature, having a half time of about 1 hour at 37C at pH 4.8. After holding a lysozyme solution at 37C (prior to addition of precipitant) for 16 hours, then cooling it back to 4C no return to the pre-warmed nucleation kinetics are observed after at least 4 weeks. Orthorhombic lysozyme crystals apparently do not undergo the flow-induced growth cessation of tetragonal lysozyme crystals. Putting the protein in the orthorhombic form does not affect the averaged face growth kinetics, only nucleation, for tetragonal crystals. This differential behaviour may be exploited to elucidate how and where flow affects the lysozyme crystal growth process. The presentation will focus on the results of these and ongoing studies in this area.

Burke, Michael↗
DOE OSTI2020

Transcriptomic and proteomic signatures of stemness and differentiation in the colon crypt

The intestinal epithelium is replaced weekly by non-quiescent stem cells with kinetics that rely on a rapid loss of stemness and choice for secretory or absorptive lineage differentiation. To determine how the cellular transcriptome and proteome changes during these transitions, we developed a new cell sorting method to purify stem cells, secretory and absorptive progenitor cells, and mature, differentiated cells. Transcriptome analyses revealed that as stem cells transition to the progenitor stage, alternative mRNA splicing and polyadenylation dominate changes in the transcriptome. In contrast, as progenitors differentiate into mature cell types, alterations in gene expression and mRNA levels drive the changes. RNA processing targets mRNAs encoding regulators of cell cycle, RNA regulators, cell adhesion, SUMOylation, and Wnt and Notch signaling. Additionally, global proteome profiling of sorted cell populations detected >2,800 proteins and revealed RNA:protein patterns of abundance and correlation. Paired together, these data highlight new potentials for autocrine and feedback regulation and provide new insights into cell state transitions in the crypt.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1985

New integration techniques for chemical kinetic rate equations. 2: Accuracy comparison

A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction, heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

Radhakrishnan, K.↗
NASA NTRS1986

New integration techniques for chemical kinetic rate equations. II - Accuracy comparison

A comparison of the accuracy of several techniques recently developed for solving stiff differential equations is presented. The techniques examined include two general purpose codes EEPISODE and LSODE developed for an arbitrary system of ordinary differential equations, and three specialized codes CHEMEQ, CREKID, and GCKP84 developed specifically to solve chemical kinetic rate equations. The accuracy comparisons are made by applying these solution procedures to two practical combustion kinetics problems. Both problems describe adiabatic, homogeneous, gas phase chemical reactions at constant pressure, and include all three combustion regimes: induction heat release, and equilibration. The comparisons show that LSODE is the most efficient code - in the sense that it requires the least computational work to attain a specified accuracy level. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for the temperature can be more accurate and efficient than computing the temperature by integrating its time derivative.

Radhakrishnan, K.↗
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 NTRS1988

Kinetic theory for electrostatic waves due to transverse velocity shears

A kinetic theory in the form of an integral equation is provided to study the electrostatic oscillations in a collisionless plasma immersed in a uniform magnetic field and a nonuniform transverse electric field. In the low temperature limit the dispersion differential equation is recovered for the transverse Kelvin-Helmholtz modes for arbitrary values of K parallel, where K parallel is the component of the wave vector in the direction of the external magnetic field assumed in the z direction. For higher temperatures the ion-cyclotron-like modes described earlier in the literature by Ganguli, Lee and Plamadesso are recovered. In this article, the integral equation is reduced to a second-order differential equation and a study is made of the kinetic Kelvin-Helmholtz and ion-cyclotron-like modes that constitute the two branches of oscillation in a magnetized plasma including a transverse inhomogeneous dc electric field.

Ganguli, G.↗
DOE OSTI2022

Kinetic Probes of the Origin of Activity in MOF-Based C–H Oxidation Catalysis

The development of homogeneous catalysis is enabled by the availability of a rich toolkit of kinetics experiments, such as the Hgdrop test, that differentiate catalytic activity at ligand-supported metal complexes from potential heterogeneous catalysts derived from decomposition of molecular species. Metal-organic frameworks (MOFs) have garnered significant attention as platforms for catalysis at site-isolated, interstitial catalyst sites. Unlike homogeneous catalysis, relatively few strategies have been advanced to evaluate the origin of catalytic activity in MOF-promoted reactions. Many of the MOFs that have been examined as potential catalysts are comprised of molecular constituents that represent viable catalysts in the absence of the extended MOF lattice, and thus interfacial sites and leached homogeneous species represent potential sources of catalyst activity. Here, we demonstrate that analysis of deuterium kinetic isotope effects (KIEs) and olefin epoxidation diastereoselectivity provides probes of the origin of catalytic activity in MOF-promoted oxidation reactions. These analyses support involvement of lattice-based Fe sites in the turnover-limiting step of C–H activation with Fe-MOF-74-based materials (i.e., the MOF functions as a bona fide catalyst) and evidence that Cu 2 -based MOF MIL-125-Cu 2 O 2 functions as a solid-state initiator for solution-phase oxidation chemistry and is not involved in the turnover limiting step (i.e., the MOF does not function as a catalyst for substrate functionalization). furthermore, we anticipate that the simple experiments described here will provide a valuable tool for clarifying the role of MOFs in C–H oxidation reactions.

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

Delayed expression of hpS2 and prolonged expression of CIP1/WAF1/SDI1 in human tumour cells irradiated with X-rays, fission neutrons or 1 GeV/nucleon Fe ions

PURPOSE: Differences in gene expression underlie the phenotypic differences between irradiated and unirradiated cells. The goal was to identify late-transcribed genes following irradiations differing in quality, and to determine the RBE of 1 GeV/n Fe ions. MATERIALS AND METHODS: Clonogenic assay was used to determine the RBE of Fe ions. Differential hybridization to cDNA target clones was used to detect differences in expression of corresponding genes in mRNA samples isolated from MCF7 cells irradiated with iso-survival doses of Fe ions (0 or 2.5 Gy) or fission neutrons (0 or 1.2 Gy) 7 days earlier. Northern analysis was used to confirm differential expression of cDNA-specific mRNA and to examine expression kinetics up to 2 weeks after irradiation. RESULTS: Fe ion RBE values were between 2.2 and 2.6 in the lines examined. Two of 17 differentially expressed cDNA clones were characterized. hpS2 mRNA was elevated from 1 to 14 days after irradiation, whereas CIP1/WAF1/SDI1 remained elevated from 3 h to 14 days after irradiation. Induction of hpS2 mRNA by irradiation was independent of p53, whereas induction of CIP1/WAF1/SDI1 was observed only in wild-type p53 lines. CONCLUSIONS: A set of coordinately regulated genes, some of which are independent of p53, is associated with change in gene expression during the first 2 weeks post-irradiation.

NASA Discipline Radiation Health↗
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.↗
NASA NTRS1984

Fast algorithms for combustion kinetics calculations: A comparison

To identify the fastest algorithm currently available for the numerical integration of chemical kinetic rate equations, several algorithms were examined. Findings to date are summarized. The algorithms examined include two general-purpose codes EPISODE and LSODE and three special-purpose (for chemical kinetic calculations) codes CHEMEQ, CRK1D, and GCKP84. In addition, an explicit Runge-Kutta-Merson differential equation solver (IMSL Routine DASCRU) is used to illustrate the problems associated with integrating chemical kinetic rate equations by a classical method. Algorithms were applied to two test problems drawn from combustion kinetics. These problems included all three combustion regimes: induction, heat release and equilibration. Variations of the temperature and species mole fraction are given with time for test problems 1 and 2, respectively. Both test problems were integrated over a time interval of 1 ms in order to obtain near-equilibration of all species and temperature. Of the codes examined in this study, only CREK1D and GCDP84 were written explicitly for integrating exothermic, non-isothermal combustion rate equations. These therefore have built-in procedures for calculating the temperature.

Radhakrishnan, K.↗
DOE OSTI2026

A projection method for particle resampling

Particle discretizations of partial differential equations are advantageous for high-dimensional kinetic models in phase-space due to their better scalability than continuum approaches with respect to dimension. Complex processes collectively referred to as particle noise hamper long time simulations with particle methods. One approach to address this problem is particle mesh adaptivity, or remapping, known as particle resampling and remeshing. Here, this work introduces a resampling method that projects particles to and from a (finite element) function space. The method is simple, using standard sparse linear algebra and finite element techniques, and it preserves all moments up to the order of a polynomial represented exactly by the continuum function space. It is distinguished from most other mesh-based methods in that new particle positions and number are decoupled from the mesh, allowing particle and continuum meshes to be adapted relatively independently. While this work is developed with structured particle and continuum phase-space grids on 1X + 1V Vlasov-Poisson models of Landau damping and two-stream instability, the method is well-suited to unstructured grids. Stable long time dynamics are demonstrated up to time T = 500. Reproducibility artifacts and data are publicly available.

Kinetic methods↗
DOE OSTI2024

Ionic Peltier effect in Li-ion electrolytes

The coupled transport of charge and heat provide fundamental insights into the microscopic thermodynamics and kinetics of materials. We describe a sensitive ac differential resistance bridge that enables measurements of the temperature difference on two sides of a coin cell with a resolution of better than 10 μK. We use this temperature difference metrology to determine the ionic Peltier coefficients of symmetric Li-ion electrochemical cells as a function of Li salt concentration, solvent composition, electrode material, and temperature. The Peltier coefficients Π are negative, i.e., heat flows in the direction opposite to the drift of Li ions in the applied electric field, large, –Π > 30 kJ mol –1 , and increase with increasing temperature at T > 300 K. The Peltier coefficient is approximately constant on time scales that span the characteristic time for mass diffusion across the thickness of the electrolyte, suggesting that heat of transport plays a minor role in comparison to the changes in partial molar entropy of Li at the interface between the electrode and electrolyte. Here, our work demonstrates a new platform for studying the non-equilibrium thermodynamics of electrochemical cells and provides a window into the transport properties of electrochemical materials through measurements of temperature differences and heat currents that complement traditional measurements of voltages and charge currents.

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

Combustion of hydrogen-air jets in local chemical equilibrium: A guide to the CHARNAL computer program

A guide to a computer program, written in FORTRAN 4, for predicting the flow properties of turbulent mixing with combustion of a circular jet of hydrogen into a co-flowing stream of air is presented. The program, which is based upon the Imperial College group's PASSA series, solves differential equations for diffusion and dissipation of turbulent kinetic energy and also of the R.M.S. fluctuation of hydrogen concentration. The effective turbulent viscosity for use in the shear stress equation is computed. Chemical equilibrium is assumed throughout the flow.

Spalding, D. B.↗
NASA NTRS1977

Neutrally stable atmospheric flow over a two-dimensional rectangular block

The phenomena of atmospheric flow over a two dimensional surface obstruction such as a building modeled as a rectangular block are analyzed by an approach using the Navier-Stokes equations with a two equation model of turbulence. The partial differential equations for the vorticity, stream function, turbulence kinetic energy, and turbulence length scale are solved by a finite difference technique. The predicted results are in agreement with the limited experimental data available. Current computed results show that the separation bubble originates from the upper front corner of the block and extends approximately 11.5 block heights behind the block. The decay of the mean velocity along the wake center line coincides almost perfectly with the experimental data. The vertical profiles of the mean velocity defect are also in reasonable agreement with wind tunnel results. Velocity profiles in the mixing region are shown to agree with the error function profile typically found in the shear layer. Details of the behavior of the turbulence kinetic energy and the turbulence length scale are also discussed.

Shieh, C. F.↗
NASA NTRS1979

Origin of cosmic magnetic fields

Reasons for describing cosmic hydromagnetic phenomena by introducing electric currents are summarized demonstrating that change from the traditional magnetic field description to a current description provides new aspects of cosmic electrodynamics. It is shown that the kink-instability of electric currents is the basic physical phenomenon responsible for magnetic flux generation in cosmic physics. A model which starts with a poloidal magnetic field is discussed; differential motions change the field configuration so that kinetic energy is transferred into magnetic energy of induced toroidal fields; the currents associated with these become unstable when the toroidal energy exceeds the poloidal energy resulting in an amplification of the original poloidal field.

Alfven, H.↗
NASA NTRS1981

Numerical solutions of atmospheric flow over semielliptical simulated hills

Atmospheric motion over obstacles on plane surfaces to compute simulated wind fields over terrain features was studied. Semielliptical, two dimensional geometry and numerical simulation of flow over rectangular geometries is also discussed. The partial differential equations for the vorticity, stream function, turbulence kinetic energy, and turbulence length scale were solved by a finite difference technique. The mechanism of flow separation induced by a semiellipse is the same as flow over a gradually sloping surface for which the flow separation is caused by the interaction between the viscous force, the pressure force, and the turbulence level. For flow over bluff bodies, a downstream recirculation bubble is created which increases the aspect ratio and/or the turbulence level results in flow reattachment close behind the obstacle.

Shieh, C. F.↗