Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “storage complex”

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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Catalytic Resonance Theory: Experimental and Kinetic Interpretation of Programmable Catalysis

The complexity of programmable catalysts that change in the physical or electronic state on the time scale of a catalytic turnover results in uncertainty in the kinetic outcome of common measurements with dynamic catalyst experiments. To provide an interpretation of utility to experimentalists, a general programmable catalytic reaction exhibiting high turnover efficiency was simulated to understand the kinetic behavior of catalysts with varying frequency, binding energy, oscillation amplitude, and minimum oscillation binding energy at different temperatures and reactant gas pressures. The time-averaged catalytic rate for varying temperature in the form of an Arrhenius plot exhibited three distinct kinetic regimes, with the intermediate temperatures or applied frequencies exhibiting zero slope, indicative of barrier-free kinetic control. Transitions in experimentally measurable kinetic regimes with both temperature and applied catalyst oscillation frequency were associated with transitions in the sensitivity of reaction rate constants and degrees of rate control. Furthermore, these distinct kinetic macroscopic features specific to programmable catalysts were identified for observation in experimental dynamic kinetic measurements.

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

Role of Sr doping and external strain on relieving bottleneck of oxygen diffusion in La 2- x Sr x CuO 4-δ

In many complex oxides, the oxygen vacancy formation is a promising route to modify the material properties such as a superconductivity and an oxygen diffusivity. Cation substitutions and external strain have been utilized to control the concentration and diffusion of oxygen vacancies, but the mechanisms behind the controls are not fully understood. Using first-principles calculations, we find how Sr doping and external strain greatly enhances the diffusivity of oxygen vacancies in La 2-x Sr x CuO 4-δ (LSCO) in the atomic level. In hole-doped case (2x > δ), the formation energy of an apical vacancy in the LaO layer is larger than its equatorial counterpart by 0.2 eV that the bottleneck of diffusion process is for oxygen vacancies to escape equatorial sites. Such an energy difference can be reduced and even reversed by either small strain (< 1.5%) or short-range attraction between Sr and oxygen vacancy, and in turn, the oxygen diffusivity is greatly enhanced. For fully compensated hole case (2x ≦ δ), the formation energy of an apical vacancy becomes too high that most oxygen vacancies cannot move but would be trapped at equatorial sites. From our electronic structure analysis, we found that the contrasting change in the formation energy by Sr doping and external strain is originated from the different localization natures of electron carrier from both types of oxygen vacancies.

25 ENERGY STORAGE↗
DOE OSTI2024

Uncertainty quantification and propagation in lithium-ion battery electrodes using bayesian convolutional neural networks

The complex nature of manufacturing processes stipulates electrodes to possess high variability with increased heterogeneity during production. X-ray computed tomography imaging has proved to be critical in visualizing the complicated stochastic particle distribution of as-manufactured electrodes in lithium-ion batteries. However, accurate prediction of their electrochemical performance necessitates precise evaluation of kinetic and transport properties from real electrodes. Image segmentation that characterizes voxels to particle/pore phase is often meticulous and fraught with subjectivity owing to a myriad of unconstrained choices and filter algorithms. Here we utilize a Bayesian convolutional neural network to tackle segmentation subjectivity and quantify its pertinent uncertainties. Otsu inter-variance and Blind/Referenceless Imaging Spatial Quality Evaluator are used to assess the relative image quality of grayscale tomograms, thus evaluating the uncertainty in the derived microstructural attributes. We analyze how image uncertainty is correlated with the uncertainties and magnitude of kinetic and transport properties of an electrode, further identifying pathways of uncertainty propagation within microstructural attributes. The coupled effect of spatial heterogeneity and microstructural anisotropy on the uncertainty quantification of transport parameters is also understood. This work demonstrates a novel methodology to extract microstructural descriptors from real electrode images through quantification of associated uncertainties and discerning the relative strength of their propagation, thus facilitating feedback to manufacturing processes from accurate image based electrochemical simulations.

25 ENERGY STORAGE↗
DOE OSTI2026

A fresh perspective on the role of band bending, and related contributors, in light-driven production of electricity and chemicals

It is widely known that semiconductor-based solar energy conversion could power our planet. This is in part because high-quality semiconductor structures are unrivalled in their ability to separate photogenerated electrons and holes. One effective approach to achieving this photoinduced charge separation relies on a phenomenon known as “band bending”. But details to justify why band bending results in photoinduced charge separation are more complex than often appreciated. This underappreciation is an impediment to the rational, hypothesis-driven design of next-generation approaches to solar energy conversion. Herein we show, by means of derivations rooted in physical chemistry, that several phenomena – not just band bending – can facilitate photoinduced charge separation, and that each is influenced by nonequilibrium species concentration and a parameter, such as diffusion coefficient or rate coefficient, that introduces dynamics. To help visualize the impact of each phenomenon, we introduce plots that depict their contributions as free energy, force, flux, force constant, and rate. We reveal that spatial dopant distributions that define band bending are predictors of initial photogenerated species transport rates. But charge separation alone does not guarantee high-efficiency operation. A photogenerated change in energy that is freely available to do useful work is also essential, and is strongly dependent on semiconductor optical properties and reaction kinetics. Notably, this information reveals that specificity of interfacial chemical reactions – even when they are not preceded by charge separation elsewhere – can result in efficient solar energy conversion. We expect that this tutorial will guide researchers in their pursuit to uncover new mechanisms for light to perform useful work.

14 SOLAR ENERGY↗
DOE OSTI2020

Performance Evaluation of AG-1 FC HEPA Filters and Medium in Nuclear Complex Facilities - 20434

High Efficiency Particulate Air (HEPA) filters are credited as the final barrier against the release of radioactive aerosol contamination in nearly every operating U. S. Department of Energy (DOE) and National Nuclear Security Agency (NNSA) nuclear facility. The Institute for Clean Energy Technology (ICET) at Mississippi State University maintains a research program that studies various aspects of these components of containment systems and seeks to answer key questions from across the industry. This programmatic overview will include results from recent and ongoing studies, including a study on filter design and performance envelope, degradation due to aging or fatigue, fire event impact on HEPA filters, and accelerated aging of medium and component parts. ICET has completed a study helping to define the threshold for combined elevated temperature and relative humidity resistance of separator style and separator-less style HEPA filters that are used in DOE complexes. The study provides experimental data that considers elevated temperature, elevated relative humidity (RH), and target differential pressure (dP) in an attempt to gain a better understanding or more comprehensive insight into the operating envelope of different configurations of ASME (American Society for Mechanical Engineers) AG-1 Section FC filters. Test variables including air temperature ranges of 48.88 deg. C, 54.44 deg. C, or 60 deg. C (120, 130, or 140 deg. F); air RH ranges including 60-70%, 80%, or 90+% with initial filter dP of either 497.68 Pa, 746.52 Pa, or 995.36 Pa (2 inches water column (in. w. c.), 3 in. w. c., or 4 in. w. c.) included for full bracketing of each set of conditions. A study examining the effects of aging and fatigue on nuclear grade HEPA filters and medium elucidates the physical properties of media along with testing of new and aged ASME AG-1 HEPA filters. In addition to performing autopsies of the tested, aged filters to help better understand the service life of their performance, newly designed accelerated aging chambers allow new medium to undergo accelerated aging and exposure treatments. Evaluation of degradation in physical properties and functionality requires a high population of aged filters of different ages, designs, manufacturers, and operational histories. The limited availability of this population and the impropriety of looking solely at formulations and components from past years requires experimental design that will provide foresight into future filter performance. This undertaking therefore involves not only properties analysis of aged medium from different operational histories, ages, designs, and manufacturers, but also prescience in analytical goals via accelerated aging studies of newly manufactured media. A more complete understanding of the mechanisms of degradation from past formulations as well as current and future formulations is possible. Accelerated Aging of newly manufactured media utilizes exposure treatments based on the Arrhenius equation to artificially age medium and therefore glimpse into the future. Evaluation of current media properties This testing allows for the comparison of performance and durability of new filters under upset or design basis conditions with aged filters that were in service under ambient conditions and other aged filters retained in storage. Susceptibility of HEPA filters to the effects of combustion byproducts, heat, and water is also being studied in an attempt to gain understanding on the prevention of filter failures during fire events in nuclear facilities. ICET seeks to determine the effect of filter performance due to smoke loading using characterized smoke from variable fuel compositions, geometries, and loads; the smoke capacity of various filters with different burning conditions, heat release rates and transport, as well as mass transport from the air stream will be studied in low flow containment systems. (authors)

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

CyberMESA: Evaluation Procedures and Metrics for Charging Infrastructure Cybersecurity

A report detailing the cybersecurity evaluation and assessment procedures for EV charging infrastructure. This report is authored by three national lab contributors. INL, PNNL, NRL. The grid integration of Electric Vehicle Charging Infrastructure (EVCI) is a large and complex system of systems (SoS) made up of components from many manufacturers and integrated by various companies responsible for their operation and maintenance. This complexity makes it challenging to implement security standards like IEC 62443 across the entire infrastructure. This document aims to aid in the vulnerability assessments of individual devices, such as Electric Vehicle Supply Equipment (EVSE), to support vendors in adopting and implementing security standards such as UL 2900, NIST SP 800-53 and Common Criteria (CC) (ISO/IEC 15408). This document serves as a guide to help evaluate the cybersecurity posture of EVCI assets. It aims to provide a consistent method for evaluating EVCI, enabling cybersecurity research teams to systematically test a wide range of electric vehicle (EV) assets. The target audiences include cybersecurity research teams, vendor development teams, automotive OEMs, and third-party evaluators.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Multi-phase characterization of pitch-carbon coated nano-silicon anodes for lithium-ion batteries

Silicon (Si) is a leading next-generation Li-ion battery anode candidate that meets rigorous performance demands for portable power including enhanced power and energy density with robust cycling performance. However, a series of complex and interrelated reactions lead to reduced calendar life in Si-containing systems and therefore challenge practical adoption. In the present work, we probe the mechanisms underlying observed performance improvements by adding a pitch-carbon coating onto nano-Si material. We pair solid-phase (X-ray photoemission spectroscopy, Fourier-transform infrared), semi-volatile phase (solid-phase microextraction-gas chromatography-mass spectrometry), and gas-phase (gas chromatography-flame-ionization detector) characterization signals to comprehensively evaluate the impact of pitch-carbon coating on the evolution of the Si solid-electrolyte interphase (SEI) and the associated impacts on electrode/electrolyte reactivity. The pitch-carbon is found to serve as a physicochemical barrier, reducing the electro-active surface area for Si/electrolyte reactivity and preventing Si oxidation. Further, the pitch-carbon coating promotes the evolution of a more-favorable SEI by subsuming substantial functionality typically associated with the fluoroethylene carbonate (FEC) electrolyte additive - such as alkoxide scavenging and suppression of transesterification pathways - and by shifting the competitive electrolyte degradation pathways' favorability. The multi-phase characterization approach enables holistic end-products evaluation from complex (electro)chemical interfacial reactions, which informs a robust interpretation of the carbon coating's role in electrochemical performance improvements. The present mechanistic evaluation aids the rational design for improved nano-Si materials.

25 ENERGY STORAGE↗
DOE OSTI2020

CyDER: A Cyber Physical Co-simulation Platform for Distributed Energy Resources in Smartgrids

The CyDER project aimed at developing an open-source, modular and scalable co-simulation platform for power grids with large shares of Distributed Energy Resources (DERs). The project partners are the Lawrence Berkeley National Lab (LBNL), Lawrence Livermore National Lab (LLNL), PG&E, SolarCity, and ChargePoint. The prime recipient is LBNL; SolarCity and ChargePoint were partners for the project’s first two years. Increased DER integration introduces a number of challenges in power grid operation including a more dynamic interaction between the transmission grid and distribution grids, and increased modeling complexity. Although specialized software exists to precisely model different components of the power system, it is far from trivial to integrate all various models and perform a holistic simulation. Instead of replicating all models in a common simulation program, a commonly accepted approach to tackle this model diversity is to couple third-party simulators and models through a co-simulation platform that coordinates information exchange among the various components. Following this line of research, this project’s objective was to develop a co-simulation platform based on a widely accepted industrial standard called Functional Mockup Interface (FMI). Within this process, the project developed models compliant with the FMI standard, called Functional Mockup Units (FMUs), and used them to perform various operational and planning power system analyses. Relying and building upon an industrial standard is the main differentiation of this project compared with previous or parallel efforts in the co-simulation area. Particular emphasis was put on delivering software utilities to facilitate setting up and running co-simulations by end-users. Furthermore, a strong aspect of this project is demonstrating that co-simulation techniques can be used to perform Hardware-in-the-Loop (HIL) simulations that couple software components (e.g., simulated models) with hardware components (e.g., real devices such PV systems and batteries). The long-term goal of CyDER project is to help establish FMI as a powerful standard for co-simulation and promote adoption by electric utilities and other interested stakeholders. The main accomplishments of the project include the development of several FMUs including distribution and transmission grid models, PV inverters with Volt/Var/Watt controllers, batteries, and predictive optimal controllers. Additionally, a unique software package was developed, called SimulatorToFMU, which is capable of exporting any Python-driven simulator or Python script as an FMU. This is an important contribution towards establishing FMI as one of the main co-simulation standards, because more and more third-party programs for sub-system modeling and simulation are delivered with Python APIs. The CyDER platform was used to perform PV hosting capacity analyses in real utility feeders with and without smart inverter controls, battery storage, and EV charging. Smart inverter controls include conventional Volt/Var/Watt controls for reactive power support and active power curtailment, but also predictive controls that optimize the charging and discharging profile of the battery connected on the DC side in order to minimize the customer’s economic benefit. Finally, an important result of this project is delivering an experimental setup that consists of residential-scale PV inverters with battery storage, a real-time grid simulator with an ideal voltage source as grid emulator, and micro Phasor Measurement Units (PMUs). All these components and additional software modules are coupled to one another using the FMI standard and can be co-simulated with the CyDER platform.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

An experimentally validated electro-thermal EV battery pack model incorporating cycle-life aging and cell-to-cell variations

Lithium-ion batteries are used in a wide variety of applications. To meet the power and energy demands of these applications battery packs are composed of hundreds to thousands of cells. The electrical and thermal interactions between cells introduce additional complexity in the pack dynamics. To capture these effects, a battery pack model composed of 192 cells based on a first-generation (2012) Nissan Leaf battery pack is developed in MATLAB/Simulink/Simscape. Here, with this model, we simulate the electrical dynamics (using a first-order equivalent-circuit model), the thermal dynamics (using a first-order lumped-parameter thermal model), and the aging dynamics (using a semi-empirical severity factor-based model) of every cell in the pack and we also create a pack thermal model that explicitly captures the heat exchange between the modules, and the cells contained within, during operation. The models are calibrated and validated, both at the cell and pack level, with experimental data. Two different case studies of this pack model are investigated. In the first case study, an initial, normally-distributed, cell-to-cell capacity variation is introduced and its effect on the pack voltage and module temperatures is studied. In the second case study, we deliberately insert cells with lower than nominal capacity into the pack and we investigate how this type of initial cell-to-cell capacity variation affects the pack’s ability to deliver energy over time. Finally, we also study how parallel-connected cells can reduce the effects of cell-to-cell variations at the expense of increased aging of the pack overall.

25 ENERGY STORAGE↗
DOE OSTI2022

Exchange-Mediated Transport in Battery Electrolytes: Ultrafast or Ultraslow?

Understanding the mechanisms of charge transport in batteries is important for the rational design of new electrolyte formulations. Persistent questions about ion transport mechanisms in battery electrolytes are often framed in terms of vehicular diffusion by persistent ion–solvent complexes versus structural diffusion through the breaking and reformation of ion–solvent contacts, i.e., solvent exchange events. Ultrafast two-dimensional (2D) IR spectroscopy can probe exchange processes directly via the evolution of the cross-peaks on picosecond time scales. However, vibrational energy transfer in the absence of solvent exchange gives rise to the same spectral signatures, hiding the desired processes. Here we employ 2D IR on solvent resonances of a mixture of acetonitrile isotopologues to differentiate chemical exchange and energy-transfer dynamics in a comprehensive series of Li + , Mg 2+ , Zn 2+ , Ca 2+ , and Ba 2+ bis(trifluoromethylsulfonyl)imide electrolytes from the dilute to the superconcentrated regime. No exchange phenomena occur within at least 100 ps, regardless of the ion identity, salt concentration, and presence of water. All of the observed spectral dynamics originate from the intermolecular energy transfer. These results place the lower experimental boundary on the ion–solvent residence times to several hundred picoseconds, much slower than previously suggested. With the help of MD simulations and conductivity measurements on the Li + and Zn 2+ systems, we discuss these results as a continuum of vehicular and structural modalities that vary with concentration and emphasize the importance of collective electrolyte motions to ion transport. These results hold broadly applicable to many battery-relevant ions and solvents.

25 ENERGY STORAGE↗
DOE OSTI2024

Consequences of plane-strain and plane-stress assumptions in fully coupled chemo-mechanical Li-ion battery models

In Li-ion battery research, it is common to simulate chemo-mechanical phenomena in reduced dimensions (e.g., 2-D) as opposed to fully resolve these complex physics in 3-D. It is common to assume either (1) the out-of-plane strain is negligible (commonly referred to as plane-strain), or (2) the out-of-plane stress is negligible (commonly referred to as plane-stress). However, there is typically little consideration as to the quantitative consequences of these approximations. Furthermore, the influence of these out-of-plane assumptions can be compounded and convoluted when chemo-mechanics models implement so-called “fully coupled” formulations, where the local species concentrations influence the stress-state and the stress-state influences the local species fluxes. Here, the present manuscript explores the implications of using plane-stress and plane-strain assumptions in 2-D as compared to simulating a full 3-D electrode particle. This comparative study includes simulating both isotropic and anisotropic particle intercalation where the particles can be surrounded by either a liquid or solid electrolyte. Additionally, common Li-ion battery-model metrics such as the state-of-stress, intercalation fraction distribution, and specific capacity are compared, while also considering the effects of particle size and C-rate. As alternatives to the pure plane-strain and plane stress approximations, two modified plane-strain assumptions are found to better approximate the fully coupled chemo-mechanical 3-D behavior.

25 ENERGY STORAGE↗
DOE OSTI2026

4D insights into lithium-ion battery sidewall rupture during thermal runaway

Thermal runaway (TR), characterized by rapid exothermic reactions, presents a serious safety risk in lithium-ion batteries (LiBs). External triggers such as high temperatures, mechanical abuse, or internal short circuits (ISCs) can initiate TR, often resulting in sidewall rupture, which may escalate to catastrophic battery pack failure. In this study, we developed and applied high-speed synchrotron imaging techniques to investigate sidewall rupture mechanisms in LiBs subjected to different triggering scenarios. Using in situ 4D tomographic imaging, we visualized the dynamic evolution of sidewall rupture with high spatial and temporal resolution. The results revealed distinct failure behaviors linked to each trigger, underscoring the complex and condition-specific nature of sidewall breach and battery failure. These insights highlight the critical importance of implementing tailored safety strategies across diverse applications. Our findings demonstrate the powerful potential of synchrotron high-speed tomography as a diagnostic tool for advanced safety testing and cell qualification.

25 ENERGY STORAGE↗
DOE OSTI2026

Hidden features in the OH-stretching spectra of amino acid decorated air–water interfaces

Chemical reactivity at the air–water interface is governed by the interfacial solvation of reactive species. For instance, during aqueous amino acid-based CO 2 absorption, water reorganizes around the reactive sites and couples dynamically with reaction pathways, facilitating the reaction. In this context, surface-sensitive vibrational sum-frequency generation (vSFG) spectroscopy can probe the OH stretch vibrations of interfacial water and determine the solvation structures around reactants and products, thereby furthering our understanding of the role of interfacial solvation. However, vSFG spectra of the air–water interface in the presence of charged species can be remarkably complex; key species-bound local water structures with distinct orientations may be hidden beneath prominent vSFG peaks arising from water–water hydrogen bonds and remain difficult to resolve. Here, we measure and compute vSFG spectra of the water OH stretch at air–water interfaces decorated with amino acids in their zwitterionic and anionic forms, as well as equimolar mixtures of these forms with bicarbonate. The latter represents post-CO 2 -absorption conditions. We find that computing depth- and frequency-dependent spectral densities—decomposed into contributions from water molecules hydrogen-bonded exclusively to other water molecules, exclusively to amines, exclusively to carboxylates, or shared between these polar/charged groups—is indispensable for accurate interpretation of the vSFG spectra. Key findings include orientational flip-flop in water sub-layers, strong carboxylate-water H-bonding, and water orientational ordering extending into the bulk aqueous phase induced by anionic amino acids. Here, this study provides a computational spectroscopic platform for improved understanding of interfacial solvation relevant to interfacial reactivity.

Air-water interface↗
DOE OSTI2022

Optimization of Water-Alternating-CO2 Injection Field Operations Using a Machine-Learning-Assisted Workflow

Summary This paper will present a robust workflow to address multiobjective optimization (MOO) of carbon dioxide (CO2)-enhanced oil recovery (EOR)-sequestration projects with a large number of operational control parameters. Farnsworth unit (FWU) field, a mature oil reservoir undergoing CO2 alternating water injection (CO2-WAG) EOR, will be used as a field case to validate the proposed optimization protocol. The expected outcome of this work would be a repository of Pareto-optimal solutions of multiple objective functions, including oil recovery, carbon storage volume, and project economics. FWU’s numerical model is used to demonstrate the proposed optimization workflow. Because using MOO requires computationally intensive procedures, machine-learning-based proxies are introduced to substitute for the high-fidelity model, thus reducing the total computation overhead. The vector machine regression combined with the Gaussian kernel (Gaussian-SVR) is used to construct proxies. An iterative self-adjusting process prepares the training knowledge base to develop robust proxies and minimizes computational time. The proxies’ hyperparameters will be optimally designed using Bayesian optimization to achieve better generalization performance. Trained proxies will be coupled with multiobjective particle swarm Optimization (MOPSO) protocol to construct the Pareto-front solution repository. The outcomes of this workflow will be a repository containing Pareto-optimal solutions of multiple objectives considered in the CO2-WAG project. The proposed optimization workflow will be compared with another established methodology using a multilayer neural network (MLNN) to validate its feasibility in handling MOO with a large number of parameters to control. Optimization parameters used include operational variables that might be used to control the CO2-WAG process, such as the duration of the water/gas injection period, producer bottomhole pressure (BHP) control, and water injection rate of each well included in the numerical model. It is proved that the workflow coupling Gaussian-SVR proxies and the iterative self-adjusting protocol is more computationally efficient. The MOO process is made more rapid by squeezing the size of the required training knowledge base while maintaining the high accuracy of the optimized results. The outcomes of the optimization study show promising results in successfully establishing the solution repository considering multiple objective functions. Results are also verified by validating the Pareto fronts with simulation results using obtained optimized control parameters. The outcome from this work could provide field operators an opportunity to design a CO2-WAG project using as many inputs as possible from the reservoir models. The proposed work introduces a novel concept that couples Gaussian-SVR proxies with a self-adjusting protocol to increase the computational efficiency of the proposed workflow and to guarantee the high accuracy of the obtained optimized results. More importantly, the workflow can optimize a large number of control parameters used in a complex CO2-WAG process, which greatly extends its utility in solving large-scale MOO problems in various projects with similar desired outcomes.

Energy & Fuels↗
DOE OSTI2021

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗
DOE OSTI2020

Surface SiO 2 Thickness Controls Uniform-to-Localized Transition in Lithiation of Silicon Anodes for Lithium-Ion Batteries

Silicon is a promising anode material for lithium-ion batteries because of its high capacity, but its widespread adoption has been hampered by a low cycle life arising from mechanical failure and the absence of a stable solid–electrolyte interphase (SEI). Understanding SEI formation and its impact on cycle life is made more complex by the oxidation of silicon materials in air or during synthesis, which leads to SiO x coatings of varying thicknesses that form the true surface of the electrode. Here, the lithiation of SiO 2 -coated Si is studied in a controlled manner using SiO 2 coatings of different thicknesses grown on Si wafers via thermal oxidation. SiO 2 thickness has a profound effect on lithiation: below 2 nm, SEI formation followed by uniform lithiation occurs at positive voltages versus Li/Li + . Si lithiation is reversible, and SiO 2 lithiation is largely irreversible. Above 2 nm SiO 2 , voltammetric currents decrease exponentially with SiO 2 thickness. For 2–3 nm SiO 2 , SEI formation above 0.1 V is suppressed, but a hold at low or negative voltages can initiate charge transfer whereupon SEI formation and uniform lithiation occur. Cycling of Si anodes with an SiO 2 coating thinner than 3 nm occurs at high Coulombic efficiency (CE). If an SiO 2 coating is thicker than 3–4 nm, the behavior is totally different: lithiation at positive voltages is strongly inhibited, and lithiation occurs at poor CE and is highly localized at pinholes which grow over time. As they grow, lithiation becomes more facile and the CE increases. Pinhole growth is proposed to occur via rapid transport of Li along the SiO 2 /Si interface radially outward from an existing pinhole, followed by the lithiation of SiO 2 from the interface outward.

25 ENERGY STORAGE↗
DOE OSTI2022

A compendium and meta-analysis of flatband potentials for TiO2, ZnO, and SnO2 semiconductors in aqueous media

Semiconductor/electrolyte interfaces are of great interest to numerous scientific fields including renewable energy, (photo)electrochemistry, and energy storage. The semiconductor flatband potential is a key parameter in locating the conduction band minimum or valence band maximum of the semiconductor material in electrolyte. Despite its importance for quantifying the energetic location of the semiconductor bands, literature reports for the same material demonstrate significant variability in the flatband potential. In this compendium and meta-analysis, reported flatband potentials of the common semiconductor materials TiO2, SnO2, and ZnO in aqueous electrolyte were compiled and assessed to quantify the spread in literature flatband potentials as well as determine the factors that lead to the significant spread. For TiO2, SnO2, and ZnO, literature flatband potentials referenced to the reversible hydrogen electrode span a range of nearly 2 V each. Flatband potential tabulations were separated by variables such as the solution pH, the crystalline polymorph, the crystal facet, the morphology, and the dimensions or combinations of these variables to assess the factors that contribute to the observed spread. Important and surprising findings from these categorizations are summarized: (1) Even for the narrowest categorizations, the spread in flatband potential is still large. (2) Flatband potentials of TiO2 and SnO2 follow the expected Nernstian dependence with solution pH. ZnO materials deviate from this Nernstian dependence. (3) In the aggregate, there is no statistically significant difference in the reported flatband potentials of anatase and rutile TiO2. Single crystal tabulations were the only distributions to have statistically significant differences in the flatband potential between anatase and rutile TiO2. (4) Anatase TiO2 materials with a nanotube morphology appear to have a +400 mV difference in mean flatband potential compared to all other morphologies, but we argue that this is likely due to widespread misuse of the Mott–Schottky analysis. Other interesting findings are revealed within the spread of literature flatband potentials, and possible explanations are provided to generate discussion. We also briefly review and discuss common techniques that were used to determine the flatband potential and the pitfalls/criticisms of these techniques. Last, we discuss some ways in which future research on the determination of the flatband potential can be performed to improve the reliability of reported values and the quality of the work. In total, the results from this meta-analysis suggest multiple factors can affect the measured flatband potential and that an abundance of caution should be applied when attempting to quantify the flatband potential of complex or nanostructured systems.

Patel, Milan Y. (ORCID:0000000291948141)↗