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 145 records · Page 8

Less can be more: Insights on the role of electrode microstructure in redox flow batteries from two-dimensional direct numerical simulations

Understanding how to structure a porous electrode to facilitate fluid, mass, and charge transport is key to enhancing the performance of electrochemical devices, such as fuel cells, electrolyzers, and redox flow batteries (RFBs). Here, using a parallel computational framework, direct numerical simulations are carried out on idealized porous electrode microstructures for RFBs. Strategies to improve an electrode design starting from a regular lattice are explored. By introducing vacancies in the ordered arrangement, it is possible to achieve higher voltage efficiency at a given current density, thanks to improved mixing of reactive species, despite reducing the total reactive surface. Careful engineering of the location of vacancies, resulting in a density gradient, outperforms disordered configurations. Our simulation framework is a new tool to explore transport phenomena in RFBs, and our findings suggest new ways to design performant electrodes.

25 ENERGY STORAGE↗

MOSCATO Development and Integration in Fiscal Year 2024

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 evolution for molten salt equipment. In FY24, significant developments were made to the code package, enhancing its capabilities in many aspects. The improvements and advancements can be summarized as follows: 1. Implementation of tritium transport capabilities and validation with experimental data: To enable modeling of tritium and other fission gases within MSRs, we implemented gas transport within MOSCATO via inclusion of couple mass transport equations within the salt and structural alloys. Comparisons to experimental data from literature showed good agreement with respect to tritium release rates. 2. Preliminary implementation of two-phase flow models in MOSCATO: To model tritium and other gases above their solubility limits, we implemented preliminary two-phase flow models within MOSCATO to account for bubble transport. The first model adopted was the Level-Set approach, which can handle the high void fraction regime, but with a requirement for high mesh resolution thus high computational expense. In this report, we present a verification of the Level-Set method using a simple benchmark case. We also performed a demonstration of the code as applied to an experimental case involving cover gas flow through salt in an experimental vessel. The second model adopted was the Eulerian-Eulerian dispersed flow model, which is computationally cheaper but limited to low void fraction regimes, such as bubbly flow. Validation and verification have not yet been performed for the Eulerian-Eulerian approach, but a preliminary implementation was completed. 3. Validation with static corrosion experiments: Static corrosion experimental data for stainless steel coupons within molten salts was used to further validate the corrosion model in MOSCATO. To do so, we leveraged the existing models in MOSCATO and simulated the sample mass loss and mass gain phenomena. Several ion species, including Cr 2+ , Fe 2+ and H + , were simulated in salt using the PNP solver, while Cr 0 and Fe 0 were simulated with a diffusion solver in stainless steel. The mass loss of the samples was compared with experimental data, and good agreement was achieved. These combined activities served to further expand the capabilities of MOSCATO and make it more generally applicable to the full range of phenomena that can control chemistry and corrosion in molten salt reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Insights into the Electrochemical Oxidation and Reduction of Nickel Oxide Surfaces

Surface oxidation/reduction processes, driven by varying electrochemical potentials, can substantially impact catalyst effectiveness and, consequently, electrolyzer performance. Here, this study combines theoretical and experimental approaches to explore the surface redox behavior of nickel oxides, which are cost-effective and efficient catalysts for many electrochemical reactions. Surface Pourbaix diagrams for three different phases of nickel oxides, i.e., nickel hydroxide (Ni(OH) 2 ), nickel oxyhydroxide (NiOOH), and nickel dioxide (NiO 2 ), were constructed using density functional theory-based simulations. Various experimental methods, including cyclic voltammetry, in situ Raman spectroscopy, and electrochemical titration, were employed to probe the surface redox processes of nickel oxide thin films. Our findings indicate that the ABAB stacking sequence of Ni(OH) 2 lacks stability under oxidizing conditions to host the surface oxidation (deprotonation) events, while the AABBCC stacking sequence of NiOOH is energetically favorable due to the presence of interlayer hydrogen bonding. Rapid charge transfer facilitated by interlayer hydrogen bonding accounts for the higher reactivity of partially oxidized/reduced NiOOH (001) surfaces compared to Ni(OH) 2 (001) and NiO 2 (001) surfaces with the same stoichiometry, where interlayer hydrogen bonding is absent. Insights presented in this work can offer guidelines for optimizing operational conditions and tailoring the surface structures and oxidation states of nickel oxides to enhance performance in applications such as electrocatalysis and supercapacitors.

36 MATERIALS SCIENCE↗

Siloxane Deposition on the Ni-YSZ Solid Oxide Fuel Cell Anode Exposed to Bio-Syngas

Siloxane, a common contaminant present in biogas, is known for adverse effects on cogeneration prime movers. In this study, the siloxane deposition products and mechanism in the solid oxide fuel cell nickel-yttria stabilized zirconia (Ni-YSZ) anode are investigated analytically and experimentally. An SOFC with Ni-YSZ anode and pure Ni/YSZ pellets were exposed to a simulated biogas-reformate fuel with octamethylcyclotetrasiloxane (D4) contamination at 750 °C. The electrochemical characterization results show that the SOFCs performance degradation caused by D4 contamination is irreversible. Morphology and XRD results illustrate that silicon and carbon deposition can both be detected in the anode and pellets. Graphite, SiC and SiO2 are all possible products based on the results of XRD test. Here, according to the formation of graphite and SiC, the new mechanism suggests that carbon is also an essential factor in siloxane contamination of Ni-YSZ anodes besides silicon, which can be explained by the catalytic and electrochemical analysis.

25 ENERGY STORAGE↗

Unlocking High Capacity and Reversible Alkaline Iron Redox Using Silicate-Sodium Hydroxide Hybrid Electrolytes

Alkaline iron (Fe) batteries are attractive due to the high abundance, low cost, and multiple valent states of Fe but show limited columbic efficiency and storage capacity when forming electrochemically inert Fe 3 O 4 on discharging and parasitic H 2 on charging. Herein, sodium silicate is found to promote Fe(OH) 2 /FeOOH against Fe(OH) 2 /Fe 3 O 4 conversions. Electrochemical experiments, operando X-ray characterization, and atomistic simulations reveal that improved Fe(OH) 2 /FeOOH conversion originates from (i) strong interaction between sodium silicate and iron oxide and (ii) silicate-induced strengthening of hydrogen-bond networks in electrolytes that inhibits water transport. Furthermore, the silicate additive suppresses hydrogen evolution by impairing energetics of water dissociation and hydroxyl de-sorption on iron surfaces. In conclusion, this new silicate-assisted redox chemistry mitigates H 2 and Fe 3 O 4 formation, improving storage capacity (199 mAh g -1 in half-cells) and coulombic efficiency (94 % after 400 full-cell cycles), paving a path to realizing green battery systems built from earth-abundant materials.

36 MATERIALS SCIENCE↗

Molecular-level solvation and selectivity behavior of Na + , K + , and Li + within glycerol-derived solvents

Ether-based acyclic solvents have a tendency to confine positively charged atoms, forming structures called podands. Due to the stability of these structures, podands have been identified as promising candidates for several different electrochemical applications. In this work, we use molecular dynamics simulations to investigate the solvation and selectivity of Na + , K + , and Li + by different glycerol-derived solvents. The structural properties suggest that the solvent-cation association is stable over a very wide temperature range. Further, a strong correlation is observed between the more electronegative oxygen sites and preferential interaction with the cations, which suggests that the electrostatic oxygen-ion interactions dominate the coordination between these species. Significantly more intense interactions are observed between the solvents and Li + , versus Na + or K + . These molecular-level thermodynamic insights provide fundamental support for further exploring the performance of glycerol-derived podand complexes in Li + capture applications and/or electrolyte formulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Parametric analysis of anodic degradation mechanisms for fast charging lithium batteries with graphite anode

In this work, we report the impact of the temperature-driven synergistically-coupled anodic degradation mechanisms on the electrochemical performance of lithium batteries with graphite anode over multiple cycles. Temperature dependence of electrochemical reactions and damage mechanisms, such as solid electrolyte interface (SEI) growth, lithium plating/stripping, dead lithium storage/dissolution, and film cracking are incorporated into the degradation model. Results of a parametric analysis are presented, evaluating the effects of charging rates (1–6 C), operating temperatures (-15 – 45°C) and electrode design parameters, on the relative performance fade in the lithium-ion battery. Thermo-electrochemical process maps are developed to provide insights into the relationship between electrode performance and failure mechanisms. The simulation results predict a severe capacity loss due to lithium plating at low temperatures, which is further aggravated at high charging rates. A common strategy for mitigating lithium plating, through charging at high temperatures, also results in rapid capacity loss due to accelerated SEI formation. Simulation results are used to identify the combination of operating conditions and electrode design parameters that improve the electrochemical performance of the battery. These results demonstrate an opportunity to design safe and high-performance lithium-ion batteries, guided by anodic degradation models.

Battery design↗

Electric-Field-Induced Assembly of an Ionic Liquid–Water Interphase Enables Efficient Heavy Metal Electrosorption

Controlling ion desolvation, transport, and charge transfer at the electrode-electrolyte interface (EEI) is critical to enabling the rational design of efficient and selective separation of targeted heavy metals and decontamination of industrial waste water. One challenge is to sufficiently resolve and interrogate the intermediate transformation steps between solvated metal cations and bare cations for electroreduction at the EEI and establish pathways to modulate these steps to achieve efficient energy transfer for targeted reactive separations. Herein, we obtained a predictive understanding of the effect of adding a hydrophobic ionic liquid, 1-ethyl-3-methylimidazolium chloride (EMIMCl), to an aqueous electrolyte on modulating the desolvation and electrosorption of Pb 2+ cations using a combination of experimental and theoretical techniques. Here, we reveal formation of a compact interphase layer consisting of EMIMCl-Pb complex clusters under an applied electric field using operando electrochemical Raman spectroscopy, atomic force microscopy, and electrochemical impedance spectroscopy measurements combined with classical molecular dynamics simulations. The application of a larger negative potential is shown to result in formation of a well oriented layer with the positive imidazolium ring of the EMIMCl lying on the electrode and the long hydrophobic alkyl chain extending into the bulk electrolyte. This oriented layer is demonstrated to facilitate desolvation of incoming solvated Pb 2+ cations and decrease the charge transfer resistance for Pb electrodeposition, which has important implications for selective removal of Pb from contaminated mixtures. Overall, our studies open up new opportunities to modulate ion desolvation using hydrophobic ionic liquids in aqueous electrolytes for efficient heavy metal separation.

36 MATERIALS SCIENCE↗

Electrochemical C–N Bond Formation within Boron Imidazolate Cages Featuring Single Copper Sites

Electrocatalysis expands the ability to generate industrially relevant chemicals locally and on-demand with intermittent renewable energy, thereby improving grid resiliency and reducing supply logistics. Herein, we report the feasibility of using molecular copper boron-imidazolate cages, BIF-29(Cu), to enable coupling between the electroreduction reaction of CO 2 (CO 2 RR) with NO 3 – reduction (NO 3 RR) to produce urea with high selectivity of 68.5% and activity of 424 μA cm –2 . Remarkably, BIF-29(Cu) is among the most selective systems for this multistep C–N coupling to-date, despite possessing isolated single-metal sites. The mechanism for C–N bond formation was probed with a combination of electrochemical analysis, in situ spectroscopy, and atomic-scale simulations. We found that NO 3 RR and CO 2 RR occur in tandem at separate copper sites with the most favorable C–N coupling pathway following the condensation between *CO and NH 2 OH to produce urea. In conclusion, this work highlights the utility of supramolecular metal–organic cages with atomically discrete active sites to enable highly efficient coupling reactions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Methylation enables the use of fluorine-free ether electrolytes in high-voltage lithium metal batteries

Lithium metal batteries represent a promising technology for next-generation energy storage, but they still suffer from poor cycle life due to lithium dendrite formation and cathode cracking. Fluorinated solvents can improve battery longevity by improving LiF content in the solid-electrolyte interphase; however, the high cost and environmental concerns of fluorinated solvents limit battery viability. Here, for this work, we designed a series of fluorine-free solvents through the methylation of 1,2-dimethoxyethane, which promotes inorganic LiF-rich interphase formation through anion reduction and achieves high oxidation stability. The anion-derived LiF interphases suppress lithium dendrite growth on the lithium anode and minimize cathode cracking under high-voltage operation. The Li + -solvent structure is investigated through in situ techniques and simulations to draw correlations between the interphase compositions and electrochemical performances. The methylation strategy provides an alternative pathway for electrolyte engineering towards high-voltage electrolytes while reducing dependence on expensive fluorinated solvents.

25 ENERGY STORAGE↗

Disentangling cation effects on ion mobility and structure in ionic liquid electrolytes

Ionic liquids (ILs) are low-temperature molten salts, where ion transport is primarily governed by ion–ion interactions. Yet, the impact of organic IL cations on critical electrolyte properties such as ion dissociation and overall transport behavior in lithium-salt-doped ILs remains poorly understood. Moreover, despite their critical role in designing IL-based electrolytes for energy storage applications, ion–ion interactions and ion-specific transport under an applied electrical potential are seldom quantified, largely due to the unique experimental and computational challenges involved. Herein, we compare transport properties obtained using 1 H, 7 Li, and 19 F pulsed-field gradient nuclear magnetic resonance (NMR) and electrophoretic NMR (eNMR) with those measured by electrochemical impedance spectroscopy. Non-equilibrium molecular dynamics (MD) simulations and eNMR confirm the presence of negatively charged [Li(TFSI) n ] (1−n) aggregates that migrate towards the positive electrode, resulting in negative lithium transference numbers. Equilibrium MD simulations reveal a vehicular Li ion transport mechanism facilitated by long-lived aggregates with Li + cations strongly bound to multiple TFSI − anions. Finally, we observe an inverse relationship between the apparent charge of the TFSI − anion in the neat IL, which is dictated by the IL cation, and Li + transport in the salt-doped systems. This highlights the opportunity to tune electrolyte performance by tailoring cation chemistry.

Li-ion batteries↗

Renewable Ammonia as An Energy Fuel for Ocean Exploration and Transportation

Renewable power generated from wave energy has faced technological and cost barriers to entry into utility-scale electricity markets. As an alternative, the production of chemical fuels, such as ammonia (NH3) which has high energy density (11.5 MJ/L) and facile storage properties, may open wave energy to new markets including ocean exploration and transportation. The electrochemical method has been studied to synthesize NH3 from air and water at ambient conditions. Based on some recent work on the electrochemical synthesis of NH3, it is possible to achieve an overall conversion efficiency of 10% from wave energy to NH3 through an electrochemical reaction between air and water. If all the recoverable wave energy in the United States (1170 TWh/yr) is used to produce renewable NH3 fuel replacing hydrocarbon fuels, this can help reduce over 300 million tons of CO2 emission every year. Several potential application scenarios at sea have been proposed for renewable NH3 fuel including production and storage for marine shipping and seasonal energy storage for Arctic exploration. Liquefied NH3 has much higher energy density, both gravimetric and volumetric, than a variety of batteries but the energy efficiency of NH3 is lower than modern batteries such as Li-ion. The Levelized cost of storing NH3 prepared using electric energy is less than $0.2/kWh and the storage time can exceed 10,000 hours which indicates that NH3 can be a promising energy storage solution to make use of abundant wave energy. However, there are some safety and environmental concerns involved in the usage of NH3 at sea. The challenges in the electrochemical catalyst for the NH3 synthesis and how molecular simulation may help to screen electrocatalyst with high efficiency and selectivity were also briefly discussed.

Liu, Jian↗

High-Fidelity Accelerated Design of High-performance Electrochemical Systems

Large-scale electrification is vital to addressing the climate crisis, but several scientific and technological challenges remain to fully electrify both the chemical industry and transportation. In both of these areas, new electrochemical materials will be critical, but their development currently relies heavily on human-time-intensive experimental trial and error and computationally expensive first-principles, meso-scale and continuum simulations. To accelerate this process, our team has developed the AutoMat platform. AutoMat can accelerate development of new electrochemical materials along two avenues: first, automated input generation and management of simulations at multiple lengthscales as well as “handoff” of outputs from one lengthscale as inputs to the next; and second, replacement of the most computationally intensive simulation processes with machine-learned surrogate models. The crux of our team’s effort was not “reinventing the wheel” by developing entirely new techniques, but rather building a “superhighway” that allows existing state-of-the-art techniques to run faster and more smoothly than before. AutoMat can utilize tools spanning from first-principles quantum chemistry computations to automated robotic experimentation, and is driven by design space search techniques to reduce the number of iterations through the full simulation loop by rapidly targeting promising regions of design spaces such as single-atom alloy catalysts or blends of liquid electrolytes.

25 ENERGY STORAGE↗

pH Regulates Ion Dynamics in Carboxylated Mixed Conductors

Coupled ionic and electronic transport underpins processes as diverse as electrochemical energy conversion, biological signaling, and soft adaptive electronics. Yet, how chemical environments such as pH modulate this coupling at the molecular scale remains poorly understood. Here, we show that the protonation state of carboxylated polythiophenes provides precise chemical control over ion dynamics, doping efficiency, solvent uptake, and mechanical response. Using a suite of multimodal operando techniques, supported by simulations, we reveal that pH dictates the balance of cation/anion uptake during electrochemical doping. Mapping across pH uncovers a quasi-nonswelling regime (≈pH 3–3.5) where charge compensation proceeds with minimal volumetric change yet pronounced stiffening. These findings establish molecular acidity as a general strategy to program ionic preference and mechanical stability, offering design principles for pH-responsive mixed conductors and soft electronic materials that couple ionic, electronic, and mechanical functionality.

PH↗

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↗

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↗

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↗

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↗