Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A combined kinetic Monte Carlo and phase field approach to model thermally activated dislocation motion

In this work, a kinetic Monte Carlo (KMC) scheme is integrated into a phase field dislocation dynamics (PFDD) model to account for thermal activation. While the implementation is general and could be applied to many different thermally activated processes, we apply the model to simulate a dislocation transmitting through an interface. Two cases are analyzed: a dislocation transmitting through a strong interface, and a dislocation transmitting through a more realistic grain boundary in tungsten. In conclusion, results and trends for both cases are discussed in detail, and ultimately show that the combined KMC-PFDD algorithm is effective in capturing thermally activated dislocation processes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A kinetic model for simulating non-equilibrium mass transport in oxides applied to hematite growth under irradiation

We propose an approach to simulating the dynamic evolution and transport of charged point defects within and through the oxide scales that form during corrosion. The method follows the cluster dynamics formalism widely adopted for radiation damage in solids, which can apply in both quasi-static and far from equilibrium conditions, such as in radiation environments. By treating each charged defect as a cluster of an atomic defect and associated unit charges, the proposed model flexibly allows charge state transitions and shifts in the Fermi level by absorption and emission of charge carriers from point defects in a reaction network with rates constrained by local equilibrium. Applying this model to hematite predicts changes in self-diffusion and oxidation kinetics in irradiation environments, surprisingly showing reduced oxidation rates in many conditions, despite enhanced self-diffusion. In conclusion, the origin of this effect lies in a change in Fermi level induced by excess vacancies formed under irradiation, which in turn suppresses the transport of cation interstitials which facilitate hematite growth.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Data-driven chemical kinetic reaction mechanism for F-24 jet fuel ignition

A data-driven chemical kinetic mechanism for the military version of Jet A, F-24, is developed for numerical simulations of the ignition process. The main purpose of this study is to obtain a practical F-24 mechanism across wide temperature and equivalence ratio ranges, with a particular focus on the negative temperature coefficient and low temperature regions. The new mechanism (ARLMech-HC-F24) is based on the HyChem model of a similar fuel and optimized using a micro-genetic algorithm against an experimental ignition delay data set of the target fuel. The development and optimization processes include reaction selection, population creation, shuffled tournament implementation based on a merit function, and child-individual creation for the next generation. Several techniques and parameters are proposed to generate an accurate mechanism through an efficient process. The newly introduced data-driven mechanism based on these techniques shows better merit value convergence and represents the ignition behavior more accurately than that without the techniques. This practical mechanism is suitable for the numerical simulations of the F-24 or Jet A ignition problem, and the suggested strategies can be employed in similar problems of rate coefficients estimation.

42 ENGINEERING↗
DOE OSTI2024

A chemical kinetic analysis of knock propensity of methanol-to-gasoline fuel

Production of low carbon gasoline-like fuels such as methanol-to-gasoline (MTG) is a promising approach to achieve rapid greenhouse gas emission reduction of the transportation sector. Despite the fact that gasoline that meets the ASTM D4814 standard for automotive spark-ignition engine fuel can be readily produced from these processes, it is unclear how the composition of MTG may affect engine performance and emissions. Here, in this paper, a surrogate for an MTG is used to numerically study the effects of gasoline composition on knock propensity and on the sensitivity of knock to thermal and fuel stratification, to oxygen dilution and to nitric oxide from exhaust gas recirculation of residual gases. Simulations were performed in ANSYS CHEMKIN-PRO using a comprehensive chemical kinetic mechanism for gasoline surrogates, and results of the MTG surrogate were compared against those of a petroleum-based regular E10 gasoline, termed PACE-20. A premium-grade MTG fuel was also formulated by adding ethanol to the MTG surrogate, and results were compared against those of four premium-grade, gasoline-like fuels representative of future alternative gasoline formulations. Surrogates and mechanism were evaluated by comparison against experimental engine data, and the model showed high accuracy at stoichiometric conditions (mean absolute error of ignition timing equal to 1.46 crank angle degrees) but larger deviations at lean conditions (mean absolute error of ignition timing equal to 5.52 crank angle degrees). Despite the fact that the MTG surrogate has a RON 1.1 units higher than that of PACE-20, it may show higher knock propensity at medium temperature conditions due to a less intense NTC behavior. MTG autoignition was more temperature- and equivalence ratio-sensitive than that of PACE20, suggesting that MTG can benefit more from naturally-occurring thermal stratification or from induced fuel stratification of the end gas to mitigate knock intensity. The sensitivity of autoignition reactivity to oxygen dilution and to NO concentration was higher for MTG than for regular gasoline at medium loads, but the opposite trend was observed at high loads due to the effect of pressure on the low-temperature chemistry of regular gasoline. Approximately 14 % vol ethanol content was required to upgrade the octane rating of MTG from regular grade to premium grade. Adding 13.6 % vol ethanol made the fuel autoignition less sensitive to both oxygen dilution and NO content (ignition time varies approx. 17 % and 50 % less with oxygen dilution and NO addition, respectively, when adding ethanol at high engine loads).

02 PETROLEUM↗
DOE OSTI2021

Mechanisms and kinetics of citrate-promoted dissolution of a uranyl phosphate mineral

The formation of uranyl phosphate precipitates is a remediation strategy to limit the mobility of uranium in contaminated soils. However, exposure to organic ligands, like the plant exudate citrate, can remobilize the uranium. The purpose of this study is to provide a more thorough comprehension of citrate-promoted dissolution of a uranyl phosphate mineral, chernikovite [(H 3 O)(UO 2 )(PO 4 )·3H 2 O], by determining the extent of uranium release from chernikovite at a wide range of citrate concentrations. Here we have quantified the kinetics of dissolution and proposed potential mechanisms of chernikovite dissolution to gain a better understanding of the fate of uranyl phosphate precipitates in the environment. Batch dissolution and continuously stirred tank reactor (CSTR) experiments indicate that increasing citrate concentrations from 0.1 mM to 50 mM increases the concentration of dissolved uranium. However, at citrate concentrations of 10 mM and greater, the effectiveness of the ligand to enhance the dissolution decreases. Density functional theory modeling and Raman spectroscopy indicate a strong interaction between citrate and the uranyl ions at the surface of the uranyl phosphate. The interaction between citrate and uranyl suggests the formation of an alteration layer of uranyl-citrate surface complexes, which could impede uranium release at citrate concentrations of 10 mM or higher. Together, these results indicate that the citrate-promoted dissolution of chernikovite is a fast reaction that is hindered by a combination of surface saturation and a secondary-phase precipitation reaction at higher concentrations of citrate and reaction times greater than or equal to 72 min based on the hydraulic residence time of CSTR experiments. These findings improve our understanding of the stability of uranyl phosphate for use in sequestration of uranium groundwater plumes at contaminated sites.

58 GEOSCIENCES↗
DOE OSTI2023

Mean field model parameterization to recover coverage-dependent kinetics

Lateral interactions between adsorbates introduce coverage dependence into adsorption energies and activation barriers of surface reactions. Lattice-based kinetic Monte Carlo (kMC) simulations can capture these interactions quantitatively but are laborious to parameterize and solve. Mean field models are more tractable, but protocols to construct and parameterize them are unclear. Here we explore the ability of a coverage-aware mean-field model to map to a lattice-kMC model of a generic two-step reaction network, including quasi-equilibrated adsorption and rate-limiting dissociation steps. We derive expressions for mean-field and coverage-dependent adsorption energies and dissociation barriers and parameterize against lattice-kMC predictions. We show that the parameterized mean-field rates correlate with ground truth lattice-kMC results across a wide range of reaction conditions and identify regions where the mean field fails. The mean field model similarly captures kMC-derived rate-order, Arrhenius and Sabatier plots at a greatly reduced computational cost. Further, the results provide guidance for parameterizing mean-field models, benchmarked against explicit lattice-based approaches for incorporating the influence of coverage effects.

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

Mechanistic insights into nitrogen activation on atomic Ru clusters in self-pillared pentasil using operando atomistic models and experimental kinetics

Alternative catalysts to the industrial Haber Bosch process have been of significant interest in the field of heterogeneous catalysis, yet realizing ammonia synthesis under mild conditions (e.g., 300 °C and 10 bar) is challenging due to the low per-pass conversion. One strategy is to promote the associative ammonia synthesis mechanism which eschews direct N-N bond cleavage. Here, in this work, we use self-pillared pentasil, a self-pillared hierarchical zeolite built by thin MFI zeolite nanosheets, as a support for subnanometric Ru clusters to synthesize ammonia. We show that Ru remains well-dispersed during reaction and further demonstrate that ammonia synthesis rates are higher than Cs-Ru/MgO. Reaction kinetics show a positive order in H 2 providing evidence for the associative mechanism, which then becomes negative in H 2 if Ru is allowed to aggregate into nanoparticles. Operando Density Functional Theory models for Ru speciation in SPP, free energy diagrams, and microkinetic modeling were then applied to develop a reaction mechanism that involves sequential hydrogenation of N 2 from metallic Ru clusters. For this hydrogenation to occur, there are site requirements for N 2 to adopt a bridge-bound configuration that facilitates sequential hydrogenation on single sites and metal clusters. These site requirements in turn inform the design of improved zeolite-supported ammonia synthesis catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A slip-spring framework to study relaxation dynamics of entangled wormlike micelles with kinetic Monte Carlo algorithm

Hypothesis: Wormlike micelles (WLMs) formed due to the self-assembly of amphiphiles in aqueous solution have similar viscoelastic properties as polymers. Owing to this similarity, in this work, it is postulated that kinetic Monte Carlo (kMC) sampling of slip-springs dynamics, which is able to model the rheology of polymers, can also be extended to capture the relaxation dynamics of WLMs. Theory: The proposed modeling framework considers the following relaxation mechanisms: reptation, union-scission, and constraint release. Specifically, each of these relaxation mechanisms is simulated as separate kMC events that capture the relaxation dynamics while considering the living nature of WLMs within the slip-spring framework. As a case study, the model is implemented to a system of sodium oleate and sodium chloride to predict the linear rheology and the characteristic relaxation times associated with the individual relaxation mechanisms at different pH and salt concentrations. Findings: Linear rheology predictions were found to be in good agreement with experimental data. Furthermore, the calculated relaxation times highlighted that reptation contributed to a continuous increase in viscosity while union-scission contributed to the decrease in viscosity of WLM solutions at a higher salinity and pH. Finally, this manifests the proposed model’s capability to provide insights into the key processes governing WLM’s rheology.

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

Structure-preserving neural networks for the regularized entropy-based closure of a linear, kinetic, radiative transport equation

The main challenge of large-scale numerical simulation of radiation transport is the high memory and computation time requirements of discretization methods for kinetic equations. In this work, we derive and investigate a neural network-based approximation to the entropy-based closure method to accurately compute the solution of the multi-dimensional moment system with a low memory footprint and competitive computational time. We extend methods developed for the standard entropy-based closure to the regularized entropy-based closures. The main idea is to interpret structure-preserving neural network approximations of the regularized entropy-based closure as a two-stage approximation to the original entropy-based closure. We conduct a numerical analysis of this approximation and investigate optimal parameter choices. Our numerical experiments demonstrate that the method has a much lower memory footprint than traditional methods with competitive computation times and simulation accuracy. The code and all trained networks are provided on GitHub.

entropy closure↗
DOE OSTI2023

Allometric-kinetic model predictions of concentration ratios for anthropogenic radionuclides across animal classes and food selection

Protection of the environment from radiation fundamentally relies on dose assessments for non-human biota. These dose assessments use measured or predicted concentrations of radionuclides in soil or water combined with Concentration Ratios (CRs) to estimate tissue concentrations in animals and plants, yet there is a paucity of CR data relative to the vast number of potential taxa and radioactive contaminants in the environment and their taxon-specific ecosystems. Because there are many taxa each having very different behaviors and biology and many possible bioavailable radionuclides, CRs have the potential to vary by orders-of-magnitude, as often seen in published data. Given the diversity of taxa, the International Commission on Radiological Protection (ICRP) has categorized non-human biota into 12 representative organisms (ROs), while the U.S. Department of Energy (DOE) uses the non-taxon specific categories of terrestrial, riparian, and aquatic animals. The question we examine here is: do these categorizations adequately “represent” the broad diversity of animals? To explore this question, here we analyze an Allometric-Kinetic (A-K) model to calculate radionuclide-specific CRs for common animal classes which are then further subcategorized into herbivores, omnivores, carnivores, and invertebrate detritivores. Comparisons in CRs among animal classes exhibited only small differences, but there was order of magnitude differences between herbivores, carnivores, and especially detritivores, for many radionuclides of interest. These findings suggest that the ICRP ROs and the DOE categories are reasonable, but their accuracy could be improved by including sub-categories related to animal dietary ecology. Finally, comparisons of A-K model predicted CR values to published CRs show order-of-magnitude variations, providing justification for additional studies of assimilation across radionuclides and animal classes.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2023

Synthesis, sintering, and grain growth kinetics of Hf 6 Ta 2 O 17

Here Aasystematic study of the solid-state synthesis, pressureless sintering, and grain growth kinetics of Hf 6 Ta 2 O 17 is presented. The ideal conditions for solids-state synthesis of Hf 6 Ta 2 O 17 powder with minimal particle necking was 1250 °C for 2 h in air. The resultant powder has an average particle size of 210 ± 70 nm. The combined synthesis and ball-milling procedure produces highly sinterable Hf 6 Ta 2 O 17 powder, achieving > 97 % of theoretical density after pressureless sintering at 1600 °C for 2 h in air. The grain growth mechanism was sensitive to processing conditions, appearing to be primarily driven by surface diffusion below 1600 °C and grain boundary diffusion above 1650 °C. The respective activation energies for grain growth were found to be Q S = 659 ± 79 kJ mol -1 and Q GB = 478 ± 63 kJ mol -1 .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Metamagnetism and kinetic arrest in a long-range ferromagnetically ordered multicaloric double perovskite Y 2 CoMnO 6

A systematic magnetic study of the double perovskite oxide Y 2 CoMnO 6 (YCMO) has been performed. A monoclinic P21/n phase of YCMO was synthesized using a sol-gel method. Neutron diffraction (ND) measurements evidence the onset of long-range ferromagnetic (FM) ordering at TC ~ 76 K, which persists down to 5 K. The presence of 25% antisite disorder, estimated from the ND data, leads to the appearance of short-range antiferromagnetic (AFM) interactions. The existence of thermal hysteresis due to competing interactions between FM and AFM phases is observed in the M vs. T measurements. The pinning of magnetic domain walls at the Co/Mn antiphase boundaries results in a metamagnetic-like behavior. The field dependence of thermomagnetic irreversibility and the nature of virgin curves indicate the occurrence of a kinetic arrest phenomenon, which is further verified via cooling and heating of the system in an unequal fields (CHUF) protocol. In full agreement with the ND data, a detailed analysis of critical exponents near the paramagnetic (PM)-FM phase transition also establishes YCMO as a long-range interacting mean-field system. A close examination of the temperature- and field-dependent magnetic entropy change yields an in-depth understanding of coexisting magnetically ordered and disordered phases in YCMO, leading to a comprehensive magnetic phase diagram of this multifunctional double perovskite system.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Conversion kinetics during melting of simulated nuclear waste glass feeds measured by dissolution of silica

In this work, we experimentally investigated the batch-to-glass conversion kinetics of two low-activity and two high-level nuclear waste melter feeds that represent widely different behavior during melting; tests comprised thermogravimetric analysis, the feed expansion test, evolved gas analysis, and x-ray diffraction. Our study indicates that the conversion processes, which include gas-evolving reactions, dissolution of crystalline phases, and primary foam expansion and collapse, are interconnected and influenced by the feed’s thermal history. The fraction of dissolved silica was adopted as a measure of the extent of conversion inside the cold cap—the floating layer of feed material on the top surface of the melt in an electric melter. The interface between the cold cap and the free-flowing melt below was associated with a constant fraction of dissolved silica. The estimated viscosity of glass-forming melt at this interface lies within a narrow range for the feeds studied.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Kinetics and transport of hydrogen in graphite at high temperature and the effects of oxidation, irradiation and isotopics

The kinetics of uptake and desorption impact the performance of graphite as a vector for tritium in high-temperature fission reactors and in the blanket of fusion reactors. Graphite components in these reactors are exposed to temperatures > 500 °C and H 2 partial pressures of few Pa and desorption temperatures are limited to < 1600 °C; limited data is available at these conditions. Here we review the mechanisms for uptake in, transport and desorption of hydrogen from graphite at high temperature, compiling data on uptake rates, diffusion coefficients and activation energies and providing a discussion of the impact of irradiation, pre-oxidation and isotope. At FHR conditions, trapping impacts uptake rates, leading to a reduction in apparent diffusivity by 35 to 80% compared to higher partial-pressure uptake. Timelines for desorption are not clearly defined; extrapolating from available data, at 1150 °C desorbing 80% of tritium uptaken at FHR conditions may take from 100 to 10,000 h.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bulk and surface diffusion of neodymium in alpha-uranium: Ab initio calculations and kinetic Monte Carlo simulations

A fundamental understanding of lanthanide transport in metallic fuels is critical for high fidelity modeling of the fuel-cladding chemical interaction (FCCI) phenomenon, which can lead to the formation of brittle intermetallic compounds and premature failure of the cladding. Here we report a combined ab initio density functional theory (DFT) and kinetic Monte Carlo (KMC) study of the bulk diffusivity of Nd in α-U, fully taking into account the effect of radiation enhanced diffusion. In this work, the vacancy mechanism is considered to be the dominant mechanism for the bulk diffusion of Nd since a Nd interstitial is found to be intrinsically unstable in α-U. The surface diffusivity of a Nd adatom on α-U (001) surface has been further predicted using KMC simulations parameterized by DFT calculations. The present study suggests that Nd transport via the surface diffusion mechanism can be many orders of magnitude faster than bulk diffusion. Furthermore, the results from the present lower length scale study can be used to inform mesoscale phase-field simulations to determine the effective diffusion coefficient of Nd through α-U with a porous microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multi-fidelity kinetic theory-based approach for the prediction of particle attrition: Application to jet cup attrition system

The timescale difference between the “fast” flow dynamics of fluidized bed reactors and the relatively “slow” rate of particle degradation makes the direct computational prediction of attrition challenging. An approach to this challenge is a multi-fidelity strategy where a high fidelity model for the flow dynamics is coupled with a lower fidelity model for the long-time resolution of the bulk attrition of the reactor inventory. We implement this approach using high-fidelity kinetic theory simulations to calculate the flow dynamics which are post-processed to calculate the frequency and intensity of the particle-particle and particle-wall collisions (e.g. collision energy spectra). This is combined with the particle breakage properties to construct the coefficients for a low-fidelity model [e.g. Monazam et al., 2018, Powder Technology 340, p. 528-536]. Simulations are performed of a jet cup attrition system containing Canadian hematite (Monazam et al. 2018). Furthermore, these are first analyzed using the collision energy spectra. Quantitative predictions of the mass loss are made using a low-fidelity model derived from the collision-spectra and a calibrated material breakage coefficient. The results are found to compare favorably with the experimental measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Chemical insights into ethyl acetate flames from experiment and kinetic modeling: Laminar burning velocity, speciation and NO$_x$ emission

Oxygenated fuels, such as alcohols, ethers, and esters, are promising alternatives to conventional fuels. These fuels can help reduce detrimental emissions like carbon monoxide and unburned hydrocarbons and enhance octane ratings. Among these oxygenates, ethyl acetate (EA), a small alkyl ester sourced from biomass, emerges as a clean, promising energy carrier. It serves as a surrogate fuel to facilitate investigations into the combustion behaviours of biodiesel. Despite its importance, the literature knowledge of EA combustion characteristics is limited. Therefore, this study aims to broaden the knowledge of the combustion behaviour of this type of oxygenated fuel compound. In this study, we measured the laminar burning velocities of EA by employing a heat flux burner and a closed combustion vessel over the equivalence ratios of 0.7 – 1.7, pressures of 1 – 10 bar and temperatures ranging from 353 – 423 K. Further, we also measured the NOx emissions in exhaust gas of the premixed flames fueled by EA/air for the first time over the equivalence ratio of 0.8 – 1.2. Additionally, we employed a non-premixed counterflow flame setup for extensive characterisation of species and their concentration under diverse conditions encompassing various strain rates and oxygen concentrations. Finally, we utilized these newly measured data to construct and validate a detailed kinetic model developed as part of this work. The newly developed model will help characterize the combustion properties of EA.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Multi-species electrolyte migration kinetics with directed flow

Here, this paper investigates an electrochemical process designed to control mixing kinetics in liquid electrolytes. The process utilizes an electrochemical cell with anode and cathode chambers connected by a channel without a membrane. This configuration depletes multiple electroactive species in the cathode chamber while selectively replenishing them in the anode chamber. A key application is the removal of impurity species from the cathode chamber electrolyte and the selective replenishment of primary electroactive species in the anode chamber. Aqueous cell experiments were conducted to validate the theoretical model, and simulations were performed for molten salt electrolytes. This approach minimizes waste by reducing the need for additional electrolyte supplies and extends the utilization of anion species, contributing to environmentally sustainable electrochemical materials processing.

Electrochemical Process↗