Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusivity”

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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel

Chemical Interaction of Palladium in Uranium Oxycarbide Nuclear Fuel Kernel Jana Howard1,3, Guang Yang3, Haiyan Zhao1*, Patrick Warren4, Tiankai Yao3, Steven Cavazos4 Elizabeth Sooby Wood4*, Ching-heng Shiau2 1 University of Idaho, Environmental Science, Idaho Falls, Idaho, USA 2 Boise State University, Microscopy and Characterization Suite, Idaho Falls, Idaho, USA 3 Idaho National Laboratory, Idaho Falls, Idaho, USA 4 University of Texas San Antonio, Department of Physics and Astronomy, San Antonio, Texas, USA *Corresponding author: haiyanz@uidaho.edu; elizabeth.soobywood@utsa.edu Tri-structural Isotropic (TRISO) particle fuel is the preferred choice for the newest and next-generation high-temperature nuclear reactors due to its robust construction [1]. The fuel design effectively contains fission products inside the particle at extremely high temperatures [2]. However, the effect of transition metal fission products on fuel performance and structural integrity remains largely unknown, creating uncertainties in predicting fuel performance and potential failure mechanisms [2]. For example, palladium (Pd) has a strong tendency to form intermetallic phases with uranium (U). These new phases tend to be hard and brittle which could lead to fracture of fuel kernel during irradiation [3,4]. Pd is not normally produced in high concentrations in normal fission reactions however, driving the production of Pd into the locality of a Uranium Oxycarbide (UCO) nuclear fuel provides an opportunity to closely observe the diffusion and chemical interactions. This study focuses on detailed transmission electron microscopy characterization of the intermetallic phases formed by the interaction between a UCO kernel and a Pd bar after annealing at 1100 °C for 100 hours. Figure 1 provides an overview of the UCO kernel in contact with the Pd bar. The UCO-Pd sample surface was examined using a Focused Ion Beam (FIB) Quanta 3D FEG FEI in VCD detector mode at 10kV, 83pA to reveal surface features. Several key surface features in the interaction region were observed including: (1) a small area with dendritic microstructure, (2) color contrast near the UCO kernel and Pd boundary, and (3) darker and lighter regions throughout the sample surface. To analyze diffusion behavior, FIB was used to prepare lamellae from five different areas of the UCO-Pd sample, as well as one for the as-received sample. These lamellae were examined using a Scanning Transmission Electron Microscope ThermoFisher Spectra 300 (STEM-Spectra). Energy Dispersive Spectroscopy results confirmed the diffusion between the UCO fuel and Pd bar. Figure 2 shows the identified phases of UC, UO2, and Pd dispersed throughout the interaction region. It was also found that the Pd concentration decreases as the radial distance from the UCO-Pd interface increases. The Pd concentration was 36.93% atomic fraction 11µm from the interaction region then decreased to 7.51% atomic fraction at 257 µm away from the interaction region. Figure 3 shows how concentrations of U and Pd vary across the interaction zone, highlighting the extent of diffusion. This study confirms that Pd interacts with surrounding material to form U-Pd phase and diffusion zones. These diffusion zones vary in composition depending on its radial distance from the point of contact with the Pd bar. These findings contribute to a deeper understanding of the fission product behavior in high temperature nuclear fuels, aiding in the prediction and mitigation of potential fuel degradation mechanisms. A B C Fig. 1. SEM micrographs of UCO fuel kernel in contact with solid Pd bar and the lift out locations. Figure 1A shows UCO fuel kernel in contact with solid Pd bar. Figure 1B shows a higher magnification image of the UCO fuel-Pd interaction zone. Figure 1C shows lift out sites for the lamellae. A B C Fig 2. EDS maps reveal the distribution of U, O, and Pd in location 3. Figure 2A shows the distribution of uranium. Figure 2B shows the distribution of oxygen. Figure 2C shows the distribution of palladium. A B C Fig 3. Palladium and uranium concentration across interaction zone in location 2. Figure 3A shows EDS map of lift out number 2. Figure 3B shows a zoomed in EDS map taken from the interaction zone. Figure 3C shows a line graph of uranium and palladium concentrations across the interaction zone. References: 1. B.E. Wells, N.R. Phillips, K.J. Geelhood. Pacific West Laboratory. (2021). TRISO Fuel: Properties and Failure Modes. https://www.nrc.gov/docs/ML2117/ML21175A152.pdf (Accessed January 16, 2025). 2. American Nuclear Society. TRISO Fuel Development Progresses in INL, ORNL. https://www.gen-4.org/gif/upload/docs/application/pdf/2014-03/nov13nn_fuel_reprint.pdf (Accessed January 16, 2025) 3. Clark R.A., M.A. Con

36 - MATERIALS SCIENCE↗
DOE OSTI2021

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nd, SbNd and Sb 3 Nd 4 and their interactions with the cladding alloy HT9

Lanthanide fission products, such as neodymium, formed during the irradiation of metallic fuels are known to cause deleterious effects from chemical interactions occurring at the fuel-cladding interface; a phenomenon known as fuel-cladding chemical interaction (FCCI). The use of fuel-based additives that bind with the lanthanide elements within the fuel meat, alleviating their interactions at the fuel-cladding interface, is one potential method proposed to mitigate the FCCI phenomenon and extend the burnup potential of such metallic fuel systems. In this study, antimony (Sb) is evaluated as one such additive, and neodymium (Nd) is used to represent the lanthanides. A Sb-Nd alloy is fabricated, which consists of two intermetallic phases, SbNd and Sb 3 Nd 4 . Isothermal diffusion couple experiments are carried out at 675 °C for 24 h between Nd and the Fe-12Cr based HT9 cladding alloy. The results are compared against similar diffusion couple experiments carried out between the Sb-Nd alloy and also HT9. Inter-diffusion between Nd/HT9 diffusion couple is characterized and the phases Fe 17 Nd 5 and Fe 17 Nd 2 were found to form whereas no interactions were observed in the Sb-Nd/HT9 diffusion couple. The lack of compound forming tendencies between Fe, the primary alloying constituent of HT9, and SbNd is elucidated through density functional theory (DFT) calculations on enthalpy of mixing, aligning well with the experimental observations. Lastly, the strong binding of Nd with Sb appears to be favorable, alleviating Nd interactions with HT9 constituent elements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Mitigation and Diagnosis of Pin-Hole Formation in Polymer Electrolyte Membrane Fuel Cells

In this work, we implement a calendering technique to flatten stray fibers within the gas diffusion media, thereby mitigating pin-hole formation in the hot-pressed MEAs. We have investigated the influence of calendering on the long-term durability for several types of gas diffusion electrodes (GDEs) using a combined chemical and mechanical accelerated stress test (AST). The calendered MEAs demonstrate an average AST lifetime improvement of 77% relative to the as-fabricated MEAs.

AMR↗
DOE OSTI2020

Evaluation of Amorphous Oxide Coatings for High-Voltage Li-Ion Battery Applications Using a First-Principles Framework

Cathode surface coatings are widely used industrially as a means to suppress degradation and improve electrochemical performance of lithium-ion batteries. However, developing an optimal coating is challenging, as different coating materials may enhance one aspect of performance while hindering another. Furthermore, to elucidate the fundamental thermodynamic and transport properties of amorphous cathode coating materials, here, we present a framework for calculating and analyzing the Li + and O 2- transport and the stability against delithiation in such materials. Our framework includes systematic workflows of ab-initio molecular dynamics calculations to obtain amorphous structures and diffusion trajectories coupled with an analysis of critical changes of the active-ion local environment during diffusion. Based on these data, we provide an estimate of room-temperature diffusivities, including statistical error bars, and the evaluation of the coating suitability in terms of its ability to facilitate Li + transport while blocking O 2- transport. Finally, we add the thermodynamic stability analysis of the coating chemistry within the operating voltage of common Li-ion cathodes. We apply this framework to two commonly used amorphous coating materials, Al 2 O 3 and ZnO. We find that (1) in general, a higher Li + content increases both Li + and O 2- diffusivities in both Al 2 O 3 and ZnO. Also, Li + and O 2- diffuse much faster in ZnO than in Al 2 O 3 . (2) However, neither Al 2 O 3 nor ZnO is expected to retain a significant concentration of Li + at high charge. (3) ZnO performs much more poorly in terms of O 2- blocking, and hence, Al 2 O 3 is preferred for high-voltage cathode applications. These results will help to quantitatively evaluate amorphous materials, such as metal oxides and fluorides, for different performance metrics and facilitate the development of optimal cathode coatings.

25 ENERGY STORAGE↗
DOE OSTI2020

Polysulfide entrapment and retardation in gel electrolyte Li–S batteries: experiments and modeling

Gelling the electrolyte in lithium–sulfur (Li–S) batteries provides advantages in terms of safety and mitigating polysulfide shuttling leading to reduced side reactions and dendrite growth on the anode side. In the current study, organic and ceramic crosslinkers have been used to improve the entrapment of polysulfides and prevent them from shuttling to the anode side. However, the performance of gel electrolytes is still far from being on par with liquid electrolytes in terms of raw capacity. The discharge curves exhibit an extended first plateau for gel electrolyte batteries compared to liquid electrolytes, showing that polysulfides are better trapped in the gel system. Yet, in the conversion of soluble polysulfides into insoluble polysulfides, the second plateau was delayed. Qualitative and quantitative analysis of the polysulfide diffusion in the liquid electrolyte and gel electrolyte systems were also conducted to predict the relative diffusivity ratio between the liquid electrolyte and gel electrolyte for simulation. To probe what happens at the end of the discharge, where the mass transport in the gel electrolyte appears to be slow, a modified discharge procedure was used where the cell is first paused after discharge and then discharged a second time before recharging normally. Coin cell tests for two gel electrolytes in addition to liquid electrolyte show up to 7% capacity recovery after the pause, which indicates that the diffusion of species during the pause is responsible for the failure recovery. Additionally, the recovery is higher for gel electrolytes compared to liquid, for higher C-rates, and for longer pause times. To better understand this behavior, a Li–S numerical model with mass transport limited reactions was used to examine the polysulfide diffusion coefficients expected in different electrolyte systems. The model is able to reproduce the trends seen in experiments and yields a higher recovery for smaller Li 2 S 4 diffusion coefficients, suggesting that insufficient Li 2 S 4 mass transport is responsible for failure at the end of discharge. By identifying the effect that gel electrolytes have on the state of the battery during discharge using both experimental and simulation methods, this work is a step towards leveraging the advantages of gel electrolytes and mitigating their weaknesses.

25 ENERGY STORAGE↗
DOE OSTI2022

Simultaneously Enhancing Exciton/Charge Transport in Organic Solar Cells by an Organoboron Additive

Efficient exciton diffusion and charge transport play a vital role in advancing the power conversion efficiency (PCE) of organic solar cells (OSCs). Here, a facile strategy is presented to simultaneously enhance exciton/charge transport of the widely studied PM6:Y6-based OSCs by employing highly emissive trans-bis(dimesitylboron)stilbene (BBS) as a solid additive. BBS transforms the emissive sites from a more H-type aggregate into a more J-type aggregate, which benefits the resonance energy transfer for PM6 exciton diffusion and energy transfer from PM6 to Y6. Transient gated photoluminescence spectroscopy measurements indicate that addition of BBS improves the exciton diffusion coefficient of PM6 and the dissociation of PM6 excitons in the PM6:Y6:BBS film. Transient absorption spectroscopy measurements confirm faster charge generation in PM6:Y6:BBS. Moreover, BBS helps improve Y6 crystallization, and current-sensing atomic force microscopy characterization reveals an improved charge-carrier diffusion length in PM6:Y6:BBS. Finally, owing to the enhanced exciton diffusion, exciton dissociation, charge generation, and charge transport, as well as reduced charge recombination and energy loss, a higher PCE of 17.6% with simultaneously improved open-circuit voltage, short-circuit current density, and fill factor is achieved for the PM6:Y6:BBS devices compared to the devices without BBS (16.2%).

14 SOLAR ENERGY↗
DOE OSTI2024

Atomistic simulations to reveal HIP-bonding mechanisms of Al6061/Al6061

Molecular dynamics simulations were employed to understand the diffusion bonding process during hot isostatic pressing (HIP) of Al6061/Al6061 alloy. Simulations of the HIP process reveal atomistic phenomena that are difficult or unlikely to be observed experimentally and provide useful insights into the mechanism of diffusion and bonding. Here, the results reveal that at the start of the HIP process, a massive incursion of oxygen atoms occurs from the pre-existing γ-Al 2 O 3 to the 6061 region across the interphase interface. These oxygen atoms interact with the enriched Mg atom layer present at the existing γ-Al 2 O 3 and 6061 matrix to form a secondary complex Mg 2 Al 2 O 5 phase. Diffusion calculations also show that transport of atoms due to the applied pressure is 4–5 orders of magnitude higher than would occur in the absence of HIP conditions. The Mg 2 Al 2 O 5 phase also provides efficient pathways for the rapid transport of Mg atoms. Because of the higher diffusion coefficients observed for Mg within the phase, Mg atoms can move more swiftly compared to their diffusion within other phases such as γ-Al 2 O 3 . This accelerated mobility facilitates the rapid movement of Mg atoms across the interface, leading to changes in the local composition and the potential growth of the Mg 2 Al 2 O 5 phase.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Physics-Informed Machine Learning Models for Predicting the Progress of Reactive-Mixing

This paper presents a physics-informed machine learning (ML) framework to construct reduced-order models (ROMs) for reactive-transport quantities of interest (QoIs) based on high-fidelity numerical simu-lations. QoIs include species decay, product yield, and degree of mixing. The ROMs for QoIs are applied to quantify and understand how the chemical species evolve over time. First, high-resolution datasets for constructing ROMs are generated by solving anisotropic reaction-di?usion equations using a non-negative finite element formulation for di?erent input parameters. The reactive-mixing model input parameters are: time-scale associated with flipping of velocity, spatial-scale controlling small/large vortex structures of velocity, perturbation parameter of the vortex-based velocity, anisotropic dispersion strength/contrast, and molecular diffusion. Second, random forests, F-test, and mutual information criterion are used to evaluate the importance of model inputs/features with respect to QoIs. We observed that anisotropic dispersion strength/contrast is the most important feature and time-scale associated with flipping of velocity is the least important feature. Third, Support Vector Machines (SVM) and Support Vector Regression (SVR) are used to construct ROMs based on the model inputs. The constructed SVR-ROMs are then used to predict scaling of QoIs. We also present estimates and inequalities on the QoIs, which inform that the species decay, mix, and produce in an exponential fashion. These inequalities also inform that a radial basis function is the most suitable kernel for the SVM/SVR models for QoIs. It is observed that R2-score for SVR-ROMs on unseen data is greater than 0.9, implying that the SVR-ROMs are able to predict the reaction-diffusion system state reasonably well. Finally, in terms of the computational cost, the proposed SVM-ROMs are O(107) times faster than running a high-fidelity finite element simulation for evaluating QoIs. This makes the proposed ML-based ROMs attractive for reactive-transport sensing and real-time monitoring applications as they are significantly faster yet reasonably accurate.

Mudunuru, Maruti K.↗
DOE OSTI2024

Integration and validation of some modules for modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures

Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.

42 ENGINEERING↗
DOE OSTI2023

Methane mass transfer in mesoporous silica saturated with liquid hydrocarbons

Mass transfer across gas/liquid interfaces plays a central role in many industrial applications. In particular, gas dissolution and diffusion in liquid hydrocarbon mixtures, confined in nanometer-sized pores, is an essential mechanism during enhanced oil recovery (EOR) from unconventional formations. In this work, we have measured methane (C 1 ) diffusion in n-decane (C 10 ), n-hexadecane (C 16 ), and mixtures of C 10 + C 16 in a mesoporous material with an average pore size of 4 nm at 50 °C and ~ 8 MPa of gas pressure. A key conclusion of this work is that the diffusivities, measured in the bulk phase, for the relevant binary systems, are sufficient to predict the diffusion behavior of the corresponding multicomponent systems in the porous medium by using Wilke’s equation for the evaluation of effective component diffusivities, combined with an accurate equation of state (EOS) representation of the (bulk) phase behavior. This observation can facilitate accurate prediction of recovery processes in unconventional formations and, potentially, guide other applications that entail gas-liquid interface mass transfer in mesoporous materials.

42 ENGINEERING↗
DOE OSTI2021

Density functional theory calculations of the thermodynamic and kinetic properties of point defects in β -U

Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2 x 2 x 2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. Finally, the diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Co-transport of water and p-xylene through carbon molecular sieve membranes

Carbon molecular sieve (CMS) materials are a potential candidate for scalable and high-performance reverse osmosis membranes due to their impressive chemical and thermal stabilities. Moreover, they have the potential to enable impressive rejections of small neutral solutes from water based on their known ability to separate small organic molecules. CMS has been extensively examined for gas and organic solvent separations, but the transport of organic and aqueous mixtures through CMS microstructures is poorly understood. In this work, we investigated the sorption, diffusion, and permeation behavior of organic compounds and water in poly(vinylidene fluoride)(PVDF)-derived CMS (PVDF-CMS). Experimental observations of diffusion, sorption, and permeation shows how the properties of penetrants such as polarity and molecular size affect the transport rates and selectivity. These basic transport and sorption parameters are utilized in sorption-diffusion models to permeation rates of water-organic mixtures in CMS membranes. The transport of water and p-xylene in CMS was experimentally confirmed to follow the sorption-diffusion mechanism. The sorption-diffusion model ideal permselectivity indicates that the CMS is p-xylene selective over water. Water/p-xylene mixture permeation experiments revealed an increased selectivity of p-xylene over water, thus providing tentative evidence for a competitive sorption-selective separation mechanism. Furthermore, this work suggests that CMS membranes exhibit organic-permeable separation properties in water/organic separations. The results presented here highlight the potential for the removal of dilute organics in water via CMS pervaporation membranes.

02 PETROLEUM↗
DOE OSTI2022

Impact of carbonation on leaching of constituents from a cementitious waste form for treatment of low activity waste at the DOE Hanford site

Carbonation can be a major aging process during disposal of alkaline cementitious waste forms and can impact constituent leaching by changes in material alkalinity, pore structure, and controlling mineral phases. The effect of carbonation on the leaching of major and trace constituents from Cast Stone, a cementitious waste form developed to treat high salt content low activity waste, was studied through a combination of leaching experiments and reactive transport simulations. Diffusive transport of constituents in the waste form was evaluated using reactive transport modeling of diffusion-controlled leaching test results and a geochemical speciation model derived from pH-dependent leaching. Comparisons between Cast Stone materials aged under nitrogen, air, and 2% carbon dioxide in nitrogen showed that carbonation impacts solubility, physical retention and observed diffusivity of major and trace constituents. Carbonation under 2% CO 2 decreased the diffusion-controlled leaching of chromium by two orders of magnitude. Finally, modeling results suggest that carbonation may also decrease solubility of technetium while changes to microstructure by carbonation increases effective diffusivity of constituents in Cast Stone.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Oxygen and Proton Transport in Flooded Graphene Pores with N-Dopants and Defects

Reactant transport is an important consideration in the design of ideal electrode structures. For the oxygen reduction reaction catalyzed by Pt/C in proton exchange membrane fuel cell cathodes, O 2 and H + delivery to Pt surfaces and H 2 O transport away are required. Some Pt nanoparticles may only be accessible via micropores that are too small for ionomer molecules to enter, possibly requiring flooding for H + transport. Here, to test if these “buried” Pt particles can play a role in activity through this proposed transport mechanism, we have performed atomic-scale simulations based on reactive force field molecular dynamics with an emphasis on confinement below 20 Å. Diffusion coefficients as a function of the molar concentration and local environment have been evaluated in water domains confined in two-dimensional graphene nanochannels of various channel heights representing a morphological model for micropores in proton exchange fuel cell cathodes. Our study shows that local atomic-scale structures can strongly modify H + , O 2 , and H 2 O transport rates in flooded micropores less than 20 Å in size. We find that there is a critical crossover in diffusion behavior around the 20 Å spacing with larger pores having bulk-like diffusion properties and confinement below 20 Å monotonically decreases diffusion rates. As pore size decreases, we observe locally dispersed water regions that ultimately strand reactants from long-distance transport. These findings suggest that flooded micropores may in fact be viable transport pathways for relevant reactants and products if the pore walls on opposite sides remain separated by ≥10 Å separation. Furthermore, the confinement effect is so strong that N-doping and C-vacancy defects in the C pore wall have only a minimal impact on diffusion rates and their effects are counterintuitively more apparent at larger spacings. These findings provide valuable insight regarding cathode performance and the role “stranded” catalyst particles may play in fuel cell cathodes.

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

Kinetic pathways of fast lithium transport in solid electrolyte interphases with discrete inorganic components

The transport of lithium ions in the solid electrolyte interphase (SEI) has been previously accepted to proceed in two steps: a fast pore diffusion through the outer, porous organic layer followed by a slow knock-off or vacancy diffusion in the inner, dense inorganic layer. The second step is believed to be the rate-limiting step during fast-charging. In this study, we have intentionally constructed a thicker SEI (SEI-rich) structure on the surface of monoclinic Nb 2 O 5 (H-Nb 2 O 5 ) by adding LiNO 3 into a conventional ethylene carbonate based electrolyte. The electrochemical performance of two electrodes, one SEI-rich and one with few SEI (SEI-lean), was found to be almost the same, including their fast-charging capability and cycling stability, despite the significant difference in their SEI structure. Importantly, analysis using cryogenic scanning/transmission electron microscopy showed the discrete decoration of individual inorganic particles (e.g., Li 2 O) and amorphous species (LiN x O y /organic components) over the surface of H-Nb 2 O 5 . These discrete inorganic particles are in contradiction to the formation of dense inner inorganic layer, which has been commonly postulated. Based on these findings, we propose a new mechanism for Li ion transport through the SEI: one-step pore diffusion, without the second step slow diffusion. This one-step pore diffusion process provides an extremely fast Li ion transport, and effectively removes the kinetic limitation of Li ion transport in the SEI for fast charging. Finally, these results strongly suggest that the influence of SEI structure on the transport kinetics of lithium ions is much less significant than previously accepted. These results offer a new understanding of possible lithium ion transport pathway within SEI and may have implications for the future designs of fast-charging battery materials.

25 ENERGY STORAGE↗
DOE OSTI2020

Mixtures of octanol and an ionic liquid: Structure and transport

Ionic liquids (ILs) with long alkyl substituents are amphiphilic, which leads to a bicontinuous liquid structure. The strongly interacting anionic and cationic head groups form a long range charge network, with the hydrocarbon tails forming a nonpolar domain. Such nonpolar domains have been shown to dissolve a variety of neutral organic solvents. In mixtures of ILs with solvents the neutral organic molecules residing in the nonpolar domains experience different environments and friction from the charged cations and anions. Thus, the neutral molecules diffuse much faster than predicted by hydrodynamic scaling using the average viscosity of the mixture. Here, we report studies on the structure and transport properties of mixtures of 1-octanol with the IL trihexyltetradecylphosphonium bis(trifluoromethylsulfonyl)imide (P 6,6,6,14 + /NTf 2 – ). The majority of the atom fraction in the P 6,6,6,14 + cation comprises four hydrocarbon substituents. The unique amphiphilic nature of ILs with the P 6,6,6,14 + cation makes 1-octanol fully miscible with the IL at ambient temperatures. X-ray scattering experiments show that the IL structure persists in the mixtures for 1-octanol mole fractions as large as x oct = 0.90. The self-diffusion coefficients of the three molecular species in the mixtures were measured by NMR experiments. The self-diffusion of the P 6,6,6,14 + cation is well described by the Stokes–Einstein equation, while the diffusivity of the NTf 2 – anion is slightly lower than the hydrodynamic prediction. The measured diffusivities of octanol in these mixtures are 1.3–4 times higher than the hydrodynamic predictions.

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

Molecular dynamics simulation studies of 1,3-dimethyl imidazolium nitrate ionic liquid with water

The fundamental understanding of intermolecular interactions of ionic liquids (ILs) with water is essential in predicting IL–water thermodynamic properties. In this study, intermolecular or noncovalent interactions were studied for 1,3-dimethyl imidazolium [DMIM]+ cation and nitrate [NO 3 ] – anion with water, employing quantum mechanics and molecular dynamics simulations. Molecular dynamics simulations were performed using a revised multipolar polarizable force field. The effect of water on ionic liquids was evaluated in terms of thermodynamic and dynamic properties. Thermodynamic properties included liquid densities ρ, excess molar volumes ΔV E , and liquid structures g(r)⁠. Dynamic properties included self-diffusion coefficients D of mixture constituents as a function of water concentration. The density of ionic liquid–water mixtures monotonically decrease with increasing concentration of water. A negative excess volume was obtained for low and high water concentrations, demonstrating strong intermolecular interactions of water with ionic liquid components. Liquid structures of ionic liquid–water mixtures revealed a tendency for anions to interact with cations at shorter intermolecular distances when the water concentration is increased. Diffusion rates were found to increase for all mixture components with increase in the fraction of water. A significant change in the diffusion rate was found at ~0.3 weight fraction of water. However, the water self-diffusion coefficient was dominant at all concentrations. Furthermore, the ratio of water/anion and anion/cation self-diffusion coefficients was found to decrease linearly with increasing concentration of water molecules.

36 MATERIALS SCIENCE↗