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 775 records · Page 43

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Steady spherical hypercritical accretion onto neutron stars

The present study of hypercritical accretion flows onto neutron stars considers steady-state, spherically symmetric flows whose accretion rate range is characterized by the carrying away of gravitational-accretion energy by neutrinos. The models used encompass pair production, radiation diffusion, and general relativistic effects. While pair pressure dominates throughout the accretion envelope when accretion rates above about 100 solar masses/yr, radiation diffusion becomes important when the accretion rate falls below about 0.001 solar masses/yr. At the lower accretion rates, free fall toward the neutron-star surface stops and an extended, quasi-static, radiation pressure-supported envelope emerges which is probably dynamically unstable.

Houck, John C.↗

NASA NTRS1992

The compositional structure and the effects of exothermicity in an unpremixed planar jet flame

Results are presented of direct numerical simulations (DNS) of a randomly perturbed compressible, spatially developing planar jet under the influence of a finite rate Arrhenius chemical reaction of the type F + O yields Product + Heat. The objectives of the simulations are to assess the compositional structure of the flame and to determine the influence of reaction exothermicity by means of statistical sampling of the data generated by DNS. It is shown that even with this idealized kinetics model, the simulated results exhibit features in accordance with experimental data. These results indicate that the Damkoehler number is an important parameter in determining the statistical composition of the reacting field, and that the results are insensitive to the mechanism by which this parameter is varied. It is demonstrated that as the intensity of mixing is increased and the effect of finite rate chemistry is more pronounced, the magnitudes of the ensemble mean and rms of the product mass fraction decrease, and those of the reactants mass fraction increase. Also, at higher mixing rates the joint probability density functions of the reactants' mass fractions shift towards higher values within the composition domain indicating a lower reactedness. These trends are consistent with those observed experimentally and are very useful in portraying the statistical structure of non-equilibrium diffusion flames. The DNS generated data are also utilized to examine the applicability of the 'laminar diffusion flamelet model' in predicting the rate of the reactant conversion with finite rate chemistry. This examination indicates that the bounds of the product formation scale reasonably well with those obtained by the flamelet model. However, in the range of the relatively low values of the Damkoehler numbers considered, the scatter of the results is substantially more than that to be completely determined by the closure.

Steinberger, C. J.↗

DOE OSTI2022

Fast estimation of reaction rates in spherical and non-spherical porous catalysts

Here we present a methodology for modeling multi-step reaction rates in porous catalyst particles for use in CFD-DEM and two fluid models. Single-step effectiveness factors based on a Thiele modulus, while useful, cannot accurately capture the cascading reaction systems common in high temperature vapor-phase chemical reactors like fluidized catalytic cracking units and catalytic biomass fast pyrolysis systems. Instead, multi-step effectiveness vectors derived from steady-state solutions to the governing reaction-diffusion equations are needed. Solutions for various catalyst shapes are presented, including spheres, cylinders, and prisms. Computational challenges inherent in repeated evaluation of reaction rates with diffusion limitations are discussed, and an efficient implementation based on pre-computed lookup tables is proposed and demonstrated on a simulation of a fluidized bed reactor. Open-source code is provided for the compilation of reaction rate tables for use in ODE, DEM, and two-fluid models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

NASA NTRS2007

Effect of EMIC Wave Normal Angle Distribution on Relativistic Electron Scattering in Outer RB

We present the equatorial and bounce average pitch angle diffusion coefficients for scattering of relativistic electrons by the H+ mode of EMIC waves. Both the model (prescribed) and self consistent distributions over the wave normal angle are considered. The main results of our calculation can be summarized as follows: First, in comparison with field aligned waves, the intermediate and highly oblique waves reduce the pitch angle range subject to diffusion, and strongly suppress the scattering rate for low energy electrons (E less than 2 MeV). Second, for electron energies greater than 5 MeV, the |n| = 1 resonances operate only in a narrow region at large pitch-angles, and despite their greatest contribution in case of field aligned waves, cannot cause electron diffusion into the loss cone. For those energies, oblique waves at |n| greater than 1 resonances are more effective, extending the range of pitch angle diffusion down to the loss cone boundary, and increasing diffusion at small pitch angles by orders of magnitude.

Khazanov, G. V.↗

NASA NTRS1984

The reactions of cobalt, iron and nickel in SO-2 atmospheres: Similarities and differences

The reactions of cobalt, iron and nickel in SO2 atmospheres are reviewed and compared. A mixed oxide-sulfide product layer is observed in all cases. Cobalt and nickel exhibits similar behavior. The observed rates are near the sulfidation rates, and the reaction rate is strongly influenced by the outward diffusion of metal through an interconnected sulfide network. A continuous interconnected sulfide is not observed in the oxide-sulfide scales formed on iron, and the reaction rates are more difficult to summarize. The differences and similarities among the three metals are explained in terms of the absence of scale-gas equilibrium and the ratio of the metal diffusivity in the corresponding oxide and sulfide.

Jacobson, N. S.↗

NASA NTRS1985

The reactions of cobalt, iron and nickel in SO2 atmospheres Similarities and differences

The reactions of cobalt, iron and nickel in SO2 atmospheres are reviewed and compared. A mixed oxide-sulfide product layer is observed in all cases. Cobalt and nickel exhibit similar behavior. The observed rates are near the sulfidation rates, and the reaction rate is strongly influenced by the outward diffusion of metal through an interconnected sulfide network. A continuous interconnected sulfide is not observed in the oxide-sulfide scales formed on iron, and the reaction rates are more difficult to summarize. The differences and similarities among the three metals are explained in terms of the absence of scale-gas equilibrium and the ratio of the metal diffusivity in the corresponding oxide and sulfide.

Jacobson, N. S.↗

NASA NTRS2001

Mathematical Modeling of Ni/H2 and Li-Ion Batteries

The modelling effort outlined in this viewgraph presentation encompasses the following topics: 1) Electrochemical Deposition of Nickel Hydroxide; 2) Deposition rates of thin films; 3) Impregnation of porous electrodes; 4) Experimental Characterization of Nickel Hydroxide; 5) Diffusion coefficients of protons; 6) Self-discharge rates (i.e., oxygen-evolution kinetics); 7) Hysteresis between charge and discharge; 8) Capacity loss on cycling; 9) Experimental Verification of the Ni/H2 Battery Model; 10) Mathematical Modeling Li-Ion Batteries; 11) Experimental Verification of the Li-Ion Battery Model; 11) Integrated Power System Models for Satellites; and 12) Experimental Verification of Integrated-Systems Model.

Weidner, John W.↗

DOE OSTI2023

Evaluating the Feasibility of Using Strain Measured During Sinusoidal Rate Pumping Tests to Characterize an Aquifer

Pumping tests with sinusoidal variation in pumping rate have been proposed as a method for improving aquifer characterization. These tests can interrogate a larger aquifer volume than slug tests and they can be more sensitive to small variations in drawdown. Current methods of using sinusoidal variations of rate are based on measuring pressure signals from the reservoir or aquifer, which requires access to monitoring wells. An alternative approach has been developed that measures the strain in the vadose zone instead of pressure in the reservoir. An instrument has been developed at Clemson University that can measure small strains using optical fiber sensors, and we propose that it can be used to measure strain during sinusoidal pumping tests. The objective of this project is to evaluate the feasibility of measuring and interpreting the vertical strain in the vadose zone generated during a sinusoidal rate pumping tests to estimate aquifer and/or overburden properties. Pumping tests were performed at a site with a water table at 9 m and saturated saprolite from 9 to 30m depth. Sinusoidal pumping tests in the saprolite used different pumping amplitudes and periods while pressure was measured at three piezometers and vertical strain was measured in the vadose zone from 1 to 7 m depth. The delay time of the pressure increased linearly with distance, and it increased with the square root of the period, while the amplitude of the pressure decreased with distance and increased with amplitude. Vertical strain in the vadose zone varied with the same period as the pumping rate and the delay times were similar to those for pressures measured at similar radial distances. Delay times for the strain increased with an increasing propagation of the strain field. The hydraulic diffusivity was estimated using delay times for pressure and strain, and the results were compared with the values estimated from a constant-rate pumping test. The hydraulic diffusivity from the constant rate test was 0.039 m2/sec, whereas it was 0.053 (±0.05) m2/sec and 0.035 (±0.03) m2/sec from the pressure data and 0.08 (±0.022) m2/sec and 0.035 (±0.001) m2/sec from the strain data from sinusoidal tests of various periods and amplitudes. There is no statistical difference among these values. This study indicates that strain measurements in the vadose zone can be used to estimate aquifer properties during sinusoidal pumping tests.

Blais, Riley↗

NASA NTRS1995

Momentum deposition on Wolf-Rayet winds: Nonisotropic diffusion with effective gray opacity

We derive the velocity and mass-loss rate of a steady state Wolf-Rayet (WR) wind, using a nonisotropic diffusion approximation applied to the transfer between strongly overlapping spectral lines. Following the approach of Friend & Castor (1983), the line list is assumed to approximate a statistically parameterized Poisson distribution in frequency, so that photon transport is controlled by an angle-dependent, effectively gray opacity. We show the nonisotropic diffusion approximation yields good agreement with more accurate numerical treatments of the radiative transfer, while providing analytic insight into wind driving by multiple scattering. We illustrate, in particular, that multiple radiative momentum deposition does not require that potons be repeatedly reflected across substantial distances within the spherical envelope, but indeed is greatest when photons undergo a nearly local diffusion, e.g., through scattering by many lines closely spaced in frequency. Our results reiterate the view that the so-called 'momentum problem' of Wolf-Rayet winds is better characterized as an 'opacity problem' of simply identfying enough lines. One way of increasing the number of thick lines in Wolf-Rayet winds is to transfer opacity from saturated to unsaturated lines, yielding a steeper opacity distribution than that found in OB winds. We discuss the implications of this perspective for extending our approach to W-R wind models that incorporate a more fundamental treatment of the ionization and excitation processes that determine the line opacity. In particular, we argue that developing statistical descriptions of the lines to allow an improved effective opacity for the line ensemble would offer several advantages for deriving such more fundamental W-R wind models.

Gayley, Kenneth G.↗

NASA NTRS1989

PAHs molecules and heating of the interstellar gas

Until now it has remained difficult to account for the rather high temperatures seen in many diffuse interstellar clouds. Various heating mechanisms have been considered: photoionization of minor species, ionization of H by cosmic rays, and photoelectric effect on small grains. Yet all these processes are either too weak or efficient under too restricting conditions to balance the observed cooling rates. A major heat source is thus still missing in the thermal balance of the diffuse gas. Using photoionization cross sections measured in the lab, it was shown that in order to balance the observed cooling rates in cold diffuse clouds (T approx. 80 K) the PAHs would have to contain 15 percent of the cosmic abundance of carbon. This value does not contradict the former estimation of 6 percent deduced from the IR emission bands since this latter is to be taken as a lower limit. Further, it was estimated that the contribution to the heating rate due to PAH's in a warm HI cloud, assuming the same PAH abundance as for a cold HI cloud, would represent a significant fraction of the value required to keep the medium in thermal balance. Thus, photoionization of PAHs might well be a major heat source for the cold and warm HI media.

Verstraete, Laurent↗

DOE OSTI2024

HTO and selenate diffusion through compacted Na-, Na–Ca-, and Ca-montmorillonite

Radionuclide transport in smectite clay barrier systems used for nuclear waste disposal is controlled by diffusion, with adsorption significantly retarding transport rates. While a relatively minor component of spent nuclear fuel, 79 Se is a major driver of the safety case for spent fuel disposal due to its long half-life (3.3×10 5 yr) and its low adsorption to clay (K D < 10 L/kg), thus a thorough understanding of Se diffusion through clay is critical for understanding the long-term safety of spent fuel disposal systems. Through-diffusion experiments with tritiated water (HTO, conservative tracer) and Se(VI) were conducted with a well-characterized, purified montmorillonite source clay (SWy-2) under a constant ionic strength (0.1 M) and three different electrolyte compositions: Na + , Ca 2+ , and a Na + -Ca 2+ mixture at pH 6.5 in order to probe the effects of electrolyte composition and interlayer cation composition on clay microstructure, Se(VI) aqueous speciation, and ultimately diffusion. Further, the results were modeled using a reactive transport modeling approach to determine values of porosity (ε), D e (effective diffusion coefficient), and K D (distribution coefficient for adsorption). HTO diffusive flux was higher in Ca-montmorillonite (D e =1.68×10 -10 m 2 s -1 ) compared to Na-montmorillonite (De=7.83×10 -11 m 2 s -1 ). This increase in flux is likely due to a greater degree of clay layer stacking in the presence of Ca 2+ compared to Na + , which leads to larger inter-particle pores. Overall, the Se(VI) flux was much lower than the HTO flux due to anion exclusion, with Se(VI) flux following the order Ca (D e = 1.03×10 -11 m 2 s -1 ) > Na–Ca (D e = 2.12×10 -12 m 2 s -1 ) > Na (D e = 1.28×10 -12 m 2 s -1 ). These differences in Se(VI) flux are due to a combination of factors, including (1) larger accessible porosity in Ca-montmorillonite due to clay layer stacking and smaller electrostatic effects compared to Na-montmorillonite, (2) larger accessible porosity for neutral-charge CaSeO4 species which makes up 32% of aqueous Se(VI) in the pure Ca system, and (3) possibly higher Se(VI) adsorption for Ca-montmorillonite. Through a combination of experimental and modeling work, this study highlights the compounding effects that electrolyte and counterion compositions can have on radionuclide transport through clay. Diffusion models that neglect these effects are not transferable from laboratory experimental conditions to in situ repository conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗

NASA NTRS1978

Subophitic basalts from Mare Crisium - Cooling rates

Subophitic basalt is the most common rock type at Mare Crisium. The cooling rate of a sample of this rock was determined by: (1) an olivine cooling speedometer and (2) Zr partitioning between ilmenite and ulvospinel. The kinetic modeling of the olivine cooling speedometer starts with a calculation of the compositional profile of olivine (the 'as-solidified profile') and proceeds to re-equilibration by diffusion as a function of cooling rate. The estimated cooling rate for subophitic basalts from the Luna 24 site is in the range of 2 C/day (about 0.1 C/hr), which is reasonably well corroborated by dynamic crystallization studies of Grove (1978).

Taylor, L. A.↗

NASA NTRS1967

A Study to Analyze the Permeation of High Density Gases and Propellant Vapors Through Single Layer Teflon or Teflon Structure Materials and Laminations

This report contains the results of a fifteen month analytical and experimental study of the leakage rate of the pressurant gases (N2, He) and the propellant vapors (N2O4,N2H4) through bladder structures consisting of two layers of Teflon separated by a metallic foil diffusion barrier containing microscopic or larger holes. Results were obtained for the steady state leakage rate through circular holes and long rectangular openings in the barrier for arbitrary thicknesses of the two Teflon layers. The effect of hole shape and relative hole position on the leakage rate were studied. The transient problem was analyzed and it was shown that steady state calculations are adequate for estimating the leakage rate. A computer program entitled "Diffusion Analyzer Program" was developed to calculate the leakage rate, both transient and steady state. Finally, the analytical results were compared to experimentally determined values of the leakage rate through a model laminated bladder structure. The results of the analysis are in good agreement with experiment. The experimental effort (Part II of the Bladder Permeation Program) measured the solubility, diffusion coefficient and permeability of helium, nitrogen and nitrogen tetroxide vapor through Teflon TFE and FEP membranes. Data were obtained in the temperature range of 25 to 100 C at pressures ranging from near vacuum to about 20 atmospheres. Results of the experimental effort were compared with the limited data previously reported. As a verification to the applicability of results to actual bladder systems, counter diffusion tests were performed with a laminated sample containing aluminum foil with a selected group of holes.

Johnson, R. L.↗

DOE OSTI2022

Towards extreme fast charging of 4.6V LiCoO 2 via mitigating high-voltage kinetic hindrance

High-voltage LiCoO 2 (LCO) is an attractive cathode for ultra-high energy density lithium-ion batteries (LIBs) in the 3C markets. However, the sluggish lithium-ion diffusion at high voltage significantly hampers its rate capability. Herein, combining experiments with density functional theory (DFT) calculations, we demonstrate that the kinetic limitations can be mitigated by a facial Mg 2+ +Gd 3+ co-doping method. The as-prepared LCO shows significantly enhanced Li-ion diffusion mobility at high voltage, making more homogenous Li-ion de/intercalation at a high-rate charge/discharge process. The homogeneity enables the structural stability of LCO at a high-rate current density, inhibiting stress accumulation and irreversible phase transition. When used in combination with a Li metal anode, the doped LCO shows an extreme fast charging (XFC) capability, with a superior high capacity of 193.1 mAhg -1 even at the current density of 20C and high-rate capacity retention of 91.3% after 100 cycles at 5C. In conclusion, this work provides a new insight to prepare XFC high-voltage LCO cathode materials.

25 ENERGY STORAGE↗

DOE OSTI2021

Enhancing high-order harmonic generation by controlling the diffusion of the electron wave packet

We experimentally study the enhancement of high-order harmonic generation (HHG) driven by synthesized ω <#comment/> − <#comment/> 3 ω <#comment/> laser fields, where we control whether the ionization rate or the electron wave packet’s diffusion is the dominant enhancement mechanism. When minimizing the electron wave packet’s diffusion, the excursion times of the corresponding electron trajectories are reduced by a factor of 2 or more. This result is important for imaging techniques that use the returning electron wave packet to probe the remaining ion. Furthermore, we achieve a 10 × <#comment/> to 3800 × <#comment/> enhancement of the harmonic yield driven by the bichromatic fields relative to that of an optimized single-color field, showing that the bichromatic fields improve HHG’s capability as a light source. We also measure that the two-color field’s harmonics have half the divergence angle compared to their single-color counterpart, suggesting that the “short” electron trajectories play a more prominent role compared to their “long” trajectory counterparts, thus improving the wavefront of the emerging harmonic beam.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1979

Rates of mass, momentum, and energy transfer at the magnetopause

Empirical estimates of the global rates of transfer of solar wind mass, tangential momentum, and energy at the Earth's magnetopause are presented for comparison against model estimates based on the four principal mechanisms that have been proposed to explain such transfer. The comparisons, although not quite conclusive, strongly favor a model that incorporates some combination of direct magnetic connection and anomalous cross field diffusion. An additional global constraint, the rate at which magnetic flux is cycled through the magnetospheric convection system, strongly suggests that direct magnetic connection plays a significant if not dominant role in the solar wind/magnetosphere interaction.

Hill, T. W.↗