Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion Couple Experiments”

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

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Evaluation of Engineered Barrier Systems FY20 Report

This report describes research and development (R&D) activities conducted during fiscal year 2020 (FY20) specifically related to the Engineered Barrier System (EBS) R&D Work Package in the Spent Fuel and Waste Science and Technology (SFWST) Campaign supported by the United States (U.S.) Department of Energy (DOE). The R&D activities focus on understanding EBS component evolution and interactions within the EBS, as well as interactions between the host media and the EBS. A primary goal is to advance the development of process models that can be implemented directly within the Generic Disposal System Analysis (GDSA) platform or that can contribute to the safety case in some manner such as building confidence, providing further insight into the processes being modeled, establishing better constraints on barrier performance, etc. The FY20 EBS activities involved not only modeling and analysis work, but experimental work as well. Despite delays to some planned activities due to COVID-19 precautions, progress was made during FY20 in multiple research areas and documented in this report as follows: (1) EBS Task Force: Task 9/FEBEX Modeling Final Report: Thermo-Hydrological Modeling with PFLOTRAN, (2) preliminary sensitivity analysis for the FEBEX in-situ heater test, (3) cement-carbonate rock interaction under saturated conditions: from laboratory to modeling, (4) hydrothermal experiments, (5) progress on investigating the high temperature behavior of the uranyl-carbonate complexes, (6) in-situ and electrochemical work for model validation, (7) investigation of the impact of high temperature on EBS bentonite with THMC modeling, (8) sorption and diffusion experiments on bentonite, (9) chemical controls on montmorillonite structure and swelling pressure, (10) microscopic origins of coupled transport processes in bentonite, (11) understanding the THMC evolution of bentonite in FEBEX-DP—coupled THMC modeling, (12) modeling in support of HotBENT, an experiment studying the effects of high temperatures on clay buffers/near-field, and (13) high temperature heating and hydration column test on bentonite.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Developing the Science Basis for Understanding Polymer Encapsulant Degradation Mechanisms: DuraMAT 2.0 Final Project Report

Polymeric encapsulants are essential materials in photovoltaic modules, protecting sensitive electronics from the environment while providing mechanical integrity to the multilayered assembly. However, these polymeric materials are susceptible to degradation processes driven by the ingress of environmental species, ultraviolet radiation, thermal stresses, and mechanical loading. In this study, we employ a combined atomistic simulation and accelerated aging experimental approach to study the molecular-scale mechanisms of encapsulant degradation. Classical molecular dynamics simulations quantify the diffusion of environmental and degradation species through the polymer matrix, producing composition-specific diffusion coefficients. Reactive simulations characterize activation energy barriers and reaction rate constants for key chemical pathways. In parallel, thermal-desorption analyses coupled with mass spectrometry monitor the emergence and concentration profiles of degradation products under controlled stressor conditions. By integrating simulation and experiment, we establish quantitative correlations between polymer composition, species diffusivity, and chemical reactivity. We anticipate that these relations and quantitative values could serve as high-fidelity inputs to reaction-diffusion models, enabling physics-informed lifetime predictions and guiding the design of more durable encapsulant materials for solar energy applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multi-scale time-resolved electron diffraction: A case study in moiré materials

Ultrafast-optical-pump — structural-probe measurements, including ultrafast electron and x-ray scattering, provide direct experimental access to the fundamental timescales of atomic motion, and are thus foundational techniques for studying matter out of equilibrium. High-performance detectors are needed in scattering experiments to obtain maximum scientific value from every probe particle. Here we deploy a hybrid pixel array direct electron detector to perform ultrafast electron diffraction experiments on a WSe 2 /MoSe 2 2D heterobilayer, resolving the weak features of diffuse scattering and moiré superlattice structure without saturating the zero order peak. Enabled by the detector’s high frame rate, we show that a chopping technique provides diffraction difference images with signal-to-noise at the shot noise limit. Finally, we demonstrate that a fast detector frame rate coupled with a high repetition rate probe can provide continuous time resolution from femtoseconds to seconds, enabling us to perform a scanning ultrafast electron diffraction experiment that maps thermal transport in WSe 2 /MoSe 2 and resolves distinct diffusion mechanisms in space and time.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Suppressed incommensurate order in swedenborgite Ca 0.5 Y 0.5 BaCo 4 O 7

Swedenborgite Ca 0.5 Y 0.5 BaCo 4 O 7 (CYBCO) with geometrically frustrated kagome and triangular lattices exhibits a disordered ground state in spite of strong antiferromagnetic couplings between the Co spins. The character of the disordered state has been debated due to ambiguous Co states in the triangular layers. Here, we perform single-crystal diffuse neutron scattering experiments to fully characterize the short-range spin correlations in CYBCO. Through reverse Monte Carlo and self-consistent Gaussian approximation analyses, we confirm that the Co ions in both the kagome and triangular layers are magnetic, and the interlayer couplings, together with the lattice distortion, promote an incommensurate magnetic order that is suppressed possibly by exchange disorder. Our work clarifies the short-range spin correlations in CYBCO and establishes a general correspondence between unequal interactions and incommensurate spin correlations on the swedenborgite lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Europium diffusion in IG-110 nuclear graphite

Europium diffusion in graphite has been recognized to be of interest in high-temperature gas reactor safety analysis, particularly as an indicator of strontium diffusion. However, no measurements of europium diffusion coefficients have been reported in the literature. In this work, the effective diffusion coefficient of europium was measured using a time-release method. Natural europium was loaded into pre-milled unirradiated IG-110 graphite spheres using a pressurized acid digestion vessel. The time-release experiments were performed in the temperature range 1823 K – 1973 K using a SiC diffusion cell connected to an inductively-coupled plasma mass spectrometer (ICP-MS) via a He gas line. Finally, the results of this work are: D Eu, IG-110 = (1.5 x 10 -3 m 2 /s) exp $(\frac{-2.87\mathrm{x}10^{5} J/mol}{RT})$ This effective diffusion coefficient can be used to aid in predictive modelling of europium transport in HTGRs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Large-Scale Molecular Dynamics Simulation of the Dehydration of a Suspension of Smectite Clay Nanoparticles

Fine-grained sediments and sedimentary rocks play important roles in a variety of modern energy technologies from petroleum geology to geological carbon sequestration and radioactive waste management. However, despite their utility and ubiquity, many of their properties remain poorly understood. In particular, the ability to predict the permeability and mechanics of these media remains a persistent fundamental challenge in the geosciences. In the present work, we show how large-scale classical molecular dynamics (MD) simulations can help interpret the properties of fine-grained sedimentary material. All-atom simulations containing 30 discrete clay particles are utilized to understand the evolution of a clay nanoparticle suspension during its progressive dehydration. Microstructural (pore size distribution, tortuosity, anisotropy), thermodynamic (enthalpy and free energy of hydration, anion exclusion), mechanical (total suction), and transport properties (diffusion coefficient tensors of water and sodium) are calculated and compared to the experiment. Overall, our results provide new insight into the coupled chemistry, mechanics, and transport properties of disordered nanoparticle assemblages and shed light upon the important role of water films in controlling these properties.

02 PETROLEUM↗
DOE OSTI2023

Investigation of desorption of hydrogen gas from polymer matrix using thermal desorption analysis and finite element modeling

Polymer materials are used widely as sealing materials and liners in hydrogen storage and transportation infrastructure. During operating conditions, polymers are exposed to high pressure hydrogen gas periodically. When in contact, the hydrogen gas diffuses through the polymer material with ease due to its smaller sized molecules compared to the polymer molecules. In the event of rapid decompression of outside hydrogen gas, the hydrogen molecules within the polymer diffuse out. This process is known as the desorption of polymer and can cause blisters or cracking causing permanent damage to the polymer. Here, we used a combined experimental and modeling approach to investigate the desorption of the hydrogen gas in the polymer. Thermal desorption analysis (TDA) experiments were performed on the polymer sample exposed to the high-pressure hydrogen gas to measure the diffusion coefficient and equilibrium hydrogen concentration. Experiments also provided the quantity of hydrogen within the polymer during the desorption process which followed Fick’s law. However, for higher pressurization values, some anomaly was observed in the quantity of hydrogen remaining in the polymer sample. A continuum mechanics based coupled diffusion–deformation-damage model was deployed to simulate the gas desorption process in the material using the finite element method to understand its effect on the material. It was found that the anomaly in TDA results were due to the blister formation on the polymer surface which resulted in longer hydrogen retention. Simulations also showed that blister formation can be affected by multiple parameters such as the addition of filler particles, and pre-existing cavity size.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Selective One-Dimensional 13 C– 13 C Spin-Diffusion Solid-State Nuclear Magnetic Resonance Methods to Probe Spatial Arrangements in Biopolymers Including Plant Cell Walls, Peptides, and Spider Silk

Two-dimensional (2D) and 3D through-space 13 C– 13 C homonuclear spin-diffusion techniques are powerful solid-state nuclear magnetic resonance (NMR) tools for extracting structural information from 13 C-enriched biomolecules, but necessarily long acquisition times restrict their applications. In this work, we explore the broad utility and underutilized power of a chemical shift-selective one-dimensional (1D) version of a 2D 13 C– 13 C spin-diffusion solid-state NMR technique. The method, which is called 1D dipolar-assisted rotational resonance (DARR) difference, is applied to a variety of biomaterials including lignocellulosic plant cell walls, microcrystalline peptide fMLF, and black widow dragline spider silk. 1D 13 C– 13 C spin-diffusion methods described here apply in select cases in which the 1D 13 C solid-state NMR spectrum displays chemical shift-resolved moieties. Furthermore, this is analogous to the selective 1D nuclear Overhauser effect spectroscopy (NOESY) experiment utilized in liquid-state NMR as a faster (1D instead of 2D) and often less ambiguous (direct sampling of the time domain data, coupled with increased signal averaging) alternative to 2D NOESY. Selective 1D 13 C– 13 C spin-diffusion methods are more time-efficient than their 2D counterparts such as proton-driven spin diffusion (PDSD) and dipolar-assisted rotational resonance. The additional time gained enables measurements of 13 C– 13 C spin-diffusion buildup curves and extraction of spin-diffusion time constants TSD, yielding detailed structural information. Specifically, selective 1D DARR difference buildup curves applied to 13 C-enriched hybrid poplar woody stems confirm strong spatial interaction between lignin and acetylated xylan polymers within poplar plant secondary cell walls, and an interpolymer distance of ~0.45–0.5 nm was estimated. Additionally, Tyr/Gly long-range correlations were observed on isotopically enriched black widow spider dragline silks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Unraveling Interdiffusion Phenomena and the Role of Nanoscale Diffusion Barriers in the Copper–Gold System

Diffusion is one of the most fundamental concepts in materials science, playing a pivotal role in materials synthesis, forming, and degradation. Of particular importance is solid state interdiffusion of metals which defines the usable parameter space for material combinations in the form of alloys. This parameter space can be explored on the macroscopic scale by using diffusion couples. However, this method reaches its limit when going to low temperatures, small scales, and when testing ultrathin diffusion barriers. Therefore, this work transfers the principle of the diffusion couples to small scales by using core–shell nanowires and in situ heating. This allows us to delve into the interdiffusion dynamics of copper and gold, revealing the interplay between diffusion and the disorder–order phase transition. Our in situ TEM experiments in combination with chemical mapping reveal the interdiffusion coefficients of Cu and Au at low temperatures and highlight the impact of ordering processes on the diffusion behavior. The formation of ordered domains within the solid-solution is examined using high-resolution imaging and nanodiffraction including strain mapping. In addition, we examine the effectiveness of ultrathin Al 2 O 3 barrier layers to control interdiffusion of the diffusion couple. Our findings indicate that a 5 nm thick layer serves as an efficient diffusion barrier. Furthermore, this research provides valuable insights into the interdiffusion behavior of Cu and Au on the nanoscale, offering potential applications in the development of miniaturized integrated circuits and nanodevices.

alloys↗
DOE OSTI2021

Particle Size Distribution Dynamics Can Help Constrain the Phase State of Secondary Organic Aerosol

Particle phase state is a property of atmospheric aerosols that has important implications for the formation, evolution, and gas/particle partitioning of secondary organic aerosol (SOA). In this work, we use a size-resolved chemistry and microphysics model (SOM-TOMAS), updated to include an explicit treatment of particle phase state, to constrain the bulk diffusion coefficient (Db) of SOA produced from ??- pinene ozonolysis. By leveraging data from laboratory experiments performed in the absence of a seed and under dry conditions, we find that the Db for SOA can be constrained (1-7 ×10-15 cm2 s-1 in these experiments) by simultaneously reproducing the time-varying SOA mass concentrations and the evolution of the particle size distribution. Another version of our model that used the predicted SOA composition to calculate the glass transition temperature, viscosity, and, ultimately, Db (~10-15 cm2 s-1) of the SOA was able to reproduce the mass and size distribution measurements when we included oligomer formation (oligomers accounted for about a fifth of the SOA mass). Our work highlights the potential of a size resolved SOA model to constrain the particle phase state of SOA by utilizing historical measurements of the evolution of the particle size distribution.

He, Yicong↗
DOE OSTI2022

Integration of Electromagnetic Geophysics Forward Simulation in Coupled Flow and Geomechanics for Monitoring a Gas Hydrate Deposit Located in the Ulleung Basin, East Sea, Korea

We investigate the feasibility of electromagnetic (EM) geophysics methods to detect the dissociation of gas hydrate specifically from a gas hydrate deposit located in the Ulleung Basin, East Sea, Korea via an integrated flow-geomechanics-EM geophysics simulation. To this end, coupled flow and geomechanics simulation is first performed with the multiple porosity model employed, where a mixed formulation with the finite volume (FV) and finite element (FE) methods are taken for the flow and geomechanics, respectively. From the saturation and porosity fields obtained from the coupled flow and geomechanics, the electrical conductivity model is established for the EM simulation. Solving the partial differential equation of electrical diffusion which is linearized using the 3D finite element method (FEM), the EM fields are then computed. For numerical experiments, particularly two approaches in the configuration for the EM methods are compared in this contribution: the surface-to-surface and the surface-to-borehole methods. When the surface-to-surface EM method is employed, the EM is found to be less sensitive, implying low detectability. Especially for the short term of production, the low detectability is attributed to the similarity of electrical resistivity between the dissociated gas (CH4) and hydrate as well as the specific dissociation pattern within the intercalated composites of the field. On the other hand, when the surface-to-borehole EM method is employed, its sensitivity to capture the produced gas flow is improved, confirming its detectability in monitoring gas flow. Hence, the EM geophysics simulation integrated with coupled flow and geomechanics can be a potential tool for monitoring gas hydrate deposits.

depressurization↗
DOE OSTI2025

Unveiling the Origin of Morphological Instability in Topologically Complex Electrocatalytic Nanostructures

Coarsening and degradation phenomena in metals have largely focused on thermally driven processes, such as bulk and surface diffusion. However, dramatic coarsening has been reported in high-surface-area, nanometer-sized Pt-based catalysts during potential cycling in an electrolyte at room temperature─a temperature too low for the process to be explained purely by surface mobility values measured in both vacuum and electrolytes (∼10 –22 and ∼10 –18 cm 2 /s, respectively). This morphological evolution must be due to a different mechanism for mass transport that is sensitive to electrochemical conditions (e.g., electrolyte composition, potential limits, and scan rate). However, there have been no notable studies of electrochemically induced coarsening in nanometer-sized electrocatalysts. Here, we unveil the origins of coarsening in an electrolyte through coupled in situ experiments and atomistic kinetic Monte Carlo (kMC) simulations. Our work demonstrates electrochemical coarsening is driven by two concurrent mechanisms that can be explained at the atomistic level: (i) dissolution/redeposition during the reduction of an oxidized species and (ii) rapid surface diffusion of undercoordinated atoms.

gold↗
DOE OSTI2020

Identification of the molecular-weight growth reaction network in counterflow flames of the C 3 H 4 isomers allene and propyne

The reaction networks responsible for aromatics formation in counterflow flames of the C 3 H 4 isomers allene and propyne are identified through a combined experimental and modeling study. Mole fraction profiles of near-atmospheric pressure (933 mbar) diffusion flames fueled by the C 3 H 4 isomers are analyzed by means of a newly assembled, chemically detailed kinetic mechanism. The experiment consists of a counterflow burner system that is coupled to a high-resolution time-of-flight molecular-beam mass spectrometer with single-photon ionization via synchrotron-generated vacuum-ultraviolet photons. Flame-sampled, mass-specific photoionization efficiency curves are used to identify the presence of aliphatically substituted aromatic species in addition to the commonly considered pericondensed ring structures. The new mechanism describes the formation and growth of aromatics through repetitive sequences of radical–radical and radical–molecule reactions that include C 1 6 intermediates. Higher concentrations of aromatic species are observed in the allene flame and the new mechanism captures the observed experimental trends very accurately. The results indicate the importance of the aliphatically substituted aromatics and of ring-enlargement reactions for the growth reactions. According to the model simulations, radical+radical recombination and PAH-radical+molecule reactions play an important role in PAH growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Electrochemical properties of poly(ethylene oxide) electrolytes above the entanglement threshold

The ion transport in electrolytes depends on three transport coefficients, conductivity (κ), salt diffusion coefficient (D), and the cation transference number with respect to the solvent velocity ($t_+^0$), and the thermodynamic factor ($T_f$). Current methods for determining these parameters involve four separate experiments, and the coupled nature of the equations used to determine them generally results in large experimental uncertainty. We present data obtained from 64 independent polymer electrolytes comprising poly(ethylene oxide) (PEO) and lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) salt. The molecular weights of PEO ranged from 5 to 275 kg mol -1 ; these samples are all above the entanglement threshold. We minimize the experimental uncertainty in transport and thermodynamic measurements by exploiting the fact that ion transport in entangled polymer electrolytes should be independent of molecular weight. The dependence of $κ, D, t_+^0$, and $T_f$ as a function of salt concentration in the range 0.035 ≤ r ≤ 0.30 are presented with a 95% confidence interval, where r is the molar ratio of lithium ions to ethylene oxide monomer units. While κ, D, and $T_f$ are all positive as required by thermodynamic constraints, there is no constraint on the sign of $t_+^0$. We find that $t_+^0$ is negative in the salt concentration range of 0.093 ≤ r ≤ 0.189.

25 ENERGY STORAGE↗
DOE OSTI2022

CEFR simulation using diffusion code system RAST-F

In this study, the RAST-F code system, which is based on the two-step approach that couples a multi-group cross-section generation Monte-Carlo (MC) code and a multi-group nodal diffusion solver, was used for the neutronic simulation of the CEFR start-up experiments. The numerical solution of the RAST-F system was verified against the full core MC solution MCS at all control rods fully inserted and withdrawn states. The RAST-F solution of the selected experimental simulations was compared against the measurement data. A good agreement between RAST-F and MCS solutions was observed with less than 120 pcm discrepancies and 1.2% root-mean-square error in terms of k{sub eff} and power distribution, respectively. Meanwhile, the RAST-F result agreed well with the experimental data within two-sigma of experimental uncertainty. The good agreement of these results indicates that RAST-F can be used to neutronic steady-state simulations for small core-size SFR, which was challenged to deterministic code systems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Magnetic excitations affected by spin-lattice coupling in the S=3/2 triangular lattice antiferromagnet Ag 2 CrO 2

Ag 2 CrO 2 is an S=3/2 triangular lattice antiferromagnet in which long-range magnetic order with a fivefold magnetic unit cell appears below the transition temperature (T N ) of 24 K. The long-range magnetic order is accompanied by a structural transition. Inelastic neutron scattering experiments were performed to investigate the magnetic interactions, anisotropy, and correlations in this material. Above T N , diffuse scattering originating from the short-ranged magnetic correlations consistent with a 120° structure was observed. Here, the magnetic excitations below T N were reproduced reasonably well using a linear spin-wave model with further-neighbor interactions and easy-axis anisotropy, which are predicted to be due to a previously proposed spin-lattice coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Computational and Experimental Mechanistic Insights into the Ethanol-to-Butanol Upgrading Reaction over MgO

The mechanism of ethanol upgrading to higher products is still under debate, especially regarding intermediate species and hydrogenation and dehydrogenation steps. In this work, we conducted a combined theoretical and experimental approach to contribute to this discussion. For such, detailed electronic structure density functional theory calculations (aiming at probing density of states, infrared spectra, geometric parameters, charge densities, and reaction energetics) and diffuse reflectance infrared Fourier transform spectroscopy experiments were carried out revealing the relevance of an appropriate combination of reactive surface sites to support the formation of several intermediates that are formed in the C-C coupling over MgO. The roles of Mg and O sites were also studied under an electronic perspective and different geometrical arrangements. We found that a kink configuration was the most adequate for ethanol to 1-butanol upgrading. Our calculations also gave us arguments to propose distinct reaction routes, whose mutual predominance would depend upon reaction temperature. At temperatures up to 573 K, the so-called β-route, which goes through scission of a Cβ-H bond and formation of an oxametallacycle-like intermediate, would dominate the coupling, whereas at higher temperatures, up to 673 K, a more usual Guerbet mechanism, via an aldol coupling step and then consecutive hydrogenations, would be expected. Here, the theoretical conclusions were followed by a careful experimental strategy using sequential experimental planning techniques in order to estimate accurate parameters with the lowest possible experimental load. Information from these different sources were coupled to develop a mathematical model for the rate of the ethanol upgrading reaction, using a Langmuir-Hinshelwood-Hougen-Watson approach. The developed and statistically validated model adequately described the experimental data at 673 K and 1.1 bar total pressure for ethanol partial pressures in the range from 0 to 20 kPa.

09 BIOMASS FUELS↗
DOE OSTI2025

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗