Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational analysis”

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 1,027 records · Page 57

Innovations in underground hydrogen storage with multiphysics simulations, optimization, and monitoring: A review

Underground Hydrogen Storage (UHS) is a promising solution for large-scale energy storage and a critical component in advancing low-carbon energy system. Ensuring the safety and efficiency of UHS necessitates a comprehensive understanding of multiphysical interactions driven by cyclic pore fluid pressure fluctuations and coupled physicochemical processes. Here, this review examines the key geomechanical responses in UHS, including rock property variations under cyclic loading, fracture evolution and propagation, reservoir stress sensitivity, and fault stability. It also explores the impact of geochemical and microbial reactions on geomechanical characteristics. We provide an in-depth analysis of Thermal-Hydraulic-Mechanical-Chemical (THMC) coupled numerical simulations, highlighting their potential for future multi-scale modeling. Limitations of current machine learning (ML) approaches in addressing UHS challenges are highlighted, emphasizing the need for innovative ML-based methodologies. Operational strategies for hydrogen injection and production are reviewed, focusing on safety, efficiency, and economic viability. The necessity for multi-objective optimization (MOO) to balance storage efficiency, risk mitigation, and cost-effectiveness is also discussed. Current monitoring technologies are evaluated to ensure safe and efficient UHS operations. Finally, this review identifies critical knowledge gaps and underscores the importance of advancing geomechanical understanding under multiphysics-coupling. We highlight the need for ML-driven multiphysics theories, enhanced modeling techniques, and robust optimization strategies to improve UHS performance. This study serves as a comprehensive reference for future research and the large-scale implementation of UHS systems.

25 ENERGY STORAGE↗

Benzo[$\mathcal{a}$]pyrene toxicokinetics in humans following dietary supplementation with 3,3'-diindolylmethane (DIM) or Brussels sprouts

Utilizing the atto-zeptomole sensitivity of UPLC-accelerator mass spectrometry (UPLC-AMS), we previously demonstrated significant first-pass metabolism following escalating (25–250 ng) oral micro-dosing in humans of [ 14 C]-benzo[a]pyrene ([ 14 C]-BaP). The present study examines the potential for supplementation with Brussels sprouts (BS) or 3,3' -diindolylmethane (DIM) to alter plasma levels of [ 14 C]-BaP and metabolites over a 48-h period following micro-dosing with 50 ng (5.4 nCi) [ 14 C]-BaP. Volunteers were dosed with [ 14 C]-BaP following fourteen days on a cruciferous vegetable restricted diet, or the same diet supplemented for seven days with 50 g of BS or 300 mg of BR-DIM ® prior to dosing. BS or DIM reduced total [ 14 C] recovered from plasma by 56–67% relative to nonintervention. Dietary supplementation with DIM markedly increased T max and reduced C max for [ 14 C]-BaP indicative of slower absorption. Further, both dietary treatments significantly reduced C max values of four downstream BaP metabolites, consistent with delaying BaP absorption. Dietary treatments also appeared to reduce the T 1/2 and the plasma AUC( 0,$\infty$ ) for Unknown Metabolite C, indicating some effect in accelerating clearance of this metabolite. Toxicokinetic constants for other metabolites followed the pattern for [ 14 C]-BaP (metabolite profiles remained relatively consistent) and non-compartmental analysis did not indicate other significant alterations. Significant amounts of metabolites in plasma were at the bay region of [ 14 C]-BaP irrespective of treatment. Although the number of subjects and large interindividual variation are limitations of this study, it represents the first human trial showing dietary intervention altering toxicokinetics of a defined dose of a known human carcinogen.

3,3'-diindolylmethane↗

Reconstructed lung doses for the million person study cohort of 26,650 Tennessee Eastman corporation workers employed between 1942 and 1947

Abstract Tennessee Eastman Corporation workers were exposed to uranium dust resulting in high-linear energy transfer (LET) irradiation to lung tissue. In this work, radiation lung doses were reconstructed for 26 650 men and women working at the plant between 1942 and 1947. Site air monitoring data of uranium concentrations and payroll records were used to determine the daily inhaled activities and annualized lung doses. Variations in the activity median aerodynamic diameter of the uranium dust, the solubility of particulate matter in the lungs and the sex-specific breathing rate were investigated as part of a sensitivity analysis. Male and female mean lung doses of 18.9 and 32.7 mGy, respectively, from high-LET alpha irradiation, and there was general agreement with evaluations from previously published epidemiological studies. Annual lung dose estimates and sensitivity analysis for the 26 650 workers in the TEC cohort have been archived on the United States Department of Energy Comprehensive Epidemiologic Data Resource.

61 RADIATION PROTECTION AND DOSIMETRY↗

Investigation of magnetic and electrical transport properties of YMB 4 (M = Cr, Fe, and Co) compounds

We have studied the structural, magnetic, and electrical transport properties of the YMB 4 series of compounds (M= Cr, Fe, and Co). These materials exhibit high stability and possess notable refractory and thermoelectric properties. Furthermore, recent theoretical predictions suggest that some of these compounds may display quantum magnetism and dimer formation among the magnetic ions. We have synthesized almost single phase YMB 4 compounds in the orthorhombic crystal structure (space group Pbam) using arc-melting and annealing. The Rietveld analysis of the room-temperature X-ray diffraction patterns reveals a clear change in the lattice parameter c, which correlates with the atomic radii of Cr, Fe, and Co. The temperature variation of resistivity measured between 2 K and 300 K reveals a metallic electron transport in all three YMB 4 compounds. However, no abrupt change in resistivity due to structural or magnetic phase transitions is observed. Additionally, a weak positive magnetoresistance of 1%–2% has been measured at 2.5 K, the contribution from the metallic impurity cannot be ruled out. The carrier concentration, on the order of 10 21 cm −3 , has been determined using conventional Hall measurements. The thermomagnetic curves recorded between 2 K and 350 K reveal the non-magnetic behavior of all these compounds. However, a small magnetic moment from paramagnetic/ferromagnetic impurity phases is detected at low temperatures. Above room temperature, no magnetic transition associated with the breaking of a dimer is observed.

borides↗

Regional infrasonic observations from surface explosions – influence of atmospheric variations and realistic terrain

SUMMARY A pair of 1-ton, conventional surface explosions were conducted at the Nevada National Security Site in the fall of 2020 producing seismoacoustic signatures observable hundreds of kilometres from the source location. Regional infrasonic observations include tropospheric ducting at large distances to the south, a wide stratospheric waveguide with signals observed more than 700 km to the east, and anomalous arrivals in the stratospheric shadow zone. Notable differences in propagation between the events are identified despite the explosions being conducted just two days apart due to a sharp temporal shift in the tropospheric winds as well as structural changes in the stratospheric winds. Propagation simulations of the two events have been completed using a combination of ray tracing and parabolic equation (PE) methods. Simulations have been conducted to quantify the impact of the temporal variations in the atmosphere as well as the influence of terrain on propagation. Temporal variations in reflection locations are found to produce notable changes in downrange propagation due to spatially varying terrain features. Finite frequency effects modelled by the PE are found to predict ensonification not included in corresponding 2D ray tracing simulations. Notable variations in predicted signal amplitude are found due to focusing by along-path and cross-path terrain gradients; though, the later of these is only modelled using fully 3D ray tracing analysis.

58 GEOSCIENCES↗

Effect of Hydrogen on Steady-State and Transient Combustion Instability Characteristics

This paper examines the effects of steady-state and transient hydrogen enrichment on thermoacoustic instability in a model gas turbine combustor. Combustion instability, a feedback loop between flame heat release rate oscillations and combustor acoustics, is characterized in a swirl-stabilized flame operated at a range of hydrogen–natural gas fuel blends and heat rates. Measurements of combustor chamber pressure fluctuations and CH* chemiluminescence imaging are used to characterize instability at a range of operating conditions. Steady-state tests show that both mixture heat rate and hydrogen content affect system stability. At a given heat rate, higher levels of hydrogen result in unstable combustion. As heat rate increases, instability occurs at lower concentrations of hydrogen in the fuel. Transient operation was tested in two directions—instability onset and decay—and two hydrogen-addition times—a short time of 1 ms and a longer time of 4 s. Results show that instability onset processes, through the transient addition of hydrogen, are highly repeatable regardless of the timescale of hydrogen addition. Certain instability decay processes are less repeatable, resulting in cases that do not fully transition from unstable to stable combustion despite similar changes in hydrogen fuel flow rate. Flame behavior before, during, and after the transient is characterized using high-speed CH* chemiluminescence imaging. Analysis of the high-speed images shows changes in flame stabilization and dynamics during the onset and decay processes. Finally, the results of this study can have implications for systems that experience variations in fuel composition, particularly in light of growing interest in hydrogen as a renewable fuel.

42 ENGINEERING↗

Fracture Analysis of Cohesive Zone Models for Modeling Residual Stress Induced Delamination in Composite Structures

A fracture study of coupon-scale composite cylinders with embedded defects was conducted with an objective to assess and validate a modeling approach using two available cohesive material models. The study included experimental and simulation evaluations of initiation of crack growth and progression. Interrupted thermal experiments used acoustic emissions monitoring to identify the onset of crack progression during each cooling interval and ultrasonic scanning provided images of defect growth. Verification, validation, and uncertainty quantification (VVUQ) processes were performed in the assessment of the simulation predicted temperature at which crack propagation begins (quantity of interest). The Sobol sensitivity analysis identified the hoop direction elastic modulus in the carbon fiber reinforced polymer (CFRP) plies as the most influential parameter for simulations using both cohesive models, accounting for at least 70% of the variation in the temperature at crack propagation. The UQ temperature range for the Tvergaard-Hutchinson model was higher (more conservative) than the experimental acoustic measurement indicators of crack progression, while the temperature range for the Thouless-Parmigiani model enveloped the experimental data points for the primary defect size of 0.75 x 1 in. The simulations could not capture the stable crack growth indicated in the experiments. This is likely due to the models’ inability to represent anisotropic fracture toughness attributed to the structure of the orthotropic fiber weave in a woven composite laminate.

42 ENGINEERING↗

Electromagnetic Analysis of ITER Electron Cyclotron Emission Components and Support Structures for the Preliminary Design Review

The ITER electron cyclotron emission (ECE) diagnostic system is located at Design Shield Module (DSM) 2, Equatorial Port (EP) 9 to measure electron temperature profile and electron temperature fluctuations and also to assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two couples of mirrors for two different optical views. This ECE diagnostic system shall be exposed to significant power due to unabsorbed electron cyclotron heating power in the plasma. It shall also receive large electromagnetic (EM) loads up to 100 MN/m 3 force density due to the eddy currents generated by the short 16-ms transient plasma disruption. The global EM models of EP12 and EP11 with the worst case of plasma disruption Major Disruption Downward Exponential 16 ms Category III (MD_DW_EXP16MS_CATIII) have been provided by ITER. However, the support structure for these critical components, including four mirrors, is still under the design and development phase, especially during the period of preliminary design reviews (PDRs). For those small in-vessel components that are not modeled explicitly for the PDRs, we can extract the magnetic field (B) data and flux variations (dB/dt) at those locations from the global model and then use the B and dB/dt method to calculate the eddy currents, EM forces, and torques on those components in the local submodels. In conclusion, this paper will present the detailed PDR EM analysis results of the ECE components and support structures.

Electromagnetic↗

Sifting for a Stream: The Morphology of the $300S$ Stellar Stream

Stellar streams are sensitive laboratories for understanding the small-scale structure in our Galaxy’s gravitational field. Here, we analyze the morphology of the $300S$ stellar stream, which has an eccentric, retrograde orbit and thus could be an especially powerful probe of both baryonic and dark substructures within the Milky Way. Due to extensive background contamination from the Sagittarius stream (Sgr), we perform an analysis combining Dark Energy Camera Legacy Survey photometry, $Gaia$ DR3 proper motions, and spectroscopy from the Southern Stellar Stream Spectroscopic Survey (S 5 ). We redetermine the stream coordinate system and distance gradient, then apply two approaches to describe $300S$ ’s morphology. In the first, we analyze stars from $Gaia$ using proper motions to remove Sgr. In the second, we generate a simultaneous model of $300S$ and Sgr based purely on photometric information. Both approaches agree within their respective domains and describe the stream over a region spanning 33° . Overall, $300S$ has three well-defined density peaks and smooth variations in stream width. Furthermore, $300S$ has a possible gap of ~4.7 and a kink. Dynamical modeling of the kink implies that $300S$ was dramatically influenced by the Large Magellanic Cloud. This is the first model of $300S$ ’s morphology across its entire known footprint, opening the door for deeper analysis to constrain the structures of the Milky Way.

Milky Way galaxy↗

MAL33 drives natural variation in maltose metabolism in Saccharomyces eubayanus

Maltose is one of the most abundant sugars in brewer’s wort, and its efficient utilization is critical for successful fermentation. However, maltose consumption varies naturally among Saccharomyces eubayanus strains isolated from different host trees, such as Quercus and Nothofagus. To identify the genetic determinants underlying these phenotypic differences, we performed bulk segregant analysis (BSA) and quantitative trait loci (QTL) mapping using an F 2 offspring derived from QC18 (Quercus-associated) and CL467.1 (Nothofagus-associated) strains. QTL mapping identified two significant genomic regions on subtelomeric loci of chromosomes V-R and XVI-L, each containing complete MAL loci composed of MAL32 (encoding maltase), MAL31 (transporter), and MAL33 (transcriptional activator) genes. Comparative polymorphism analyses identified mutations in MAL32 and MAL33 of QC18, including frameshift mutations resulting in premature stop codons. Functional validation demonstrated that the heterologous expression of MAL33 ChrV from CL467.1 fully restored maltose utilization in QC18, indicating the functional presence of MAL33 cis-regulatory sequences and MAL32 and MAL31 genes in QC18. While structural protein predictions identified truncation and impaired functionality in the maltose-responsive activation domain of Mal33p from QC18, overexpression of QC18’s own MAL33 ChrV allele also improved maltose metabolism, suggesting dosage-dependent transcriptional limitations rather than complete functional loss. These results indicate that allelic variations in the maltose-responsive activation domain of Mal33p result in differences in maltose consumption between strains. Here, we hypothesized that reduced maltose metabolism in QC18 is an adaptive response to the distinct sugar composition in Quercus robur bark, contrasting with the starch-rich environment of Nothofagus pumilio. These findings highlight subtelomeric MAL gene diversity as a reservoir of genetic variation, representing a key evolutionary mechanism that influences maltose adaptation among natural Saccharomyces isolates.

evolutionary plasticity↗

Leveraging microfluidic dielectrophoresis to distinguish compositional variations of lipopolysaccharide in E. coli

Lipopolysaccharide (LPS) is the unique feature that composes the outer leaflet of the Gram-negative bacterial cell envelope. Variations in LPS structures affect a number of physiological processes, including outer membrane permeability, antimicrobial resistance, recognition by the host immune system, biofilm formation, and interbacterial competition. Rapid characterization of LPS properties is crucial for studying the relationship between these LPS structural changes and bacterial physiology. However, current assessments of LPS structures require LPS extraction and purification followed by cumbersome proteomic analysis. This paper demonstrates one of the first high-throughput and non-invasive strategies to directly distinguish Escherichia coli with different LPS structures. Using a combination of three-dimensional insulator-based dielectrophoresis (3DiDEP) and cell tracking in a linear electrokinetics assay, we elucidate the effect of structural changes in E. coli LPS oligosaccharides on electrokinetic mobility and polarizability. We show that our platform is sufficiently sensitive to detect LPS structural variations at the molecular level. To correlate electrokinetic properties of LPS with the outer membrane permeability, we further examined effects of LPS structural variations on bacterial susceptibility to colistin, an antibiotic known to disrupt the outer membrane by targeting LPS. Our results suggest that microfluidic electrokinetic platforms employing 3DiDEP can be a useful tool for isolating and selecting bacteria based on their LPS glycoforms. Future iterations of these platforms could be leveraged for rapid profiling of pathogens based on their surface LPS structural identity.

59 BASIC BIOLOGICAL SCIENCES↗

Deterministic modelling of the SPERT IV reactor transients using the multi-physics capability of WIMS

The ANSWERS WIMS reactor physics code is being developed for whole core multi-physics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding. An integrated thermal hydraulic solver permits temperature and density variations to feedback to the neutronics calculation. This capability can be applied to both steady state and time dependent transient problems. Nuclear reactor design and safety case development requires assessment of a range of reactor transients, to inform both normal operation limitations and accident scenario analysis. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver for the simulation of reactor transients in whole core models. To support the validation of the multi-physics capability of WIMS, this capability has been applied to the reactivity insertion transient experiments performed in the plate fuelled SPERT-IV reactor. This study employs WIMS using the whole core solver MERLIN, which calculates the time-dependent flux distribution and magnitude, coupled to ARTHUR, which solves for the time-dependent thermal hydraulics solution, and provides thermal feedback via the cross sections generated by GEOM, which performs resonance shielding calculations and generates cross section data for the plate geometry. This spatial kinetics model, with dynamic cross section generation, allows for variations in the temperature and neutron flux profile with time.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Dynamical Large Deviations of Two-Dimensional Kinetically Constrained Models Using a Neural-Network State Ansatz

We use a neural-network ansatz originally designed for the variational optimization of quantum systems to study dynamical large deviations in classical ones. We use recurrent neural networks to describe the large deviations of the dynamical activity of model glasses, kinetically constrained models in two dimensions. We present the first finite size-scaling analysis of the large-deviation functions of the two-dimensional Fredrickson-Andersen model, and explore the spatial structure of the high-activity sector of the South-or-East model. These results provide a new route to the study of dynamical large-deviation functions, and highlight the broad applicability of the neural-network state ansatz across domains in physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First M87 Event Horizon Telescope Results. IX. Detection of Near-horizon Circular Polarization

Event Horizon Telescope (EHT) observations have revealed a bright ring of emission around the supermassive black hole at the center of the M87 galaxy. EHT images in linear polarization have further identified a coherent spiral pattern around the black hole, produced from ordered magnetic fields threading the emitting plasma. Here we present the first analysis of circular polarization using EHT data, acquired in 2017, which can potentially provide additional insights into the magnetic fields and plasma composition near the black hole. Interferometric closure quantities provide convincing evidence for the presence of circularly polarized emission on event-horizon scales. We produce images of the circular polarization using both traditional and newly developed methods. All methods find a moderate level of resolved circular polarization across the image ( $\langle$| v |$\rangle$ < 3.7%), consistent with the low image-integrated circular polarization fraction measured by the Atacama Large Millimeter/submillimeter Array (| v int |< 1%). Despite this broad agreement, the methods show substantial variation in the morphology of the circularly polarized emission, indicating that our conclusions are strongly dependent on the imaging assumptions because of the limited baseline coverage, uncertain telescope gain calibration, and weakly polarized signal. We include this upper limit in an updated comparison to general relativistic magnetohydrodynamic simulation models. This analysis reinforces the previously reported preference for magnetically arrested accretion flow models. We find that most simulations naturally produce a low level of circular polarization consistent with our upper limit and that Faraday conversion is likely the dominant production mechanism for circular polarization at 230 GHz in M87*.

79 ASTRONOMY AND ASTROPHYSICS↗

Targeted Dy intercalation under graphene/SiC for tuning its electronic band structure

Metal intercalation of graphene is a promising method to tune its electronic band structure and generate novel electronic and topological phases. The tuning depends critically on the ability to bond the intercalated atoms at predesigned, subsurface interlayer locations because the emerging band structure depends on metal location. In this work, we have studied Dy intercalation under single-layer graphene (SLG) on SiC using spot profile analysis–low-energy electron diffraction and scanning tunneling microscopy (STM). The experimental work is complemented with density-functional theory (DFT) analysis. Because different diffraction spots originate from different subsurface interlayer regions, it is possible to identify changes in the intercalation location by monitoring the spot intensity as a function of growth conditions. DFT calculations of the chemical potential as a function of intercalated Dy coverage support the variation of the stability of the intercalated phase at different intercalated locations. The preferred location is confirmed from STM studies showing the removal of the 6 × 6 moiré corrugation at the preferred location, observed at higher Dy coverage.

2-dimensional systems↗

Considering interplay between multiple physical phenomena to elucidate single crystal-like texture, phase transformations, and mechanical behavior of directed energy deposited SS316L

A widespread implementation of large scale additive manufacturing (AM) processes, such as wire arc-directed energy deposition (WA-DED) AM can transform the current manufacturing supply chain networks. Naturally, such implementation requires control of the microstructural attributes, such as texture and phase evolution in the processed alloys. Currently, the texture evolution in fusion-based AM (F-BAM) processes is majorly rationalized by the phenomena occurring only during solidification. However, such rationalization is insufficient for understanding the evolution of primary and secondary crystallographic orientations, and consequently, fails to offer a comprehensive understanding and control of overall texture in F-BAM processed alloys. To this end, we report a single crystal (SX)-like texture in WA-DED processed SS316L for the first time. Furthermore, we assess the physical phenomena that may lead to such unique microstructural evolution during WA-DED AM. Subsequently, using microstructural characterization spanning the build height and thermomechanical simulations we investigate the effect of competitive growth and epitaxial growth occurring during solidification and thermally induced plastic deformation occurring post solidification on the overall texture of WA-DED processed SS316L. A spatial variation in solidification pathway is also established and correlated with variation in undercoolings across the build. Tensile tests revealed a strong orientation dependence of deformation mechanisms with over 110% elongation to failure of specimens deformed along <011>. Such anisotropy is rationalized using Schmid's analysis of dislocation slip and deformation twinning. Importantly, overall, the mechanisms outlined in this work will facilitate an enhanced understanding and subsequent control of texture evolution, solidification behavior and mechanical behavior of WA-DED processed steels.

36 MATERIALS SCIENCE↗

Mitigating Injection‐Induced Seismicity Along Basement Faults by Extraction: Application to 2016–2018 Pohang Earthquakes

Abstract The risk of induced earthquake is considered as the major obstacle for sustainable subsurface energy activities. As a mitigation approach, mechanistic risk assessment can be essential to minimize the potential of seismic hazard. Coupled hydro‐mechanical modeling approach is integrated with geomechanical analysis to better assess the occurrence of seismicity induced by injection and/or extraction, which will provide a basis to update operational conditions (e.g., number and location of mitigating wells, rate and duration of operations) for mitigating earthquake nucleation. The Coulomb stability analysis and seismicity rate estimates from this generic study indicate that simultaneous extraction can relieve injection‐induced pore‐pressure buildup, but also transfer and accumulate elastic energy along the faults. Field‐scale simulation with variation in well design for the Pohang EGS site, South Korea, show that multiple wells aligned across a preexisting fault form gradients of pressure and stress fields normal to the fault, which can enhance the direct impact of diffusion process and poroelastic stressing on the fault stability perturbation. Our mechanistic study suggests multiphysics‐oriented mitigation strategies that consider (1) possible mechanical deformation, caused by hydro‐mechanical coupling effects, and (2) proper configuration of well operations based on site‐specific geological characteristics.

Chang, K. W.↗

Stabilizing lattice oxygen redox in layered sodium transition metal oxide through spin singlet state

Abstract Reversible lattice oxygen redox reactions offer the potential to enhance energy density and lower battery cathode costs. However, their widespread adoption faces obstacles like substantial voltage hysteresis and poor stability. The current research addresses these challenges by achieving a non-hysteresis, long-term stable oxygen redox reaction in the P3-type Na 2/3 Cu 1/3 Mn 2/3 O 2 . Here we show this is accomplished by forming spin singlet states during charge and discharge. Detailed analysis, including in-situ X-ray diffraction, shows highly reversible structural changes during cycling. In addition, local CuO 6 Jahn-Teller distortions persist throughout, with dynamic Cu-O bond length variations. In-situ hard X-ray absorption and ex-situ soft X-ray absorption study, along with density function theory calculations, reveal two distinct charge compensation mechanisms at approximately 3.66 V and 3.99 V plateaus. Notably, we observe a Zhang-Rice-like singlet state during 3.99 V charging, offering an alternative charge compensation mechanism to stabilize the active oxygen redox reaction.

25 ENERGY STORAGE↗