Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic modeling”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

N 2 vibrational excitation in atmospheric pressure ns pulse and RF plasma jets

Time-resolved N 2 vibrational temperature and translational–rotational temperature in quasi-two-dimensional atmospheric pressure plasma jets sustained by ns pulse and RF discharges in nitrogen/noble gas mixtures are measured by the broadband vibrational Coherent Anti-Stokes Raman Scattering (CARS) . The results indicate a much stronger vibrational excitation in the RF plasma jet, due to the lower reduced electric field and higher discharge power. In a ns pulse discharge in N 2 /He, N 2 vibrational temperature is significantly lower compared to that in N 2 /Ar, due to the more rapid vibration–translation (V–T) relaxation of nitrogen by helium atoms. In the RF plasma jets in N 2 /Ne and N 2 /Ar, the vibrational excitation increases considerably as the nitrogen fraction in the mixture is reduced. The experimental data in the RF plasma jet in N 2 /Ar jet are compared with the kinetic modeling predictions. The results indicate that nitrogen vibrational excitation in N 2 /Ar plasma jets with a small N 2 fraction in the mixture (several percent) is controlled primarily by electron impact, anharmonic vibration–vibration (V–V) pumping, and V–T relaxation by N atoms. In comparison, V–V energy transfer from the vibrationally excited molecules in the first excited electronic state, N 2 (A 3 Σ u + , v), which are generated primarily by the energy transfer from the metastable Ar atoms, has a minor effect on the vibrational populations of the ground electronic state, N 2 (X 1 Σ g + , v). Although the discharge energy fraction going to electronic excitation is significant, the predicted quasi-steady-state N 2 (A 3 Σ u + ) number density, controlled by the energy pooling and quenching by N atoms, remains relatively low. Because of this, the net rate of N 2 (X 1 Σ g + ) vibrational excitation by the V–V energy transfer from N 2 (A 3 Σ u + ) is much lower compared to that by the direct electron impact. Furthermore, the results show that atmospheric pressure RF plasma jets can be used as sources of highly vibrationally excited N 2 molecules and N atoms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
DOE OSTI2024

Thermodynamic control on the decomposition of organic matter across different electron acceptors

The increasing availability of high-resolution characterization of natural organic matter (OM) data has shifted the paradigm of lumped descriptions of OM components and potential microbial activities. Our recent development of a substrate-explicit thermodynamic model uniquely enables incorporating complex OM pools to formulate biogeochemical reaction models based on their elemental compositions. While this previous work facilitates prediction of aerobic respiration of complex OM, it is equally imperative to consider the role of non-oxygenic electron acceptors in regulating OM turnover and the fate of carbon. In this study, we significantly expand our previous model by flexibly incorporating both detailed OM chemistry and electron acceptors other than oxygen. Here, our modeling analysis has revealed substantial variations in the energy status of OM molecules across different soils, which drive the co-occurrence of different electron-accepting processes. We demonstrated the effectiveness of the proposed model using a consistency check with experimental data. Through systematic evaluation of the impact of diverse chemical inputs (both electron donors and acceptors) on OM decomposition, the new model also revealed how key microbial growth parameters such as carbon use efficiency (CUE) and reaction rates vary across different electron-accepting processes. Our model provides a unified framework integrating thermodynamic and kinetic constraints on microbial metabolic activities. It complements traditional kinetic models, which are often designed solely to capture mass fluxes. We conclude that thermodynamic modeling emerges as a powerful tool for describing the mechanisms underlying the interplay between microbial growth and OM chemistry and cycling across different electron acceptors, enhancing our ability to project complex ecosystem behaviors in dynamic environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

High-pressure CO 2 dissociation with nanosecond pulsed discharges

Herein we investigate the conversion of CO 2 into CO and O 2 with nanosecond repetitively pulsed (NRP) discharges in a high-pressure batch reactor. Stable discharges are obtained at up to 12 bar. By-products are measured with gas chromatography. The energy efficiency is determined for a range of processing times, pulse energy, and fill pressures. It is only weakly sensitive to the plasma operating parameters, i.e the extent of CO 2 conversion is almost linearly-dependent on the specific energy invested. A conversion rate as high as 14% is achieved with an energy efficiency of 23%. For long processing times, saturation in the yield and a drop in efficiency are observed, due to the increasing role of three-body recombination reactions, as described by zero-dimensional detailed kinetic modeling. The modeling reveals the presence of three-stage kinetics between NRP pulses, controlled by electron-impact CO 2 dissociation, vibrational relaxation, and neutral elementary kinetics. Transport effects are shown to be important for CO 2 conversion at high pressures. For fill pressures beyond 10 bar, CO 2 may locally transit into supercritical states. The supercritical plasma kinetics may bypass atomic oxygen pathways and directly convert CO 2 into O 2 . This work provides a detailed analysis of plasma-based high-pressure CO 2 conversion, which is of great relevance to future large-scale sustainable carbon capture, utilization, and storage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

LOKI

LOKI is a simulation code that models weakly-coupled plasma evolution using Eulerian kinetic models in four phase space dimensions. The code is useful to understand basic plasma nonlinear plasma physics phenomena - in particular, wave-particle and wave-wave resonances that occur in laser plasma interactions. LOKI can simulate fundamental physical process in low-density plasmas on a scale of hundreds of light wavelengths. It does so by solving the relativistic Vlasov and Valsov-Fokker-Planck differential transport equations describing the collision less and low-collisionality evolution of charged particles in electromagnetic fields in two configuration space and two velocity space dimensions. Two principle electromagnetic models are provided: electrostatic (Poisson's Equation) and electromagnetic (Maxwell's equations).

DIMITS, ANDRIS↗
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 OSTI2022

2D Reconstruction of Magnetotail Electron Diffusion Region Measured by MMS

Models for collisionless magnetic reconnection in near-Earth space are distinctly characterized as 2D or 3D. In 2D kinetic models, the frozen-in law for the electron fluid is usually broken by laminar dynamics involving structures set by the electron orbit size, while in 3D models the width of the electron diffusion region is broadened by turbulent effects. We present an analysis of in situ spacecraft observations from the Earth's magnetotail of a fortuitous encounter with an active reconnection region, mapping the observations onto a 2D spatial domain. While the event likely was perturbed by low-frequency 3D dynamics, the structure of the electron diffusion region remains consistent with results from a 2D kinetic simulation. As such, the event represents a unique validation of 2D kinetic, and laminar reconnection models.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
DOE OSTI2021

Understanding copper diffusion in CuInSe 2 with first-principles based atomistic and continuum models

We investigate the diffusion of copper in CuInSe 2 using thermodynamic and kinetic models based on density functional theory calculations, attempting to reconcile large differences in reported experimental diffusivities. We find that observations of rapid chemical diffusion can be explained by large thermodynamic factors, which we calculate using a compositionally constrained model of intrinsic point defect formation. We further characterize how copper diffusion coefficients depend on material synthesis conditions and exhibit their variation across the CuInSe 2 secondary phase diagram. In doing so, we identify stable off-stoichiometries that are dominated by either vacancy- or interstitial-mediated diffusion mechanisms. These results are employed in the development of a continuum reaction–diffusion model, which we use to simulate experimental depth profiles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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.↗
DOE OSTI2021

Reaction Mechanism Generator v3.0: Advances in Automatic Mechanism Generation

In chemical kinetics research, kinetic models containing hundreds of species and tens of thousands of elementary reactions are commonly used to understand and predict the behavior of reactive chemical systems. Reaction Mechanism Generator (RMG) is a software suite developed to automatically generate such models by incorporating and extrapolating from a database of known thermochemical and kinetic parameters. Here, we present the recent version 3 release of RMG and highlight improvements since the previously published description of RMG v1.0. Most notably, RMG can now generate heterogeneous catalysis models in addition to the previously available gas- and liquid-phase capabilities. For model analysis, new methods for local and global uncertainty analysis have been implemented to supplement first-order sensitivity analysis. The RMG database of thermochemical and kinetic parameters has been significantly expanded to cover more types of chemistry. The present release includes parallelization for faster model generation and a new molecule isomorphism approach to improve computational performance. RMG has also been updated to use Python 3, ensuring compatibility with the latest cheminformatics and machine learning packages. Overall, RMG v3.0 includes many changes which improve the accuracy of the generated chemical mechanisms and allow for exploration of a wider range of chemical systems.

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

Heat Integration Optimization and Dynamic Modeling Investigation for Advancing the Coal-Direct Chemical Looping Process

The purpose of the project is to address the optimization and startup operation of a modular coal direct chemical looping (CDCL) combustion system integrated with a steam cycle for power generation to reduce the risks involved in further scale-up of the technology. The modular reactor design of the CDCL process provides flexibility in the fabrication of the reactor and in its operating capacity (i.e. turndown ratio) at the cost of a more complex heat exchange network (HEN) design and integration. To address the technology gaps and advance the efficiency and economic feasibility of the CDCL technology, the project will perform a detailed and comprehensive analysis of the integration of a modular CDCL reactor system and a steam cycle system under both static and transient conditions via HEN process performance simulations and system dynamic modeling, respectively. The scope of work consists of 1) Experimental and computational studies of the CDCL combustor reactor 2) Comprehensive static (i.e. steady-state) system HEN design analysis in CDCL 550 MWe commercial unit for power generation and 3) Dynamic modeling of site specific design of 10MWe CDCL large pilot plant. The project team has successfully developed and validated a kinetic model for the oxidation of oxygen carriers in the combustor using the unreacted shrinking core model (UCSM). The model is capable of capturing the oxidation kinetics of fully or partially reduced oxygen carrier particles. A computational fluid dynamics (CFD) model is developed to simulate the hydrodynamics, heat transfer, and chemical reaction occurring in the CDCL combustor. The model is developed in MFIX and ANSYS Fluent. Key aspects of CDCL combustor operation, including heat transfer, oxygen carrier oxidation, and the transport of oxygen carrier particles, are simulated using this CFD model. The HEN for a commercial scale 550 MWe CDCL power plant is simulated and optimized using ASPEN Plus. Practical design considerations are incorporated based on industrial experiences. The performance and cost for the commercial CDCL plant is updated based on these analyses. A dynamic model for the 10 MWe CDCL pilot plant is developed in ProTRAX simulation software. The model is based on the pilot plant design developed in project DE-FE0027654 “10 MWe CDCL Large Pilot Plang – Pre-FEED Study” and the steam cycle data obtained from Dover Light & Power plant. The transient behaviors during pilot plant load variation are simulated using the dynamic model.

01 COAL, LIGNITE, AND PEAT↗
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.↗