SEARCH · Engineering Papers

Results for “Diffusion rate”

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 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Nuclear Magnetic Resonance Dynamics of LiTFSI–Pyrazole Eutectic Solvents

Deep Eutectic Solvents (DESs) have emerged as promising candidates to replace conventional organic solvents in various technological applications due to their low vapor pressure, non-flammability, and ease of preparation at low costs. In particular, Type IV DESs, which are composed of metal salts and hydrogen bond donors, are possible replacements for lithium-ion battery electrolytes. In this study, we investigate the molecular dynamics of solvents of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) and pyrazole (PYR) at varying LiTFSI:PYR molar ratios (1:2, 1:3, 1:4, 1:5) using Nuclear Magnetic Resonance Dispersion (NMRD) and Pulsed Field Gradient (PFG) Nuclear Magnetic Resonance (NMR). PFG NMR reveals composition-dependent diffusion trends, while NMRD provides molecular-level insights into the longitudinal relaxation rate (R 1 = 1/T 1 ). Notably, the LiTFSI:PYR (1:2) sample shows distinct behavior across both techniques, exhibiting enhanced relaxation rates and lower self-diffusion for 1 H compared to the other nuclei ( 19 F and 7 Li), suggestive of stronger and more efficient Li + –pyrazole interactions, as confirmed by the modeling of the relaxation profiles. Our study advances understanding of ion dynamics in azole-based eutectic solvents, supporting their potential use in safer battery electrolytes.

FFC-NMR↗

DOE OSTI2021

Unlocking bulk and surface oxygen transport properties of mixed oxide-ion and electron conducting membranes with combined oxygen permeation cell and oxygen probe method

Surface exchange kinetics and bulk diffusion of oxygen are of paramount importance to the activity of oxygen electrocatalysis and performance of electrochemical devices such as fuel cell, metal-air batteries, and oxygen separation membranes. Conventional approaches to obtaining these transport properties are often limited to single property under a specific non-operation related condition. Here we use a combined oxygen permeation cell and oxygen probe methodology to simultaneously attain rates of oxygen surface exchange and bulk conductivity/chemical diffusivity of three representative mixed oxide-ion and electron conductors, namely SrCo 0.9 Ta 0.1 O 3-δ (SCT), La 0.6 Sr 0.4 CoO 3-δ (LSC) and La 0.6 Sr 0.4 FeO 3-δ (LSF), operated under a steady-state oxygen flux. The results explicitly show that SCT exhibit the highest oxide-ion conductivity/chemical diffusivity, fastest rates of surface oxygen exchange kinetics, thus promising to be the best oxygen electrocatalyst. Here, we have also mapped out the distribution of oxygen chemical potential gradient across the membranes and applied B-transport number concept to illustrate the rate-limiting steps in the overall oxygen permeation process.

42 ENGINEERING↗

DOE OSTI2020

Computational modeling of charge hopping dynamics along a disordered one-dimensional wire with energy gradients in quantum environments

This computational study investigates the effects of energy gradients on charge hopping dynamics along a one-dimensional chain of discrete sites coupled to quantum bath, which is modeled at the level of Pauli master equation (PME). This study also assesses the performance of different approximations for the hopping rates. Three different methods for solving the PME, a fourth order Runge–Kutta method, numerical diagonalization of the rate matrix followed by analytic propagation, and kinetic Monte Carlo simulation method, are tested and confirmed to produce virtually identical values of time dependent mean square displacement, diffusion constant, and mobility. Five different rate expressions, exact numerical evaluation of Fermi’s Golden Rule (FGR) rate, stationary phase interpolation (SPI) approximation, semiclassical approximation, classical Marcus rate, and Miller–Abrahams rate, are tested to help understand the effects of approximations in representing quantum environments in the presence of energy gradients. The results based on direct numerical evaluation of FGR rate exhibit transition from diffusive to non-diffusive behavior with the increase in the gradient and show that the charge transport in the quantum bath is more sensitive to the magnitude of the gradient and the disorder than in the classical bath. Among all the four approximations for the hopping rates, the SPI approximation is confirmed to work best overall. A comparison of two different methods to calculate the mobility identifies drift motion of the population distribution as the major source of non-diffusive behavior and provides more reliable information on the contribution of quantum bath.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2003

Finite-Rate Ablation Boundary Conditions for Carbon-Phenolic Heat-Shield

A formulation of finite-rate ablation surface boundary conditions, including oxidation, nitridation, and sublimation of carbonaceous material with pyrolysis gas injection, has been developed based on surface species mass conservation. These surface boundary conditions are discretized and integrated with a Navier-Stokes solver. This numerical procedure can predict aerothermal heating, chemical species concentration, and carbonaceous material ablation rate over the heatshield surface of re-entry space vehicles. In this study, the gas-gas and gas-surface interactions are established for air flow over a carbon-phenolic heatshield. Two finite-rate gas-surface interaction models are considered in the present study. The first model is based on the work of Park, and the second model includes the kinetics suggested by Zhluktov and Abe. Nineteen gas phase chemical reactions and four gas-surface interactions are considered in the present model. There is a total of fourteen gas phase chemical species, including five species for air and nine species for ablation products. Three test cases are studied in this paper. The first case is a graphite test model in the arc-jet stream; the second is a light weight Phenolic Impregnated Carbon Ablator at the Stardust re-entry peak heating conditions, and the third is a fully dense carbon-phenolic heatshield at the peak heating point of a proposed Mars Sample Return Earth Entry Vehicle. Predictions based on both finite-rate gas- surface interaction models are compared with those obtained using B' tables, which were created based on the chemical equilibrium assumption. Stagnation point convective heat fluxes predicted using Park's finite-rate model are far below those obtained from chemical equilibrium B' tables and Zhluktov's model. Recession predictions from Zhluktov's model are generally lower than those obtained from Park's model and chemical equilibrium B' tables. The effect of species mass diffusion on predicted ablation rate is also examined.

Chen, Y.-K.↗

NASA NTRS2004

Non-Additive Voltametric Currents From a Mixture of Two, Three and Four Redox-Active Compounds and Electroanalytical Implications

We have published recently the effect of dissimilar diffusion coefficients on the size of the voltammetric waves from a mixture of two redox-active compounds. Similarly, at the potential range where three redox-active species, decamethylferrocene (dMeFc), ferrocene (Fc) and N-methylphenothiazine (MePTZ), are oxidized simultaneously with rates controlled by linear diffusion, electrogenerated radicals diffusing outwards from the electrode react with the original species diffusing towards the electrode from the bulk; thus, Fc(+) reacts with dMeFc producing Fc and dMeFc(+), while MePTZ(+) reacts both with dMeFc producing MePTZ and dMeFc(+), and with Fc producing MePTZ and Fc(+). These reactions replace dMeFc with Fc at the second plateau, and both dMeFc and Fc with MePTZ at the third plateau. Since the diffusion coefficients of the three species are not equal, the mass-transfer limited currents of the second and the third oxidation wave plateaus change by approx. 10%. Numerical simulations of the experimental voltamograms support this mechanism. Similar results were also obtained for a mixture of four redoxactive compounds. The implications of this non-additive nature of currents on: (a) the use of internal voltammetric standards for quantitative analysis of a mixture of redox-active compounds; and, (b) the half wave potentials (E1/2) of the 2nd, 3rd and 4th waves for qualitative analysis, will be discussed.

Dass, Amala↗

DOE OSTI2020

Transport Processes in a Li-ion Cell during an Internal Short-Circuit

Internal short-circuit in a lithium-ion cell causes an abrupt increase in cell temperature and triggers subsequent thermal runaway. In this work, we present a detailed electrochemical-thermal model to investigate the physical behavior during an internal short-circuit. Simulations at wide range of heat transfer coefficients and short-circuit resistances are conducted to illustrate electrochemical and thermal behavior under a wide range of conditions. The Joule heating at the shorted region promotes electrochemical reactions nearby, causing in-plane non-uniformity of electrolyte and active material transport. Furthermore, it is found that diffusion in solid active materials plays a significant role at very high shorting currents (~20 C), as electrochemical reactions rate are being controlled progressively more by availability of Li + at the interface, due to limitations in diffusion through the active material with increasing discharge rates. This diffusion limitation causes a drop in available energy, and subsequently a decrease in cell equilibrium potential and the heat generation rate at the location of the short. On the other hand, rapid depletion of lithium concentration in the electrolyte and accumulation on the electrode surface results in highly non-uniform transport properties resulting in higher heat generation rates. Hence, the heating regime shifts from "local heating" to "global heating". Based on the findings, important design parameters for battery safety are discussed.

25 ENERGY STORAGE↗

DOE OSTI2021

Analytical modeling for redox flow battery design

Deeper market penetration of redox flow batteries requires optimization of the cell performance. Though important for optimization, detailed analytical solutions have not been developed for electrolyte flow, mass and charge transport, and reaction kinetics within redox flow batteries. To this end, here we present analytical solutions to active species concentration and over-potential based on advection-diffusion transport for ions and Bulter-Volmer model for interface reaction kinetics. The solutions were validated with results from a finite element model. These solutions were then applied to investigate the relationship between over-potential and state of charge, current density, reaction rate constant, flow velocity, diffusivity, total active species concentration, and electrode structure. Explicit formulas were identified for minimum activation over-potential and limiting current density as well as their dependence on electrolyte properties, operation conditions, and electrode structure. With our new mathematical formulas, this work provides a theoretical framework for flow battery design.

25 ENERGY STORAGE↗

NASA NTRS1986

A rapid and accurate method for calculation of stratospheric photolysis rates with molecular scattering

A method for calculating the photodissociation rates needed for photochemical modeling of the stratosphere, which includes the effects of molecular scattering, is described. The procedure is based on Sokolov's method of averaging functional correction. The radiation model and approximations used to calculate the radiation field are examined. The approximated diffuse fields and photolysis rates are compared with exact data. It is observed that the approximate solutions differ from the exact result by 10 percent or less at altitudes above 15 km; the photolysis rates differ from the exact rates by less than 5 percent for altitudes above 10 km and all zenith angles, and by less than 1 percent for altitudes above 15 km.

Boughner, Robert E.↗

DOE OSTI2023

Mass transfer in catalytic depolymerization: External effectiveness factors and serendipitous processivity in stagnant and stirred melts

Several heterogeneous catalysts are being developed to recycle plastics. Most operate in viscous polymer melts, where external mass transfer effects could limit the supply of co-reactants to active sites. External mass transfer can also impede the diffusion of long chain products away from the catalyst after each cut. Product egress limitations could potentially confer unintentional processivity to catalyst operation, i.e. a tendency for the catalyst to repeatedly cut the same chain after an initial encounter. We formulate reaction–diffusion equations to quantify mass transfer effects on the co-reactant transport to the catalyst and the degree of serendipitous processivity. Results are developed for catalysts in stagnant or stirred melts, with simple expressions involving Damkohler, Peclet, and Sherwood numbers, i.e. dimensionless combinations of rate constants, catalyst particle size, polymer diffusivities, and shear rates (where applicable). In conclusion, we estimate the impact of these effects for a spherical core–shell catalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Visualization and analysis of coupling between plasmas self-organization and plasma-induced fluid circulation in 1 atm DC glows with liquid anode

The physical processes prevailing when plasma contacts liquid water are poorly understood; however, it is this very interaction that is the basis of new, exciting technologies for water purification, the treatment of cancer and disease and the production of new, high value chemical products. A chief impediment to widespread application of plasmas for the aforementioned applications is scale up and cost. Here we investigate plasma induced flows which facilitate transport of reactivity from the plasma to liquid water at rates much faster than diffusion. We find that this plasma induced motion varies depending on whether or not the plasma is self organized. The insight from this work paves the way to a better understanding of transport and ultimately informing plasma treatment technologies on how to best optimize the rate at which processes such as plasma based water purification occurs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Modeling mesoscale fission gas behavior in UO2 by directly coupling the phase field method to spatially resolved cluster dynamics

Abstract Fission gas release within uranium dioxide nuclear fuel occurs as gas atoms diffuse through grains and arrive at grain boundary (GB) bubbles; these GB bubbles grow and interconnect with grain edge bubbles; and grain edge tunnels grow and connect to free surfaces. In this study, a hybrid multi-scale/multi-physics simulation approach is presented to investigate these mechanisms of fission gas release at the mesoscale. In this approach, fission gas production, diffusion, clustering to form intragranular bubbles, and re-solution within grains are included using spatially resolved cluster dynamics in the Xolotl code. GB migration and intergranular bubble growth and coalescence are included using the phase field method in the MARMOT code. This hybrid model couples Xolotl to MARMOT using the MultiApp and Transfer systems in the MOOSE framework, with Xolotl passing the arrival rate of gas atoms at GBs and intergranular bubble surfaces to MARMOT and MARMOT passing evolved GBs and bubble surface positions to Xolotl. The coupled approach performs well on the two-dimensional simulations performed in this work, producing similar results to the standard phase field model when Xolotl does not include fission gas clustering or re-solution. The hybrid model performs well computationally, with a negligible cost of coupling Xolotl and MARMOT and good parallel scalability. The hybrid model predicts that intragranular fission gas clustering and bubble formation results in up to 70% of the fission gas being trapped within grains, causing the increase in the intergranular bubble fraction to slow by a factor of six. Re-solution has a small impact on the fission gas behavior at 1800 K but it has a much larger impact at 1000 K, resulting in a twenty-times increase in the concentration of single gas atoms within grains. Due to the low diffusion rate, this increase in mobile gas atoms only results in a small acceleration in the growth of the intergranular bubble fraction. Finally, the hybrid model accounts for migrating GBs sweeping up gas atoms. This results in faster intergranular bubble growth with smaller initial grain sizes, since the additional GB migration results in more immobile gas clusters reaching GBs.

Kim, Dong-Uk↗

NASA NTRSDate not supplied

Effect of Alloying Additions on Twinning in Ni-based Superalloys

Micro-twinning is the dominant creep deformation mechanism in Ni-based superalloys at temperatures above 700 °C. We use atomistic simulations to study two mechanisms of twin nucleation and growth that are characterized by qualitatively different rate limiting processes. In case of the mechanism proposed by Kolbe, the rate limiting process is diffusion-mediated atomic reshuffling. In case of the other mechanism, we proposed recently, the rate limiting process is nucleation of Shockley partial dislocation. We demonstrate the effects of alloying additions on functionality of these mechanisms.

Valery V Borovikov↗

NASA NTRSDate not supplied

Effect of Alloying Additions on Twinning in Ni-Based Superalloys

Micro-twinning is the dominant creep deformation mechanism in Ni-based superalloys at temperatures above 700 C. We use atomistic simulations to study two mechanisms of twin nucleation and growth that are characterized by qualitatively different rate limiting processes. In case of the mechanism proposed by Kolbe, the rate limiting process is diffusion-mediated atomic reshuffling. In case of the other mechanism, we proposed recently, the rate limiting process is nucleation of Shockley partial dislocation. We demonstrate the effects of alloying additions on functionality of these mechanisms.

Ni-based superalloys↗

DOE OSTI2022

Mechanistic calculation of the effective silver diffusion coefficient in polycrystalline silicon carbide: Application to silver release in AGR-1 TRISO particles

The silicon carbide (SiC) layer in tristructural isotropic (TRISO) fuel particles serves as a barrier to prevent the escape of fission products produced and not retained in the fuel kernel. The release of silver (Ag) is a concern due to the long half-life of the 110m Ag isotope. However, accurately determining the fission gas release rate requires knowing the diffusion coefficient through the SiC layer. In this study, we leverage atomistic calculations of Ag diffusivity in SiC bulk and grain boundaries (GBs) to develop a mesoscale effective Ag diffusion coefficient (D eff ) in SiC. Since GBs serve as pathways for Ag diffusion, D eff is defined as a function of temperature and microstructure variables. In particular, the size of SiC grains in the direction perpendicular to diffusion is shown to significantly affect Ag diffusion. The prediction of the mechanistic, mesoscale approach falls within one order of magnitude of empirical values. The temperature and microstructure-dependent effective Ag diffusivity in SiC is implemented in the fuel performance code Bison with a correction factor to predict Ag release from AGR-1 TRISO fuel particles. We hereby quantify the impact of SiC grain size on Ag release and improve Bison’s predictions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2001

Suppression of Low Strain Rate Nonpremixed Flames by an Agent

The agent concentration required to achieve the suppression of low strain rate nonpremixed flames is an important consideration for fire protection in a microgravity environment such as a space platform. Currently, there is a lack of understanding of the structure and extinction of low strain rate (<20 s(exp -1)) nonpremixed flames. The exception to this statement is the study by Maruta et al., who reported measurements of low strain rate suppression of methane-air diffusion flames with N2 added to the fuel stream under microgravity conditions. They found that the nitrogen concentration required to achieve extinction increased as the strain rate decreased until a critical value was obtained. As the strain rate was further decreased, the required N2 concentration decreased. This phenomenon was termed "turning point" behavior and was attributed to radiation-induced nonpremixed flame extinction. In terms of fire safety, a critical agent concentration assuring suppression under all flow conditions represents a fundamental limit for nonpremixed flames. Counterflow flames are a convenient configuration for control of the flame strain rate. In high and moderately strained near-extinction nonpremixed flames, analysis of flame structure typically neglects radiant energy loss because the flames are nonluminous and the hot gas species are confined to a thin reaction zone. In counterflowing CH4-air flames, for example, radiative heat loss fractions ranging from 1 to 6 percent have been predicted and measured. The objective of this study is to investigate the impact of radiative emission, flame strain, agent addition, and buoyancy on the structure and extinction of low strain rate nonpremixed flames through measurements and comparison with flame simulations. The suppression effectiveness of a number of suppressants (N2, CO2, or CF3Br) was considered as they were added to either the fuel or oxidizer streams of low strain rate methane-air diffusion flames.

Hamins, A.↗

DOE OSTI2022

Hybrid quantum-classical simulations of magic angle spinning dynamic nuclear polarization in very large spin systems

Solid-state nuclear magnetic resonance can be enhanced using unpaired electron spins with a method known as dynamic nuclear polarization (DNP). Fundamentally, DNP involves ensembles of thousands of spins, a scale that is difficult to match computationally. This scale prevents us from gaining a complete understanding of the spin dynamics and applying simulations to design sample formulations. We recently developed an ab initio model capable of calculating DNP enhancements in systems of up to ~1000 nuclei; however, this scale is insufficient to accurately simulate the dependence of DNP enhancements on radical concentration or magic angle spinning (MAS) frequency. We build on this work by using ab initio simulations to train a hybrid model that makes use of a rate matrix to treat nuclear spin diffusion. We show that this model can reproduce the MAS rate and concentration dependence of DNP enhancements and build-up time constants. Finally, we then apply it to predict the DNP enhancements in core–shell metal-organic-framework nanoparticles and reveal new insights into the composition of the particles’ shells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Evaluation of Oxygen Consumption as a Sensitive Measure of Electrical Cable Polymer Insulation Degradation

Polymers used as dielectric materials in nuclear grade electrical cables, such as cross-linked polyethylene (XLPE) and ethylene-propylene rubber (EPR), are long-lived and durable. Estimating the decades-long useful service life of polymers experimentally in a feasible timeframe of weeks or months requires greatly accelerated aging in the laboratory at elevated temperatures and/or gamma dose rates. Many limitations of accelerated aging, such as due to artefacts of the process including diffusion limited oxidation and dose rate effects, have been identified that can make accurate lifetime prediction at service conditions challenging. This can be due to differences in degradation mechanisms between rapid and extended aging or to the presence of induction periods during which early aging is undetectable using common metrics of polymer insulation characterization such as tensile elongation at break. Conceptually, an ideal measure of aging would track clearly and predictably with material changes throughout the lifecycle of the evolving material state from early aging through end of useful life. With such a metric accelerated aging could be performed at much milder conditions, conditions closer to those found in long term service, and end of life calculated from the early data. Oxygen consumption has previously been proposed as a sensitive measure of insulation polymer oxidation that might meet such desirable characteristics. In this work, nuclear cable insulation samples of the two most common types, XLPE and EPR, from two of the most sourced manufacturers, RSCC and The Okonite Company, were subjected to gamma irradiation at dose rates of 100, 200, and 300 Gy/hr for up to 42 days at 26°C. The aged samples were used to evaluate oxygen consumption as a sensitive measure of aging for these materials and these conditions. For comparison and validation, the aged samples were also characterized using the more conventional tests of mass change, density, carbonyl index, and yellowness index. For both materials, oxygen consumption during irradiation was found to track consistently with applied dose—oxygen consumption rate was linear with time. Mass change and carbonyl index were also clearly proportional to exposure dose. Yellowness index, a surface sensitive measure, was observed to increase linearly with dose at low exposure levels before plateauing. Sample density was not observed to change under the conditions explored. Oxygen consumption was found to be a particularly sensitive metric to measure the degradation of polymeric cable insulation exposed to radiation. Oxygen consumption was compared to Arrhenius and time-temperature-superposition methods for activation energy calculation using literature date and found to be a useful and viable alternative for lifetime prediction. The techniques developed in this work are anticipated to be valuable in the estimation of insulation aging rates under mild conditions approaching service conditions. Direction comparison of useful life prediction by oxygen consumption with prediction via Arrhenius methodology for thermal aging and existing methods for radiation and combined thermal/radiation aging would serve to further validate this promising approach and refine our understanding of material-specific aging rates.

36 MATERIALS SCIENCE↗