SEARCH · Search NASA

Results for “Kinetic model”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS1994

A three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawn and dusk sides of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H+ fluxes at tens of keV, which are always over-estimated. A newly-invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm-time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, j(sub o)(1+Ay(exp n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗

NASA NTRS1995

Three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawnside and duskside of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H(+) fluxes at tens of keV, which are always overestimated. A newly invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, J(sub o)(1 + Ay(sup n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗

DOE OSTI2025

Transient kinetic insights into selective propene oxidation over industrial bismuth molybdate catalysts

Selective oxidation of propene to acrolein over industrial multicomponent bismuth molybdate (BMO) catalysts significantly depends on reaction conditions that include operating parameters and catalyst state. Here, this work investigates selective propene oxidation in the intrinsic kinetic paradigm of Temporal Analysis of Products (TAP) reactor by systematically varying catalyst redox state, temperature, and oxygen-to-propene feed ratio. A 93 % propene conversion with an acrolein yield of 80 % is achieved at elevated temperatures (450 °C) on oxidized catalysts under oxygen-rich conditions (O 2 :C 3 H 6 = 10). However, these conditions diminish acrolein-to-CO 2 selectivity due to enhanced total oxidation to CO 2 . In contrast, a reduced catalyst state, moderate temperature (350 °C), and lower oxygen feed (O 2 :C 3 H 6 = 1) nearly doubles the acrolein to-CO 2 selectivity, albeit at a lower acrolein yield (33 %). Transient kinetic studies together with a kinetic model that simplify the major product formation pathways in lumped non-elementary forms reveal that the availability of surface oxygen species plays a pivotal role in governing reaction pathways. Additionally, density functional theory (DFT) calculations on pure BMO catalysts inform the role of surface redox states on propene and oxygen activation barriers. Lattice oxygen at acrolein-selective sites drives both acrolein and CO 2 formation, while adsorbed oxygen at activation sites favors unselective CO 2 generation. This work establishes a critical relationship between transient product selectivity and surface oxygen availability, which is strongly influenced by catalyst redox state, feed ratio, and reaction temperature. These insights underscore the importance of dynamic reactor operation strategies and offer a foundation for designing next-generation propene oxidation processes with tunable acrolein selectivity.

42 - ENGINEERING↗

NASA NTRS1999

A Gas-Kinetic Scheme for Multimaterial Flows and Its Application in Chemical Reaction

This paper concerns the extension of the multicomponent gas-kinetic BGK-type scheme to multidimensional chemical reactive flow calculations. In the kinetic model, each component satisfies its individual gas-kinetic BGK equation and the equilibrium states of both components are coupled in space and time due to the momentum and energy exchange in the course of particle collisions. At the same time, according to the chemical reaction rule one component can be changed into another component with the release of energy, where the reactant and product could have different gamma. Many numerical test cases are included in this paper, which show the robustness and accuracy of kinetic approach in the description of multicomponent reactive flows.

Lian, Yongsheng↗

DOE OSTI2024

SHOCK TUBE IGNITION STUDIES OF RENEWABLE DIESEL FUELS FOR MEDIUM AND HEAVY-DUTY TRANSPORTATION

Currently extensive research on alternative fuels is being conducted due to their increasing demand to reduce greenhouse emissions. One renewable fuel studied in this work is dimethyl ether (DME) blended with propane(C3H8) as a potential mixture for heavy-duty engines used in semi-trucks. The blend has the potential to drastically reduce particulate and greenhouse gas emissions compared to a conventional diesel engine operating under similar conditions. To develop the use of mixture, one must conduct detailed conceptual and simulation studies before progressing to detail studies in CFD, engine modifications, and live testing. For simulations, accurate high-fidelity chemical kinetic models are necessary. However, the validity of the chemical kinetic mechanism for operating conditions of a heavy-duty mixing-controlled compression (MCCI) engine was widely unknown until recent work presented here and published. In this work, we studied the ignition of DME and propane blends in a shock tube under MCCI engine conditions. Ignition delay time (IDT) gathered behind the reflected shock for DME-propane mixtures for heavy-duty compression ignition (CI) engine parameters. Testing was conducted for undiluted varieties spanning from temperatures of 700 to 1100 K at pressures ranging from 55 to 84 bar for various blends (100% CH3OCH3, 100% C3H8, 60% CH3OCH3/ 40% C3H8) of DME and propane were combusted in synthetic air (21% O2/ 79% N2). Several experiments were conducted at higher pressures (90-120 bar) to improve the model performance and accuracy. The ignition delay times (IDTs) were compared to recent mechanisms, including Aramco3.0, NUIG, and Dames et al. A common trend among the mechanisms was overpredicted experimental IDTs. Further studies were conducted by a sensitivity analysis using the Dames et al. model, and critical reactions sensitive to IDTs of DME-propane mixture near 60 bar are outlined. Chemical analysis was conducted on the NTC region to explain chemical kinetics which is critical for developing MCCI heavy duty engines.

Mohammed, Zuhayr Pasha↗

NASA NTRS1987

Grain dynamics in zero gravity

The dynamics of granular materials has proved difficult to model, primarily because of the complications arising from inelastic losses, friction, packing, and the effect of many grains being in contact simultaneously. One interesting limit for which it was recently possible to construct a theory is that where the grain-grain interactions are dominated by binary collisions. The kinetic model of granular systems if similar to the kinetic theory of gases, except that collisional energy losses are always present in the former and must be treated explicitly. Few granular materials on Earth are describable by this limiting model, since gravity tends to collapse the grains into a high-density state where Coulombic friction effects are dominant. The planned Space Station offers an unusual opportunity to test the kinetic grain model and to explore its predictions. Without gravity, the regime of low interparticle velocities, where an elastic description of the collision is still valid, is investigated. This will allow direct interpretation by dynamical computer simulations as well as by kinetic theory.

Werner, B. T.↗

DOE OSTI2026

Investigating the Adsorption–Desorption Kinetics of a Molecular Water Oxidation Catalyst at an Electrode Interface

Probing the dynamics of molecular catalysts at electrode–electrolyte interfaces is essential for understanding catalytic mechanisms. Structure-specific spectroscopic methods are particularly powerful for examining electrocatalytic interfaces but are mostly used under steady-state conditions. Herein, we combined surface-enhanced infrared absorption spectroscopy (SEIRAS) with phase-sensitive detection (PSD) to investigate the dynamics of a molecular Ir-based water oxidation catalyst at the Au–electrolyte interface. We found that the amplitude of the absorbance of the catalyst is anticorrelated to that of interfacial water. This anticorrelation can be understood by the adsorption of the electrooxidized catalyst on the electrode and concurrent displacement of interfacial water. The infrared signals from the interface exhibit an increasing phase lag with respect to the electrode potential with an increasing scan rate of the potential. Kinetic modeling suggests that the potential-dependent adsorption–desorption kinetics of the molecular catalyst on the electrode gives rise to this phase lag. Furthermore, this study shows that PSD-SEIRAS is a powerful tool for investigating the interfacial dynamics of electrocatalytic systems.

Catalysts↗

NASA NTRS1990

Fracture-tough, corrosion-resistant bearing steels

The fundamental principles allowing design of stainless bearing steels with enhanced toughness and stress corrosion resistance has involved both investigation of basic phenomena in model alloys and evaluation of a prototype bearing steel based on a conceptual design exercise. Progress in model studies has included a scanning Auger microprobe (SAM) study of the kinetics of interfacial segregation of embrittling impurities which compete with the kinetics of alloy carbide precipitation in secondary hardening steels. These results can define minimum allowable carbide precipitation rates and/or maximum allowable free impurity contents in these ultrahigh strength steels. Characterization of the prototype bearing steel designed to combine precipitated austenite transformation toughening with secondary hardening shows good agreement between predicted and observed solution treatment response including the nature of the high temperature carbides. An approximate equilibrium constraint applied in the preliminary design calculations to maintain a high martensitic temperature proved inadequate, and the solution treated alloy remained fully austenitic down to liquid nitrogen temperature rather than transforming above 200 C. The alloy can be martensitically transformed by cryogenic deformation, and material so processed will be studied further to test predicted carbide and austenite precipitation behavior. A mechanistically-based martensitic kinetic model was developed and parameters are being evaluated from available kinetic data to allow precise control of martensitic temperatures of high alloy steels in future designs. Preliminary calculations incorporating the prototype stability results suggest that the transformation-toughened secondary-hardening martensitic-stainless design concept is still viable, but may require lowering Cr content to 9 wt. pct. and adding 0.5 to 1.0 wt. pct. Al. An alternative design approach based on strain-induced martensitic transformation during cryogenic forming, thus removing the high martensitic constraint, may permit alloy compositions offering higher fracture roughness.

Olson, Gregory B.↗

NASA NTRS1994

Direct numerical simulations of turbulent non-premixed methane-air flames modeled with reduced kinetics

Turbulent non-premixed stoichiometric methane-air flames modeled with reduced kinetics have been studied using the direct numerical simulation approach. The simulations include realistic chemical kinetics, and the molecular transport is modeled with constant Lewis numbers for individual species. The effect of turbulence on the internal flame structure and extinction characteristics of methane-air flames is evaluated. Consistent with earlier DNS with simple one-step chemistry, the flame is wrinkled and in some regions extinguished by the turbulence, while the turbulence is weakened in the vicinity of the flame due to a combination of dilatation and an increase in kinematic viscosity. Unlike previous results, reignition is observed in the present simulations. Lewis number effects are important in determining the local stoichiometry of the flame. The results presented in this work are preliminary but demonstrate the feasibility of incorporating reduced kinetics for the oxidation of methane with direct numerical simulations of homogeneous turbulence to evaluate the limitations of various levels of reduction in the kinetics and to address the formation of thermal and prompt NO(x).

Card, J. M.↗

NASA NTRS1970

A mathematical model for the thermoradiation inactivation of dry Bacillus subtilis var. Niger spores

Empirically based kinetic model describing synergistic inactivation of dry Bacillus subtilis by combined heat and gamma radiation environment

Dugan, V. L.↗

NASA NTRS1977

Bhatnagar-Gross-Krook model solution of back-scattering of outgas flow from spherical spacecraft

The BGK kinetic model for a binary mixture of gases was used to predict the backscattering return flux of molecules outgassing from a spherical spacecraft. For simplicity, the outgassing molecules were considered to be a single species, with the ambient atmospheric molecules being the second species in the binary mixture. Two basic collision mechanisms were considered: 1) collisions involving only outgas molecules (self-scattering); and 2) collisions between outgas molecules and molecules in the ambient atmosphere (ambient-scattering). Results were obtained in the form of simple equations relating outgas return flux to parameters based on outgas characteristics, spacecraft size, and the characteristics of the ambient flow.

Robertson, S. J.↗

DOE OSTI2025

Experimental and modeling study of the autoignition behavior of a saturated heterocycle: Pyrrolidine

Experiments are conducted in both rapid compression machine (RCM) and shock tube (ST) to better quantify autoignition behavior (e.g., ignition delay, heat release) and understand heteroatomic effects in heterocyclic compounds, which are important reference components for the combustion of biomass-derived liquid fuels. These tests focus on the nitrogen-containing, five-membered saturated ring, pyrrolidine, at diluted conditions covering pressures of 20 and 50 bar, temperatures of 720–1450 K and a range of stoichiometries (ϕ = 0.5–2). A chemical kinetic model is developed and coupled to an existing combustion kinetics framework describing key nitrogen containing intermediates (e.g. pyrrole, ammonia and NOx). H-abstraction reactions by OH, H, CH 3 and HO 2 , are determined using ab-initio transition state theory methods, while analogies to cyclopentane are adopted for many other reactions, such as ring-opening. Here, the autoignition measurements reveal the lack of negative temperature coefficient (NTC) behavior and low-temperature chemistry for pyrrolidine, as opposed to its saturated hydrocarbon analogue, cyclopentane. Interestingly, at the lowest temperatures (T < 750 K), the reactivity of cyclopentane is greater than pyrrolidine, while at higher temperatures, pyrrolidine becomes more reactive. Agreement between the experimental measurements and the model is good, and it is found that H-abstraction reactions by HO 2 and ensuing chemistry play key roles in controlling the reactivity of this cyclic amine. Most of the fluxes, i.e., >70 %, are predicted to move through 1- or 2-pyrroline (C 4 H 7 N) and then the cyclic C 4 H 6 N radical, at both lower and higher temperatures, to form either CH 2 CHCHCHNH via ring-opening or pyrrole via β-scission. It appears that the ring opens more easily at lower temperature whereas the C–H β-scission dominates at higher temperature and lower pressure, such that the reaction of the fuel radical intermediate carrying an unpaired electron on the nitrogen atom with HO 2 is the next most notable in promoting oxidation. When comparing pyrrolidine and cyclopentane, which exhibits distinct pathways in different temperature regimes, the pyrrolidine pathways and sensitivity analysis align more closely to the high temperature case of cyclopentane where the important role of HO 2 radicals is seen to provide chain branching through HO 2 reaction with the fuel, accompanied by H 2 O 2 formation and decomposition to OH. The formation of 5-membered diene rings and ring opening reactions are also found to be highly relevant. Of particular note, it is found that there is little influence of small molecule nitrogen-chemistry, e.g., NH 2 , HCN, NO/NO 2 on the reactivity of the pyrrolidine mixtures investigated here where no recirculated combustion gases are included.

Autoignition↗

NASA NTRS2002

Fundamentals of Plasma Sails Propulsion Concept

The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000], is based on the two-fluid plasma model and requires a 15-km frontal standoff distance (or 20-km cross-sectional diameter) in order for the magnetic bubble to absorb sufficient momentum from the SW to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Winglee et al. [2000] derived the above size requirement based on an extrapolation of their simulated results in which a system much smaller than a M2P2 was used (p. 21,074 of their study). We submit, however, that a fluid model has no validity for such a small scale size-even in the region near the plasma source! It is assumed in the MHD fluid model, normally applied to the magnetosphere, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the SW. In the case of the M2P2, however, the size of the magnetic bubble is actually less than or, at best, comparable to, the scale of these characteristic parameters and, therefore, a kinetic approach, which addresses the smallscale physical mechanisms involved, must be used. A fully three-dimensional version of the hybrid code is used in our M2P2 (Plasma Sails) studies was originally developed by Delamere et al. [1999]. The M2P2 plasma sail is an excellent application for this hybrid code. The primary advantage of this code is the seamless interface between fluid and kinetic descriptions of the ion populations. A kinetic description is not necessary for the dense inner regions of the magnetic bubble and tremendous computational savings can be realized by treating this dense, magnetized ion population with the fluid description. It is essential, however, that the outer bubble regions be treated kinetically as well as the SW protons. Comparison of full size M2P2 simulation based on 3D MHD and kinetic models show that kinetic treatment introduces much more asymmetry to the considering problem and the possibility of kinetic instabilities development.

Khazanov, G. V.↗

NASA NTRS2024

Fundamental Experimental Tests and Modeling of LOX/CH4 Engines at High Pressures

In recent years, the space travel industry has grown exponentially, resulting in the need for a low-cost, efficient rocket engine fuel. As such, there has been renewed interest in utilizing liquid natural gas (LNG) as it is less likely to soot than kerosene-based fuels and has widespread availability and low cost of use as compared to the traditional RP-X and even liquid methane (LCH4) fuels due to the reduced need for refinement. However, current literature does not fully cover varying blends of natural gas and the effects that impurities have on natural gas/methane fuels, with most studies being confined to low-pressure applications. Therefore, new experimental ignition delay time measurements at rocket engine relevant pressures are essential to ensure current chemical kinetic models capture the ignition behavior at these elevated conditions and to determine the optimal purity levels for a reliable, cost-efficient, aerospace-grade LCH4/LNG fuel. The current work explores various blends of natural gas/methane fuel with various impurities, including higher hydrocarbons as well as nitrogen-carriers. A shock tube study was carried out to study the ignition delay times and carbon monoxide time histories of these natural gas blends utilizing chemiluminescence and laser absorption spectroscopy at 20, 50 and 100 atm over temperatures ranging from 1400-1700K and at an equivalence ratio of 1. The experimental data was then compared to two different chemical kinetic mechanisms: the industry-standard GRI 3.0 and the in-house UCF 2022. The data is used to improve chemical kinetic mechanisms and modeling of LNG fuels for rocket engine applications.

Methane↗

DOE OSTI2025

Active Site Dynamics in Molybdenum-Based Silica-Supported Olefin Metathesis Catalysts: Site Renewal and Decay Beyond the Chauvin Cycle

Heterogeneous olefin metathesis catalysts exhibit low active site densities and unpredictable kinetics due to dynamic active site formation and decay processes. Here, in this study, we establish a quantitative framework that captures active site generation, renewal, and decay in olefin metathesis over silica-supported molybdenum oxide catalysts, enabling a mechanistic explanation of catalytic behavior and strategies to achieve high, stable activity. Steady-state active site titrations reveal that 2,3-dimethyl-butene isomers (4MEs) cofeeding increases active site density by up to 4.3-fold, directly correlating with enhanced metathesis rates. Spectroscopic studies demonstrate that 4MEs facilitate Mo(VI) reduction to Mo(IV) and interact strongly with surface Si–OH groups, generating labile protons that promote active site formation via a 1,2-proton shift mechanism. Kinetic modeling indicates that ethylene acts as a decay promoter, shifting kinetic control away from the Chauvin cycle and suppressing metathesis activity. Comparative studies on catalysts with varying Mo loading reveal that promotion is most effective for dispersed molybdate species, with a decline at higher Mo loadings. These findings provide a unified mechanistic framework for heterogeneous olefin metathesis, offering new strategies to enhance active site accessibility, mitigate deactivation, and optimize catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

New physics in multi-lepton tau decays

Dark particles with lepton-flavor-violating couplings to the tau lepton can induce rare neutrinoless $τ$ decays with large final state multiplicities. We study models where transitions of the type $τ^\pm\to \ell^\pm\,ϕ$, with $ϕ$ a light new particle, initiate a chain of decays in the dark sector that terminate with decays into electrons, muons, or pions. These decay cascades appear as rare five or even seven-body $τ$ decays with multiple reconstructable resonances. We survey several representative models: kinetically mixed dark photon, gauged $L_i-L_j$ models, and other more exotic charge assignments such as chiral $U(1)'$ extensions of the Standard Model. The main new ingredient is the possibility of flavor violation at very high scales. In these models, a number of channels that have not yet been searched for experimentally, such as $τ\to 5μ$, $τ\to 3μ\,2e$, $τ\to μ\,4e$, and hadronic channels like $τ\to μ\,4π$, typically dominate over the previously-considered signatures such as $τ\to 3μ$. While some of the models, such as the gauged $L_i-L_j$ ones, also contain more challenging channels with missing energy due to decays to neutrinos, they can still be searched for via fully visible channels.

Ema, Yohei [Florida U.]↗

NASA NTRS2002

Origins of Solar Systems Program

O'D. Alexander, Conel↗

NASA NTRS2003

Final Summary of Research Report to the National Aeronautics and Space Administration Origins of Solar Systems Program

The chondrites are aggregates of components (e.g. chondrules, chondrule rims and matrix) that formed in the nebula but, at present, there is no consensus on how any of these components formed or whether their formation produced or post dated the chemical fractionations between the chondrites. Chondrites are, at present, the most primitive Solar System objects available for laboratory study and the conditions under which their principle components formed would provide the most direct constraints for models of nebula formation and evolution. The conditions under which chondrules formed is of particular importance because, if their relative abundance in chondrites approximates that in the nebula, they are the products of one of the most energetic and pervasive processes that operated in the early Solar System. The goal of this proposal was to combine theoretical modeling with a comprehensive study of the elemental and isotopic compositions of the major components in unequilibrated ordinary chondrites (UOCs), with the aim of determining the conditions in the nebula at the time of their formation. The isotopes of volatile and moderately volatile elements should be particularly revealing of conditions during chondrule formation, as evaporation under most conditions would lead to isotopic mass fractionation. Modeling of chondrule and matrix formation requires the development of a kinetic model of evaporation and condensation, and calibration of this model against experiments. Cosmic spherules present an opportunity to test our evaporation models under flash heating conditions that would be difficult to simulate experimentally. However, there is surprisingly little known about the isotopic compositions of silicate cosmic spherules, and a number of questions need to be addressed. Is the range of compositions they exhibit due to evaporation? If they are, are the relative volatilities consistent with the models/experiments and are the isotopic fractionations consistent with Rayleigh conditions? For instance, do the alkalis and S evaporate prior to significant melting so that conditions did not meet the Rayleigh criteria of rapid diffusion? If so, their isotopic fractionation might be considerably suppressed. Could this mechanism of K loss apply to chondrule formation? The Fe isotopic fractionation during evaporation of silicates has not been measured, so cosmic spherules might provide a clue to whether FeO diffusion is fast enough to maintain Rayleigh conditions during evaporation. And so on.

O'D. Alexander, Conel↗