SEARCH · Search NASA

Results for “Interface kinetics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

DOE OSTI2012

Structure Determination and Functional Analysis of a Chromate Reductase from Gluconacetobacter hansenii

N/A

Jin, Hongjun↗

DOE OSTI2020

Small- and Wide-Angle X-ray Scattering Studies of Liquid–Liquid Phase Separation in Silicate Melts

Incongruent liquid-liquid phase separation is a common phenomenon in molten oxides, but there are few structural studies due to the high temperatures involved ca. >1600oC. In this work, we present a high-energy x-ray combined small/wide angle scattering (SAXS/WAXS) method that can directly probe non-equilibrium phase diagrams. The Porod exponent in the SAXS signal is extracted, which measures the roughness of the interface formed between different phases. Problematic volatilization of silica in Al 2 O 3 -, CaO-, and MgO-bearing silica-rich liquids were minimized using levitation, laser heating and rapid time-resolved measurements (lasting tens of seconds), to probe the kinetics of the phase separation process. The WAXS measurement simultaneously probes the local atomic structure of the liquid state and shows that the primary structural changes are associated with the SiO 2 network. Cation coordination numbers extracted from pair distribution functions support the model that the origin of phase separation in binary silicate systems is due to Coulombic repulsions between poorly screened cations. For CaO-SiO 2 , the liquid-liquid phase separation is found to correlate with a pause in the temperature dependent atomic disorder associated with SiO4 tetrahedra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations

Traditional Kinetic Monte Carlo (KMC) approaches, rooted in Gillespie’s stochastic simulation algorithm, become computationally demanding in systems with a large range of timescales. The goal of this work is to propose and study new adaptive lattice-KMC time integration strategies for spatially non-uniform systems. To that end, two novel adaptive tau-leaping methods and their corresponding time integration strategies are developed based on the idea of the “n-fold” direct KMC method. These strategies allow for the simultaneous execution of multiple reactions, advancing time by adaptively selected coarse increments. We present numerical experiments comparing the proposed methods with existing approaches in a catalytic surface kinetics application involving ammonia decomposition.

Bimolecular reactions↗

DOE OSTI2016

Kinetic studies of ICF implosions

Here, kinetic effects on inertial confinement fusion have been investigated. In particular, inter-ion-species diffusion and suprathermal ion distribution have been analyzed. The former drives separation of the fuel constituents in the hot reacting core and governs mix at the shell/fuel interface. The latter underlie measurements obtained with nuclear diagnostics, including the fusion yield and inferred ion burn temperatures. Basic mechanisms behind and practical consequences from these effects are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Many-body effects in transport and energy transfer

Many-body physics is common in defining responses of nanoscale to external perturbations such as bias, optical driving, temperature gradient. When applied to open systems these perturbations lead to matter and energy transfer between the system and its surroundings (baths) and/or cause intra-system dynamics. While many different theoretical techniques have been developed to describe system responses, not all of them are convenient (or even appropriate) when applied at nanoscale. Traditional theoretical techniques are either based on kinetic schemes (e.g., Redfield-Lindblad quantum master equation) or utilize quasiparticle description (e.g., single-particle scattering theory, input-output formalism, standard nonequilibrium Green’s functions technique). The former fails to account for quantum mechanics of the system-bath interface and may fail qualitatively when used in open interacting systems. The latter is inconvenient when many-body states description of the system is preferable, which is usual, e.g., in molecular optical spectroscopy. This project was devoted to development of new nonequilibrium Green’s function (NEGF) techniques and application of standard and many-body flavors of the NEGF to description of charge and energy transport, optical response, nonadiabatic molecular dynamics, and thermodynamics of nanoscale devices in general, and, in particular, of single-molecule junctions. In particular, Hubbard NEGF (new many-body flavor of the NEGF) was introduced as a new technique especially convenient in treatment of optoelectronic experiments in single-molecule junctions and useful in discussion of bias induced nuclear dynamics in junctions. The research also led to further technique developments and applications beyond the timing of this Award.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2023

Zero-Order Reaction Kinetics (Zero-RK): Enabling the Use of Detailed Chemical Kinetics in Combustion Simulations (CRADA Final Report)

This was a collaborative effort between Lawrence Livermore National Security, LLC (LLNS), as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Gamma Technologies, LLC (GT or Participante), to incorporate the ability to access LLNL chemical kinetics technologies while using GT-SUITE, GT’s market leading engine simulation software. At the end of the project, LLNL has released Zero-RK version 3.5 with zero- and one-dimensional (0-D and 1-D) solver functionality that interfaces with GT’s GT-SUITE v2023 and later releases. GT has tested its product to assure their customers that the interface can provide reduction in chemistry solution time for detailed chemistry simulations. The process has also positioned GT to easily benefit from future improvements of the Zero-RK suite of tools developed under the DOE Vehicle Technologies Office.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2026

Zero Order Reactioin Kinetics: Enabling the use of Detailed Chemical Kinetics in Combustion Simulations (Final CRADA Report)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1990

Compatibility of Fe-40Al with various fibers

Chemical reaction can occur at the fiber/matrix interface of intermetallic matrix composites, leading to a degradation of mechanical properties. Fe-40Al matrix composites were fabricated using SiC, B, W, Mo-base, and Al2O3 fibers. Composite samples were heat treated up to 1500 K to study the reaction kinetics, and reaction rates were determined from reaction zone thickness measurements. The Al2O3 and W fibers were found to be compatible with the Fe-40Al matrix, while the Mo-based fibers reacted moderately and the B and SiC fibers reacted severely. Experimental results are compared to theoretical thermodynamic predictions.

Draper, S. L.↗

NASA NTRS1999

In Situ Monitoring of Crystal Growth Using MEPHISTO, Mission STS 87-Program USMP-4: Experimental Analysis

This report summarizes the results of the In situ Monitoring of Crystal Growth Using MEPHISTO (Material por l'Etude des Phenomenes Interessant de la Solidification sur Terre et en Orbite) experiment on USMP-4. The report includes microstructural and compositional data obtained during the first year of the post flight analysis, as well as numerical simulation of the flight experiment. Additional analyses are being continued and will be reported in the near future. The experiments utilized MEPHISTO hardware to study the solidification and melting behavior of bismuth alloyed with 1 at% tin. The experiments involved repeated melting and solidification of three samples, each approximately 90 cm long and 6mm in diameter. Half of each sample also included a 2 mm. diameter growth capillary, to assist in the formation of single grain inside. One sample provided the Seebeck voltage generated during melting and freezing processes. Another one provided temperature data and Peltier pulsed demarcation of the interface shape for post flight analysis. The third sample provided resistance and velocity measurements, as well as additional thermal data. The third sample was also quenched at the end of the mission to preserve the interface composition for post flight determination. A total of more than 45cm of directionally solidified alloy were directionally solidified at the end of the flight for post mission structural and compositional characterization. Metallurgical analysis of the samples has shown that the interfacial kinetics play a key role in controlling the morphological stability of faceted alloys. Substantial differences were observed in the Seebeck signal between the ground-based experiments and the space-based experiments. The temperature gradient in the liquid for the ground-based experiments was also significantly lower than the temperature gradient in the liquid for the space-based experiments. Both of these observations indicate significant influence of liquid convection for the ground-based experiments.

Abbaschian, Reza↗

DOE OSTI2023

Stabilizing High–Nickel Cathodes with High–Voltage Electrolytes

Electrolytes connect the two electrodes in a lithium battery by providing Li + transport channels between them. Advanced electrolytes are being explored with high-nickel cathodes and the lithium-metal anode to meet the high energy density and cycle life goals, but the origin of the performance differences with different electrolytes is not fully understood. Here, the mechanisms involved in protecting the high-capacity, cobalt-free cathode LiNiO 2 with a model high-voltage electrolyte (HVE) are delineated. The kinetic barrier posed by a thick surface degradation layer with poor Li + -ion transport is found to be the major contributor to the fast capacity fade of LiNiO 2 with the conventional carbonate electrolyte. In contrast, HVE reduces the side reactions between the electrolyte and the electrodes, leading to a thinner nano-interphase layer comprised of more beneficial species. Crucially, the HVE leads to a different surface reorganization pathway involving the formation of a thinner nanoscale LiNi 2 O 4 spinel phase on the LiNiO 2 surface. Finally, with a high 3D Li + -ion and electronic conductivity, the spinel LiNi 2 O 4 reorganization nanolayer preserves fast Li + transport across the cathode–electrolyte interface, reduces reaction heterogeneity in the electrode and alleviates intergranular cracking within secondary particles, resulting in superior long-term cycle life.

25 ENERGY STORAGE↗

DOE OSTI2016

Predicting oxidation-limited lifetime of thin-walled components of NiCrW alloy 230

Using alloy 230 as an example, a generalized oxidation lifetime model for chromia-forming Ni-base wrought alloys is proposed, which captures the most important damaging oxidation effects relevant for component design: wall thickness loss, scale spallation, and the occurrence of breakaway oxidation. For deriving input parameters and for verification of the model approach, alloy 230 specimens with different thicknesses were exposed for different times at temperatures in the range 950–1050 °C in static air. The studies focused on thin specimens (0.2–0.5 mm) to obtain data for critical subscale depletion processes resulting in breakaway oxidation within reasonably achievable test times up to 3000 h. The oxidation kinetics and oxidation-induced subscale microstructural changes were determined by combining gravimetric data with results from scanning electron microscopy with energy dispersive X-ray spectroscopy. The modeling of the scale spallation and re-formation was based on the NASA cyclic oxidation spallation program, while a new model was developed to describe accelerated oxidation occurring after longer exposure times in the thinnest specimens. The calculated oxidation data were combined with the reservoir model equation, by means of which the relation between the consumption and the remaining concentration of Cr in the alloy was established as a function of temperature and specimen thickness. Based on this approach, a generalized lifetime diagram is proposed, in which wall thickness loss is plotted as a function of time, initial specimen thickness, and temperature. As a result, the time to reach a critical Cr level at the scale/alloy interface of 10 wt% is also indicated in the diagrams.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Population balance modeling of polyurethane foam formation with pressure‐dependent growth kernel

Abstract Polyurethane foams are widely used materials often chosen for their useful characteristics such as low thermal conductivity, ease of application, and high strength‐to‐weight ratios. Computational models are needed to predict the dynamics of the flow and expansion, and the resulting material properties, to improve manufacturing processes. In this paper, a model for PMDI, a water‐blown polyurethane foam, is presented. By extending a kinetics‐based approach by adding bubble‐scale information via a population balance equation (PBE) using the quadrature method of moments, we can track bubble size distributions during foaming. We present results from a three‐dimensional computational fluid dynamics model using arbitrary Lagrangian–Eulerian interface tracking implemented in finite element software. The model compares favorably with experimental data, including dynamics, bubble distributions measured by both camera and diffusion wave spectroscopy, and post‐test bubble size from scanning electron microscopy and density measurements from x‐ray computed tomography.

Ortiz, Weston↗

DOE OSTI2012

Evaluation of strain caused by coherent precipitates in an Al alloy using TEM techniques

Elastic strains, caused by GP zones in an aged Al alloy, were determined quantitatively using two techniques: Dark Field In-line Holography (DFH) and High Resolution Transmission Electron Microscopy-Geometric Phase Analysis (HRTEM-GPA). The results obtained by both techniques showed that the elastic strain was not uniform along the precipitate-matrix interface. In some areas, it was found that strain had negligible value and this was attributed to the loss of coherence between the lattices. It is suggested that a possible explanation for this fact could be a variation in the 'vacancies pump mechanism' kinetics. To obtain a better interpretation of the experimental deformation maps, a reference GP precipitate-matrix structure was built using QSTEM software. The main advantages of DFH over HRTEM-GPA were a bigger field of view and low electron dose requirements without spatial resolution loss. Another difference found was that crystalline defects such as dislocations were evidenced by HRTEM-GPA in contrast to the result obtained by DFH. However, strain measurements were affected by mask size effect in the former. - Highlights: Black-Right-Pointing-Pointer Strain around GP precipitates was not uniform and was in a more relaxed state. Black-Right-Pointing-Pointer Strain values determined by HRTEM had more artifacts compared with holography ones. Black-Right-Pointing-Pointer Strain distribution was attributed to differences in the release of vacancies. Black-Right-Pointing-Pointer Holography showed important advantages over HRTEM regarding field of view mainly.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Elucidating the Discharge Behavior of Aqueous Zinc Sulfur Batteries in the Presence of Molybdenum(IV) Chalcogenide Catalyst: The Criticality of Interfacial Electrochemistry

The aqueous zinc-sulfur battery holds promise for significant capacity and energy density with low cost and safe operation based on environmentally benign materials. However, it suffers from the sluggish kinetics of the conversion reaction. Here, we highlight the efficacy of molybdenum(IV) sulfide (MoS 2 ) to reduce the overpotential of S-ZnS conversion in aqueous electrolytes and study the discharge products formed at the solid-solid and solid-liquid interfaces using experimental and theoretical approaches. Specifically, the MoS 2 -catalyzed electrochemical conversion reaction is characterized via ex situ X-ray diffraction (XRD), transmission electron microscopy (TEM) with energy dispersive spectroscopy (EDS), Raman spectroscopy, synchrotron-based Mo K-edge X-ray absorption spectroscopy (XAS), and in situ synchrotron-based X-ray computed tomography (XCT). Additionally, operando synchrotron-based S K-edge XAS and X-ray fluorescence (XRF) maps are collected to determine the spatial evolution of sulfur-based species at the electrode-electrolyte interface. Further, coupling the operando S K-edge XAS data with the simulated spectra and fitting the data suggested a possible ZnS 2 intermediate phase.

25 ENERGY STORAGE↗

DOE OSTI2019

Edge turbulence evolution and intermittency development near the density limit on the HL-2A tokamak

The development of intermittent non-Gaussian processes is researched in the edge turbulence of ohmically heated HL-2A discharges approaching the density limit. As the density limit is approached, the E × B shear flow at the last closed flux surface (LCFS) weakens, a strong positive skewness develops in the scrape-off layer (SOL), and negative skewness develops inside the LCFS of turbulent density fluctuations. A conditional averaging analysis confirms more frequent increased amplitude positive (negative) going density fluctuation activity in the SOL (inside the LCFS) as the density limit is approached. The measured turbulent stress across the edge, LCFS, and SOL region is decomposed into diffusive and residual stress components, and the nonlinear exchange of kinetic energy between the turbulence and the low-frequency E × B shear flow is determined. Residual stress acts to amplify the E × B flow at the LCFS, while the diffusive stress acts to dissipate the flow just inside this region, at the interface between the core plasma and the LCFS. The relative strength of the flow drive associated with the residual stress weakens as the density limit is approached, while the turbulent viscosity associated with the diffusive stress increases at high density. The adiabatic parameter, k z 2 v th 2 / ω ν e , drops significantly to about 0.5 in the SOL when the density limit is approached, indicating a transition from the adiabatic regime to the hydrodynamic regime due to increased collisionality. Such changes enhance the particle transport through the nonadiabatic electron response and hence should end in a stronger edge cooling at fixed edge plasma heat flux.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Critical initial-slip scaling for the noisy complex Ginzburg–Landau equation

We employ the perturbative fieldtheoretic renormalization group method to investigate the universal critical behavior near the continuous non-equilibrium phase transition in the complex Ginzburg–Landau equation with additive white noise. This stochastic partial differential describes a remarkably wide range of physical systems: coupled nonlinear oscillators subject to external noise near a Hopf bifurcation instability; spontaneous structure formation in non-equilibrium systems, e.g., in cyclically competing populations; and driven- dissipative Bose–Einstein condensation, realized in open systems on the interface of quantum optics and many-body physics, such as cold atomic gases and exciton-polaritons in pumped semiconductor quantum wells in optical cavities. Our starting point is a noisy, dissipative Gross–Pitaevski or nonlinear Schrödinger equation, or equivalently purely relaxational kinetics originating from a complex-valued Landau–Ginzburg functional, which generalizes the standard equilibrium model A critical dynamics of a non-conserved complex order parameter field. We study the universal critical behavior of this system in the early stages of its relaxation from a Gaussian-weighted fully randomized initial state. In this critical aging regime, time translation invariance is broken, and the dynamics is characterized by the stationary static and dynamic critical exponents, as well as an independent ‘initial-slip’ exponent. We show that to first order in the dimensional expansion about the upper critical dimension, this initial-slip exponent in the complex Ginzburg–Landau equation is identical to its equilibrium model A counterpart. Here, we furthermore employ the renormalization group flow equations as well as construct a suitable complex spherical model extension to argue that this conclusion likely remains true to all orders in the perturbation expansion.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1972

The effect of alloy composition on the mechanism of stress-corrosion cracking of titanium alloys in aqueous environments

A detailed study was made of the relation between the size distribution of Ti3Al particles in a Ti-8Al alloy and the tensile properties measured in air and in saltwater. The size distribution of Ti3Al was varied by isothermal aging for various times at temperatures in the range 770 to 970 K (930 to 1290 F). The aging kinetics were found to be relatively slow. Quantitative measurements of the particle coarsening rate at 920 K (1200 F) showed good agreement with the predicted behavior for coarsening controlled by matrix diffusion, and suggested that the specific free energy of the Ti3Al alpha interface in negligible small. In all cases, the Ti3Al particles were sheared by the glide dislocations. It was concluded that there is a definite correlation between the presence of deformable Ti3Al particles and an alloy's susceptibility to aqueous stress corrosion cracking. Furthermore, the appearance of the surface slip lines and the dislocation substructure in deformed specimens suggest that the specific effect of the Ti3Al particles is to cause a nonhomogeneous planar slip character and an enhanced chemical potential of the slip bands.

Wood, R. A.↗