Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational methods Electrochemistry”

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

An immersed interface method for microstructure-scale electrochemical battery models: numerical formulation and performance portable implementation

We present the numerical formulation, verification, and performance portable implementation of an immersed interface method for microstructure scale electrochemical modeling of batteries. The innovation in this approach is the resolution of chemical species and electrostatic potential discontinuities at active interfaces without the use of interface conforming unstructured grids. A unified formulation on Cartesian grids for all domains (electrodes and electrolyte) is used with interfacial flux conditions applied using volume fraction or “color” function gradients. We have developed one dimensional and two dimensional test cases with analytic solutions for electrochemical modeling using which we verified the consistency and accuracy of our scheme. Our solver is also validated against solutions from a macroscale model and an unstructured multi-subdomain solver for a full lithium ion cell. We then demonstrated the utility of our solver on an image-based complex battery electrode microstructure at high charging rate. Our technique also exhibits good scalability on distributed memory architectures using central processing units (CPU), with problem sizes up to 1.8 billion degrees of freedom and with 5400 ranks. Initial performance studies of our open-source performance portable solver showed about 70 times speed up using a graphics processing unit (GPU) compared to single compute core for a problem with 4 million cells.

25 ENERGY STORAGE↗

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↗

Autonomous Electrochemistry Platform with Real-Time Normality Testing of Voltammetry Measurements Using ML

Electrochemistry workflows utilize various instruments and computing systems to execute workflows consisting of electrocatalyst synthesis, testing and evaluation tasks. The heterogeneity of the software and hardware of these ecosystems makes it challenging to orchestrate a complete workflow from production to characterization by automating its tasks. We propose an autonomous electrochemistry computing platform for a multi-site ecosystem that provides the services for remote experiment steering, real-time measurement transfer, and AI/ML-driven analytics. We describe the integration of a mobile robot and synthesis workstation into the ecosystem by developing custom hub-networks and software modules to support remote operations over the ecosystem’s wireless and wired networks. We describe a workflow task for generating I-V voltammetry measurements using a potentiostat, and a machine learning framework to ensure their normality by detecting abnormal conditions such as disconnected electrodes. We study a number of machine learning methods for the underlying detection problem, including smooth, non-smooth, structural and statistical methods, and their fusers. We present experimental results to illustrate the effectiveness of this platform, and also validate the proposed ML method by deriving its rigorous generalization equations.

Alnajjar, Anees↗

An Electrochemical Strategy to Synthesize Disilanes and Oligosilanes from Chlorosilanes**

Silanes are important compounds in industrial and synthetic chemistry. Here, we develop a general approach for the synthesis of disilanes as well as linear and cyclic oligosilanes via the reductive activation of readily available chlorosilanes. The efficient and selective generation of silyl anion intermediates, which are arduous to achieve by other means, allows for the synthesis of various novel oligosilanes by heterocoupling. In particular, this work presents a modular synthesis for a variety of functionalized cyclosilanes, which may give rise to materials with distinct properties from linear silanes but remain challenging synthetic targets. In comparison to the traditional Wurtz coupling, our method features milder conditions and improved chemoselectivity, broadening the functional groups that are compatible in oligosilane preparation. Furthermore, computational studies support a mechanism whereby differential activation of sterically and electronically distinct chlorosilanes are achieved in an electrochemically driven radical-polar crossover mechanism.

36 MATERIALS SCIENCE↗

Mathematical Model of a Regenerative Fuel Cell for System Optimization

This thesis developed a system-level optimization model of a regenerative fuel cell (RFC) system for long-duration, off-world energy storage applications. Prior RFC design studies have typically been limited to reduced parameter sets and simplified constraints due to computational limitations relative to the number of relevant degrees of freedom. As a result, important nonlinear interactions between subsystems have not been fully captured. This work began to address that gap by developing a higher-fidelity, nonlinear optimization framework that incorporates a broader set of design variables and coupled constraints, enabling a multidimensional model that captures the coupled behavior of RFC subsystems and demonstrates the feasibility of applying optimization to such systems. An expanded system-level optimization approach was established that captures interactions between electrochemical performance, structural requirements, and storage design. This enabled a more comprehensive evaluation of trade-offs than conventional formulations. The model integrates four coupled subsystems: a fuel cell, an electrolyzer, reactant gas, and high-pressure storage tanks, and was formulated to accommodate a wide range of mission parameters, including operational time and required output power. It incorporates constraints on available solar array power, reactant mass balance between production and consumption, and pressure-dependent storage requirements. To enable reliable convergence, the optimization problem was reformulated to reduce dimensionality and improve numerical stability, with subsystem models organized for efficient evaluation. Problem dimensionality was reduced by consolidating lower-level design variables into higher-level representative quantities, and subsystem behavior was evaluated within the optimization loop. A multi-start initialization strategy was employed to mitigate sensitivity to local minima and improve solution quality, while nonlinear relationships were solved using robust numerical methods. The results showed that convergence was achieved across a range of required output power values. Specific energy reached a maximum at a critical mission power level, where the electrolyzer power matched the available solar input and operated near its voltage and current density limits. Beyond this point, further increases in required power resulted in less mass-efficient operation, increasing total system mass and reducing overall performance. The developed model represents an advancement in RFC system-level optimization by enabling analysis of a broader and more tightly coupled design space than previous considerations. While convergence behavior and computational cost remain challenges, the methods introduced improve solvability and allow inclusion of additional design variables with minimal loss of physical fidelity. However, the numerical results should not be interpreted as definitive design recommendations, as the model includes simplifying assumptions and omits several higher-order effects. Future work should extend this framework by incorporating additional subsystems and loss mechanisms, such as thermal management, parasitic power consumption, and reactant losses, to improve fidelity and ensure more representative design conclusions.

Electrochemistry↗

Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

Implicit solvent models are a computationally efficient method of representing solid/liquid interfaces prevalent in electrocatalysis, energy storage, and materials science. However, electronic structure changes induced at the metallic surface by the dielectric continuum are not fully understood. To address this, we perform DFT calculations for the Pt(111)/water interface, in order to compare Poisson–Boltzmann continuum solvation methods with ab initio molecular dynamics (AIMD) simulations of explicit solvent. We show that the implicit solvent cavity can be parametrized in terms of the electric dipole moment change at the equilibrated explicit Pt/water interface to obtain the potential of zero charge (PZC). We also compare the accuracy of aqueous enthalpies of adsorption of phenol on Pt(111) using geometry and charge density based dielectric cavitation methods. The ability to parametrize the cavity according to individual atoms, as afforded in the geometry based approach, is key to obtaining accurate enthalpy changes of adsorption under aqueous conditions. Additionally, we show that the electronic structure changes induced by explicit solvent and our proposed implicit solvent parametrization scheme yield comparable density difference profiles and d-band projected density of states. We therefore demonstrate the capability of implicit solvent approaches to capture both the energetics of adsorption processes and the main electronic effects of aqueous solvent on the metallic surface. Therefore, this work provides a scheme for computationally efficient simulations of interfacial processes for applications in areas such as heterogeneous catalysis and electrochemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Physics-Based Modeling Framework for Prognostic Studies

Prognostics and Health Management (PHM) methodologies have emerged as one of the key enablers for achieving efficient system level maintenance as part of a busy operations schedule, and lowering overall life cycle costs. PHM is also emerging as a high-priority issue in critical applications, where the focus is on conducting fundamental research in the field of integrated systems health management. The term diagnostics relates to the ability to detect and isolate faults or failures in a system. Prognostics on the other hand is the process of predicting health condition and remaining useful life based on current state, previous conditions and future operating conditions. PHM methods combine sensing, data collection, interpretation of environmental, operational, and performance related parameters to indicate systems health under its actual application conditions. The development of prognostics methodologies for the electronics field has become more important as more electrical systems are being used to replace traditional systems in several applications in the aeronautics, maritime, and automotive fields. The development of prognostics methods for electronics presents several challenges due to the great variety of components used in a system, a continuous development of new electronics technologies, and a general lack of understanding of how electronics fail. Similarly with electric unmanned aerial vehicles, electrichybrid cars, and commercial passenger aircraft, we are witnessing a drastic increase in the usage of batteries to power vehicles. However, for battery-powered vehicles to operate at maximum efficiency and reliability, it becomes crucial to both monitor battery health and performance and to predict end of discharge (EOD) and end of useful life (EOL) events. We develop an electrochemistry-based model of Li-ion batteries that capture the significant electrochemical processes, are computationally efficient, capture the effects of aging, and are of suitable accuracy for reliable EOD prediction in a variety of usage profiles.

Li-ion Batteries↗

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis and characterization of the dinuclear cobalt(III) complex: [(C 2 F 5 ) 3 Co(μ-F)] 2 2–

Here, reaction of the versatile tris(perfluoroethyl) cobalt(III) precursor [fac-(MeCN) 3 Co(C 2 F 5 ) 3 ] with [NMe 4 ]F and [PPh 4 ]Cl in THF formed the unexpected cobalt(III) bridging fluoride dimer ([(C 2 F 5 ) 3 Co(μ-F)] 2 2– . The new fluoro-organometallic cobalt(III) complex was characterized by NMR and UV-vis spectroscopies, X-ray crystallography, cyclic voltammetry, and by computational methods. In the strongly coordinating solvent MeCN, ([(C 2 F 5 ) 3 Co(μ-F)] 2 2– exhibits dynamic processes on the NMR timescale which are consistent with solvent coordination. However, in the weakly coordinating solvent CH 2 Cl 2 , a more static structure is observed suggesting that the dimer retains its structure in solution state. An electrochemical analysis of ([(C 2 F 5 ) 3 Co(μ-F)] 2 2– was performed, and the data are compared to previously reported cobalt(III) perfluoroethyl complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Charting C–C coupling pathways in electrochemical CO 2 reduction on Cu(111) using embedded correlated wavefunction theory

The electrochemical CO 2 reduction reaction (CO 2 RR) powered by excess zero-carbon-emission electricity to produce especially multicarbon (C 2+ ) products could contribute to a carbon-neutral to carbon-negative economy. Foundational to the rational design of efficient, selective CO 2 RR electrocatalysts is mechanistic analysis of the best metal catalyst thus far identified, namely, copper (Cu), via quantum mechanical computations to complement experiments. Here, we apply embedded correlated wavefunction (ECW) theory, which regionally corrects the electron exchange-correlation error in density functional theory (DFT) approximations, to examine multiple C–C coupling steps involving adsorbed CO (*CO) and its hydrogenated derivatives on the most ubiquitous facet, Cu(111). We predict that two adsorbed hydrogenated CO species, either *COH or *CHO, are necessary precursors for C–C bond formation. The three kinetically feasible pathways involving these species yield all three possible products: *COH–CHO, *COH–*COH, and *OCH–*OCH. The most kinetically favorable path forms *COH–CHO. In contrast, standard DFT approximations arrive at qualitatively different conclusions, namely, that only *CO and *COH will prevail on the surface and their C–C coupling paths produce only *COH–*COH and *CO–*CO, with a preference for the first product. This work demonstrates the importance of applying qualitatively and quantitatively accurate quantum mechanical method to simulate electrochemistry in order ultimately to shed light on ways to enhance selectivity toward C 2+ product formation via CO 2 RR electrocatalysts.

30 DIRECT ENERGY CONVERSION↗

High entropy powering green energy: hydrogen, batteries, electronics, and catalysis

A reformation in energy is underway to replace fossil fuels with renewable sources, driven by the development of new, robust, and multi-functional materials. High-entropy materials (HEMs) have emerged as promising candidates for various green energy applications, having unusual chemistries that give rise to remarkable functionalities. This review examines recent innovations in HEMs, focusing on hydrogen generation/storage, fuel cells, batteries, semiconductors/electronics, and catalysis—where HEMs have demonstrated the ability to outperform state-of-the-art materials. We present new master plots that illustrate the superior performance of HEMs compared to conventional systems for hydrogen generation/storage and heat-to-electricity conversion. We highlight the role of computational methods, such as density functional theory and machine learning, in accelerating the discovery and optimization of HEMs. The review also presents current challenges and proposes future directions for the field. We emphasize the need for continued integration of modeling, data, and experiments to investigate and leverage the underlying mechanisms of the HEMs that are powering progress in sustainable energy.

batteries↗

Quantifying Volume Change in Porous Electrodes via the Multi-Species, Multi-Reaction Model

Automotive manufacturers are working to improve individual cell and overall pack design by increasing their performance, durability, and range, while reducing cost; and active material volume change is one of the more complex aspects that needs to be considered during this process. As the time from initial design to manufacture of electric vehicles is decreased, design work that used to rely solely on testing needs to be supplemented or replaced by virtual methods. As electrochemical engineers drive battery and system design using model-based methods, the need for coupled electrochemical/mechanical models that take into account the active material change utilizing physics based or semi-empirical approaches is necessary. In this study, we illustrated the applicability of a mechano-electrochemical coupled modeling method considering the multi-species, multi-reaction model as popularized by Verbrugge and Baker. To do this, validation tests were conducted using a computer-controlled press apparatus that can control the press displacement and press force with precision. The coupled MSMR volume change model was developed and its applicability to graphite and NMC cells was illustrated. The increased accuracy of the model considering the coupled MSMR volume change approach shows in the importance of accounting for individual gallery volume change behavior on cell level predictions.

25 ENERGY STORAGE↗

Gerischer Electrochemistry Today

Semiconductor photoelectrochemistry is a dynamic and interdisciplinary field at the forefront of research in solar fuels, energy conversion, and catalysis. Here, this Perspective captures the collective insights from the second Gerischer Electrochemistry Today Symposium, held at Colorado State University in Fort Collins, CO, in August 2024, which convened leading researchers, early-career scientists, and industry partners to define the critical next steps for the field. Through interactive sessions, technical talks, panel discussions, and training initiatives─including a Semiconductor Electrochemistry Bootcamp─the symposium emphasized three pillars of advancement: (i) facilitating the exchange of new ideas in semiconductor electrochemistry and charge separation; (ii) fostering the development of future researchers, research topics, and participation in the semiconductor workforce; and (iii) building community. This Energy Focus distills key themes from the meeting and identifies major knowledge gaps in the following areas: mechanisms of charge separation and recombination, role of defects and disorder, dynamic and operando characterization methods, interfacial chemistry and surface passivation, theoretical and modeling limitations, and standardization and benchmarking. The inclusive and collaborative structure of the symposium enabled the generation of this comprehensive report that will serve as a roadmap for fundamental and applied research in the rapidly evolving field of semiconductor electrochemistry over the next decade.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Simplified Universal Equations for Ionic Conductivity and Transference Number

Nernst-Einstein equation can provide a reasonable estimate of the ionic conductivity of dilute solutions. For concentrated solutions, alternate methods such as Green–Kubo relations and Einstein relations are more suitable to account for ion-ion interactions. Such computations can be expensive for multicomponent systems. Simplified mathematical expressions like the Nernst-Einstein equation do not exist for concentrated multicomponent mixtures. Newman’s treatment of multicomponent concentrated solutions yields a conductivity relation in terms of species concentration and Onsager phenomenological coefficients. However, the estimation of these phenomenological coefficients is not straightforward. Here, mathematical formulations that relate the phenomenological coefficients with the friction coefficients are developed, leading to simplified, ready-to-use expressions of conductivity and transference numbers that can be used for a wide range of ionic mixtures. This approach involves spectral decomposition of the matrix of Onsager phenomenological coefficients. The general analytical expressions for conductivity and transference number are simplified for binary electrolytes, and numerical solutions are provided for ternary and quaternary mixtures with ion dissociation.

Electrochemistry↗

Electro-carburization of mild steel through two-step reaction in molten-salt electrolyte

The rapid carburization of mild steel is demonstrated via the electrochemical deposition of solid carbon, followed by thermochemical iron-carbon (Fe-C) reaction using a molten lithium-potassium carbonate (LiKCO 3 ) salt electrolyte and carbon dioxide as the carbon source. Depth-dependent microstructural and compositional changes in the steel are shown to be controlled by the applied cell voltage and reaction temperatures. Reaction temperatures above the Fe-C eutectoid temperature promote faster carburization due to carbon dissolution into austenite. The carbon deposition rate is regulated by the applied cell voltage: At - 2.4 V, carbon deposition occurs much faster than its diffusion into the steel, resulting in the accumulation of a carbon-rich layer exceeding 1 mm, which does not diffuse into the steel. In contrast, at - 1.8 V, all deposited carbon effectively diffuses into the steel, leaving a minimal external deposit. In both cases, carbon ingress forms a carbide-rich surface layer and precipitates additional carbides along the grain boundaries within the steel electrode. This method offers a potential approach to increase the carbon content of mild or low-carbon steel feedstocks, such as scrap or direct reduced iron, for industrial applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Topology optimization for the design of porous electrodes

Porous electrodes are an integral part of many electrochemical devices since they have high porosity to maximize electrochemical transport and high surface area to maximize activity. Traditional porous electrode materials are typically homogeneous, stochastic collections of small-scale particles and offer few opportunities to engineer higher performance. Fortunately, recent breakthroughs in advanced and additive manufacturing are yielding new methods to structure and pattern porous electrodes across length scales. These architected electrodes are emerging as a promising new technology to continue to drive improvement; however, it is still unclear which structures to employ and few tools are available to guide their design. Here we address this gap by applying topology optimization to the design of porous electrodes. We demonstrate our framework on two applications: a porous electrode driving a steady Faradaic reaction and a transiently operated electrode in a supercapacitor. We present computationally designed electrodes that minimize energy losses in a half-cell. For low-conductivity materials, the optimization algorithm creates electrode designs with a hierarchy of length scales. Further, the designed electrodes are found to outperform undesigned, homogeneous electrodes. Finally, we present three-dimensional porous electrode designs. We thus establish a topology optimization framework for designing porous electrodes.

25 ENERGY STORAGE↗

Structural and electrochemical evolution of alloy interfacial layers in anode-free solid-state batteries

“Anode-free” solid-state batteries feature high energy density since there is no anode active material. Although the beneficial effects of interfacial layers at the anode-solid electrolyte interface have been demonstrated, the mechanisms through which they influence lithium plating/stripping are unclear. Here, in this study, we reveal the evolution of 100-nm silver and gold interfacial layers during lithium plating/stripping using electrochemical methods, electron microscopy, X-ray microcomputed tomography, and modeling. The alloy layers improve Coulombic efficiency and resistance to short circuiting, even though the alloys form solute regions or particulates that detach from the current collector during plating. In situ electrochemical impedance spectroscopy shows that the alloys return to the interface and mitigate contact loss at the end of stripping, avoiding a critical vulnerability of anode-free cells. Contact retention is driven by even Li thickness that promotes spatially uniform stripping, as well as local alloy delithiation in response to current concentrations that homogenizes current and diminishes voiding.

25 ENERGY STORAGE↗