Search NASA⌕ Search

SEARCH · Search NASA

Results for “Extended Storage and Transportation”

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.

126 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Localized high-concentration electrolytes get more localized through micelle-like structures

Liquid electrolytes in batteries are typically treated as macroscopically homogeneous ionic transport media despite having a complex chemical composition and atomistic solvation structures, leaving a knowledge gap of the microstructural characteristics. Here, we reveal a unique micelle-like structure in a localized high-concentration electrolyte, in which the solvent acts as a surfactant between an insoluble salt in a diluent. The miscibility of the solvent with the diluent and simultaneous solubility of the salt results in a micelle-like structure with a smeared interface and an increased salt concentration at the centre of the salt–solvent clusters that extends the salt solubility. These intermingling miscibility effects have temperature dependencies, wherein a typical localized high-concentration electrolyte peaks in localized cluster salt concentration near room temperature and is used to form a stable solid–electrolyte interphase on a Li metal anode. Importantly, these findings serve as a guide to predicting a stable ternary phase diagram and connecting the electrolyte microstructure with electrolyte formulation and formation protocols of solid–electrolyte interphases for enhanced battery cyclability.

25 ENERGY STORAGE↗
DOE OSTI2020

Wettability of ultra-small pores of carbon electrodes by size-asymmetric ionic fluids

Recently, we studied the phase behavior of ionic fluids under confinement using the classical density functional theory within the framework of the restricted primitive model. Additionally, the theoretical results indicate that narrowing the pore size may lead to a drastic reduction in the electric double layer capacitance, while increasing the surface electrical potential would improve the ionic accessibility of micropores. In this work, we extend the theoretical investigation to systems containing size-asymmetric electrolytes that may exhibit a vapor-liquid like phase transition in the bulk phase. The effects of pore size and surface electric potential on the phase diagram and microscopic structures of the confined electrolytes were studied over a broad range of parameters. We found that decreasing the pore size or increasing the surface potential could destabilize the liquid phase in micropores, and capillary evaporation could occur regardless of the size asymmetry between cations and anions. Compared to that in a symmetric ionic system, the vapor-liquid phase separation is more likely to take place as the size asymmetry becomes more pronounced. The phase transition would alter the “accessibility” of ions to micropores and lead to coexisting micropores with different surface charge densities as identified by Monte Carlo simulation.

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

Nonlinear electrical transport near the metal–insulator transition in V 4 O 7 thin films

Nonlinear electrical transport associated with correlated electronic states has been widely investigated in transition-metal oxides near metal–insulator transitions. Here, in this study, we investigate nonlinear transport and threshold switching in sputter-deposited V 4 O 7 thin films grown on fused silica substrates. Temperature-dependent transport measurements show a metal–insulator transition near 240 K with negligible thermal hysteresis, defining the temperature scale that governs the nonlinear electrical response. Current–voltage measurements reveal reproducible bipolar threshold switching over a broad temperature range, with the threshold voltage decreasing systematically as the transition temperature is approached. Electrothermal finite-element simulations reproduce the measured switching characteristics and show that switching occurs when localized Joule heating drives a confined region of the device toward the transition temperature. The simulated peak local temperature at threshold lies close to the intrinsic metal–insulator transition temperature of V 4 O 7 , establishing that the nonlinear electrical response originates from electrothermal feedback acting on the strongly temperature-dependent conductivity near the transition. These results identify V 4 O 7 thin films as a model correlated-oxide system in which volatile threshold switching emerges from proximity to an extended metal–insulator transition.

25 ENERGY STORAGE↗
DOE OSTI2024

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Lithium-compatible and air-stable vacancy-rich Li 9 N 2 Cl 3 for high–areal capacity, long-cycling all–solid-state lithium metal batteries

Attaining substantial areal capacity (>3 mAh/cm 2 ) and extended cycle longevity in all–solid-state lithium metal batteries necessitates the implementation of solid-state electrolytes (SSEs) capable of withstanding elevated critical current densities and capacities. In this study, we report a high-performing vacancy-rich Li 9 N 2 Cl 3 SSE demonstrating excellent lithium compatibility and atmospheric stability and enabling high–areal capacity, long-lasting all–solid-state lithium metal batteries. The Li 9 N 2 Cl 3 facilitates efficient lithium-ion transport due to its disordered lattice structure and presence of vacancies. Notably, it resists dendrite formation at 10 mA/cm 2 and 10 mAh/cm 2 due to its intrinsic lithium metal stability. Furthermore, it exhibits robust dry-air stability. Incorporating this SSE in Ni-rich LiNi 0.83 Co 0.11 Mn 0.06 O 2 cathode-based all–solid-state batteries, we achieve substantial cycling stability (90.35% capacity retention over 1500 cycles at 0.5 C) and high areal capacity (4.8 mAh/cm 2 in pouch cells). These findings pave the way for lithium metal batteries to meet electric vehicle performance demands.

25 ENERGY STORAGE↗
DOE OSTI2025

Revealing the roles of the solid–electrolyte interphase in designing stable, fast-charging, low-temperature Li-ion batteries

Designing the solid–electrolyte interphase (SEI) is critical for stable, fast-charging, low-temperature Li-ion batteries. Fostering a “fluorinated interphase,” SEI enriched with LiF, has become a popular design strategy. Although LiF possesses low Li-ion conductivity, many studies have reported favorable battery performance with fluorinated SEIs. Such a contradiction suggests that optimizing SEI must extend beyond chemical composition design to consider spatial distributions of different chemical species. In this work, we demonstrate that the impact of a fluorinated SEI on battery performance should be evaluated on a case-by-case basis. Sufficiently passivating the anode surface without impeding Li-ion transport is key. We reveal that a fluorinated SEI containing excessive and dense LiF severely impedes Li-ion transport. In contrast, a fluorinated SEI with well-dispersed LiF (i.e., small LiF aggregates well mixed with other SEI components) is advantageous, presumably due to the enhanced Li-ion transport across heterointerfaces between LiF and other SEI components. An electrolyte, 1 M LiPF 6 in 2-methyl tetrahydrofuran (2MeTHF), yields a fluorinated SEI with dispersed LiF. This electrolyte allows anodes of graphite, μSi/graphite composite, and pure Si to all deliver a stable Coulombic efficiency of 99.9% and excellent rate capability at low temperatures. Pouch cells containing layered cathodes also demonstrate impressive cycling stability over 1,000 cycles and exceptional rate capability down to −20 °C. Through experiments and theoretical modeling, we have identified a balanced SEI-based approach that achieves stable, fast-charging, low-temperature Li-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2026

Site-Selective Modification of Lanthanum Oxychloride to Modulate Halide-Ion Conduction

Design principles for solid-state halide-ion conduction remain poorly defined despite the increasing importance of halide ions as charge carriers in a variety of energy storage and electrochemical computing technologies. Here, we employ a siteselective modification strategy in which aliovalent cations are preferentially introduced at the La 3+ crystallographic site of LaOCl in the 2c Wyckoff position, enabling controlled generation of chloride vacancies and modification of lattice dynamics to enhance chloride-ion conductivity. Aliovalent substitution of La 3+ with Mg 2+ , Ca 2+ , and Sr 2+ generates charge-compensating Cl vacancies while preserving the matlockite crystal structure. X-ray excited optical luminescence measurements with Dy 3+ as a reporter chromophore evidence vacancy-derived midgap electronic states and an extended energy range of radiation-less Auger electron emission corresponding to substantial modification of electronic structure and local electrostatic potentials. Ca alloying at 8−10 at. % increases the chloride-ion conductivity by three- to 4 orders of magnitude as compared to unalloyed LaOCl, whereas comparable amounts of Sr- and Mg-alloying in LaOCl imbue less pronounced conductivity enhancements. Temperature-dependent Raman spectroscopy measurements reveal that Ca- and Sr-alloying substantially soften the La−Cl sublattice and yield a more compliant crystal lattice that can deform to accommodate Cl-ion migration. Structure solutions derived from Rietveld refinements to powder Xray diffraction reveal larger O−La−Cl bond-angle deviations and enhanced out-of-plane cation displacements for Ca- and Sr-alloyed compositions as compared to Mg-alloyed LaOCl. Such local distortions enhance chloride-ion mobility by reshaping and flattening vacancy migration energy landscapes and by modulating lattice dynamics governing anion conduction. We further illustrate that coalloying of Ca with Mg and Sr induces a nonmonotonic conductivity−defect stoichiometry relationship that can be rationalized based on cooperative interactions. Together, these results establish site-selective aliovalent alloying of LaOCl as an effective route to halide-ion solid electrolytes and provide broadly generalizable design principles for site-selective modification to induce vacancy formation and lattice softening to engender facile anion transport

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

Thermal Dehydration of Aluminum (Oxy)hydroxides on Fuel Cladding Material - 20200

The aluminum cladding of research-reactor fuel undergoes general corrosion with resulting formation of adherent aluminum (oxy)hydroxide films during in-reactor and post-discharge exposure to water under various conditions and temperatures. These (oxy)hydroxides contain chemically-bound water that poses challenges for extended dry storage due to the risk of thermal or radiolytic decomposition releasing free water and/or hydrogen and oxygen gases. This study comprised laboratory experiments of thermal drying behavior of aluminum (oxy)hydroxides to identify approaches for reducing bound water on fuel cladding prior to sealed dry storage. Effective drying strategies for adherent (oxy)hydroxides on fuel cladding will improve the safety of dry storage by mitigating potential avenues for additional corrosion and/or generation of flammable gases inside the storage canister. This work is part of a broader investigation to address knowledge gaps and technical data needs for dry storage of aluminum-clad spent nuclear fuel (ASNF), which included investigation of (oxy)hydroxide formation on aluminum alloy substrates immersed in water, characterization of service-grown films from ASNF, and radiolytic yield of hydrogen from (oxy)hydroxide powders and films. The current work comprises experimental thermal dehydration of aluminum trihydroxides (bayerite and gibbsite) characteristic of low-temperature (<80 deg. C) corrosion in water and aluminum oxyhydroxide (boehmite) characteristic of high temperature (>80 deg. C) corrosion in water. First, commercially produced (oxy)hydroxide powders (boehmite and gibbsite) were tested via thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC). The effort aimed to identify or confirm key temperature ranges/thresholds for the thermal decomposition reactions, as well as impacts of ramp rate and hold times, for drying of isolated, high-surface-area (oxy)hydroxides. The information gleaned from powder tests was used to guide subsequent drying tests on adherent (oxy)hydroxide films grown on aluminum alloy substrates. TGA was used to analyze small samples of the aluminum coupons with adherent (oxy)hydroxide films. Specimens were characterized both pre- and post-drying using X-ray diffraction (XRD) to determine the film composition and scanning electron microscopy (SEM) to determine its morphology. Drying of adherent (oxy)hydroxide films was anticipated to pose additional complications relative to drying of (oxy)hydroxide powders. Anticipated challenges include achieving sufficient drying of the (oxy)hydroxide at low-enough temperature to avoid melting or other undesirable phase changes in the aluminum substrate and overcoming transport limitations of water through the thickness of a dense film or through a tortuous pore structure to reach the outer surface of the film. Inn addition, dehydration of (oxy)hydroxides significantly changes the film density, which may lead to substantial alterations in the morphology, including potential cracking and spalling of the film. TGA/DSC tests of gibbsite powders resulted in successful conversion to boehmite or alumina, depending on the maximum temperature reached. At low ramp rates (≤5 deg. C/min), the conversion to boehmite occurred around 300 deg. C (210-340 deg. C). XRD confirmed that boehmite was the only phase detected after tests reaching 450 deg. C. The transition to alumina occurred around 510 deg. C (470-550 deg. C), with XRD detecting only alumina after tests reaching at least 600 deg. C. Boehmite powders dehydrated to alumina at about 400 deg. C (330-460 deg. C) for coarse (77-μm particle diameter) powder and about 490 deg. C (420-520 deg. C) for fine (0.7-μm particle diameter) powder. For all powders tested, the maximum percent mass loss after drying to high temperature slightly exceeded the theoretical mass loss for complete dehydration of the stoichiometric (oxy)hydroxide to alumina, which is likely attributable to physisorbed water in the powder. Drying of a thick (∼8.6 μm), predominantly bayerite adherent film displayed a dramatic change in film morphology, with the initially continuous trihydroxide film cracking into sections on the order of 50 μm square and partially delaminating. The outermost layer of the film completely faked of in some regions. The layer exposed under the spalled oxide also displayed prominent cracking, with spacing on the order of 10 μm, but this surface layer appeared to remain completely adhered to the aluminum substrate. TGA showed mass losses per unit surface area up to ∼1.0 mg/cm{sup 2}. XRD characterization of the remaining oxide was inconclusive, with no crystalline phases detected. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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 OSTI2020

Charge Transport in Zirconium-Based Metal–Organic Frameworks

Metal–organic frameworks (MOFs) are a class of crystalline porous materials characterized by inorganic nodes and multitopic organic linkers. Because of their molecular-scale porosity and periodic intraframework chemical functionality, MOFs are attractive scaffolds for supporting and/or organizing catalysts, photocatalysts, chemical-sensing elements, small enzymes, and numerous other functional-property-imparting, nanometer-scale objects. Notably, these objects can be installed after the synthesis of the MOF, eliminating the need for chemical and thermal compatibility of the objects with the synthesis milieu. Thus, postsynthetically functionalized MOFs can present three-dimensional arrays of high-density, yet well-separated, active sites. Depending on the application and corresponding morphological requirements, MOF materials can be prepared in thin-film form, pelletized form, isolated single-crystal form, polycrystalline powder form, mixed-matrix membrane form, or other forms. For certain applications, most obviously catalytic hydrolysis and electro- or photocatalytic water splitting, but also many others, an additional requirement is water stability. MOFs featuring hexa-zirconium(IV)-oxy nodes satisfy this requirement. For applications involving electrocatalysis, charge storage, photoelectrochemical energy conversion, and chemiresistive sensing, a further requirement is electrical conductivity, as embodied in electron or hole transport. As most MOFs, under most conditions, are electrically insulating, imparting controllable charge-transport behavior is both a chemically intriguing and chemically compelling challenge. Herein, we describe three strategies to render zirconium-based metal–organic frameworks (MOFs) tunably electrically conductive and, therefore, capable of transporting charge on the few nanometers (i.e., several molecular units) to few micrometers (i.e., typical dimensions for MOF microcrystallites) scale. The first strategy centers on redox-hopping between periodically arranged, chemically equivalent sites, essentially repetitive electron (or hole) self-exchange. Zirconium nodes are electrically insulating, but they can function as grafting sites for (a) redox-active inorganic clusters or (b) molecular redox couples. Alternatively, charge hopping based on linker redox properties can be exploited. Marcus’s theory of electron transfer has proven useful for understanding/predicting trends in redox-hopping based conductivity, most notably, in accounting for variations as great as 3000-fold depending on the direction of charge propagation through structurally anisotropic MOFs. In MOF environments, propagation of electronic charge via redox hopping is necessarily accompanied by movement of charge-compensating ions. Consequently, rates of redox hopping can depend on both the identity and concentration of ions permeating the MOF. In the context of electrocatalysis, an important goal is to transport electronic charge fast enough to match or exceed the inherent activity of MOF-based or MOF-immobilized catalysts. Bandlike electronic conductivity is the focus of an alternative strategy: one based on the introduction of molecular guests capable of forming donor–acceptor charge transfer complexes with the host framework. Theory again can be applied predictively to alter conductivity. A third strategy similarly emphasizes electronic conductivity, but it makes use of added bridges in the form of molecular oligomers or inorganic clusters that can then be linked to span the length of a MOF crystallite. For all strategies, retention of molecular-scale porosity is emphasized, as this property is key to many applications. Lastly, while our focus is on Zr-MOFs, the described approaches clearly are extendable to other MOF compositions, as has already been demonstrated, in part, in studies by others.

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

Localized High-Concentration Electrolytes Get More Localized

Liquid electrolytes in batteries are typically treated as macroscopically homogeneous ionic transport media, despite having complex chemical composition and atomistic solvation structures, thus leaving a knowledge gap of microstructural characteristics. A promising electrolyte system design has recently been developed with the title of a localized high-concentration electrolyte (LHCE). These electrolytes are a summation of three components – a charge-carrying salt, a solvent that solvates the salt, and a diluent that is sparingly interactive with the salt in terms of ion solvation. This creates the perception of a high-concentration electrolyte (HCE) with high anion-cation interactions that supports a salt-derived, stable solid-electrolyte interphase, while an LHCE improves bulk electrolyte properties, such as ionic conductivity and viscosity, when compared to an HCE. Much trial by error design of this LHCE concept has been supported by simulations to confirm ionic interactions, cycling data to confirm cell stability, and post-mortem analysis to provide evidence of mechanisms. However, a systematic design to extend the capabilities of LHCEs is still to be desired. Here, we analyzed a specific LHCE mixture to dictate the best viable mixture that extends the localization of solvation structures while maintaining viable cell-level capabilities. We seek to observe and understand solvent-to-salt interactions that impact electrolyte performance. These interactions range from single-ion solvation, multidentate interactions (i.e., solvation-driven clustering), to more widespread network-level interactions. With the support of both experimental and computational analysis of different electrolyte features, we have begun a methodology of formulating the most ideal LHCE mixtures for improving the cyclability of high-capacity battery chemistries.

25 ENERGY STORAGE↗
DOE OSTI2025

Elucidating ion capture and transport mechanisms of Preyssler anions in aqueous solutions using biased MACE-accelerated MD simulations

Equilibrium and biased multi-atomic cluster expansion (MACE) accelerated molecular dynamics (MD) simulations in aqueous solutions are performed to investigate the ion capture and transport mechanisms of the {P 5 W 30 } Preyssler anion (PA) as the smallest representative member of the extended polyoxometalate (POM) family with an internal cavity. The unique interatomic interactions present in the internal cavity vs the exterior of PA are carefully investigated using equilibrium MACE MD simulations for two representative Na(H 2 O)@PA and Na@PA complexes in aqueous solutions. Our careful analyses of radial distribution functions and coordination numbers show that the presence of confined water in Na(H 2 O)@PA has profound modulating effects on the nature of the interactions of the encapsulated ion with the oxygens of the PA cavity. Using well-converged MACE-accelerated multiple walker well-tempered metadynamics simulations with nanosecond timescales, two different associative (ion exchange) and dissociative (ion ejection) ion transport mechanisms were carefully investigated for Na + as one of the most abundant and representative ions present in seawater and saline solutions. By comparing systems with and without confined water, it was found that the presence of only one pre-encapsulated confined water in Na(H 2 O)@PA dramatically changes the free energy landscape of ion transport processes. It was also found that the contraction and dilation of the two windows present in PA directly influence the Na + and H 2 O transport. Furthermore, the results from this work are helpful, as they show a viable path toward tuning the ion exchange and transport phenomena in aqueous solutions of POM molecular clusters and frameworks.

25 ENERGY STORAGE↗
DOE OSTI2022

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

3D Morphology of Bimodal Porous Copper with Nano-Sized and Micron-Sized Pores to Enhance Transport Properties for Functional Applications

Multiscale porous metals with multiscale porosity from nanometer to micrometer have a high specific surface area and high effective diffusivity for ion transport, thereby enhancing functionalities and extending the applications of porous metals. Herein, the Cu–Fe–Al ternary system was selected as the precursor alloy to construct multiscale, bimodal porous copper by the chemical dealloying method. The effect of the phase composition and initial microstructure of precursor alloys (Al x Fe 75– x Cu 25 , x = 10–60) on the three-dimensional (3D) morphology of multiscale porous metals was systematically investigated, with a goal to precisely control the multiscale porous structure. The four crystal structure phases (body-centered cubic (BCC), face-centered cubic (FCC), CsCl type (B2), and monoclinic) in precursor alloys were analyzed by synchrotron X-ray diffraction refinement. The 3D morphology, feature size distribution, and tortuosity of four representative porous Cu after dealloying Al x Fe 75– x Cu 25 ( x = 10, 30, 50, and 60) precursor alloys were directly visualized and quantified via advanced synchrotron X-ray nanotomography. The relationship between the phases/crystal structures of precursor alloys and their corresponding porous morphology was established: the micron-sized pores in bimodal porous Cu are formed by dissolving the CuFeAl phase with BCC and monoclinic crystal structures, and the nano-sized pores are formed by dealloying the CuFeAl phase with FCC and B2 crystal structures. The size of the nanoporous structure depends on the ratio between the more noble and more active components in the precursor alloy, while the size of the microporous structure depends on the corresponding phase size in the precursor alloy. The tortuosity results showed that the multiscale porous structure with both nanoporosity and microporosity exhibits lower tortuosity, which will enhance transport properties for functional applications.

25 ENERGY STORAGE↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
DOE OSTI2025

Probing the Sustainable Reduction of CO 2 , N 2 and NO 3 to Fuels and Chemicals using Non-Traditional Porphyrinoid Catalysts

For the last several decades, numerous strategies have been proposed to be able to interconvert electricity and commodity fuels for transportation and other applications. Even as solar and wind energy become more widely available, these renewable energy sources are not always available where the population is most dense or during peak times of energy demand. One plausible solution to the issue of electricity storage and transport is to use electricity to drive the formation of high-energy compounds and chemical fuels. As a result, the electrochemical conversion of CO 2 into chemical fuels has received major attention as a multi-pronged approach for electricity coversion, storage, and transport. Similar strategies are also attractive for conversion of N 2 and NO x into value added compounds such as ammonia. Metalloporphyrin and metallocorrole complexes are comprised of aromatic tetrapyrrole scaffolds and have been extensively studied as homogeneous catalysts for critical catalytic processes like the electrochemical CO 2 reduction reaction (eCO 2 RR) to generate CO as feedstock for the Fischer–Tropsch process to generate a variety of hydrocarbons as chemical fuels. Despite their significance and impact, efficiently driving such catalytic reduction processes is challenging for multiple reasons. Such processes require multiple-electron transfer events, as opposed to single-electron redox chemistry that forms high-energy singly-reduced intermediates. In addition, these electron transfer events must be coupled to proton-transfer steps and avoid application of high overpotentials where proton (H + ) reduction to generate H 2 (as well as other side-reactions) can lower the selectivity and energy efficiency of the eCO 2 RR process. Certain Fe(III) aromatic tetrapyrroles (e.g., porphyrins and corroles) have been reported to overcome these challenges, but are often difficult to synthesize and modify, and have a tendency to support single-electron redox chemistry at the metal as opposed to multielectron redox involving the ligand and metal in concert. Other families of tetrapyrroles in which all four meso-carbons are reduced (i.e., sp 3 -hybridized) are also known in the literature. These non-aromatic tetrapyrroles are known as porphyrinogens and support multi-electron redox chemistry which has been elaborated by several groups, prompting researchers to consider whether these scaffolds may improve the kinetics and efficiencies of multi-electron steps attendant to activation of thermodynamically stable small-molecule substrates such as CO 2 , N 2 , NO 3 − , and NO 2 − . Although porphyrinogens support multi-electron redox properties, the four sp 3 -hybridized meso-carbons inherent to the tetrapyrrole core completely disrupt π-conjugation between pyrrolic units, which compromises the photochemical properties of porphyrinogens in comparison to that of traditional porphyrinoids. Moreover, the highly reducing nature of porphyrinogens predisposes these platforms to extreme air and water sensitivity. Accordingly, such metalated porphyrinogens are often pyrophoric and/or incompatible with many common organic solvents (CH 2 Cl 2 , CHCl 3 , CH 3 CN, EtOAc, etc.). Hence, reports of efficient small-molecule activation or catalysis supported by porphyrinogens have been limited. To overcome the instability of porphyrinogens while retaining the scaffold’s attractive multi-electron redox properties, we have directed attention to establishing less-well studied groups of non-aromatic tetrapyrroles (i.e., phlorins, biladienes, and isocorroles). These non-traditional tetrapyrroles each contain just a single sp 3 -hybridized meso-carbon, which ablates the platforms aromaticity but still provides an extended π-framework. In addition to developing innovative platforms that may be used to sustainably generate value-added chemicals, fuels, and ammonia via either photo- or electrocatalytic approaches our work also has significantly broadened our understanding of how to prepare and modulate the properties of non-aromatic tetrapyrroles for other applications. Our efforts entailed a synergistic partnership and active collaborations with Drs. David C. Grills and Mehmed Z. Ertem (both of Brookhaven National Lab) to better characterize isocorroles (and related non-aromatic tetrapyrroles) using a variety of advanced spectroscopic and theoretical methods. Our combined efforts have shown that the combination of properties that isocorroles (and related tetrapyrroles) provide results in good chemical stability paired with unique redox and photochemical characteristics that are not typically supported by simple and/or unadorned aromatic tetrapyrroles. Beside redox chemistry, the photochemistry of isocorroles and other non-aromatic tetrapyrroles containing one sp 3 -hybridized meso-carbon is also appealing due to their absorption in the long-visible to near-IR regions, which are essential for photocatalysis and related applications that benefit from direct excitation at these wavelengths. The results reported vastly improve our understanding of non-aromatic tetrapyrrole synthesis, properties and utility for activation of small molecule substrates such as O 2 and CO 2 in the presence of weak cationic organic acids.

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

Atomically Dispersed Zn/Co–N–C as ORR Electrocatalysts for Alkaline Fuel Cells

Hydrogen fuel cells have drawn increasing attention as one of the most promising next-generation power sources for future automotive transportation. Developing efficient, durable, and low-cost electrocatalysts, to accelerate the sluggish oxygen reduction reaction (ORR) kinetics, is urgently needed to advance fuel cell technologies. Herein, we report on metal–organic frameworks-derived nonprecious dual metal single-atom catalysts (SACs) (Zn/Co–N–C), consisting of Co–N 4 and Zn–N 4 local structures. These catalysts exhibited superior ORR activity with a half-wave potential (E 1/2 ) of 0.938 V versus RHE (reversible hydrogen electrode) and robust stability (ΔE 1/2 = –8.5 mV) after 50k electrochemical cycles. Moreover, this remarkable performance was validated under realistic fuel cell working conditions, achieving a record-high peak power density of ~1 W cm –2 among the reported SACs for alkaline fuel cells. Operando X-ray absorption spectroscopy was conducted to identify the active sites and reveal catalytic mechanistic insights. The results indicated that the Co atom in the Co–N 4 structure was the main catalytically active center, where one axial oxygenated species binds to form an O ads –Co–N 4 moiety during the ORR. In addition, theoretical studies, based on a potential-dependent microkinetic model and core-level shift calculations, showed good agreement with the experimental results and provided insights into the bonding of oxygen species on Co–N 4 centers during the ORR. Furthermore, this work provides a comprehensive mechanistic understanding of the active sites in the Zn/Co–N–C catalysts and will pave the way for the future design and advancement of high-performance single-site electrocatalysts for fuel cells and other energy applications.

25 ENERGY STORAGE↗
DOE OSTI2021

Topical Report – Findings on Subtask 2.7 – Wet ESP and Aerosol Testing at Coal Creek Station

Growing concerns over the impact of CO 2 emissions from combustion sources on global climate change have prompted numerous research and development projects aimed at developing cost-effective technologies for CO 2 capture. One family of technologies being demonstrated at pilot and full scale globally is postcombustion carbon capture (PCCC) systems that employ amine-based solvents. The captured CO 2 can be compressed and permanently stored underground or used for enhanced oil recovery. The proximity of North Dakota’s lignite-fired fleet of power plants to potential CO 2 storage options creates a unique atmosphere for PCCC within the state. However, the unique components present in lignite flue gas present a challenge for large-scale PCCC at North Dakota power plants by contributing to aerosol formation. Aerosols can negatively impact the long-term performance of amine-based solvents for CO 2 capture. Amine-based solvents are volatile, and flue gas particulate provides nucleation sites where amine vapors can condense as aerosols. Because aerosols cannot be easily captured at the column outlet using conventional technologies, the amine-laden aerosols escape the system and lead to accelerated solvent losses. Moreover, particulate components can chemically react with amines to form degradation products that can permanently deactivate the amine, cause fouling, and lead to hazardous emissions. Many of the elements that have been shown to catalyze solvent degradation are present in lignite coals and can exacerbate solvent replacement economics. Understanding this issue is critical to the implementation of solvent-based CO 2 capture systems as applied to lignite-fired generation systems. The Energy & Environmental Research Center (EERC) designed and carried out this project to fully characterize aerosol behavior with various control technologies installed to better optimize aerosol mitigation technology for CO 2 capture. To meet the goal of this project, the following objectives were identified: Determine the effectiveness of a wet electrostatic precipitator (WESP) on mitigating formation of problematic aerosols at Great River Energy’s Coal Creek Station, upstream of the PCCC system. Determine the effectiveness of the Mitsubishi Heavy Industries (MHI) proprietary amine emission reduction unit (AERU) as a postcapture solvent recovery system for reducing aerosol emissions and extending solvent life downstream of the PCCC system. Determine the impact of aerosols on the efficiency and degradation products of both commercial and advanced solvents within the PCCC system. Work was conducted at Coal Creek Station Unit 1 using a slipstream of flue gas from the outlet of the plant’s flue gas desulfurization (FGD) unit. Flue gas was routed through a pilot-scale FGD unit to remove SO 2 to very low levels (~1 ppm) and then through a direct contact cooler (DCC) to further cool the gas and to remove moisture. The gas exiting the DCC was then optionally routed through a WESP before passing to the CO 2 absorber columns. The MHI solvent was used to scrub CO 2 from the slipstream through a set of two absorber columns. The rich solvent was regenerated in a stripper column by heating to drive off captured CO 2 . Flue gas exiting the absorber column was routed to MHI’s proprietary AERU to recover entrained solvent. The system operated using a catch-and-release method where the CO 2 was separated to provide data on the process, but the captured CO 2 was released back into the host site stack. Particulate was measured, collected, and analyzed from multiple locations throughout the pilot-scale system. Unlike the performance observed in prior work, the inlet FGD and DCC did not remove significant particulate matter from the flue gas. This appears to be due to a difference in the nature of the particulate. The DCC seemed to increase particulate size and count, most likely owing to water condensing onto the surfaces of fly ash particles. When the WESP was operated, it achieved >95% particulate capture. Very little particulate matter or indications of solvent were detected at the AERU outlet. When operating with advanced KS 21 solvent, the particulate material at the AERU outlet was even further decreased. Solvent analysis showed that some species derived from flue gas and ash were slowly concentrating in the solvent over the duration of the test. The levels observed were reported to be within expected ranges and were not of concern to MHI. A high-level techno-economic assessment of installing CO 2 capture at Coal Creek Station suggested that, when using a standard monoethanolamine (MEA)-based solution with simple heat integration, the energy penalty to net generation would be 34%. The bulk of this was due to steam losses for regenerating solvent, followed by parasitic electrical demand for CO 2 compression and then by increased parasitic load for pushing flue gas through the absorber column. These demands could be decreased with a more advanced solvent that exhibits lower heat of regeneration and lower pressure drop than does a simple MEA solution. Further energy could be saved with more thorough heat integration to recover useful energy from the steam used for solvent regeneration. Installing a WESP was predicted to increase the cost of electricity by nearly $5/MWh. This would become cost-effective if solvent losses were roughly 10 times the baseline estimate when not using a WESP but could be returned to baseline by installing the WESP. Piping CO 2 for storage in more favorable geology could help with carbon capture and storage economics. Although storing off-site would necessitate construction of a CO 2 transport pipeline, the cost of this pipeline might be more than offset by reducing the number of wells required, the depths of the wells required, and the electrical demand for the CO 2 compressor. Additional factors that favor off-site storage costs include smaller expected CO 2 plume sizes, which translates to less monitoring and fewer landowner agreements. More detailed assessment of the specific geology in the region under and around Coal Creek Station would be needed to accurately assess the costs and benefits of different storage site options. This subtask was cofunded through the Energy & Environmental Research Center–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission.

Strege, Joshua↗