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 505 records · Page 28

Deciphering controls on metal migration within floodplains: The critical role of redox environments on metal-organic complexes (Final Technical Report)

Our project elucidates the effect of differing hydrologic regimes on nutrient and contaminant dynamics within floodplains, and we explore the impact of varying sources of hydrologic extremes on nutrient-contaminant fate and transport. We use a combination of field measurements and laboratory experiments to examine the relationships between redox conditions, functionality of dissolved organic matter, metal speciation (specifically examining metal-ligand complexes), and nutrient-contaminant fate, particularly during hydrologic extremes. Floodplain biogeochemical variation through climatic extremes and beaver-induced hydrologic variation provide a unique look at controls on nutrient-contaminant fate and transport. Through our work, we have developed a novel approach that provides an unprecedented ability to resolve aqueous metal-organic complexes that reveals metal binding preferences for natural organic matter based on metal-ligand chemistry. We have developed novel methods to separate, quantify, and characterize organic-metal complexes in natural DOM that pair with field experiments comprising floodplains having systematic variation in biogeochemical conditions. We also developed a comprehensive computational approach for molecular characterization of metal-organic species that merges LC-coupled high resolution (HR) mass spectrometry and ICP-MS analysis. Importantly, we show that beaver dams have an outsized influence on hydrologic extremes, dwarfing the impact of climatic variation, that ultimate control biogeochemical conditions and water quality. Additionally, beaver dams impose a hydrologic driver that leads to improved, river water quality by divert water through soils and the hyporheic zone. We further show that iron(III) minerals have a dominant control on carbon mineralization across floodplains independent of hydrologic conditions, and that iron and aluminum associations with organic matter are a controlling factor in carbon mineralization. Ultimately, our work is helping to advancing a robust predictive understanding of how hydrologic changes in watersheds affect water quality and carbon and metal cycling in riparian zones.

54 ENVIRONMENTAL SCIENCES↗

Unbalanced Nested Random Effects Estimation of Variance Components

To investigate the contributing factors of variance in the measurement of an iodine 127 (127I) sample, we implement a nested random effects analysis of variance (ANOVA). Historically, the reported uncertainty on a measurement of 127I has been obtained by methods of forward uncertainty propagation because there is typically only one replicate of a given sample for which to estimate the uncertainty. When samples are processed there are several types of quality control (QC) standards analyzed along-side the unknown samples with two to five replicates for each. Assuming the variance observed in these replicate QC standards is representative of that of the unknown samples, we use these data in a nested random effects ANOVA to estimate the total uncertainty of a measurement. We demonstrate this approach with two sets of measurements from Idaho National Laboratory and compare the results with the forward uncertainty propagation approach. Variance component estimates for the coarsest level of the nesting structure were most imprecise because replicates were most limited at these levels. We find that the results largely agree between forward propagation and ANOVA, and the greatest contributors of variance are due to instrument variation and chemical processing, with human processing being among the smallest contributors. This analysis provides reassurance that the reported uncertainties using forward propagation are reasonable and the process is well controlled. We propose a future designed experiment to increase replicates at the coarsest level of the hierarchy to improve estimates of these variance components.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Variation of effective filament diameter, irreversibility field, anisotropy, and pinning efficiency in Bi-2212 round wires

Abstract In recent years Bi 2 Sr 2 CaCu 2 O x (Bi-2212) received increasing attention due to its round wire multifilamentary architecture, a unique feature in high- T c superconductor. In fact, round wires are preferable for magnet designs, including solenoids for nuclear magnetic resonance (NMR) or research purpose and accelerator magnets. However, due to the narrow over-pressure heat treatment conditions necessary to obtain high J c and to the peculiar microstructure of Bi-2212 wires, a full understanding of the correlations between the different properties has not yet been developed. In this paper we investigate the effect of a vital part of Bi-2212 optimization, the maximum heat-treatment temperature T max in the range of 885 °C–896 °C, on the variations of J c , effective filament diameter d eff , anisotropy γ , INTER- and intra-grain irreversibility fields and pinning energies U 0 , all critical parameters in unravelling the complex mix of vortex pinning and connectivity that ultimately determines the critical current density. We found that d eff of the higher J c wires heat-treated at lower temperature is much smaller than for the lower J c wires. Moreover, a systematic increase of the irreversibility field and a decrease of the intrinsic Bi-2212 anisotropy underpins the higher J c . The analysis of the pinning energies reveals that there is little sample-to-sample variation in the INTER-grain pinning, whereas in all samples the intra-grain pinning has an enhancement below ∼40–45 K becoming more and more evident with increasing J c . These results suggest that the overall J c performance are not only related to the wire microstructure and connectivity, which obviously affect the INTER-grain properties, but they are also intimately related to the intrinsic and intra-grain properties such as γ and U 0 .

Physics↗

Assessing the numerical stability of physics models to equilibrium variation through database comparisons on DIII-D

High fidelity kinetic equilibria are crucial for tokamak modeling and analysis. Manual workflows for constructing kinetic equilibria are time consuming and subject to user error, motivating development of automated equilibrium reconstruction tools to provide accurate and consistent reconstructions for downstream physics analysis. These automated tools also provide access to kinetic equilibria at large database scales, which enables the quantification of general uncertainties arising from equilibrium reconstruction techniques. In this paper, we compare a large database of DIII-D kinetic equilibria generated manually by physics experts to equilibria from automated kinetic reconstruction tools, assessing the impact of reconstruction method on equilibrium parameters and resulting magnetohydrodynamic stability calculations. We find agreement among scalar parameters, whereas profile quantities, such as the bootstrap current, show larger disagreements. We analyze ideal kink and classical tearing stability with DCON and STRIDE respectively, finding that the kink stability calculation is generally more robust than the tearing index Δ' calculation. We find that in 90% of cases, both kink stability classifications are unchanged between the manual expert and automated kinetic equilibria.

CAKE↗

Impact of U-10Mo HALEU fuel element tolerances on the Massachusetts Institute of Technology reactor safety and operational performance – Thermal hydraulics

The U.S. is coordinating efforts for the conversion of six U.S. High Performance Research Reactors (USHPRR) including one critical facility from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts of Institute of Technology Reactor (MITR), and achieve similar performance, high-assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being evaluated. The impact of the fabrication specification and tolerances was assessed following the preliminary design of the MITR LEU fuel elements using the U-10Mo monolithic alloy. This research focuses on the analysis of fabrication specification impact on thermal hydraulics (TH) characteristic of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., coolant channel gap thickness, fuel plate thickness, etc.). The analyses are performed based on an all-fresh LEU fuel conversion plan identified in a preliminary safety analysis report submitted to the Nuclear Regulatory Commission. The reactor power margin to the onset of nucleate boiling (ONB) is assessed under the limiting safety system settings (LSSS), where a scram occurs, to ensure there is sufficient margin to the reactor safety limit, which is defined by the onset of flow instability that occurs after the ONB. The best estimate plus uncertainty approach is employed to analyze this TH characteristic, which yields realistic results while maintaining adequate conservatism, utilizing a statistical uncertainty propagation method with the STAT7 code. The TH characteristic is analyzed as a function of the variability of the specification parameters resulting from the fabrication process. The main findings of this study show that the MITR core can meet the TH safety and operational requirements at the all-LEU initial core startup (cycle 1), selected transition cycles (most reactive cycle and most limiting cycle: cycle 3 and 5, respectively) and equilibrium (cycle 14) cores under all limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core power margin to ONB on those specification parameters that have the most direct impact on TH performance is non-linear but monotonically decreasing within the specification tolerances. The third order polynomial fit curves are reported in detail for selected limiting cases and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion↗

Circularity Futures Workshop Series: Summary Report

The aim of this report is to synthesize key feedback received from the three-part Circularity Futures workshop series held in Spring 2024. The workshop series was conducted by the National Renewable Energy Laboratory (NREL) on behalf of U.S. Department of Energy, Office Energy Efficiency and Renewable Energy (EERE), and was broken into three workshops: Workshop 1 - Circularity Analysis Needs and Priorities; Workshop 2 - Circularity Metrics and Indicators; and Workshop 3 - Circularity Data. Together, the workshops focused on identifying the existing priorities and gaps in the circularity modeling space, understanding different stakeholders' use and interpretation of circularity metrics and indicators, identifying common data gaps and data quality challenges, and assessing the robustness of available solutions. The workshop series brought a diverse group of stakeholders - including representatives from U.S. government offices, national labs, nonprofit organizations, industry, and academia - to collect first-hand feedback on needs, priorities, challenges and opportunities in the circularity modeling and analysis space. The workshop discussions highlighted numerous common needs, priorities and challenges among the interviewed groups. Several topics were frequently discussed, including: 1) Circularity as a pathway for sustainable economic growth: While circularity is generally defined in terms of resource conservation and reducing wasteful disposal of materials, participants agreed that circular strategies should serve broader economic, environmental, and social goals. It is therefore crucial for circularity analysis to look beyond waste reduction and instead evaluate a variety of impact metrics such as cost savings, job creation, air quality, and pollutant emissions. Mutli-criteria decision-making frameworks may be useful for making sense of disparate metrics and evaluating tradeoffs between impact categories.; 2) Economic and social factors are not well understood: Underdevelopment of existing end-of-life (EOL) management infrastructure, inconsistent standardization codes and policy space in reusing recycled content, and suboptimal collection and sorting strategies collectively contribute to uncertainty about the economic potential of circular pathways. The latter observation is consistent among all technologies but more emphasized for renewable energy systems. Social impacts of circularity practices are less understood and less researched than other sustainability aspects.; 3) Inconsistent methods for assessing emerging technologies: LCA and TEA results vary widely depending on the assumptions made with regards to market adoption of new technologies. Emerging technologies suffer limited availability of data needed to conduct a robust circularity analysis. Yet, understanding projected impacts of proposed nascent technology is a key need for different stakeholder groups.; and 4) Lack of temporally and geospatially explicit data: There is a need for open data that represents variations in circularity technologies over time and location. The lack thereof leads to aggregated and potentially misrepresented results in circularity analysis. Sensitivity analyses should be included to verify whether options perceived as more sustainable align with real-world practices.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Real-Space Pseudopotential Method for the Calculation of Third-Row Elements X-ray Photoelectron Spectroscopic Signatures

X-ray photoelectron spectroscopy (XPS) is a powerful characterization technique that unveils subtle chemical environment differences via core– electron binding energy (CEBE) analysis. We extend the development of realspace pseudopotential methods to calculating 1s, 2s, and 2p 3/2 CEBEs of third-row elements (S, P, and Si) within the framework of Kohn–Sham density-functional theory (KS-DFT). The new approach systematically prevents variational collapse and simplifies core-excited orbital selection within dense energy level distributions. However, careful error cancellation analysis is required to achieve accuracy comparable to all-electron methods and experiments. Combined with real-space KS-DFT implementation, this development enables large-scale simulations with both Dirichlet boundary conditions and periodic boundary conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advancing SiC clad fuel performance model: bridging micro- and macro-scale models and experimental validation

This report presents a workflow for advancing fuel performance modeling of SiC composite cladding for light-water reactors by linking microscale, experimental data-informed finite element analysis with rod-scale fuel performance codes such as BISON. The workflow uses X-ray computed tomography (XCT) to capture the actual geometry and processing-induced defects of as-fabricated SiC composite tube specimens, particularly porosity and wall-thickness variations, and converts the segmented XCT volumes into image-based finite element meshes for high fidelity structural analysis.

Koyanagi, Takaaki [Oak Ridge National Laboratory (↗

Stability analysis of the Eulerian–Lagrangian finite volume methods for nonlinear hyperbolic equations in one space dimension

In this paper, we construct a novel Eulerian–Lagrangian finite volume (ELFV) method for nonlinear scalar hyperbolic equations in one space dimension. It is well known that the exact solutions to such problems may contain shocks though the initial conditions are smooth, and direct numerical methods may suffer from restricted time step sizes. To relieve the restriction, we propose an ELFV method, where the space-time domain was separated by the partition lines originated from the cell interfaces whose slopes are obtained following the Rakine–Hugoniot junmp condition. Unfortunately, to avoid the intersection of the partition lines, the time step sizes are still limited. To fix this gap, we detect effective troubled cells (ETCs) and carefully design the influence region of each ETC, within which the partitioned space-time regions are merged together to form a new one. Then with the new partition of the space-time domain, we theoretically prove that the proposed first-order scheme with Euler forward time discretization is total-variation-diminishing and maximum-principle-preserving with at least twice larger time step constraints than the classical first order Eulerian method for Burgers’ equation. Numerical experiments verify the optimality of the designed time step sizes.

97 MATHEMATICS AND COMPUTING↗

A comparison of the power potential for surface- and seabed-deployed tidal turbines in the San Juan Archipelago, Salish Sea, WA

The San Juan archipelago lies along the axis of tidal movement between Straits of Juan de Fuca and Strait of Georgia in the Salish Sea. The amplitude of the tidal exchange produces significant tidal currents between the islands, as well as in Rosario and Haro Straits. These are of interest as a future source of electrical power generation, given the islands’ dependence on electricity supply by a subsea cable from the mainland. Here, we evaluate the tidal current energy potential in this region through a re-analysis of measurements collected by the National Ocean Service (NOS) and a high-resolution numerical model. Given the considerable variations in water depth and vertical velocity profiles across candidate tidal energy sites, we consider the trade-offs between tidal turbines deployed from a floating platform and those anchored to the seabed. Measurement re-analysis indicates several locations that could support tidal current power generation by MW-scale turbines with an acceptable balance between turbine size, rated power, and capacity factor. Even for relatively large (30 m) turbine diameters, surface-deployed turbines would be expected to produce up to 30 % more electricity than the same turbine deployed near the seabed due to vertical shear, with this difference increasing for smaller diameter turbines. A benchmark comparison at the measurement locations with the highest power generation potential shows relatively good model fidelity, though, even with relatively high resolution, time-average power density disagreements of ± 50 % persist throughout the water column. Nonetheless, the simulation identifies several locations with power generation potential more than twice as high as locations in the measurement re-analysis. These sites were either not surveyed by NOS or excluded due to data quality issues. Altogether, these results demonstrate the complimentary of models and measurements for tidal energy site assessment and suggest that tidal currents could be an important source of electricity generation in the San Juan archipelago.

16 TIDAL AND WAVE POWER↗

Antiprotons and Elementary Particles over a Solar Cycle: Results from the Alpha Magnetic Spectrometer

We present results over an 11-year Solar cycle of cosmic antiprotons based on 1.1 × 10 6 events in the rigidity range from 1.00 to 41.9 GV. The $\bar{𝑝}$ fluxes exhibit distinct properties. The magnitude of the $\bar{𝑝}$ flux temporal variation is significantly smaller than those of 𝑝, 𝑒 − , and 𝑒 + . A hysteresis between the $\bar{𝑝}$ fluxes and the 𝑝 fluxes is observed, whereas the $\bar{𝑝}$ and 𝑒 − fluxes show a linear correlation. With a model-independent analysis, we found a universal relation between the shape of the rigidity spectrum and the magnitude of flux temporal variation over an 11-year Solar cycle for both positively and negatively charged particles. The simultaneous results on $\bar{𝑝}$ and 𝑝, 𝑒 − , and 𝑒 + provide unique information for understanding particle transport in the Solar System as a function of mass, charge, and spectral shape.

cosmic ray composition & spectra↗

ENSO-Induced Teleconnection: Process-Oriented Diagnostics to Assess Rossby Wave Sources and Ambient Flow Properties in Climate Models

Climate model fidelity in representing ENSO-induced teleconnection is assessed with process-oriented diagnostics that examine a chain of processes, from equatorial Pacific precipitation to the midlatitude circulation pattern over the Pacific–North American regions. Such processes are rarely addressed during model development. Using an upper-tropospheric divergent level, local vorticity gradient of the ambient zonal flow (∂ 2 U¯/∂y 2 ) and a restoring force for Rossby waves (β * ) are estimated, the equivalent barotropic vorticity equation is solved, and an anomalous Rossby wave source (RWS') quantified. The analysis is applied to AMIP5 and AMIP6 simulations. For a realistic circulation response representation, the hypothesis that models accurately represent the strength and location of RWS', and spatial variations in β * is tested. Compared to AMIP5, in AMIP6 there are clear improvements in representing RWS' and β * . To validate the hypothesis, the analysis identifies two metrics: spatially coherent RWS' in the subtropical North Pacific, and longitudes of negative β * over the western-central North Pacific. By projecting these metrics in two and three-dimensional views, improvements or degradations in model versions are apparent. If a model’s fidelity in representing ∂ 2 U¯/∂y 2 and RWS' are compromised, then radiated Rossby waves are reflected more equatorward, resulting in zonally elongated circulation anomalies over the central North Pacific. Thus, during climate model development, applying this analysis frequently will keep a regular check on the fidelity of the modeled response to anomalous El Niño convection in conjunction with changing model ambient flow dependencies. As a result, this analysis is intended to form a process-oriented diagnostics package, a community contribution to the NOAA Model Diagnostics Task Force.

54 ENVIRONMENTAL SCIENCES↗

Signature Analysis Utilizing a Dynamic Molten Salt Reactor Model for MC&A

Moving from traditional fixed fuel nuclear reactor systems to a mobile, dynamic fuel system that has dissolved special nuclear material in a molten salt is a paradigm shift in several respects. One major consideration is how to develop effective nuclear material controls and accounting for these novel reactor systems. The Molten Salt Reactor Experiment is one example of a critical molten salt system and provided a significant reference library based on the documented effort. But it was low thermal power that was not intended to reflect a commercial scale electricity production design. Therefore, to facilitate and assist vendors with domestic licensing considerations, research is underway to identify methods that could be used for domestic safeguards approaches for these novel reactor systems. This research presents the results of a signature analysis of data generated from three simulated scenarios using the molten salt demonstration reactor model defined in the Transient Simulation Framework of Reconfigurable Modules. Each scenario provides 1 hour isotope inventories over a 180 day period. The scenarios investigated provide test cases to examine if direct gamma-ray spectroscopy of the fuel can be used to identify changes comparing a base case (no reactivity control and fixed fission contribution) to an insertion of reactivity (10 pcm no change in fission composition) and a change in fission composition. The analysis demonstrates that monitoring the total count rate in a highly collimated high-resolution photon energy spectrum is sensitive to perturbations imposed into the reactor model. The total photon count rate changes ≈2% for the fission composition change and ≈4.5% for the reactivity insertion compared to the base case. However, both scenarios show an increase in the total photon count rate. The total photon count rate can be used to identify changes due to power (number of fissions) but not due to a change in the material undergoing fission. To distinguish between the two cases of increased power, the photon spectrum would require an intensive analysis technique. A photon peak count rate ratio analysis could be used to identify changes in the fissile material fission generation in the core through identification of a static peak that shows little variation to the source of fission and a highly varying peak. The photon peak strength will ultimately be determined by the isotope’s fission yield. A preliminary analysis investigating the coupling of the isotope’s fission yield and its concentration in the fuel salt derived from the modeling has been performed. A ratio analysis of the photon peak count rates of 140 La to 99 Mo demonstrated that the reactivity insertion creates a distinct difference in the ratio compared to the fission composition change scenario.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The scalar chemical potential in cosmological collider physics

Non-analyticity in co-moving momenta within the non-Gaussian bispectrum is a distinctive sign of on-shell particle production during inflation, presenting a unique opportunity for the “direct detection” of particles with masses as large as the inflationary Hubble scale ( H ). However, the strength of such non-analyticity ordinarily drops exponentially by a Boltzmann-like factor as masses exceed H . In this paper, we study an exception provided by a dimension-5 derivative coupling of the inflaton to heavy-particle currents, applying it specifically to the case of two real scalars. The operator has a “chemical potential” form, which harnesses the large kinetic energy scale of the inflaton, \( {\overset{\cdot }{\phi}}_0^{1/2}\approx 60H \) , to act as an efficient source of scalar particle production. Derivative couplings of inflaton ensure radiative stability of the slow-roll potential, which in turn maintains (approximate) scale-invariance of the inflationary correlations. We show that a signal not suffering Boltzmann suppression can be obtained in the bispectrum with strength f NL ~ \( \mathcal{O} \) (0 . 01–10) for an extended range of scalar masses \( \lesssim {\overset{\cdot }{\phi}}_0^{1/2} \) , potentially as high as 10 15 GeV, within the sensitivity of upcoming LSS and more futuristic 21-cm experiments. The mechanism does not invoke any particular fine-tuning of parameters or breakdown of perturbation-theoretic control. The leading contribution appears at tree-level , which makes the calculation analytically tractable and removes the loop-suppression as compared to earlier chemical potential studies of non-zero spins. The steady particle production allows us to infer the effective mass of the heavy particles and the chemical potential from the variation in bispectrum oscillations as a function of co-moving momenta. Our analysis sets the stage for generalization to heavy bosons with non-zero spin.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Myna: Connecting powder bed fusion build data to simulation tools for digital twin applications

Additive manufacturing (AM), as a digital process, can generate a detailed digital thread linking a part’s design and manufacturing to its operational performance. As AM systems advance, an increasing amount of process data is stored in manufacturing databases. In principle, this data can be utilized by simulation-based digital twin approaches, such as real-time process control and asynchronous post-processing guidance. However, few tools currently exist for systematically integrating digital thread data with computational tools. Here, in this study, we propose a software package, called Myna, for connecting data from powder bed fusion processes to simulation tools. The utility of such a platform is demonstrated using build data from the Oak Ridge National Laboratory Manufacturing Demonstration Facility “Peregrine v2023-10” public dataset to automatically configure and run 54 semi-analytical 3DThesis melt pool simulations, 78 numerical Additive FOAM melt pool simulations, and 3 ExaCA microstructure simulations. The simulated, spatially registered microstructures are then compared directly with electron backscatter diffraction characterization of the corresponding as-built part locations. The resulting simulated microstructure showed variation as a function of process parameters, particularly stripe width; however, the experimental data had little variation between the microstructure texture and grain size resulting from different processing conditions. Analysis of the discrepancies suggest that it is possible a two-phase ferritic-austenitic solidification model is needed to accurately predict grain size and texture for certain stainless steel 316L feedstock compositions under powder bed fusion conditions, providing direction for future research. As illustrated here, due to the number and complexity of the simulations involved in AM process-structure–property predictions, automated methods to connect process data and simulations will remain necessary tools for testing hypotheses and implementing digital twin applications.

Knapp, Gerald L. [Oak Ridge National Laboratory (O↗

Dynamic fracture of glass fiber-reinforced ductile polymer matrix composites and loading rate effect

Here, the dynamic fracture of S-2 glass fiber-reinforced polymer matrix composites (FRPMCs) was investigated in this study. The matrix ductility was improved by a recently developed network topology modification technique via mixing partially reacted substructures (mPRS). The composite material was manufactured and characterized by micro-CT scanning and scanning electron microscopy (SEM). Dynamic single-edge notched bending (d-SENB) experiments were performed on the composites by using a modified split-Hopkinson pressure bar. Each specimen’s fracture process was visualized by ultrafast X-ray imaging. Such in-situ radiography enabled identifying the damage initiation below 50-micron scale and inspecting its propagation through the internal structures of opaque composites, thereby accurately quantifying the composites’ mechanical properties. Furthermore, the identical d-SENB experiments were designed and the digital image correlation (DIC) was employed to monitor the stress wave propagation on the composite specimens. The force and deflection measurements were modified and correlated to the physical damage processes. Besides, quasi-static SENB experiments were conducted to identify the loading rate effects on the composites’ fracture behaviors. The force and deflection history, bending stiffness, energy dissipation, and fracture toughness at different loading rates were quantified and compared. Finally, post-fracture analysis by micro-CT scanning and SEM provided physical observations on the variation of the fracture morphology by different loading rates.

42 ENGINEERING↗

Optimization of well design and CO 2 injection strategy for risk reduction in Class VI geological carbon sequestration wells

The safety and durability of Class VI wells are critical for geological carbon sequestration (GCS). However, current GCS operations face unique challenges: unlike traditional Class II wells, Class VI CO 2 injection wells operate at rates up to 100 times higher, dramatically increasing the risk of wellbore leakage and structural compromise due to severe temperature drops and associated mechanical stresses. Despite existing guidelines on material selection, there remains a substantial gap in understanding how rapid CO 2 injection rates, low surface temperatures, and variable reservoir conditions interact to threaten long-term well integrity. This study presents a comprehensive, original workflow integrating advanced analytical and numerical models for both well flow and well integrity analysis. By systematically simulating a wide range of field-relevant scenarios—including variations in injection rate, CO 2 temperature, and reservoir pressure—this work provides the first cross-validated assessment of cooling effects on wellbore. The results reveal that extreme temperature drops, up to 60 °C, can occur under high injection rates, particularly in depleted reservoirs, significantly increasing the risk of cement failure. Building on these insights, the study proposes innovative, practical well design and operational strategies, including ductile cement formulations, pre-stressing techniques, advanced insulation coatings, and proactive management of injection rates. The safety of Class VI well extends beyond simply using CO 2 resistant materials. Cement materials should possess optimal thermo-hydraulic-mechanical-chemical properties for effective performance. This work provides a scientific basis for optimizing Class VI well designs, with direct benefits for minimizing environmental risk, lowering operational costs, and enhancing the long-term reliability of GCS.

25 ENERGY STORAGE↗

Neutron Irradiation Induced Changes in Isotopic Abundance of 6Li and 3D Nanoscale Distribution of Tritium in LiAlO2 Pellets Analyzed by Atom Probe Tomography

Tritium, a radioactive isotope of hydrogen is produced by neutron irradiation of 6Li enriched LiAlO2 pellets in Tritium producing burnable absorber rods (TPBARs). Three-dimensional nanoscale mapping of Tritium in irradiated pellets is critical to understand its spatial variation and its association with various microstructural features in the pellets. Spatially resolved analysis of 6Li isotopic enrichment in such ceramic pellets before and after irradiation can provide insights to heterogeneities in Li isotopic distribution. However, Li being a light isotope, its spatially resolved compositional or isotopic analysis is a challenging task to most analytical methods. Here we used atom probe tomography, to evaluate the Li isotopic enrichment in both the LiAlO¬2 matrix and the Li deficient secondary phase LiAl5O8 before irradiation. These results were then compared with 6Li isotopic ratio measured by APT analysis of neutron irradiated pellets to clearly demonstrate the unique capability of APT to quantitatively analyze the isotopic enrichment of 6Li as well as for other light elements present in these materials. Evidence for heterogenous nanoscale distribution of 3H and O3H within irradiated pellet microstructure is also provided, which is attributed to 3H trapping in irradiation induced vacancy clusters or voids within the LiAlO2 matrix. The characterization results given here, now prove the feasibility for using APT for analyzing nanoscale spatially resolved isotopic enrichment of Li as well as other light elements such as Tritium in TPBARs. This work also paves way to adopt APT for analyzing light elements such as Li and tritium in materials used for nuclear fusion reactor applications, geochemistry, astrophysics as well as for energy storage applications.

Devaraj, Arun↗