Search NASA⌕ Search

SEARCH · Search NASA

Results for “material response simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model Updates

This work presents the latest improvements to, and investigations performed with, the Fluoride-Cooled High-Temperature Pebble-Bed Reactor reference plant models for the United States Nuclear Regulatory Commission. These models, developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, serve as the foundation for the future development of detailed design evaluation models based on license applications. BlueCRAB is the code suite proposed for non-light-water reactor systems safety analyses, and it incorporates various simulation tools developed by the Nuclear Energy Advance Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of this work is to assess the level of readiness of BlueCRAB for modeling fluoride-cooled high-temperature pebble-bed reactors. We first developed numerical models in BlueCRAB that include the key physics for this technology, ensuring an adequate level of fidelity for modeling the core performance during accident scenarios. This was followed by simulation of transient scenarios, two loss-of-forced-cooling events (one protected and one unprotected), and two control rod withdrawal events (one delayed and one prompt supercritical reactivity insertion). The analysis includes comparisons between the 2-D thermal fluid porous media models in Pronghorn and SAM, comparisons between coupled Pronghorn-Griffin and coupled SAM-Griffin models for two loss-of-forced cooling events and one control rod withdrawal event, and comparisons between SAM single-solve and domain-overlapping approaches for multi-scale thermal fluid coupling. In addition, we performed comparisons between 3-D, 2-D, and 0-D neutronic models for the two control rod withdrawal scenarios with Pronghorn-Griffin. The results show that the BlueCRAB models led to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields such as the neutron flux, power, temperatures, densities, pressures, and velocities are all within the expected ranges, and their distributions can be explained from the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. The simulations demonstrate that there are differences between the various levels of fidelity, and it is advisable to have flexible tools that can cover the breadth and depth of needs that may arise in future technical evaluations. We believe that BlueCRAB’s capabilities represent a significant asset for confirmatory analyses aimed at resolving important safety questions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Dispatch Manager for NEML2 Constitutive Model Calculations Embedded in MOOSE

This report describes the extended capabilities of the NEML2 constitutive modeling library, including a flexible and efficient work dispatching system designed to leverage both CPU and GPU resources. This enhancement addresses one of the primary computational challenges in large-scale simulations: the ability to distribute and execute batches of material model evaluations across heterogeneous computing devices. The new dispatch system introduces a modular set of dispatcher and scheduler classes that coordinate the flow of data and execution between devices. The dispatcher is responsible for efficiently packaging work, managing device-specific memory operations, and synchronizing results. This modularity allows for extensibility, making it straightforward to integrate additional computing backends in the future. From an implementation standpoint, the dispatcher system interfaces seamlessly with NEML2's existing models. They handle device-aware tensor operations, optimize memory transfers, and support asynchronous execution when applicable. This design ensures that batches of material points can be evaluated concurrently, substantially improving throughput compared to previous single-device or serial implementations. These improvements not only enhance the raw performance of NEML2 but also improve its usability in multiscale and high-fidelity simulations, where the simultaneous evaluation of large material point batches is critical. Benchmarks included in the report demonstrate the system’s scalability, highlighting its effectiveness when leveraging modern GPU architectures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Nanoprecipitates to Enhance Radiation Tolerance in High-Entropy Alloys

The growth of advanced energy technologies for power generation is enabled by the design, development, and integration of structural materials that can withstand extreme environments, such as high temperatures, radiation damage, and corrosion. High-entropy alloys (HEAs) are a class of structural materials in which suitable chemical elements in four or more numbers are mixed to typically produce single-phase concentrated solid solution alloys (CSAs). Many of these alloys exhibit good radiation tolerance like limited void swelling and hardening up to relatively medium radiation doses (tens of displacements per atom (dpa)); however, at higher radiation damage levels (>50 dpa), some HEAs suffer from considerable void swelling limiting their near-term acceptance for advanced nuclear reactor concepts. In this study, we developed a HEA containing a high density of Cu-rich nanoprecipitates distributed in the HEA matrix. The Cu-added HEA, NiCoFeCrCu 0.12 , shows excellent void swelling resistance and negligible radiation-induced hardening upon irradiation up to high radiation doses (i.e., higher than 100 dpa). The void swelling resistance of the alloy is measured to be significantly better than NiCoFeCr CSA and austenitic stainless steels. Density functional theory simulations predict lower vacancy and interstitial formation energies at the coherent interfaces between Cu-rich nanoprecipitates and the HEA matrix. The alloy maintained a high sink strength achieved via nanoprecipitates and the coherent interface with the matrix at a high radiation dose (~50 dpa). From our experiments and simulations, the effective recombination of radiation-produced vacancies and interstitials at the coherent interfaces of the nanoprecipitates is suggested to be the critical mechanism responsible for the radiation tolerance of the alloy. Furthermore, the materials design strategy based on incorporating a high density of interfaces can be applied to high-entropy alloy systems to improve their radiation tolerance.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-purity germanium semiconductor modeling in the detector response function toolkit

In this study, we have extended the detector response function toolkit (DRiFT) to provide modeling capabilities of semiconductor sensors. DRiFT provides realistic nuclear instrumentation response by post-processing Monte-Carlo N-particle (MCNP®) radiation transport outputs. MCNP® is capable of modeling radiation transport in complex environments, but has limited detector physics and readout electronics modeling capabilities. Semiconductor detector response can be calculated with a high-fidelity for a flexible range of environments by utilizing MCNP® to simulate radiation interactions inside of detector volumes, and then using DRiFT to model charge transport and signal formation in the semiconductor, as well as the readout electronics. DRiFT models charge transport in the semiconductor, the preamplifier, shaping amplifier, pulse pile-up, and electronic noise to generate detector response. The semiconductor application in DRiFT can model a range of semiconductor materials, shapes, and sizes; and is demonstrated here for a large volume coaxial high-purity germanium (HPGe) detector. Here, we compare detector response functions of a coaxial HPGe detector with measurement of 60 Co, 133 Ba, and 137 Cs at varying count rates, and we conduct a parameter study to demonstrate the effect of changing parameters in the DRiFT simulation. The HPGe detector response function shows excellent agreement with measurements of difference sources with varying dead times and count rates.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Giant electric field-induced second harmonic generation in polar skyrmions

Electric field-induced second harmonic generation allows electrically controlling nonlinear light-matter interactions crucial for emerging integrated photonics applications. Despite its wide presence in materials, the figures-of-merit of electric field-induced second harmonic generation are yet to be elevated to enable novel device functionalities. Here, we show that the polar skyrmions, a topological phase spontaneously formed in PbTiO3/SrTiO3 ferroelectric superlattices, exhibit a high comprehensive electric field-induced second harmonic generation performance. The second-order nonlinear susceptibility and modulation depth, measured under non-resonant 800 nm excitation, reach ~54.2 pm V −1 and ~664% V −1 , respectively, and high response bandwidth (higher than 10 MHz), wide operating temperature range (up to ~400 K) and good fatigue resistance (>10 10 cycles) are also demonstrated. Through combined in-situ experiments and phase-field simulations, we establish the microscopic links between the exotic polarization configuration and field-induced transition paths of the skyrmions and their electric field-induced second harmonic generation response. Our study not only presents a highly competitive thin-film material ready for constructing on-chip devices, but opens up new avenues of utilizing topological polar structures in the fields of photonics and optoelectronics.

Nonlinear optics↗
DOE OSTI2023

A Novel Spatio-Temporal Regime Tracking Method for Impact Simulations

In this proposal, we present a novel method of tracking the rheological regimes activated during impact cratering events that will allow researchers to gain new insights into cratering mechanics. Rheology describes the stress-strain response of rocks to different conditions. Planetary impact cratering events often occur on too large of a scale to be feasibly captured in controlled experiments. Instead, these dynamic events are primarily studied using multi-physics codes equipped with complex material models that enable calculations of the impact event at scale. However, determining which physical processes are required for the problem of interest is challenging. Because the dominant rheological regimes change with space and time during crater formation, it is difficult to link numerical simulations with observable features of craters at the end of the event. The basis of this work is the implementation of numerical flags that track the activation of each rheological regime throughout impact simulations. We demonstrate this with the ‘Rock Model’ implemented in the CTH shock-physics code. We will use this rheology tracking method to ’zoom in’ on a specific region within an event and track the conditions the rock experiences over time. This work will develop community benchmarks to validate and distribute our implemented rheological models. We will focus on improving the melt models used by the planetary impact modeling community by developing an EOS-aware rheological transition from solid to melt. Through analysis of cell and tracer-particle based tracking data, we will study the effects of different rheologies on modeled outcomes, particularly on the volume and distributions of impacts melts. We will also use this method to link observable features with the rheological mechanisms responsible for them. The deliverables (peer-reviewed papers) from the proposed work are (i) tests of the implemented rheologic processes and demonstrations of the tracking flags; (ii) the first calculations of the spatio-temporal evolution of the dominant rheologies during impact cratering events; and (iii) application to delivery of impactor iron during basin-scale impacts.

58 GEOSCIENCES↗
DOE OSTI2023

Investigation of local distortion effects on X-ray absorption of ferroelectric perovskites from first principles simulations

Understanding the role of ferroelectric polarization in modulating the electronic and structural properties of crystals is critical for advancing these materials for overcoming various technological and scientific challenges. However, due to difficulties in performing experimental methods with the required resolution, or in interpreting the results of methods therein, the nanoscale morphology and response of these surfaces to external electric fields has not been properly elaborated. Here, in this work, we investigate the effect of ferroelectric polarization and local distortions in a BaTiO 3 perovskite, using two widely used computational approaches which treat the many-body nature of X-ray excitations using different philosophies, namely the many-body, delta-self-consistent-field determinant (mb-ΔSCF) and the Bethe–Salpeter equation (BSE) approaches. We show that in agreement with our experiments, both approaches consistently predict higher excitations of the main peak in the O–K edge for the surface with upward polarization. However, the mb-ΔSCF approach mostly fails to capture the L 2,3 separations at the Ti–L edge, due to the absence of spin–orbit coupling in Kohn–Sham density functional theory (KS-DFT) at the generalized gradient approximation level. On the other hand, and most promising, we show that application of the GW/BSE approach successfully reproduces the experimental XAS, both the relative peak intensities as well as the L 2,3 separations at the Ti–L edges upon ferroelectric switching. Thus simulated XAS is shown to be a powerful method for capturing the nanoscale structure of complex materials, and we underscore the need for many-body perturbation approaches, with explicit consideration of core-hole and multiplet effects, for capturing the essential physics in these systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Surface enrichment dictates block copolymer orientation

Orientation of block copolymer (BCP) morphology in thin films is critical to applications as nanostructured coatings. Despite being well-studied, the ability to control BCP orientation across all possible block constituents remains challenging. Here, in this study, we deploy coarse-grained molecular dynamics simulations to study diblock copolymer ordering in thin films, focusing on chain makeup, substrate surface energy, and surface tension disparity between the two constituent blocks. We explore the multi-dimensional parameter space of ordering using a machine-learning approach, where an autonomous loop using a Gaussian process (GP) control algorithm iteratively selects high-value simulations to compute. The GP kernel was engineered to capture known symmetries. The trained GP model serves as both a complete map of system response, and a robust means of extracting material knowledge. We demonstrate that the vertical orientation of BCP phases depends on several counter-balancing energetic contributions, including entropic and enthalpic material enrichment at interfaces, distortion of morphological objects through the film depth, and of course interfacial energies. BCP lamellae are found more resistant to these effects, and thus more robustly form vertical orientations across a broad range of conditions; while BCP cylinders are found to be highly sensitive to surface tension disparity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

ELM sputter erosion modeling of a tungsten coated small angle slot divertor in DIII-D

We modeled plasma edge localized mode (ELM) sputter erosion for a Small Angle Slot divertor with a tungsten coated region (SAS-VW), designed for experiments in the DIII-D tokamak, and proposed for use in future advanced tokamaks. The simulations use a free-streaming, 1000 eV, C +6 and D +1 ELM impingement model, with SOLPS-ITER, ITMC-DYN, and REDEP/WBC code packages for background plasma, material response, and erosion/redeposition respectively. The results show ELM'ing plasma gross and net tungsten erosion fluxes of the mixed-material C/W surface peaking at the slot entrance region, and an order of magnitude higher than for non-ELMs. The per-pulse erosion, however, remains low, of order 0.5 nm, due to expected moderate ELM frequencies and duration in DIII-D. Here, the ELMs result in a ~25x higher peak sputtered W current leaving the divertor slot region, towards the core plasma, compared to the ELM-free plasma case. The time-integrated escape current, however, may not significantly affect core plasma high-Z contamination concerns, for a 1% ELM duty factor, but may be an issue for higher frequency ELMs. In general, the modeling results appear favorable for effective testing of the SAS-VW divertor in DIII-D, and extrapolation to innovative divertor designs in future ITER-like and DEMO fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Practical Analysis Framework for Assessment of Printed Circuit Heat Exchangers in High-Temperature Nuclear Service

Abstract Printed circuit heat exchangers (PCHEs) have high thermal efficiency because of the numerous minuscule channels. These minuscule channels result in a high thermal exchange area per unit volume, making PCHE a top contender for an intermediate heat exchanger (IHX) in high-temperature reactors. Thousands of minuscule channels make finite element analysis of the PCHE computationally infeasible. A two-dimensional analysis is usually performed for the PCHE core, which cannot simulate the local channel level responses reasonably because of the absence of global constraint influence. At present, there is no analysis technique available in the ASME Code or literature that is computationally efficient and suitable for engineers to estimate PCHE local responses. A novel but practical two-step analysis framework is proposed for performing PCHE analysis. In the first step, the channeled core is replaced by orthotropic solids with similar stiffness to simulate the global thermomechanical elastic responses of the PCHE. In the second step, local submodel analysis with detailed channel geometry and loading is performed using the elastic-perfectly plastic (EPP) material model. The proposed two-step analysis technique provides a unique capability to estimate the channel corner responses to be used for PCHE performance assessment. This study first developed a methodology for calculating the elastic orthotropic properties of the PCHE core. Next, the two-step analysis is performed for a realistic size PCHE core, and different issues observed in the results are scrutinized and resolved. Finally, a practical finite element analysis framework for PCHEs in high-temperature nuclear service is recommended.

Engineering↗
DOE OSTI2022

Single Photon Detection with On-Chip Number Resolving Capability

Single photon detection (SPD) plays an important role in many forefront areas of fundamental science and advanced engineering applications. In recent years, rapid developments in superconducting quantum computation, quantum key distribution, and quantum sensing call for SPD in the microwave frequency range. We have explored in this LDRD project a new approach to SPD in an effort to provide deterministic photon-number-resolving capability by using topological Josephson junction structures. In this SAND report, we will present results from our experimental studies of microwave response and theoretical simulations of microwave photon number resolving detector in topological Dirac semimetal Cd 3 As 2 . These results are promising for SPD at the microwave frequencies using topological quantum materials.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

The effect of differential mineral shrinkage on crack formation and network geometry

Rock, concrete, and other engineered materials are often composed of several minerals that change volumetrically in response to variations in the moisture content of the local environment. Such differential shrinkage is caused by varying shrinkage rates between mineral compositions during dehydration. Using both 3D X-ray imaging of geo-architected samples and peridynamic (PD) numerical simulations, we show that the spatial distribution of the clay affects the crack network geometry with distributed clay particles yielding the most complex crack networks and percent damage (99.56%), along with a 60% reduction in material strength. We also demonstrate that crack formation, growth, coalescence, and distribution during dehydration, are controlled by the differential shrinkage rates between a highly shrinkable clay and a homogeneous mortar matrix. Sensitivity tests performed with the PD models show a clay shrinkage parameter of 0.4 yields considerable damage, and reductions in the parameter can result in a significant reduction in fracturing and an increase in material strength. Additionally, isolated clay inclusions induced localized fracturing predominantly due to debonding between the clay and matrix. These insights indicate differential shrinkage is a source of potential failure in natural and engineered barriers used to sequester anthropogenic waste.

58 GEOSCIENCES↗
DOE OSTI2025

Quantitative Nonlinear Optical Polarimetry with High Spatial Resolution

Nonlinear optical microscopy such as in the optical second-harmonic generation (SHG) modality has become a popular tool today for probing materials in the physical and biological sciences. While imaging and spectroscopy are widely used in the microscopy mode, nonlinear polarimetry, which can shed light on materials’ symmetry and microstructure, is relatively underdeveloped. This is partly because quantitative analytical modeling of the optical SHG response for anisotropic crystals and films largely assumes low-numerical aperture (NA) focusing of light, where the plane-wave approximation is sufficient. Tight focusing provides unique benefits in revealing out-of-plane polarization responses, which cannot be detected by near-plane-wave illumination at normal incidence. Here, we outline a method for quantitatively analyzing SHG polarimetry measurements obtained under high-NA focusing within a microscope geometry. Experiments and simulations of a variety of standard samples, from single crystals to thin films, are in good agreement, including measured and simulated spatial SHG maps of ferroelectric domains. A solution to the inverse problem is demonstrated, where the spatial distribution of an SHG tensor with unknown tensor coefficient magnitudes is determined by experimentally measured polarimetry. The ability to extract the out-of-plane component of the nonlinear polarization in normal incidence is demonstrated, which can be valuable for high-resolution polarimetry of 2D materials, thin films, heterostructures, and uniaxial crystals with a strong out-of-plane response.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Structure and Flow-Viscosity of Filled-Polymer-Based 3D Printing Ink: Exploration through Coarse-Grained Molecular Dynamics

The addition of nanofiller particles to a polymer matrix has long been known to enhance or modify the composite’s mechanical and rheological properties. However, quantitatively capturing such changes with molecular level simulations remains computationally challenging. Toward that goal, we performed coarse-grained molecular dynamics of a nanocomposite system at a fixed (25 vol %) filler loading under nonspecific, weak polymer–filler interactions representative of a broad class of technologically important materials. We report several interesting results, including: (1) the equilibrium chain-configuration remains Gaussian-like as in an unfilled melt; (2) smaller filler particles display a stronger tendency to cluster; (3) larger fillers act as plasticizers by reducing the entanglement density and accelerating the chain mobility; and (4) fillers enhance the tensile response modulus, with the effect being stronger for larger particles. We also simulate cluster breakup, yielding, and elongational flow under an applied time-linear tensile strain and study the flow viscosity as a function of filler-size and chain-length.

Materials science↗
DOE OSTI2025

Directed self-assembly of chiral liquid crystals into biomimetic bouligand structures in thin film

The Bouligand structure, renowned for its helicoidal arrangement and enhanced mechanical properties, has attracted significant research interest for its ability to impart enhanced strength to intrinsically soft materials. Biomimetic approaches have centered on fibrous structures in bulk materials, but translating this architecture into thin-film regime for miniaturized-wearable devices with programmable functions remains challenging. Here, we direct the self-assembly of cholesteric liquid crystals (CLCs) into hierarchical helical structures using chemically patterned surfaces. Alternating surface anchoring regions align uniform lying-down helices at the nanoscale, guiding a secondary microscale helical structure exhibiting both left- and right-handed twists. This mimetic Bouligand structure in CLCs enables optical modulation under applied field and strain with enhanced mechanical response. Simulations reveal the structural evolution from initial Bouligand configuration in LC layers to alternating twist helices. This research provides a basis for designing and manufacturing miniaturized or wearable devices with nanometer-scale precision in regulating electro-optical and mechanical properties. Bouligand structures, which offer strength in natural materials, are of interest but difficult to obtain. Here, the authors report the development of such structures by directed self-assembly of cholesteric single crystals into hierarchical helical structures, with the secondary structure having right and left-handed twists.

Bouligand structure↗
DOE OSTI2020

Exploring the relationship between solvent-assisted ball milling, particle size, and sintering temperature in garnet-type solid electrolytes

Garnet-type solid electrolytes, such as Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 (LLZTO), are promising materials for solid-state batteries, but processing remains a challenge, in part due to the high sintering temperature required for densification. Additionally, this temperature can be lowered by decreasing the initial particle size via solvent-assisted ball milling, but the relationship between solvent choice, particle properties, sintering behavior, and ionic conductivity is not well understood. In this work, we systematically explore these parameters, showing that milling in commonly used protic solvents, such as alcohols, effectively decreases the particle size but results in lithium loss (through Li + /H + exchange) that leads to poor sintering. By contrast, milling in aprotic solvents with surfactant reduces the particle size to ~220 nm without lithium loss, enabling the fabrication of dense samples (5.1 g/cm 3 ) with good ionic conductivity (0.43 mS/cm at 25 °C) at a lower sintering temperature (1000 °C). We compare ionic conductivities and activation energies for samples prepared with different particle sizes and sintering temperatures and use multiphase-field simulations to identify the mass transport and microstructural mechanisms responsible for the observed sintering dependence on particle size. These results further clarify the relationship between processing parameters and performance and represent important progress toward overcoming fabrication challenges for these materials.

25 ENERGY STORAGE↗
DOE OSTI2022

Monitoring of mechanical changes in a pipe assembly with complex geometry using multi-mode acoustic signals

Guided wave ultrasonic testing has been widely used for nondestructive and long-range evaluation of diverse mechanical structures, such as pipes and plates. In general, a large number of guided wave modes exist even in the simplest geometries, and they exhibit highly dispersive characteristics. Thus, sophisticated excitation methods and elaborate modal analyses need to be implemented to excite, acquire, and process well-defined pure modes as part of such a measurement method. If the interrogated structure has complex geometry, for example, with cross-sectional variations or multiple joints, the utilization of pure modes becomes significantly more difficult. To overcome this difficulty, this paper proposes the use of multi-mode, multi-frequency guided waves for long-range inspection of structures with arbitrarily complex geometries. Furthermore, using broadband chirp signals and multi-mode time–frequency analysis in a prototypical complex pipe assembly that includes multiple pipe diameters, elbows, and flange connections, we experimentally show sensitive detection of mechanical changes (both material addition and loss). In addition, we demonstrate that the complexity of the structure can be modeled in the digital domain with sufficient detail to allow accurate numerical simulation of the acoustic response using this method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Damage Modeling of Solid Oxide Fuel Cells Accounting for Redox Effects

This work applies a multiscale mechanistic damage model developed for brittle ceramics and implemented in commercial finite element (FE) packages via user subroutines to study progressive damage in solid oxide fuel cells (SOFC) subjected to thermomechanical loading under normal operating and shutdown conditions including redox effects. The damage model captures the micromechanics of stiffness reduction due to material porosity change and microcracking and integrates the as-obtained stiffness reduction law into a continuum damage mechanics (CDM) formulation for the evolution of microcracks up to fracture. The volumetric “swelling” that occurs during redox is treated in constitutive modeling similarly to thermal expansion, but swelling strains are irreversible. Furthermore, this damage model was first validated through predictions of strength and stress-strain response for the SOFC electrode materials. Next, it has been applied to predict the potential for degradation in a generic planar SOFC stack with large active area cells. Multicell stack models were simulated in both co-flow and counter-flow configurations. In addition, a constant temperature redox cycle was also simulated to capture overall cell electrode damage due to volumetric swelling of the nickel (Ni)-based anode in the anode-supported cells.

25 ENERGY STORAGE↗