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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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.↗
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↗
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.↗
NASA NTRS1984

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

A comparison of the efficiency of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations is presented. The methods examined include two general-purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D, and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes are applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code currently available for the integration of combustion kinetic rate equations. An important finding is that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient than evaluating the temperature by integrating its time-derivative.

Radhakrishnan, K.↗
NASA NTRS1994

ASTM E 1559 method for measuring material outgassing/deposition kinetics has applications to aerospace, electronics, and semiconductor industries

The American Society for Testing and Materials has published a new standard test method for characterizing time and temperature-dependence of material outgassing kinetics and the deposition kinetics of outgassed species on surfaces at various temperatures. This new ASTM standard, E 1559(1), uses the quartz crystal microbalance (QCM) collection measurement approach. The test method was originally developed under a program sponsored by the United States Air Force Materials Laboratory (AFML) to create a standard test method for obtaining outgassing and deposition kinetics data for spacecraft materials. Standardization by ASTM recognizes that the method has applications beyond aerospace. In particular, the method will provide data of use to the electronics, semiconductor, and high vacuum industries. In ASTM E 1559 the material sample is held in vacuum in a temperature-controlled effusion cell, while its outgassing flux impinges on several QCM's which view the orifice of the effusion cell. Sample isothermal total mass loss (TML) is measured as a function of time from the mass collected on one of the QCM's which is cooled by liquid nitrogen, and the view factor from this QCM to the cell. The amount of outgassed volatile condensable material (VCM) on surfaces at higher temperatures is measured as a function of time during the isothermal outgassing test by controlling the temperatures of the remaining QCM's to selected values. The VCM on surfaces at temperatures in between those of the collector QCM's is determined at the end of the isothermal test by heating the QCM's at a controlled rate and measuring the mass loss from the end of the QCM's as a function of time and temperature. This reevaporation of the deposit collected on the QCM's is referred to as QCM thermogravimetric analysis. Isothermal outgassing and deposition rates can be determined by differentiating the isothermal TML and VCM data, respectively, while the evaporation rates of the species can be obtained as a function of temperature by differentiating the QCM thermogravimetric analysis data.

Garrett, J. W.↗
NASA NTRS1993

An upwind, kinetic flux-vector splitting method for flows in chemical and thermal non-equilibrium

We have developed new upwind kinetic difference schemes for flows with non-equilibrium thermodynamics and chemistry. These schemes are derived from the Boltzmann equation with the resulting Euler schemes developed as moments of the discretized Boltzmann scheme with a locally Maxwellian velocity distribution. Splitting the velocity distribution at the Boltzmann level is seen to result in a flux-split Euler scheme and is called Kinetic Flux Vector Splitting (KFVS). Extensions to flows with finite-rate chemistry and vibrational relaxation is accomplished utilizing nonequilibrium kinetic theory. Computational examples are presented comparing KFVS with the schemes of Van Leer and Roe for a quasi-one-dimensional flow through a supersonic diffuser, inviscid flow through two-dimensional inlet, and viscous flow over a cone at zero angle-of-attack. Calculations are also shown for the transonic flow over a bump in a channel and the transonic flow over an NACA 0012 airfoil. The results show that even though the KFVS scheme is a Riemann solver at the kinetic level, its behavior at the Euler level is more similar to the existing flux-vector splitting algorithms than to the flux-difference splitting scheme of Roe.

Eppard, W. M.↗
NASA NTRS2004

Effects of Kinetic Processes in Shaping Io's Global Plasma Environment: A 3D Hybrid Model

The global dynamics of the ionized and neutral components in the environment of Io plays an important role in the interaction of Jupiter's corotating magnetospheric plasma with Io. The stationary simulation of this problem was done in the MHD and the electrodynamics approaches. One of the main significant results from the simplified two-fluid model simulations was a production of the structure of the double-peak in the magnetic field signature of the I0 flyby that could not be explained by standard MHD models. In this paper, we develop a method of kinetic ion simulation. This method employs the fluid description for electrons and neutrals whereas for ions multilevel, drift-kinetic and particle, approaches are used. We also take into account charge-exchange and photoionization processes. Our model provides much more accurate description for ion dynamics and allows us to take into account the realistic anisotropic ion distribution that cannot be done in fluid simulations. The first results of such simulation of the dynamics of ions in the Io's environment are discussed in this paper.

Lipatov, Alexander S.↗
DOE OSTI2025

Structure-preserving neural networks for the regularized entropy-based closure of a linear, kinetic, radiative transport equation

The main challenge of large-scale numerical simulation of radiation transport is the high memory and computation time requirements of discretization methods for kinetic equations. In this work, we derive and investigate a neural network-based approximation to the entropy-based closure method to accurately compute the solution of the multi-dimensional moment system with a low memory footprint and competitive computational time. We extend methods developed for the standard entropy-based closure to the regularized entropy-based closures. The main idea is to interpret structure-preserving neural network approximations of the regularized entropy-based closure as a two-stage approximation to the original entropy-based closure. We conduct a numerical analysis of this approximation and investigate optimal parameter choices. Our numerical experiments demonstrate that the method has a much lower memory footprint than traditional methods with competitive computation times and simulation accuracy. The code and all trained networks are provided on GitHub.

entropy closure↗
DOE OSTI2026

Catalytic resonance theory for the kinetics of photon-promoted catalysis

The illumination of catalytic surfaces with a continuous or pulsed stream of photons dynamically modulates surface chemistry for faster rates, non-equilibrium conversion, or product selectivity control. To establish fundamental principles of dynamic photon-modulated catalysis, the photocatalytic conversion of a generic surface reaction was simulated using the kinetic Monte Carlo method to understand the kinetic implications of an independent stream of photons that promotes surface product desorption. The time-averaged photocatalytic rate at differential conditions for varying photon flux and temperatures indicated three kinetic regimes described by product thermal desorption control, surface reaction control, and an intermediate kinetic regime with a zero slope Arrhenius plot, consistent with a degree of rate control dominated by the photon arrival frequency (i.e., per-site photon flux). Here, the maximum photocatalytic rate occurred orders of magnitude above the Sabatier limit at the resonance frequency, identified as the photon arrival frequency matching the surface reaction rate constant.

Canavan, Jesse R. [University of Minnesota, Minnea↗
NASA NTRS1988

Simplified half-life methods for the analysis of kinetic data

The analysis of reaction rate data has as its goal the determination of the order rate constant which characterize the data. Chemical reactions with one reactant and present simplified methods for accomplishing this goal are considered. The approaches presented involve the use of half lives or other fractional lives. These methods are particularly useful for the more elementary discussions of kinetics found in general and physical chemistry courses.

Eberhart, J. G.↗
DOE OSTI2025

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Identification of local modes using modal mass distributions

An approach for local modal identification using a modal mass spatial distribution matrix is proposed. This modal mass matrix is obtained by overlay multiplying the modal matrix to the product of the mass and the modal matrices of the structure. Each column in the spatial modal mass matrix can be interpreted as a measure of the kinetic energy distribution over the physical degrees of freedom for the corresponding mode. Local modes are identified due to the highly concentrated kinetic energy distribution in the localized regions. This spatial modal mass approach is more sensitive than the modal amplitude ratio method, since kinetic energy is in the quadratic form of modal amplitudes. Numerical examples involving identification of local modes using the proposed modal mass spatial distribution matrix approach are presented and discussed.

Chou, Chaur-Ming↗
NASA NTRS1995

Grain formation around carbon stars. 1: Stationary outflow models

Asymptotic giant branch (AGB) stars are known to be sites of dust formation and undergo significant mass loss. The outflow is believed to be driven by radiation pressure on grains and momentum coupling between the grains and gas. While the physics of shell dynamics and grain formation are closely coupled, most previous models of circumstellar shells have treated the problem separately. Studies of shell dynamics typically assume the existence of grains needed to drive the outflow, while most grain formation models assume a constant veolcity wind in which grains form. Furthermore, models of grain formation have relied primarily on classical nucleation theory instead of using a more realistic approach based on chemical kinetics. To model grain formation in carbon-rich AGB stars, we have coupled the kinetic equations governing small cluster growth to moment equations which determine the growth of large particles. Phenomenological models assuming stationary outflow are presented to demonstrate the differences between the classical nucleation approach and the kinetic equation method. It is found that classical nucleation theory predicts nucleation at a lower supersaturation ratio than is predicted by the kinetic equations, resulting in significant differences in grain properties. Coagulation of clusters larger than monomers is unimportant for grain formation in high mass-loss models but becomes more important to grain growth in low mass-loss situations. The properties of the dust grains are altered considerably if differential drift velocities are ignored in modeling grain formation. The effect of stellar temperature, stellar luminosity, and different outflow velocities are investigated. The models indicate that changing the stellar temperature while keeping the stellar luminosity constant has little effect on the physical parameters of the dust shell formed. Increasing the stellar luminosity while keeping the stellar temperature constant results in large differences in grain properties. For small outflow velocities, grains form at lower supersaturation ratios and close to the stellar photosphere, resulting in larger but fewer grains. The reverse is true when grains form under high outflow velocities, i.e., they form at higher supersaturation ratios, farther from the star, and are much smaller but at larger quantities.

Egan, Michael P.↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗