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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Determination of the characteristic magnetic pre-sheath length at divertor surfaces using micro-engineered targets on DiMES at DIII-D

The magnetic pre-sheath (MPS) width, L MPS , is a critical parameter to define the sheath potential, which controls the ion trajectory of low-Z species (D, T, He, and C), as well as the prompt re-deposition of high-Z species. In this work, to determine L MPS , we fabricated micro-trenches (30×30×4 µm) via focused ion beam (FIB) milling on a silicon surface and exposed them to L-mode deuterium plasmas in DIII-D via the Divertor Material Evaluation System (DiMES) removable sample exposure probe. The areal distribution of impurity depositions, mainly consisting of carbon, was measured by energy-dispersive X-ray spectroscopy (EDS) to reveal the deuterium ion shadowing effect on the trench floors. The carbon deposition profiles showed that the erosion was maximized for the azimuthal direction of φ = -40° (referenced to the toroidal magnetic field direction) as well as the polar angle of θ = 80°. A Monte Carlo equation-of-motion model, based on a collisionless MPS, was used to calculate the azimuthal and polar deuterium ion angle distributions (IADs) for a range of L MPS = k × ρ i , where ρ i is the ion gyro radius and k = 0.5-4. Then, gross erosion profiles were calculated by a Monte Carlo micro-patterning and roughness (MPR) code for ion sputtering using as input the calculated azimuthal and polar IADs for each value of k . Good agreement with the experimental C deposition profiles was obtained for the case k = 2.5-3.5. This result is consistent with a previous kinetic modeling prediction of k ~ 3, as well as previous analytical investigations that predicted the L MPS to be several ion gyro radii. A validation of theoretical sheath models supports its applicability to ITER and pilot plant divertors to successfully predict plasma-materials interactions.

Erosion↗
DOE OSTI2026

Development of a deep potential model for F and CF 2 etching of Si and SiO 2

An understanding of plasma-surface interactions at increasingly smaller scales is invaluable for the development of novel technologies and processing techniques. Molecular dynamics (MD) simulations can provide insights into atomic-scale interactions, though they are restricted by the availability of interatomic potentials. Machine learning methods, such as Deep Potential Molecular Dynamics (DeepMD), provide a systematic framework for the development of accurate and flexible ab initio-based models. In this work, we develop DeepMD models for the ion-enhanced etching of Si and SiO 2 by F and CF 2 radicals. We employ an active learning process to expand the data set on which the model is trained and demonstrate its effect on the model accuracy. The DeepMD results are compared to data from classical MD simulations and experiments. Physical sputtering yields of SiO 2 by Ar + ions show good agreement with previous simulation results using conventional interatomic potentials, though the predicted depth profiles are different. Etching yields are calculated as a function of ion energy and neutral to ion flux ratio for the Ar + ion-enhanced etching of SiO 2 and Si by F atoms, as well as for etching of SiO 2 by CF 2 radicals, showing reasonable agreement with experimental data. Finally, an ion-enhanced surface kinetic model is fitted to the DeepMD etch yields, and the fitted parameters are compared to quantities computed directly from DeepMD simulations. This study illustrates how molecular dynamics simulations using machine learning potentials can provide an accurate model of etching processes relevant to device manufacturing.

Kounis-Melas, Andreas [Princeton Univ., NJ (United↗
DOE OSTI2022

The role of filamentation in activation and DNA sequence specificity of the sequence-specific endonuclease SgrAI

Filament formation by metabolic, biosynthetic, and other enzymes has recently come into focus as a mechanism to fine-tune enzyme activity in the cell. Filamentation is key to the function of SgrAI, a sequence-specific DNA endonuclease that has served as a model system to provide some of the deepest insights into the biophysical characteristics of filamentation and its functional consequences. Structure-function analyses reveal that, in the filamentous state, SgrAI stabilizes an activated enzyme conformation that leads to accelerated DNA cleavage activity and expanded DNA sequence specificity. The latter is thought to be mediated by sequence-specific DNA structure, protein–DNA interactions, and a disorder-to-order transition in the protein, which collectively affect the relative stabilities of the inactive, non-filamentous conformation and the active, filamentous conformation of SgrAI bound to DNA. Full global kinetic modeling of the DNA cleavage pathway reveals a slow, rate-limiting, second-order association rate constant for filament assembly, and simulations of in vivo activity predict that filamentation is superior to non-filamenting mechanisms in ensuring rapid activation and sequestration of SgrAI's DNA cleavage activity on phage DNA and away from the host chromosome. In vivo studies demonstrate the critical requirement for accelerated DNA cleavage by SgrAI in its biological role to safeguard the bacterial host. Collectively, these data have advanced our understanding of how filamentation can regulate enzyme structure and function, while the experimental strategies used for SgrAI can be applied to other enzymatic systems to identify novel functional roles for filamentation.

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

SAM Code Enhancements for Fission Product Tracking of Noble Gases and Metals in MSRs

This report documents fiscal year 2026 enhancements to the System Analysis Module (SAM) for modeling fission product transport in liquid-fueled molten salt reactors (MSRs). The work advances three principal areas: noble gas transport, noble metal deposition, and user interface improvements. The noble gas transport capability integrates drift-flux gas transport, Henry’s law two-film interphase mass transfer with pressure-based nucleation suppression, Knudsen-regime pore diffusion into porous graphite with a conjugate salt-graphite interface constraint, built-in material properties, five Sherwood-number mass transfer correlations including three derived from high-fidelity NekRS simulations, and xenon-135 reactivity feedback through SAM’s point-kinetics model. This work also presents a comprehensive verification test suite, including new analytically verified cases for pressure-dependent onset of interphase gas transfer in a stagnant vertical pipe, a postulated FLiBe-graphite Xe extraction permeator, a gravity riser with a fission-product source, and a descending pipe with gas redissolution driven by hydrostatic pressure. A machine learning framework for bubble rise velocity prediction in molten salt systems is developed and benchmarked on molten-salt and diverse aqueous bubble datasets. The best-performing fine-tuned transfer-learning networks achieve an 82% reduction in RMSE relative to the Clift correlation, and is implemented directly in SAM. The noble metal transport capability is developed, including a liquid-wall deposition model and a gas-surface flotation mechanism that transfers insoluble particles entrained by sparging gas to wetted structures. Verification tests and demonstration cases cover the surface deposition, flotation efflux, and flotation shedding. Finally, a new [SpeciesTransport] input structure replaces positional global vectors with selfcontained, order-independent, named species blocks, simplifies the specification of multiphase species and decay chains, and remains fully compatible with existing SAM input files. Together, these developments improve the physical fidelity, verification basis, and usability of SAM for system-level analyses of fissionproduct behavior in MSRs.

Mui, Travis (ORCID:0000000303736470)↗
DOE OSTI2025

Ab Initio Molecular Dynamics Insights into Stress Corrosion Cracking and Dissolution of Metal Oxides

Oxide phases such as α-Fe2O3 (hematite) and α-Al2O3 (corundum) are highly insoluble in water; however, subcritical crack growth has been observed in humidity nonetheless. Chemically induced bond breaking at the crack tip appears unlikely due to sterically hindered molecular transport. The molecular mechanics of a crack in corundum with a reactive force field reveal minimal lattice trapping, leading to bond breaking before sufficient space opens for water transport. To address this, we model a pre-built blunt crack with space for H2O molecule adsorption at the tip and show that it reduces fracture toughness by lowering the critical J-integral. Then, we explore stress-enhanced dissolution to understand the mechanism of crack tip blunting in the oxide/water system. Density functional theory combined with metadynamics was employed to describe atomic dissolution from flat hematite and corundum surfaces in pure water. Strain accelerates dissolution, stabilizing intermediate states with broken bonds before full atom detachment, while the free energy profile of unstrained surfaces is almost monotonic. The atomistic calculations provided input for a kinetic model, predicting the shape evolution of a blunt crack tip, which displays three distinct regimes: (i) dissolution primarily away from the tip, (ii) enhanced blunting near but not at the apex, and (iii) sharpening near the apex. The transition between regimes occurs at a low strain, highlighting the critical role of water in the subcritical crack growth of oxide scales, with dissolution as the fundamental microscopic mechanism behind this process.

Chemistry↗
DOE OSTI2021

Temperature-responsive smart tracers for field-measurement of inter-well thermal evolution: Heterogeneous kinetics and field demonstration

Temperature-responsive smart tracers enable advanced warning of a “premature thermal breakthrough” and may therefore improve reservoir management and reduce financial risk. A successful calculation of inter-well fluid temperatures requires that temperature-dependent kinetics are known, which typically result from homogeneous batch reactor experiments. However, recent meso-scale field experiments at the Altona Field Laboratory involving fluid flow in a discretely-fractured reservoir suggest that silica-fluid interfaces may accelerate hydrolysis kinetics. Here, the Arrhenius parameters of phenyl acetate hydrolysis are investigated under heterogeneous reaction conditions in packed-bed column experiments. The breakthrough curve of the reaction product is compared to an inert reference tracer (carbon-cored nanoparticles). Temperatures experienced during field testing ranging from 10 to 40 °C were studied using a phosphate buffer solution to maintain a near-neutral pH, as found in the reservoir. The empirically-determined pre-exponential factor and activation energy values are subsequently used in a first-order kinetic model to improve calculations of effective reservoir temperatures for these meso-scale field tests. Here, the results suggest that hydrolysis reactions can directly monitor inter-well reservoir temperatures, but only if the influence of solid-fluid interactions are carefully accounted for.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Single-Molecule Kinetics of Styrene Hydrogenation on Silica-Supported Vanadium: The Role of Disorder for Single-Atom Catalysts

We report a theoretical approach for the study of supported atom catalysis is developed based on recent advances in the study of single-molecule kinetics. This view is particularly useful in exhibiting the role of disorder in single-atom and single-site catalysts on amorphous supports. The distribution of passage times (or waiting times) through a complex catalytic network originating from a set of coupled active sites is described by a probability distribution function (PDF), f(t), that reflects the local environment of the reaction center. An efficient algorithm is developed based on the linear algebra of the Markov transition matrix that produces f(t) or its moments. The kinetics of the hydrogenation reaction of styrene on an organovanadium(III) catalyst supported on amorphous silica is studied. A kinetic model consisting of three intertwined catalytic cycles emanating from three chemically distinct active sites is proposed to describe the chemistry. Density functional theory (DFT) calculations are employed to determine the free energy barriers of the reactions, which are used to construct the rate coefficient matrix. The disorder induced by the amorphous support material is divided into a low-dimensional short-range component reflecting the covalent structures near the reaction center and a weaker long-range component modeling the bulk randomness. The results are computed and analyzed for a wide range of concentration values and disorder scenarios. The unusual structure in the f(t) PDF is found to occur for certain cases that reveal the contribution of multiple catalytic pathways acting in concert.

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

FAGE measurements of tropospheric HO with measurements and model of interferences

Ambient HO measurements by low-pressure laser-excited fluorescence with chemical modulation, and supporting ozone and water-vapor data, are presented for periods in May and August 1987. The observed peak daytime ambient HO concentrations are in the range (2.5 to 8) x 10 exp 6 molecules/cu cm and show small negative offsets due to photochemical interference. Direct measurements of the interference at fixed (O3) give the dependence on ambient (H2O) and on the modulating reagent (isobutane). At ambient (O3) = 30 ppb and 10 torr H2O, with excitation and detection at a total pressure of 4 torr, the net interference is equal to (HO) = -1.3 x 10 exp 6 molecules/cu cm. Production of HO by the reaction of isobutane with O(1D) accounts for the negative interference. Quenching of HO fluorescence by the modulating reagent contributes a smaller positive term to the interference; kinetic measurements of the quenching rate coefficient are reported. The experimental interference results are compared with a detailed kinetic model of HO production, excitation, relaxation, and detection; reasonable agreement is found.

Hard, T. M.↗
DOE OSTI2020

Novel Application of Hybrid Anion Exchange Resin for Phosphate Desorption Kinetics in Soils: Minimizing Re-Adsorption of Desorbed Ions

The process of phosphate desorption from soils is difficult to measure using stirred batch techniques because of the accumulation of desorbed ions in a bathing solution. To accurately measure the apparent rate coefficient of phosphate desorption from soils, it is necessary to remove the desorbed ions. In this study, a novel hybrid (i.e., iron oxide coated) anion exchange resin was used as a sink to study long-term (seven days) P desorption kinetics in intensively managed agricultural soils in the Midwestern U.S. (total phosphorus (TP): 196–419 mg/kg). The phosphate desorption kinetics in the hybrid anion exchange resin method were compared with those in the other conventional batch desorption method with pure anion exchange resins or without any sink. The extent of P desorption in the hybrid resin methods was >50% of total desorbed phosphate in the other methods. The initial kinetic rate estimated in the pseudo-second-order kinetic model was also highest (3.03–31.35 mg/(g·hr)) in the hybrid resin method when the same soil system was compared. This is because adsorbed P in the hybrid resins was nearly irreversible. The hybrid anion exchange resin might be a new and ideal sink in measuring the P desorption process in soils and sediments

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

Plume Characterization of a Laboratory Model 22 N GPIM Thruster via High-Frequency Raman Spectroscopy

The Green Propellant Infusion Mission (GPIM) will demonstrate the capability of a green propulsion system, specifically, one using the monopropellant, AF-M315E. One of the risks identified for GPIM is potential contamination of sensitive areas of the spacecraft from the effluents in the plumes of AF-M315E thrusters. Plume characterization of a laboratory-model 22 N thruster via optical diagnostics was conducted at NASA GRC in a space-simulated environment. A high-frequency pulsed laser was coupled with an electron-multiplied ICCD camera to perform Raman spectroscopy in the near-field, low-pressure plume. The Raman data yielded plume constituents and temperatures over a range of thruster chamber pressures and as a function of thruster (catalyst) operating time. Schlieren images of the near-field plume enabled calculation of plume velocities and revealed general plume structure of the otherwise invisible plume. The measured velocities are compared to those predicted by a two-dimensional, kinetic model. Trends in data and numerical results are presented from catalyst mid-life to end-of-life. The results of this investigation were coupled with the Raman and Schlieren data to provide an anchor for plume impingement analysis presented in a companion paper. The results of both analyses will be used to improve understanding of the nature of AF-M315E plumes and their impacts to GPIM and other future missions.

Raman Spectroscopy↗
NASA NTRS1980

Solar-pumped gas laser development

A survey of gas properties through detailed kinetic models led to the identification of critical gas parameters for use in choosing appropriate gas combinations for solar pumped lasers. Broadband photoabsorption in the visible or near UV range is required to excite large volumes of gas and to insure good solar absorption efficiency. The photoexcitation density is independent of the absorption bandwidth. The state excited must be a metastable state which is not quenched by the parent gas. The emission bandwidth must be less than 10 A to insure lasing threshold over reasonable gain lengths. The system should show a high degree of chemical reversibility and an insensitivity to increasing temperature. Other properties such as good quantum efficiency and kinetic efficiency are also implied. Although photoexcitation of electronic vibrational transitions is considered as a possible system if the emission bands sufficiently narrow, it appears that photodissociation into atomic metastables is more likely to result in a successful solar pumped laser system.

Wilson, J. W.↗
DOE OSTI2025

Pyridine Catalysis of Anhydride Hydrolysis within Carbodiimide‐Driven Reaction Networks

Carbodiimide-fueled reaction networks offer a versatile platform for nonequilibrium chemical systems. Typically, the carbodiimide converts a carboxylic acid to its anhydride, called “activation,” which subsequently undergoes hydrolysis, called “deactivation.” Here, we investigate pyridines with variable nucleophilicity as catalysts to control deactivation (pyridine, 4-methylpyridine, 4-methoxypyridine, and 4-dimethylaminopyridine). Reactions have been monitored by NMR spectroscopy. Although this reaction network is simple, determination of well-defined rate constants from kinetic modeling is challenging because of correlation between the parameters. This issue can be addressed by analyzing the anhydride hydrolysis independently. The rate of attack of the pyridines on the anhydride follows expected nucleophilicity trends, although this is offset by increased protonation of more-nucleophilic pyridines at typical pH's. The optimized parameters can be used to model the full carbodiimide-driven process, although the presence of the common carbodiimide EDC has unanticipated effects on the anhydride hydrolysis rate. The results offer context for controlling carbodiimide-fueled reaction networks through the choice of suitable catalysts and pH.

Anhydrides↗
NASA NTRS2010

Modeling Nucleation and Grain Growth in the Solar Nebula: Initial Progress Report

The primitive solar nebula was a violent and chaotic environment where high energy collisions, lightning, shocks and magnetic re-connection events rapidly vaporized some fraction of nebular dust, melted larger particles while leaving the largest grains virtually undisturbed. At the same time, some tiny grains containing very easily disturbed noble gas signatures (e.g., small, pre-solar graphite or SiC particles) never experienced this violence, yet can be found directly adjacent to much larger meteoritic components (chondrules or CAIs) that did. Additional components in the matrix of the most primitive carbonaceous chondrites and in some chondritic porous interplanetary dust particles include tiny nebular condensates, aggregates of condensates and partially annealed aggregates. Grains formed in violent transient events in the solar nebula did not come to equilibrium with their surroundings. To understand the formation and textures of these materials as well as their nebular abundances we must rely on Nucleation Theory and kinetic models of grain growth, coagulation and annealing. Such models have been very uncertain in the past: we will discuss the steps we are taking to increase their reliability.

Nuth, Joseph A.↗
DOE OSTI2024

Ultrafast Exciton Dynamics of CH 3 NH 3 PbBr 3 Perovskite Nanoclusters

Exciton dynamics of perovskite nanoclusters has been investigated for the first time using femtosecond transient absorption (TA) and time-resolved photoluminescence (TRPL) spectroscopy. The TA results show two photoinduced absorption signals at 420 and 461 nm and a photoinduced bleach (PB) signal at 448 nm. The analysis of the PB recovery kinetic decay and kinetic model uncovered multiple processes contributing to electron-hole recombination. The fast component (~8 ps) is attributed to vibrational relaxation within the initial excited state, and the medium component (~60 ps) is attributed to shallow carrier trapping. The slow component is attributed to deep carrier trapping from the initial conduction band edge (~666 ps) and the shallow trap state (~40 ps). The TRPL reveals longer time dynamics, with modeled lifetimes of 6.6 and 93 ns attributed to recombination through the deep trap state and direct band edge recombination, respectively. In conclusion, the significant role of exciton trapping processes in the dynamics indicates that these highly confined nanoclusters have defect-rich surfaces.

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

Conformer-Dependent Chemistry: Experimental Product Branching of the Vinyl Alcohol + OH + O 2 Reaction

The concentration of formic acid in Earth's tropo-sphere is underestimated by detailed chemical models compared to field observations. Phototautomerization of acetaldehyde to its less stable tautomer vinyl alcohol, followed by the OH-initiated oxidation of vinyl alcohol, has been proposed as a missing source of formic acid that improves the agreement between models and field measurements. Theoretical investigations of the OH + vinyl alcohol reaction in excess O 2 conclude that OH addition to the alpha carbon of vinyl alcohol produces formaldehyde + formic acid + OH, whereas OH addition to the β site leads to glycoaldehyde + HO 2 . Furthermore, these studies predict that the conformeric structure of vinyl alcohol controls the reaction pathway, with the anti-conformer of vinyl alcohol promoting α OH addition, whereas the syn-conformer promotes β addition. However, the two theoretical studies reach different conclusions regarding which set of products dominate. We studied this reaction using time-resolved multiplexed photoionization mass spectrometry to quantify the product branching fractions. Our results, supported by a detailed kinetic model, conclude that the glycoaldehyde product channel (arising mostly from syn-vinyl alcohol) dominates over formic acid production with a 3.6:1.0 branching ratio. Additionally, this result supports the conclusion of Lei et al. that conformer-dependent hydrogen bonding at the transition state for OH-addition controls the reaction outcome. As a result, tropospheric oxidation of vinyl alcohol creates less formic acid than recently thought, increasing again the discrepancy between models and field observations of Earth's formic acid budget.

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

A thermodynamic approach to obtain materials properties for engineering applications

With the ever increases in the capabilities of computers for numerical computations, we are on the verge of using these tools to model manufacturing processes for improving the efficiency of these processes as well as the quality of the products. One such process is casting for the production of metals. However, in order to model metal casting processes in a meaningful way it is essential to have the basic properties of these materials in their molten state, solid state as well as in the mixed state of solid and liquid. Some of the properties needed may be considered as intrinsic such as the density, heat capacity or enthalpy of freezing of a pure metal, while others are not. For instance, the enthalpy of solidification of an alloy is not a defined thermodynamic quantity. Its value depends on the micro-segregation of the phases during the course of solidification. The objective of the present study is to present a thermodynamic approach to obtain some of the intrinsic properties and combining thermodynamics with kinetic models to estimate such quantities as the enthalpy of solidification of an alloy.

Chang, Y. Austin↗
DOE OSTI2022

CaMn 0.9 Ti 0.1 O 3 based redox catalysts for chemical looping – Oxidative dehydrogenation of ethane: Effects of Na 2 MoO 4 promoter and degree of reduction on the reaction kinetics

Reduction kinetics and stability of 20 wt% Na 2 MoO 4 -promoted CaMn 0.9 Ti 0.1 O 3 were investigated for its applications in Chemical Looping – Oxidative Dehydrogenation (CL-ODH) of ethane, a potential alternative for ethylene production with higher efficiency and lower emissions. Here, the present work reports a kinetics model and parameters for a Na 2 MoO 4 -promoted, Ti-doped CaMnO 3 (CaMn 0.9 Ti 0.1 O 3 ) redox catalyst under H 2 and C 2 H 4 . A first-order reaction model provides the best fit for the reduction of Na 2 MoO 4 /CaMn 0.9 Ti 0.1 O 3 under H 2 , while the C 2 H 4 reduction is well described by an Avrami–Erofe’ev model. The activation energy for C 2 H 4 oxidation is approximately three times higher than that for H 2 conversion, showing that the activation of C 2 H 4 is significantly more difficult on the surface of the redox catalyst. The reduction rate of Na 2 MoO 4 /CaMn 0.9 Ti 0.1 O 3 under H 2 at 750 °C is more than two orders of magnitude greater than that under C 2 H 4 , while the reduction rate of unpromoted CaMn 0.9 Ti 0.1 O 3 is comparable under H 2 and C 2 H 4 , showing that the addition of Na 2 MoO 4 effectively suppresses C 2 H 4 combustion relative to H 2 oxidation. The kinetics results for Na 2 MoO 4 /CaMn 0.9 Ti 0.1 O 3 confirm its excellent selectivity towards hydrogen combustion, making it a promising candidate under CL-ODH. Additionally, the stability of the CaMn 0.9 Ti 0.1 O 3 @ Na 2 MoO 4 core-shell structure, which was the underlying reason for the excellent selectivity, was examined under both shallow and deep reductions. It was determined that deep reduction of the redox catalyst, e.g. higher than 80% solid conversion, would lead to loss of sodium and hence to decreased selectivity for hydrogen combustion. In contrast, the core-shell structure was well-maintained, exhibiting excellent performance after 50 redox cycles when deep reduction of the redox catalyst was avoided. This study offers a basis for both the CL-ODH reactor design and redox catalyst optimizations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗