Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electrochemical 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mechanistic Insights on Permeation of Water over Iron Cations in Nanoporous Silicon Oxide Films for Selective H 2 and O 2 Evolution

Electrocatalysts encapsulated by an ultrathin and semipermeable oxide layer offer a promising avenue for efficient, selective, and cost-effective production of hydrogen through photoelectrochemical water splitting. This architecture is especially attractive for Z-scheme water splitting, for which a nanoporous oxide film can be leveraged to mitigate undesired, yet kinetically facile, reactions involving redox shuttles, such as aqueous iron cations, by limiting transport of these species to catalytically active sites. In this work, molecular dynamics simulations were combined with electrochemical measurements to provide a mechanistic understanding of permeation of water and Fe(III)/Fe(II) redox shuttles through nanoporous SiO 2 films. It is shown that even for SiO 2 pores with a width as small as 0.8 nm, water does not experience any energy barrier for permeating into the pores due to a favorable interaction with hydrophilic silanol groups on the oxide surface. In contrast, permeation of Fe(III) and Fe(II) into microporous SiO 2 pores is limited due to high energy barriers, which stem from a combination of distortion and dehydration of the second and third ion solvation shells. Furthermore, our simulations and experimental results show that SiO 2 coatings can effectively mitigate undesired Fe(III)/Fe(II) redox reactions at underlying electrodes by attenuating permeation of iron cations, while allowing water to permeate and thus participate in water splitting reactions. In a broader context, our study demonstrates that selectivity of solvated cations can be manipulated by controlling the pore size and surface chemistry of oxide films.

14 SOLAR ENERGY↗
DOE OSTI2024

In situ formed catalysts for active, durable, and thermally stable ammonia protonic ceramic fuel cells at 550 °C

Ammonia protonic ceramic fuel cells (NH 3 -PCFCs) are promising and attractive energy-conversion devices owing to their high energy density, zero-carbon emission, and safety. The development of NH 3 -PCFCs, however, depends largely on the insufficient activity and poor durability of typical Ni-based anodes for ammonia decomposition, especially at low temperatures such as 550 °C. Herein, we report a self-assembled heterostructured Ru 0.95 Cu 0.05 Ni x (RCN) catalyst obtained through an in situ reaction between the surface-decorated Ru 0.95 Cu 0.05 nanoparticles and the Ni grain in the anode under typical processing conditions. At 550 °C, Ni–BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3 anode-supported PCFCs with RCN catalysts exhibit a high peak power density of 0.732 W cm -2 and a significantly enhanced durability of 100 h in NH 3 . Moreover, the cells demonstrate improved thermal stability compared with the bare cell during a 31-cycle thermal cycling test in NH 3 between 550 and 700 °C. In conclusion, the enhanced performance is likely attributed to the synergistic effects of Ru and Cu in RCN for NH 3 decomposition, resulting in a more vital interaction of NH 3 than that of the bare anode surfaces, as confirmed by NH 3 thermal conversion, electrochemical performance, and theoretical simulations.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Exploring Li-Ion Transport Properties of Li 3 TiCl 6 : A Machine Learning Molecular Dynamics Study

We performed large-scale molecular dynamics simulations based on a machine-learning force field (MLFF) to investigate the Li-ion transport mechanism in cation-disordered Li 3 TiCl 6 cathode at six different temperatures, ranging from 25°C to 100°C. In this work, deep neural network method and data generated by ab − initio molecular dynamics (AIMD) simulations were deployed to build a high-fidelity MLFF. Radial distribution functions, Li-ion mean square displacements (MSD), diffusion coefficients, ionic conductivity, activation energy, and crystallographic direction-dependent migration barriers were calculated and compared with corresponding AIMD and experimental data to benchmark the accuracy of the MLFF. From MSD analysis, we captured both the self and distinct parts of Li-ion dynamics. The latter reveals that the Li-ions are involved in anti-correlation motion that was rarely reported for solid-state materials. Similarly, the self and distinct parts of Li-ion dynamics were used to determine Haven’s ratio to describe the Li-ion transport mechanism in Li 3 TiCl 6 . Obtained trajectory from molecular dynamics infers that the Li-ion transportation is mainly through interstitial hopping which was confirmed by intra- and inter-layer Li-ion displacement with respect to simulation time. Ionic conductivity (1.06 mS/cm) and activation energy (0.29eV) calculated by our simulation are highly comparable with that of experimental values. Overall, the combination of machine-learning methods and AIMD simulations explains the intricate electrochemical properties of the Li 3 TiCl 6 cathode with remarkably reduced computational time. Thus, our work strongly suggests that the deep neural network-based MLFF could be a promising method for large-scale complex materials.

Selvaraj, Selva Chandrasekaran (ORCID:000000029023↗
DOE OSTI2025

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

Novel Superionic Na Conductors for Solid‐State Na batteries: Materials and Interface

This project aims to advance the fundamental understanding on how the cation or anion doping influence on the crystal structure, phase stability and ion transport of novel Na superionic conductors. Specifically, we focus on the Na 3-x A x SbS 4-y X y family solid electrolytes (A: cation dopant, X: anion dopant), a class of promising materials for all-solid-state sodium batteries owing to their high ionic conductivity and compositional tunability. By integrating materials synthesis, neutron scattering, and electrochemical characterization and theoretical simulations, this project sees to elucidate 1) how does the synthesis approaches (solid-state reaction, low-temperature reaction) on the phase purity and doping process; 2) how chemical doping (A: cation dopants such as Sn; X: anion dopants such as Se, F) modulate lattice dynamics and Na+ migration as well as interface stability toward Na metal. Specifically, this project’s proposed technology involves three objectives: (1) synthesize novel Na-ion conductors (anion or cation doped Na 3 SbS 4 ) with high phase purity and examine their structural features (e.g., phase stability); (2) understand Na-ion conductive properties of new Na-ion conductors through a combination of experimental characterizations and theoretical simulations; (3) investigate the interface stability of novel Na superionic conductors with electrode materials in solid-state Na batteries. These objectives will be achieved through combined experimental synthesis/characterizations and computational approaches. The obtained knowledge will provide a comprehensive synthesis-structure-property correlation, enabling the rational design of novel solid-state conductors with excellent conductivity and stability.

25 ENERGY STORAGE↗
DOE OSTI2026

Multi-physics Preconditioning for Thermally Activated Batteries

Thermal batteries, also known as molten-salt batteries, are single-use reserve power systems activated by pyrotechnic heat generation, which transitions the solid electrolyte into a molten state. The simulation of these batteries relies on multiphysics modeling to evaluate performance and behavior under various conditions. This paper presents advancements in scalable preconditioning strategies for the Thermally Activated Battery Simulator (TABS) tool, enabling efficient solutions to the coupled electrochemical systems that dominate computational costs in thermal battery simulations. We propose a hierarchical block Gauss-Seidel preconditioner implemented through the Teko package in Trilinos, which effectively addresses the challenges posed by tightly coupled physics, including charge transport, porous flow, and species diffusion. The preconditioner leverages scalable subblock solvers, including smoothed aggregation algebraic multigrid (SA-AMG) methods and domain-decomposition techniques, to achieve robust convergence and parallel scalability. Strong and weak scaling studies demonstrate the solver’s ability to handle problem sizes up to 51.3 million degrees of freedom on 2048 processors, achieving near sub-second setup and solve times for the end-to-end electrochemical solve. These advancements significantly improve the computational efficiency and turnaround time of thermal battery simulations, paving the way for higher-resolution models and enabling the transition from 2D axisymmetric to full 3D simulations.

25 ENERGY STORAGE↗
DOE OSTI2025

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

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

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

Molecular Dynamics Simulations of Supercritical Carbon Dioxide and Water using TraPPE and SWM4-NDP Force Fields

The increased levels of carbon dioxide (CO 2 ) emissions due to the combustion of fossil fuels and the consequential impact on global climate change have made CO 2 capture, storage, and utilization a significant area of focus for current research. In most electrochemical CO 2 applications, water is used as a proton donor due to its high availability and mobility and use as a polar solvent. Additionally, supercritical CO 2 is a promising avenue for electrochemical applications due to its unique chemical and physical properties. Consequently, understanding the interactions between water and supercritical CO 2 is of great importance for future electrochemical applications. Molecular dynamics (MD) simulation is a powerful tool that enables atomistic-resolution dynamics of molecular systems, which can complement and guide future experimental investigations. This study employed atomistic MD to study the cosolubilities, codiffusivities, and structure of supercritical CO 2 and water systems, with a polarizable water model (SWM4-NDP) and a nonpolarizable CO 2 model (TraPPE). Additionally, ab initio MD simulations were used to better understand how atomistic polarizable/nonpolarizable models compare to explicit modeling of electron densities. The polarizable water model exhibited substantial improvement in water-associated properties. In conclusion, we anticipate the development of a compatible polarizable CO 2 model to yield similar improvement, providing a pathway for realizing novel high-pressure electrochemical systems.

25 ENERGY STORAGE↗
DOE OSTI2025

C-C Coupling Mechanism on Cu(100) A Molecular Dynamics Study at 298K

The electrochemical reduction of carbon dioxide (CO2) into valuable fuels such as C1 (syngas, methane) and C2 (ethylene, ethanol) products is a key strategy for achieving a carbon-neutral economy. Computational studies of C-C coupling, a critical step in CO2 reduction, are essential for designing more efficient catalysts. However, simulating these processes under realistic electrochemical conditions, including temperature and solvent effects, is computationally demanding. In this work, we develop a machine learning-based atomistic potential to study CO2 reduction on Cu(100) surfaces, accounting for temperature and explicit water solvent effects. We compute thermodynamic free energies of the possible C-C coupling pathways, CO*+CO*→OCCO*, CO*+CHO*→OCCHO*, CO*+COH*→OCCOH*, CHO*+COH*→OHCCOH*, COH*-COH*→HOCCOH*, and CHO*-CHO*→OHCCHO*. Our results quantify the thermodynamic tendencies of these reactions and reveal that, in addition to the well-established CO* + CO* → OCCO* pathway, CHO* is a critical intermediate in the formation of C2 products on Cu(100). Furthermore, we demonstrate that the machine learning approach offers a cost-efficient framework for studying CO2 reduction on diverse catalysts under realistic electrochemical conditions.

machine learning↗
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 OSTI2026

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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 OSTI2025

Unveiling the Mechanism of Mn Dissolution Through a Dynamic Cathode‐Electrolyte Interphase on LiMn2O4

Abstract Understanding the formation and evolution of the cathode‐electrolyte interphase (CEI), which forms at the interface between the cathode and electrolyte, is crucial for revealing degradation mechanisms in cathode materials, especially for developing strategies to stabilize the interphase in the strongly oxidizing conditions that evolve at high operating voltages in next‐generation Li‐ion batteries. However, The present understanding of the CEI is challenged by its complex and dynamic nature. In this work, near‐edge X‐ray absorption fine structure spectroscopy, electrochemical characterization, and reactive molecular dynamics simulations are combined to reveal a mechanism for CEI formation and evolution above model LiMn 2 O 4 (LMO) thin‐film electrodes in contact with conventional carbonate‐based electrolytes. It is found that Mn dissolution from LMO can be understood in terms of repetitive Mn 3 O 4 formation and dissolution behavior during cycling, which is closely connected to electrolyte decomposition and a key aspect of the CEI formation and growth. The behavior of the CEI in this model system offers detailed insight into the dynamic chemistry of the interphase, underscoring the important role of electrolyte composition and cathode surface structure in interphase degradation.

Ou, Wenhan↗
DOE OSTI2023

Controlling swelling in mixed transport polymers through alkyl side-chain physical cross-linking

Semiconducting conjugated polymers bearing glycol side chains can simultaneously transport both electronic and ionic charges with high charge mobilities, making them ideal electrode materials for a range of bioelectronic devices. However, heavily glycolated conjugated polymer films have been observed to swell irreversibly when subjected to an electrochemical bias in an aqueous electrolyte. The excessive swelling can lead to the degradation of their microstructure, and subsequently reduced device performance. An effective strategy to control polymer film swelling is to copolymerize glycolated repeat units with a fraction of monomers bearing alkyl side chains, although the microscopic mechanism that constrains swelling is unknown. Here we investigate, experimentally and computationally, a series of archetypal mixed transporting copolymers with varying ratios of glycolated and alkylated repeat units. Experimentally we observe that exchanging 10% of the glycol side chains for alkyl leads to significantly reduced film swelling and an increase in electrochemical stability. Through molecular dynamics simulation of the amorphous phase of the materials, we observe the formation of polymer networks mediated by alkyl side-chain interactions. When in the presence of water, the network becomes increasingly connected, counteracting the volumetric expansion of the polymer film.

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

Defluorination Mechanisms and Real-Time Dynamics of Per- and Polyfluoroalkyl Substances on Electrified Surfaces

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants found in groundwater sources and a wide variety of consumer products. In recent years, electrochemical approaches for the degradation of these harmful contaminants have garnered a significant amount of attention due to their efficiency and chemical-free modular nature. However, these electrochemical processes occur in open, highly non-equilibrium systems, and a detailed understanding of PFAS degradation mechanisms in these promising technologies is still in its infancy. To shed mechanistic insight into these complex processes, we present the first constant-electrode potential (CEP) quantum calculations of PFAS degradation on electrified surfaces. These advanced CEP calculations provide new mechanistic details about the intricate electronic processes that occur during PFAS degradation in the presence of an electrochemical bias, which cannot be gleaned from conventional density functional theory calculations. We complement our CEP calculations with large-scale ab initio molecular dynamics simulations in the presence of an electrochemical bias to provide time scales for PFAS degradation on electrified surfaces. Taken together, our CEP-based quantum calculations provide critical reaction mechanisms for PFAS degradation in open electrochemical systems, which can be used to prescreen candidate material surfaces and optimal electrochemical conditions for remediating PFAS and other environmental contaminants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗