Engineering Papers⌕ Search

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.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

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)↗
DOE OSTI2020

Modeling strategy for progressive failure prediction in lithium-ion batteries under mechanical abuse

The prediction for the internal failure of lithium-ion batteries (LIBs) under external mechanical abuse loading remains a challenge for safe design. This paper systematically studies the modeling approach for progressive failure simulation and short-circuit prediction. Hemispherical indentation experiments are conducted on a battery specimen containing 5 representative sandwich (RS) layers to examine the deformation and progressive failure behavior of battery components. Three different models, viz, a high-fidelity detailed model, an intermediate homogenized model and a fully homogenized model, are developed and employed to simulate the progressive failure in a multi-layer battery specimen. The simulation results obtained using the three models all correlate well with experimental phenomena, capturing the overall stress-strain response and ultimate failure of the battery structure. The numerical results can help us understand the failure behavior and the mechanism for mechanical abuse-induced short circuit in LIBs. The feasibility of the three models for use in progressive failure prediction is compared by analyzing the capability of the models in predicting the critical failure events, the computational efficiency of the models, and the feasibility to implement coupled mechanical-electrochemical-thermal simulations. The results of this work provide useful insights on practical choices for the modeling strategy and safety design of LIBs under mechanical abuse conditions.

25 ENERGY STORAGE↗
DOE OSTI2020

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↗
DOE OSTI2021

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↗
DOE OSTI2022

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↗
DOE OSTI2022

Effect of charge and anisotropic diffusivity on ion exchange kinetics in nuclear waste form materials

Motifs of radionuclide absorption particles (absorbers or fillers) for the nuclear waste treatment usually have crystal structures with nano-sized tunnels for fast ion diffusion. On the other hand, the accumulation of ions generates an electric field which affects the ion diffusion. Therefore, the strong anisotropic kinetic properties and the electric field affect ion exchange kinetics, especially in an assembly of these particles which may block the fast diffusion tunnels for each other. In this work, we developed a mesoscale model to describe the effect of anisotropic kinetic properties and electric field on ion exchange kinetics. With the model we simulated the effect of anisotropic diffusivity, electrochemical potential, particle morphology, size and aggregation on ion exchange kinetics during batch experiments. Here, tt is found that 1) ion exchange results in charge accumulation inside the particle because different ions have different diffusivities and different diffusion driving forces. The electric neutrality is not preserved; 2) the electric field speeds up the slower diffusion ions (either uptake or release) and slows down the faster diffusion ions to reduce charge accumulation and maintains charge neutrality; 3) the ion exchange becomes very slow at the later stages because a) diffusion driving force (electrochemical gradient, Cs + concentration in the solution and Na + inside the particle) continuously decreases, as observed in batch experiments; and b) the diffusivity decreases due to the Cs + concentration dependence of ion diffusivity. This can explain the observed experimental phenomena that the radionuclide capacity usually cannot be reached; 4) for needle-like particle and their clusters, both surface area and diffusion distance along the direction with the fastest diffusivities have significant impact on the ion exchange kinetics.

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

Heat Generation versus Heat Dissipation Maps for Batteries (Battery Safety Maps) [SWR-21-27]

NREL's Battery Safety Map (Heat Generation versus Heat Dissipation Maps for Batteries) generates a response surface that is grouped into regions where the battery is deemed safe from the propagation of thermal runaway, those where an engineering solution would render the battery safe, and finally instances where thermal runaway is inevitable. The demarcation for the different regions is provided by specifications or constraints imposed by the end user. NREL's Battery Safety Map is a user-defined subroutine that sweeps through a set of parameters describing the operating conditions and short-circuit scenarios for battery safety simulations. For each run, the electrochemical-thermal model implementing reaction kinetics and gas generation are run to evaluate the amount of heat generated. A thermal model then captures the amount of heat dissipation based on boundary conditions defined (e.g., forced convection, constant temperature at the surface, etc.). Based on the rate of heat transfer and the amount of heat generated, each run is classified into "safe," "mitigation possible," or "unsafe" zones. These categories are defined based on constraints or tolerance specifications provided by the end user.

Santhanagopalan, Shriram↗
NASA NTRS1976

Mathematical model of one-man air revitalization system

A mathematical model was developed for simulating the steady state performance in electrochemical CO2 concentrators which utilize (NMe4)2 CO3 (aq.) electrolyte. This electrolyte, which accommodates a wide range of air relative humidity, is most suitable for one-man air revitalization systems. The model is based on the solution of coupled nonlinear ordinary differential equations derived from mass transport and rate equations for the processes which take place in the cell. The boundary conditions are obtained by solving the mass and energy transport equations. A shooting method is used to solve the differential equations.

Source record↗
DOE OSTI2022

Probing the Role of Multi-scale Heterogeneity in Graphite Electrodes for Extreme Fast Charging

Electrode-scale heterogeneity can combine with complex electrochemical interactions to impede lithium-ion battery performance, particularly during fast charging. This research investigates the influence of electrode heterogeneity at different scales on the lithium-ion battery electrochemical performance under operational extremes. We employ image-based mesoscale simulation in conjunction with a three-dimensional electrochemical model to predict performance variability in 14 graphite electrode X-ray computed tomography data sets. Our analysis reveals that the tortuous anisotropy stemming from the variable particle morphology has a dominating influence on the overall cell performance. Cells with platelet morphology achieve lower capacity, higher heat generation rates, and severe plating under extreme fast charge conditions. On the contrary, the heterogeneity due to the active material clustering alone has minimal impact. Our work suggests that manufacturing electrodes with more homogeneous and isotropic particle morphology will improve electrochemical performance and improve safety, enabling electromobility.

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

Dendrite Growth Morphology Modeling in Liquid and Solid Electrolytes

The main goal of this project is to develop a multi-scale modeling approach that connects micron-scale phase-field models and atomic-scale density functional theory (DFT)-based simulations via parameter- and relationship-passing in order to predict Li-metal dendrite morphology evolution, in both liquid and solid electrolytes. The key hypothesis of the DFT-informed phase-field multiscale modeling approach is that it can capture the electrochemical-mechanical driving forces and incorporate the roles of nano-meter-thin solid electrolyte interphase (SEI) in liquid electrolytes as well as of the microstructures of micro-meter-thick solid electrolytes (SEs) for all-solid-state batteries. In this project, we have formulated and implemented phase-field models to incorporate the electrochemical driving forces in liquid electrolytes and then incorporate mechanical driving forces to simulate dendrite growth in solid electrolytes with resolved microstructures. We have implemented two treatments for the SEI: an explicit model to include the microstructure of the SE or SEI in the phase field model and an implicit model to simulate the impact of nano-meter thick SEI in liquid electrolytes by varying the electrode/electrolyte interfacial properties. The key interfacial properties, including the electronic and ionic transport properties, the charge transfer reaction kinetics, and mechanical properties, were computed by DFT-based calculations. At the DFT-based model, one key advancement is to directly predict the charge transfer reaction kinetics at a complex Li/SEI/electrolyte interface by linking DFT with density functional tight binding (DFTB) calculations. As the main accomplishments, we have demonstrated two successful predictions in both solid electrolyte and liquid electrolyte based on this multiscale approach. The predicted intergranular Li dendrite growth in LLZO revealed the importance of trapped electrons at internal interfaces in the microstructure of LLZO. The predicted electroplating morphology of mossy Li and faceted Mg agreed well with experiments. The insights provided by the multiscale model and the model enabled electrolyte and SEI design will accelerate the development of Li-metal electrode for high energy density batteries, that meet DOE’s target on cell density (>350 Wh/kg) and cost below $100/kWhuse for EV applications.

25 ENERGY STORAGE↗
DOE OSTI2023

Understanding the Effect of Electrochemical Properties and Microstructure on the Microgalvanic Corrosion of Mg Alloys via Phase-Field Simulations

In this work, we apply a phase-field model to gain insights into the corrosion behavior of Mg alloys. Specifically, we study the effect of electrochemical properties of the materials and their environment and the spatial distribution of second phases in a magnesium alloy. To this end, we perform sensitivity analyses in which we separately vary the exchange current density and corrosion potential of the second phase, as well as the electrolyte conductivity. We describe the effects of these parameters on the corrosion current and the resulting corrosion-front morphology. We also confirm that the ratio between the Wagner length, calculated for the anodic phase, and the initial width of the anodic domain is a good predictor of the morphological characteristic of the corroding surface. In addition, we examine the effect of the microstructure on the corrosion behavior in both 2D and 3D systems by varying the shape, size, and distribution of the second-phase regions. Our findings provide insights into the influence of the two-phase microstructure on the corrosion behavior in a magnesium alloy.

Electrochemistry↗
DOE OSTI2020

Transport Processes in a Li-ion Cell during an Internal Short-Circuit

Internal short-circuit in a lithium-ion cell causes an abrupt increase in cell temperature and triggers subsequent thermal runaway. In this work, we present a detailed electrochemical-thermal model to investigate the physical behavior during an internal short-circuit. Simulations at wide range of heat transfer coefficients and short-circuit resistances are conducted to illustrate electrochemical and thermal behavior under a wide range of conditions. The Joule heating at the shorted region promotes electrochemical reactions nearby, causing in-plane non-uniformity of electrolyte and active material transport. Furthermore, it is found that diffusion in solid active materials plays a significant role at very high shorting currents (~20 C), as electrochemical reactions rate are being controlled progressively more by availability of Li + at the interface, due to limitations in diffusion through the active material with increasing discharge rates. This diffusion limitation causes a drop in available energy, and subsequently a decrease in cell equilibrium potential and the heat generation rate at the location of the short. On the other hand, rapid depletion of lithium concentration in the electrolyte and accumulation on the electrode surface results in highly non-uniform transport properties resulting in higher heat generation rates. Hence, the heating regime shifts from "local heating" to "global heating". Based on the findings, important design parameters for battery safety are discussed.

25 ENERGY STORAGE↗
DOE OSTI2020

Microstructure Generation via Generative Adversarial Network for Heterogeneous, Topologically Complex 3D Materials

Using a large-scale, experimentally captured 3D microstructure dataset, we implement the generative adversarial network (GAN) framework to learn and generate 3D microstructures of solid oxide fuel cell electrodes. The generated microstructures are visually, statistically, and topologically realistic, with distributions of microstructural parameters, including volume fraction, particle size, surface area, tortuosity, and triple phase boundary density, being highly similar to those of the original microstructure. These results are compared and contrasted with those from an established, grain-based generation algorithm (DREAM.3D). Importantly, simulations of electrochemical performance, using a locally resolved finite element model, demonstrate that the GAN generated microstructures closely match the performance distribution of the original, while DREAM.3D leads to significant differences. Finally, the ability of the generative machine learning model to recreate microstructures with high fidelity suggests that the essence of complex microstructures may be captured and represented in a compact and manipulatable form.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Understanding and Strategies for Controlled Interfacial Phenomena in Lithium-Ion Batteries and Beyond

Electrolyte chemistry and properties and electrode structure and chemical properties particularly at the interfaces are crucial for the development of advanced battery components. This project focused on elucidating the roles of the electrolyte and electrode on the formation and evolution of the SEI layer and cell electrochemical performance on silicon and on Li metal anodes. For Si nano/microstructures we evaluated lithiation, volume changes, reactivity, and chemo-mechanical transformations as functions of nanoparticle size, shape, presence of coatings, electrolyte composition, electron leakage to the electrolyte, and cycling. For Li metal anodes we provided detailed analyses of the electrode and electrolyte effects on interfacial reactivity and as a function of microstructural evolution and underlying stochasticity as dendrites nucleate and grow. We thoroughly investigated and simulated how the chemistry of the various components of the electrolyte and the electrode architecture may affect the electrochemical reactions as well as cell degradation. A multiscale modeling approach was utilized, where atomistic simulations informed about the microscopic behavior of the system, and intermediate time and length scales were investigated with mesoscopic models. Results were also tested against selected experiments from our collaborators. The combined theoretical-experimental strategy saves costs by utilizing predictions from computational analyses to guide experimentation.

25 ENERGY STORAGE↗
NASA NTRS1988

Nickel cadmium battery performance simulation - A means toward optimum design in satellite applications

It is noted that electrochemically based cell models can advance state-of-the-art battery prediction by allowing variations in the parameters that constitute wear-out and predict end-of-life performance. This modeling method is used as the foundation for the simulation of Ni-Cd battery performance. This simulation tool consists of a comprehensive computer model to be developed which would predict the voltage behavior of a battery under mission conditions.

Clark, K.↗
DOE OSTI2024

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, in this work, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Electric Field Effects on Water and Ion Structure and Diffusion at the Orthoclase (001)–Water Interface

Understanding the electrochemical properties of mineral–water interfaces tends to rely upon electrical double layer (EDL) models, but these models are based on the assumption that electrostatic equilibrium is constantly maintained. In reality, interfacial reactions, ion diffusion, and their electrochemical signatures are based in nonequilibrium conditions of locally or globally imbalanced electrical fields where current EDL models have limited purview. In this work, we performed molecular dynamics (MD) simulations of the orthoclase (001) surface in contact with a 1 M NaCl aqueous solution under various electric fields, to explore the interplay between EDL structure and dynamics when perturbed by electric fields of different direction and strength, by confinement, and by different distributions of structural surface charge. The simulations showed that confinement between two opposing (001) surfaces led to the development of an induced field when the applied field was perpendicular to the surfaces and, as a result, to ionic diffusion coefficients that were independent of electric field strength. In contrast, when the applied field was parallel to the surfaces, confinement resulted in ionic diffusion coefficients that were more strongly dependent on the magnitude of the electric field than in bulk water. Differences in the density and distribution of aluminol groups on the two surfaces had a significant impact on how the interfacial structure and dynamics varied in the presence of an electric field. Notably, these differences resulted in an electro-osmotic flow with opposite directions at the two surfaces under parallel applied electric field. Overall, the MD simulations highlighted the importance of considering atomic-level structure and heterogeneities when developing models of the electrochemical properties of mineral–water interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗