Search NASASearch

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 253 records · Page 14

A Study of the Kinetics of the CO Oxidation Catalyst in a Human Spaceflight Fire Cartridges as a Method to Understand and Predict Performance

Fire Cartridges (FC) are an integral part of fire response on ISS, Orion and Gateway. When a pair of FCs are integrated with an Emergency Mask, they can be used to provide safe breathing air on these vehicles in a post-fire environment. The FCs have two primary mechanisms for removing contaminants from the air, an activated carbon (AC) bed which adsorbs the majority of fire byproducts, and a catalyst bed consisting of gold nanoparticles supported on iron oxide which oxidized CO to CO2 under near ambient conditions. While the FC catalyst has proven its ability to convert CO under relevant conditions, this paper attempts to address the apparent inconsistent performance of the catalyst seen in testing. Numerous variables effect the performance of the catalyst, including but not limited to: reaction environment, reactant gas composition, inert gas composition, catalyst age, and catalyst lot. This paper discusses the importance of these variables and attempts a fundamental mechanistic explanation for the effects of each variable. The fundamental mechanisms are then used to establish the magnitude of expected performance variability, which were validated against real test data of flight quality FCs.

Fire Cartridge

A Study of the Kinetics of the CO Oxidation Catalyst in a Human Spaceflight Fire Cartridges as a Method to Understand and Predict Performance

Fire Cartridges (FC) are an integral part of fire response on ISS, Orion and Gateway. When a pair of FCs are integrated with an Emergency Mask, they can be used to provide safe breathing air on these vehicles in a post-fire environment. The FCs have two primary mechanisms for removing contaminants from the air, an activated carbon (AC) bed which adsorbs the majority of fire byproducts, and a catalyst bed consisting of gold nanoparticles supported on iron oxide which oxidized CO to CO2 under near ambient conditions. While the FC catalyst has proven its ability to convert CO under relevant conditions, this paper attempts to address the apparent inconsistent performance of the catalyst seen in testing. Numerous variables effect the performance of the catalyst, including but not limited to: reaction environment, reactant gas composition, inert gas composition, catalyst age, and catalyst lot. This paper discusses the importance of these variables and attempts a fundamental mechanistic explanation for the effects of each variable. The fundamental mechanisms are then used to establish the magnitude of expected performance variability, which were validated against real test data of flight quality FCs.

Fire Cartridge

Improved Modeling of Finite-Rate Turbulent Combustion Processes in Research Combustors

The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.

VanOverbeke, Thomas J.

Isothermal Analysis of the Crystallization Kinetics in Lithium Disilicate Glass using Trans Temp Furnace

Crystallization kinetics for lithium disilicate, Li2O2SiO2, (LS2) glass has been studied extensively by nonisothermal methods, but only a few studies on the isothermal crystallization kinetics of LS2 are available. In the present research, isothermal crystallization experiments or the LS2 glass were conducted in a Trans Temp furnace between 600 and 635 C, and selected properties such as the activation energy for crystallization (E), crystal growth index or Avrami parameter (n), the concentration of quenched-in nuclei in the starting glass (Ni) and the crystal nucleation rate (I) were measured. The crystal nucleation rate (I) was measured at only one selected temperature of 452 C, at this time. This commercial furnace has a 13 cm long isothermal heating zone (+/- 1 C) that allows precise heat treatment of relatively large samples. By placing a thermocouple within approx. 2 mm of the sample, it was possible to detect the heat of crystallization in the form of an isothermal crystallization exotherm during isothermal heat treatment of the sample. The values of E (318 plus or minus 10 kJ/mol), n (3.6 plus or minus 0.l), and N(sub i) (1.6 x 10(exp l2) m(sup -3)) calculated by analyzing these isotherms using the standard Johnson-Mehl-Avrami (JMA) equation were reproducible and in agreement with the literature values. The value of I, 1.9 x 10(exp 10) m(sup -3) s(sup -1) at 452 C, is an order of magnitude higher than the reported value for LS2.

Fuss, T.

Exponential-fitted methods for integrating stiff systems of ordinary differential equations: Applications to homogeneous gas-phase chemical kinetics

Conventional algorithms for the numerical integration of ordinary differential equations (ODEs) are based on the use of polynomial functions as interpolants. However, the exact solutions of stiff ODEs behave like decaying exponential functions, which are poorly approximated by polynomials. An obvious choice of interpolant are the exponential functions themselves, or their low-order diagonal Pade (rational function) approximants. A number of explicit, A-stable, integration algorithms were derived from the use of a three-parameter exponential function as interpolant, and their relationship to low-order, polynomial-based and rational-function-based implicit and explicit methods were shown by examining their low-order diagonal Pade approximants. A robust implicit formula was derived by exponential fitting the trapezoidal rule. Application of these algorithms to integration of the ODEs governing homogenous, gas-phase chemical kinetics was demonstrated in a developmental code CREK1D, which compares favorably with the Gear-Hindmarsh code LSODE in spite of the use of a primitive stepsize control strategy.

Pratt, D. T.

Numerical schemes for 3-wave kinetic equations: A complete treatment of the collision operator

In our previous work Walton and Tran (2023), numerical schemes for a simplified version of 3-wave kinetic equations, in which only the simple forward-cascade terms of the collision operators are kept, have been successfully designed, especially to capture the long time dynamics of the equation given the multiple blow-up time phenomenon. In this second work in the series, we propose numerical treatments for the complete 3-wave kinetic equations, in which the complete, much more complicated collision operators are fully considered based on a novel conservative form of the equation. Here we then derive an implicit finite volume scheme to solve the equation. The new discretization uses an adaptive time-stepping method which allows for the simulations to be carried to very long times. Our computed solutions are compared with previously derived long-time asymptotic estimates for the decay rate of total energy of time-dependent solutions of 3-wave kinetic equations and found to be in excellent agreement.

97 MATHEMATICS AND COMPUTING

Atomic oxygen testing with thermal atom systems - A critical evaluation

The use of thermal atom (kinetic energy near 0.04 eV) test methods as a materials selection and screening technique for LEO spacecraft is critically evaluated in this paper. The physics and chemistry of the thermal atom environments are shown to produce specific mass loss rates (mg/sq cm per min) and reaction efficiencies (Re) radically different from those produced in the LEO environment. A response surface study shows that specific mass loss rates change rapidly with plasma-asher parameters and seldom agree with flight data. FEP Teflon is shown to react by a different mechanism than Kapton, polyethylene, or graphite. The Re (Re = volume of material removed/oxygen atom) of Kapton, polyethylene, Mylar, Tedlar, FEP Teflon, and graphite measured in a flowing afterglow apparatus are 0.001 to 0.0001 those measured with high-energy atoms (kinetic energy 1.5 eV or greater) in beam systems or in LEO. The effect of sample temperature and atom impact energy on Re is discussed. A simple kinetic model describing the reaction of atomic oxygen with polymer surfaces is developed. Guidelines and recommendations for thermal atom testing and interpretation of test results are presented.

Koontz, Steven L.

A Bottom-Up Approach to Rational Design of Crystalline Materials: Investigation of Vibronic Coherences Underlying Exciton Dynamics in Semiconductors

In this project we uncovered structure-function relationships of donor-acceptor co-crystals used to develop next-generation optoelectronic devices. Unraveling the photodynamics of molecular crystalline materials poses many challenges for spectroscopy due to broad, overlapping features representing numerous underlying dynamical processes. This leads researchers to make many assumptions about the dynamics of a system in choosing an appropriate kinetic fitting model. Computationally, electronic structure methods are either prohibitively expensive or underdeveloped for computing the excited state structure of molecular materials, especially states that exhibit charge transfer. Researchers must therefore perform calculations of excited electronic states using truncated models of molecular materials. Here we present a joint experimental-theoretical approach to bridging the gap between the photodynamics of a molecular material and its constituent molecules. We focus our efforts on quantifying the timescales and mechanisms of photoexcitation in donor-acceptor co-crystals and donor-acceptor dimers where the lowest-lying excited state is characterized by charge transfer from the donor to the acceptor. We employ ultrafast UV pump, UV-Vis probe transient absorption spectroscopy to unravel the time-resolved spectroscopic signatures of the photodynamics in both the crystalline material and donor-acceptor dimers in solution. We perform electronic structure and excited state dynamics calculations of the dimers to inform kinetic fitting models and assign the spectral features. The photodynamics of the crystal vs. dimer systems have many similarities, enabling unprecedented insights into the formation and evolution of charge transfer excitons in the crystalline systems.

36 MATERIALS SCIENCE

Studies in electron phenomena in MOS structures: The pulsed C-V method

The pulse hysteresis capacitance voltage (C-V) provides a straight forward technique for measuring the change of various charges in MOS structures and a tool for investigating the kinetics of various electron phenomena is developed and described. The method can be used for measuring the energy distribution and kinetics of surface states with the resolution of about 1/5 x 10 to the -9 power cm eV. Some transients in an MOS structure, particularly, the thermal generation of minority charge carriers via surface states and the relaxation of minority charge carriers supplied from the inversion layer outside the MOS structure are theoretically investigated. Analytical expressions which clearly present the physics of those electron phenomena are derived.

Kaplan, G.

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics

Closed Loop Terminal Guidance Navigation for a Kinetic Impactor Spacecraft

A kinetic impactor spacecraft is a viable method to deflect an asteroid which poses a threat to the Earth. The technology to perform such a deflection has been demonstrated by the Deep Impact (DI) mission, which successfully collided with comet Tempel 1 in July 2005 using an onboard autonomous navigation system, called AutoNav, for the terminal phase of the mission. In this paper, we evaluate the ability of AutoNav to impact a wide range of scenarios that an deflection mission could encounter, varying parameters such as the approach velocity, phase angle, size of the asteroid, and the determination of spacecraft attitude. Using realistic Monte Carlo simulations, we tabulated the probability of success of the deflection as a function of these parameters, and the highest sensitivity to be due the spacecraft attitude determination mode. In addition, we also specifically analyzed the impact probability for a proposed mission which would send an impactor to the asteroid 1999RQ36. We conclude with some recommendations for future work.

closed loop guidance

Curing Composite Materials Using Lower-Energy Electron Beams

In an improved method of fabricating composite-material structures by laying up prepreg tapes (tapes of fiber reinforcement impregnated by uncured matrix materials) and then curing them, one cures the layups by use of beams of electrons having kinetic energies in the range of 200 to 300 keV. In contrast, in a prior method, one used electron beams characterized by kinetic energies up to 20 MeV. The improved method was first suggested by an Italian group in 1993, but had not been demonstrated until recently. With respect to both the prior method and the present improved method, the impetus for the use of electron- beam curing is a desire to avoid the high costs of autoclaves large enough to effect thermal curing of large composite-material structures. Unfortunately, in the prior method, the advantages of electron-beam curing are offset by the need for special walls and ceilings on curing chambers to shield personnel from x rays generated by impacts of energetic electrons. These shields must be thick [typically 2 to 3 ft (about 0.6 to 0.9 m) if made of concrete] and are therefore expensive. They also make it difficult to bring large structures into and out of the curing chambers. Currently, all major companies that fabricate composite-material spacecraft and aircraft structures form their layups by use of automated tape placement (ATP) machines. In the present improved method, an electron-beam gun is attached to an ATP head and used to irradiate the tape as it is pressed onto the workpiece. The electron kinetic energy between 200 and 300 keV is sufficient for penetration of the ply being laid plus one or two of the plies underneath it. Provided that the electron-beam gun is properly positioned, it is possible to administer the required electron dose and, at the same time, to protect personnel with less shielding than is needed in the prior method. Adequate shielding can be provided by concrete walls 6 ft (approximately equal to 1.8 m) high and 16 in. (approximately equal to 41 cm) thick, without a ceiling. The success of the present method depends on the use of a cationic epoxy as the matrix material in the prepreg tape, heating the prepreg tape to a temperature of 50 C immediately prior to layup, and exposing the workpiece to an electron-beam dose of approximately 2 Mrad. Experiments have shown that structures fabricated by the present method have the same mechanical properties as those of nominally identical structures fabricated by the prior method with electron beams of 3 to 4 MeV.

Byrne, Catherine A.

Charge Transport in Solvated Donor–Acceptor Functionalized Peptoids: Molecular Dynamics and Rate Theory

Scalable solar-energy conversion requires photoactive materials that combine the efficiency of natural photosynthetic systems with the stability and processability needed for practical applications. Achieving reliable charge transport in soft, self-assembled organic materials remains challenging, as structural fluctuations and environmental effects strongly influence charge-transfer (CT) rates. Here, we present a broadly applicable computational framework for evaluating CT rates in the condensed phase, combining Fermi’s golden rule rate theory with inputs from all-atom molecular dynamics (MD) simulations and first-principles electronic-structure calculations. The approach does not rely on system-specific parametrization and is applicable to a wide range of soft and disordered materials. We demonstrate the applicability and usefulness of the framework on redox-active peptoids functionalized with iron–porphyrin (Fe–P) complexes, a bioinspired platform with programmable donor–acceptor units and tunable three-dimensional organization. The calculated CT rates exhibit strong sensitivity to molecular conformation, with variations spanning several orders of magnitude. This dependence is shown to arise from the pronounced variation in diabatic electronic coupling with the relative orientations and separations of the Fe–P complexes across the conformational ensemble. The framework provides a consistent route for connecting atomistic structure to CT kinetics in the condensed phase and enables analysis of structure–rate relations in organic semiconducting systems.

Charge transfer

Gyrokinetic electromagnetic particle simulations in triangular meshes with C1 finite elements

The triangular mesh-based gyrokinetic scheme enables comprehensive axis-to-edge studies across the entire plasma volume. Our approach employs triangular finite elements with first-derivative continuity (C1), building on previous work to facilitate gyrokinetic simulations. Additionally, we have adopted the mixed variable/pullback scheme for gyrokinetic electromagnetic particle simulations. The filter-free treatment in the poloidal cross-section with triangular meshes introduces unique features and challenges compared to previous treatments using structured meshes. Our implementation has been validated through benchmarks using ITPA-toroidicity-induced Alfvén eigenmode parameters, showing its capability in moderate to small electron skin depth regimes. Additional examinations using experimental parameters confirm its applicability to realistic plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Closure and transport theory for plasmas with multiple ion species

This report summarizes a DOE-funded research project on closure and transport theory for plasmas with multiple ion species. The work developed accurate moment-based closure models across the full range of collisionality, incorporating both kinetic and collisional effects. Key contributions include high-fidelity closure relations, multi-temperature models, and transport coefficients applicable to fusion, space, and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Multimodal Ligand Non-Innocence for Reduction Catalysis (Final Technical Report)

The goal of this project was to develop multimodal ligand non-innocence in late transition metal complexes of aminophenolate ligands and to apply this ligand-centered reactivity to small-molecule reduction reactions relevant to energy-conversion chemistry. Specific aims were, 1) Prepare new coordination complexes of aminophenolate ligands and characterize them using a battery of structural, spectroscopic, and computational methods. 2) Benchmark the thermodynamics and kinetics of e - , H + , H • , and H - transfer reactivity to establish the baseline reactivity of new metal complexes containing aminophenolate ligands. 3) Apply new aminophenolate complexes to catalytic reduction reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Production and analysis of thermal decomposition products from polymeric materials

A description is presented of a strategy for analyzing the combustion process and the degradation products which are formed. One of three primary objectives in the study of polymer degradation is related to the characterization of the material to be studied and the investigation of the thermal behavior of the material. Another objective is concerned with the definition of the nature of the decomposition process by identification and quantitation of the degradation products. The third objective involves the determination of the mechanism and kinetics of the decomposition process. The methods of sample degradation include pyrolysis, oxidative degradation, flaming combustion, and the use of large-scale combustion chambers. Methods of chemical separation and identification are considered, taking into account low-boiling volatiles, high-boiling volatiles, and ancillary techniques.

Chatfield, D. A.