Mechanistic electrochemical modeling of simulated spent nuclear fuel (SIMFUEL) degradation
Poster summarizing modeling of UO2 (SIMFUEL) corrosion using the Zuzax code suite
Jove-Colon, Carlos F. [Sandia National Laboratorie↗
SEARCH · Engineering Papers
Results for “electrochemical modeling and simulation”
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.
Time-Resolved X-ray Operando Observations of Lithiation Gradients across the Cathode Matrix and Individual Oxide Particles during Fast Cycling of a Li-Ion Cell
Lithiated transition metal oxides serve as active materials in the positive electrode (cathode) of lithium-ion cells. During electrochemical cycling, lithium ions intercalate and deintercalate into these oxide particles. This behavior causes two types of lithiation gradients to emerge: (i) a bulk gradient across the depth of the cathode matrix (averaged over individual oxide particles) and (ii) a microscopic gradient across the particles themselves, which also depends on their location in the electrode. Here we show how both gradients can be studied using operando X-ray diffraction during 4C charge and 4C discharge. The oxide (de)lithiation is estimated from the unit cell parameters by indexing the X-ray diffraction spectra. By fitting the lithiation profiles with orthogonal polynomials, the bulk gradients across the electrode thickness are quantified. These gradients develop as the current flows through the cell and dissipate during open-circuit and potentiostatic-hold periods. Further details of lithiation dynamics can be obtained through shape analysis of the Bragg peaks. In particular, from electrochemical model simulations, we show that the width and skewness of the (003) peak track (de)lithiation fronts moving across the individual oxide particles.
(003) peak↗
Opportunities and challenges in process modeling and simulation of electrochemical systems
Electrochemical technologies have garnered intense interest in both academic research and industrial applications due to their potential to increase energy efficiency and reduce carbon footprint. However, electrochemical process fundamentals have been absent from the chemical process simulators available to millions of chemical engineers worldwide. To expedite process research and development of electrochemical technologies in the chemical industry, it is imperative that essential electrochemical process fundamentals be incorporated into process simulators to support modeling and simulation of electrochemical processes. Here, this study examines three process fundamentals key to the research and development of electrochemical processes: electrochemical reaction kinetics, electrolyte thermodynamics, and heat and mass transfer. It further illustrates application of these process fundamentals with a case study modeling an electrochemical process for the conversion of acrylonitrile to adiponitrile. The modeling example highlights the roles of applied voltage on reaction rates, and electron flow rates on the performance of electrochemical processes. It suggests that the applied voltage and electron flow rates are two unique concepts that should be included in simulators for electrochemical processes.
09 BIOMASS FUELS↗
Modeling Isotope Separation in Electrochemical Lithium Deposition
Naturally occurring Li consists of two stable isotopes, 6 Li with an abundance of about 7.5%, and 7 Li making up the remainder with 92.5%. The development of a 6 Li enrichment technique, in terms of technical reliability and environmental safety to reach 6 Li future requirements, represents a key step in the roadmap for nuclear fusion energy supply worldwide. This paper uses finite element analysis-based models to simulate electrochemical Li isotope separation, which is an attractive method in terms of simplicity, safety, and scalability. In the model, we quantitatively analyze how different electrochemical factors including thermodynamics, charge-transfer kinetics, and diffusivities affect the separation process (separation factor), together with cell parameters, such as cell length and current density. The maximum separation factor of 1.128 could be obtained with the cell under the optimal thermodynamic, kinetic, and diffusive conditions, which is among the highest separation factors ever reported. Furthermore, these results will assist in designing the actual isotope separation setup with large separation factor and appropriate timing for sample collection.
07 ISOTOPE AND RADIATION SOURCES↗
High-Throughput Laser Processing and Acoustic Diagnostics for Enhanced Battery Performance and Manufacturing: Cooperative Research and Development Final Report, CRADA Number CRD-20-17191
The National Laboratory of the Rockies (NLR) will conduct electrochemical testing and characterization on the partners lithium ion cells as well as provide and validate laser-ablated electrode materials for enhancing cell rate performance.
36 MATERIALS SCIENCE↗
Deriving effective electrode–ion interactions from free-energy profiles at electrochemical interfaces
Understanding ion adsorption at electrified metal–electrolyte interfaces is essential for accurate modeling of electrochemical systems. Here, in this study, we systematically investigate the free energy profiles of Na + , Cl − , and F − ions at the Au(111)–water interface using enhanced sampling molecular dynamics with both classical force fields and machine-learned interatomic potentials (MLIPs). Our classical metadynamics results reveal a strong dependence of predicted ion adsorption on the Lennard-Jones parameters, highlighting that—without due care—standard mixing rules can lead to qualitatively incorrect descriptions of ion–metal interactions. We present a systematic methodology for tuning the cross term LJ parameters to control adsorption energetics in agreement with more accurate models. As a surrogate for an ab initio model, we employed the recently released Universal Models for Atoms MLIP, which validates classical trends and displays strong specific adsorption for chloride, weak adsorption for fluoride, and no specific adsorption for sodium, in agreement with experimental and theoretical expectations. By integrating molecular-level adsorption free energies into continuum models of the electric double layer, we show that specific ion adsorption substantially alters the interfacial ion population, the potential of zero charge, and the differential capacitance of the system. Our results underscore the critical importance of force field parameterization and advanced interatomic potentials for the predictive modeling of ion-specific effects at electrified interfaces and provide a robust framework for bridging molecular simulations and continuum electrochemical models.
Roncoroni, Fabrice [Lawrence Berkeley National Lab↗
A Diffuse-Interface Model for Predicting the Evolution of Metallic Negative Electrodes and Interfacial Voids in Solid-State Batteries with Homogeneous and Polycrystalline Solid Electrolyte Separators
Here, this paper presents a novel diffuse-interface electrochemical model that simultaneously simulates the evolution of the metallic negative electrode and interfacial voids during the stripping and plating processes in solid-state batteries. The utility and validity of this model are demonstrated for the first time on a cell with a sodium (Na) negative electrode and a Na-β″-alumina ceramic solid electrolyte (SE) separator. Three examples are simulated. First, stripping and plating with a perfect electrode/electrolyte interface; second, stripping and plating with a single interfacial void at the electrode/electrolyte interface; third, stripping with multiple interfacial voids. Both homogeneous and polycrystalline SEs with low and high-conductivity grain boundaries (GBs) are considered for all three examples. Heterogeneous GB conductivity marginally impacts the thickness of the Na electrode in cases with a perfect electrode/electrolyte interface. Moreover, it results in local changes to void growth due to the interactions between the void edge and the GBs. The void growth rate is a linear function of the flux of Na atoms at the void edge, which in turn depends on the applied current density. We also show that the void coalescence rate increases with applied current density and can be marginally influenced by GB conductivity.
batteries↗
Electrode Mesoscale as a Collection of Particles: Coupled Electrochemical and Mechanical Analysis of NMC Cathodes
Battery electrodes are composed of polydisperse particles and a porous, composite binder domain. These materials are arranged into a complex mesostructure whose morphology impacts both electrochemical performance and mechanical response. We present image-based, particle-resolved, mesoscale finite element model simulations of coupled electrochemical-mechanical performance on a representative NMC electrode domain. Beyond predicting macroscale quantities such as half-cell voltage and evolving electrical conductivity, studying behaviors on a per-particle and per-surface basis enables performance and material design insights previously unachievable. Voltage losses are primarily attributable to a complex interplay between interfacial charge transfer kinetics, lithium diffusion, and, locally, electrical conductivity. Mesoscale heterogeneities arise from particle polydispersity and lead to material underutilization at high current densities. Particle-particle contacts, however, reduce heterogeneities by enabling lithium diffusion between connected particle groups. While the porous composite binder domain (CBD) may have slower ionic transport and less available area for electrochemical reactions, its high electrical conductivity makes it the preferred reaction site late in electrode discharge. Mesoscale results are favorably compared to both experimental data and macrohomogeneous models. This work enables improvements in materials design by providing a tool for optimization of particle sizes, CBD morphology, and manufacturing conditions.
42 ENGINEERING↗
Mechanical, Electrochemical & Thermal Modeling of Electric Vehicle Batteries for Crash Simulation (CRADA CRD-19-00811 Final Report)
Under the proposed effort in partnership with Hyundai Motor Company (HMC), the National Laboratory of the Rockies (NLR) will cooperate with HMC to develop mathematical models for battery cells and modules for simulating abuse response in batteries subject to the type of mechanical crushing that can occur in a full motor vehicle crash.
33 ADVANCED PROPULSION SYSTEMS↗
Modeling the Nucleation and Growth of Lead Sulfate Particles on Lead Electrodes
Lead-acid batteries (LABs) play a pivotal role in the energy storage sector with applications spanning from starting-lighting-ignition batteries to grid energy storage. Passivation of lead negative electrodes by PbSO 4 particles is a fundamental mechanism limiting the performance of LAB. In this regard, an electrochemical model is developed that simulates the nucleation and growth (N&G) dynamics of PbSO 4 particles on a flat lead electrode, responsible for its passivation. The model considers the electrochemical reactions between lead electrode and sulfuric acid, N&G of PbSO 4 particles, passivation of the lead surface, and the ternary transport of PbSO 4 (aq), bisulfate, and protons in H 2 SO 4 electrolyte. The model is validated with a dataset of cyclic voltammetry (CV) responses collected at several scan rates and H 2 SO 4 concentrations. The model shows remarkable qualitative and quantitative agreement with the experimental data including CV peak features, discharge capacity, and particle size. The model was employed to explore key N&G quantities, such as supersaturation, nucleation rate, particle count, growth rate, particle size, and surface coverage, and to examine how their interactions influence electrode utilization. Parametric studies were also conducted to evaluate how scan rates and acid concentrations influence the previously mentioned N&G quantities and, subsequently, the utilization of the electrode.
25 ENERGY STORAGE↗
Chapter 6: Abuse Response of Batteries Subjected to Mechanical Impact
Electrochemical and thermal models to simulate nominal performance and abuse response of lithium-ion cells and batteries have been reported widely in the literature. Studies on mechanical failure of cell components and how such events interact with the electrochemical and thermal response are relatively less common. This chapter outlines a framework developed under the Computer Aided Engineering for Batteries program to couple failure modes resulting from external mechanical loading to the onset and propagation of electrochemical and thermal events that follow. Starting with a scalable approach to implement failure criteria based on thermal, mechanical, and electrochemical thresholds, we highlight the practical importance of these models using case studies at the cell and module level. The chapter also highlights a few gaps in our understanding of the comprehensive response of batteries subjected to mechanical crash events, the stochastic nature of some of these failure events, and our approach to build safety maps that help improve robustness of battery design by capturing the sensitivity of some key design parameters to heat generation rates under different mitigation strategies.
abuse simulations↗
A comprehensive modeling method for proton exchange membrane electrolyzer development
Hydrogen attracts significant interests as an effective energy carrier that can be derived from renewable sources. Hydrogen production using a proton-exchange membrane (PEM) electrolyzer can efficiently convert renewable power via water splitting in wide scales—from large, centralized generation to on-site production. Mathematical models with multiple scales and fidelities facilitate the continuing improvements of PEM electrolyzer development to improve performance, cost, and reliability. The model scopes and methods are presented in this paper, which also introduces a comprehensive PEM electrolysis modeling tool based on computational fluid dynamics (CFD) software, ANSYS/Fluent. The modeling tool incorporates electrochemical model of a PEM electrolysis cell to simulate the performance of coupled thermal-fluid, species transport, and electrochemical processes in a product-scale cell or stack by leveraging the powerful meshing generation and CFD solver of ANSYS/Fluent. The thermal-fluid modeling includes liquid water/gas two-phase flow and simulates a PEM electrolysis cell by using Fluent user-defined functions as add-on modules accounting for PEM-specific species transport and electrochemical processes. The modeling outcomes expediate PEM electrolyzer scaling up from basic material development and laboratory testing.
08 HYDROGEN↗
Charge state-dependent ion condensation near conjugated polymer backbones
Despite the technological appeal of polymeric organic mixed ionic/electronic conductors (OMIECs) for diverse applications, a deep understanding of the fundamentals of mixed charge transport in these materials, especially regarding the complex interplay between polymer, ion and solvent structure in determining transport, is lacking. Herein, extensive molecular dynamics (MD) simulations of a model OMIEC representing various electrochemically gated states are reported that reveal charge state-dependent counterion condensation. X-ray diffraction simulations based on the MD data predict a measurable change in the scattering intensity at the counterion absorption edge, indicative of counterion repositioning with charging. We leverage an operando resonant X-ray scattering technique to experimentally corroborate the simulated scattering and report excellent agreement between predicted and experimental data, confirming that counterions preferentially reside in the lamellar mid-plane of crystallites at low doping, and near the polymer backbone at higher doping. Driving forces for ion type-dependent spatial repositioning and implications thereof are discussed.
36 MATERIALS SCIENCE↗
Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments
Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.
42 ENGINEERING↗
Electrochemical Modeling of Iodide Oxidation in Metal-Halide Molten Salts
The oxidation of iodide in NaI-AlBr 3 , NaI-AlCl 3 , and NaI-GaCl 3 molten salts was analyzed using simulation software to extract relevant kinetic parameters. The experimental oxidation potentials were ordered AlCl 3 < AlBr 3 < GaCl 3 , with higher oxidation potentials correlating with softer Lewis acidity of the metal halide. An iodide oxidation and metal halide speciation model was developed and simulated to fit the electrochemical response, enabling determination of electrochemical charge transfer parameters and chemical equilibrium constants. NaI-AlBr 3 displayed the fastest electron transfer rates yet showed the lowest current densities. All salts revealed smaller than expected current densities, explained by equilibrium between various species, where some are not electrochemically active at the studied potentials. These equilibrium reactions are due to the various metal halide species, controlling the reactant concentration of iodide and the resultant current. We hypothesize the electrochemically active iodide species, present as a metal halide monomer (MX 3 I − ), is decreased dramatically from the expected concentration, sequestered as a more stable metal halide dimer species (M 2 X 6 I − ) with a higher oxidation potential. Traditional Tafel analysis of the experimental data supports the validity of the simulations. These results increase understanding of iodide oxidation in low-temperature Lewis acidic molten salts and inform task-specific molten salt design.
Lee, Rose Y. (ORCID:0000000226587220)↗
Extended Finite Element Based Approach in Additive Manufacturing Modeling for Optimizing Highly Complex Manifold in Protonic Ceramic Electrochemical Cells
The objective of this project is to develop a novel numerical approach based on the extended finite element method (XFEM) for modeling moving boundaries of material deposited in additive manufacturing (AM) process with improved accuracy and reduced computational cost. One of the major challenges in simulating AM processes is that the boundaries of the computational domain need to change as material is deposited. Previous approaches have achieved this by activating new finite elements, but this requires a high level of mesh refinement to capture small movement of the boundary. XFEM allows the solution boundary to move independently of the mesh, permitting smooth representation of the evolution of the boundary as material is deposited. As a result, this new method will significantly improve the simulation accuracy while reducing the computational cost compared with existing approaches. Furthermore, this project will bring significant improvements for the design and AM process optimization in an iterative fashion among numerical simuation, parameter optimization, and experimental validation. Specifically, this project will support the protonic ceramic electrochemical cells (PCEC) stack development at INL by providing insights into the density, residual stress, thermal and mechanical properties of the PCEC manifold and interconnect, which is critical in enhancing the system lifetime and reducing the overall cost. Upon the success of proposed development and validation of the proposed approach, simulation will be applied to determine an optimal set of AM process parameters for the PCEC manifold production. This project will consolidate Idaho National Laboratory (INL)’s simulation capabilities provided by the MOOSE framework and the Valhalla AM simulation application, with encouraging expansion to emerging PCEC applications and AM technologies at INL and beyond.
97 MATHEMATICS AND COMPUTING↗
Investigating the Impact of Thickness, Calendering and Channel Structures of Printed Electrodes on the Energy Density of LIBs - 3D Simulation and Validation
Current lithium ion batteries (LIBs) are expensive and bulky, limited by relatively low charging rates. To increase the rate of charging and reduce weight, thin electrodes with high energy density are required. The increase in energy density can be achieved by several techniques including boosting electrolyte transport, high loading/utilization of active material, employing high conductive electrolytes and electrodes with advanced architectures, and increasing cell temperature. In this paper, a 3D physics-based electrochemical model of LIBs is developed in COMSOL simulation software for different thickness, calendering steps as well as channel structures (conical, cylindrical) to optimize the electrode design and in turn maximize volumetric energy density. The simulation results demonstrated that calendering the electrodes with high initial porosity increases the volumetric energy density of the cell. In addition, cylindrical channel structures with relatively lower edge-to-edge distance also results in increased volumetric energy density. The simulation results of the 3D model was validated by comparing it with experimental results.
improving volumetric energy density↗
Numerical Simulation and Modeling of Hydrogen Gas Evolution on Planar and Microwire Array Electrodes
The impact of gas evolution on the electrochemical characteristics of planar electrodes and microwire array electrodes has been analyzed using modeling and simulation. The impacts can mainly be broken into three phenomena: a) a shift in the local reversible hydrogen electrode potential; b) hyperpolarization; and c) an increase in the solution resistance of the electrolyte. The local reversible hydrogen electrode potential shift was found to play the most important role, constituting >40% of the total potential drop between the cathode and reference electrode, following correction for cell resistance. Compared to planar electrodes, a microwire array structure reduces the impact of bubbles on the solution conductance, but the shift in the local reversible hydrogen electrode potential varies with distance from the actual electrode surface.
25 ENERGY STORAGE↗