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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Accelerated Aging of Lead-Free Propellant

Following higher than expected 2-NDPA depletion rates in a lead-free doublebase formulation (RPD-422), an accelerated aging study was conducted to verify the depletion rates. A test plan was prepared to compare the aging characteristics of lead-free propellant and NOSIH-AA2. The study was also designed to determine which lead-free ballistic modifiers accelerated 2-NDPA depletion. The increased depletion rate occurred in propellants containing monobasic copper salicylate. Four lead-free propellants were then formulated to improved aging characteristics over previous lead-free propellant formulations. The new formulations reduced or replaced the monobasic copper salicylate. The new formulations had improved aging characteristics. Their burn rates, however, were unacceptable for use in a 2.75 inch rocket. To compare aging characteristics, stabilizer depletion rates of RPD-422, AA2, M28, and RLC 470/6A were measured or taken from the literature. The data were fit to a kinetic model. The model contained first and zero order terms which allowed the stabilizer concentration to go to zero. In the model, only the concentration of the primary stabilizer was considered. Derivatives beyond the first nitrated or nitroso derivative of 2-NPDA were not considered. The rate constants were fit to the Arrhenius equation and extrapolated to lower temperatures. The time to complete stabilizer depletion was estimated using the kinetic model. The four propellants were compared and the RPD-422 depleted faster at 45 C than both A22 and M28. These types of predictions depend on the validity of the model and on confidence in the Arrhenius relationship holding at lower temperatures. At 45 C, the zero order portion of the model dominates the depletion rate.

Furrow, Keith W.↗
NASA NTRS1985

Direct solar-pumped iodine laser amplifier

This semiannual progress report covers the period from April 1, 1985 to Sept. 30, 1985 under NASA grant NAS1-441 entitled direct solar pumped iodine laser amplifier. During this period the parametric studies of the iodine laser oscillator pumped by a Vortek simulator was carried out before the amplifier studies. The amplifier studies are postponed to the extended period following completion of the parametric studies. In addition, the kinetic modeling of a solar pumped iodine laser amplifier, and the experimental work for a solar pumped dye laser amplifier are in progress. This report contains three parts: (1) the radiation characteristics of solar simulator and the parametric characteristics of photodissociation iodine laser continuously pumped by a Vortek solar simulator; (2) kinetic modeling of a solar pumped iodine laser amplifier; and (3) the study of the dye laser amplifier pumped by a Tamarack solar simulator.

Han, K. S.↗
NASA NTRS2009

Nitric Oxide and Oxygen Air-Contamination Effects on Extinction Limits of Non-Premixed Hydrocarbon-Air Flames for a HIFiRE Scramjet

Unique nitric oxide (NO) and oxygen air-contamination effects on the extinction Flame Strength (FS) of non-premixed hydrocarbon (HC) vs. air flames are characterized for 7 gaseous HCs, using a new idealized 9.3 mm straight-tube Opposed Jet Burner (OJB) at 1 atm. FS represents a laminar strain-induced extinction limit based on cross-section-average air jet velocity, Uair, that sustains combustion of a counter jet of gaseous fuel just before extinction. Besides ethane, propane, butane, and propylene, the HCs include ethylene, methane, and a 64 mole-% ethylene / 36 % methane mixture, the writer s previously recommended gaseous surrogate fuel for HIFiRE scramjet tests. The HC vs. clean air part of the work is an extension of a May 2008 JANNAF paper that characterized surrogates for the HIFiRE project that should mimic the flameholding of reformed (thermally- or catalytically-cracked) endothermic JP-like fuels. The new FS data for 7 HCs vs. clean air are thus consolidated with the previously validated data, normalized to absolute (local) axial-input strain rates, and co-plotted on a dual kinetically dominated reactivity scale. Excellent agreement with the prior data is obtained for all 7 fuels. Detailed comparisons are also made with recently published (Univ. Va) numerical results for ethylene extinction. A 2009-revised ethylene kinetic model (Univ. Southern Cal) led to predicted limits within approx. 5 % (compared to 45 %, earlier) of this writer s 2008 (and present) ethylene FSs, and also with recent independent data (Univ. Va) obtained on a new OJB system. These +/- 5 % agreements, and a hoped-for "near-identically-performing" reduced kinetics model, would greatly enhance the capability for accurate numerical simulations of surrogate HC flameholding in scramjets. The measured air-contamination effects on normalized FS extinction limits are projected to assess ongoing Arc-Heater-induced "facility test effects" of NO production (e.g., 3 mole-%) and resultant oxygen depletion (from 21 to 19.5 %), for testing the "64/36" surrogate fuel in Langley s Arc-Heated Scramjet Test Facility for HIFiRE engine designs. The FS results show a generally small (< 4 %) "nitric oxide enhancement" effect, relative to clean air, for up to 3 % NO (freestream Mach number up to 7 in Arc Jet testing). However, a progressively large "oxygendeficiency weakening" effect develops. For 3 % NO, a net weakening of 26 % in FS is derived for the "64/36" fuel vs. air. The corresponding net weakening for pure ethylene is 20 %. A number of practical recommendations regarding facility test effects are offered.

Pellett, Gerald L.↗
NASA NTRS2011

Experimental and Modeling Studies of the Combustion Characteristics of Conventional and Alternative Jet Fuels. Final Report

The objectives of this project have been to develop a comprehensive set of fundamental data regarding the combustion behavior of jet fuels and appropriately associated model fuels. Based on the fundamental study results, an auxiliary objective was to identify differentiating characteristics of molecular fuel components that can be used to explain different fuel behavior and that may ultimately be used in the planning and design of optimal fuel-production processes. The fuels studied in this project were Fischer-Tropsch (F-T) fuels and biomass-derived jet fuels that meet certain specifications of currently used jet propulsion applications. Prior to this project, there were no systematic experimental flame data available for such fuels. One of the key goals has been to generate such data, and to use this data in developing and verifying effective kinetic models. The models have then been reduced through automated means to enable multidimensional simulation of the combustion characteristics of such fuels in real combustors. Such reliable kinetic models, validated against fundamental data derived from laminar flames using idealized flow models, are key to the development and design of optimal combustors and fuels. The models provide direct information about the relative contribution of different molecular constituents to the fuel performance and can be used to assess both combustion and emissions characteristics.

Meeks, Ellen↗
DOE OSTI2022

Oxyanion Surface Complexes Control the Kinetics and Pathway of Ferrihydrite Transformation to Goethite and Hematite

The rate and pathway of ferrihydrite (Fh) transformation at oxic conditions to more stable products is controlled largely by temperature, pH, and the presence of other ions in the system such as nitrate (NO 3 - ), sulfate (SO 4 2- ), and arsenate (AsO 4 3- ). Although the mechanism of Fh transformation and oxyanion complexation have been separately studied, the effect of surface complex type and strength on the rate and pathway remains only partly understood. We have developed a kinetic model that describes the effects of surface complex type and strength on Fh transformation to goethite (Gt) and hematite (Hm). Two sets of oxyanion-adsorbed Fh samples were prepared, nonbuffered and buffered, aged at 70 ± 1.5 °C, and then characterized using synchrotron X-ray scattering methods and wet chemical analysis. Kinetic modeling showed a significant decrease in the rate of Fh transformation for oxyanion surface complexes dominated by strong inner-sphere (SO 4 2- and AsO 4 3- ) versus weak outer-sphere (NO 3 - ) bonding and the control. The results also showed that the Fh transformation pathway is influenced by the type of surface complex such that with increasing strength of bonding, a smaller fraction of Gt forms compared with Hm. In conclusion, these findings are important for understanding and predicting the role of Fh in controlling the transport and fate of metal and metalloid oxyanions in natural and applied systems.

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

Modeling electron beam irradiation of methane

We built a chemical kinetic model for the electron beam irradiation of methane and report simulation results with dose rates ranging from 2 to 200 kGy s -1 at 298 K and 1 atm pressure. The model contains 92 species and 778 gas-phase reactions. We report the time evolution of concentrations of stable products (including H 2 , C 2 H 6 , C 3 H 8 , C 2 H 4 , C 3 H 6 , C 4 H 10 , C 2 H 2 , and C 3 H 4 ) as well as some radicals and ions. Here we discuss how yields of products depend on dose rate and how they change during the course of irradiation. To understand the processes controlling product formation, we quantify the extent to which various reactions produce and destroy major products and the key radicals and ions. Normalized sensitivity coefficients are computed for the production of stable products to gain further insight into factors controlling product yields. We were able to reduce the mechanism to 63 species and 331 reactions while incurring less than 5% error in the final concentrations of stable products.

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

Thermophysics Characterization of Kerosene Combustion

A one-formula surrogate fuel formulation and its quasi-global combustion kinetics model are developed to support the design of injectors and thrust chambers of kerosene-fueled rocket engines. This surrogate fuel model depicts a fuel blend that properly represents the general physical and chemical properties of kerosene. The accompanying gaseous-phase thermodynamics of the surrogate fuel is anchored with the heat of formation of kerosene and verified by comparing a series of one-dimensional rocket thrust chamber calculations. The quasi-global combustion kinetics model consists of several global steps for parent fuel decomposition, soot formation, and soot oxidation, and a detailed wet-CO mechanism. The final thermophysics formulations are incorporated with a computational fluid dynamics model for prediction of the combustor efficiency of an uni-element, tri-propellant combustor and the radiation of a kerosene-fueled thruster plume. The model predictions agreed reasonably well with those of the tests.

Wang, Ten-See↗
NASA NTRS2001

Thermophysics Characterization of Kerosene Combustion

A one-formula surrogate fuel formulation and its quasi-global combustion kinetics model are developed to support the design of injectors and thrust chambers of kerosene-fueled rocket engines. This surrogate fuel model depicts a fuel blend that properly represents the general physical and chemical properties of kerosene. The accompanying gaseous-phase thermodynamics of the surrogate fuel is anchored with the heat of formation of kerosene and verified by comparing a series of one-dimensional rocket thrust chamber calculations. The quasi-global combustion kinetics model consists of several global steps for parent fuel decomposition, soot formation, and soot oxidation and a detailed wet-CO mechanism to complete the combustion process. The final thermophysics formulations are incorporated with a computational fluid dynamics model for prediction of the combustion efficiency of an unielement, tripropellant combustor and the radiation of a kerosene-fueled thruster plume. The model predictions agreed reasonably well with those of the tests.

Wang, Ten-See↗
DOE OSTI2022

Time-resolved CO 2 , CO, and N 2 vibrational population measurements in Ns pulse discharge plasmas

Abstract Time-resolved CO 2 and N 2 vibrational populations and translational-rotational temperature are measured in a CO 2 –N 2 plasma sustained by a ns pulse discharge burst in plane-to-plane geometry. Time-resolved, absolute number density of CO generated in the plasma is also inferred from the experimental data. CO 2 and CO vibrational populations are measured by mid-IR, tunable quantum cascade laser absorption spectroscopy, and N 2 vibrational populations are measured by the ns broadband vibrational CARS. Transient excitation of N 2 and CO 2 asymmetric stretch vibrational energy modes is detected during the discharge burst. The time-resolved rate of CO generation does not correlate with N 2 or CO 2 ( ν 3 ) vibrational temperatures, indicating that CO 2 dissociation via the vibrational excitation is insignificant at the present conditions. The rate of CO generation decreases gradually during the discharge burst. The estimated specific energy cost of the CO product is close to that of N atoms in pure nitrogen, measured previously at similar operating conditions. Comparison of the experimental data with the kinetic modeling analysis indicates that CO 2 dissociation in collisions with electronically excited N 2 molecules is the dominant channel of CO generation at the present conditions, although the inferred CO yield in these processes is significantly lower than 1. The effect of vibrational energy transfer between N 2 and CO 2 on the plasma chemical processes is insignificant. The kinetic model underpredicts a rapid reduction of the N 2 and CO 2 ( ν 3 ) vibrational temperatures during the later half of the discharge burst and in the afterglow. V–T relaxation of N 2 by N and O atoms generated in the ns pulse discharge plasma does not affect the vibrational relaxation rate in a significant way. However, rapid V–T relaxation of CO 2 by O atoms has a significant effect on the relaxation rate. The difference between the experimental data and the modeling predictions may be due to the unknown scaling of the CO 2 –O V–T rates with the vibrational quantum number.

Physics↗
DOE OSTI2022

The influence of cooling rate on condensation of iron, aluminum, and uranium oxide nanoparticles

Fundamental observations of particle size distributions are needed to develop models that predict the fate and transport of radioactive materials in the atmosphere following a nuclear incident. The extent of material transport is influenced by the time scales of particle formation processes (e.g., condensation, coagulation). In this study, we investigated the influence of cooling time scales on size distributions of uranium, aluminum, and iron oxide particles that are synthesized separately under identical run conditions inside the controlled environment of an argon plasma flow reactor. Two distinct temperature distributions are imposed along the flow reactor by varying the argon flow rate downstream of the plasma torch. The vaporized reactants of uranium, aluminum, and iron are cooled from about 5000K to 1000K before they are collected on silicon wafers for ex situ scanning electron microscope analysis. The microscope images show that the sizes of the largest aluminum and iron oxide particles heavily depend on the cooling time scales, whereas significant size variation with cooling rate is not observed for uranium oxide particles. In addition, the size distribution of aluminum oxide particles exhibits the broadest range among all three metal oxides studied. We performed simulations of particle size distributions using a kinetic model that couples gas phase oxidation chemistry with particle formation processes, including nucleation, condensation, and coagulation. The model results demonstrate the strong sensitivity of particle size distribution to different cooling histories (i.e., temperature vs residence time) along the flow reactor. In conclusion, the kinetic model also helps identify directions for future research to improve the predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine learning force field model for kinetic Monte Carlo simulations of itinerant Ising magnets

Here, we present a scalable machine learning (ML) framework for large-scale kinetic Monte Carlo (kMC) simulations of itinerant electron Ising systems. As the effective interactions between Ising spins in such itinerant magnets are mediated by conducting electrons, the calculation of energy change due to a local spin update requires solving an electronic structure problem. Such repeated electronic structure calculations could be overwhelmingly prohibitive for large systems. Assuming the locality principle, a convolutional neural network (CNN) model is developed to directly predict the effective local field and the corresponding energy change associated with a given spin update based on Ising configuration in a finite neighborhood. As the kernel size of the CNN is fixed at a constant, the model can be directly scalable to kMC simulations of large lattices. Our approach is reminiscent of the ML force field models widely used in first-principles molecular dynamics simulations. Applying our ML framework to a square-lattice double-exchange Ising model, we uncover unusual coarsening of ferromagnetic domains at low temperatures. Our work highlights the potential of ML methods for large-scale modeling of similar itinerant systems with discrete dynamical variables.

machine learning↗
DOE OSTI2020

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

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

Non-premixed counterflow methane flames in DC/AC/NS electric fields

The response of counterflow diffusion flames to sub-breakdown DC and AC electric fields, as well as their superposition with ns pulse discharge waveforms, is studied in the plane-to-plane electrode geometry. Sub-breakdown DC and low-frequency AC electric fields cause the flame displacement toward the grounded electrode, in the direction of the applied field, indicating that the body force on the positive ions exceeds that on the electrons and negative ions. As the AC frequency increases, the flame response becomes less pronounced, due to the reduction of the electrohydrodynamic (EHD) body force impulse over the AC half-period. The electric field in the electrode gap is determined by ps Electric Field Induced Second Harmonic (E-FISH) generation, with absolute calibration using sub-breakdown ns pulses overlapped with the measured electric field waveform. The results show that the electric field distribution across the flame in the current saturation regime follows the Laplacian field. This indicates that the space charge density in the gap is too low to distort the applied DC or AC field, consistent with the kinetic modeling predictions. Combining a nanosecond pulse discharge with a sub-breakdown DC field generates a diffuse plasma across the entire gap. Time-resolved and spatially resolved measurements of the electric field in the discharge indicate the ionization wave propagation between the electrodes. The present results do not exhibit a detectable flame displacement enhancement by ns discharge pulses combined with a sub-breakdown field, observed previously. Kinetic modeling calculations show that the absence of this effect in the plane-to-plane geometry is due to the rapid plasma self-shielding. This indicates that alternative electrode geometries limiting the self-shielding would be more effective for the plasma / electric field enhanced flameholding and flame stabilization applications.

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

A Model for Radiolysis in a Flowing-Water Target during High-Intensity Proton Irradiation

At the Facility for Rare Isotope Beams (FRIB), interactions between heavy-ion beams and beam-dump water will create a wide variety of radionuclides which can be accessed by a technique known as “isotope harvesting”. However, irradiation of water is always accompanied by the creation of numerous radical, ionic, and molecular radiolysis products. Some of the radiolysis products have sufficiently long lifetimes to accumulate in the irradiated water and affect the harvesting chemistry. Here we investigate the formation of hydrogen peroxide, molecular hydrogen, and molecular oxygen during a high-intensity proton irradiation of a flowing-water isotope-harvesting target and compare the experimental results to simulations. The simulations kinetically model the chemical reactions occurring in the homogeneous phase of radiolysis in flowing water and establish an “effective yield”. In both the experiment and simulations, the bulk quantities of H 2 , H 2 O 2 , and O 2 are considerably lower than predicted by primary radiolysis yields (escape yields), meaning that in the high beam intensity regime the homogeneous phase reactions have a considerable impact on the overall chemical composition of the water. Further, it could be shown that for radiation which is characterized by a limited linear energy transfer, such as the here applied protons, the bulk outcome of the microscopic kinetic modeling could be estimated by a simplified steady-state model.

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

Aldehyde cool-flame chemistry explains a missing source of organic acids

Combustion emission is a significant source of organic acids, impacting atmospheric chemistry and climate. Their formation mechanisms, however, remain poorly understood, leading to underestimation in kinetic models. We investigate the cool-flame oxidation of key combustion intermediates—C 1 –C 4 aldehydes and benzaldehyde. Using in-situ synchrotron vacuum ultraviolet photoionization mass spectrometry, we observe the direct conversion of aldehydes to organic acids, a process enhanced by HO 2 radicals. Quantum chemistry calculations reveal that the reaction of RC(O)O 2 with HO 2 on the singlet potential energy surface contributes to organic acids. Incorporating this pathway into a kinetic model significantly improves organic acid prediction. Despite the high-temperature nature of engine combustion, significant spatial and temporal inhomogeneities (e.g., near-wall regions and crevice volumes) lead to localized cool-flame conditions, facilitating organic acid formation and emission. Elucidating the acid formation under cool-flame conditions provides a critical mechanism for accurately modelling anthropogenic organic acid emissions and developing mitigation strategies.

SVUV-PIMS↗
DOE OSTI2021

Revisiting trends in the exchange current for hydrogen evolution

Nørskov and collaborators proposed a simple kinetic model to explain the volcano relation for the hydrogen evolution reaction on transition metal surfaces such that j 0 = k 0 f(ΔG H ) where j 0 is the exchange current density, f(ΔG H ) is a function of the hydrogen adsorption free energy ΔG H as computed from density functional theory, and k 0 is a universal rate constant. Herein, focusing on the hydrogen evolution reaction in acidic medium, we revisit the original experimental data and find that the fidelity of this kinetic model can be significantly improved by invoking metal-dependence on k 0 such that the logarithm of k 0 linearly depends on the absolute value of ΔG H . Here, we further confirm this relationship using additional experimental data points obtained from a critical review of the available literature. Our analyses show that the new model decreases the discrepancy between calculated and experimental exchange current density values by up to four orders of magnitude. Furthermore, we show the model can be further improved using machine learning and statistical inference methods that integrate additional material properties.

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

Volume Diffusion Growth Kinetics and Step Geometry in Crystal Growth

The role of step geometry in two-dimensional stationary volume diff4sion process used in crystal growth kinetics models is investigated. Three different interface shapes: a) a planar interface, b) an equidistant hemispherical bumps train tAx interface, and c) a train of right angled steps, are used in this comparative study. The ratio of the super-saturation to the diffusive flux at the step position is used as a control parameter. The value of this parameter can vary as much as 50% for different geometries. An approximate analytical formula is derived for the right angled steps geometry. In addition to the kinetic models, this formula can be utilized in macrostep growth models. Finally, numerical modeling of the diffusive and convective transport for equidistant steps is conducted. In particular, the role of fluid flow resulting from the advancement of steps and its contribution to the transport of species to the steps is investigated.

Mazuruk, Konstantin↗
NASA NTRS2005

Heat Pipe Reactor Dynamic Response Tests: SAFE-100 Reactor Core Prototype

The SAFE-I00a test article at the NASA Marshall Space Flight Center was used to simulate a variety of potential reactor transients; the SAFEl00a is a resistively heated, stainless-steel heat-pipe (HP)-reactor core segment, coupled to a gas-flow heat exchanger (HX). For these transients the core power was controlled by a point kinetics model with reactivity feedback based on core average temperature; the neutron generation time and the temperature feedback coefficient are provided as model inputs. This type of non-nuclear test is expected to provide reasonable approximation of reactor transient behavior because reactivity feedback is very simple in a compact fast reactor (simple, negative, and relatively monotonic temperature feedback, caused mostly by thermal expansion) and calculations show there are no significant reactivity effects associated with fluid in the HP (the worth of the entire inventory of Na in the core is .<el, so fluid movement and temperature changes will cause very minor effects). In previous SAFE-100 tests, the point kinetics model was based on core thermal expansion via deflection measurements. It was found that core deflection was a strung function of how the SAFE-100 modules were fabricated and assembled (in terms of straightness, gaps, and other tolerances). To remove the added variable of how this particular core expands as compared to a different concept, it was decided to use a temperature based feedback model (based on several thermocouples placed throughout the core).

Bragg-Sitton, Shannon M.↗