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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Accelerated Steam Methane Reforming by Dynamically Applied Charges

Catalyst design has traditionally focused on tuning active site properties to optimally bind reaction intermediates and balance the kinetic requirements of multiple competing chemical processes, as necessitated by the Sabatier principle. It has recently been proposed that for reactions following certain potential energy landscapes, the activity limit imposed by the Sabatier principle may be overcome by using programmed oscillations of surface electron density at the timescales of surface reactions (i.e., “catalytic resonance”). Here, we use a combination of density functional theory (DFT) simulations and transient kinetic models (TKMs) to simulate the kinetics of steam methane reforming (SMR) on Ru(211) surfaces under statically and dynamically applied charges. DFT-calculated binding energies of SMR intermediates and transition states exhibit strong sensitivity to positively applied charges and follow unique scaling relationships that deviate from linear periodic trends across transition metals. Our simulations demonstrate that applying a small positive charge to Ru dramatically enhances the steady-state turnover frequency (TOF) of SMR by up to 5 orders of magnitude above the TOF observed over neutral Ru. Thus, statically charging Ru catalysts may be an effective strategy to lower the temperature requirements for SMR. Dynamic square-wave oscillations in charge resulted in SMR catalytic resonance with an onset frequency f ∼ 106 Hz and the corresponding average TOFs exceeding the statically charged Ru surface by an additional 15%. Here, based on sensitivity analyses performed for the two end points of oscillation, we propose that dynamic TOF improvement beyond the Sabatier maximum can be expected when the system oscillates between two kinetic regimes that are uniquely controlled by distinct elementary steps.

Catalysts↗
DOE OSTI2019

Analysis and Simulation of One-Dimensional Transport Models for Lithium Symmetric Cells

The main objective of this paper is to analyze transport models for lithium symmetric cells and arrive at an efficient model and code to simulate the same. Two one-dimensional models are considered. The first model uses the dilute solution Nernst–Planck equations in conjunction with the electroneutrality assumption (EN-NP model). For binary electrolytes, an analytical solution for electrolyte concentration is derived and compared with numerical solutions by an approximate finite volume method. The second approach relaxes the electroneutrality assumption by way of Poisson’s equation for the electrostatic potential (PNP model). The computational difficulties of the PNP model are tackled using the approximate finite volume method, with demonstrated convergence characteristics even with bulk dimensions several orders of magnitude larger than the characteristic double layer size. A robust code is developed for the PNP model. The computationally efficient transport models can facilitate simulations, physical understanding, and analysis. This is illustrated by a case study in which these models are coupled with modified kinetic models, which are then parameterized with experimental voltage response data using a systems-level approach. The estimated parameters provide further insight into the electrochemical phenomena underpinning voltage transitions in symmetric cells.

25 ENERGY STORAGE↗
DOE OSTI2022

Integrated Computational Materials and Mechanical Modeling for Additive Manufacturing of Alloys with Graded Structure Used in Fossil Fuel Power Plants

Wire-arc additive manufacturing (WAAM) has demonstrated its unique capability of producing large-size alloy components with a significantly reduced fabrication time and enhanced geometry design freedom. In this project, the team has developed an ICME (Integrated Computational Materials Engineering) modeling framework, which supports the WAAM of the AUSC (Advanced Ultra-Supercritical) power plant components. The manufacturing design has been applied to Inconel 740H, steel P91, as well as the dissimilar alloy components between steel P91 and Inconel 740H. The ICME model framework is developed by considering two types of modeling. First, mechanistic modeling has been applied to control the printing quality and understand the sequence of the dissimilar printing of the wall structure. The following models have been included in the developed ICME framework: finite element thermal model, grain structure model, residual stress simulation, crystal plasticity model, CALPHAD-based precipitation kinetic model, phase stability prediction, thermal expansion predictive model, and heuristic creep model. Secondary, a physics-based machine learning model has also been developed based on the ICME model structure. The machine learning model development is based on the ICME model prediction with calibration of the experiments. In addition, the WAAM has been utilized as a high-throughput experimental tool rapidly generating a gradient of alloy composition to facilitate experimental database generation for process-structure-property relationships. Such a database directly supported the ICME-enhanced machine learning, which further assisted in intermediate composition block design between P91 and 740H. A high-throughput screening study of the oxidation resistance has been performed based on such high-throughput experimentation. Based on the computational design, several dissimilar alloy manufacturing with post-heat treatment have been performed with a comprehensive evaluation of mechanical performance, including hardness mapping, yield strength, creep resistance. In this project, the single component of P91 and 740H processed by WAAM after heat treatment designed by ICME has demonstrated higher performance in yield strength and creep resistance than the wrought materials. The P91 sample prepared by WAAM with ICME-designed heat treatment performs better than P92 in creep resistance. The designed graded alloy printing with intermediate block shows a promising performance that exceeds the traditional welding. Moreover, the current research indicates the high need for location-specific design analysis with uncertainty quantification, an important topic that deserves more dedicated research. The achievement of this project demonstrated the promising future of WAAM in structural alloy manufacturing for energy power plant development. Successful printing requires synergetic efforts made by manufacturing, mechanical, and materials sciences.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1999

SiC and Si3N4 Recession Due to SiO2 Scale Volatility Under Combustor Conditions

Silicon carbide (SiC) and Si3N4 materials were tested in various turbine engine combustion environments chosen to represent either conventional fuel-lean or fuel-rich mixtures proposed for high-speed aircraft. Representative chemical vapor-deposited (CVD), sintered, and composite materials were evaluated by furnace and high-pressure burner rig exposures. Although protective SiO2 scales formed in all cases, the evidence presented supports a model based on paralinear growth kinetics (i.e., parabolic growth moderated simultaneously by linear volatilization). The volatility rate is dependent on temperature, moisture content, system pressure, and gas velocity. The burner tests were thus used to map SiO2 volatility (and SiC recession) over a range of temperatures, pressures, and velocities. The functional dependency of material recession (volatility) that emerged followed the form A[exp(-Q / RT)](P(sup x)v(sup y). These empirical relations were compared with rates predicted from the thermodynamics of volatile SiO and SiOxHy reaction products and a kinetic model of diffusion through a moving boundary layer. For typical combustion conditions, recession of 0.2 to 2 micrometers/hr is predicted at 1200 to 1400 C, far in excess of acceptable long-term limits.

Smialek, James L.↗
DOE OSTI2020

Biological evaluation of Keggin‐type polyoxometalates on tyrosinase: Kinetics and molecular modeling

Abstract Abnormal overexpression of tyrosinase activity can lead to the production of hyperpigmentation in human skin and enzymatic browning in fruits and vegetables. Herein, the inhibition and mechanism of the H 3 PMo 12 O 40 and two transition metal‐substituted Keggin‐type polyoxometalates (Na 7 PMo 11 CoO 40 and Na 7 PMo 11 ZnO 40 ) on tyrosinase were studied by kinetics and molecular modeling. Kinetic studies indicated that all compounds had more potent inhibitory activities than standard arbutin, and H 3 PMo 12 O 40 (IC 50 = 0.443 ± 0.006 m m ) is ~15‐fold stronger inhibition than arbutin. Additionally, all compounds inhibited tyrosinase in a reversible competitive manner. Intriguingly, molecular modeling elucidated that three compounds competitively bind to tyrosinase mainly through more interactions with Cu 2+ ions and the amino acid residue capable of forming cation‐π and hydrogen bonding, forming a reversible non‐covalent complex. Molecular simulation study correlated well with the biological activity of three compounds in vitro. This work provided new insights into design and synthesis of polyoxometalates as tyrosinase inhibitors in the field of medicine, cosmetic, and food.

Chi, Guoxiang↗
DOE OSTI2021

Measurements of Intermediate Species in Fuel-Rich Oxidation of Ethylene, Toluene, and n -Decane

Polycyclic aromatic hydrocarbons (PAHs) are important precursors to the formation of soot and are also pollutant emissions from combustion devices, including internal combustion engines. To understand the formation of PAHs, the low-molecular-weight hydrocarbon chemistry leading to their formation needs to be understood. Toward this goal, fuel-rich oxidation of three different hydrocarbons, i.e., ethylene, toluene, and n-decane, has been investigated in an atmospheric pressure flow reactor at varying temperatures (1000–1350 K), equivalence ratios (φ = 3.0–12.0), and residence times of 0.25–1.5 s. The major C 1 to C 7 intermediates, such as methane, acetylene, ethylene, allene, propyne, propylene, vinylacetylene, 1,3-butadiene, cyclopentadiene, benzene, and toluene, were quantified using a gas chromatograph equipped with a flame ionization detector. The experimental tendency of intermediate species formation, e.g., the dependence of temperature, equivalence ratio, and residence time, was similar in ethylene oxidation and n-decane oxidation. The concentrations of intermediate species up to C 4 were higher in ethylene/n-decane oxidation than toluene oxidation, while a nearly equal or larger amount of C5–C7 species was produced in toluene oxidation. The experimental data were compared with modeling results using a detailed chemical kinetic mechanism. The calculated data using the kinetic model were in agreement with the experimental results. Additionally, a comprehensive kinetic analysis on the reaction pathways of each species was conducted to assess the differences in the oxidation chemistry with the change in the structure of hydrocarbons. In particular, in-depth analysis of benzene formation was performed, elucidating that two benzene formation pathways were important in ethylene and n-decane oxidation: (1) H elimination of the 2,4-cyclohexadienyl radical produced from the isomerization of the 2,4-cyclopentadienylmethyl radical, and (2) dehydrogenation of cyclohexadiene produced from reactions of vinyl radical with 1,3-butadiene. In toluene oxidation, it was found that benzene was primarily produced from toluene through the replacement of the methyl group with hydrogen.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

A purely kinetic description of the evaporation of water droplets

The process of water evaporation, although deeply studied, does not enjoy a kinetic description that captures known physics and can be integrated with other detailed processes such as drying of catalytic membranes embedded in vapor-fed devices and chemical reactions in aerosol whose volumes are changing dynamically. In this work, we present a simple, three-step kinetic model for water evaporation that is based on theory and validated by using well-established thermodynamic models of droplet size as a function of time, temperature, and relative humidity as well as data from time-resolved measurements of evaporating droplet size. Additionally, the kinetic mechanism for evaporation is a combination of two limiting processes occurring in the highly dynamic liquid-vapor interfacial region: direct first order desorption of a single water molecule and desorption resulting from a local fluctuation, described using third order kinetics. The model reproduces data over a range of relative humidities and temperatures only if the interface that separates bulk water from gas phase water has a finite width, consistent with previous experimental and theoretical studies. The influence of droplet cooling during rapid evaporation on the kinetics is discussed; discrepancies between the various models point to the need for additional experimental data to identify their origin.

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

Dehydroxylation kinetics of kaolinite and montmorillonite examined using isoconversional methods

The use of calcined clays as supplementary cementitious materials (SCMs) in concrete is a promising strategy towards decarbonizing the cement and concrete industry. This is especially relevant considering the ever-increasing demand for concrete. Comprehensive understanding of the kinetics of calcination is essential towards maximizing the potential reactivity of clay minerals while ensuring energy efficiency. In this study, the kinetics of the dehydroxylation of kaolinite and montmorillonite are investigated under non-isothermal conditions at constant heating rate. Activation energies ( E a ) are determined via Friedman differential and advanced Vyazovkin incremental methods over the isoconversional range; these are devoid of computational approximations, thus allowing kinetic analysis without assuming a specific reaction model. Kinetic equations—in the differential form as well as a combination of differential and integral forms are compared against the experimentally determined reaction models to identify the most probable dehydroxylation mechanism for kaolinite and montmorillonite. A reaction order mechanism is established for dehydroxylation of kaolinite, while montmorillonite is noted to undergo dehydroxylation via a single-step reversible diffusion-controlled process. Kinetic triplet—comprising activation energy, reaction model and pre-exponential factor—is used to predict isothermal calcination conditions, which is further verified using analytical techniques. Heat release rates of clay-portlandite blends from isothermal calorimetry are used within a thermodynamic framework to quantify reactivity of the calcined clays. Here, the study demonstrates a general approach based on isoconversional methods to predict calcination conditions for different clays that can be used in efficient and optimized production of blended cements or SCMs.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Performance of Polymer Electrolyte Fuel Cell Electrodes with Atomically Dispersed (AD) Fe-C-N ORR Catalyst

Activity of (AD)Fe-N-C catalyst with low Fe content is investigated in differential cells prepared by hot pressing anode gas diffusion electrodes (0.2 mgPt/cm2) to N211 membrane and brush painting the cathode catalyst ink. Polarization curves obtained in H2/O2 show double Tafel slopes which, in conjunction with the redox potential observed in cyclic voltammetry traces, forms the basis for a distributed ORR (oxygen reduction reaction) kinetic model with potential-dependent available sites. Application of this model to polarization data in H2/air provides the basis for formulating an oxygen transport model and leads to 7.7 ± 0.6 mA/cm2 catalyst activity at 0.9 V, 31.5-34.3 mA/cm2 cell performance at 0.8 V, and 0.8-2 s/cm O2 transport resistance. A coupled kinetic, O2 transport and proton transport illustrates the projected improvements needed in catalyst activity and electrode structure to approach the automotive target of 1000 mW/cm2 stack power density while meeting the 1.45 kW/oC heat rejection constraint at 1.5 cathode stoichiometry. The model indicates that we need twelve-fold higher mass activity for reducing the kinetic losses, doubling of active site density and an engineered electrode structure for 50% lower proton transport resistance, as well as a 50% reduction in electrode thickness to limit the O2 transport losses

ACTIVE-SITES↗
DOE OSTI2024

Flame kinetics at scramjet-engine-relevant conditions: Role of prompt dissociation of weakly-bound radicals

Combustion in high-speed ram-based propulsion engines occurs under distinct thermodynamic conditions of high reactant temperatures (greater than 1000 K) and relatively low pressures (<5 atm). There is a lack of fundamental flame measurements at such conditions that result in adiabatic flame temperatures (T ad ) exceeding 2500 K. In this work, we have measured laminar flame speeds of oxygen-enriched CH 4 /oxidizer mixtures at sub-atmospheric conditions to probe kinetics at high T ad using the isobaric spherically expanding flame approach. Simulations with recent kinetic models revealed increasing differences between data and model predictions with increasing T ad , reaching up to 25 %. Kinetic analyses reveal that at the thermodynamic conditions in these O 2 -enriched flames, i.e., lower pressures and higher T ad , the effects of HCO prompt dissociation are accentuated. In addition to HCO, the prompt dissociations of CH 2 OH and C 2 H 5 are also considered. Here, the prompt dissociations of all three radicals were evaluated and their effects considered in flame speed simulations. Reaction path analysis for the present flames revealed that approximately half of the reaction flux for HCO formation undergoes prompt dissociation to H + CO. Furthermore, these analyses also revealed that the pathways and sensitive reactions are similar between oxygen-enriched fuel/oxidizer mixtures and preheated fuel/air mixtures, if both have similar T ad . Thus, flames of oxygen-enriched mixtures could be a surrogate to probe the flame chemistry of highly preheated mixtures at relatively low pressures that are often encountered in ram-based propulsion engine combustors.

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

An initiation-promotion model of tumour prevalence from high-charge and energy radiations

A repair/misrepair kinetic model for multiple radiation-induced lesions (mutation inactivation) is coupled to a two-mutation model of initiation-promotion in tissue to provide a parametric description of tumour prevalence in the mouse Harderian gland from high-energy and charge radiations. Track-structure effects are considered using an action-cross section model. Dose-response curves are described for gamma rays and relativistic ions, and good agreement with experiment is found. The effects of nuclear fragmentation are also considered for high-energy proton and alpha-particle exposures. The model described provides a parametric description of age-dependent cancer induction for a wide range of radiation fields. Radiosensitivity parameters found in the model for an initiation mutation (sigma 0 = 7.6 x 10(-10) cm2 and D0 = 148.0 Gy) are somewhat different than previously observed for neoplastic transformation of C3H10T1/2 cell cultures (sigma 0 = 0.7 x 10(-10) cm2 and D0 = 117.0 Gy). We consider the two hypotheses that radiation acts solely as an initiator or as both initiator and promoter and make model calculations for fractionation exposures from gamma rays and relativistic Fe ions. For fractionated Fe exposures, an inverse-dose-rate effect is provided by a promotion hypothesis with an increase of 30% or more, dependent on the dose level and fractionation schedule, using a mutation rate for promotion similar to that of single-gene mutations.

NASA Discipline Radiation Health↗
DOE OSTI2022

Experimental and modeling study of water time histories during H 2 S-N 2 O combustion in a shock tube

Here, with the interest in directly burning sour gas in gas turbines, and the fact that even small amounts of H 2 S or its combustion products can alter combustion characteristics, many research studies have been performed to better understand the combustion chemistry of H 2 S. In the present study, the water formation was followed by laser absorption with N 2 O as an oxidant, instead of O 2 . Nitrous Oxide being essentially consumed via N 2 O (+M) ⇌ N 2 + O (+M), the water formation via the H 2 S + O route can then be probed to further validate the models. Three H 2 S/N 2 O mixtures diluted in 98% Ar were studied to cover the following range of equivalence ratios: 0.5, 1.0, and 2.0, over a wide range of temperatures (1580–1940 K) around atmospheric pressure. A chemical kinetic model was then developed first by validating the base N 2 O kinetics mechanism, then by investigating the experimental results presented herein. The N 2 O kinetics mechanism was updated with recent work on the low- and high-pressure limits for N 2 O decomposition (N 2 O (+M) ⇌ N 2 + O (+M)) as well as a review of rate coefficients for N 2 O + H ⇌ N 2 + OH from the literature. Good agreement is shown for H 2 /N 2 O mixtures. Updates were then made to the H 2 S kinetics mechanism, specifically, an update from the literature on SO 2 + H (+M) ⇌ SO + OH (+M) and an adjustment to SO + SH ⇌ S 2 + OH. Additionally, reactions between SH and N 2 O were determined using W1BD data and transition state theory which required the addition of an NNS sub-mechanism. With these updates, the mechanism provides good agreement with the H 2 S/N 2 O experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

A theoretical explanation of the effect of oxygen poisoning on industrial Haber-Bosch catalysts

The Haber-Bosch process has been investigated extensively, however a low-temperature, low-pressure process remains elusive. As has been demonstrated many times, this stems in part from the difficulty of breaking the N-N triple bond. In this work, we highlight an additional reason for the lack of a low-temperature ammonia synthesis process: the effect of oxygen poisoning at low temperature. Using density functional theory (DFT), we have created a new model for the active site of industrial Haber-Bosch catalysts which explicitly includes the potassium promoter. Furthermore, we present a new micro-kinetic model for ammonia synthesis that includes the effect of oxygen poisoning due to trace water content in the input gas stream. Our model agrees well with previous experiments and shows that devising a strategy to avoid oxygen poisoning is crucial to creating a low-temperature Haber-Bosch process. Additionally, the model suggests that using a weaker-binding catalyst is one way to avoid oxygen poisoning if it is impractical to remove all water from the reactor.

08 HYDROGEN↗
DOE OSTI2019

Metabolic model guided strain design of cyanobacteria

Cyanobacteria are oxygenic photoautotrophs that serve as potential platforms for the production of biochemicals from cheap and renewable raw materials – sunlight, water, and carbon dioxide. Systems level analysis of the metabolic network of these organisms could enable the successful engineering of these organisms for the enhanced production of target chemicals. Metabolic modeling techniques including both stoichiometric and kinetic modeling with a genome-wide coverage enable a global assessment of metabolic capabilities. Recent studies guided by such modeling techniques have engineered strains for the enhanced production of valuable chemicals such as ethanol, n-butanol, 1,3-propanediol, glycerol, limonene, and isoprene from CO 2 .

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Ultra-low concentration terbium (Tb) adsorption on garlic peels biosorbent and its application for $\mathrm{Nd}$-$\mathrm{Fe}$-$\mathrm{B}$ scraps recovery

Motivated by the strategic value of middle-heavy rare earth elements (MHREEs), we proposed a bio-adsorption method to recover ultra-low concentration terbium (Tb(Ⅲ)) from industrial wastewater that originated from the Nd-Fe-B scraps recovery process. It could be found that Tb(Ⅲ) ions as low as 5.34 ng·mL -1 (ppb) could be adsorbed onto Ca(Ⅱ)-modified garlic peels (Ca-GP) within 10 min with the adsorption efficiency of 99.2% at pH 3.5. The adsorption behavior of Tb(Ⅲ) onto Ca-GP conformed the Langmuir isotherm model and pseudo-second-order kinetic model. Attenuated total reflection Fourier Transform Infrared Spectroscopy (ATR-FTIR), X-ray photoelectron spectroscopy (XPS) and density-functional theory (DFT) calculations results showed that Tb(Ⅲ) ions could be ion-exchanged with cations in Ca-GP such as Ca(Ⅱ) and -COOH. In application, column experiment results showed that the maximum bed adsorption capacity for Tb, praseodymium (Pr), neodymium (Nd), dysprosium (Dy), and Total rare earth elements (REEs) ions calculated by the Thomas model was 0.06, 3.50, 9.65, 7.34, and 27.37 µg·g -1 , respectively, with the initial concentration of 0.37 ng·mL -1 Tb, 20 ng·mL -1 Pr, 44 ng·mL -1 Nd, 67 ng·mL -1 Dy and 170 ng·mL -1 total REEs ions of actual solutions (below the solubility of rare-earth hydroxide). Here this manuscript thus provided a novel cost-effective and efficient approach to the enrichment and recovery of MHREEs from ultra-low concentration REEs ions-containing solutions.

54 ENVIRONMENTAL SCIENCES↗