Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “utility scale”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dual Mode Intermediate Temperature Fuel Cell: Liquid Fuels and Electricity - CRADA 351 (Abstract)

FuelCell Energy, Inc. (FCE) in collaboration with Pacific Northwest National Laboratory (PNNL), The University of Connecticut (UCONN), The Energy and Environmental Research Center at The University of North Dakota (EERC-UND), and Massachusetts Institute of Technology (MIT), has entered into a co-operative agreement with the U.S. Department of Energy Advanced Research Project Agency – Energy (ARPA-E), agreement number DE-AR0000503, to develop an electrochemical device for conversion of methane gas to a usable liquid product, such as methanol or formaldehyde. The objective of this CRADA is to develop the anode electro-catalyst for the aforementioned electrochemical cell during Q1-Q8. PNNL will utilize thermodynamic modeling, batch-mode reactor screening and button-scale electrochemical testing to develop this anode."

09 BIOMASS FUELS↗
DOE OSTI2025

CABLE Stage 3 – Ultra‐High Strength/Highly Conductive Al Alloys

This work is to be conducted in support of the CABLE Conductor Manufacturing Prize Program awarded by NREL. The intent is to connect competitor teams with national laboratories that can help accelerate the development of innovative solutions and products. Teams who receive CABLE Prize awards are eligible to utilize vouchers at national laboratories to advance their ideas. NanoAL LLC has been awarded a voucher to utilize at Ames National Laboratory. Ames National Laboratory will perform a scale‐up trial of Nano 6000‐T9 to evaluate and optimize the heat treatment and other processing variables to validate the lab‐scale work done at NanoAL. The final processing schedule will be used to fabricate the submission samples for the Stage 3 contest. All work will be supported by mechanical and electrical conductivity testing be performed by Ames.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Onshore U.S. Carbon Pipeline Deployment: Siting, Safety, and Regulation

Carbon capture, utilization, and storage technology has significant potential to reduce greenhouse gas emissions and mitigate the impact of climate change, particularly in hard-to-decarbonize industrial and commercial sectors. Unless utilization or storage occurs at the same location as carbon capture, which is rare, carbon must be transported from a point source to a utilization or storage site. Pipelines offer significant advantages for large-scale transportation of carbon dioxide over other methods. While studies show that reaching net-zero carbon emissions in the United States by 2050 will require between 29,000 and 66,000 miles of carbon pipelines, the U.S. had deployed fewer than 6,000 miles of carbon pipelines by 2022. Although closing this gap is important to achieving low-carbon goals, carbon pipelines operate in a complex and uncertain local, state and federal regulatory landscape and face public concerns about safety and siting. This report covers numerous regulatory issues surrounding carbon pipeline development, including the current narrow federal definition of carbon dioxide and the considerable variation in state and local governments’ laws and regulations. The report serves as a primer for regulators and stakeholders who seek to better understand the regulatory challenges and opportunities facing this critical infrastructure

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Bridging the Gap between Cosmological Simulations with Graph Neural Networks and Domain Adaptation

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross dataset tasks (up to 28% better relative error and up to almost an order of magnitude better χ 2 ). Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Domain Adaptive Graph Neural Networks for Constraining Cosmological Parameters Across Multiple Data Sets

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross-dataset tasks (up to $28\%$ better relative error and up to almost an order of magnitude better $\chi^2$). Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2025

A visco-plastic constitutive model for accurate densification and shape predictions in powder metallurgy hot isostatic pressing

Powder metallurgy hot isostatic pressing (PM-HIP) is an advanced manufacturing process that produces near net shape parts with high material utilization and uniform microstructures. Despite being used frequently to produce small-scale components, the application of PM-HIP to large-scale components is limited due to inadequate understanding of its complex mechanisms that cause unpredictable post-HIP shape distortions. A computational model can provide necessary information about the intermediate and final stages of the HIP process that can help understand it better and make accurate predictions. Generally, two types of computational models are employed for PM-HIP of metal powders, namely, plastic and visco-plastic models. Between these, the plastic model is preferred due to its cheaper calibration approach requiring less experimental data. However, the plastic model sometimes produces incorrect predictions when slight variations of the HIP conditions are encountered in practical situations. Therefore, this work presents a visco-plastic model that addresses these limitations of the plastic model. A novel modified calibration approach is employed for the visco-plastic model that utilizes less experimental data than existing approaches. With the new approach, the data requirement is same for both plastic and visco-plastic models. This also enables a quantitative comparison of plastic and visco-plastic models, which have been only qualitatively compared in the past. When calibrated with the same experimental data, both the models are found to produce similar results. In conclusion, the calibrated visco-plastic model is applied to several complex geometries, and the predictions are found to be in good agreement with experimental observations.

Hot isostatic pressing↗
DOE OSTI2021

Improved forward voltage and external quantum efficiency scaling in multi-active region III-nitride LEDs

Ultra-low voltage drop tunnel junctions (TJs) were utilized to enable multi-active region blue light emitting diodes (LEDs) with up to three active regions in a single device. The multi-active region blue LEDs were grown monolithically by metal-organic chemical vapor deposition (MOCVD) without growth interruption. This is the first demonstration of a MOCVD grown triple-junction LED. Optimized TJ design enabled near-ideal voltage and EQE scaling close to the number of junctions. Furthermore, this work demonstrates that with proper TJ design, improvements in wall-plug efficiency at high output power operation are possible by cascading multiple III-nitride based LEDs.

42 ENGINEERING↗
DOE OSTI2022

Field Study of Grid-connected Heat Pump Water Heaters in the Southeast U.S.: The Next Right Thing

The integration of grid-connected functionality and advanced control algorithms into heat pump water heaters (HPWH) offers the capability to shift load with minimal customer impact. This capability provides a flexible grid resource to utilities, while the increased energy efficiency of HPWHs offers customers a lower electric bill. The combined utility/customer value from grid-connected HPWHs is compelling in the Southeast U.S. where residential electric water heating is prevalent and utility load management is common. This paper presents the results of a HPWH load shifting study conducted in Central Florida using the CTA-2045 standard. Building upon previous research in the Pacific Northwest, this study consisted of approximately 45 occupied homes equipped with HPWHs undergoing load shifting strategies weekly for over a year. Curtailment durations ranged from three to five hours in the morning, and four to five hours in the evening to coincide with high-value periods for utility coincident load for system-wide electric demand reductions. During the morning and afternoon, a one- or two-hour load-up event preceded curtailment. Baseline data were collected across varied Florida weather in which no load shifting events were implemented. Results from traditional load shifting strategies were analyzed across seasons and used to devise unique load shifting approaches to increase renewable energy use during periods of high solar energy generation. Regional impacts are forecasted for large-scale implementation of strategies. Lessons learned and recommendations are also provided for how utilities, manufacturers, and regional planners can maximize load shifting benefits from grid-connected HPWHs.

Heat pump water heaters, residential building, loa↗
DOE OSTI2023

Hierarchical off-diagonal low-rank approximation of Hessians in inverse problems, with application to ice sheet model initialization

Obtaining lightweight and accurate approximations of discretized objective functional Hessians in inverse problems governed by partial differential equations (PDEs) is essential to make both deterministic and Bayesian statistical large-scale inverse problems computationally tractable. The cubic computational complexity of dense linear algebraic tasks, such as Cholesky factorization, that provide a means to sample Gaussian distributions and determine solutions of Newton linear systems is a computational bottleneck at large-scale. These tasks can be reduced to log-linear complexity by utilizing hierarchical off-diagonal low-rank (HODLR) matrix approximations. In this work, we show that a class of Hessians that arise from inverse problems governed by PDEs are well approximated by the HODLR matrix format. In particular, we study inverse problems governed by PDEs that model the instantaneous viscous flow of ice sheets. In these problems, we seek a spatially distributed basal sliding parameter field such that the flow predicted by the ice sheet model is consistent with ice sheet surface velocity observations. Here, we demonstrate the use of HODLR Hessian approximation to efficiently sample the Laplace approximation of the posterior distribution with covariance further approximated by HODLR matrix compression. Computational studies are performed which illustrate ice sheet problem regimes for which the Gauss–Newton data-misfit Hessian is more efficiently approximated by the HODLR matrix format than the low-rank (LR) format. We then demonstrate that HODLR approximations can be favorable, when compared to global LR approximations, for large-scale problems by studying the data-misfit Hessian associated with inverse problems governed by the first-order Stokes flow model on the Humboldt glacier and Greenland ice sheet.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Noble-Metal-Free, Nickel-Based Dual Functional Materials for Improved Methane Production from In Situ Carbon Dioxide Capture and Conversion

Promoters for dual functional materials have not been well explored, but promoters could improve the efficiency of the process by improving the selectivity of the CO 2 methanation process. Utilizing integrated capture and conversion, where CO 2 is captured and converted to useful products, would allow for a useful avenue to control CO 2 emissions. One such way to accomplish this would be to utilize materials that can both capture and convert CO 2 to useful products. However, these materials are often based on costly noble metals, like ruthenium and platinum, decreasing their viability on an industrial scale. Less expensive metals, for example, nickel, would allow for dual functional materials to be more readily utilized in industrial settings. Nickel-based dual functional materials often do not react with the captured CO 2 and merely desorb the CO 2 rather than form a useful product. However, promoters have not been well explored for these types of materials to improve the catalytic properties, which would be beneficial to improve nickel-based materials. Herein, we report the addition of ytterbium on a nickel-based dual functional material and the improvements to the production of methane from captured CO 2 with the incorporated ytterbium promoter. The ytterbium promoter improves the selectivity of the catalysts for the hydrogenation of captured CO 2 to methane and increases the ability for the material to capture CO 2 due to additional basic sites being formed on the surface of alumina. As a result, the 12%Ni/4%Yb/6%Na 2 O/Al 2 O 3 catalyst was utilized to capture carbon dioxide and then convert the captured CO 2 to methane over five cycles, where both the amount captured and the amount converted remained stable, indicating the stability of the material over long-term use.

Catalysts↗
DOE OSTI2022

Elevated temperature nanoscratch of Inconel 617 Superalloy

Inconel 617 superalloy is a main candidate to be used for mechanical and tribo-components in high temperature helium-cooled reactors. Recent findings show that it grows a unique surface oxide, especially under high temperature helium with distinct wear, friction, and contact properties. Here, this study reports the elevated temperature nanoscratch behavior of Inconel 617 and further utilizes it to understand the effect of temperature on contact friction constituent contributors, adhesion and plowing at small scales. Inconel 617 is aged in high temperature helium, and consequently, the total kinetic friction coefficient of the alloy surface oxide is obtained in temperatures ranging from 25 °C to 400 °C. A finite element model is developed and validated based on the experimental results. The model is then utilized along with previously established techniques to determine the adhesion and plowing components of the friction coefficient. At small scale, the experimental results show that with increasing temperature the friction coefficient increases. It was inferred that this increase is mainly due to the increased contribution of plowing friction at high levels of deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Generalized Framework for Service Restoration in a Resilient Power Distribution System

An electric power grid is one of the complex infrastructures, and because of its complex nature, there is simply no way that outages can be completely avoided. Thus, a modern society that depends on reliable electric supply requires a resilient electric system that can recover from disruptions while integrating emerging smart grid technologies. Here, this article presents a novel approach for service restoration in a modern power distribution system with controllable switches and distributed generation (DG) resources for any kind of outage. The proposed framework supports both the traditional service restoration using feeder reconfiguration and the grid-forming DG-assisted intentional islands that are dynamically sized using algorithms based on the fault scenario, available resources, and priority of loads. The problem is formulated as a mixed-integer linear program that incorporates critical system connectivity and operating constraints. Simulations are performed to demonstrate the effectiveness of the proposed approach using a large-scale four-feeder 1069-bus three-phase unbalanced distribution test system. It is demonstrated that the framework is effective in utilizing all available resources in quickly restoring the power supply to improve resiliency during extreme events and is scalable for a large-scale unbalanced power distribution system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Analyzing the Economic and Technical Potential of H2@Scale

H2@Scale is one of EERE's flagship initiatives to enable clean and affordable hydrogen production, storage, distribution, and utilization across multiple sectors in the economy. The January H2IQ hour will discuss results from a recent multi-laboratory analysis report quantifying the technical and economic potential of H2@Scale, given the use of hydrogen as an energy intermediate to integrate sectors in the U.S. energy system. Mark Ruth of the National Renewable Energy Laboratory, the lead author of the report, will discuss specific examples of the role hydrogen can play in industry and transportation sectors, such as steelmaking, synthetic fuels, and blending into natural gas pipelines. He will also explain how hydrogen production can be integrated with the electric grid to monetize low-cost energy and incentivize growth in renewable energy supply. The H2IQ Hour will dive into estimates for future hydrogen consumption in current and emerging sectors, given R&D advances and varying prices of natural gas and electricity.

39 EE - Hydrogen and Fuel Cell Technologies (EE-3F↗
DOE OSTI2022

The Imperial Valley Dark Fiber Project: Toward Seismic Studies Using DAS and Telecom Infrastructure for Geothermal Applications

We report that the Imperial Valley is a seismically active basin occupying the southern end of the Salton trough, an area of rapid extension, high heat flow, and abundant geothermal resources. This report describes an ongoing large-scale distributed acoustic sensing (DAS) recording study acquiring high-density seismic data on an array between Calipatria and Imperial, California. This 27 km array, operating on dark fiber since 9 November 2020, has recorded a wealth of local seismic events as well as ambient noise. The goal of the broader Imperial Valley Dark Fiber project is to evaluate passive DAS as a tool for geothermal exploration and monitoring. This report is intended to provide installation information, noise characteristics, and metadata for future studies utilizing the data set. Because of the relatively small number of basin-scale DAS studies that have been conducted to date, we also provide a range of lessons learned during the deployment to assist future researchers exploring this acquisition strategy.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Bridging confined phase behavior of CH 4 -CO 2 binary systems across scales

Phase behavior of confined fluids may deviate significantly from that of the bulk fluid due to the fluid-wall interactions being a significant portion of all intermolecular interactions under confinement. Despite recent advancements in understanding confined phase behavior of pure fluids, confined phase behavior of mixtures remains an understudied topic. In this work, we examine the confined phase behavior of a CH 4 -CO 2 binary system by combining Monte Carlo (MC) simulations, a cubic equation of state (EoS), and the lattice Boltzmann method (LBM). First, the effects of confinement on density and phase distribution in nano-size pores are established using Gibbs Ensemble MC calculations, which produce precise results of liquid and vapor confined pressures and account for the modification of the phase change location. By comparing the phase envelopes of bulk and confined mixtures at a fixed temperature, here it is observed that the phase envelopes shrink with reductions in pore size. Based on this observation, we extend a modified Peng-Robinson EoS, which was originally developed for pure fluids under confinement, to mixtures via van-der-Waals-type mixing rules and by accounting for shifts in the critical properties of confined CH 4 -CO 2 . The resulting phase envelopes are in good agreement with the MC data. In addition, a local density model is used in combination with the confined EoS to calculate adsorption isotherms of CH 4 -CO 2 mixtures and to characterize the behavior of confined matter in nanopores. Finally, we incorporate this EoS in a multicomponent multiphase LBM that uses a pseudopotential model to represent intermolecular forces. This workflow utilizes multiscale simulation techniques to bridge the behavior of multicomponent systems across scales and to shed light on the confined phase behavior of CH 4 -CO 2 binary systems.

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

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

When and Where Lithium Plating Occurs, Its Correlation with Microstructure Heterogeneity, and the Mechanisms That Initiate and Self-Regulate Electrochemical Heterogeneity (A02-0444)

A microstructure scale electrochemical LIB model was used to investigate lithium plating onset, material non-uniform utilization, and in-plane heterogeneities for an NMC-graphite full cell. Model predicts active material particle surface roughness and size distribution (respectively, non-uniform curvature within and between particles) initiate in-plane heterogeneity, and that particle size heterogeneity at the separator interface controls the lithium plating preferential deposition ("Where"). These in-plane heterogeneities are then exacerbated by through-plane heterogeneities induced at fast charge as electrolyte depletion occurs and concentrates intercalation reaction near the anode-separator interface. Also, magnitude and occurrence of lithium plating is controlled by effective, or macroscale, microstructure parameters ("When"). As local states of charge start to diverge between nearby active material regions, overpotential differences induced by OCP difference kick in and contribute to reduce these SOC local heterogeneities. However, for staged materials such as graphite, with OCP profile alternating between plateaus and varying regions, this balancing mechanism is, respectively, inactive and active. This leads to a dynamic, non-monotonic, in-plane heterogeneity time evolution for state of charge and Faraday current density, for which their respective in-plane heterogeneity magnitude alternates. Such behavior has been modeled both for the whole electrode at the microstructure scale and at the particle scale. In-plane heterogeneities are usually considered to be detrimental, as they result in material non-uniform utilization (i.e., under and over stressed regions) and earlier degradations. However, this work provides a more granular approach as it discriminates between a harmful in-plane heterogeneity (non-uniform curvature) that triggers SOC in-plane heterogeneity, and a beneficial in-plane heterogeneity (Faraday current density) that contributes to reduce SOC in-plane heterogeneity. This work comprehensively explains the mechanisms that initiate, exacerbate, and regulate heterogeneity at the microstructure scale, while providing some design suggestions to reduce both in-plane and through-plane heterogeneities, as summarized in the graphical abstract.

ADVANCED PROPULSION SYSTEMS↗