Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion processes”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Linking pore-scale and basin-scale effects on diffusive methane transport in hydrate bearing environments through multi-scale reservoir simulations

We explore the gas hydrate-generating capacity of diffusive methane transport induced by solubility gradients due to pore size contrasts in lithologically heterogeneous marine sediments. Through the use of 1D, 2D, and 3D reactive transport simulations, we investigate scale-dependent processes in diffusion-dominated gas hydrate systems. These simulations all track a sand body, or series of sands, surrounded by clays as they are buried through the gas hydrate stability zone. Methane is sourced by microbial methanogenesis in the clays surrounding the sand layers. In 1D, simulations performed in a Lagrangian reference frame demonstrate that gas hydrate in thin sands (3.6 m thick) can occur in high saturations (upward of 70%) at the edges of sand bodies within the upper 400 meters below the seafloor. Diffusion of methane toward the center of the sand layer depends on the concentration gradient within the sand: broader sand pore size distributions with smaller median pore sizes enhance diffusive action toward the sand’s center. Incorporating downhole log- and laboratory-derived sand pore size distributions, gas hydrate saturations in the center of the sand can reach 20% of the hydrate saturations at the sand’s edges. Furthermore, we show that hydrate-free zones exist immediately above and below the sand and are approximately 5 m thick, depending on the sand-clay solubility contrast. A moving reference frame is also adopted in 2D, and the angle of gravity is rotated relative to the grid system to simulate a dipping sand layer. This is important to minimize diffusive edge effects or numerical diffusion that might be associated with a dipping sand in an Eulerian grid system oriented orthogonal to gravity. Two-dimensional simulations demonstrate the tendency for gas hydrate to accumulate downdip in a sand body because of greater methane transport at depth due to larger sand-clay solubility contrasts. In 3D, basin-scale simulations illuminate how convergent sand layers in a multilayered system can compete for diffusion from clays between them, resulting in relatively low hydrate saturations. All simulations suggest that when hydrate present in clays dissociates with burial, the additional dissolved methane is soaked up by nearby sands preserving high hydrate saturations.

Nole, Michael↗
DOE OSTI2024

Unraveling Interdiffusion Phenomena and the Role of Nanoscale Diffusion Barriers in the Copper–Gold System

Diffusion is one of the most fundamental concepts in materials science, playing a pivotal role in materials synthesis, forming, and degradation. Of particular importance is solid state interdiffusion of metals which defines the usable parameter space for material combinations in the form of alloys. This parameter space can be explored on the macroscopic scale by using diffusion couples. However, this method reaches its limit when going to low temperatures, small scales, and when testing ultrathin diffusion barriers. Therefore, this work transfers the principle of the diffusion couples to small scales by using core–shell nanowires and in situ heating. This allows us to delve into the interdiffusion dynamics of copper and gold, revealing the interplay between diffusion and the disorder–order phase transition. Our in situ TEM experiments in combination with chemical mapping reveal the interdiffusion coefficients of Cu and Au at low temperatures and highlight the impact of ordering processes on the diffusion behavior. The formation of ordered domains within the solid-solution is examined using high-resolution imaging and nanodiffraction including strain mapping. In addition, we examine the effectiveness of ultrathin Al 2 O 3 barrier layers to control interdiffusion of the diffusion couple. Our findings indicate that a 5 nm thick layer serves as an efficient diffusion barrier. Furthermore, this research provides valuable insights into the interdiffusion behavior of Cu and Au on the nanoscale, offering potential applications in the development of miniaturized integrated circuits and nanodevices.

alloys↗
DOE OSTI2015

Structure and ionic conductivity of well-aligned polycrystalline sodium titanogallate grown by reactive diffusion

We prepared the b-axis-oriented polycrystalline Na{sub 0.85}Ti{sub 0.51}Ga{sub 4.37}O{sub 8} (NTGO) embedded in Ga{sub 2}O{sub 3}-doped Na{sub 2}Ti{sub 4}O{sub 9} matrix using the reactive diffusion technique. When the sandwich-type Ga{sub 2}TiO{sub 5}/NaGaO{sub 2}/Ga{sub 2}TiO{sub 5} diffusion couple was heated at 1323 K for 24 h, the NTGO polycrystal was readily formed in the presence of a liquid phase. The resulting polycrystalline material was characterized by X-ray diffractometry, electron microscopy and impedance spectroscopy. We mechanically processed the annealed diffusion couple and obtained the thin-plate electrolyte consisting mostly of the grain-aligned NTGO polycrystal. The ionic conductivity (σ) of the electrolyte along the common b-axis direction steadily increased from 1.3×10{sup −4} to 7.3×10{sup −3} S/cm as the temperature increased from 573 to 1073 K. There was a slope change at ca. 792 K for the Arrhenius plot of σ; the activation energies were 0.39 eV above this temperature and 0.57 eV below it. The NTGO showed the crystal structure (space group C2/m) with substantial positional disordering of one of the two Ga sites. The Na{sup +} ions occupied ca. 43% of the Wyckoff position 4i site, the deficiency of which would contribute to the relatively high ionic conductivity along the b-axis. The reactive diffusion could be widely applicable as the novel technique to the preparation of grain-aligned ceramics of multi-component systems. - Graphical abstract: We have prepared the b-axis-oriented Na{sub 0.85}Ti{sub 0.51}Ga{sub 4.37}O{sub 8} polycrystal embedded in Ga{sub 2}O{sub 3}-doped Na{sub 2}Ti{sub 4}O{sub 9} matrix by the heat treatment of sandwich-type diffusion couple of Ga{sub 2}TiO{sub 5}/NaGaO{sub 2}/Ga{sub 2}TiO{sub 5}. The resulting Na{sub 0.85}Ti{sub 0.51}Ga{sub 4.37}O{sub 8} electrolyte showed the ionic conductivity ranging from 1.3×10{sup −4} S/cm at 573 K to 7.3×10{sup −3} S/cm at 1073 K. - Highlights: • The b-axis-oriented polycrystalline Na{sub 0.85}Ti{sub 0.51}Ga{sub 4.37}O{sub 8} is successfully prepared. • Crystal structure of Na{sub 0.85}Ti{sub 0.51}Ga{sub 4.37}O{sub 8} is determined by single-crystal XRD. • The polycrystal shows relatively high Na{sup +} ion conductivity along the common b-axis. • Reactive diffusion is successfully used for the preparation of grain-aligned ceramics.

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

Laboratory Experiments to Evaluate Matrix Diffusion of Dissolved Organic Carbon Carbon-14 in Southern Nevada Fractured-rock Aquifers

Dissolved inorganic carbon (DIC) carbon-14 ( 14 C) is used to estimate groundwater ages by comparing the DIC 14 C content in groundwater in the recharge area to the DIC 14 C content in the downgradient sampling point. However, because of chemical reactions and physical processes between groundwater and aquifer rocks, the amount of DIC 14 C in groundwater can change and result in 14 C loss that is not because of radioactive decay. This loss of DIC 14 C results in groundwater ages that are older than the actual groundwater ages. Alternatively, dissolved organic carbon (DOC) 14 C in groundwater does not react chemically with aquifer rocks, so DOC 14 C ages are generally younger than DIC 14 C ages. In addition to chemical reactions, 14 C ages may also be altered by the physical process of matrix diffusion. The net effect of a continuous loss of 14 C to the aquifer matrix by matrix diffusion and then radioactive decay is that groundwater appears to be older than it actually is. Laboratory experiments were conducted to measure matrix diffusion coefficients for DOC 14 C in volcanic and carbonate aquifer rocks from southern Nevada. Experiments were conducted using bromide (Br-) as a conservative tracer and 14 C-labeled trimesic acid (TMA) as a surrogate for groundwater DOC. Outcrop samples from six volcanic aquifers and five carbonate aquifers in southern Nevada were used. The average DOC 14 C matrix diffusion coefficient for volcanic rocks was 2.9 x 10 -7 cm 2 /s, whereas the average for carbonate rocks was approximately the same at 1.7 x 10 -7 cm 2 /s. The average Br- matrix diffusion coefficient for volcanic rocks was 10.4 x 10 -7 cm 2 /s, whereas the average for carbonate rocks was less at 6.5 x 10 -7 cm 2 /s. Carbonate rocks exhibited greater variability in DOC 14 C and Br- matrix diffusion coefficients than volcanic rocks. These results confirmed, at the laboratory scale, that the diffusion of DOC 14 C into southern Nevada volcanic and carbonate aquifers is slower than DIC 14 C. Because of the apparent sorption of 14 C-labeled TMA in the experiments, matrix diffusion coefficients are likely even lower. The reasons for the higher than expected Br-/ 14 C-labeled TMA are unknown. Because the molecular size of TMA is on the low end of the range in molecular size for typical humic substances, the matrix diffusion coefficients for the 14 C-labeled TMA likely represent close to the maximum diffusion rates for DOC 14 C in the volcanic and carbonate aquifers in southern Nevada.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Hydrogen uptake in graphite matrix at high temperature

Tritium management is a critical challenge for the next generation of nuclear reactors, such as Fluoride Salt Cooled High Temperature Reactors (FHRs) and High Temperature Gas-cooled Reactors (HTGRs), due to the higher production rate (up to 10,000 times) than conventional Light Water Reactors (LWRs). Graphitic materials employed as moderator, reflector, and fuel pebbles offer a potential pathway for tritium recovery by serving as a sink for tritium. Prediction of uptake capacity under reactor relevant conditions remains a challenge due to a lack of low partial pressure data and significant inter-grade variability of graphite. This study addresses these gaps by providing a comprehensive characterization of hydrogen (as a tritium surrogate) uptake and release behavior in the A3-3 graphite matrix (GM) used in fuel pebbles. Uptake measurements are performed at reactor relevant temperatures of 600- 800 °C, 1-200 Torr hydrogen pressure, and 15-120 min equilibration time, followed by thermal desorption spectroscopy up to 1100 °C. Uptake experiments at different equilibration times demonstrate the role of kinetics in hydrogen uptake, which can be modeled as a diffusion-with-trapping process. In the thermodynamics limit, the Sips adsorption model is shown to capture the uptake in A3-3 GM well. Our campaign provides a set of new results for hydrogen uptake in A3-3, including limiting uptake capacity at 600 °C, apparent diffusion coefficient at 600 °C, and the first estimates of the FHR/HTGR relevant (600 °C, 20 Pa partial pressure) equilibrium uptake capacity and time to saturation. Desorption data highlights a new site for hydrogen uptake, not observed in nuclear graphite, which we attribute to the non-graphitized binder. Using the Kissinger method, we estimate activation energy for release from the desorption peaks, confirming the activation energy for release from the basal planes and providing the first estimate for the activation energy of release from the binder.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

GRIZZLY MODEL OF FULLY COUPLED HEAT TRANSFER, MOISTURE DIFFUSION, ALKALI-SILICA REACTION AND FRACTURING PROCESSES IN CONCRETE

Idaho National Laboratory (INL) has been developing a fully coupled thermo-hydro-mechanical-chemical (THMC) concrete model, referred as Grizzly, with the ultimate goal of reliably simulating and predicting the long-term performance of concrete structures subjected to various aging and degradation mechanisms. In this paper, we present a fully coupled concrete model implemented in our Grizzly code for simulating coupled moisture/heat transfer, stress-dependent ASR swelling and elastic/inelastic mechanics. Example simulation results and preliminary model validation results against experimental data reported in the literature will be provided also in this paper. The close matches between the simulation results and experiments clearly reveal the potential of the fully coupled solution strategy for solving tightly coupled THMC concrete models. We also present simulation results of using a rate-dependent damaging mechanics model for simulating the ASR swelling induced degradation of concrete structures and its feedback to moisture and heat transport processes.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2019

Degradation of Mg-doped zinc oxide buffer layers in thin film CdTe solar cells

Cadmium Sulphide is the conventional n-type buffer layer used in thin film Cadmium Telluride solar cells. It is well known that Cadmium Sulphide causes optical losses and sulphur diffuses into the absorber during high temperature activation. Sputter-deposited Mg-doped ZnO (MZO) has been shown to be an attractive buffer layer for Cadmium Telluride solar cells due to its transparency and tuneable band gap. It is also stable to high temperature processing and avoids diffusion of elements into the cadmium telluride absorber during the cadmium chloride activation treatment. However, degradation is observed in solar cells incorporating MZO buffer layers. Analysis of the MZO film surface potential has revealed significant fluctuations in the thin film work function once the layer is exposed to the atmosphere following deposition. These fluctuations are attributed to the high reactivity to water vapour of the MgO contained in the MZO films. This has been analysed using X-ray Photoelectron Spectroscopy to determine corresponding changes in the surface chemistry. The Zinc Oxide component is relatively stable, but the analysis shows that MgO forms a Mg(OH) 2 layer on the MZO surface which forms a secondary barrier at the MZO/CdTe interface and/or at the interface between MZO and the Fluorine-doped SnO 2 . This affects the Fill Factor and as a consequence it degrades the conversion efficiency.

14 SOLAR ENERGY↗
DOE OSTI2014

Model of evolution of surface grain structure under ion bombardment

Diffusion and chemical reactions in multicomponent systems play an important role in numerous technology applications. For example, surface treatment of materials and coatings by particle beam leads to chemical composition and grain structure change. To investigate the thermal-diffusion and chemical processes affecting the evolution of surface structure, the mathematical modeling is efficient addition to experiment. In this paper two-dimensional model is discussed to describe the evolution of titanium nitride coating on the iron substrate under implantation of boron and carbon. The equation for diffusion fluxes and reaction rate are obtained using Gibbs energy expansion into series with respect to concentration and their gradients.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Computing Nonequilibrium Responses with Score-Shifted Stochastic Differential Equations

Using equilibrium fluctuations to understand the response of a physical system to an externally imposed perturbation is the basis for linear response theory, which is widely used to interpret experiments and shed light on microscopic dynamics. For nonequilibrium systems, perturbations cannot be interpreted simply by monitoring fluctuations in a conjugate observable and general response results rely on path ensemble averaging. Furthermore, these techniques do not apply to perturbations that affect the diffusion tensor in a stochastic system. Here, we introduce an “effective” physical process that represents the diffusion perturbed dynamics and enables accurate calculations of responses to a change in the diffusion. Interestingly, the effective dynamics contain an additional drift involving the instantaneous “score” of the system, and we leverage score matching algorithms to carry out nonequilibrium response calculations on systems for which the exact stationary distribution is unknown.

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

Voltammetric Study of Antitumor Efficient Ethylenediamine-Type of Ligands

The electrochemical behavior of four dialkyl esters (methyl, ethyl, n-propyl and n-butyl) of (S,S)-α,α ′-(1,2-ethanediyldiimino)- bis[cyclohexanepropanoic acid] with moderate antitumor activity has been studied by cyclic voltammetry. I-E curves recorded at boron-doped diamond electrode in -2.00 < Ep < 0.00 V range show a single oxidative peak located in the -0.95 < Ep,a < -0.52 V region in a form of one-electron wave. The electrochemical processes appear to be irreversible for all compounds except for the compound with methyl substituent which demonstrates well-defined reversible electron transfer. Derivatives with bulkier moieties are easier to oxidize. For all investigated compounds an amine-centered redox process is dominant. Further investigation indicates that electrochemical reaction is diffusion controlled, chemically irreversible process. The results are discussed in terms of linking their electrochemical activity to antitumor effect. © 2018 The Electrochemical Society.

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

Quantification of the Crack Evolution Process by Extracting Relevant Signal Components from Wave Propagation and Diffusive Transport Front Measurements

Wave propagation and diffusive transport phenomena in a geological rock sample undergoing crack evolution process are expected to interact with the mechanical discontinuities in the medium. The measurements of the signals associated with these phenomena can be used to assess and monitor the crack-driven micromechanical alterations in the rock. Different wave/diffusion phenomena, such as sonic propagation, pressure diffusion, and acoustic emission (AE), are sensitive to different elements of the mechanical discontinuities generated during the evolution of the crack clusters from initiation to coalescence. Sonic propagation, AE, and pressure diffusion monitoring have the potential to map the crack evolution because the transmitter-receiver arrays can be designed, arranged and tuned to (1) achieve maximum recovery of the scattered waveforms and travel times, (2) capture the later arrivals and multiple reflections, and (3) illuminate large rock volume. However, the structural/topological complexities of the mechanical discontinuities, complex distribution of the stress fields, complex mechanical alterations in media, and fluid redistribution in the crack system pose serious challenges for the detection and modeling of the crack evolution process (from here on, we will use the term ‘crack evolution process’ to mean that the crack evolution occurred under shallow crustal conditions). For purposes of accurately accounting such complexities and heterogeneities in the absence of reliable physical laws, simulation methods, and signal processing techniques, my early-career research proposal will develop and apply novel data-driven machine learning methods to: (1) extract signal components relevant to the various phases of crack evolution and (2) generate a 2D visual map of the crack evolution process.

58 GEOSCIENCES↗
DOE OSTI2021

Kinetics of solid–liquid interface motion in molecular dynamics and phase-field models: crystallization of chromium and silicon

The results of molecular dynamics (MD) simulations of the crystallization process in one-component materials and solid solution alloys reveal a complex temperature dependence of the velocity of the crystal–liquid interface featuring an increase up to a maximum at 10–30% undercooling below the equilibrium melting temperature followed by a gradual decrease of the velocity at deeper levels of undercooling. At the qualitative level, such non-monotonous behaviour of the crystallization front velocity is consistent with the diffusion-controlled crystallization process described by the Wilson–Frenkel model, where the almost linear increase of the interface velocity in the vicinity of melting temperature is defined by the growth of the thermodynamic driving force for the phase transformation, while the decrease in atomic mobility with further increase of the undercooling drives the velocity through the maximum and into a gradual decrease at lower temperatures. At the quantitative level, however, the diffusional model fails to describe the results of MD simulations in the whole range of temperatures with a single set of parameters for some of the model materials. The limited ability of the existing theoretical models to adequately describe the MD results is illustrated in the present work for two materials, chromium and silicon. It is also demonstrated that the MD results can be well described by the solution following from the hodograph equation, previously found from the kinetic phase-field model (kinetic PFM) in the sharp interface limit. The ability of the hodograph equation to describe the predictions of MD simulation in the whole range of temperatures is related to the introduction of slow (phase field) and fast (gradient flow) variables into the original kinetic PFM from which the hodograph equation is obtained. The slow phase-field variable is responsible for the description of data at small undercoolings and the fast gradient flow variable accounts for local non-equilibrium effects at high undercoolings. The introduction of these two types of variables makes the solution of the hodograph equation sufficiently flexible for a reliable description of all nonlinearities of the kinetic curves predicted in MD simulations of Cr and Si. This article is part of the theme issue ‘Transport phenomena in complex systems (part 1)’.

Science & Technology - Other Topics↗
DOE OSTI2023

Rubidium and potassium isotopic variations in chondrites and Mars: Accretion signatures and planetary overprints

As moderately volatile elements, isotopes of Rb and K can trace volatilization processes in planetary bodies. Rubidium isotopic data are however very scarce, especially for non-carbonaceous meteorites. Here, in this study, we report combined Rb and K isotopic data (δ 87/85 Rb and δ 41/39 Κ) for 7 ordinary, 6 enstatite, and 4 Martian meteorite falls to understand the causes for the variations in volatile abundances and isotopic compositions. Bulk Rb and K isotopic compositions of planetary bodies are estimated to be (Table 1): Mars +0.10 ± 0.03 ‰ for Rb and -0.26 ± 0.05 ‰ for K, bulk OCs $-0.12^{+0.15}_{-0.24}$ ‰ for Rb and $-0.72^{+0.28}_{-0.41}$ ‰ for K, bulk ECs $-0.02^{+0.29}_{-0.26}$ ‰ for Rb and $-0.33^{+0.37}_{-0.23}$ ‰ for K. The bulk K isotopic compositions of subgroup OCs are estimated to be $-0.72^{+0.26}_{-0.55}$ ‰ for H chondrites, $-0.71^{+0.23}_{-0.39}$ ‰ for L chondrites, and $-0.77^{+0.63}_{-0.30}$ ‰ for LL chondrites. A broad correlation between the Rb and K isotopic compositions of planetary bodies is observed. The correlation follows a slope that is consistent with kinetic evaporation and condensation processes, suggesting volatility-controlled mass-dependent isotope fractionation (as opposed to nucleosynthetic anomalies). Individual ordinary and enstatite chondrites show large Rb and K isotopic variations (-1.02 to +0.29 ‰ for Rb and -0.91 to -0.15 ‰ for K). Samples of lower metamorphic grades display correlated elemental and isotopic fractionations between Rb and K, while samples of higher metamorphic grades show great scatter, suggesting that chondrite parent-body processes have decoupled the two elements and their isotopes at the sample scale. Several processes could have contributed to the observed isotopic variations of Rb and K, including (i) chondrule “nugget effect”, (ii) volatilization during parent-body thermal metamorphism (heat-induced vaporization and gas transport within parent bodies), (iii) thermal diffusion during parent-body metamorphism, and (iv) impact/shock heating. Quantitative modeling of the first two processes suggests that neither of them could produce isotopic variations large enough to explain the observed isotopic variations. Volatilization during parent-body thermal metamorphism [the scenario (ii)], which has been commonly invoked to explain the isotopic variations of volatile elements, is gas transport-limited and its effect on isotopic fractionations of moderately volatile elements should be negligible. Modeling of diffusion processes suggests that (iii) could produce K isotopic variation comparable to the observed variation. The large isotopic variations in non-carbonaceous meteorites are thus most likely due to diffusive redistribution of K and Rb during metamorphism and/or shock-induced heating and vaporization.

58 GEOSCIENCES↗
DOE OSTI2016

Experimental and Computational Investigation of Lepidocrocite Anodes for Sodium-Ion Batteries

In this work, we investigated several titanates with lepidocrocite-type structures (general formula A x Ti1- y M y O 4 , with A = Na and M = Li or Mg), having potential utility as anode materials for sodium-ion batteries. First-principles calculations were used to determine key battery metrics, including potential profiles, structural changes during sodiation, and sodium diffusion energy barriers for several compositions, and were compared to experimental results. Site limitations were found to be critical determinants of the gravimetric capacities, which are also affected both by the stacking arrangement of the corrugated layers and the identity of M (Li or Mg). To explain the experimentally observed lattice parameter changes observed as a function of the state of charge, it was necessary to assume the participation of water/solvent during the sodium intercalation process. Here, sodium diffusion barriers were also found to vary as a function of state of charge and diffusion direction, with a spread of 0.06-1.3 eV at low sodium contents, narrowing to 0.3-0.5 eV at higher sodium contents. Based on these results, strategies for selecting and improving the performance of these electrode materials are suggested.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Live-monitoring of Te inclusions laser-induced thermo-diffusion and annealing in CdZnTe crystals

The presence of Te inclusions is one of the main factors limiting performances of CdZnTe crystals as X-ray detectors. We show that by means of infrared laser radiation it is possible to move and anneal tellurium inclusions exploiting a thermo-diffusion mechanism. The process is studied live during irradiation by means of an optical microscope equipment. Experimental conditions, and, in particular, energy laser fluence, for annealing inclusions of different dimensions are determined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗