Search NASA⌕ Search

SEARCH · Search NASA

Results for “DIFFUSER”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

MODIS Solar Diffuser Attenuation Screen Modeling Results

On-orbit calibration of the reflected solar bands on the EOS Moderate Resolution Imaging Spectroradiometer (MODIS) is accomplished by have the instrument view a high reflectance diffuse surface illuminated by the sun. For some of the spectral bands this proves to be much too bright a signal that results in the saturation of detectors designed for measuring low reflectance (ocean) surfaces signals. A mechanical attenuation device in the form of a pin hole screen is used to reduce the signals to calibrate these bands. The sensor response to solar illumination of the SD with and without the attenuation screen in place will be presented. The MODIS detector response to the solar diffuser is smooth when the attenuation screen is absent, but has structures up to a few percent when the attenuation screen is present. This structure corresponds to non-uniform illumination from the solar diffuser. Each pin hole produces a pin-hole image of the sun on the solar diffuser, and there are very many pin hole images of the sun on the solar diffuser for each MODIS detector. Even though there are very many pin-hole images of the sun on the solar diffuser, it is no longer perfectly uniformly illuminated. This non-uniformly illuminated solar diffuser produces intensity variation on the focal planes. The results of a very detailed simulation will be discussed which show how the illumination of the focal plane changes as a result of the attenuation, and the impacts on the calibration will be discussed.

Waluschka, Eugene↗
NASA NTRS2004

The Performance of a Subsonic Diffuser Designed for High Speed Turbojet-Propelled Flight

An initial-phase subsonic diffuser has been designed for the turbojet flowpath of the hypersonic x43B flight demonstrator vehicle. The diffuser fit into a proposed mixed-compression supersonic inlet system and featured a cross-sectional shape transitioning flowpath (high aspect ratio rectangular throat-to-circular engine face) and a centerline offset. This subsonic diffuser has been fabricated and tested at the W1B internal flow facility at NASA Glenn Research Center. At an operating throat Mach number of 0.79, baseline Pitot pressure recovery was found to be just under 0.9, and DH distortion intensity was about 0.4 percent. The diffuser internal flow stagnated, but did not separate on the offset surface of this initial-phase subsonic diffuser. Small improvements in recovery (+0.4 percent) and DH distortion (-32 percent) were obtained from using vane vortex generator flow control applied just downstream of the diffuser throat. The optimum vortex generator array patterns produced inflow boundary layer divergence (local downwash) on the offset surface centerline of the diffuser, and an inflow boundary layer convergence (local upwash) on the centerline of the opposite surface.

Biesiadny, Thomas J.↗
NASA NTRS1950

High-Subsonic Performance Characteristics and Boundary-Layer Investigations of a 12 10-Inch-Inlet-Diameter Conical Diffuser

Performance and boundary-layer data were taken in a 12 degree 10-inch inlet-diameter conical diffuser of 2:1 exit- to inlet-area ratio. These data were taken for two inlet-boundary-layer conditions. The first condition was that of a thinner inlet boundary later (boundary-layer displacement thickness, delta* approximately equal to 0.034) produced by an inlet section approximately 1 inlet diameter in length between the entrance bell and the diffuser. The second condition was a thicker inlet boundary layer (delta* approximately equal to 0.120) produced by an additional inlet section length of approximately 6 diameters. Longitudinal static-pressure distributions were measured fro wall static orifices. Transverse total- and static-pressure surveys were made at the inlet and exit stations. Boundary-layer velocity distributions were measured at seven stations between the inlet and exit. These data were obtained for a Reynolds number (based on inlet diameter) range of 1 x 10(exp 6) to 3.9 x 10(exp 6). The corresponding Mach number range was from M = 0.2 to choking. At the maximum-power-available condition supersonic flow was obtained as far as 4.5 inches downstream from the diffuser inlet with a maximum Mach number of M approximately equal to 1.5. The total-pressure loss through the diffuser in percentage of inlet dynamic pressure was approximately 2.5 percent for the thinner inlet boundary later and 5.5 percent for the thicker inlet boundary later over the lower subsonic range. These valued increased with increasing flow rate- the values for the thicker inlet boundary later more than those for the thinner inlet boundary layer. The diffuser effectiveness, expressed as the ratio of the actual static-pressure rise to the ideal static-pressure rise, was about 85 percent for the thinner inlet boundary layer and about 67 percent for the thicker inlet boundary later in the lower subsonic range. These values decrease with increasing flow rate. Separated flow was observed for both inlet-boundary-layer conditions in the region of adverse pressure gradient just downstream of the transition curvature from inlet section to diffuser. The flow for the thinner-inlet-boundary-layer condition did not fully re-establish itself along the diffuser walls. The thicker inlet-boundary-layer flow, while not completely re-establishing the normal flow pattern downstream of the separated region, did re-establish more successfully than the thinner inlet boundary layer.

Little, B. H., Jr.↗
NASA NTRS1952

Amplitude of Supersonic Diffuser Flow Pulsations

A theoretical method for evaluating the stability characteristics and the amplitude and the frequency of pulsation of ram-jet engines without heat addition is presented herein. Experimental verification of the theoretical results are included where data were available. Theory and experiment show that the pulsation amplitude of a high mass-flow-ratio diffuser having no cone surface flow separation increases with decreasing mass flow. The theoretical trends for changes in amplitude, frequency, and mean-pressure recovery with changes in plenum-chamber volume were experimentally confirmed. For perforated convergent-divergent-type diffusers, a stability hysteresis loop was predicted on the pressure-recovery mass-flow-ratio curve. At a given mean mass-flow ratio, the higher.value of mean pressure recovery corresponded to oscillatory flow in the diffuser while the lower branch was stable. This hysteresis has been observed experimentally. The theory indicates that for a ram-jet engine of given diameter, the amplitude of pulsation of a supersonic diffuser is increased by decreasing the relative size of the plenum chamber with respect to the diffuser volume down to a critical value at which oscillations cease. In the region of these critical values, the stable mass-flow range of the diffuser may be increased either by decreasing the combustion chamber volume or by increasing the length of the diffuser.

Sterbentz, William H.↗
NASA NTRS2010

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗
NASA NTRS2013

Kinetic Monte Carlo Simulation of Cation Diffusion in Low-K Ceramics

Low thermal conductivity (low-K) ceramic materials are of interest to the aerospace community for use as the thermal barrier component of coating systems for turbine engine components. In particular, zirconia-based materials exhibit both low thermal conductivity and structural stability at high temperature, making them suitable for such applications. Because creep is one of the potential failure modes, and because diffusion is a mechanism by which creep takes place, we have performed computer simulations of cation diffusion in a variety of zirconia-based low-K materials. The kinetic Monte Carlo simulation method is an alternative to the more widely known molecular dynamics (MD) method. It is designed to study "infrequent-event" processes, such as diffusion, for which MD simulation can be highly inefficient. We describe the results of kinetic Monte Carlo computer simulations of cation diffusion in several zirconia-based materials, specifically, zirconia doped with Y, Gd, Nb and Yb. Diffusion paths are identified, and migration energy barriers are obtained from density functional calculations and from the literature. We present results on the temperature dependence of the diffusivity, and on the effects of the presence of oxygen vacancies in cation diffusion barrier complexes as well.

Good, Brian↗
NASA NTRS2016

Passive Rocket Diffuser Theory: A Re-Examination of Minimum Second Throat Size

Second-throat diffusers serve to isolate rocket engines from the effects of ambient back pressure during testing without using active control systems. Among the most critical design parameters is the relative area of the diffuser throat to that of the nozzle throat. A smaller second throat is generally desirable because it decreases the stagnation-to-ambient pressure ratio the diffuser requires for nominal operation. There is a limit, however. Below a certain size, the second throat can cause pressure buildup within the diffuser and prevent it from reaching the start condition that protects the nozzle from side-load damage. This paper presents a method for improved estimation of the minimum second throat area which enables diffuser start. The new 3-zone model uses traditional quasi-one-dimensional compressible flow theory to approximate the structure of two distinct diffuser flow fields observed in Computational Fluid Dynamics (CFD) simulations and combines them to provide a less-conservative estimate of the second throat size limit. It is unique among second throat sizing methods in that it accounts for all major conical nozzle and second throat diffuser design parameters within its limits of application. The performance of the 3-zone method is compared to the historical normal shock and force balance methods, and verified against a large number of CFD simulations at specific heat ratios of 1.4 and 1.25. Validation is left as future work, and the model is currently intended to function only as a first-order design tool.

Jones, Daniel R.↗
NASA NTRS2016

Generating a Simulated Fluid Flow over a Surface Using Anisotropic Diffusion

A fluid-flow simulation over a computer-generated surface is generated using a diffusion technique. The surface is comprised of a surface mesh of polygons. A boundary-layer fluid property is obtained for a subset of the polygons of the surface mesh. A gradient vector is determined for a selected polygon, the selected polygon belonging to the surface mesh but not one of the subset of polygons. A maximum and minimum diffusion rate is determined along directions determined using the gradient vector corresponding to the selected polygon. A diffusion-path vector is defined between a point in the selected polygon and a neighboring point in a neighboring polygon. An updated fluid property is determined for the selected polygon using a variable diffusion rate, the variable diffusion rate based on the minimum diffusion rate, maximum diffusion rate, and the gradient vector.

Rodriguez, David L.↗
NASA NTRSDate not supplied

Zircon (U-Th)/He Impact Crater Thermochronometry and the Effects of Shock Microstructures on Helium Diffusion Kinetics

Accurate and precise age determination of impact cratering events remains challenging and often contentious; less than half of all known craters are regarded as accurately and precisely dated. Zircon (U-Th)/He (ZHe) dating of impactites can be employed to date medium to large impact structures as ZHe ages can be fully reset in minutes at T >1000°C, a plausible scenario in the central melt pool. In contrast, complete resetting of ZHe at 200-300°C, encountered near the crater margins or due to post-impact hydrothermal overprinting, may take >103-6 years. To test the reliability of ZHe impact dating, we have quantified the effects of shock-induced microstructures on helium diffusion kinetics in well-characterized variably shocked zircon. We investigated samples from two impact structures, the Chicxulub multi-ring crater and Ries complex crater, to compare diffusion kinetics from structures with different size, age, and hydrothermal system longevity. Shock microstructures were characterized by backscattered-electron imaging prior to determining the He diffusion kinetics by prograde and retrograde fractional-release experiments via light-bulb furnace with incremental step-heating (10°C) from 300°C to 600°C. Next, we examine the internal interconnectivity and sizes of the diffusion domains by EBSD. While we found that zircon with few shock microstructures exhibited no marked deviation from helium diffusion kinetics of undamaged zircon, zircon grains with planar microstructures and granular textures are characterized by a dramatic decrease in helium retentivity due to the reduction in the effective domain size and the introduction of interconnected fast diffusion pathways. A subset of grains were ZHe dated and showed that less deformed grains yielded a weighted mean age within error of the accepted impact ages, while the grains with planar microstructures or granular textures gave systematically younger ages. These new diffusion data and ZHe ages demonstrate that highly shocked grains are unsuitable for ZHe impact crater dating. Therefore, detailed characterization of impact-induced microstructures is critical for determining accurate ZHe impact ages and offers the possibility of investigating post-impact hydrothermal circulation.

Zircon↗
NASA NTRSDate not supplied

Zircon (U-Th)/He Impact Crater Thermochronometry and the Effects of Shock Microstructures on Helium Diffusion Kinetics

Absolute age determination of impact cratering events remains difficult and often controversial; a challenge that has resulted in < 50% of known impact craters regarded as accurately and precisely dated. Besides conventional 40Ar/39Ar and U-Pb methods, zircon (U-Th)/He (ZHe) dating of impactites has been applied to large- to medium-sized impact structures. ZHe dates can be fully reset in minutes at 1000°C, which is commonly reached in central sections of the melt sheet, whereas resetting of ZHe at <300°C, which might be encountered near the crater margins or persist in post-impact hydrothermal systems, may take >103-6 years. There is a critical need to quantify the effects of shock-induced microstructures and impact metamorphism on helium diffusion kinetics in well-characterized, variably shocked zircon to further establish the reliability of (U-Th)/He for dating impacts. For this purpose, we investigated suevite and impact melt samples from two impact structures, the Chicxulub multi-ring basin and the Ries complex crater, which enables us to compare zircon helium diffusion kinetics from impact structures with differing sizes, ages, and hydrothermal system longevities. Shock microstructures were characterized by backscattered-electron (BSE) imaging prior to diffusion step-heating fractional release experiments using light-bulb furnace with prograde and retrograde incremental 10°C steps from 250°C to 600°C. Afterward, we characterize the diffusion domain sizes and their interconnectivity within the shocked zircon grains using electron backscatter diffraction (EBSD). We find that zircon with few shock microstructures exhibit no significant deviation from helium diffusion kinetics of undamaged zircon. In contrast, zircon grains with planar deformation features (PDFs) and granular textures classified by BSE and EBSD are characterized by a dramatic decrease in helium retentivity, similar to radiation damaged grains, due to a reduction in the effective domain size and the introduction of interconnected fast diffusion pathways created by shock microstructures. A subset of grains were dated by ZHe after the external morphology of the grains was determined by BSE imaging. The euhedral grains yielded a weighted mean age within the uncertainty of the accepted impact ages, whereas the grains with PDFs or granular textures exhibited younger ages. Thus, these diffusion experiments and ZHe dates suggest that the dramatic decrease in domain size likely renders shocked grains more susceptible to impact-induced hydrothermal resetting and subsequent overprinting. Hence, characterization of shock microstructures is critical for determining accurate impact ages using ZHe methods especially when applied to previously unconstrained craters. The thermochronometer also offers the opportunity to determine the magnitude and duration of post-impact hydrothermal circulation.

zircon↗
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
DOE OSTI2024

A Finite Difference informed Random Walk solver for simulating radiation defect evolution in polycrystalline structures with strongly inhomogeneous diffusivity

Diffusivity of species and defects on grain boundaries is usually several orders of magnitude larger than that inside grains. Such strongly inhomogeneous diffusivity requires prohibitively high computational demands for modeling microstructural evolution. Here, this paper presents a highly-efficient numerical solver, combining the Finite Difference method and Random Walk model, designed for accurately modeling strongly inhomogeneous diffusion within polycrystalline structures. The proposed solver, termed Finite Difference informed Random Walk (FDiRW), integrates a customized Finite Difference (cFD) scheme tailored for fast diffusion along thin grain boundaries represented by a single-layer of nodes. Numerical experiments demonstrate that the FDiRW solver achieves an impressive efficiency gain of 1560x compared to traditional Finite Difference methods while maintaining accuracy, making it feasible for personal computer machines to handle diffusional systems with strongly inhomogeneous diffusivity across static polycrystalline microstructures. The model has been successfully applied to simulate radiation defect evolution, showcasing its scalability to engineering scales in both length and time dimensions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of lithium diffusion into Ga 2 O 3 thin films

The integration of lithium based compounds (e.g., Li:NiO/Ga 2 O 3 , LiGa 5 O 8 /Ga 2 O 3 ) in Ga 2 O 3 based pn-heterojunctions raises concerns about interface stability, i.e., Li diffusion effects on Ga 2 O 3 properties. In this work the ex-situ diffusion of Li is investigated in three different Ga 2 O 3 epilayers systems [(001) κ-Ga 2 O 3 and (−201) β-Ga2O3 heteroepitaxy on (001) α-Al 2 O 3 , and (010) β-Ga 2 O 3 homoepitaxy] at relevant temperatures for the synthesis / processing of Li-based epilayers. It is here experimentally demonstrated and quantified the Li diffusion in all the investigated Ga 2 O 3 epilayers systems and Li bulk (D Li,bulk ) and 2D defects (D Li,2D ) diffusion coefficients are provided. In the case of the (010) β-Ga 2 O 3 homoepitaxial layer (nominally free of structural defects), hybrid functional theory calculations foresee a diffusion mechanism mediated by Ga vacancies (VGa). Moreover, in the (010) β-Ga 2 O 3 homo-layer a significant effect on its functional properties (e.g., additional Raman vibrational modes, induced conductivity in an otherwise insulating sample) upon the Li-diffusion process is experimentally highlighted and tentatively related to the passivation of acceptor defects (i.e., formation of V Ga -nLi complexes).

Defects↗
DOE OSTI2025

Effect of Mg and Ni impurities on tritium diffusion in lithium ceramics through cluster dynamics simulations

This study investigates the impact of Mg and Ni doping on tritium diffusion in LiAlO 2 and LiAl 5 O 8 ceramics, that are used in tritium-producing burnable absorber rods (TPBARs). Utilizing Centipede simulations across a broad temperature range (500 K to 1250 K), we explore the interplay between defect dynamics, cluster formation, and tritium mobility. In LiAlO 2 , Mg doping significantly enhances tritium diffusivity by increasing tritium interstitial concentrations and diffusion coefficients of key species, thereby doubling the overall tritium diffusivity. Ni doping, while shifting the dominant defect to Li vacancies, maintains high tritium mobility due to the low binding energy of Li vacancy-tritium complexes, which ensures effective tritium migration. In LiAl 5 O 8 , Mg and Ni doping results in a slight reduction in the diffusion coefficients of key species, yet the dramatic increase in tritium interstitial concentrations compensates, leading to a net small increase in tritium diffusivity. In conclusion, the findings highlight the critical role of defects in tritium transport and the effect of Mg and Ni defects on the performance of these ceramics in demanding nuclear environments.

Cluster dynamics↗
DOE OSTI2025

First-principles investigation of cerium and neodymium diffusion in BCC chromium and vanadium via vacancy-mediated transport

Lanthanide transport plays a crucial role in the performance and longevity of metallic nuclear fuels. This study examines the diffusion behavior of Ce and Nd—two major fission products—in body-centered cubic (BCC) Cr and V, which are potential liner or coating materials for mitigating fuel-cladding chemical interactions (FCCI). Using density functional theory (DFT) calculations and self-consistent mean-field (SCMF) analysis, the vacancy-mediated diffusion coefficients are evaluated. Our findings reveal that Ce and Nd act as oversized solutes and are strongly bound to vacancies in BCC Cr and V, with diffusivities in Cr significantly lower than in V and in hexagonal closed-packed (HCP) Zr, as investigated in our previous work. The activation energies for Ce and Nd diffusion are 3.39 and 3.32 eV, respectively, in BCC Cr, and 2.56 and 2.33 eV, respectively, in BCC V. Analysis of vacancy drag and partial diffusion coefficient ratios indicates a strong tendency for lanthanide enrichment at vacancy sinks in BCC Cr, and to a lesser extent in BCC V, with this effect persisting up to the melting point in Cr and remaining substantial for Nd in V at high temperatures. Under irradiation, the increase in vacancy concentration is expected to enhance lanthanide transport, potentially accelerating interactions at liner-cladding interfaces. Although BCC Cr exhibits relatively low lanthanide diffusivities under equilibrium conditions, the expected segregation tendencies under irradiation suggest that Zr liners may be a more favorable option. Further investigations using rate theory, cluster dynamics, and phase-field modeling are required to quantitatively assess the performance of these materials in reactor environments.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

Adsorption and Solvation Modulate Rhodamine B Diffusion in Ethanol/Water-Filled Anodic Alumina Nanopores

Confinement of solvents and solutes within nanoporous materials frequently leads to the emergence of unique mass transport behaviors that, once fully understood, may lead to improved chemical separations. Here, the diffusion of Rhodamine B (RhB) dye within 10 and 20 nm diameter anodic aluminum oxide (AAO) nanopores filled with binary ethanol/water mixtures is investigated. Mixture compositions spanning from pure ethanol to pure water are employed. The results of confocal fluorescence correlation spectroscopy studies reveal that RhB diffusion occurs by a two-component mechanism comprising composition-dependent fast and slow motions, characterized by diffusion coefficients D f and D s . The results are consistent with those of previous studies performed under more limited conditions [J. Phys. Chem. C, 2023, 127, 411-420]. The fast component scales with mixture viscosity and is assigned to hindered bulk-like diffusion in central pore regions. Slow diffusion likely involves adsorption of RhB to the pore surface and may be described by a desorption mediated mechanism. The occurrence of RhB adsorption to the AAO surface is verified at the single-molecule level by wide-field fluorescence imaging of membrane cross-sectional surfaces. Unique composition dependent trends in the autocorrelation amplitude and in D s that mimic bulk RhB solubility are revealed. D s is found to be smallest in pure ethanol and pure water and largest in intermediate mixtures. These results suggest that RhB surface adsorption is strongest in the pure liquids and weakest in mixtures of intermediate composition, where the dye is least soluble, and most soluble, respectively. As a result, molecular dynamics simulations reveal that a water layer appears on the pore surface under most conditions, while RhB is solvated primarily by ethanol. The composition dependence of RhB diffusion is concluded to reflect its solvation dependent interactions with the pore walls.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Experimental Characterization of Hydrogen Diffusion in Shale Rocks for Geologic Storage Applications

As global energy systems undergo a transition to cleaner alternatives, geologic hydrogen storage has emerged as a promising solution for large-scale energy storage. A critical factor in determining the feasibility of this approach is the effectiveness of caprock formations, such as shale, in preventing hydrogen migration. This study investigates the diffusion behavior of hydrogen through shale to assess its suitability as a caprock for geologic hydrogen storage. Using a novel double-seal core holder design and a through-diffusion apparatus, hydrogen diffusion was measured through shale rock from the Eagle Ford and Wolfcamp Formations under dry conditions. These measurements were complemented by microstructural and mineralogical analyses using low-pressure nitrogen adsorption and X-ray diffraction. The effective diffusion coefficient of hydrogen in these shale caprocks ranged from 2.51 × 10 –8 to 9.85 × 10 –8 m 2 /s. Notably, we observed that the diffusion behavior was more related to the pore network structure and could not be attributed to differences in the total pore volume between shale types alone. Here, to further understand the role of pore network complexity, a fractal pore model was developed to correlate tortuosity with the fractal dimension of the pore structure (a measure of pore network complexity). The proposed model closely matched tortuosity values obtained from diffusion experiments, outperforming existing theoretical tortuosity–porosity correlations. These findings provide key quantitative parameters needed to assess the feasibility of geologic hydrogen storage as well as insights that can be applied to hydrogen storage in a range of geologic formations.

08 HYDROGEN↗
DOE OSTI2025

Defect-Mediated Diffusion Pathways in Spodumene Accelerate Lithium Transport

Lithium extraction from naturally occurring α-spodumene is hindered by poor lithium diffusivity, necessitating high-temperature phase transformation to a low-density β polymorph. Although β spodumene exhibits up to 5 orders of magnitude higher lithium-ion diffusivity, both phases have diffusion activation energies between 0.8 and 1 eV, indicating that polymorph density is not the controlling factor over diffusivity. We show that aluminum vacancies facilitate lithium-ion diffusion in α-spodumene by reducing the migration barrier from 2.4 to 0.9 eV. Bond valence site energy and nudged elastic band calculations show a new lithium local minimum site which promotes a one-dimensional percolation network by reducing the lithium intersite distance from 4.5 Å to 2.9 Å. However, aluminum vacancies are energetically unfavorable to percolate through the whole structure, resulting in very low net lithium diffusivity and highlighting the critical role of nonstoichiometric defects in facilitating lithium transport in rigid aluminosilicate structures.

Chemical structure↗