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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Exploring Li-Air Batteries for High Specific Energy and High Power Applications: A Simulation Study

Commercialization of lithium-air batteries face many challenges, such as electrolyte decomposition, short cycle life, low energy efficiency, low power density, etc. However, commercialization of Li-air batteries for mass sensitive applications such as electric vehicles, portable power source, and drones is more challenging due to additional constraints of safety, electrolyte evaporation, high specific energy requirements, and reliable discharge times. In this presentation, we will present our finite element simulation results comparing Li-O2 and Li-air batteries using power density, energy density, specific power, and discharge times as metrics to evaluate different electrolytes and electrode geometry to reduce total mass and maximize discharge current. We use a finite element model and a discharge product model developed in which is based on porous electrode and concentrated electrolyte theories and the discharge product is modeled using quantum tunneling; for reaction kinetics and oxygen diffusion, an improved model was used. The electrolyte properties such as ion conductivity and ion diffusion were obtained from Molecular Dynamics (MD) simulations while the other parameters for the finite element model were calibrated to match experiments at high discharge current densities (>1.5 mA/cm2). The mass densities of different electrolytes were computed using MD simulations as well. For this presentation, we examine the practical electrochemical mass of a system at different current ratings, the sensitivity of mass to the use of ambient air as compared to pure oxygen as well as the electrolyte, which affects maximum current density and the total mass associated with the electrolyte (which includes the mass of additional components), and, the optimization of battery geometry for total discharge time, average discharge voltage, maximum discharge current, and minimum electrochemical mass.

Mehta, M.↗
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 OSTI2025

Unraveling membrane electrode assembly design for electrochemical conversion of carbon dioxide to formate/formic acid

This work presents a one-dimensional continuum modeling approach to investigate various cell architectures used for electrochemical conversion of CO 2 to formate/formic acid. Ion transport is simulated by a system of generalized modified Poisson–Nernst–Planck (GMPNP) equations that reflect the reactive transport phenomena including steric effects as the electrolyte solutions become concentrated. In the cathode catalyst layer, ionic current contributions from both the supporting electrolyte and solid-state ionomer are considered. Voltage and CO 2 utilization breakdowns are utilized to deconvolute the impacts of the cell architecture. The origins of (bi)carbonate formation in the cathode are explored, as the subsequent decrease in CO 2 availability is a key reason for low faradaic efficiencies to formate/formic acid. In addition, the role of a supporting electrolyte (KOH) is investigated to understand its tradeoffs: while the K + ions can improve both conductivity and electrochemically active surface area in the cathode, the presence of OH − ions raises the pH and leads to deleterious formation of (bi)carbonates. To this end, we also present parametric studies on the concentration and flow rate of supplied KOH to the cell, to establish a path towards eliminating the need for a supporting electrolyte.

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

Onset of instabilities in rotating flows by direct numerical simulation

A rotating disk is the canonical experiment for measuring surface reaction rates in geochemical and electrochemical systems. Using the similarity solution for laminar flow around an infinite disk, the mass transfer coefficient can be simply related to the intrinsic reaction rate at the surface. However, measurements of mass transfer rates use a finite-size disk within a larger container of solution; here the flow is no longer strictly laminar, but there must always be some recirculation. Our interest was initially in the assumption of a uniform radial concentration field, how this breaks down near the perimeter of the disk, and what effect that might have on the measured mass transfer rates. However, our numerical simulations suggest that the flow around a finite-size disk becomes time dependent at Reynolds number ($\textit{Re}$) below 1000, which is much smaller than the typical values in mass-transfer measurements ($\textit{Re}$~10 4 ). In this work, we observe the formation of coherent structures in the flow, which suggest a non-uniform mass transfer at the disk surface. The rotating-disk flow follows a similar sequence of instabilities to the Taylor–Couette flow: a centrifugal instability leading an axisymmetric, time-invariant flow, followed by a Hopf bifurcation to a time-periodic flow. To minimise the possibility that our results are a numerical artefact, we have also simulated the instability in the stationary boundary layer of a rotor–stator flow, comparing with self-similar solutions at low $\textit{Re}$ and with spectral methods near the critical Reynolds number.

42 ENGINEERING↗
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 OSTI2024

Functionalization of Cathode–Electrolyte Interface with Ionic Liquids for High-Performance Quasi-Solid-State Lithium–Sulfur Batteries: A Low-Sulfur Loading Study

We introduce a quasi-solid-state electrolyte lithium-sulfur (Li–S) battery (QSSEB) based on a novel Li-argyrodite solid-state electrolyte (SSE), Super P–Sulfur cathode, and Li-anode. The cathode was prepared using a water-based carboxymethyl cellulose (CMC) solution and styrene butadiene rubber (SBR) as the binder while Li6PS5F0.5Cl0.5 SSE was synthesized using a solvent-based process, via the introduction of LiF into the argyrodite crystal structure, which enhances both the ionic conductivity and interface-stabilizing properties of the SSE. Ionic liquids (IL) were prepared using lithium bis(trifluoromethyl sulfonyl)imide (LiTFSI) as the salt, with pre-mixed pyrrolidinium bis(trifluoromethyl sulfonyl)imide (PYR) as solvent and 1,3-dioxolane (DOL) as diluent, and they were used to wet the SSE–electrode interfaces. The effect of IL dilution, the co-solvent amount, the LiTFSI concentration, the C rate at which the batteries are tested and the effect of the introduction of SSE in the cathode, were systematically studied and optimized to develop a QSSEB with higher capacity retention and cyclability. Interfacial reactions occurring at the cathode–SSE interface during cycling were also investigated using electrochemical impedance spectroscopy, cyclic voltammetry, and X-ray photoelectron spectroscopy supported by ab initio molecular dynamics simulations. This work offers a new insight into the intimate interfacial contacts between the SSE and carbon–sulfur cathodes, which are critical for improving the electrochemical performance of quasi-solid-state lithium–sulfur batteries.

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 OSTI2025

Print-and-Plate Architected Electrodes for Electrochemical Transformations Under Flow

Flow cell electrodes are typically composed of porous carbon materials, such as papers, felts, and cloths. However, their random architecture hinders the fundamental characterization of electrode structure-performance relationships during in situ operation of porous electrochemical flow systems. This work describes a “print-and-plate” method that combines direct ink writing of micro-periodic lattices with a two-step metal plating process that converts them into highly conductive (sheet resistance 40 mΩ sq -1 ) electrodes. Their operando performance is assessed in an anthraquinone disulfonic acid half-cell using widefield electrochemical fluorescence microscopy, where output current and fluorescence intensity are in excellent agreement. The pressure drop associated with flow through three electrode designs is determined via simulations from which the most efficient design is identified and manufactured via print-and-plate. Confocal fluorescence microscopy is then used to create a 3D map of the state of charge (SOC) inside this print-and-plate electrode. The experimental state of the charge map is in good agreement with computational predictions. The rapid design, simulation, and fabrication of print-and-plate electrodes enable fundamental investigations of how architected porosity affects electrochemical performance under flow.

3D-printing↗
DOE OSTI2025

Validating corrosion models: A comparison of governing equations

Experimental validation of Finite Element Method (FEM) models varying electrochemical governing equations, inclusion of chemical reactions, and time on the resultant damage profile for two galvanic couples is explored. Two anode materials (Magnesium AZ31 and Carbon Steel) in contact with a cathode (Stainless Steel 304 L) were modeled in/exposed to NaCl (1 and 0.1 M respectively for the anode materials) for up to one week. The physics approach, inclusion of chemical reactions, and the boundary conditions required to accurately represent the damage profile in FEM models depended on the galvanic couple materials and, ultimately, the corrosion rate. For high rates of corrosion (i.e., magnesium anode), the Nernst-Planck equation with Electroneutrality was sufficient to describe the damage, while, for low rates of corrosion (i.e., carbon steel anode), the Laplace equation was sufficient. In all cases, the most complete governing equation (Nernst-Planck-Poisson Equation) was not necessary to accurately describe the damage. Precipitation reactions in solution also played a critical role in the predicted damage profile, especially for high corrosion rate systems. Finally, for short time periods (< 6 h), the choice of governing equations does not significantly influence damage profile results. Overall, the choice of physics to reduce error in simulations relies on the boundary conditions, geometry, conductivity of the solution, electrochemical potential differences, and time of exposure. The above results are discussed with regard to accuracy and computational savings.

Carbon steel↗
DOE OSTI2023

Purification of Lithium-Ion Battery Black Mass through Tailored Alkaline Corrosion

Obtaining high-purity material outputs is crucial to the viability of novel process aimed at direct recycling of lithium-ion batteries. Metallic impurities in recycled cathodes have been shown to inhibit performance, thereby threatening mainstream acceptance of recycled battery products. Thus, shredded black mass (BM) must be purified to remove metallic contaminants, and specifically Al and Cu originating from the electrode current collectors. We herein explore a process to ionize solid copper and aluminum to ionic form based on tailored alkaline chemistry, without incurring damage to the target cathode material (Li(NixMnyCo1-x-y)O2; NMC). Al and Cu corrosion may be enhanced through the addition of chloride salt, elevated temperatures, and the use of ultrasonication - all of which disrupt the formation of passivating films on the metallic surface, and thereby increase corrosion rate. We demonstrate optimized parameters for Al and Cu corrosion both from a kinetic and overall process cost perspective. Further, we analyze the impact of these conditions on the structural (XRD, SEM), chemical (EDS, ICP), and electrochemical (impedance, cycling, dQ/dV) properties of NMC, and suggest that the present purification method does not significantly disrupt NMC performance. Finally, we present preliminary results from a promising bench-scale demonstration of this purification process applied to a simulated black mass.

ADVANCED PROPULSION SYSTEMS,ENERGY STORAGE↗
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 OSTI2025

The Effects of Morphology and Hydration on Anion Transport in Self-Assembled Nanoporous Membranes

Ordered nanoporous polymer membranes offer opportunities for systematically probing the mechanisms of ion transport under confinement and for realizing useful materials for electrochemical devices. Here, we examine the impact of morphology and ion hydration on the transport of hydroxide and bromide anions in nanostructured polymer membranes with 1 nm scale pores. We use aqueous lyotropic self-assembly of an amphiphilic monomer, with a polymerizable surfactant to create direct hexagonal (H I ) and gyroid mesophases. UV-induced cross-linking leads to the formation of nanoporous polymers with water continuous channels. The membranes are mechanically robust and chemically durable, resisting degradation during extended exposure to 1 M NaOH solutions. We use a combination of electrochemical impedance spectroscopy, pulsed-field gradient NMR spectroscopy, and molecular simulations to elucidate anion and water transport. The as-prepared hexagonal systems display higher conductivity and lower activation energies for both anions relative to the gyroid system. When compared at equivalent hydration, however, gyroid and hexagonal membranes show similar activation energies, with nearly identical conductivities at ambient temperatures. Both ionic conductivity and water diffusivity increase with increasing hydration. Here, the water uptake as a function of relative humidity for the hexagonal and gyroid mesophases ultimately dictates the water diffusion and magnitude of the ionic conductivity, with the hexagonal system showing overall higher capacity for hydration and thus faster ion transport. The durability of these materials under aggressive alkaline conditions and their relatively high hydroxide ion conductivity suggest that these nanostructured polymers could be of interest as membranes for alkaline fuel cells.

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

Quantitative Voltammetry Measurements of High-Concentration Actinides in Molten Chloride Fuel Salts

Molten chloride mixtures containing high concentrations of UCl 3 (i.e., >70 wt%) are candidates for next-generation molten salt reactor (MSR) fuel salts that enable new plant optimization and fuel cycle strategies. MSR licensing, however, is contingent on the development of effective material control and accounting (MC&A) practices for these salts. To close this gap, we have developed an electroanalytical approach that enables accurate in situ measurements of uranium concentrations in these highly loaded fuel salts. This methodology uses a multielectrode array voltammetry sensor combined with digital simulations that account for non-idealities such as uncompensated resistance effects. The theory-based numerical simulations of voltammetry responses yielded correction factors that are agnostic of electrochemical cell geometry and produced consistent voltammetry peak current results when correcting for varying amounts of ohmic resistance. In doing so, we have demonstrated accurate concentration measurements for UCl 3 in LiCl-KCl-UCl 3 across a range of mass fractions between 55 and 72.3 wt% at 550 °C.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Solvent effects determine the sign of the charges of maximum entropy and capacitance at silver electrodes

Fully harnessing electrochemical interfaces for reactions requires a detailed understanding of solvent effects in the electrochemical double layer. Predicting the significant impact of solvents on entropic and electronic properties of electrochemical interfaces has remained an open challenge of computational electrochemistry. Using molecular dynamics simulations of silver–water and silver–acetonitrile interfaces, we show that switching the solvent changes the signs for both the charge of maximum capacitance (CMC) and charge of maximum entropy (CME). Contrasting the capacitance and CME behavior of these two interfaces, we demonstrate that the preferred orientation of the solvent molecule and the corresponding charge density determine the sign of the CMC and CME and, hence, the qualitatively different charge asymmetry of the electrochemical interface.

Chemistry↗
DOE OSTI2025

Unveiling the lithium-ion transport mechanism in Li{sub 2}ZrCl{sub 6} Solid-State Electrolyte {ital via} deep learning-accelerated molecular dynamics simulations.

Lithium zirconium chlorides (LZCs) present a promising class of cost-effective solid electrolytes for next-generation all-solid-state batteries. The unique crystal structure of LZCs plays a crucial role in facilitating lithium-ion mobility, which further affects the electrochemical performance. To understand the underlying mechanism governing ion transport, we employed deep learning-accelerated molecular dynamics simulation on Li2ZrCl6 (trigonal alpha- and monoclinic beta-LZC), focusing specifically on the zirconium coordination environment. Our results reveal that disordered alpha-LZC exhibits the highest ionic conductivity, while beta-LZC demonstrates significantly lower conductivity, closely aligning with experimental findings. The study confirms that across all phases, lithium migration proceeds via the site-to-site hopping mechanism, where variations in site residence times critically impact the overall ionic conductivity. In alpha-LZCs, lithium ions prefer to anisotropically diffuse across interlayers as the result of a lower energy barrier, driven primarily by collective diffusion. In contrast, lithium ions in beta-LZC primarily isotropically diffuse within the intralayer, hindered by higher energy barriers and determined by individual diffusion. The variation in ZrCl6 2- octahedral unit softening, induced by the specific layered arrangement of zirconium atoms, emerges as a critical determinant of the energy barriers across the LZC phases. These atomic-scale insights into the transport processes provide valuable guidance for the rational design and optimization of LZCs-based electrolytes, accelerating their practical application in advanced energy storage technologies.

Guo, Hanzeng↗
DOE OSTI2024

Three-chamber electrochemical reactor for selective lithium extraction from brine

Efficient lithium recovery from geothermal brines is crucial for the battery industry. Current electrochemical separation methods struggle with the simultaneous presence of Na + , K + , Mg 2+ , and Ca 2+ because these cations are similar to Li + , making it challenging to separate effectively. We address these challenges with a three-chamber reactor featuring a polymer porous solid electrolyte in the middle layer. This design improves the transference number of Li + (t Li+ ) by 2.1 times compared to the two-chamber reactor and also reduces the chlorine evolution reaction, a common side reaction in electrochemical lithium extraction, to only 6.4% in Faradaic Efficiency. Employing a lithium-ion conductive glass ceramic (LICGC) membrane, the reactor achieved high t Li+ of 97.5% in LiOH production from simulated brine, while the concentrations of Na + K + , Mg 2+ , and Ca 2+ are below the detection limit. Electrochemical experiments and surface analysis elucidated the cation transport mechanism, highlighting the impact of Na + on Li + migration at the LICGC interface.

Science & Technology - Other Topics↗
DOE OSTI2021

Probing the electronic properties of the electrified silicon/water interface by combining simulations and experiments.

Silicon (Si) is broadly used in electrochemical and photoelectrochemical devices, where the capacitive and Faradaic reactions at the Si/water interfaces are critical for signal transduction or noise generation. However, probing the electrified Si/water interface at the microscopic level remains a challenging task. Here we focus on hydrogenated Si surfaces in contact with water, relevant to transient electronics and photoelectrochemical modulation of biological cells and tissues. We show that by carrying out first principles molecular dynamics simulations of the Si(100)/water interface in the presence of an electric field we can realistically correlate the computed flat-band potential and tunneling current images at the interface with experimentally measured capacitive and Faradaic currents. Specifically, we validate our simulations in the presence of bias by performing pulsed chronoamperometry measurements on Si wafers in solution. Consistent with prior experiments, our measurements and simulations indicate the presence of voltage-dependent capacitive currents at the interface. We also find that Faradaic currents are weakly dependent on the applied bias, which we relate to surface defects present in newly prepared samples.

Ye, Zifan↗
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↗