Engineering PapersSearch

SEARCH · Engineering Papers

Results for “electrode geometry”

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

Reaction Diffusion Modelling of 3D Pillar Electrodes in Single-Catalyst CO 2 Reduction Cascades

Effective electrochemical CO 2 reduction to liquid fuels requires that the local catalytic environment facilitates the desired reactivity, yet a microscopic understanding of this environment is difficult to achieve from experiment alone. In this work, a 3D reaction-diffusion model was developed to explore the effects of electrode surface area and local geometry on the performance of a heterogeneous catalyst that performs a two-step CO 2 reduction cascade reaction to CO and then CH 3 OH under aqueous conditions. Kinetic parameters for the model were inspired by experimental results using a cobalt phthalocyanine (CoPc) catalyst. Three-dimensional architectures composed of arrays of square pillars with varying dimensions and either smooth or periodically modulated surfaces were tested, revealing the extent to which geometry modulates the performance of the cascade reactions. Although structural variations modulate local concentration gradients, we find that electrochemically active surface area predominantly governs the overall cascade reaction. Moreover, the results suggest that supersaturation of CO, with concentrations up to ten-fold higher than the equilibrium solubility limit, might be critical for more efficient conversion to CH 3 OH. For any given geometry, the spatially averaged ratio of [CO] to [CO 2 ] is dictated by the electrochemically active surface area and determines the yield of CH 3 OH. For a fixed surface area, geometries that spatially confine the electrolyte yield moderate local [CO] to [CO 2 ] ratios within small volumes. In contrast, less confining geometries result in a broader distribution of local ratios spread over larger volumes, with both configurations yielding the same spatially averaged [CO] to [CO 2 ] ratio. These insights provide valuable design principles—highlighting the critical importance of surface area and possibly CO supersaturation—for engineering advanced electrode architectures that leverage intermediate trapping and CO supersaturation to enhance overall performance in tandem CO 2 reduction systems.

COMSOL

In Situ Liquid Electron Microscope Cells Strongly Attenuate Electrochemical Behavior

In situ electrochemical Scanning Transmission Electron Microscopy (ec-S/TEM) probes the dynamics of electrochemical processes in real time and at high spatial resolution. This method utilizes microfabricated electrochemical cells, with enclosed dimensions of 800 μm wide, 800 μm long, and 1–2 μm tall, with electron-transparent silicon nitride windows. Working, counter, and reference electrodes in the configuration studied are deposited Pt with a planar surface area of about 1 μm 2 and a height of 50 nm. The cell confines the electrolyte within a 1 μm tall channel, and the microfabricated electrodes have a non-standard geometry, are coplanar, and are not front-facing. As such, standard assumptions of electrochemical experiments do not apply. COMSOL modeling was used to compare and contrast differences in electrochemical behavior between conventional experimental setups and the in situ ec-S/TEM cells. Cell height strongly affected voltammogram measurements: peak magnitudes and shifts correlate directly to cell height. Reactant concentration above the electrode quickly decreases during deposition, leading to a planar diffusion-dominated regime. The concentration gradient prioritizes particle growth along the thin edges of the electrode and not on the planar part of the electrode. Experiments done in the ec-S/TEM ex situ of a TEM verify the model’s accuracy.

Wittman, Reed M. [Oak Ridge National Laboratory (O

Membrane Degradation in PEM Fuel Cells: Part I. Modeling Gas Crossover and the Pt Band

Understanding chemical degradation of the proton-exchange membrane in fuel cells is crucial for extending their lifetimes. Herein, various degradation reactions reported in literature are organized and analyzed, including direct radical generation and an indirect (Fenton) pathway. To understand the transport of dissolved H 2 and O 2 crossover gases as they relate to membrane degradation, an agglomerate-scale model is introduced, treating gas, ionomer, and catalyst as discrete phases. The model reveals a key phenomenon: at working potentials, dissolved gases are mostly consumed at the interface between the catalyst layer and the membrane, leaving little gas to cross the membrane. Under open-circuit conditions, dissolved gases are not consumed and can then cross the membrane. This explains high H 2 O 2 concentrations and degradation rates seen in experiments but not captured in previous models. Following mixed-potential theory, crossover gases supply the hydrogen-oxidation and oxygen-reduction reactions, which occur simultaneously on individual Pt particles comprising the Pt band in the membrane, forming reactive species (H 2 O 2 , OH·). Results show crossover gas almost entirely reacts on the Pt band, allowing little to reach the opposite electrode. Furthermore, the micro-scale geometry of the catalyst-layer/membrane interface impacts the gas crossover at working potentials, indicating that cell construction affects membrane durability.

Johnson, Evan F. [Lawrence Berkeley National Labor

Random 3D Interpenetrating Electrode Design for Energy Storage Applications

We introduce the concept of random interpenetrating electrode structures, inspired by spinodal decomposition, as a transformative approach for energy storage applications. While ordered architectures, such as double gyroids and Kelvin cell structures, have long been favored for their structural uniformity and surface area utilization, we demonstrate that random structures offer distinct advantages including shorter diffusion pathways, higher packing density, and enhanced reaction kinetics, particularly at low temperatures. Using a combination of computational and electrochemical analysis, we show that random architecture outperforms their ordered counterparts in scenarios where diffusion dominates. This work represents the first systematic exploration of random structures in energy storage, highlighting their potential to overcome limitations of traditional designs. However, challenges such as structural reproducibility and scalability remain, necessitating further investigation into fabrication techniques and material selection. By establishing a framework to optimize random electrode architectures, we provide critical insights into the interplay between geometry, diffusion, and reaction kinetics. This study not only introduces a new design principle for 3D electrodes but also opens pathways for next-generation energy storage systems that demand both high performance and adaptability to a range of operating conditions.

Materials science

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Feasibility Assessment of Using Electrical Impedance Tomography for Damage Localization in Graphite Microreactor Components

Microreactors have great potential to decrease capital costs and construction timelines, reducing the barriers to implementing advanced nuclear reactor technologies. However, the lower power output of these microreactors introduces economic challenges if the operation and maintenance costs cannot be reduced sufficiently. Many microreactor concepts use graphite materials for in-core moderator and structural purposes, which will require periodic inspection or, ideally, in situ structural health monitoring. This work describes an initial evaluation of the feasibility of leveraging graphite’s semiconducting properties to perform electrical impedance tomography (EIT) for defect localization. First, a study was performed to identify the best methods for bonding electrical contacts to minimize the effects of contact resistance that interfere with impedance measurements of the graphite. After determining the best bonding approach, electrodes were attached to multiple graphite components with varying geometries (e.g., simple block and more representative microreactor hexagonal block). In parallel, finite element analysis approaches were developed and implemented to analyze the impact of defects and contact resistances and to inform an iterative inverse model for reconstructing the conductivity distribution. The results of this study show that defect localization in graphite components is possible using EIT if a sufficient number of electrodes (to improve spatial resolution) can be bonded using a technique with very low contact resistance and those contact resistances remain stable during reactor operation. If the reactor core and vessel design do not allow electrodes to be bonded during reactor operation, then it may be possible to detach/reattach electrodes between reactor operation and outages if a suitable mechanical connection that offers repeatable contact resistance can be identified. Alternatively, future work could focus on identifying contacts that can survive the harsh reactor operating environment so that contacts do not need to be removed.

22 GENERAL STUDIES OF NUCLEAR REACTORS

Investigating instabilities in magnetized low-pressure capacitively coupled RF plasma using particle-in-cell (PIC) simulations

The effect of a uniform magnetic field on particle transport in low-pressure radio frequency (RF) capacitively coupled plasma (CCP) has been studied using a particle-in-cell model. Three distinct regimes of plasma behavior can be identified as a function of the magnetic field. In the first regime at low magnetic fields, asymmetric plasma profiles are observed within the CCP chamber due to the effect of $\overrightarrow{E}$ × $\overrightarrow{B}$ drift. As the magnetic field increases, instabilities develop and form self-organized spoke-shaped structures that are distinctly seen within the bulk plasma closer to the sheath. In this second regime, the spoke-shaped coherent structures rotate inside the plasma chamber in the -$\overrightarrow{E}$ × $\overrightarrow{B}$ direction, where $\overrightarrow{E}$ × $\overrightarrow{B}$ are the DC electric and magnetic field vectors, respectively, and the DC electric field exists in the sheath and pre-sheath regions. The spoke rotation frequency is in the megahertz range. As the magnetic field strength increases further, the rotating coherent spokes continue to exist near the sheath. The coherent structures are, however, accompanied by new small-scale incoherent structures originating and moving within the bulk plasma region away from the sheath. This is the third regime of plasma behavior. The threshold values of the magnetic field between these regimes were found not to vary with changing plasma reactor geometry (e.g., area ratio between ground and powered electrodes) or the use of an external capacitor between the RF-powered electrode and the RF source. The threshold values of the magnetic field between these regimes shift toward higher values with increasing gas pressure. Analysis of the results indicates that the rotating structures are due to the lower hybrid instability driven by density gradients and electron-neutral collisions. This paper provides guidance on the upper limit of the magnetic field for instability-free operation in low-pressure CCP-based semiconductor deposition and etch systems that use the external magnetic field for plasma uniformity control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Dilute Regeneration-Driven Membrane Capacitive Deionization of Synthetic Seawater Using Nanopatterned Membranes and Prussian Blue Analog Electrodes

Membrane capacitive deionization (MCDI) offers energy-efficient seawater desalination but is limited at high salinity by membrane resistance and incomplete electrode regeneration. Nanopatterned ion-exchange membranes, dilute regeneration protocols, and Prussian blue analog (PBA)-functionalized electrodes are combined in a flow-by-MCDI cell. Nanopatterned ion-exchange membranes (hexagonal, octagonal, double-ring, rectangular) enhance interfacial ion transport, with hexagonal geometry delivering ≈12.5% greater surface area and the best performance. PBA-functionalized electrodes increase salt adsorption and charge-transfer kinetic rates. The integrated system lowers the area-specific resistance by 45 Ω cm2, resulting in a 500 mV reduction in the cell voltage for a current density of 2 mA cm−2 for a 35 000 ppm NaCl feed. This improves the energy-normalized salt adsorption six fold (64–382 mmol J−1). Low salinity (2000 ppm) and mixed-salt regeneration sustains a ≈39% water recovery and stable performance for at least seven cycles. Overall, combining nanopatterned membranes, which promote confinement-enhanced ion mobility, and PBA electrodes, which enhance salt adsorption, improved the energy efficiency of MCDI.

Hasan, Mahmudul

Viability of Additively Manufactured Electrodes for Lithium-Ion Batteries

As the global economy becomes increasingly electrified, the demand for batteries and energy storage is expected to rise significantly, particularly in the transportation and electricity sectors. Lithium-ion batteries (LIBs) are currently the most advanced and widely used technology in this field. Traditionally, LIBs are manufactured using simple 2D planar geometries to maximize production efficiency and minimize costs. However, this approach limits energy density due to the restricted design flexibility of the electrodes. Additive manufacturing (AM) offers a promising solution to enhance the energy density and efficiency of LIBs by enabling the design of architectures that reduce diffusive losses and allow for a greater amount of active material to be incorporated within the same device footprint, thereby minimizing the use of inactive materials. Different AM techniques come with their own set of limitations, including printing speed, material compatibility, and scale, which must be considered when designing electrodes. Scalable and cost-effective methods are particularly important for electric vehicle batteries, while achieving higher energy densities in microbatteries is crucial for the miniaturization of wearable electronics and medical devices. Here, in this study, we simulate various 3D porous electrode designs for LIBs using graphite and nickel manganese cobalt oxide (NMC) electrodes. These designs are selected to represent structures that could be produced using different AM techniques, such as direct ink writing, fused deposition modeling, and stereolithography. Our results indicate that at higher charging rates and increased areal mass loading, 3D structures can outperform traditional 2D electrodes, although the benefits may diminish with more complex designs that are harder to manufacture. The observed gains in energy density are attributed to improved electrode utilization and reduced diffusive energy losses. This comprehensive analysis of structure–performance relationships will provide valuable insights to guide future research on 3D designs, material selection, and AM techniques for additively manufactured battery electrodes.

25 ENERGY STORAGE

Design considerations for optimizing the transient CHI injector on QUEST

Transient coaxial helicity injection (transient CHI) on the QUEST (Q-shu University Experiment with Steady-State ST) spherical tokamak (ST) has recently validated the floating biased electrode configuration for solenoid-free plasma startup. In support of a significant divertor upgrade on QUEST, the details of the transient CHI injector geometry on QUEST have been examined more closely using tokamak simulation code (TSC) simulations. QUEST uses a HIT-II-like (Helicity Injected Torus) injector configuration in which the injector region is comprised of coaxial injector electrodes located in the lower part of the machine. TSC simulations indicate high injector flux operation may benefit from an increased gap between the coaxial electrodes. This is the first study to examine the differences between an open and closed electrode configuration for solenoid-free plasma startup using transient CHI. Results show that both configurations can generate similar levels of closed flux, but the closed configuration may be easier to implement in some ST designs as a much smaller portion of the injector needs to be close to the injector flux coil. Results show that increasing the electrode gap width from the present 10.8 cm to about 15–20 cm would increase the closed flux fraction by about 40%. The results presented in this paper are generally applicable to the CHI design for other STs.

coaxial helicity injection

Spatially Patterned Architectures to Modulate CO 2 Reduction Cascade Catalysis Kinetics

Electrochemical CO 2 reduction using renewable sources of electrical energy holds promise for converting CO 2 into fuels and chemicals. The complex interactions among chemical/electrochemical reactions and mass transport make it difficult to analyze the effect of an individual process on electrode performance based only on experimental methods. Here, we developed a generalized steady-state simulation to describe an electrode surface in which sequential cascade catalysts are patterned in a periodic trench design. If appropriately constructed, this trench geometry is hypothesized to be able to yield a higher net current density for a CO 2 reduction (CO 2 R) cascade reaction. We have used realistic experimental reaction kinetics to investigate the role of trench geometry in mass transport, local microenvironments, and selectivity for a model CO 2 R cascade reaction. The model considers local concentration gradients of bicarbonate species at quasi-equilibrium and catalytic surface reactions based on concentration-dependent Butler–Volmer kinetics. Our results suggest that varying the spatial distribution of active sites plays a significant role in facilitating effective mass transport between active sites, modulating selectivity for the cascade reaction, and enhancing the yield of desirable cascade products. Moreover, we observe that this trench geometry significantly alters the cascade reaction rate by affecting the local pH, which can cause inadvertent depletion of available aqueous CO 2 to limit the CO 2 R cascade kinetics and modest suppression of the hydrogen evolution reaction (HER). The results highlight the trade-offs between mass transport, pH, and reaction kinetics that become apparent only when considering the coupled physics of all processes at the electrode surface. Here, this model can thus serve as a primary tool to build more selective and efficient patterned architectures for the CO 2 R cascade catalysis.

CO2 reduction

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS

Spin-filter tunneling detection of antiferromagnetic resonance with electrically tunable damping

Antiferromagnetic spintronics offers the potential for higher-frequency operations and improved insensitivity to magnetic fields compared to ferromagnetic spintronics. However, previous electrical techniques to detect antiferromagnetic dynamics have utilized large, millimeter-scale bulk crystals. In this work, we demonstrate direct electrical detection of antiferromagnetic resonance in structures on the few-micrometer scale using spin-filter tunneling in platinum ditelluride (PtTe 2 )/bilayer chromium sulfide bromide (CrSBr)/graphite junctions in which the tunnel barrier is the van der Waals antiferromagnet CrSBr. This sample geometry allows not only efficient detection but also electrical control of the antiferromagnetic resonance through spin-orbit torque from the PtTe 2 electrode. The ability to efficiently detect and control antiferromagnetic resonance enables detailed studies of the physics governing these high-frequency dynamics.

Cham, Thow Min Jerald [Cornell University, Ithaca,

Concavity-based local erosion and sphere-size-based local dilation applied to lithium-ion battery electrode microstructures for particle identification

Performance metrics of lithium-ion batteries can be extracted from the analysis of electrode microstructures nanoscale imaging. The characterization workflow can involve a challenging particle identification, or instance segmentation, step. In this work, we propose a new identification method based on an original transformation: a sphere-size-based local dilation followed by a concavity-based local erosion, that is local morphology closing. The new transformation is much more efficient than the global morphology closing, with correct identification achieved with only 1.7 % dilation volume and 2.6 % erosion volume on a test geometry, compared to 39.2 % and more than 50 %, respectively, with its global counterpart. The new method has been then benchmarked versus other identification algorithms (watershed and pseudo coulomb repulsive field) on a real electrode microstructure with equal or better segmentation achieved.

25 ENERGY STORAGE

Rheo–Electric Foundations and Engineering of Carbon Slurry Flow Electrodes

Flowable carbon slurry electrodes promise continuous, scalable electrochemical desalination and energy storage when their rheology (flow and stability) and electronic connectivity are jointly optimized. This report integrates six complementary thrusts that we executed under the Laboratory Directed Research and Development (LDRD) project: (i) a fast, ex situ, centrifuge-compatible method to quantify pellet (packed slurry) electronic conductivity, with a validated cell constant (“shape factor”) and quantitative structure–property insights for carbon mixtures; (ii) a controlled circulation study showing that peristaltic pumping can irreversibly alter particle size and morphology over laboratory-relevant durations; (iii) development of a miniature jet pump that preserves particle integrity while achieving entrainment ratios of 2–4 at practical flow rates; (iv) a comprehensive mapping of how a representative non-ionic surfactant (Tween 20) collapses the carbon gel network, separating network-percolation losses from intrinsic particle contact resistance; (v) in-situ visualization revealing how gravitational settling and flow-driven mixing compete to control charge transport, with dynamic bed layers preventing saturation while achieving 30-fold local concentration enhancement; and (vi) a simple, widely adoptable fabrication route to structured, gas-nitrided titanium nitride (TiN) titanium current collectors that dramatically lower interfacial losses and enable geometry-driven performance gains. Together, these results provide a coherent toolchain and design rules to accelerate formulation screening, de-risk device prototyping, and standardize characterization for carbon slurry electrodes in flow-electrode capacitive deionization (FCDI), electrochemical flow capacitor (EFC), and related systems.

25 ENERGY STORAGE

Design and Preliminary Experimental Results on a Uniform-Field Test Fixture for Power Flow Experiments

Electrodes that produce a uniform electric field across their surfaces are desired for a variety of high-voltage systems including spark-gap switches and transversely excited atmospheric lasers. Here, the implementation of uniform-field electrodes to a parallel-plate load is described for the study of vacuum power flow under pulsed-power applications. We begin with an overview of general uniform-field geometries, and provide an open-source program that generates these geometries as described by Ernst. Following that, a modification to an existing power-flow experimental platform, that adopts a uniform-field geometry, is introduced and its utility is described. Lastly, preliminary experimental results are presented and future steps are outlined.

Capacitors

Modeling Electrodeposition in 3D Porous Architectures for Solid-State Li-Metal Batteries

Li-metal storage in three-dimensional (3D) electrodes is considered a potential dendrite-mitigation strategy. The large surface area and high porosity of these electrodes result in reduced local Li-plating current densities. The porous topology provides a scaffold for Li-deposition and stripping, maintaining both mechanical integrity and Li accessibility. The goal of this study is to understand how characteristics, such as geometry and material properties, affect the current distribution and deposition pattern. To this end, we developed a computational method to track material growth driven by electrodeposition within a complex geometry. This method ensures that the finite-element discretization remains conforming to the moving boundary while preserving an adequate mesh quality, and thus maintains solution accuracy. Using this new computational tool, we analyze the conditions under which porous anode architectures effectively expand the surface area of the charge-transfer interface, and self-regulate current density and dendrite growth.

3D electrode architectures

EchemAMR (electro-chemical microsctructure scale models with adaptive meshing) [SWR-23-111]

A 3D microstructure resolving electrochemical transport and interfacial chemistry solver. Electrode microstructure plays an important role in determining the performance of an electrochemical system, e.g. lithium ion battery. EchemAMR is a microstructure scale model that solves the governing equations for ion transport, electrical current continuity, interfacial chemistry and structural mechanics. Complex microstructure geometries from imaging can be directly imported into EchemAMR. A volume fraction based description of the geometry on Cartesian grid with an immersed interface formulation enables simplified meshing and large-scale simulations with millions of degrees of freedom. EchemAMR has been tested against systems with analytic solutions for numerical convergence and highly resolved lithium ion battery microstructures. EchemAMR demonstrates excellent mass conversation and efficient scaling on heterogenous High-Performance Computing (HPC) with central and graphics processing units.

Sitaraman, Hariswaran