Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion couples”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Integrated geological, economic, and risk assessment of underground hydrogen storage

Underground hydrogen storage (UHS) is a promising option to buffer variable renewable power and support the hydrogen economy. Yet this technology is early-stage, and key uncertainties remain about how coupled processes such as geochemical and geomechanical impacts affect long-term storage security. In this paper, we seek to integrate the results from a large, multi-scale research program to quantify core feasibility metrics and situate these studies within ongoing efforts. Through this work, we evaluate hydrogen recoverability during storage by assessing diffusive losses, losses to geochemical interactions and their subsequent impact on geomechanical properties, and losses to residual saturation during injection and withdrawal. Molecular and laboratory studies show that lithologic heterogeneity, pore geometry, and mineral surface chemistry govern hydrogen–rock interactions, controlling loss pathways. High-pressure coreflooding is used to estimate relative permeabilities and measure parameters needed for reservoir-scale sensitivity analyses. These core-scale experiments found early hydrogen breakthrough at low saturations driven by capillary and viscous fingering, which limits pore-space utilization at early times, while repeated injection and withdrawal cycles improve deliverability over time. Technoeconomic factors and purity requirements for end uses are also considered as part of our comprehensive feasibility analysis. Finally, we underscore the importance of geophysical monitoring and tailored injection strategies to maintain integrity and efficiency. Together, these results provide a quantitative foundation for safe, scalable deployment of UHS and highlight priorities for future integrated studies.

08 HYDROGEN↗
DOE OSTI2023

Local inversion-symmetry breaking in a bismuthate high-Tc superconductor

Abstract The doped perovskite BaBiO 3 exhibits a maximum superconducting transition temperature ( T c ) of 34 K and was the first high- T c oxide to be discovered, yet pivotal questions regarding the nature of both the metallic and superconducting states remain unresolved. Although it is generally thought that superconductivity in the bismuthates is of the conventional s -wave type, the pairing mechanism is still debated, with strong electron-phonon coupling and bismuth valence or bond disproportionation possibly playing a role. Here we use diffuse x-ray scattering and Monte Carlo modeling to study the local structure of Ba 1-x K x BiO 3 across its insulator-metal boundary. We find no evidence for either long- or short-range disproportionation, which resolves a major conundrum, as disproportionation and the related polaronic effects are likely not relevant for the metallic and superconducting states. Instead, we uncover nanoscale structural correlations that break inversion symmetry, with far-reaching implications for the electronic physics. This unexpected finding furthermore establishes that the bismuthates belong to the broader classes of materials with hidden spin-orbit coupling and a tendency towards inversion-breaking displacements.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Simulation study of the effects of polymer network dynamics and mesh confinement on the diffusion and structural relaxation of penetrants

The diffusion of small molecular penetrants through polymeric materials represents an important fundamental problem, relevant to the design of materials for applications such as coatings and membranes. Polymer networks hold promise in these applications because dramatic differences in molecular diffusion can result from subtle changes in the network structure. In this paper, we use molecular simulation to understand the role that cross-linked network polymers have in governing the molecular motion of penetrants. By considering the local, activated alpha relaxation time of the penetrant and its long-time diffusive dynamics, we can determine the relative importance of activated glassy dynamics on penetrants at the segmental scale vs entropic mesh confinement on penetrant diffusion. We vary several parameters, such as the cross-linking density, temperature, and penetrant size, to show that cross-links primarily affect molecular diffusion through the modification of the matrix glass transition, with local penetrant hopping at least partially coupled to the segmental relaxation of the polymer network. This coupling is very sensitive to the local activated segmental dynamics of the surrounding matrix, and we also show that penetrant transport is affected by dynamic heterogeneity at low temperatures. To contrast, only at high temperatures and for large penetrants or when the dynamic heterogeneity effect is weak, does the effect of mesh confinement become significant, even though penetrant diffusion more broadly empirically follows similar trends as established models of mesh confinement-based transport.

Chemistry↗
DOE OSTI2023

Pore connectivity influences mass transport in natural rocks: Pore structure, gas diffusion and batch sorption studies

For this work, six rocks (one granodiorite, one limestone, two chalks, one mudstone, and one dolostone) with different extents of heterogeneity at six different particle sizes (from 75 to 8000 μm) were studied to describe the effects of pore connectivity on mass transport. The methods applied were (i) porosity measurement of granular rocks, (ii) analyses of gas-phase diffusive transport in a bed of packed particles, along with a solid quartz method at these six particle sizes being developed to identify the contribution of intraparticle diffusion, and (iii) batch sorption tests of multiple ions (anions and cations) with subsequent analyses of inductively coupled plasma-mass spectrometry. Granular porosity measurement results reveal that with decreasing particle sizes, the effective porosities for the “heterogenous” group of rocks (Grimsel granodiorite and Edwards limestone) increase, whereas the porosities of another “homogeneous” group (two Israel chalk samples, Japan mudstone, and Wyoming dolostone) remain constant. Gas diffusion results show that the intraparticle gas diffusion coefficient among these two sample groups, varying in the magnitude of 10 -8 to 10 -6 m 2 /s, are not directly correlated to the porosity differences. Moreover, the batch sorption work displays a different affinity of rocks for various tracers. For Grimsel granodiorite, Japan mudstone, and Wyoming dolostone, the adsorption capacity of Sm 3+ and Eu 3+ increases as the particle size decreases. In general, this integrated research of grain size distribution, granular rock porosity, intraparticle diffusivity, and ionic sorption capacity gives insights into the pore connectivity effect on both physical and chemical transport behaviors for different lithologies and/or different particle sizes.

58 GEOSCIENCES↗
DOE OSTI2025

Developing the Science Basis for Understanding Polymer Encapsulant Degradation Mechanisms: DuraMAT 2.0 Final Project Report

Polymeric encapsulants are essential materials in photovoltaic modules, protecting sensitive electronics from the environment while providing mechanical integrity to the multilayered assembly. However, these polymeric materials are susceptible to degradation processes driven by the ingress of environmental species, ultraviolet radiation, thermal stresses, and mechanical loading. In this study, we employ a combined atomistic simulation and accelerated aging experimental approach to study the molecular-scale mechanisms of encapsulant degradation. Classical molecular dynamics simulations quantify the diffusion of environmental and degradation species through the polymer matrix, producing composition-specific diffusion coefficients. Reactive simulations characterize activation energy barriers and reaction rate constants for key chemical pathways. In parallel, thermal-desorption analyses coupled with mass spectrometry monitor the emergence and concentration profiles of degradation products under controlled stressor conditions. By integrating simulation and experiment, we establish quantitative correlations between polymer composition, species diffusivity, and chemical reactivity. We anticipate that these relations and quantitative values could serve as high-fidelity inputs to reaction-diffusion models, enabling physics-informed lifetime predictions and guiding the design of more durable encapsulant materials for solar energy applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Anharmonic phonons, thermal expansion, and nuclear quantum effects in Zn

The anharmonic behavior of phonons and thermal expansion of hexagonal zinc were studied from 15 to 690 K by inelastic neutron scattering (INS) and ab initio simulations. Phonon spectra were measured for Q-points covering the full Brillouin zone, giving the phonon density of states (DOS) and dispersions along high-symmetry directions. The dispersions were sharp at 15 K, but diffuse intensity was observed at energies above them. The dispersions broadened with temperature T and the diffuse intensity became stronger. This diffuse intensity appeared in all INS measurements and simulations, except for classical molecular dynamics at 15 K. The temperature-dependent effective potential (TDEP) method, which included the nuclear quantum effect from zero-point vibrational dynamics, was used to calculate the free energy and thermal expansion. For T < 100 K nuclear quantum effects were important for obtaining the correct negative thermal expansion, and path integral molecular dynamics (PIMD) was particularly effective for obtaining the negative thermal expansion in the basal plane. A Heisenberg-Langevin model for interacting phonons coupled to a thermal bath was able to reproduce the shape and intensity of the diffuse spectral features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Exploring interlayer coupling in the twisted bilayer PtTe 2

We have investigated interlayer interactions in the bilayer PtTe 2 system, which influence the electronic energy bands near the Fermi level. Our diffusion Monte Carlo (DMC) calculations for the high-symmetry bilayer stackings (AA, AB, AC) manifest distinct interlayer binding characteristics among the stacking modes by revealing significantly different interlayer separations depending on the stacking, which is critical to understanding the interlayer coupling of the twisted bilayers consisting of various local stacking arrangements. Furthermore, a comparison between the interlayer separations obtained from DMC and density functional theory (DFT) shows that meta-generalized gradient approximation (GGA)-based van der Waals–DFT results agree with DMC for different layer stackings, including twisted bilayers, but only the ground-state AA stacking matches well with GGA-based DFT predictions. This underscores the importance of accurate exchange-correlation potentials even for capturing the stacking-dependent interlayer binding properties. We further show that the variability in DFT-predicted interlayer separations is responsible for the large discrepancy of band structures in the 21 . 79 ∘ twisted bilayer PtTe 2 , affecting its classification as metallic or insulating. These results demonstrate the importance of obtaining a correct description of stacking-dependent interlayer coupling in modeling delicate bilayer systems at finite twists. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2020

Solution-Based, Anion-Doping of Li 4 Ti 5 O 12 Nanoflowers for Lithium-Ion Battery Applications

Solution–based, anionic doping represents a convenient strategy with which to improve upon the conductivity of candidate anode materials such as Li 4 Ti 5 O 12 (LTO). As such, novel synthetic hydrothermally–inspired protocols have primarily been devised herein, aimed at the large–scale production of unique halogen–doped, micron–scale, three–dimensional, hierarchical LTO flower–like motifs. Although fluorine (F) doping has been explored, the use of chlorine (Cl) dopants is the primary focus here. Several experimental variables, such as dopant amount, lithium hydroxide concentration, and titanium butoxide purity, were probed and perfected. Furthermore, the Cl doping process did not damage the intrinsic LTO morphology. The analysis, based on interpreting a compilation of SEM, XRD, XPS, and TEM–EDS results, was used to determine an optimized dopant concentration of Cl. Electrochemical tests demonstrated an increased capacity via cycling of 12 % for a Cl–doped sample as compared with pristine LTO. Moreover, the Cl–doped LTO sample described in this study exhibited the highest discharge capacity yet reported at an observed rate of 2C for this material at 143mAh g –1 . Overall, these data suggest that the Cl dopant likely enhances not only the ion transport capabilities, but also the overall electrical conductivity of our as–prepared structures. Furthermore, to help explain these favorable findings, theoretical DFT calculations were used to postulate that the electronic conductivity and Li diffusion were likely improved by the presence of increased Ti 3+ ion concentration coupled with widening of the Li migration channel.

25 ENERGY STORAGE↗
DOE OSTI2022

Effects of EMIC Wave-Driven Proton Precipitation on the Ionosphere

In this report we investigate the proton precipitation caused by electromagnetic ion cyclotron (EMIC) waves and its impact on the ionosphere with the self-consistent ring current-atmosphere interactions model (RAM-SCBE) coupled with an ionospheric particle transport model Global Airglow (GLOW). The EMIC wave diffusion process causes significant precipitation of tens of keV protons, particularly in the afternoon to midnight sector. These precipitating energetic protons further impact the ionosphere-thermosphere and contribute to the ionization in the E/F regions. The integrated auroral conductance is significantly enhanced in the dusk-to-midnight sector. Although the EMIC waves do not directly interact with the ring current electrons, after the EMIC wave scattering included in the model, remarkable changes are found in the global distribution of precipitating electron flux. This means that the addition of proton precipitation driven by EMIC waves results in feedback effects on the ring current electron dynamics through the circulation system. Validation is also conducted by comparing the simulated precipitation flux and ionospheric electron density with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

High-power Mg batteries enabled by heterogeneous enolization redox chemistry and weakly coordinating electrolytes

Magnesium batteries have long been pursued as potentially low-cost, high-energy and safe alternatives to Li-ion batteries. However, Mg 2+ interacts strongly with electrolyte solutions and cathode materials, leading to sluggish ion dissociation and diffusion, and consequently low power output. Here we report a heterogeneous enolization chemistry involving carbonyl reduction (C=O ↔ C–O-), which bypasses the dissociation and diffusion difficulties, enabling fast and reversible redox processes. This kinetically favoured cathode is coupled with a tailored, weakly coordinating boron cluster-based electrolyte that allows for dendrite-free Mg plating/stripping at a current density of 20 mA cm -2 . The combination affords a Mg battery that delivers a specific power of up to 30.4 kW kg -1 , nearly two orders of magnitude higher than that of state-of-the-art Mg batteries. Here, the cathode and electrolyte chemistries elucidated here propel the development of magnesium batteries and would accelerate the adoption of this low-cost and safe battery technology.

25 ENERGY STORAGE↗
DOE OSTI2020

Calculating RF current condensation with consistent ray-tracing and island heating

Finally, by exploiting the nonlinear amplification of the power deposition of RF waves, current condensation promises new pathways to the stabilization of magnetic islands. We present a numerical analysis of current condensation, coupling a geometrical optics treatment of wave propagation and damping to a thermal diffusion equation solver in the island. Taking into account the island geometry and relativistic damping, previous analytical theory can be made more precise and specific scenarios can be realistically predicted. With this more precise description, bifurcations and associated hysteresis effects could be obtained in an ITER-like scenario at realistic parameter values. Moreover, it is shown that dynamically varying the RF wave launching angles can lead to hysteresis and help to avoid the nonlinear shadowing effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Resolving the dynamic correlated disorder in KTa 1- x Nb x O 3

Understanding the complex temporal and spatial correlations of ions in disordered perovskite oxides is critical to rationalize their functional properties. Here, in this study, we provide new insights into the longstanding controversy regarding the off-centering of transition metal ions in the archetypal ferroelectric alloy KTa 1-x Nb x O 3 (KTN). By mapping the full energy (E) and wavevector (Q) dependence of the dynamical structure factor S(Q, E) using neutron scattering, and rationalizing our observations with atomistic simulations leveraging machine learning, we fully resolve the static vs dynamic nature of diffuse scattering sheets, as well as their composition (x) and temperature dependence. Our first-principles simulations, extended with machine-learning molecular dynamics, reproduce both inelastic neutron spectra and diffuse features, and establish how dynamically-correlated transition metal off-centerings couple to phonons, unifying local and collective viewpoints. This study sheds new light into an exemplary ferroelectric systems and shows the importance of mapping the full S(Q, E) to reveal critical spatio-temporal correlations of atomic disorder from which functional properties emerge.

42 ENGINEERING↗
DOE OSTI2023

Diagnosing Near-Surface Model Errors with Candidate Physics Parameterization Schemes for the Multiphysics Rapid Refresh Forecast System (RRFS) Ensemble during Winter over the Northeastern United States and Southern Great Plains

Abstract During the winter of 2020/21 an ensemble of FV3-LAM forecasts was produced over the contiguous United States for the Winter Weather Experiment using five physics suites. These forecasts are evaluated with the goal of optimizing physics parameterizations within the future operational Rapid Refresh Forecast System (RRFS) in the Unified Forecast System (UFS) realm and for selecting suitable physics suites for a multiphysics RRFS ensemble. The five physics suites have different combinations of land surface models (LSMs), planetary boundary layer (PBL) parameterizations, and surface layer schemes, chosen from those used in current and possible future operational systems and likely to be supported in the operational UFS. Full-season evaluation reveals a persistent near-surface cold bias in the U.S. Northeast from one suite and a nighttime warm bias in the southern Great Plains in another suite, while other suites have smaller biases. A representative case is chosen to diagnose the cause for each of these biases using sensitivity simulations with different physics combinations or modified parameters and verified with additional mesonet observations. The cold bias in the Northeast is attributed to aspects of the Noah-MP LSM over snow cover, where Noah-MP simulates lower soil water content, and thus lower thermal conductivity than other LSMs, leading to less upward ground heat flux during nighttime and consequently lower surface temperature. The nighttime warm bias found in the southern Great Plains is attributed to overestimation of vertical mixing in the K -profile-based eddy-diffusivity mass-flux (K-EDMF) PBL scheme and insufficient land–atmospheric coupling from the GFS surface layer scheme over short vegetation. A few key parameters driving these systematic biases are identified.

Meteorology & Atmospheric Sciences↗
DOE OSTI2020

Selective One-Dimensional 13 C– 13 C Spin-Diffusion Solid-State Nuclear Magnetic Resonance Methods to Probe Spatial Arrangements in Biopolymers Including Plant Cell Walls, Peptides, and Spider Silk

Two-dimensional (2D) and 3D through-space 13 C– 13 C homonuclear spin-diffusion techniques are powerful solid-state nuclear magnetic resonance (NMR) tools for extracting structural information from 13 C-enriched biomolecules, but necessarily long acquisition times restrict their applications. In this work, we explore the broad utility and underutilized power of a chemical shift-selective one-dimensional (1D) version of a 2D 13 C– 13 C spin-diffusion solid-state NMR technique. The method, which is called 1D dipolar-assisted rotational resonance (DARR) difference, is applied to a variety of biomaterials including lignocellulosic plant cell walls, microcrystalline peptide fMLF, and black widow dragline spider silk. 1D 13 C– 13 C spin-diffusion methods described here apply in select cases in which the 1D 13 C solid-state NMR spectrum displays chemical shift-resolved moieties. Furthermore, this is analogous to the selective 1D nuclear Overhauser effect spectroscopy (NOESY) experiment utilized in liquid-state NMR as a faster (1D instead of 2D) and often less ambiguous (direct sampling of the time domain data, coupled with increased signal averaging) alternative to 2D NOESY. Selective 1D 13 C– 13 C spin-diffusion methods are more time-efficient than their 2D counterparts such as proton-driven spin diffusion (PDSD) and dipolar-assisted rotational resonance. The additional time gained enables measurements of 13 C– 13 C spin-diffusion buildup curves and extraction of spin-diffusion time constants TSD, yielding detailed structural information. Specifically, selective 1D DARR difference buildup curves applied to 13 C-enriched hybrid poplar woody stems confirm strong spatial interaction between lignin and acetylated xylan polymers within poplar plant secondary cell walls, and an interpolymer distance of ~0.45–0.5 nm was estimated. Additionally, Tyr/Gly long-range correlations were observed on isotopically enriched black widow spider dragline silks.

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

The role of high-order electron correlation effects in a model system for non-valence correlation-bound anions

The diffusion Monte Carlo (DMC), auxiliary field quantum Monte Carlo (AFQMC), and equation-of-motion coupled cluster (EOM-CC) methods are used to calculate the electron binding energy (EBE) of the non-valence anion state of a model (H2O)4 cluster. Two geometries are considered, one at which the anion is unbound and the other at which it is bound in the Hartree–Fock (HF) approximation. It is demonstrated that DMC calculations can recover from the use of a HF trial wave function that has collapsed onto a discretized continuum solution, although larger EBEs are obtained when using a trial wave function for the anion that provides a more realistic description of the charge distribution and, hence, of the nodal surface. For the geometry at which the cluster has a non-valence correlation-bound anion, both the inclusion of triples in the EOM-CC method and the inclusion of supplemental diffuse d functions in the basis set are important. DMC calculations with suitable trial wave functions give EBE values in good agreement with our best estimate EOM-CC result. AFQMC using a trial wave function for the anion with a realistic electron density gives a value of the EBE nearly identical to the EOM-CC result when using the same basis set. For the geometry at which the anion is bound in the HF approximation, the inclusion of triple excitations in the EOM-CC calculations is much less important. The best estimate EOM-CC EBE value is in good agreement with the results of DMC calculations with appropriate trial wave functions.

Upadhyay, Shiv (ORCID:0000000285010501)↗
DOE OSTI2022

Evaluation of subsurface transport processes of delayed gas signatures applicable to underground nuclear explosions

Radioactive gas signatures from underground nuclear explosions (UNEs) result from gas-migration processes occurring in the subsurface. The processes considered in this study either drive or retard upward migration of gases from the detonation cavity. The relative importance of these processes is evaluated by simulating subsurface transport in a dual-permeability medium for the multi-tracer Noble Gas Migration Experiment (NGME) originally intended to study some aspects of transport from a UNE. For this experiment, relevant driving processes include weak two-phase convection driven by the geothermal gradient, over pressuring of the detonation cavity, and barometric pumping while gas sorption, dissolution, radioactive decay, and usually diffusion represent retarding processes. From deterministic simulations we found that over-pressuring of the post-detonation chimney coupled with barometric pumping produced a synergistic effect amplifying the tracer-gas reaching the surface. Bounding simulations indicated that the sorption and dissolution of gases, tending to retard transport, were much smaller than anticipated by earlier laboratory studies. The NGME observations themselves show that differences in gas diffusivity have a larger effect on influencing upward transport than do the combined effects of tracer-gas sorption and dissolution, which is consistent with a Sobol’ sensitivity analysis. Both deterministic simulations and those considering parametric uncertainties of transport-related properties predict that the excess in concentration of SF 6 compared to 127 Xe as might be captured in small volumetric samples should be much smaller than the order-of-magnitude contrast found in the large-volume gas samples taken at the site. While extraction of large-volume subsurface gas samples is shown to be capable of distorting in situ gas compositions, the highly variable injection rate of SF 6 into the detonation cavity relative to that of 127 Xe at the start of the field experiment is the most likely explanation for the large difference in observed concentrations.

58 GEOSCIENCES↗
DOE OSTI2023

The Effect of Monomer Size on Fusion and Coupling in Colloidal Quantum Dot Molecules

The fusion step in the formation of colloidal quantum dot molecules, constructed from two core/shell quantum dots, dictates the coupling strength and hence their properties and enriched functionalities compared to monomers. Herein, studying the monomer size effect on fusion and coupling, we observe a linear relation of the fusion temperature with the inverse nanocrystal radius. This trend, similar to that in nanocrystal melting, emphasizes the role of the surface energy. The suggested fusion mechanism involves intraparticle ripening where atoms diffuse to the reactive connecting neck region. Moreover, the effect of monomer size and neck filling on the degree of electronic coupling is studied by combined atomistic-pseudopotential calculations and optical measurements, uncovering strong coupling effects in small QD dimers, leading to significant optical changes. Understanding and controlling the fusion and hence coupling effect allows tailoring the optical properties of these nanoscale structures, with potential applications in photonic and quantum technologies.

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

Facile electrocatalytic proton reduction by a [Fe–Fe]-hydrogenase bio-inspired synthetic model bearing a terminal CN — ligand

An azadithiolate bridged CN — bound pentacarbonyl bis-iron complex, mimicking the active site of [Fe–Fe] H 2 ase is synthesized. The geometric and electronic structure of this complex is elucidated using a combination of EXAFS analysis, infrared and Mössbauer spectroscopy and DFT calculations. The electrochemical investigations show that complex 1 effectively reduces H + to H 2 between pH 0–3 at diffusion-controlled rates (10 11 M —1 s —1 ) i.e. 10 8 s —1 at pH 3 with an overpotential of 140 mV. Electrochemical analysis and DFT calculations suggests that a CN— ligand increases the p K a of the cluster enabling hydrogen production from its Fe(I)–Fe(0) state at pHs much higher and overpotential much lower than its precursor bis-iron hexacarbonyl model which is active in its Fe(0)–Fe(0) state. The formation of a terminal Fe–H species, evidenced by spectroelectrochemistry in organic solvent, via a rate determining proton coupled electron transfer step and protonation of the adjacent azadithiolate, lowers the kinetic barrier leading to diffusion controlled rates of H 2 evolution. The stereo-electronic factors enhance its catalytic rate by 3 order of magnitude relative to a bis-iron hexacarbonyl precursor at the same pH and potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗