Search NASA⌕ Search

SEARCH · Search NASA

Results for “species diffusion”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The influence of LiH and TiH 2 on hydrogen storage in MgB 2 II. XPS study of surface and near-surface phenomena

We report that Mg(BH 4 ) 2 is a promising solid-state hydrogen storage material, releasing 14.9 wt% hydrogen upon conversion to MgB 2 . The rehydrogenation of MgB 2 is particularly challenging, requiring prolonged exposure to high pressures of hydrogen at high temperature. Here we report an XPS study probing the influence of LiH and TiH 2 on the hydrogen storage properties of MgB 2 in the surface and near-surface regions, as a complementary investigation to a preceding study of the bulk properties. Surface and near-surface properties are important considerations for nanoscale and bulk hydrogen storage materials. If there are reactions occurring at the surface that modify the chemical composition in the near-surface region, species diffusion can alter the chemical composition even deep into the bulk of the material. For LiH/MgB 2 , metastable LiH–B and LiH–Mg species are produced that are more reactive than Bulk MgB 2 . With prolonged glovebox storage, the LiH/MgB 2 material shows increased reactivity towards O and C and enriched levels of Li and B in the near-surface region. In addition, Li induces the growth of Li 2 CO 3 in the surface and near surface regions. Exposing LiH/MgB 2 to hydrogen at 700 bar and 280 °C for 24 h produces borohydride at a temperature 100 °C below the threshold for bulk MgB 2 hydrogenation. In a specifically surface process with macroscopic implications, the hydrogenation conditions also cause Li 2 CO 3 to react with boron hydroxide in the sample to form a Li-deficient glassy lithium borate melt at the interfaces of the particles, bonding them together. Subsequent heating to 380 °C dehydrogenates the borohydride and eliminates the Li-deficient glassy lithium borate. The LiH/MgB 2 material is not reversible because desorption does not lead back to LiH/MgB 2 , but rather to elemental B and Mg metal in the near-surface region. In contrast to LiH, TiH 2 does not react with MgB 2 , despite the favorable thermodynamics for destabilization via TiB 2 formation. Furthermore, high pressure hydrogenation yields only unreacted TiH 2 and MgB 2 in the surface and near-surface regions. Thus, added TiH 2 provides no benefit to MgB 2 hydrogenation, in agreement with the findings of the preceding bulk study.

08 HYDROGEN↗
NASA NTRS1985

Study of the kinetics and mechanism of the thermal nitridation of SiO2

X-ray photoelectron spectroscopy (XPS) has been used to study the nitridation time and temperature dependence of the nitrogen distribution in thermally nitrided SiO2 films. The XPS data show that the maximum nitrogen concentration near the (SiO(x)N(y)/Si interface is initially at the interface, but moves 20-25 A away from the interface with increasing nitridation time. Computer modeling of the kinetic processes involved is carried out and reveals a mechanism in which diffusing species, initially consisting primarily of nitrogen, react with the substrate, followed by formation of the oxygen-rich oxynitride due to reaction of the diffusing oxygen displaced by the slower nitridation of the SiO2. The data are consistent with this mechanism provided the influence of the interfacial strain on the nitridation and oxidation kinetics is explicitly accounted for.

Vasquez, R. P.↗
DOE OSTI2026

Diffusion of NiH on Ni(111) and of Ni with Adsorbed H on Ni(111) in the Context of Ni Coarsening in Solid Oxide Cells

In the Ni-based hydrogen electrode of a solid oxide cell (SOC), Ni coarsening is an important degradation mechanism and could be enhanced by NiH diffusing on the Ni particle surfaces. Here, in this work, the average lifetime and diffusion distance of NiH on Ni(111) are computed using density-functional theory and kinetic Monte Carlo methods. It is found that NiH is extremely short-lived and, thus, cannot promote coarsening in the SOC. Also, the diffusion of Ni on Ni(111) is shown to be at most slightly accelerated by H along NiH dissociation paths involving the movement of Ni with a nearby H. Based on this result, coarsening is not likely to be dramatically accelerated by the diffusion of Ni with H on Ni(111). However, support is provided for the experimental procedure of measuring the product of surface coverage and single-species diffusivity of Ni on Ni(111) in a hydrogen-rich atmosphere.

Ni coarsening↗
DOE OSTI2021

Microstructure-dependent rate theory model of defect segregation and phase stability in irradiated polycrystalline LiAlO 2

We report gamma lithium aluminate (LiAlO 2 ) is a breeder material for tritium and is one of key components in a tritium-producing burnable absorber rod (TPBAR). Dissolution and precipitation of second phases such as LiAl 5 O 8 and voids are observed in irradiated LiAlO 2 . Such microstructure changes cause the degradation of thermomechanical properties of LiAlO 2 and affect tritium retention and release kinetics, and hence, the TPBAR performance. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed for investigating the accumulation of species and phase stability in polycrystalline LiAlO 2 structures under irradiation. Three sublattices (i.e. [Li, Al, V] I [O, V o ] II [Li i , Al i , O i , V i ] III ), and concentrations of six diffusive species (i.e. Li; vacancy of Li or Al at [Li, Al, V] I sublattice, O vacancy at [O, V o ] II sublattice, and Li, Al and O interstitials at [Li i , Al i , O i , V i ] III interstitial sublattices; are used to describe spatial and temporal distributions of defects and chemistry. Microstructure-dependent thermodynamic and kinetic properties including the generation, reaction, and chemical potentials of defects and defect mobility are taken into account in the model. The parametric studies demonstrated the capability of the developed RIS model to assess the effect of thermodynamic and kinetic properties of defects on the segregation and depletion of species in polycrystalline structures and to explain the phase stability observed in irradiated LiAlO 2 samples. The developed RIS model will be extended to study the precipitation of LiAl 5 O 8 and voids and tritium retention by integrating the phase-field method.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding of Low‐Porosity Sulfur Electrode for High‐Energy Lithium–Sulfur Batteries

Abstract The lithium–sulfur (Li–S) battery is a promising technology for large‐scale energy storage and vehicle electrification due to its high theoretical energy density and low cost. Reducing the sulfur cathode porosity has been identified recently as a viable strategy for improving the cell practical energy density and minimizing pore‐filling electrolytes to extend cell life at lean electrolyte conditions. Direct use of a low‐porosity cathode for Li–S battery results in poor electrode wetting, nonuniform electrode reactions, and thus early cell failure. To understand and mitigate the barriers associated with the use of low‐porosity electrodes, multiscale modeling is performed to predict electrode wetting, electrolyte diffusion, and their impacts on sulfur reactions in Li–S cells by explicitly considering the electrode wettability impacts and electrode morphologies. The study elucidates the critical impact of low tortuosity and large channel pore design for promoting electrode wetting and species diffusion. It is suggested that the secondary particle size should be comparable with the electrode thickness to effectively promote electrolyte wettability and sulfur reactivity. This study provides new insights into the low‐porosity electrode material and designs and is expected to accelerate the development of practical high‐energy Li–S batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Cluster dynamics simulations of tritium and helium diffusion in lithium ceramics

Tritium (T) and He diffusion in LiAlO 2 and LiAl 5 O 8 phases influences the performance of tritium producing burnable absorber rods (TPBARs) by affecting the gas release, swelling and thermal conductivity of Li-bearing ceramic pellets. Frenkel pair defects and clusters created by irradiation can attract T and He interstitials and form clusters of the type He i x Li, He i x Al, He i x O, T i x Li, T i x Al, and T i x O, 1 ≤ x ≤ 4 in a Li, Al or O vacancy site (notation denotes x He or T atoms in a 1 Li, 1 Al or 1 O vacant site). The concentration and mobility of each of these clusters collectively contribute to the diffusion of the He and T gases in LiAlO 2 and LiAl 5 O 8 . Here, in this work, free energy cluster dynamics simulations implemented in the Centipede code, are used to obtain the concentration and diffusivities of these clusters which are then used to calculate the total diffusivity of T and He gases in LiAlO 2 and LiAl 5 O 8 . The results show that diffusivity of T is at least one order of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 whereas He diffusion is 2–13 orders of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 . There is a higher concentration of highly diffusive species (T interstitials and T i 03 Li for the case of tritium and He i 01 Li, He i 02 Li, and He i 03 Li for the case of He) in LiAlO 2 than in LiAl 5 O 8 which increase the total diffusion of T and He in LiAlO 2 .

36 MATERIALS SCIENCE↗
NASA NTRS1979

Distribution of energetic positive ion species above a diffuse midnight aurora

The present paper deals with species-identifying distribution function measurements of auroral primary particles, made during a magnetically quiet presubstorm period above a hydrogen-associated diffuse aurora. Only ions identified as H(+) and He(++) were detected. In the mass spectrum data, the He(++) was not clearly above background. At energies between 2.5 and 12 keV, the He(++)/H(+) intensity ratio had an upper limit of 2 to 4 percent. The same upper limit applies to all other ions, such as He(+) and O(+). Though this would suggest a solar wind source for these ions, an admixture including an appreciable fraction of polar wind protons is not precluded. This situation contrasts sharply with a number of recently reported observations of large intensities of precipitating O(+) ions during magnetic storms, and may be characteristic of undisturbed periods.

Moore, T. E.↗
DOE OSTI2020

A direct numerical simulation of Jet A flame kernel quenching

The safe operation of aeronautical engines requires an understanding of flame ignition, propagation and extinction. In this study, direct numerical simulations are performed using a 29 species reduced chemical mechanism for jet fuel surrogate Jet A to understand the flame quenching process. Here, initially laminar spherical flames of varying sizes and equivalence ratios are subject to an identical periodic domain of decaying and isotropic high intensity turbulence with a turbulent Reynolds number of 2400. All cases become quenched, except for the larger kernel with lower Karlovitz number. An analysis of the flame structure shows broadened preheat zone, flame shortening on the product side, differential species diffusion and partial fuel pyrolysis in the fresh mixture. Two extinction mechanisms are identified arising from flame shortening and high flame stretch. Flame shortening occurs due to turbulence-chemistry interactions that resemble the flame–flame interaction in a laminar counterflow reactant-to-reactant configuration, which contorts and breaks up the ignition kernel. Flame stretch is a local effect that attenuates the heat release rate and causes the flame to retreat towards the product mixtures, similar to what has been observed for reactant-to-product laminar counterflow flames. Chemical explosive mode analysis was also performed to quantify the flame structure and local combustion mode. The diffusion–reaction balance in pinched-off flame islands favors extinction of these smaller structures, while auto-ignition modes are observed within the flame kernel after fresh mixture is engulfed and preheated in the product kernel. Statistics of the density-weighted displacement speed conditional on local combustion mode indicates strong correlation between the local extinction mode and negative displacement speed. The local balance between diffusion and reaction ultimately determines the propensity for local extinction in both laminar and turbulent flames, the extent of which has an impact on global flame propagation.

42 ENGINEERING↗
DOE OSTI2021

Parabolic oxidation kinetics of a plutonium alloy at room temperature

Auger electron spectroscopy (AES) was used to characterize the oxide scale formed on an ~0.5 wt. % Ga plutonium metal alloy. Analysis of changes to the characteristic Auger transitions for Pu and O during sputter depth profiles were used in combination with previously measured sputter yields to determine the thickness of the oxide scale and provide insight into its chemistry. From this data, the parabolic kinetics for the Pu-Ga alloy were determined as a function of relative humidity, RH = 14–92%, and oxygen partial pressure, p(O 2 ) = 31–157 Torr, at 25°C. For each relative humidity investigated, the parabolic rate constant was found to be independent of oxygen partial pressure. Furthermore, power-law fits to the relative humidity effect on the rate constant was found to have a power of n = 1.265, indicating the diffusion species which rate-limits the reaction isn’t molecular water, which subsequent dissociates at the oxide-metal interface, or hydroxyl from dissociation at the gas-surface interface.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ab initio molecular dynamics (AIMD) simulations of NaCl, UCl 3 molten salts

Ab initio molecular dynamics (AIMD) simulations are used to calculate select thermophysical and thermodynamic properties of NaCl, UCl and NaCl-UCl 3 molten salts. Following established approaches, the AIMD simulations include a model for Van der Waals interactions. The Langreth & Lundqvist (vDW-DF2), DFT-D3, and density-dependent energy correction (DFT-dDsC) dispersion models are first tested for molten NaCl in order to assess their accuracy for density and heat capacity predictions across a range of temperatures. Based on the NaCl results, the vdW-DF2 and DFT-dDsC methods are extended to the UCl 3 system and compared to available experimental data. Next, mixtures of NaCl-UCl 3 are investigated and analyzed with respect to the deviation from ideal solution behavior. For the DFT-dDsC simulations, density deviates by up to 2% from an ideal mixture, with the maximum occurring close to the eutectic composition. The mixing energy also deviates from an ideal solution and exhibits a minimum of about -0.074 to eV per formula unit, again close to the eutectic composition. Finally, the compressibility and species diffusivity of the pure and mixed salt systems are calculated. The diffusivities are slightly reduced in the mixed compared to the pure systems and the compressibilities loosely follow a linear correlation as a function of the UCl 3 composition. The trends observed for mixing properties are rationalized by correlating them to the coordination chemistry, which emphasizes the importance of maintaining the extended network formed by U and Cl ions as the NaCl concentration increases. The concentration at which breakdown of the extended network occurs, roughly coincides with the minimum of both the mixing energy and the deviation from ideal solution behavior for density.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Detailed chemistry modeling of rotating detonations with dilute n -heptane sprays and preheated air

Utilization of liquid fuels is crucial to enabling commercialization of rotating detonation engines (RDE) in the near future. In this study, Eulerian-Lagrangian simulations are conducted for rotating detonative combustion with dilute n-heptane sprays and preheated air. Two-dimensional flattened configuration is used and a skeletal chemical mechanism with 44 species and 112 elementary reactions for n-heptane combustion is adopted. The flow structure, droplet distribution, and thermochemical parameters in the refill zone are first analyzed. It is shown that the mixture in the refill zone is heterogeneous, including evaporating droplets, vapor, and air. When the total temperature is below 950 K, the average equivalence ratio increases with the total temperature. When it is higher than 950 K, the average equivalence ratio is almost constant. Subsequently, the chemical explosive mode analysis is applied to identify the controlling reactions and dominant combustion modes in the fuel refill zone and reaction fronts. Results demonstrate that the initiation reaction (R104: n-C 7 H 16 + O 2 2-C 7 H 15 + HO 2 ) and low-temperature reaction (R107: RO 2 R'O 2 H) are dominant in the upstream and downstream of the refill zone, respectively. The intermediate species from low-temperature chemistry, R'O 2 H, is found to be important for the chemical explosive mode in the undetonated mixture. The influence of species diffusion and dispersed droplets is further analyzed. Results show that vapor autoignition facilitated by droplet evaporation occurs in the refill zone. Finally, the effects of the air total temperature on the detonation propagation speed and RDE propulsion performance are investigated. Further, it is found that the detonation propagation speed and specific impulse increase with air total temperature. The total pressure ratio first increases and then decreases as the air total temperature increases. Moreover, when the total temperature of the preheated air is above 1,300 K, the effects of low-temperature chemistry are negligible.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Magnetohydrodynamic Control of Interfacial Degradation in Lithium-Ion Batteries for Fast Charging Applications

Interfacial anodic degradation in graphitic materials under fast charging conditions causes severe performance loss and safety hazard in lithium ion batteries. We present a novel method for minimizing the growth of these aging mechanism by application of an external magnetic field. Under magnetic field, paramagnetic lithium ions experience a magnetohydrodynamic force, which rotates the perpendicularly diffusing species and homogenizes the ionic transport. This phenomenon minimizes the overpotential hotspots at the anode/separator interface, consequently reducing SEI growth, lithium plating, and interfacial fracture. In situ electrochemical measurements indicate an improvement in capacity for lithium cobalt oxide/graphite pouch cell (20 mAh) charged from 1–5 C under an applied field of 1.8 kG, with a maximum capacity gain of 22% at 5C. Post-mortem FE-SEM and EDS mapping shows that samples charged with magnetic field have a reduced lithium deposition at 3C and a complete suppression of interfacial fracture at 5C. At 5C, a 24% reduction in the lithium content is observed by performing XPS on the anodic interfacial film. Finally, fast charging performance under variable magnetic field strengths indicate a saturation behavior in capacity at high fields (>2 kG), thereby limiting the field and consequent energy requirements to obtain maximum capacity gain under extreme conditions.

electromagnetism↗
DOE OSTI2024

Anion and Cation Migration at 2D/3D Halide Perovskite Interfaces

Here, this study explores the ionic dynamics in 2D/3Dperovskite solar cells, which are known for their improved efficiency and stability. The focus is on the impact of halide choice in 3D perovskites treated with phenethylammonium halide salts (PEAX, X = Br and I). Our findings reveal that light and heat drive ionic migration in these structures, with PEA + species diffusing into the 3D film in PEABr-treated samples. Mixed-halide 3D perovskites show halide interdiffusion, with bromine migrating to the surface and iodine diffusing into the film. Cathodoluminescence microscopy reveals localized 2Dphases on the 3D perovskite, which become more evenly distributed after thermal treatment. Both PEAX salts enhance the performance of photovoltaic devices. This improvement is attributed to the passivation capabilities of the salts themselves and their respective Ruddlesden–Popper (RP) phases. Annealed PEAI-treated devices show a better balance between efficiency and statistical distribution of photovoltaic parameters.

25 ENERGY STORAGE↗
DOE OSTI2025

Comment on “Brownian motion of droplets induced by thermal noise”

We simulate phase separated fluids using the Cahn-Hillard fluctuating hydrodynamic (CH-FHD) model and measure the statistical properties of capillary waves generated by thermal fluctuations. Our measurements are in good agreement with stochastic lubrication theory and molecular dynamics simulations but differ significantly from recent CH-FHD results by Zhang et al. [Phys. Rev. E 109, 024208 (2024)2470-004510.1103/PhysRevE.109.024208]. Specifically, we find that capillary wave statistics at thermodynamic equilibrium are independent of transport properties, namely viscosity and species diffusion.

Bell, J B↗
DOE OSTI2026

Mass Spectrometer Transient Analysis

This software implements a complete preprocessing pipeline for transient mass spectrometry (MS) data collected during TAP (Temporal Analysis of Products) experiments. It is designed to extract chemically meaningful fluxes from overlapping ion signals by applying a calibrated defragmentation matrix and solving the resulting linear system using non-negative least squares (NNLS) regression. The core script, preprocess_mass_spec.py, performs the following operations: Gain correction: Applies amplifier gain scalars derived from inert-packed calibration pulses to normalize signal intensities across AMUs and acquisition settings. Background subtraction: Removes experiment baselines using user-defined time windows, ensuring compatibility with slow-diffusing species and preventing negative values that would interfere with NNLS. Options to subtract before and after defragmentation. Defragmentation: Constructs a fragmentation matrix A from zeroth moments of calibration pulses (equal molar gas:inert mixtures) and solves Ax=b at each time point, where b is the raw MS signal and x is the estimated species flux. The matrix is normalized to inert signals and accounts for instrument-specific fragmentation behavior. Pulse-mode handling: Supports both averaged and individual pulse modes, enabling statistical treatment of fluxes and calculation of standard deviations. Integration and output: Computes zeroth moments (integrated fluxes) and exports time-resolved and integrated data in CSV format, suitable for downstream kinetic modeling. The software is validated using both virtual TAP simulations (VTAP) and experimental data from propane dehydrogenation (PDH) on CrOx/Al2O3 catalysts. It preserves temporal resolution by applying NNLS point-by-point across the pulse duration (typically 6,000+ time slices per pulse), leveraging the linear superposition principle to reconstruct full flux profiles. The defragmented outputs are compatible with kinetic extraction methods such as the G and Y procedures, which are used to derive rate–concentration relationships from TAP data. The details of these validations are discussed in detail in the supporting manuscript and supporting information. Example data and output files are also included. The methodology is robust to experimental noise and drift, with calibration protocols that account for pulse size effects, MS aging, and inert gas normalization. The software is modular, reproducible, and tailored for high-throughput TAP-MS workflows in catalysis research.

Kristy, Stephen [Idaho National Laboratory (INL), ↗
DOE OSTI2022

Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel

Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Aagesen, Larry K. (ORCID:000000034936676X)↗
DOE OSTI2022

Engineering & Microstructure Scale PIE Report on EBR-II X441A Metallic Fuel Pins for the MORPH Experiment

The objectives of this project are to increase fundamental understanding of irradiation induced metallic U-Pu-Zr fuel behavior and to obtain data needed for the development of irradiation models for metallic fuels in MARMOT. The requested metallic fuel pins are from the X441A experiment irradiated in Experimental Breeder Reactor (EBR)-II and will be provided by the PIs for the duration of the experiment. The purpose of the pins irradiated in the X441A assembly was to vary Zr composition and fuel slug diameter in order to provide data for the metallic fuel performance code “LIFEMETAL” [ANL-IFR-125]. The fuels of interest to this project are metallic fuels with several Zr compositions (in wt.%): U-19Pu-6Zr, U-19Pu-10Zr, and U-19Pu-14Zr, which were irradiated to a peak burnup of approximately 11 at.%. The cladding was the same for all three fuel pins (austenitic stainless steel D9), which allows investigation of fuel-cladding interaction (FCI) phenomena. Varying Zr content in fuel pins enables investigation of the effect of Zr on fuel restructuring and fuel-cladding compatibility. Engineering and microstructure scale PIE activities will be focused on investigation of fundamental aspect of fuel performance such as species diffusion and migration, fission product behavior, and constituent redistribution. Obtained microstructural information will be used as the basis for the development of MARMOT models of U-Pu-Zr fuel performance at the mesoscale. Experimental data obtained through this proposal will be used to provide this fundamental understanding which will serve as the foundation of the development of radiation models for U-Pu-Zr in MARMOT. It will also provide a starting point for the design of new experiments to provide data for the validation of these MARMOT models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Multi-physics Preconditioning for Thermally Activated Batteries

Thermal batteries, also known as molten-salt batteries, are single-use reserve power systems activated by pyrotechnic heat generation, which transitions the solid electrolyte into a molten state. The simulation of these batteries relies on multiphysics modeling to evaluate performance and behavior under various conditions. This paper presents advancements in scalable preconditioning strategies for the Thermally Activated Battery Simulator (TABS) tool, enabling efficient solutions to the coupled electrochemical systems that dominate computational costs in thermal battery simulations. We propose a hierarchical block Gauss-Seidel preconditioner implemented through the Teko package in Trilinos, which effectively addresses the challenges posed by tightly coupled physics, including charge transport, porous flow, and species diffusion. The preconditioner leverages scalable subblock solvers, including smoothed aggregation algebraic multigrid (SA-AMG) methods and domain-decomposition techniques, to achieve robust convergence and parallel scalability. Strong and weak scaling studies demonstrate the solver’s ability to handle problem sizes up to 51.3 million degrees of freedom on 2048 processors, achieving near sub-second setup and solve times for the end-to-end electrochemical solve. These advancements significantly improve the computational efficiency and turnaround time of thermal battery simulations, paving the way for higher-resolution models and enabling the transition from 2D axisymmetric to full 3D simulations.

25 ENERGY STORAGE↗