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

Multi-layer anodes for high-current applications

In this work, we present a combination of experiments and modeling of a two-layer anode structure designed by EnPower Inc. for high-energy and fast-charge capabilities. The anode consists of lower porosity near the current collector and higher porosity near the separator with comparable active materials in both regions. A pseudo-two-dimensional electrochemical–thermal model was designed to represent the performance of this electrode. Simulations, consistent with experiments, show superior ion transport and lithiation in the multi-layer anode (MLA) compared to a conventional single-layer anode (SLA). Surprisingly, this improved transport in MLA manifests as enhanced cathode performance during high-rate discharge and this, in turn, provides higher energy density for MLA. Similarly, during fast charge, less irreversible lithium is deposited due to this improved transport in MLA, and hence MLA exhibits less capacity fade compared to SLA. Polarization analysis demonstrates marginally lower cumulative overpotential for MLA in different case studies; however, MLA cells maintain a significantly higher capacity in the same conditions and have more than double the cycle life. Further, this means despite the apparently limited polarization benefit provided by MLA; the MLA structure can more reliably be employed in cell designs. Additional design changes are also analyzed by means of the model.

25 ENERGY STORAGE↗
DOE OSTI2022

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. Here, a vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

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

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Analytical-based simulation approach for an anion exchange membrane fuel cell

An analytical and empirical-based 1-D, non-isothermal, steady-state model for anion exchange membrane fuel cell capable of capturing two-phase phenomena is presented in this study. Coupled multi-physics including mass and charge transport, electrochemical reactions, heat transfer, and two-phase water transport are considered in the model and the simulated results are compared to experimental data. To better represent actual material properties and localized conditions, the model applies multilayer discretization in the gas diffusion electrode to enhance prediction accuracy. The model successfully predicts the baseline performance at 70 °C, 131 kPa abs., 92% RH with pure H 2 /O 2 gas as well as the limiting current at 10% H 2 . The robust simulation approach allows for simplistic and accurate estimation of cell performance without the complications of applying two-phase parameters and expensive computational need for numerical models. In addition, the results from the sensitivity studies of material properties and operating conditions provide valuable insights on water management strategies and optimal component design for advancing anion exchange membrane fuel cell technology.

1-D model↗
DOE OSTI2024

SEI Formation and Lithium-Ion Electrodeposition Dynamics in Lithium Metal Batteries via First-Principles Kinetic Monte Carlo Modeling

The stabilization and enhanced performance of lithium metal batteries (LMBs) depend on the formation and evolution of the Solid Electrolyte Interphase (SEI) layer as a critical component for regulating the Li metal electrodeposition processes. This study employs a first-principles kinetic Monte Carlo (kMC) model to simulate the SEI formation and Li + electrodeposition processes on a lithium metal anode, integrating both the electrochemical electrolyte reduction reactions and the diffusion events giving place to the SEI aggregation processes during battery charge and discharge processes. The model replicates the competitive interactions between organic and inorganic SEI components, emphasizing the influence of the cycling regime. Results indicate that grain boundaries within the SEI facilitate faster lithium-ion transport compared to crystalline regions, crucial for improving the performance and stability of LMBs. The findings underscore the importance of dynamic SEI modeling for further development of next-generation high-energy-density batteries.

25 ENERGY STORAGE↗