Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Lowering the barrier to access information-rich transient kinetic data for machine learning methods

Transient kinetic data contain a wealth of information about intrinsic features of a catalyst as well as the reaction mechanism. Currently, high volume transient data is underutilized, and data science methods could both increase the value of information that can be extracted from this data, integrate experimental with theoretical data sources, and accelerate the pace of catalyst technology advancement. Transient kinetic characterizations with simple probe molecules exhibiting reversible adsorption, irreversible adsorption and bulk-surface diffusion are presented as training components for similar experiments with more complex surface reactions. In conclusion, by increasing the availability and accessibility of transient kinetic data through details of its structure and acquisition, we aim to decrease the barrier for data scientists to apply machine learning methods to this valuable data source.

Catalysis↗
DOE OSTI2025

Ensemble cure kinetics network (ECK-Net): A method to derive cure kinetics of thermosetting resin

This paper introduces an Ensemble Cure Kinetics Network (ECK-Net), a neural network (NN)–based framework for modeling the cure kinetics of thermosetting resins within a phenomenological context. ECK-Net replaces traditional analytic models, which require extensive chemical insight and multiple isothermal/non-isothermal experiments, with a data-driven surrogate that maps nonlinear relationships between temperature, degree of cure, and reaction rate from differential scanning calorimetry data. The proposed approach predicts input-dependent kinetic coefficients of a generalized nth-order reaction equation rather than reaction rates directly, enabling a single unified model to represent various epoxy systems without relying on iso-conversional analysis or predefined functional forms. To ensure robustness, multiple independently trained networks under different random initializations are blended through an ensemble strategy, effectively mitigating the stochastic variability inherent to neural networks. The framework is validated using experimental datasets from multiple resin systems, including aerospace-grade materials (Toray 3900-2, Cycom 5320-1, and Hexcel 8552) and a windmill-grade resin (RIMR 035c). The model accurately reproduces the temporal evolution of the degree of cure under manufacturers’ recommended cure cycles across all tested resins systems, yielding Pearson’s correlation coefficients of 0.992, 0.994, 0.993, 0.997, respectively. To demonstrate process-level applicability, the trained network was implemented within the Abaqus environment to simulate out-of-autoclave (OOA) curing process of the CFRP panel composed of Toray T830H-6K/3900-2D prepreg. The simulation results showed excellent agreement with experimental temperature response (maximum peak temperature, simulation: 189.6 °C, experiment: 188.5 °C) and the final degree of cure (simulation: 0.948, experiment: 0.960 ± 0.013), confirming ECK-Net’s capability as a reliable alternative to conventional cure kinetics modeling methods.

Composite curing↗
DOE OSTI2024

A hybrid Monte Carlo, discontinuous Galerkin method for linear kinetic transport equations

Here we present a hybrid method for time-dependent particle transport problems that combines Monte Carlo (MC) estimation with deterministic solutions based on discrete ordinates. For spatial discretizations, the MC algorithm computes a piecewise constant solution and the discrete ordinates use bilinear discontinuous finite elements. From the hybridization of the problem, the resulting problem solved by Monte Carlo is scattering free, resulting in a simple, efficient solution procedure. Between time steps, we use a projection approach to “relabel” collided particles as uncollided particles. In conclusion, from a series of standard 2-D Cartesian test problems we observe that our hybrid method has improved accuracy and reduction in computational complexity of approximately an order of magnitude relative to standard discrete ordinates solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Effect of natural carbon filler on thermo-oxidative degradation of thermoplastic-based composites

Here thermo-oxidative degradation of thermoplastic-based composites filled with two types of carbon filler including bituminous (Pittsburgh No.8, P8) and sub-bituminous (Powder River Basin, PRB) coals were investigated using differential scanning calorimetry (DSC) and thermogravimetric analysis. Oxidation induction time (OIT) and activation energy (AE) for coal plastic composites (CPCs), HDPE, and commercially available wood plastic composites (WPCs) were determined using isothermal and isoconversional kinetic methods. OIT values for CPC increased with coal content, indicating higher thermo-oxidative stability. AE values obtained using isothermal method were 52–116 kJ/mol for all CPC formulations, 32 kJ/mol for HDPE, and 33–114 kJ/mol for WPCs. AEs for CPCs obtained using isoconversional methods were not constant but rather dependent on the degree of degradation with CPC/ P8 coal possessing higher overall AE. Incorporating coal into HDPE increased thermo-oxidative stability of the resulting composite, indicating coal potentially acts as a natural primary and secondary antioxidant for polymer materials.

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

Kinetic theory based new upwind methods for inviscid compressible flows

Two new upwind methods called the Kinetic Numerical Method (KNM) and the Kinetic Flux Vector Splitting (KFVS) method for the solution of the Euler equations have been presented. Both of these methods can be regarded as some suitable moments of an upwind scheme for the solution of the Boltzmann equation provided the distribution function is Maxwellian. This moment-method strategy leads to a unification of the Riemann approach and the pseudo-particle approach used earlier in the development of upwind methods for the Euler equations. A very important aspect of the moment-method strategy is that the new upwind methods satisfy the entropy condition because of the Boltzmann H-Theorem and suggest a possible way of extending the Total Variation Diminishing (TVD) principle within the framework of the H-Theorem. The ability of these methods in obtaining accurate wiggle-free solution is demonstrated by applying them to two test problems.

Deshpande, S. M.↗
NASA NTRS1973

A compilation of studies from Atmospheric Variability Experiment (AVE)

Five methods for obtaining the sign of vertical motion at various levels in the atmosphere were investigated to determine which gave the best explanation for areas of rain and no-rain in the southeastern United States during the period February 19-22, 1964. The methods used were the terrain-induced vertical motion, the kinematic method including the terrain effect, the adiabatic method, the omega equation, and the vorticity equation combined with Dines' Compensation Principle. Stability and moisture availability were considered but not as extensively as vertical motion. Values of vertical motion obtained by the kinetic method, including orographic lifting produced the best agreement with the observed areas. When terms in the omega equation were added through use of constant multipliers, results comparable to the adiabatic method were obtained. Without this addition large uncertainties occurred. Maps showing areas where terms of the omega equation would indicate positive vertical motion did not correlate well with the observed rainfall patterns.

Scoggins, J. R.↗
NASA NTRS1981

Theory of strong turbulence by renormalization

The hydrodynamical equations of turbulent motions are inhomogeneous and nonlinear in their inertia and force terms and will generate a hierarchy. A kinetic method was developed to transform the hydrodynamic equations into a master equation governing the velocity distribution, as a function of the time, the position and the velocity as an independent variable. The master equation presents the advantage of being homogeneous and having fewer nonlinear terms and is therefore simpler for the investigation of closure. After the closure by means of a cascade scaling procedure, the kinetic equation is derived and possesses a memory which represents the nonMarkovian character of turbulence. The kinetic equation is transformed back to the hydrodynamical form to yield an energy balance in the cascade form. Normal and anomalous transports are analyzed. The theory is described for incompressible, compressible and plasma turbulence. Applications of the method to problems relating to sound generation and the propagation of light in a nonfrozen turbulence are considered.

Tchen, C. M.↗
NASA NTRS1972

Analysis of a boron-carbide-drum-controlled critical reactor experiment

In order to validate methods and cross sections used in the neutronic design of compact fast-spectrum reactors for generating electric power in space, an analysis of a boron-carbide-drum-controlled critical reactor was made. For this reactor the transport analysis gave generally satisfactory results. The calculated multiplication factor for the most detailed calculation was only 0.7-percent Delta k too high. Calculated reactivity worth of the control drums was $11.61 compared to measurements of $11.58 by the inverse kinetics methods and $11.98 by the inverse counting method. Calculated radial and axial power distributions were in good agreement with experiment.

Mayo, W. T.↗
NASA NTRS1972

Approximation method for the kinetic Boltzmann equation

The further development of a method for approximating the Boltzmann equation is considered and a case of pseudo-Maxwellian molecules is treated in detail. A method of approximating the collision frequency is discussed along with a method for approximating the moments of the Boltzmann collision integral. Since the return collisions integral and the collision frequency are expressed through the distribution function moments, use of the proposed methods make it possible to reduce the Boltzmann equation to a series of approximating equations.

Shakhov, Y. M.↗
NASA NTRS2011

Wind and Temperature Spectrometry of the Upper Atmosphere in Low-Earth Orbit

Wind and Temperature Spectrometry (WATS) is a new approach to measure the full wind vector, temperature, and relative densities of major neutral species in the Earth's thermosphere. The method uses an energy-angle spectrometer moving through the tenuous upper atmosphere to measure directly the angular and energy distributions of the air stream that enters the spectrometer. The angular distribution gives the direction of the total velocity of the air entering the spectrometer, and the energy distribution gives the magnitude of the total velocity. The wind velocity vector is uniquely determined since the measured total velocity depends on the wind vector and the orbiting velocity vector. The orbiting spectrometer moves supersonically, Mach 8 or greater, through the air and must point within a few degrees of its orbital velocity vector (the ram direction). Pointing knowledge is critical; for example, pointing errors 0.1 lead to errors of about 10 m/s in the wind. The WATS method may also be applied without modification to measure the ion-drift vector, ion temperature, and relative ion densities of major ionic species in the ionosphere. In such an application it may be called IDTS: Ion-Drift Temperature Spectrometry. A spectrometer-based coordinate system with one axis instantaneously pointing along the ram direction makes it possible to transform the Maxwellian velocity distribution of the air molecules to a Maxwellian energy-angle distribution for the molecular flux entering the spectrometer. This implementation of WATS is called the gas kinetic method (GKM) because it is applied to the case of the Maxwellian distribution. The WATS method follows from the recognition that in a supersonic platform moving at 8,000 m/s, the measurement of small wind velocities in the air on the order of a few 100 m/s and less requires precise knowledge of the angle of incidence of the neutral atoms and molecules. The same is true for the case of ion-drift measurements. WATS also provides a general approach that can obtain non-equilibrium distributions as may exist in the upper regions of the thermosphere, above 500 km and into the exosphere. Finally, WATS serves as a mass spectrometer, with very low mass resolution of roughly 1 part in 3, but easily separating atomic oxygen from molecular nitrogen.

Herrero, Federico↗
DOE OSTI2023

Compressible Rayleigh–Taylor instability subject to isochoric initial background state

The effect of compressibility on the single-mode Rayleigh–Taylor instability is examined using two (2D) and three-dimensional (3D) direct numerical simulations. To isolate compressibility from background stratification effects, this work employs a constant density profile on each side of the interface. The numerical simulations are performed at various Reynolds numbers using the gas kinetic method for static Mach numbers up to M = 0.4. The most important finding is that compressibility acting in isolation enhances the instability and perturbations grows faster with increasing Mach number, unlike previous results with background isothermal state, which show suppression of the instability at higher static Mach numbers. In addition, compressibility is also shown to increase the bubble-spike asymmetry. While the instability grows faster for the 3D case, the findings are qualitatively similar in 2D and 3D. Furthermore, the dynamical reasons underlying the effect of compressibility are elucidated by examining the evolution of vorticity and turbulent kinetic energy transport equations.

42 ENGINEERING↗
DOE OSTI2021

Non-Boltzmann Effects in Chain Branching and Pathway Branching for Diethyl Ether Oxidation

Low-temperature (LT) engine applications have several potential benefits, including reduced emissions and increased efficiency. Attaining these benefits requires accurate kinetic modeling of LT chain branching, which depends heavily on ketohydroperoxide (KHP) decomposition. For diethyl ether (DEE), a promising biofuel, current estimates of the KHP decomposition rate constant are largely based on empirical fits to data. In this study, we investigate the most important KHP isomer in DEE LT oxidation by applying variable reaction coordinate transition state theory to the main pathway for KHP decomposition: OO bond fission to produce •OH and a keto-alkoxy radical, •OQ'O. We also use ab initio kinetics methods to investigate the decomposition of •OQ'O, where we find dominant branching to acetic acid, with the remaining flux going to CH 3 C(O)OCHO. Additionally, new time-resolved measurements of DEE and acetic acid concentrations during LT (450–600 K) DEE oxidation are obtained in a laser photolysis flow reactor coupled with multiplexed photoionization mass spectrometry. These new experimental data, along with jet-stirred reactor data in the literature, are compared with the predictions of a recent DEE mechanism (Tran et al. Proc. Comb. Inst. 2019, 37, 511-519) that was modified with the newly calculated ab initio rate constants for KHP and •OQ'O decomposition. The predictions of the modified mechanism are quite poor when compared to the experimental data; this is primarily due to the new KHP ⇌ •OQ'O + •OH rate constant, which is 1–2 orders of magnitude slower than empirical values employed in recent mechanisms. To reconcile the new KHP rate constant and the experimental data, we explore and quantify the possible role of non-Boltzmann (nB) reaction sequences. The nB reactions have a substantial effect on both the overall mechanism reactivity and the •OQ'O branching to acetic acid. We also provide guidance on the proper implementation of nB reactions in kinetic mechanisms.

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

Zinc Speciation and Propane Dehydrogenation in Zn/H-ZSM-5 Catalysts

Zn/H-ZSM-5 catalysts have been frequently investigated for propane dehydrogenation (PDH); however, the active site remains unresolved due to the complexity of the system. We employed in situ FTIR spectroscopy and a kinetics method to correlate the Zn speciation and PDH activity in Zn/H-ZSM-5 with two Si/Al ratios (15 and 39) and a range of Zn/Al ratios (0–1.7). Incremental additions of zinc show that Zn 2+ sites are preferentially formed on H-ZSM-5 over a fraction of paired Al sites followed by [Zn-O-Zn] 2+ and [ZnOH] + sites and then ZnO x clusters. The [Zn-OH] + and [Zn-O-Zn] 2+ sites in H-ZSM-5 are more active and selective than isolated Zn 2+ for PDH. [Zn-OH] + species sublimate over time on stream, leading to catalyst deactivation, while [Zn-O-Zn] 2+ species are stable even after high-temperature reduction (750 °C for 60 min). Three distinct Zn sites ([Zn-O-Zn] 2+ , Zn 2+ , and [ZnOH] + ) show a similar propane reaction order (close to 1) and H 2 reaction order (close to 0). Combined with the lack of Zn hydride when propane flows over the catalyst at 550 °C, it is concluded that propane adsorption and dissociation is a rate-determining step and H 2 desorption is fast. This work indicates that the preparation of H-ZSM-5 with abundant Al pairs may be a strategy to form stable and selective Zn/H-ZSM-5 catalysts for propane dehydrogenation. Furthermore, it is also highlighted that examining the effect of both metal/Al ratios and Al distribution of the zeolite is crucial in identifying the metal cations in metal–zeolite systems.

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

Proton affinities of pertechnetate (TcO 4 - ) and perrhenate (ReO 4 - )

The anions pertechnetate, TcO 4 - , and perrhenate, ReO 4 - , exhibit very similar chemical and physical properties. Revealing and understanding disparities between them enhances fundamental understanding of both. Electrospray ionization generated the gas-phase proton bound dimer (TcO 4 - )(H + )(ReO 4 - ). Collision induced dissociation of the dimer yielded predominantly HTcO 4 - and ReO 4 - , which according to Cooks’ kinetic method indicates that the proton affinity (PA) of TcO 4 - is greater than that of Re O 4 - . Density functional theory computations agree with the experimental observation, providing PA[TcO 4 - ] = 300.1 kcal mol -1 and PA[ReO 4 - ] = 297.2 kcal mol -1 . Attempts to rationalize these relative PAs based on elementary molecular parameters such as atomic charges indicate that the entirety of bond formation and concomitant bond disruption needs to be considered to understand the energies associated with such protonation processes. Although in both the gas and solution phases, TcO 4 - is a stronger base than ReO 4 - , it is noted that the significance of even such qualitative accordance is tempered by the very different natures of the underlying phenomena.

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

Ten-moment fluid model for low-temperature magnetized plasmas

In this paper, a one-dimensional 10-moment multi-fluid plasma model is developed and applied to low-temperature magnetized plasmas. The 10-moment model solves for six anisotropic pressure terms, in addition to density and three components of fluid momentum, which allows the model to capture finite kinetic effects. The results are benchmarked against a 5-moment model, which assumes that the gas constituents follow a Maxwellian velocity distribution function (VDF), and a particle-in-cell Monte Carlo collision model that allows for arbitrary non-Maxwellian VDFs. The models are compared in a one-dimensional, low-temperature, partially magnetized plasma test case. The 10-moment results accurately reproduce the anisotropic temperature profile in low-temperature magnetized plasmas, where shear gradients exist due to the E×B drift. We discuss the mechanisms by which the anisotropic pressure can be generated in low-temperature magnetized plasmas. In addition, the importance of a self-consistent heat flux closure to the 10-moment model is studied, showing consistency with other models only when the assumptions of the underlying model are met. The 10-moment model allows for study of electron inertia effects and non-Maxwellian VDFs without the need for kinetic methods that are more computationally expensive.

Kuldinow, Derek Amur (ORCID:0000000319730196)↗
DOE OSTI2022

Neural Phase Simulation

The Neural Phase Simulation (NPS) is a package of codes for simulating microstructure evolution and accelerated molecular dynamics with deep neural-networks based surrogate models. NPS is designed to offer quantitatively accurate and computationally efficient simulation capabilities by leveraging modern machine-learning techniques. The primary intended use cases of NPS are training neural network surrogate models, though performing simulations on a single node is also supported. The NPS surrogate models can be trained from ground truth simulation methods, which are supposed to be accurate but expensive, such as molecular dynamics, phase field methods, kinetic Monte Carlo and discrete dislocation dynamics.

Zhou, Fei↗