Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Temperature Dependence of Diffusivities in Liquid Elements (LMD)

This research was to advance the understanding of diffusion mechanisms in liquid metals and alloys through accurate diffusivity measurements over a wide range of temperatures, including the proximity of the materials melting points. Specifically, it was driven towards developing a methodology (and subsequent flight hardware) to enable several diffusion coefficient measurements (i.e., at several different temperatures) to be performed using a single sample. The Liquid Metal Diffusion (LMD) was funded as a Flight Definition Project in February 1993 in response to NRA 91-OSSA-20 (Microgravity Science and Applications Division). The Science Concept Review for LAID was held during April 1994. In January 1995 we were informed that we had failed this review and the project was change to ground-based activities only. A new proposal was submitted for the next NRA addressing the panels concerns. As part of NASA's Risk Mitigation program, a scaled-down version of the hardware was funded in July of 1995 for a flight opportunity utilizing experiment on the Microgravity Isolation Mount. This experiment was to determine the self-diffusivity of indium at 185 C. The LMD was transferred to the Mir Space Station in STS-81 and returned on STS-84 (January - May 1997). Three, out of five, self-diffusion data sets were returned. A description of this experiment/hardware is included below. This summary is only intended to give the reader an overview of the results obtained for the tasks outlined in the original proposal. Research that was not published is explained in more detail. At the end of this report is a list of refereed publications and invited talks that were given as a result of this work. The reader is directed to these for further details. Attachment: Real-time diffusivity measurements in liquids at several temperatures with one sample, On the insensitivity of liquid diffusivity measurements to deviations from 1D transport, and Numerical simulations of the convective contamination of diffusivity measurements in liquids.

Banish, R. Michael↗
NASA NTRS2002

Flow/Soot-Formation Interactions in Nonbuoyant Laminar Diffusion Flames

This is the final report of a research program considering interactions between flow and soot properties within laminar diffusion flames. Laminar diffusion flames were considered because they provide model flame systems that are far more tractable for theoretical and experimental studies than more practical turbulent diffusion flames. In particular, understanding the transport and chemical reaction processes of laminar flames is a necessary precursor to understanding these processes in practical turbulent flames and many aspects of laminar diffusion flames have direct relevance to turbulent diffusion flames through application of the widely recognized laminar flamelet concept of turbulent diffusion flames. The investigation was divided into three phases, considering the shapes of nonbuoyant round laminar jet diffusion flames in still air, the shapes of nonbuoyant round laminar jet diffusion flames in coflowing air, and the hydrodynamic suppression of soot formation in laminar diffusion flames.

Dai, Z.↗
NASA NTRS2007

Diffusion Flame Stabilization

Diffusion flames are commonly used for industrial burners in furnaces and flares. Oxygen/fuel burners are usually diffusion burners, primarily for safety reasons, to prevent flashback and explosion in a potentially dangerous system. Furthermore, in most fires, condensed materials pyrolyze, vaporize, and burn in air as diffusion flames. As a result of the interaction of a diffusion flame with burner or condensed-fuel surfaces, a quenched space is formed, thus leaving a diffusion flame edge, which plays an important role in flame holding in combustion systems and fire spread through condensed fuels. Despite a long history of jet diffusion flame studies, lifting/blowoff mechanisms have not yet been fully understood, compared to those of premixed flames. In this study, the structure and stability of diffusion flames of gaseous hydrocarbon fuels in coflowing air at normal earth gravity have been investigated experimentally and computationally. Measurements of the critical mean jet velocity (U(sub jc)) of methane, ethane, or propane at lifting or blowoff were made as a function of the coflowing air velocity (U(sub a)) using a tube burner (i.d.: 2.87 mm) (Fig. 1, left). By using a computational fluid dynamics code with 33 species and 112 elementary reaction steps, the internal chemical-kinetic structures of the stabilizing region of methane and propane flames were investigated (Fig. 1, right). A peak reactivity spot, i.e., reaction kernel, is formed in the flame stabilizing region due to back-diffusion of heat and radical species against an oxygen-rich incoming flow, thus holding the trailing diffusion flame. The simulated flame base moved downstream under flow conditions close to the measured stability limit.

Takahashi, Fumiaki↗
NASA NTRS2008

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. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗
NASA NTRS2006

Diffusion Flame Stabilization

Diffusion flames are commonly used for industrial burners in furnaces and flares. Oxygen/fuel burners are usually diffusion burners, primarily for safety reasons, to prevent flashback and explosion in a potentially dangerous system. Furthermore, in most fires, condensed materials pyrolyze, vaporize, and burn in air as diffusion flames. As a result of the interaction of a diffusion flame with burner or condensed-fuel surfaces, a quenched space is formed, thus leaving a diffusion flame edge, which plays an important role in flame holding in combustion systems and fire spread through condensed fuels. Despite a long history of jet diffusion flame studies, lifting/blowoff mechanisms have not yet been fully understood, compared to those of premixed flames. In this study, the structure and stability of diffusion flames of gaseous hydrocarbon fuels in coflowing air at normal earth gravity have been investigated experimentally and computationally. Measurements of the critical mean jet velocity (U(sub jc)) of methane, ethane, or propane at lifting or blowoff were made as a function of the coflowing air velocity (U(sub a)) using a tube burner (i.d.: 2.87 mm). By using a computational fluid dynamics code with 33 species and 112 elementary reaction steps, the internal chemical-kinetic structures of the stabilizing region of methane and propane flames were investigated. A peak reactivity spot, i.e., reaction kernel, is formed in the flame stabilizing region due to back-diffusion of heat and radical species against an oxygen-rich incoming flow, thus holding the trailing diffusion flame. The simulated flame base moved downstream under flow conditions close to the measured stability limit.

Takahashi, Fumiaki↗
DOE OSTI2019

A Prediction of Uranium(VI) Diffusion in Montmorillonite at Various Chemical Solution Conditions

Most nuclear waste disposal options include compacted bentonite, predominately consisting of montmorillonite clay, as part of a barrier system to minimize contaminant mobility. Uranium (U) is the primary element in spent nuclear fuel, and a potential contaminant of water resources. Uranium sorption onto clay and its slow diffusive transport away from waste canisters are expected to limit U(VI) mobility, and possibly also control nuclear fuel degradation rates based on mass action considerations. A prediction of U(VI) diffusion in these systems is complicated by two main factors: (1) the complex mineralogical structure of montmorillonite clay, leading to two types of clay porosities and surface sites, and (2) the dependence of U(VI) sorption on U(VI) solution speciation and the specific chemical system conditions. For instance, depending on pH and ionic strength, U(VI) may sorb onto montmorillonite due to weak ion exchange reactions or the formation of stable surface complexes. Furthermore, in contrast to conventional porous media, U(VI) diffusion in montmorillonite is not only limited by contaminant sorption, but also by a (potential) lack of access to the full clay porosity for some U(VI) solution species. For instance, at the high degrees of clay compaction expected in future barrier systems, the partial or full exclusion of U(VI) from negatively-charged clay interlayer spaces may substantially decrease U(VI) fluxes. In this modeling study, we specifically investigated how the charge of U(VI) solution species affects U(VI) sorption and diffusion behavior. For this purpose, we first simplified complex chemical solution speciation diagrams to charge speciation diagrams. Then, we coupled an existing U(VI) surface complexation model (Tournassat et al., Geochim. Cosmochim. Acta, 2018) with an analytical solution of the diffusion equation to simulate U(VI) mobility in lab-scale through-diffusion experiments for various chemical conditions. Our results suggest that there are two important parameters that govern U(VI) diffusion in these systems. First, U(VI) surface complexation on clay edge surfaces leads to contaminant retardation. On the field scale, this will determine how early a plume may cross site boundaries and pose risks to drinking water resources. Second, the degree of access of U(VI) solution species to diffuse layers in clay interlayer spaces largely controls U(VI) diffusive fluxes. In the field, this parameter will govern the level of U(VI) concentrations that can be expected in contaminated waters.

Pistorino, Jonathan C.↗
DOE OSTI2018

Neoclassical Diffusion of Radiation-Belt Electrons Across Very Low L -shells

In the presence of drift-shell splitting intrinsic to the IGRF magnetic field model, pitch-angle scattering from Coulomb collisions experienced by radiation-belt electrons in the upper atmosphere and ionosphere produces extra radial diffusion, a form of neoclassical diffusion. The strength of the neoclassical radial diffusion at L < 1.2 exceeds that expected there from radial-diffusion mechanisms traditionally considered, and decreases with increasing L-shell. In this study we construct a numerical model for this coupled (radial and pitch-angle) collisional diffusion process and apply it to simulate raw count-rate data observed aboard the Gemini spacecraft for several years after the 1962 Starfish nuclear detonation. The data show apparent lifetimes 10-100 times as long as would have been expected from collisional pitch-angle diffusion and Coulomb drag alone. Our model reproduces apparent lifetimes for >0.5-MeV electrons in the region 1.14 < L < 1.26 to within a factor of two (comparable to the uncertainty quoted for the observations). We conclude that neoclassical radial diffusion (resulting from drift-shell splitting intrinsic to IGRF's azimuthal asymmetries) mitigates the decay expected from collisional pitch-angle diffusion and inelastic energy loss alone and thus contributes importantly to the long apparent lifetimes observed at these low L-shells.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Probing the Impact of Vacancy Diffusion on Void Dynamics at the Lithium Metal–Solid Electrolyte Interface

Lithium (Li) metal-based solid-state batteries (SSBs) are considered promising candidates for next-generation energy storage due to their superior energy density and enhanced safety compared to conventional Li-ion systems. However, their practical application is limited by challenges such as void formation at the Li-solid electrolyte (SE) interface, which disrupts ion transport and accelerates interfacial degradation. This work investigates how the coupled effects of electro-dissolution kinetics and surface diffusion at the Li metal surface govern the evolution of interfacial morphology during stripping. This work examines the influence of three distinct surface diffusion modes, which are terrace diffusion, step diffusion, and interlayer diffusion, on maintaining interfacial stability. In addition, how the dominant surface diffusion mechanism can overcome the contact loss due to high reaction kinetics is explored. Furthermore, the roughness of the Li metal anode surface is quantified, and the influence of different diffusion mechanisms on the evolution of the dynamic solid–solid interface is examined. The critical role of temperature in enhancing Li surface diffusivity and expanding the regime of stable contact is highlighted. By identifying distinct regimes of interface stability, this study analyzes how non-uniform electrochemical dynamics dictate void morphology evolution and interfacial contact. These insights offer guiding principles for engineering robust Li–SE interfaces in SSBs.

lithium metal anode↗
DOE OSTI2024

Cluster dynamics simulations of tritium and helium diffusion in lithium ceramics

Tritium (T) and He diffusion in LiAlO 2 and LiAl 5 O 8 phases influences the performance of tritium producing burnable absorber rods (TPBARs) by affecting the gas release, swelling and thermal conductivity of Li-bearing ceramic pellets. Frenkel pair defects and clusters created by irradiation can attract T and He interstitials and form clusters of the type He i x Li, He i x Al, He i x O, T i x Li, T i x Al, and T i x O, 1 ≤ x ≤ 4 in a Li, Al or O vacancy site (notation denotes x He or T atoms in a 1 Li, 1 Al or 1 O vacant site). The concentration and mobility of each of these clusters collectively contribute to the diffusion of the He and T gases in LiAlO 2 and LiAl 5 O 8 . Here, in this work, free energy cluster dynamics simulations implemented in the Centipede code, are used to obtain the concentration and diffusivities of these clusters which are then used to calculate the total diffusivity of T and He gases in LiAlO 2 and LiAl 5 O 8 . The results show that diffusivity of T is at least one order of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 whereas He diffusion is 2–13 orders of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 . There is a higher concentration of highly diffusive species (T interstitials and T i 03 Li for the case of tritium and He i 01 Li, He i 02 Li, and He i 03 Li for the case of He) in LiAlO 2 than in LiAl 5 O 8 which increase the total diffusion of T and He in LiAlO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

Diffusion-limited hydrolysis in polymeric materials

Polymer degradation under aggressive environmental stressors often develops heterogeneities due to diffusion-limited reaction phenomena. This is well established for diffusion-limited oxidation (DLO), which is known to occur for most polymeric materials at elevated temperatures but has been less summarized for the conditions of diffusion-limited hydrolysis (DLH). Here, an overview of hydrolysis for several materials and a computational model, analogous to the underlying equations for DLO, able to define this diffusion-reaction system is presented. A systematic study of the influence of various parameters, such as water diffusivity, reaction rate and order, and a more in-depth focus on residual isocyanate hydrolysis in a polymeric methylene diphenyl diisocyanate (pMDI) based polyurethane (PU) foam is given. For this system, we present experimental data for model ‘input’ parameters and discuss predictions for different conditions. We conceptually compare the behavior of diffusion-limited oxidation to that of diffusion-limited hydrolysis (DLH). With the mathematical framework and key material properties presented herein, any DLH phenomena following Fickian diffusion behavior can be understood, modeled, and predicted.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Study of Tritium Diffusivity in Pure and Sn-Defective Zr: A First-Principles Density Functional Theory Approach

Zirconium alloys (e.g., zircaloy-4) are used as tritium ( 3 H) getter materials in tritium-producing burnable absorber rods (TPBARs) owing to their ability to capture 3 H and chemically convert them into metal hydrides. Understanding of 3 H diffusion mechanisms in zircaloy is crucial for the optimal design of material performance in nuclear technology. Here, for this work, we perform first-principles density functional theory calculations to study the 3 H diffusion mechanism in pure and impure Zr with a low concentration of tin (Sn) atoms to determine the impact of the presence of Sn on the movement of 3 H atoms through the material. First, we calculated the diffusion barriers for 3 H in pure Zr by taking different migration pathways. We then introduced a low concentration of Sn impurity and systematically explored the impurity effect on the diffusion barriers for 3 H. Using our calculated diffusion energy barriers, we further obtained the diffusion coefficients and analyzed the results by comparing them with the experimental and previously calculated values. A diffusion coefficient of the order of 10 –8 m 2 /s is predicted. We also found that the presence of a Sn impurity could reduce the diffusivity up to 4 orders of magnitude. In conclusion, our results could serve as guidelines for further experimental investigations.

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

Electrode and Microstructure Dependence of Oxygen Diffusion in Ferroelectric Hafnium Zirconium Oxide Thin Films

Hafnia-based ferroelectrics hold promise to reduce energy demand for computing by enabling compute-in-memory and as non-volatile memories. The ferroelectric phase in this material system is, in part, stabilized by oxygen vacancies. While oxygen vacancies may be a necessity for phase stability, they limit device endurance through diffusion and accumulation into conducting channels. Herein, it is shown that oxygen diffusion is spatially variable within individual grains of ferroelectric hafnium zirconium oxide (HZO). Using 18 O tracers and finite difference modeling, it is shown that grain boundaries and regions near electrode interfaces allow for relatively rapid oxygen diffusion, with values as much as 10 4 larger than the grain cores. Further, the selection of electrode material affects the diffusion coefficients across all microstructural regions. HZO films in contact with TiN electrodes result in more oxygen-deficient HZO films and higher oxygen diffusion coefficients. Tungsten electrodes result in fewer vacancies and lower diffusion coefficients. Diffusion activation energy differences between the HZO with the two electrodes is reconciled by differing populations of charged and uncharged oxygen vacancies. This insight into the local vacancy populations and diffusion pathways provides a platform for designing hafnia-based films, deposition processes, and integration strategies to reduce vacancy gradients and improve performance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quantifying Anionic Diffusion in 2D Halide Perovskite Lateral Heterostructures

Anionic diffusion strongly impacts the stability of halide perovskite materials, but it is still not well understood. Here, a quantitative investigation of in-plane thermally driven anionic inter-diffusion in a series of novel 2D and quasi-2D halide perovskites lateral heterostructures is reported. The calculated diffusion coefficients (D) reveal the inhibition of Br–I inter-diffusion with bulky π-conjugated organic cations compared with short-chain aliphatic organic cations. Furthermore, halide diffusion is found to be faster in quasi-2D (n > 1) than 2D perovskites (n = 1). The increment becomes less apparent as the “n” number increases, akin to the quantum confinement effect observed for band gaps. These trends are rationalized by molecular dynamics simulations of free energy barriers for halide diffusion that reveal mechanisms for suppressing diffusion. We report this work provides important fundamental insights on the anionic migration and diffusion process in halide perovskite materials.

2D materials↗
DOE OSTI2017

Hydrophilic domain structure in polymer exchange membranes: Simulations of NMR spin diffusion experiments to address ability for model discrimination

Here, we detail the development of a flexible simulation program (NMR_DIFFSIM) that solves the nuclear magnetic resonance (NMR) spin diffusion equation for arbitrary polymer architectures. The program was used to explore the proton ( 1 H) NMR spin diffusion behavior predicted for a range of geometrical models describing polymer exchange membranes. These results were also directly compared with the NMR spin diffusion behavior predicted for more complex domain structures obtained from molecular dynamics (MD) simulations. The numerical implementation and capabilities of NMR_DIFFSIM were demonstrated by evaluating the experimental NMR spin diffusion behavior for the hydrophilic domain structure in sulfonated Diels-Alder Poly(Phenylene) (SDAPP) polymer membranes. The impact of morphology variations as a function of sulfonation and hydration level on the resulting NMR spin diffusion behavior were determined. These simulations allowed us to critically address the ability of NMR spin diffusion to discriminate between different structural models, and to highlight the extremely high fidelity experimental data required to accomplish this. A direct comparison of experimental double-quantum-filtered 1 H NMR spin diffusion in SDAPP membranes to the spin diffusion behavior predicted for MD-proposed morphologies revealed excellent agreement, providing experimental support for the MD structures at low to moderate hydration levels.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Time-Dependent Diffusion Coefficient of Fe in MgB{sub 2} Superconductors

The iron (Fe) diffusion in superconducting MgB{sub 2} bulk samples has been studied for sintering time durations of 15 min, 30 min, 1 h, 2 h, and 4 h at 900{sup ∘}C. Fe coating bulk polycrstalline superconducting MgB{sub 2} samples for Fe coating were prepared by pelletizing and used in the diffusion experiments with initial sintering at 800{sup ∘}C for 1 h. A thin layer of Fe was coated on MgB{sub 2} pellets by evaporation in vacuum. Effects of Fe diffusion on the structural, electrical, and superconducting properties of MgB{sub 2} have been investigated by X-ray diffraction (XRD), scanning electron microscopy (SEM), infrared spectroscopy (IR), energy-dispersive X-ray spectroscopy (EDS), and resistivity measurements. Fe diffused samples have slightly increased critical transition temperatures and have larger lattice parameter c values, in comparison with bare samples. Fe diffusion coefficients were calculated from depth profiles of c parameter and room temperature resistivity values. Depth profiles were obtained by successive removal of thin layers from Fe diffused surfaces of the samples. Our results have shown that the Fe diffusion coefficient decreases with increasing sintering time and resistivity measurements can be utilized for determination of diffusion coefficient.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Mutual Diffusion Coefficients and Refractive Index Increments of K 2 SO 4 (aq) at 298.15 K from Rayleigh Interferometry

Here, isothermal mutual diffusion coefficients (interdiffusion coefficients) were measured for K 2 SO 4 (aq) at 298.15 ± 0.005 K, at numerous concentrations ranging from dilute solutions to near saturation (0.59648 mol∙dm -3 ; 0.61349 mol∙kg -1 ) under free diffusion boundary conditions, using high precision Rayleigh interferometry. Under the experimental conditions these diffusion coefficients are on the volume-fixed reference frame D v . Two series of experiments were performed, the first using the traditional experimental approach with a mercury lamp source and with the interference patterns being recorded on glass photographic plates, and the second with a He-Ne laser light source and computerized data acquisition using a photodiode array. The results from both series of experiments are in excellent agreement, and generally yield diffusion coefficients precise to at least 0.002 x 10 -9 m 2 ∙s 1 (0.15% to 0.19%) and in most cases to 0.001 x 10 -9 m 2 ∙s -1 (0.07% to 0.1%). These experiments also yield accurate values of the refractive index differences for the solution pairs used in the diffusion experiments. The new diffusion coefficients are compared to two sets of published values of diffusion coefficients for K 2 SO 4 (aq) which are somewhat discrepant from each other. This study extends and complements our earlier work on the diffusion coefficients of the most common brine salts: NaCl(aq), KCl(aq), MgCl 2 (aq), CaCl 2 (aq), Na 2 SO 4 (aq), ] MgSO 4 (aq), NaHCO 3 (aq), and KHCO 3 (aq) at 298.15 K.

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

Leveraging inter-firm influence in the diffusion of energy efficiency technologies: An agent-based model

Energy efficiency technologies (EETs) are crucial for saving energy and reducing carbon dioxide emissions. However, the diffusion of EETs in small and medium-sized enterprises is rather slow. Literature shows the interactions between innovation adopters and potential adopters have significant impacts on innovation diffusion. Enterprises lack the motivation to share information, and EETs usually lack observability, which suppress the inter-firm influence. Thus, an information platform, together with proper policies encouraging or forcing enterprises to disclose EET-related information, should help harness inter-firm influence to accelerate EETs' diffusion. To explore whether and how such an information platform affects EETs' diffusion in small and medium-sized enterprises, this report builds an agent-based model to mimic EET diffusion processes. Based on a series of controlled numerical experiments, some counter-intuitive phenomena are discovered and explained. The results show that the information platform is a double-edged sword that notably accelerates EETs' diffusion by approximately 47% but may also boost negative information to diffuse even faster and delay massive adoption of EETs. Increasing network density and the intensity of inter-firm influence are effective to speed EET diffusion, but their impacts diminish drastically after reaching some critical values (0.05 and 0.15 respectively) and eventually harm the stability of the system. Ultimately, the findings implicate that EET suppliers should carefully launch their promising but immature products; policies that can reduce the perceived risk by enterprises and the effort to maintain an informative rather than judgmental information platform can prominently mitigate the negative side effects brought by high fluidity of information.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗