Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “current distribution mapping”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Developing novel electrodes with ultralow catalyst loading for high-efficiency hydrogen production in proton exchange membrane electrolyzer cells

Hydrogen plays more crucial roles for decarbonizing the planets and meeting the climate challenges because of its high energy density and zero-emission. It can be produced with proton exchange membrane electrolyzer cells (PEMECs) driven by sustainable and renewable energy resources. Although PEMECs have a number of advantages, including high purity production, quick response, and the ability to operate at high pressure facilitating the gas delivering, their performance and cost greatly hinder their commercial-scale applications. To achieve high-efficiency and cost-reduced hydrogen production in PEMECs, we proposed thin engineered liquid/gas diffusion layers (LGDLs) and associated electrodes, i.e., catalyst-coated LGDLs (CCLGDLs), over conventional porous transport layers (PTLs) and catalyst-coated membranes (CCMs). The research approaches in this project are based on material synthesis, in-situ and ex-situ characterizations, component design and treatment, numerical modeling, and cost analysis. The thin and tunable LGDLs (TT-LGDLs) and CCLGDLs were successfully developed with great performance improvement as demonstrated in lab-scale, bench-scale, and system-scale electrolyzer tests. The electrode thickness was reduced from 370 µm to less than 100 µm with simplified fabrication processes. With the catalytically enhanced Ir-based catalyst coating, the as-developed CCLGDLs with a catalyst loading of 0.34 mg Ir /cm 2 achieved a cell performance of 1.77 V at 2 A cm -2 , exhibiting the catalyst mass activity enhanced by >20 times with significant catalyst saving over conventional catalyst cell design. In-situ PEMEC characterizations, including the current distribution mapping and high-speed and multiscale visualizations, were conducted for a deeper understanding of mass transport and electrochemical reactions within an electrolyzer with LGDLs and CCLGDLs. A 2D cell model was developed and validated for the enhanced performance on TT-LGDL through reducing ohmic losses due to nonuniform hydration and water transport. Further, the cost analysis results have shown a path to move beyond equivalency and surpass costs associated with the project baseline. In this project, the design and fabrication of TT-LGDLs and CCLGDLs will contribute to the performance enhancement, manufacturing simplification, and cost reduction for PEMECs and other energy conversion devices, thus shortening their pathways towards commercialization. This project also provides a good foundation for furthering the in-situ reaction interface research.

08 HYDROGEN↗
DOE OSTI2024

EWTN: Quantifying Mass Transport to Enable Water Electrolyzer Architectures with Low Flow-Rate Sensitivity

To develop cost-effective and high-performing polymer electrolyte water electrolyzers (PEWEs) for gigawatt-scale applications, researchers have focused on reducing precious metal catalyst loadings and optimizing porous transport layers. However, the performance of PEWEs is also affected by system architecture. Mass transport losses are dependent on localized architecture and material interactions. In-situ measurements, such as current density distribution maps have demonstrated advantages in understanding the intricate characteristics and influence of two-phase flow within PEWEs. This study proposes the parameter of effective water transport number (EWTN) as a quantitative tool to investigate such current density distribution (CDD) measurements for PEWEs. Results show that higher flow-rates have EWTN values of 0.95 and above, indicating no mass transport limitations; while lower flow-rates with large gradients CDD have EWTN values between 0.6-0.8, indicating mass transport limited conditions. The new analysis also identified a correlation between mass transport losses due to bubble accumulation, membrane hydration, and ohmic overpotentials. To address these limitations, an unitized pin-type LGDL/flow-field design was developed, which effectively prevents local gas phase accumulation, resulting in improved mass transport characteristics. The results of this work show reduced flow-rate sensitivity with the pin-type architecture and ∼13% increased performance at 0.24 ml / min / cm 2 .

Electrochemistry↗
DOE OSTI2022

Reflection Optical Imaging to Study Oxygen Evolution Reactions

Discovery of high-efficiency water-splitting catalysts is critical for clean hydrogen fuel generation. Quantitative and reliable chemical imaging methods can play a crucial role in providing critical catalyst efficiency information for materials development and screening. Here, we developed a reflection optical imaging method to map the performance of transition metal oxide catalysts for oxygen evolution reaction (OER). The generated high concentration of oxygen in the electrolyte will decrease the refractive index of the solutions near the electrode and change the reflectivity of optical images. We successfully imaged and measured the OER on cobalt oxide, and the localized current density and overpotential distribution maps have been extracted. The oxidation and reduction of cobalt oxide can be imaged with our reflection optical imaging method as well. These maps showed that the catalytic activities are quite different at different locations, and some regions showed much higher activities. This is due to the non-uniformity created by the catalyst synthesis process and demonstrates the importance of using the optical imaging method to map the localized catalytic efficiency. In addition, the NiFe alloy’s OER catalytic activities have been studied as well. The proposed reflection optical imaging method could find many applications in catalyst activity measurement and screening.

Electrochemistry↗
DOE OSTI2025

Unimodal Imaging of Monovalent Metal-Chelator Complexes and Lipids by MALDI Imaging Mass Spectrometry

Careful regulation of monovalent metal ions (M + ) is necessary to maintain a functional cellular system. Of these ions, appropriate sodium (Na + ) and potassium (K + ) concentrations are particularly integral for electrochemical signaling, as well as the secondary transport of nutrients and waste. Dysregulation of M + homeostasis can disrupt these mechanisms, potentially influencing the metabolism of downstream biomolecules such as lipids. Thus, the relationship between M + abundances and related biomolecular distributions must be elucidated to better understand the physiology of healthy and disordered tissues. Traditional techniques for imaging biological metal distributions include SIMS, LA-ICP-MS, and XRF; however, these capabilities are limited to elemental analysis or the analysis of molecular fragments and must be paired with other modalities to visualize distributions of more complex biomolecules within the same or similar samples. Conversely, matrix-assisted laser desorption/ionization imaging mass spectrometry (MALDI IMS) is a powerful tool often used for mapping such biomolecular distributions, but current methods are unable to detect metals within tissue. This study illustrates a novel methodology that adds metal detection to the MALDI IMS repertoire through which the simultaneous detection of M + metals and lipids is achievable. Using a robotic sprayer for homogeneous application, on-tissue deposition of the chelator deferiprone (DEF) enables subsequent detection of the ionizable metal-chelator complex by MALDI without hindering lipid detection. Our work provides proof-of-concept data for the simultaneous detection of K + , Na + , and intact lipids using MALDI IMS.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A Scalable Gaussian Process Approach to Shear Mapping with MuyGPs

Analysis of cosmic shear is an integral part of understanding structure growth across cosmic time, which in turn provides us with information about the nature of dark energy. Conventional methods generate shear maps from which we can infer the matter distribution in the universe. Current methods (e.g., Kaiser–Squires inversion) for generating these maps, however, are tricky to implement and can introduce bias. Recent alternatives construct a spatial process prior for the lensing potential, which allows for inference of the convergence and shear parameters given lensing shear measurements. Realizing these spatial processes, however, scales cubically in the number of observations—an unacceptable expense as near-term surveys expect billions of correlated measurements. Therefore, we present a linearly scaling shear map construction alternative using a scalable Gaussian process prior called MuyGPs. MuyGPs avoids cubic scaling by conditioning interpolation on only nearest neighbors and fits hyperparameters using batched leave-one-out cross-validation. This work is the first step toward a full, scalable mass mapping method. We work in a simplified regime where we validate our method by interpolating and analyzing maps given noisy point-estimate data from all three shear fields, taken from a suite of N -body ray-tracing simulations. We also show that we can perform these operations at the scale of billions of galaxies on high-performance computing platforms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dynamic Heat Flow and Current Distribution Analysis in the Bottom Anode of an Electric Arc Furnace Using Fiber-Optic Sensors

A reliable method for monitoring bottom anode wear during DC Electric Arc Furnace (DC-EAF) operation is of critical importance for safe and efficient steel production. Underestimation of bottom wear poses a serious safety risk that must be avoided, while overestimation of bottom wear also poses challenges, as premature anode replacement is expensive and affects EAF productivity. Previously, we demonstrated that fiber-optic sensors can be successfully deployed to create a spatially distributed temperature map to monitor the health of the anode. The present work explores the heat flow and current density distribution in bottom anode pins to predict bottom wear, steel penetration events, and monitor refractory erosion. Small dynamic variations in pin temperature induced by joule heating during arcing also provide a means to observe local current flows in each pin. When mapped, these measurements provide a real-time view of the non-uniform and dynamic current flow in the bottom anode during EAF operation that can affect bottom wear.

Bottom Anode↗
DOE OSTI2021

Updated distribution maps of predominant Culex mosquitoes across the Americas

Abstract Background Estimates of the geographical distribution of Culex mosquitoes in the Americas have been limited to state and provincial levels in the United States and Canada and based on data from the 1980s. Since these estimates were made, there have been many more documented observations of mosquitoes and new methods have been developed for species distribution modeling. Moreover, mosquito distributions are affected by environmental conditions, which have changed since the 1980s. This calls for updated estimates of these distributions to understand the risk of emerging and re-emerging mosquito-borne diseases. Methods We used contemporary mosquito data, environmental drivers, and a machine learning ecological niche model to create updated estimates of the geographical range of seven predominant Culex species across North America and South America: Culex erraticus , Culex nigripalpus , Culex pipiens , Culex quinquefasciatus , Culex restuans , Culex salinarius , and Culex tarsalis . Results We found that Culex mosquito species differ in their geographical range. Each Culex species is sensitive to both natural and human-influenced environmental factors, especially climate and land cover type. Some prefer urban environments instead of rural ones, and some are limited to tropical or humid areas. Many are found throughout the Central Plains of the USA. Conclusions Our updated contemporary Culex distribution maps may be used to assess mosquito-borne disease risk. It is critical to understand the current geographical distributions of these important disease vectors and the key environmental predictors structuring their distributions not only to assess current risk, but also to understand how they will respond to climate change. Since the environmental predictors structuring the geographical distribution of mosquito species varied, we hypothesize that each species may have a different response to climate change. Graphical abstract

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Current and Projected Future Distribution of Plant Community Types for the Southern Seward Peninsula, Alaska

Landscape scale maps of plant community distribution at 5 m resolution were developed for southern Seward Peninsula. A Random Forest based environmental niche model was developed using topographic and climate data sets. Models were trained for the 2010-2019 period and applied to develop plant community distributions for contemporary (2010-2019) period. Trained model was applied to project the plant community distributions under future climate. Using downscaled projections for RCP8.5 scenario from five climate models (CCSM4, GFDL-CM3, GISS-E2-R, IPSL-CM5A-LR, MRI-CGCM3) and multi-model mean, a six member ensemble of plant community distribution were developed for current decade (2010-2019) and for future decades (2020-2029, 2030-2039, 2040-2049, 2050-2059). Contains *.tif and *.sh files. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulating water dynamics related to pedogenesis across space and time: Implications for four-dimensional digital soil mapping

Digital soil mapping (DSM) relies on machine-learning and geostatistics to represent soil property observations across space. DSM techniques are powerful but often empirical, being limited to the quality and density of point samples. Water dynamics are closely related to soil variability, and the physics that govern water movement are well known. Hydrological properties can hence be simulated by physical models through space and time, unveiling key characteristics about soils. We propose the use of hydrologic models to map soils across the surface (2D), depth (1D), and time (1D)–which provides a 4D approach to digital soil mapping (4DSM). The Distributed Hydrology Soil Vegetation Model (DHSVM) was applied to a watershed currently under pasture. Moisture sensors and wells were installed at different depths in the watershed on summit, sideslope and toeslope positions to validate the model. DHSVM simulations of soil moisture distribution and depth to saturation were performed during the hydrological year (October 2008-September 2009). Clusters of similar pixels based on soil moisture values were determined using Dynamic Time Warping (DTW) to align temporal data and K-means. Clustering was performed both seasonally and for the entire year. Temporal patterns simulated by DHSVM matched measurements given by moisture sensors and wells. Seasonal clusters differed from the annual cluster. Distinct clusters were observed for each season and with depth, showing that spatiotemporal soil variability is lost when statically assessing soils. Spatiotemporal clusters corroborated field observations of fragipan occurrence not explicitly spatially mapped by Soil Survey Geographic Database (SSURGO). If a connection can be made between water and soils, static and dynamic soil variability can be predicted using physically based hydrologic models. Hydrologic models can benefit soil mapping by enabling reliable 4D simulation of water dynamics, which are fundamental to soil variability and soil classification and directly relate to biological, physical and chemical soil processes not captured by typical soil sampling protocols.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Improved Vertical Carrier Transport for Green III-Nitride LEDs Using ( In , Ga ) N Alloy Quantum Barriers

We report on experimental and simulation-based results using ( In , Ga ) N alloy quantum barriers in c -plane green light-emitting diode (LED) structures as a means to improve vertical carrier transport and reduce forward voltage ( V F ) . Three-dimensional device simulations that include random alloy fluctuations are used to understand carrier behavior in a disordered potential. The simulated current density–voltage ( J - V ) characteristics and modified electron-hole overlap | F mod | 2 indicate that increasing the indium fraction in the ( In , Ga ) N quantum barriers leads to a reduced polarization discontinuity at the interface between the quantum barrier and quantum well, thereby reducing V F and improving | F mod | 2 . Maps of electron and hole current through the device show a relatively homogenous distribution in the X Y plane for structures using Ga N quantum barriers; in contrast, preferential pathways for vertical transport are identified in structures with ( In , Ga ) N barriers as regions of high and low current. A positive correlation between hole (electron) current in the p -side ( n -side) barrier and indium fraction reveals that preferential pathways exist in regions of high indium content. Furthermore, a negative correlation between the strain ε z z and indium fraction shows that high indium content regions have reduced strain-induced piezoelectric polarization in the Z direction due to the mechanical constraint of the surrounding lower indium content regions. Experimentally, multiple quantum well green LEDs with ( In , Ga ) N quantum barriers exhibit lower V F and blue-shifted wavelengths relative to LEDs with Ga N quantum barriers, consistent with simulation data. These results can be used to inform heterostructure design of low V F , long-wavelength LEDs and provide important insight into the nature of carrier transport in III-nitride alloy materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Tree Canopies Reflect Mycorrhizal Composition

Mycorrhizae alter global patterns of CO 2 fertilization, carbon storage, and elemental cycling, yet knowledge of their global distributions is currently limited by the availability of forest inventory data. Here, we show that maps of tree-mycorrhizal associations (hereafter “mycorrhizal maps”) can be improved by the novel technology of imaging spectroscopy because mycorrhizal signatures propagate up from plant roots to impact forest canopy chemistry. We analyzed measurements from 143 airborne imaging spectroscopy surveys over 112,975 individual trees collected across 13 years. Results show remarkable accuracy in capturing ground truth observations of mycorrhizal associations from canopy signals across disparate landscapes (R 2 = 0.92, p < 0.01). Additionally, upcoming imaging spectroscopy satellite missions can reveal new insights into landscape-scale variations in water, nitrogen, phosphorus, carotenoid/anthocyanin, and cellulose/lignin composition. Applied globally, this approach could improve the spatial precision of mycorrhizal distributions by a factor of roughly 104 and facilitate the incorporation of dynamic shifts in forest composition into Earth system models.

58 GEOSCIENCES↗
DOE OSTI2024

Rapid mapping of electrochemical processes in energy-conversion devices

Electrochemical impedance spectroscopy (EIS) is ubiquitously applied to identify physicochemical processes governing the performance of energy-conversion devices. However, deconvolution and interpretation of impedance phenomena are limited by measurement throughput and a dearth of scalable analysis methods. Here, we demonstrate an approach to quickly collect and coherently analyze large volumes of electrochemical data. In this study, we accelerate impedance characterization by combining rapid measurements in time and frequency domains, which are interpretably transformed using the distribution of relaxation times (DRT) and a new distribution of phasances (DOP) model. This method provides excellent agreement with EIS and decreases measurement time by an order of magnitude. High-throughput spectra are then distilled into detailed electrochemical maps. This approach is applied to a Li-ion battery and a protonic ceramic electrochemical cell as practical case studies, demonstrating how mapping can richly characterize physicochemical relationships that are difficult to decipher with conventional measurement and analysis methods.

25 ENERGY STORAGE↗
DOE OSTI2020

Sample preparation with sucrose cryoprotection dramatically alters Zn distribution in the rodent hippocampus, as revealed by elemental mapping

Transition metal ions (Fe, Mn, Cu, Zn) are essential for healthy brain function, but altered concentration, distribution, or chemical form of the metal ions has been implicated in numerous brain pathologies. Currently, it is not possible to image the cellular or sub-cellular distribution of metal ions in vivo and therefore, studying brain-metal homeostasis largely relies on ex vivo in situ elemental mapping. Sample preparation methods that accurately preserve the in vivo elemental distribution are essential if one wishes to translate the knowledge of elemental distributions measured ex vivo toward increased understanding of chemical and physiological pathways of brain disease. Furthermore, the choice of sample preparation is particularly important for metal ions that exist in a labile or mobile form, for which the in vivo distribution could be easily distorted by inappropriate sample preparation. One of the most widely studied brain structures, the hippocampus, contains a large pool of labile and mobile Zn. Herein, we describe how sucrose cryoprotection, the gold standard method of preparing tissues for immuno-histochemistry or immuno-fluorescence, which is also often used as a sample preparation method for elemental mapping studies, drastically alters hippocampal Zn distribution. Based on the results of this study, in combination with a comparison against the strong body of published literature that has used either rapid plunge freezing of brain tissue, or sucrose cryo-protection, we strongly urge investigators in the future to cease using sucrose cryoprotection as a method of sample preparation for elemental mapping, especially if Zn is an analyte of interest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Single-shot ionization-based monitor for pulsed electron beams

In this work, we present an experimental demonstration of a single-shot, nondestructive electron beam diagnostic based on the ionization of a low-density pulsed gas jet. In our study, 7 MeV electron bunches from a radio-frequency photoinjector, carrying up to 100 pC of charge, traversed a localized distribution of nitrogen gas (N 2 ). The interaction of the electron bunches with the N 2 gas generated a correlated signature in the ionized particle distribution, which was spatially magnified using a series of electrostatic lenses and recorded with a microchannel-plate detector. Various modalities, including point-to-point imaging and velocity mapping, are investigated. A temporal trace of the detector current enabled the identification of single- and double-ionization events. The characteristics of the ionization distribution, dependence on gas density, total bunch charge, and other parameters, are described. Approaches to scaling to higher electron bunch density and energy are suggested. Additionally, the instrument proves useful for comprehensive studies of the ionization process itself.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Weak lensing in DES Y3: redshift distributions, shape catalogue, and mass mapping

The size of galaxy surveys have rapidly escalated with time, growing from thousands of galaxies of the first pioneering surveys to current and future surveys (DES, KiDS, HSC, LSST, Euclid, WFIRST), which map (or plan to map) a few hundreds millions of galaxies. These ambitious programs have the potential to answer many open questions in Cosmology (or, depending on their results, to open new exciting scenarios, involving new physics beyond the standard cosmological model and the standard model of particle physics). In particular, this thesis focuses on the first three years of data from the Dark Energy Survey (DES Y3), which comprises 100 million galaxies observed in 5 different optical bands over 5000 square degrees of the southern hemisphere. DES and the data set used in this thesis are discussed in detail in Chapter 3, which together with Chapter 1 and 2 constitutes Part I of this thesis. The scientifically positive outcome of DES and other photometric galaxy surveys goes h and in h and with a supreme understanding of the data, of the observational systematics affecting the measurements, and of the limits of the theoretical modelling of the observables. One of the main challenge of galaxy photometric surveys like DES is to correctly estimate the distance to the galaxies. Photometric surveys observe the galaxies in a limited number of broad optical bands, and the estimate of their \textit{redshifts} is somewhat uncertain, which constitutes one of the main limiting factor of their cosmological analyses. Part II and its Chapter 4 are devoted to a new auxiliary technique recently developed that allows to infer (or calibrate) redshift distributions exploiting cross-correlation with small reference samples with secure redshifts. The so called ``clustering-redshift'' method is discussed in depth, as well as its role in the overall DES Y3 redshift calibration strategy. Part III and Chapter 5 are devoted to another key aspect of photometric surveys, that enables the stu dy of weak gravitational lensing: the shape catalogue. The official DES Y3 shape catalogue is presented, together with all the systematic tests performed to validate it such that it can be considered science-ready. Last, part IV focuses on how to use weak gravitational lensing and the DES Y3 shape catalogue to create weak lensing mass maps, i.e., maps of the mass distribution of the Universe. Chapter 6 discusses different map making techniques, whereas Chapter 7 presents a simulated cosmological analysis using the high order moments of the mass maps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Microscopic analysis of copper current collectors and mechanisms of fragmentation under compressive forces

Extensive fragmentation of copper current collectors was observed after spherical indentation on prismatic and large-format pouch Li-ion cells by 3D X-ray computed tomography (XCT). Microscopic analysis including scanning electron microscopy (SEM), scanning transmission electron microscopy (STEM) and x-ray photoelectron microscopy (XPS) was carried out on copper current collectors from used commercial cells and pristine anodes. The copper-graphite cross-section images showed rough interface areas affected by reactions and diffusion in the used cell. Electron probe micro-analyzer (EPMA) element mapping showed the interface area was rich in oxygen and phosphorus. A detectable amount of phosphorus was also uniformly distributed inside the current collector. The same oxygen and phosphorus distributions were confirmed by STEM/EDS analysis. XPS depth profiles on multiple elements revealed the interface area of the aged anode was rich in Li, F, P, O and C and diffused at least 50 nm into the copper. In comparison, the pristine anode showed a very smooth C/Cu interface. No other elements were detected. For commercial cells, the reactions in the interface area and diffusion of multiple elements into the lattice and grain boundaries were responsible for the embrittlement of the copper current collectors. Finally, permanent cell capacity loss was observed in electrochemical performance of the indented cells.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mapping valence electron distributions with multipole density formalism using 4D-STEM

Recent advancement in aberration correction and detector technology opened a door to various applications using 4D-STEM, which yields a diffraction pattern for each scanning position within a crystal unit-cell in scanning transmission electron microscopy (STEM) and generates incredible amounts of data in momentum space. Currently 4D-STEM analysis relies on the center-of-mass of the diffraction patterns in electric field and charge density mapping. It only derives the total projected charge density and is limited to phase objects, e.g. extremely thin samples. Here, we propose a new analytical method to accurately map aspherical valence electron distributions with atom-centered multipolar functions formalism using the whole 4D-STEM dataset. We demonstrate that, with the full dynamical calculations for various sample thicknesses, the method is sensitive not only to the miniscule charge transfer, but also to the atomic site symmetry and aspherical electron orbitals. The process of the refinement is much more robust and reliable than quantitative convergent beam electron diffraction.

4D-STEM↗
DOE OSTI2025

Multiscale Modeling Framework for Lithium Nucleation in 3D Porous Carbon Anodes

Porous carbon scaffolds offer a promising route for mitigating non-uniform lithium (Li) plating to enhance the safety and longevity of Li metal batteries. However, the influence of microstructural morphology on Li nucleation is not well understood. Here, we present a multiscale modeling framework to investigate how the porous microstructure of carbon materials affects Li nucleation behavior. Ab initio molecular dynamics simulations quantify the nucleation energy barriers of Li on graphene as a function of Li content, surface curvature, and applied potential, providing key parameters for a classical nucleation theory (CNT) model. From macroscale half-cell simulations, we obtained Li concentration and electrical potential profiles to define boundary conditions for mesoscopic simulations. At the mesoscale, three distinct synthetic 3D microstructures with different porosities and characteristic feature sizes are generated to resolve local distributions of Li flux, current density, and mechanical stress. These outputs are integrated into the CNT model to map spatial variation in nucleation rates. Our findings reveal trade-offs between suppressing nucleation rates and achieving spatial uniformity, offering design guidelines for optimizing porous carbon anodes to balance nucleation control and mechanical integrity.

Materials science↗