Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Combustion modeling in internal combustion engines

The fundamental assumptions of the Blizard and Keck combustion model for internal combustion engines are examined and a generalization of that model is derived. The most significant feature of the model is that it permits the occurrence of unburned hydrocarbons in the thermodynamic-kinetic modeling of exhaust gases. The general formulas are evaluated in two specific cases that are likely to be significant in the applications of the model.

Zeleznik, F. J.↗
NASA NTRS1983

Solar-simulator-pumped atomic iodine laser kinetics

The literature contains broad ranges of disagreement in kinetic data for the atomic iodine laser. A kinetic model of a solar-simulator-pumped iodine laser is used to select those kinetic data consistent with recent laser experiments at the Langley Research Center. Analysis of the solar-simulator-pumped laser experiments resulted in the following estimates of rate coefficients: for alkyl radical (n-C3F7) and atomic iodine (I) recombination, 4.3 x 10 to the 11th power (1.9) + or - cu cm/s; for n-C3F7I stabilized atomic iodine recombination (I + I) 3.7 x 10 to the -32nd power (2.3) + or -1 cm to the 6th power/s; and for molecular iodine (I2) quenching, 3.1 x 10 to the -11th power (1.6) + or - 1 cu cm/s. These rates are consistent with the recent measurements.

Wilson, H. W.↗
NASA NTRS1989

Group-kinetic theory and modeling of atmospheric turbulence

A group kinetic method is developed for analyzing eddy transport properties and relaxation to equilibrium. The purpose is to derive the spectral structure of turbulence in incompressible and compressible media. Of particular interest are: direct and inverse cascade, boundary layer turbulence, Rossby wave turbulence, two phase turbulence; compressible turbulence, and soliton turbulence. Soliton turbulence can be found in large scale turbulence, turbulence connected with surface gravity waves and nonlinear propagation of acoustical and optical waves. By letting the pressure gradient represent the elementary interaction among fluid elements and by raising the Navier-Stokes equation to higher dimensionality, the master equation was obtained for the description of the microdynamical state of turbulence.

Tchen, C. M.↗
NASA NTRS1994

Quasiglobal reaction model for ethylene combustion

The objective of this study is to develop a reduced mechanism for ethylene oxidation. The authors are interested in a model with a minimum number of species and reactions that still models the chemistry with reasonable accuracy for the expected combustor conditions. The model will be validated by comparing the results to those calculated with a detailed kinetic model that has been validated against the experimental data.

Singh, D. J.↗
NASA NTRS2004

Computational and Experimental Study of Energetic Materials in a Counterflow Microgravity Environment

Counterflow diffusion flames are studied for various fuels flowing against decomposition products from solid ammonium perchlorate (AP) pellets in order to obtain fundamental understanding of composite propellant flame structure and chemistry. We illustrate this approach through a combined experimental and numerical study of a fuel mixture consisting of C2H4 CO + H2, and C2H2 + C2H4 flowing against solid AP. For these particular AP-fuel systems, the resulting flame zone simulates the various flame structures that are ex+ to exist between reaction products from Ap crystals and a hydrocarbon binder. As in all our experimental studies, quantitative species and temperature profiles have been measured between the fuel exit and AP surface. Species measured included CN, NH, NO, OH, N2, CO2, CO, H2, CO, HCl, and H2O. Temperature was measured using a thermocouple at the exit, spontaneous Raman scattering measurements throughout the flame, OH rotational population distributions, and NO vibrational population distributions. The burning rate of AP was also measured as a function of strain rate, given by the separation distance between the AP surface and the gaseous hydrocarbon fuel tube exit plane. This distance was nominally set at 5 mm, although studies have been performed for variations in separation distance. The measured 12 scalars are compared with predictions from a detailed gas-phase kinetics model consisting of 86 species and 531 reactions. Model predictions are found to be in good agreement with experiment and illustrate the type of kinetic features that may be expected to occur in propellants when AP particle size distributions are varied. Furthermore, the results constitute the continued development of a necessary database and validation of a comprehensive model for studying more complex AP-solid fuel systems in microgravity. Exploratory studies have also been performed with liquid and solid fuels at normal gravity. Because of melting (and hence dripping) and deep thermal wave penetration into the liquid, these experiments were found feasible, but not used for obtaining quantitative data. Microgravity experiments are needed to eliminate the dripping and boiling phenomena of these systems at normal gravity. Microgravity tests in the NASA Glenn 2.2 second drop tower were performed (1) to demonstrate the feasibility of performing propellant experiments using the NASA Glenn microgravity facilities, (2) to develop the operational procedures for safe handing of the energetic materials and disposal of their toxic combustion by-products and (3) to obtain initial measurements of the AP burning rate and flame structure under microgravity conditions. Experiments were conducted on the CH4/AP system previously studied at normal gravity using a modified design of the counterflow burner and a NASA Glenn Pig Rig, i.e., one of the existing drop rigs for general-purpose usage. In these experiments, the AP burning rate was measured directly with a linear variable differential transducer (LVDT) and video imaging of the flame structure was recorded ignition was achieved by hot wires stretched across the AP surfaces. Initial drop tower combustion data show that with the same burner separation distance and flow conditions of the normal gravity experiments, the AP burning rate is approximately a factor of two lower. This difference is likely a result of radiation effects, but further tests with longer test times need to be conducted to verify that steady state conditions were achieved under microgravity conditions.

Takahashi, Fumiaki↗
NASA NTRS2012

AdapChem

AdapChem software enables high efficiency, low computational cost, and enhanced accuracy on computational fluid dynamics (CFD) numerical simulations used for combustion studies. The software dynamically allocates smaller, reduced chemical models instead of the larger, full chemistry models to evolve the calculation while ensuring the same accuracy to be obtained for steady-state CFD reacting flow simulations. The software enables detailed chemical kinetic modeling in combustion CFD simulations. AdapChem adapts the reaction mechanism used in the CFD to the local reaction conditions. Instead of a single, comprehensive reaction mechanism throughout the computation, a dynamic distribution of smaller, reduced models is used to capture accurately the chemical kinetics at a fraction of the cost of the traditional single-mechanism approach.

Oluwole, Oluwayemisi O.↗
DOE OSTI2023

Electron heating in kinetic-Alfvén-wave turbulence

We report analytical and numerical investigations of subion-scale turbulence in low-beta plasmas using a rigorous reduced kinetic model. We show that efficient electron heating occurs and is primarily due to Landau damping of kinetic Alfvén waves, as opposed to Ohmic dissipation. This collisionless damping is facilitated by the local weakening of advective nonlinearities and the ensuing unimpeded phase mixing near intermittent current sheets, where free energy concentrates. The linearly damped energy of electromagnetic fluctuations at each scale explains the steepening of their energy spectrum with respect to a fluid model where such damping is excluded (i.e., a model that imposes an isothermal electron closure). In conclusion, the use of a Hermite polynomial representation to express the velocity-space dependence of the electron distribution function enables us to obtain an analytical, lowest-order solution for the Hermite moments of the distribution, which is borne out by numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Ignition delay time and H 2 O measurements during methanol oxidation behind reflected shock waves

To improve detailed chemical kinetics models, the oxidation of methanol was investigated behind reflected shock waves in shock tubes. Ignition delay times of methanol–air mixtures, with Ar as diluent, were studied between 940 and 1540 K in a heated shock tube, for pressures up to 14.9 atm and for equivalence ratios of 0.5, 1.0, and 2.0. Water profiles were measured by utilizing a laser absorption technique in the 1350-to-1600-K temperature range, at an average pressure of 1.3 atm and for similar equivalence ratios. The present study shows here the ignition delay times of methanol to be in very good agreement with results from the literature (Fieweger et al., 1997), whereas the other conditions have never been investigated before. The ignition delay time data are also in good agreement with modern detailed kinetics mechanisms such as the AramcoMech 3.0 model. The water time-history profiles were modeled using well-known literature mechanisms. Discrepancies were observed between these kinetics mechanisms, and poor predictions were observed for the lower temperatures investigated. Sensitivity and rate-of-production analyses were performed using 3 literature mechanisms (namely, AramcoMech 3.0, Princeton, and JetSurfII). Discrepancies were found among the models when predicting important reactions dominating the oxidation of methanol as well as the rate-of-production of H 2 O.

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

Experimental and computational investigation of the influence of iso -butanol on autoignition of n -decane and n -heptane in non-premixed flows

Experimental and computational investigation is carried out to elucidate the influence of iso-butanol on critical conditions of autoignition of n-decane and n-heptane. The counterflow configuration is employed. In this configuration an axisymmetric stream of air is directed over the surface of an evaporating pool of a liquid fuel. A stagnation plane is established. The temperature of the air is increased until autoignition is achieved in the mixing layer in the vicinity of the stagnation plane. The temperature of the air stream at autoignition, T ig , is measured at various values of the strain rate, which is defined as the axial gradient of the axial component of the flow velocity at the stagnation plane. Fuels tested are n-decane, n-heptane, iso-butanol and various mixtures of n-decane/iso-butanol and n-heptane/iso-butanol. Kinetic modeling is carried out using the recently updated version of the comprehensive CRECK chemical-kinetic mechanism. Critical conditions of autoignition are predicted for all fuels and fuel mixtures and the results are compared with the measurements. Computations show that low-temperature chemistry plays a significant role in promoting autoignition of n-decane and n-heptane, and the influence of low-temperature chemistry decreases with increasing strain rates because there is insufficient residence for the low-temperature reactions to take place. Experimental data and numerical simulations show that addition of even small amounts of iso-butanol to n-decane or n-heptane increases the value of T ig at low strain rates indicating that iso-butanol strongly inhibits the low-temperature chemistry of n-decane and n-heptane. Furthermore, the simulations show that when iso-butanol is present in the liquid fuel, only a small amount of n-decane is available in the gas phase and its low concentration cannot effectively sustain the low temperature degenerate branching oxidation process. Here, predicted flame structures show that the peak values of mole fraction of ketohydroperoxide are significantly reduced when iso-butanol is added to n-decane indicating that the kinetic pathway to low temperature ignition is blocked. This observation is confirmed from sensitivity analysis.

42 ENGINEERING↗
DOE OSTI2023

RAM-SCB Runs on Request at CCMC

The Ring current–Atmosphere interactions Model (RAM) with Self-Consistent magnetic (B) field (SCB) combines a large-scale kinetic model of ring current plasma with a three-dimensional (3-D) force-balanced model of the terrestrial magnetic field. RAM-SCB simulates the evolution of major ion species (H + , O + , and He + by default) and electrons as a function of azimuth, radial distance, energy, and pitch angle. Simulation outputs include the Dst index, and pressure and differential flux of the modeled species, thus providing benefit as a science code and to inform surface charging hazard within the model domain. Version 2.2 of this open-source simulation code has now—as of 8 August 2023—been made available for Runs-On-Request via the Community Coordinated Modeling Center.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Enhanced reactivity at the oil–water interface accelerates the synthesis of zymonic acid in microemulsions

Chemical reactions in microscale compartments, such as aerosols and emulsions, can exhibit significantly faster reaction rates relative to macroscale containers. This enhancement in chemistry is often due to the elevated importance of surfaces as reaction vessels are reduced to picoliter volumes. While most studies have focused on the air-water interface of droplets, there are comparably fewer studies of reactions in micron-scale aqueous solutions encapsulated by oil. Here we investigate the condensation reaction of pyruvic acid (PA) to form zymonic acid (ZA) and water. Using microfluidics and optical trapping, chemical kinetics are measured in monodisperse micron-sized emulsion droplets in situ via Raman spectroscopy. Relative to a macroscopic bulk solution, which exhibits little to no reaction over many days, we find efficient production of ZA over the same time period. A kinetic model is developed to elucidate the role of the interface in accelerating the microdroplet reaction kinetics. After quantifying the surface partitioning of PA from interfacial tension measurements, the rate coefficient for the condensation reaction at the oil-water interface is determined to be 1.8 × 10-2 M-1 s-1. This rate coefficient is estimated to be 105 larger than the reaction rate in bulk aqueous solutions. Compared to previous studies of accelerated ZA formation at the air-water interface on nanodroplets, we find that the reaction at the oil-water interface is 20 times more efficient. Despite this difference, the overall ZA formation rate in emulsions is significantly slower than in the same-sized aerosols, which arises from the weaker partitioning of PA to the oil-aqueous relative to air-water interface. These results highlight the interplay between interfacial partitioning and reactivity in accelerating chemistry in microcompartments and provides new insights into how interfacial composition influences condensation reactions.

Wallace, Brandon J↗
DOE OSTI2021

Modeling of Glycolate Destruction in the Recycle Collection Tank

The Savannah River Site’s DWPF is being upgraded with the introduction of the NG flowsheet. Glycolic acid has been shown to be a superior alternative to formic acid for sludge processing. The new flowsheet improves or maintains necessary parameters such as 1) reduction of mercury, 2) adjustment of feed rheology, 3) pH stability, and 4) adjustment of melter oxidation/reduction potential. Further, the use of glycolic acid virtually eliminates the potential for catalytic hydrogen generation in DWPF processing DWPF process condensates are collected and returned to the Savannah River Site (SRS) CSTF. The RCT collects off-gas condensate during chemical processing, vitrification, and other unit operations performed in DWPF and is the singular return vessel delivering recycle effluent back to CSTF. Each batch of recycle will have a small amount of glycolate from chemical processing and melter off-gas condensates. To avoid potential flammability issues due to thermolysis of glycolate in the CSTF, Savannah River National Laboratory (SRNL) provided to Savannah River Remediation (SRR) at their request a Task Technical and Quality Assurance Plan (TTQAP) to quantify and mitigate glycolate returns via DWPF’s recycle stream. The request included testing of a process to oxidize glycolate and other organic species that are responsible for hydrogen generation from thermolysis. Following that work SRR provided a Task Technical Request (TTR) that requested process modeling. In 2021 a TTQAP was issued to cover the modeling work. Modeling draws data from laboratory scale studies using chemical simulants and radioactive waste samples. Chemical kinetic modeling was performed to evaluate the feasibility of using sodium permanganate to destroy glycolate in the RCT. The results from the laboratory studies were summarized in a series of reports. Reference 9 is a report of lab scale processing of actual DWPF Slurry Mix Evaporate Condensate Tank (SMECT) and Offgas Condensate Tank (OGCT) samples in the SRNL Shielded Cells. Tests at caustic conditions demonstrated sodium permanganate was effective in converting glycolate to oxalate, and permanganate (Mn 7+ ) is reduced to manganate (Mn 6+ ) with no significant formation of carbon dioxide or carbonate. Equation (1) was found to best describe the observed reaction of glycolate with permanganate under nominal (60 to 145 mg/L in RCT) glycolate entrainment conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1987

Numerical calculations of turbulent reacting flow in a gas-turbine combustor

A numerical study for confined, axisymmetrical, turbulent diffusion flames is presented. Local mean gas properties are predicted by solving the appropriate conservation equations in the finite difference form with the corresponding boundary conditions. The k-epsilon two-equation turbulence model is employed to describe the turbulent nature of the flow. A two-step kinetic model is assumed to govern the reaction mechanism. The finite reaction rate is the smaller of an Arrhenius type of reaction rate and a modified version of eddy-breakup model. Reasonable agreement is observed between calculations and measurements, but to obtain better agreement, more work is needed on improvements of the above mathematical models. However, the present numerical study offers an improvement in the analysis and design of the gas turbine combustors.

Lin, Chin-Shun↗
DOE OSTI2019

Impact of Ethanol additions on Autoignition Characteristics of a Full Boiling Range Gasoline and its Surrogates at Advanced Engine Conditions

Ethanol is widely blended into commercial gasolines, as an oxygenate to mitigate pollutant formation, as well as an octane enhancer to improve knock resistance. Ethanol can also facilitate the use of lower-grade, or less-refined gasoline in modern internal combustion (IC) engines. However, the chemical kinetic interactions of ethanol with petroleum-based gasoline are not well understood, especially towards autoignition propensity over a range of conditions. These interactions are critical for improving knock resistance, as well as towards performance in low temperature combustion schemes, such as gasoline compression ignition. This study investigates the chemical kinetic interactions of ethanol with a full boiling range gasoline, where blending ratios of 0 to 30 vol./vol. are covered. A mid-octane range (Anti Knock Index 91.5) Fuels for Advanced Combustion Engines gasoline is employed, while a variety of multi-component surrogate blends are also studied. Autoignition data are acquired within a rapid compression machine at temperature, pressure, and fuel loading conditions representative of advanced IC engines. Changes in low-temperature behavior, including first-stage ignition times and extents of low temperature heat release, as well as main ignition times are quantified for the neat fuels and the ethanol-blended fuels. A recently-updated, detailed chemical kinetic model for the gasoline surrogates is used to interpret the measurements, and better understand the kinetic pathways.

autoignition↗
DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
DOE OSTI2026

Drift-kinetic effects of tungsten on plasma response to RMP in ITER

Here, effects of high- Z ( Z is the particle charge number) tungsten impurity ions on the plasma response to the resonant magnetic perturbation (RMP) field are numerically investigated for the ITER 15 MA baseline scenario, where the tungsten contribution to the plasma response is computed with a drift-kinetic model while the bulk thermal particle contributions follow the fluid approximation. The study yields three highlights: (i) the drift-kinetic contribution of the tungsten impurity exerts minor influence on the plasma response compared to that computed by the pure fluid model without tungsten; (ii) a new figure of merit, based on the resonant spectrum perturbation at the plasma boundary surface, results in different optimal coil phasing compared to that previously obtained by maximizing the edge-peeling plasma response; (iii) the optimal $n = 3$ RMP (for edge localized mode (ELM) control, $n$ is the toroidal mode number) is found to induce a large tungsten particle influx near the plasma edge associated with the neoclassical toroidal viscosity. The study thus provides useful data on the compatibility of the full tungsten wall with RMP ELM control in ITER.

ITER↗
DOE OSTI2026

Rate constants and third-body collision efficiencies for recombination of Na with OH and O2: Implications for flame inhibition by alkali salts

Flame inhibition by alkali metals has implications for solid fuel combustion, fire safety, and a number of industrial processes. While the mechanism of inhibition is fairly well understood, details related to thermodynamic properties and rate constants are still in question. In the present work, the recombination of Na with OH (R5) and O2 (R7), respectively, was characterized theoretically, and the implications for modeling laminar premixed hydrogen flames doped with sodium species were examined. Third-body collision efficiencies and low-pressure-limit rate constants were obtained using newly fitted ab initio-based potential energy surfaces, classical trajectories, and one-dimensional master equation calculations. The results are consistent with available experimental results and aid in the present modeling study by providing rate information for bath gases and conditions that remain unexplored experimentally. Chemical kinetic modeling of relative Na and absolute H and OH profiles in H2-fueled laminar, premixed flames doped with a sodium salt shows that the most important radical removal cycle is the sequence Na + OH (+M) → NaOH (+M) (R5), NaOH + H → Na + H2O (R12), even under oxidizing conditions. A secondary cycle, Na + O2 (+M) ⇄ NaO2 (+M) (R7), NaO2 + OH → NaOH + O2 (R14), is less important due to the low thermal stability of NaO2. In most flames, reaction R7 is partially equilibrated, and reaction R14 becomes rate-limiting for the second cycle. The flame analysis supports a lower value of k14 than indicated by recent work on KO2 + OH, but more work is required to confirm this.

Jasper, Ahren W.↗
DOE OSTI2021

Stochasticity and positive feedback enable enzyme kinetics at the membrane to sense reaction size

Here, we present detailed kinetic analyses of a panel of soluble lipid kinases and phosphatases, as well as Ras activating proteins, acting on their respective membrane surface substrates. The results reveal that the mean catalytic rate of such interfacial enzymes can exhibit a strong dependence on the size of the reaction system—in this case membrane area. Experimental measurements and kinetic modeling reveal how stochastic effects stemming from low molecular copy numbers of the enzymes alter reaction kinetics based on mechanistic characteristics of the enzyme, such as positive feedback. For the competitive enzymatic cycles studied here, the final product—consisting of a specific lipid composition or Ras activity state—depends on the size of the reaction system. Furthermore, we demonstrate how these reaction size dependencies can be controlled by engineering feedback mechanisms into the enzymes.

59 BASIC BIOLOGICAL SCIENCES↗