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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Mesoscale Effects of Composition and Calendering in Lithium-Ion Battery Composite Electrodes

Macrohomogeneous battery models are widely used to predict battery performance, necessarily relying on effective electrode properties, such as specific surface area, tortuosity, and electrical conductivity. While these properties are typically estimated using ideal effective medium theories, in practice they exhibit highly non-ideal behaviors arising from their complex mesostructures. In this paper, we computationally reconstruct electrodes from X-ray computed tomography of 16 nickel–manganese–cobalt-oxide electrodes, manufactured using various material recipes and calendering pressures. Due to imaging limitations, a synthetic conductive binder domain (CBD) consisting of binder and conductive carbon is added to the reconstructions using a binder bridge algorithm. Reconstructed particle surface areas are significantly smaller than standard approximations predicted, as the majority of the particle surface area is covered by CBD, affecting electrochemical reaction availability. Finite element effective property simulations are performed on 320 large electrode subdomains to analyze trends and heterogeneity across the electrodes. Significant anisotropy of up to 27% in tortuosity and 47% in effective conductivity is observed. Electrical conductivity increases up to 7.5× with particle lithiation. We compare the results to traditional Bruggeman approximations and offer improved alternatives for use in cell-scale modeling, with Bruggeman exponents ranging from 1.62 to 1.72 rather than the theoretical value of 1.5. We also conclude that the CBD phase alone, rather than the entire solid phase, should be used to estimate effective electronic conductivity. This study provides insight into mesoscale transport phenomena and results in improved effective property approximations founded on realistic, image-based morphologies.

25 ENERGY STORAGE↗
DOE OSTI2023

Understanding the Surprising Ionic Conductivity Maximum in Zn(TFSI) 2 Water/Acetonitrile Mixture Electrolytes

Aqueous electrolytes composed of 0.1 M zinc bis-(trifluoromethyl-sulfonyl)-imide (Zn-(TFSI) 2 ) and acetonitrile (ACN) were studied using combined experimental and simulation techniques. The electrolyte was found to be electrochemically stable when the ACN V% is higher than 74.4. In addition, it was found that the ionic conductivity of the mixed solvent electrolytes changes as a function of ACN composition, and a maximum was observed at 91.7 V% of ACN although the salt concentration is the same. This behavior was qualitatively reproduced by molecular dynamics (MD) simulations. Detailed analyses based on experiments and MD simulations show that at high ACN composition the water network existing in the high water composition solutions breaks. As a result, the screening effect of the solvent weakens and the correlation among ions increases, which causes a decrease in ionic conductivity at high ACN V%. Furthermore, this study provides a fundamental understanding of this complex mixed solvent electrolyte system.

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

Electrification and decarbonization of spent Li-ion batteries purification by using an electrochemical membrane reactor

The expanding electric vehicle market brings with it exponential growth in the use of lithium (Li)-ion batteries (LIB) for which a wave of spent LIB is expected to come within the next 5 to 10 years. Due to the economic and strategic value imbedded within the metals contained in LIB, different recycling technologies, including hydrometallurgy, pyrometallurgy and direct recycling, are under development. Being different from previous hydrometallurgical methods, which may have high chemical consumption and negative environmental impact, an electrochemical membrane reactor is designed and validated for the first time, to electrify and decarbonize the impurity removal process. This reactor electroplates copper (Cu) and electrochemically precipitates aluminum (Al) and iron (Fe) from simulated spent LIB leachates, by consuming only air, water, and electricity, and the impurities are reduced to <1 ppm. The purified leachate maintains 99.5 % of the nickel (Ni), 95.4 % of the cobalt (Co) and 99.14 % of manganese (Mn) from the original leachate solution, and then can be directly applied for cathode precursor synthesis. Additionally, the purification process doesn’t introduce extra impurity, and the reactor restoration process generates valuable by-product hydro sulfate (H2SO4). This electrochemical process can reduce the cost, because of the much less chemical consumption and the valuable by-product generation, and mitigates the waste emissions, because of no extra impurity introduced and no greenhouse gas (GHG) produced. In conclusion, the chemical precipitation method uses significant amount of NaOH, which induced GHG emission during the manufacturing process.

25 ENERGY STORAGE↗
DOE OSTI2022

Numerical Simulation and Modeling of Hydrogen Gas Evolution on Planar and Microwire Array Electrodes

The impact of gas evolution on the electrochemical characteristics of planar electrodes and microwire array electrodes has been analyzed using modeling and simulation. The impacts can mainly be broken into three phenomena: a) a shift in the local reversible hydrogen electrode potential; b) hyperpolarization; and c) an increase in the solution resistance of the electrolyte. The local reversible hydrogen electrode potential shift was found to play the most important role, constituting >40% of the total potential drop between the cathode and reference electrode, following correction for cell resistance. Compared to planar electrodes, a microwire array structure reduces the impact of bubbles on the solution conductance, but the shift in the local reversible hydrogen electrode potential varies with distance from the actual electrode surface.

25 ENERGY STORAGE↗
DOE OSTI2024

CO 2 Conversion to Butene via a Tandem Photovoltaic–Electrochemical/Photothermocatalytic Process: A Co-design Approach to Coupled Microenvironments

Here, we developed a tandem, unassisted, solar-driven electrochemical and photothermocatalytic process for the single-pass conversion of CO 2 to butene using only simulated solar irradiation as the energetic input. The two-step process involves electrochemical CO 2 reduction (CO 2 R) to ethylene followed by ethylene dimerization to butene. We assessed two unassisted electrochemical setups to concentrate ethylene in the CO 2 R reactor, achieving concentrations up to 5.4 vol.% with 1.8% average solar-to-ethylene conversion and 5.6% average CO 2 -to-ethylene single-pass conversion under 1-sun illumination. When passed through the photothermocatalytic ethylene oligomerization reactor, we generated 600 ppm of butene under 3-sun illumination. Through analysis of this process, we identified that the presence of H 2 , CO, and H 2 O leads to rapid deactivation of the Ni-based ethylene oligomerization catalyst.

14 SOLAR ENERGY↗
DOE OSTI2025

Bubble Transport through a Porous Lattice with an Applied Inlet Flow

Within gas-evolving electrochemical systems, bubbles negatively impact performance by covering electrode active sites for reactions, blocking electric field lines, and obstructing liquid electrolyte flow causing pressure buildup. Recent additive manufacturing advances have enabled tuned porous electrode microstructures to be created, but producing systems that maximize electrochemical throughput and minimize bubble impact remains challenging. Thus, improved physical understanding of and modeling capabilities for bubble behavior are critical to improve electrolyzer design. To address this need, this study examines rising stage bubbles within a lattice with an applied liquid flow—an underexplored regime that strongly influences an electrochemical bubble’s fate. Notably, theoretical predictions and resolved bubble simulations are complemented by experiments from a 3D-printed visualization cell that matches the simulation geometry. The minimum threshold flow rate to achieve bubble breakthrough is found to be larger for higher porosities and for smaller bubbles. Different-sized bubbles decrease expected electrochemical performance in different ways; smaller bubbles tend to stay stuck but cover less solid surface, while larger bubbles more readily break through but cover more surface while in the lattice. The bubble trajectory, deformation, and contact area provide insight into these different behaviors. These findings provide design guidelines toward creating more effective electrolyzers.

Guo, Jack [Lawrence Livermore National Laboratory ↗
NASA NTRS1988

Nickel cadmium battery performance simulation - A means toward optimum design in satellite applications

It is noted that electrochemically based cell models can advance state-of-the-art battery prediction by allowing variations in the parameters that constitute wear-out and predict end-of-life performance. This modeling method is used as the foundation for the simulation of Ni-Cd battery performance. This simulation tool consists of a comprehensive computer model to be developed which would predict the voltage behavior of a battery under mission conditions.

Clark, K.↗
NASA NTRS1981

Computer simulation and experimental validation of isoelectric focusing in ampholine-free systems

An abbreviated version of a mathematical model of the steady state in isoelectric focusing is presented. Details of the mathematical transformations leading to a model suitable for computer implementation and numerical solution are given in the Appendix. The model describes the structure of concentration, pH, conductivity and potential gradients arising from the focusing of electrochemically defined ampholytes. Some of the results of computer simulations of two and three component systems are of particular interest to the experimentalist and are presented together with experimental validation of this model.

Bier, M.↗
DOE OSTI2020

Transport Processes in a Li-ion Cell during an Internal Short-Circuit

Internal short-circuit in a lithium-ion cell causes an abrupt increase in cell temperature and triggers subsequent thermal runaway. In this work, we present a detailed electrochemical-thermal model to investigate the physical behavior during an internal short-circuit. Simulations at wide range of heat transfer coefficients and short-circuit resistances are conducted to illustrate electrochemical and thermal behavior under a wide range of conditions. The Joule heating at the shorted region promotes electrochemical reactions nearby, causing in-plane non-uniformity of electrolyte and active material transport. Furthermore, it is found that diffusion in solid active materials plays a significant role at very high shorting currents (~20 C), as electrochemical reactions rate are being controlled progressively more by availability of Li + at the interface, due to limitations in diffusion through the active material with increasing discharge rates. This diffusion limitation causes a drop in available energy, and subsequently a decrease in cell equilibrium potential and the heat generation rate at the location of the short. On the other hand, rapid depletion of lithium concentration in the electrolyte and accumulation on the electrode surface results in highly non-uniform transport properties resulting in higher heat generation rates. Hence, the heating regime shifts from "local heating" to "global heating". Based on the findings, important design parameters for battery safety are discussed.

25 ENERGY STORAGE↗
DOE OSTI2022

Review—Development of Highly Active and Stable Catalyst Supports and Platinum–Free Catalysts for PEM Fuel Cell

Metal-free, metal-containing, and template-assisted metal-containing nitrogen-modified carbon-based catalysts carbon composite supports, and highly active and electrochemically stable hybrid cathode catalysts for oxygen reduction reaction are reviewed in this manuscript. Here, novel procedures are developed for the synthesis of (i) highly catalytically active PGM-free catalyst and stable carbon composite catalyst supports. The carbon composite catalyst supports are engineered with optimized BET surface area and pore size distribution and with a well-defined kinetic and mass transfer region during the reaction and their performance is dicussed in detail in this review. The surface activation results in increasing carbon graphitization and inclusion of non-metallic active sites on the support surface. The USC catalysts exhibit the initial performance of 0.91 V and a maximum power density of 177 mW cm –2 , with well-defined kinetic and mass-transfer regions and ~2.5% H 2 O 2 production. Activated carbon composite support (ACCS) is modified to optimize its kinetic activity and its electrochemical stability shows excellent thermal stability and support stability under simulated start-up/shut-down operating conditions. The stabilities of various supports developed in this study are compared with those of a commercial Pt/C catalyst. The active sites for the ORR identified with electrochemical and physicochemical methods are pyridinic-N and quaternary-N.

25 ENERGY STORAGE↗
DOE OSTI2020

In-situ/operando X-ray absorption spectroscopic investigation of the electrode/electrolyte interface on the molecular scale

A method for performing X-ray absorption spectroscopy (XAS) in the soft X-ray region was developed to investigate the surface-sensitive electron signal of electrode/electrolyte interfaces in common multivalent based organic solvents under in situ/operando conditions. Our approach enables us to probe the molecular-scale structure of electrode interfaces by measuring total electron yield and is suitable for redox systems exhibiting low intrinsic electrochemical current. In-situ F K-edge XAS measurements in a 0.5 M Magnesium bis(trifluoromethanesulfonimide)/2-Methyltetrahydrofuran (Mg(TFSI) 2 /2-MeTHF) electrolyte were carried out to determine the evolution of interfacial species during the electrochemical charging/discharging process. Time-dependent density-functional theory (TD-DFT) simulation indicate that the F K-edge evolution is the result of interfacial chemical environment change driven by electrochemical potential. In addition, we performed the “operando XAS” which runs the CV and collects spectra simultaneously. Furthermore, by using this method, the non-equilibrium state of condensed matter interfaces and interfacial dynamical transient process can be revealed.

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

Removal of As(V) from wastewaters using magnetic iron oxides formed by zero-valent iron electrocoagulation

Electrocoagulation of zero-valent iron has been widely applied to the removal of dissolved arsenic, but the solid-liquid separation of arsenic-containing precipitates remains technically challenging. In this work, zero-valent iron was electrochemically oxidized to magnetic iron oxides for the removal of As(V) from simulated and actual mining wastewaters. The results indicated that lepidocrocite was formed when zero-valent iron was oxidized by dissolved oxygen, but ferrihydrite and green rust were first formed and then transformed to magnetic iron oxides (mainly magnetite and maghemite) in the electrochemical oxidation from 0 to 0.9 V (vs. SCE), which facilitates the adsorption of As(V) and subsequent solid-liquid separation under a magnetic field. In simulated As(V)-containing solution with initial pH 7.0, zero-valent iron was electrochemically oxidized to magnetite and maghemite at 0.6 V (vs. SCE) for 2 h. The As(V) concentration first decreased from 5127.5 to 26.8 μg L -1 with a removal ratio of 99.5%. In actual mining wastewaters, zero-valent iron was electrochemically oxidized to maghemite at 0.6 V (vs. SCE) for 24 h, and the As(V) concentration decreased from 5486.4 to 3.6 μg L -1 with a removal ratio of 99.9%. The removal ratio of As(V) increased slightly with increasing potential, and increased first and then decreased with increasing initial pH. Compared with that of SO 4 2- and NO 3 - , the presence of Cl - significantly enhanced the removal of As(V). This work provides a highly efficient, facile and low-cost technique for the treatment of arsenic-containing wastewaters.

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

A Segregated Approach for Modeling the Electrochemistry in the 3-D Microstructure of Li-Ion Batteries and Its Acceleration Using Block Preconditioners

Abstract Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the microstructure scale. A fast-iterative solver with an appropriate preconditioner is then required to simulate large representative volume in a reasonable time. In this work, a finite element electrochemical model is developed to resolve the concentration and potential within the electrode active materials and the electrolyte domains at the microstructure scale, with an emphasis on numerical stability and scaling performances. The block Gauss-Seidel (BGS) numerical method is implemented because the system of equations within the electrodes is coupled only through the nonlinear Butler–Volmer equation, which governs the electrochemical reaction at the interface between the domains. The best solution strategy found in this work consists of splitting the system into two blocks—one for the concentration and one for the potential field—and then performing block generalized minimal residual preconditioned with algebraic multigrid, using the FEniCS and the Portable, Extensible Toolkit for Scientific Computation libraries. Significant improvements in terms of time to solution (six times faster) and memory usage (halving) are achieved compared with the MUltifrontal Massively Parallel sparse direct Solver. Additionally, BGS experiences decent strong parallel scaling within the electrode domains. Last, the system of equations is modified to specifically address numerical instability induced by electrolyte depletion, which is particularly valuable for simulating fast-charge scenarios relevant for automotive application.

25 ENERGY STORAGE↗
DOE OSTI2025

Cooperative Effects Associated with High Electrolyte Concentrations in Driving the Conversion of CO2 to C2H4 on Copper

Compared to a conventional electrolyte concentration of 1 M HCOOK, the use of a highly concentrated 7.1 M HCOOK electrolyte increases the Faradaic efficiency (FE) ratio of C2H4/CO from 2.2 +- 0.3 to 18.3 +- 4.8 at -1.08 V vs. reversible hydrogen electrode (RHE) on a Cu gas-diffusion electrode. Based on electrochemical analysis and ab initio molecular dynamics (AIMD) simulation, the identity and concentration of the cation and anion play more important roles in controlling the CO2R reaction pathway than the bulk CO2 solubility and the bulk pH of electrolytes. In situ attenuated reflectance surface enhanced infrared absorption spectroscopy (ATR-SEIRAS) suggests that, unlike 1 M HCOOK, the *CO-bridge-binding mode on Cu is dominant in 7.1 M HCOOK electrolyte, which potentially results in less CO release and higher yield of C2H4. This study demonstrates that although we can tailor the electrolyte composition to shift product selectivity, the factors that control the product selectivity are numerous and cannot be distilled down into one correlated property-reactivity relationship.

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

Spontaneous sodium ion storage behaviors of reduced graphene oxide anodes exceeding 100% Coulombic efficiency by modulated ion solvation

Rechargeable batteries are essential energy storage devices that power portable devices and electrical vehicles throughout the world. In general, it is thought that the electrochemical performance of rechargeable batteries is mostly determined by the electrodes within them and that the electrolyte plays a relatively passive role. However, ion transport and storage can be greatly influenced by the electrolyte solution structure, specifically, ion solvation within the bulk and ion desolvation across the electrode/electrolyte interfaces. Herein, we studied the role of the electrolyte as an active component of electrochemical energy storage devices. We found that with an appropriate electrolyte formulation, ion storage in disordered carbonaceous anode materials can occur spontaneously without externally supplied electrical energy. Reduced graphene oxide (RGO) in an ether-based electrolyte demonstrates ‘spontaneous' ion storage behaviors of adsorbing and inserting the solvated ions utilizing facilitated permeability and wettability of RGO, which results in Coulombic efficiency of ~145% due to additional charging capacity of ~180 mAh g -1 during electrochemical processes. The unexpected spontaneous ion storage behavior was extensively investigated using a combination of electrochemical analyses and diagnostics, advanced characterizations, and computational simulation. In conclusion, we believe the spontaneous ion storage behavior offers a new way to further improve the energy efficiency of practical rechargeable batteries.

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

C–C Bond Formation during Electrochemical CO 2 Reduction on Pristine Cu(100) Unlikely to Involve Adsorbed CO at Any Potential

Formation of hydrocarbons containing two or more carbon atoms (C 2+ ) during heterogeneous electrochemical CO and CO 2 reduction (ECOR and ECO 2 R) only occurs, among pure metals, on Cu electrodes. Moreover, the activity and selectivity is facet dependent, with Cu(100) generally preferentially forming ethylene over methane. Previously, we found via quantum-mechanics-based modeling that, unlike standard density functional theory, more accurate correlated wavefunction methods predict that non-electroactive coupling pathways involving two adsorbed COs (*CO) or a *CO and a *COH to form C–C bonds on Cu(100) are kinetically inhibited, with the former also thermodynamically unfavorable. Here, we extend that embedded complete active space second order perturbation theory (ECASPT2) study, further showing that electrochemical coupling of two *COs to form an anionic dimer [OC*–*CO] (1+δ)– , followed by protonation to form [OC*–*COH] δ− , is not kinetically competitive with the reduction of *CO to *COH at relevant ECO/CO 2 R potentials. Our simulations therefore suggest that the ability of Cu(100) to electrochemically synthesize C 2+ molecules from CO and CO 2 is unlikely to be via *CO, at least on pristine Cu(100). Instead, hydrogenated CO species (*COH, *CH x OH, or *CH x ) are most likely to be the key intermediates in C–C bond formation.

Martirez, John Mark P. [Princeton Plasma Physics L↗
DOE OSTI2024

Rational Design of High-Performing Electrodes in Energy Storage Devices

Our approach in this project is to control conductivities of different ions by modulating the pore sizes in MOFs as sulfur hosts in cathodes. To validate our hypothesis that controlling the pore size of MOFs can effectively prevent LiPSs shuttling while allowing sufficient Li ion transport, which will ensure electrochemical cycling performance, we performed molecular dynamics (MD) simulations using a series of UiO MOFs with different linker sizes including LiPSs and electrolyte components (1M LiTFSI in DME/DOL 1:1 volume) in pores. Specifically, UiO-66, UiO-67, and UiO-68 are Zr-based MOFs with a linker length of 5.67 Å, 10.07 Å, and 14.35 Å, respectively, leading to different pore diameters of 8.4 Å, 12.8 Å, and 17.8 Å. We calculated the mean-squared displacements of Li and LiPS from MD simulations at 300 K for 50 ns or over to assess their mobility in MOFs.

25 ENERGY STORAGE↗
DOE OSTI2026

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗