Search NASA⌕ Search

SEARCH · Search NASA

Results for “Microstructures”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Invertible Temper Modeling using Normalizing Flows and the Effects of Structure Preserving Loss

Advanced manufacturing research and development is typically small-scale, owing to costly experiments associated with these novel processes. Deep learning techniques could help accelerate this development cycle but frequently struggle in small-data regimes like the advanced manufacturing space. While prior work has applied deep learning to modeling visually plausible advanced manufacturing microstructures, little work has been done on data-driven modeling of how microstructures are affected by heat treatment, or assessing the degree to which synthetic microstructures are able to support existing workflows. We propose to address this gap by using invertible neural networks (normalizing flows) to model the effects of heat treatment, e.g., tempering. The model is developed using scanning electron microscope imagery from samples produced using shear-assisted processing and extrusion (ShAPE) manufacturing. This approach not only produces visually and topologically plausible samples, but also captures information related to a sample’s material properties or experimental process parameters. We also demonstrate that topological data analysis, used in prior work to characterize microstructures, can also be used to stabilize model training, preserve structure, and improve downstream results. We assess directions for future work and identify our approach as an important step towards end-to-end deep learning system for accelerating advanced manufacturing research and development.

Howland, Sylvia↗
DOE OSTI2025

Impacts of Biomass Feedstock Pre-Processing on Heat and Mass Transfer During Pyrolysis Using X-Ray Computed Tomography and Multiscale Modeling

Knowledge of the transport properties of biomass particles such as porosity, tortuosity, and permeability is paramount for high-fidelity modeling of biomass pyrolysis due to the heat and mass transfer limitations imposed by particle microstructure. X-ray computed tomography (XCT) is a non-destructive imaging method that enables full 3D reconstructions of the biomass particle microstructure with high resolution, permitting direct calculation of porosity, tortuosity, and permeability from real particle geometries. In this study, XCT imaging revealed the 3D microstructures of particles and chars from pyrolytic conversion of cylindrically cut or milled/pelletized loblolly pine samples. The porosity, tortuosity, and permeability were calculated directly from the XCT geometries via open-source microstructural analysis tool MATBOX+TauFactor (https://github.com/NREL/MATBOX_Microstructure_analysis_toolbox) and computational fluid dynamics (CFD) simulations using our solver, Mesoflow (https://github.com/NREL/mesoflow). These properties were used in a reactor scale model developed in COMSOL of the single particle reactor at NREL to investigate the impact of feedstock pre-processing on biomass conversion during pyrolysis with rigorous experimental validation.

biomass↗
DOE OSTI2025

PROCESS-STRUCTURE-PROPERTY RELATIONSHIPS IN LASER POWDER BED FUSION PRODUCED 17-4 PH STEEL

Laser powder bed fusion (LPBF) is a metal additive manufacturing method that produces non-traditional microstructures as a result of the rapid solidification and thermal cycling inherent to the process. When using LPBF-produced material in application, these unique microstructures challenge the applicability of well developed mechanical property databases achieved by conventional heat treatments. For wider adoption of this technology, a more holistic understanding is necessary on how process attributes develop material structure, which dictate mechanical properties. This dissertation explores the process– structure–property relationships in LPBF 17-4 PH steel through systematic evaluation of atmospheric processing and heat treatment effects on microstructure and mechanical performance. Specimens were fabricated under controlled build environments, subjected to a range of solutionizing, homogenizing, and aging treatments, and characterized using optical microscopy, electron back scatter diffraction (EBSD), and X-ray diffraction (XRD) to quantify phase evolution. Tensile testing was performed to directly link heat treatment pathway and nitrogen absorption to mechanical performance. This work demonstrates where conventional heat treatment standards are applicable to LPBF 17-4 PH steel and where modifications are required. By directly correlating phase stability, nitrogen effects, and tensile response, this work provides practical guidelines for tailoring post-processing strategies. These findings underscore that successful application of LPBF 17-4 PH steel requires explicit consideration of both build environment and post-processing. By linking processing conditions to microstructure and performance, this work advances understanding of critical variables that govern reliability of additively manufactured precipitation-hardened stainless steels in demanding applications.

Brown, Benjamin [Kansas City National Security Cam↗
DOE OSTI2025

In situ irradiation of spent nuclear fuels

To improve the economics of commercial nuclear reactors, nuclear vendors and utilities are seeking approval for increased burnup usage of the existing nuclear fleet. This is especially critical for meeting the clean energy initiative by the United States Government, Department of Energy (DOE). However, one of the key challenges the nuclear industry faces in this regard is that the regions exposed to high burnup and low temperatures exhibit a fine-grained microstructure with large bubbles known as high-burnup structure (HBS) [1]. The formation of HBS has been correlated to the diminished performance of the reactor, as well as fuel fragmentation and pulverization during transient and accidental conditions [2]. Therefore, it is paramount to understand the mechanisms for HBS formation along with its impact on the properties and performance of nuclear fuels. While existing programs, such as Nuclear Energy Advanced Modeling and Simulation (NEAMS) and Advanced Fuel Campaign (AFC) are focusing on evaluating the performance impact of HBS, the physical mechanisms contributing to HBS formation are still not fully understood. In addition, having predictive capabilities and sound understanding of the microstructural evolution of nuclear fuel is essential for accelerated development, qualification, and deployment of new nuclear materials and novel reactor designs for advanced nuclear reactors. There is a lack of consensus among researchers regarding the mechanisms leading to such restructuring observed in HBS. Grain subdivision due to polygonization versus recrystallization, continuous versus discrete recrystallization occurring in tandem or conjunction, etc., have been proposed and debated. In general, it is hypothesized that defect accumulation and dislocation interaction within the grains cause the realignment of dislocations into grain boundaries, leading to the new subgrain formation, which over time transforms into new grains. However, due to the lack of transient data, the importance of fission rate, irradiation, thermal, and stress history of the fuel on the restructuring could not be assessed. In situ microstructural evolution under various irradiation conditions is desired to bridge this gap. Alternatively, phase-field-based models have been developed to capture HBS formation via discrete recrystallization utilizing the classical nucleation approach [3–5]. However, in these models, the grain nucleation criteria are often defined based on empirical relations for burnup and fission gas density leading to dislocation density change. A mechanistic approach to capture the dislocation interaction with the microstructural features leading to grain subdivision is lacking.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Structure–property relationships of reduced graphene oxide membranes intercalated with polycyclic aromatics

Graphene oxide (GO) membranes intercalated with various organic moieties have shown excellent potential for a range of water processing applications. However, microstructure–functional property relationships in these structurally disordered membranes are not well understood. We demonstrate a practical methodology for developing such relationships for GO membranes intercalated with molecular species, with polycyclic aromatic toluidine blue O (TBO) as an example functional intercalant. We use solid-state UV–vis absorbance and fluorescence measurements to quantitatively track the arrangements of TBO in a series of TBO-loaded reduced GO (rGO) membranes. This study reveals the evolution of diverse arrangements including TBO monomers, lateral and stacked dimers, and other aggregates, as a function of overall TBO loading. These microstructures are then correlated to changes in overall properties such as interlayer d-spacings, permeate fluxes, and solute rejections. The characterization of these different intercalant microstructures explains non-intuitive flux and rejection trends, which can circumvent flux and solute rejection trade-offs.

02 PETROLEUM↗
DOE OSTI2024

Hypereutectic Al-Ce-X (X=Mn, Cr, V, Mo, W) alloys fabricated by laser powder-bed fusion

Here, we characterize the microstructures and high-temperature mechanical properties of Al-2Ce and ternary Al-2Ce-1X (at.%) alloys fabricated by laser powder-bed fusion (LPBF), where X = Mn, Cr, V, Mo, and W are slow-diffusing transition metals. All ternary alloys show a hypereutectic microstructure in the as-LPBF state, containing an interconnected network of eutectic Al 11 Ce 3 phases (~10 vol.%) and an additional population of submicron, equiaxed Al 20 CeX 2 primary precipitates (~10 vol.%) which are isomorphous among these five alloys. Similar microstructures are present in arc-melted rods and atomized powders but are coarser due to the slower cooling rates in these processes. The hardness of the as-LPBF ternary Al-Ce-X alloys (1300–1400 MPa) is higher than that of the binary Al-Ce alloy (~1100 MPa) due to the higher volume fraction of strengthening phases. Furthermore, during exposure at 400 °C for up to three months, greater hardness retention is achieved in the ternary Al-Ce-X alloys (65–75%) than in the binary Al-Ce alloy (~55%), which is attributed to the extreme coarsening resistance of the Al 20 CeX 2 precipitates imparted by the very slow-diffusing ternary solute. These coarsening-resistant Al 20 CeX 2 precipitates also substantially improve alloy creep resistance, increasing the threshold stress for dislocation creep at 300°C from ~32 MPa for the binary Al-Ce alloy to ~77–100 MPa for the ternary Al-Ce-X alloys, and at 400°C from <10 MPa for the binary Al-Ce alloy to >40 MPa for the ternary Al-Ce-V alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

ExaCA v2.0: A versatile, scalable, and performance portable cellular automata application for additive manufacturing solidification

The previously established ExaCA software for performance portable alloy grain structure simulation has been updated to better represent the solidification behavior during complex alloy processing conditions, such as those encountered during metal additive manufacturing (AM), and for improved performance and scalability. Here, an extension to the time–temperature history input data format and the core ExaCA algorithm to include an arbitrary number of melting and solidification events yielded improved prediction of texture for various melt pool geometries, expanding the range of AM-relevant conditions that can be accurately simulated. Improved heat transport process simulation coupling, including the creation of large raster datasets from single track time–temperature history data and in-memory coupling with the new, performance portable finite difference code Finch, were also demonstrated in example studies on the effect of multilayer AM microstructure predictions on hatch spacing and cell size, respectively. Additional new features are detailed and demonstrated, including the ability to perform simulations using various interfacial response function forms, execute simulations on state-of-the-art hardware, improved usability through post-processing versatility, and improved strong and weak scaling performance. The performance, physics, and versatility improvements demonstrated here will further enable large-scale studies on AM process–microstructure relationships that were not previously possible. Furthermore, the usability improvements and ability to run coupled AM process–microstructure simulations using the Finch-ExaCA workflow will facilitate broader use of this open-source software by the computational materials community.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Influence of rolling reduction and annealing on recrystallization and grain structure in Ta-2.5W alloys

The microstructural evolution of wrought Tantalum - 2.5 wt% Tungsten (Ta-2.5W) alloys during thermomechanical processing is critical for optimizing their mechanical reliability in demanding applications such as aerospace, chemical processing, and nuclear technology. Despite the widespread use of Ta-W alloys, a comprehensive understanding of how rolling reduction, annealing temperature, and elemental inhomogeneity interact to determine recrystallization behavior and grain refinement remains incomplete. Here, in this study, we systematically investigate the effects of cold rolling and subsequent annealing on the microstructure of Ta-2.5W, with particular attention to grain orientation, stored energy, and elemental banding. Our results demonstrate that higher rolling reduction rates lower the onset and completion temperatures for recrystallization, resulting in finer and more homogeneous grain structures. Electron backscatter diffraction (EBSD) analysis reveals that grains with a 〈111〉 parallel to the plate normal possess higher stored energy and nucleate recrystallization more readily than grains with a 〈001〉 parallel to the plate normal. Elemental mapping shows that tungsten inhomogeneity leads to localized bands of accelerated recrystallization and hardness variation. These findings provide new insights into the mechanisms of microstructure refinement in Ta-2.5W alloys, offering guidance for tailoring processing routes to achieve superior performance in demanding engineering environments.

Annealing↗
DOE OSTI2025

Suppression of irradiation hardening in tungsten-coated ferritic steel for fusion reactor blanket applications

W-coated reduced activation ferritic steels have been developed for use as plasma facing components in fusion reactor blankets, offering excellent sputtering resistance and structural strength. Previous high-temperature coating methods, such as diffusion bonding and brazing, caused interfacial deterioration due to thermal stress from mismatched thermal expansion between W and reduced activation ferritic steel. To address this, underwater explosive welding was introduced as a high-velocity cold process that joins dissimilar materials while maintaining a strong, thin interface without the thermal issues associated with traditional methods. In this study, the effects of neutron irradiation on the hardness and microstructure in W-coated F82H reduced activation ferritic steel (W/F82H) joined by underwater explosive welding are investigated. Following neutron irradiation at 290 °C, irradiation hardening is suppressed in W, F82H, and their interface within the W/F82H material. Furthermore, microstructural observations indicate that the recovery of work hardening and relaxation of elastic strain introduced during coating significantly contribute to the suppression of irradiation hardening in W/F82H. In conclusion, W/F82H exhibits significantly suppressed irradiation hardening compared with those in stand-alone materials. This suppression is explained by residual stress from thermal expansion mismatch and the unique microstructure at the interface. These results provide valuable insights for the development of more durable materials in nuclear fusion applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Laser ablation of high-loading Li-ion battery electrodes improves accessible capacity and cycle life for Behind-the-Meter Storage

Adoption of Behind-the-Meter Storage (BTMS) requires design of batteries that enable high safety, long cycle life, and low cost at the system level. Pairing Li 4 Ti 5 O 12 (LTO) with LiMn 2 O 4 (LMO) achieves targets related to safety and cycle life, but these materials' low energy densities contribute to higher cost at the system scale. Increasing electrode loading is a simple approach to improve energy density, but comes with a trade-off in electrode utilization due to long, tortuous Li + diffusion pathways. Here, laser ablation is used to microstructure (pattern) high-loading electrodes to enhance electrode performance through improved Li + diffusion pathways. Four cell types, comprising combinations of standard or patterned anode and cathode, were prepared to evaluate the effects of laser ablation at each electrode. A rate test shows that patterning electrodes enhances active material utilization at ≳1C rates. Patterning the cathode yields the most benefit, as cells with a patterned cathode demonstrate a ~20% higher accessible capacity than those without at 1.4C. Additionally, 1C capacity retention of cells with patterned cathode (91% through 3000 cycles) is significantly improved over cells with only the anode patterned (64%) and non-patterned electrodes (50%). Characterization of post-mortem cells before and after refreshing their electrolyte suggests that 1C capacity retention is improved by mitigation of electrode "dry-out". We hypothesize that the microstructure acts as a reservoir of additional electrolyte, or a path for gas to escape, so that active material remains wetted throughout long-term cycling, and/or the microstructure may reduce localized, gas-forming overpotentials in the high-loading electrode.

25 ENERGY STORAGE↗
DOE OSTI2025

Additive manufacturing of metal matrix composites

Although Metal matrix composites (MMCs) are superior to most sought-after metallic alloys, their challenging fabricability has limited their widespread use in bulk-form applications. Among the many advanced fabrication techniques, Additive Manufacturing (AM), owing to its unique capabilities to produce near-net shapes, has drawn significant traction in the past two decades, especially for materials that are difficult to process using traditional methods. However, unlike pure metal/alloy systems, MMCs are highly sensitive to the processing conditions prevailing in AM techniques due to factors such as the high melting point of reinforcement particles and the potential for in-situ reactions. Therefore, it may be a while before metal matrix composites are commercially produced via AM. This review will discuss the current state-of-the-art design, fabricability, and performance of various additively manufactured MMCs. A particular focus will be on microstructural evolution and microstructure-property relationships. The most employed AM techniques, such as directed energy deposition, powder bed fusion, binder jetting, sheet lamination, and solid-state friction stir processing, are fundamentally different in terms of thermo-kinetics, forming the perspective for this review. A detailed comparison of microstructural evolution and process parameter optimization, including feedstock preparation methods and the role of machine learning and modeling among the different AM processes, is also presented. Finally, a critical evaluation of emerging AM technologies for MMCs is also provided, highlighting their potential advantages and challenges.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

A mechanistic interpretation of Nelson curves for PVP failures under high temperature hydrogen attack

As an empirically established design criterion, Nelson curves that relate the service temperature and the allowable hydrogen partial pressure have been developed and utilized for more than sixty years in pressure vessels and piping (PVP) safety design. Despite a relatively clear thermodynamic understanding of the high-temperature-hydrogen-attack (HTHA) problem, the detailed fracture process on the microstructural length scales, however, remains elusive, and a quantitative assessment of the PVP lifetime under HTHA from the available creep fracture dataset is still not possible. This work develops a microstructure-informed and micromechanics-based model by incorporating a synergy between hydrogen transport and intergranular-cavity-based fracture process. Based on the available creep lifetime data of C-0.5Mo steels, we are able to calibrate material constitutive parameters, and then conduct nonlinear finite element simulations that reveal a real-time stress-induced hydrogen diffusional transport along grain boundaries, coupled with a microstructure-explicit failure process, from which Nelson curves can be computed. Such failure analyses allow us to delineate two distinct regimes on the Nelson curves, i.e., dislocation-creep-controlled or grain boundary diffusion-assisted cavity growth. More importantly, we found that a small change of the pipe thickness and applied stresses can significantly shift these lifetime curves. However, these two parameters are usually not provided in Nelson curves, thus limiting their usage in material selection and safety design. In conclusion, this discrepancy can clearly be mitigated by extensive parametric studies from our micromechanical modeling/simulation framework.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Recrystallization driven softening and heating rate dependencies of FeCrAl nuclear fuel cladding during accident transients

A refined understanding of FeCrAl cladding behavior during rapid transients is critical for its potential deployment in light-water reactors. Current assessments focus on transient burst testing metrics such as balloon geometry, burst temperature, and hoop stress, often used as proxies for simpler conventional tensile properties. However, directly correlating isothermal tensile and creep data with accident transient scenarios remains a challenge, although it is essential for high fidelity model development. Recent modeling based on tensile tests up to 800 °C, conducted with both immediate loading and a 10-minute soak, showed that immediate loading better predicts experimental burst temperatures, indicating a thermal softening effect. Building upon this observation, the current study connects transient performance, microstructural evolution, and high-temperature tensile properties by leveraging results from C26M claddings burst tests performed at heating rates of 1–50 °C/s and hoop stresses from 25 to 100 MPa. At 25 MPa, rupture temperatures varied by only 6 °C, but at 100 MPa, the difference reached 116 °C, with faster heating yielding higher burst temperatures. Microstructural analysis identified recrystallization as the primary cause of heating rate-dependent softening, eliminating prior cold-working. In-situ thermomechanical data linked ballooning onset to localized instabilities, similar to ultimate tensile strength behavior in conventional tensile tests. High heating rates correlated with immediate loading tensile data, while lower rates matched soaked data. Furthermore, by linking burst performance to microstructural evolution and tensile properties, this work provides a foundation for more accurate modeling of FeCrAl claddings and potentially other Fe-based materials under accident conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

X-ray Micro-Computed Tomography for Structural Analysis of All-Solid-State Battery at Pouch Cell Level

Characterizing the microstructure of all-solid-state batteries (ASSBs) during fabrication and operation is vital for their advancement, particularly as scaling to pouch cell levels introduces challenges in probing large-scale microstructural evolution. This work highlights the potential of synchrotron X-ray micro-computed tomography (sXCT) as a nondestructive, rapid (<30 min), and high-resolution technique for visualizing and quantifying key microstructural features, including overhang, porosity, contact loss, active surface area, and tortuosity, in all-solid-state pouch cells. The large field of view (up to millimeters) of sXCT enables detailed analysis at an industry-relevant scale, bridging the gap between laboratory research and commercial applications. Furthermore, integrating realistic sXCT-derived 3D models into multiphysics simulations could provide insights into chemo-mechanical degradation, particularly at the edges of the pouch cells, offering a pathway for designing robust, high-performance ASSBs. This perspective establishes sXCT as an indispensable tool for advancing both the understanding and the engineering of next-generation energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2024

The effect of neutron irradiation on conventional and nanocrystalline nickel: NSUF Work Package UA-22PN080501 Milestone M3UA-22PN0805012 (Final Report)

Even though nanocrystalline materials (20-100 nm) present an unprecedented potential, scientific knowledge related to the effect of neutron irradiation on the mechanical properties and microstructure is still scarce. Most of the past studies were conducted using ion irradiation which may not have the same effect as neutron irradiation because of the smaller irradiation volume and the higher dose rate. To reach a firm conclusion on the potential of nanocrystalline materials for nuclear reactor applications, extensive study of model metals with different stacking fault energy (SFE) is required to elucidate their behavior in radiation environments. Nanocrystalline copper and nickel are typically chosen because they are commonly used as model FCC metals in studies of radiation effects. Nickel is an FCC metal with a high stacking-fault energy (~125 mJ/m 2 ) compared to copper (~45 mJ/m 2 ). Hence, microcrystalline and nanocrystalline nickel samples were irradiated in the INL’s Advanced Test Reactor (ATR), as a part of FY08 North Carolina State University NSUF Irradiation Experiment# 96 to evaluate the irradiation behavior of these materials. The objective of this FY20 NSUF project# 19122 is to perform PIE (at PNNL NSUF facility) on previously ATR-neutron irradiated (1.2 and 2.6 dpa; 80-89°C) nanocrystalline and microcrystalline nickel samples to investigate the changes in mechanical properties and microstructures and evaluate whether nanocrystalline nickel is relatively more radiation resistant compared to conventional microcrystalline nickel. To perform PIE at PNNL, sixteen neutron irradiated specimens (microcrystalline and nanocrystalline) were transferred from the NSUF Nuclear Fuels and Materials Library at INL. Experimental techniques such as SEM/EBSD, XRD, TEM, Vickers microhardness and tensile testing were employed to characterize the effect of neutron irradiation on the microstructure and mechanical properties of nanocrystalline nickel and compared them with corresponding characteristics of microcrystalline nickel.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Highly Ordered Eutectic Mesostructures via Template‐Directed Solidification within Thermally Engineered Templates

Template-directed self-assembly of solidifying eutectics results in emergence of unique microstructures due to diffusion constraints and thermal gradients imposed by the template. Here, the importance of selecting the template material based on its conductivity to control heat transfer between the template and the solidifying eutectic, and thus the thermal gradients near the solidification front, is demonstrated. Simulations elucidate the relationship between the thermal properties of the eutectic and template and the resultant microstructure. The overarching finding is that templates with low thermal conductivities are generally advantageous for forming highly organized microstructures. When electrochemically porosified silicon pillars (thermal conductivity < 0.3 Wm −1 K −1 ) are used as the template into which an AgCl-KCl eutectic is solidified, 99% of the unit cells in the solidified structure exhibit the same pattern. In contrast, when higher thermal conductivity crystalline silicon pillars (≈100 Wm −1 K −1 ) are utilized, the expected pattern is only present in 50% of the unit cells. The thermally engineered template results in mesostructures with tunable optical properties and reflectances nearly identical to the simulated reflectances of perfect structures, indicating highly ordered patterns are formed over large areas. This work highlights the importance of controlling heat flows in template-directed self-assembly of eutectics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Two–Photon Printing of Base–Catalyzed Resins Involving Thiol–Michael Network Formation

Two-photon printing is accomplished using a photobase generator, 2-(2-nitrophenyl)-propyloxycarbony tetramethyl guanidine, to locally catalyze thiol-ene coupling between multifunctional monomers while mitigating chain-growth polymerization with free radical scavengers. Microstructures printed from base-catalyzed resins exhibit higher resolution (linewidths < 400 nm) and lower print error (<3%) than analogous microstructures printed using a photosensitive, free radical initiator. Further, Raman spectroscopy reveals that resins polymerized using the photobase generator in the presence of free radical scavengers exhibit higher selectivity of thiol-ene coupling over radical polymerization, resulting in stiffer, more uniform polymer networks. The base-catalyzed resins are capable of producing 3D microstructures printed with high accuracy and minimal post-processing defects. As a result, a direct comparison between free radical and base-initiated resins highlights the need for mindful consideration of how chemical reaction pathway influences printability and network end-properties when designing resins.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method

Abstract An interface-modified reproducing kernel particle method (IM-RKPM) is introduced in this work to allow for a direct model construction from image pixels of heterogeneous polycrystalline Li-ion battery microstructures. The interface-modified reproducing kernel (IM-RK) approximation is constructed through scaling of a kernel function by a regularized distance function in conjunction with strategic placement of interface node locations. This leads to RK shape functions with either weak or strong discontinuities across material interfaces, suitable for modeling various interface mechanics. With the placement of a triple junction node and distance-based scaling of kernel functions, the resulting IM-RK shape function also possesses proper discontinuities at the triple junctions. This IM-RK approximation effectively remedies the well-known Gibb’s oscillation in the smooth approximation of discontinuities. Different from the conventional meshfree approaches for interface discontinuities, this IM-RK approach is done without additional degrees of freedom associated with the enrichment functions, and it is formulated with the standard procedures in the RK shape function construction. This work focuses on identifying the accuracy and convergence properties of IM-RKPM for modeling the coupled electro-chemo-mechanical system. A linear patch test is formulated and numerically tested for the electro-chemo-mechanical coupled problem with a Butler–Volmer boundary condition representing the physical conditions in Li-ion battery microstructures. This is followed by verification of the optimal rates of convergence of IM-RKPM for solving the coupled problem with higher order solutions. The image-based modeling of Li-ion battery microstructures in the numerical examples demonstrates the applicability of the proposed method to realistic Li-ion battery materials modeling.

25 ENERGY STORAGE↗