Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “load shapes”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Investigation of nonlinear mechanical behavior of two superfine-grained graphites with DIC assisted disc splitting test

Using standardized uniaxial tensile test specimens is not practical due to the limited volumes in irradiation capsules, molten salt degradation facilities, or oxidation apparatus. The ASTM International Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite (ASTM D8289) was developed to provide a convenient way for estimating tensile strength. Unlike the traditional uniaxial tensile test in ASTM C749 (International Standard Test Method for Tensile Stress–Strain of Carbon and Graphite), which uses dog-bone shaped specimens larger than 12.95 mm × 120.65 mm, the ASTM D8289 standard uses smaller discs with diameters of 6–12.7 mm. A digital image correlation (DIC) system, which uses a full-field noncontact surface displacement measurement technique, was applied along with the ASTM D8289 disc splitting test on IG-110 samples and compared with previous measurements from Mersen 2114 samples. Results confirmed that the DIC technique can measure the surface displacement/strain on these small (Ø6 mm × 3 mm) graphite specimens with good repeatability. However, Mersen 2114 and IG-110 samples exhibited strain discrepancies when DIC measurements were compared with analytical and finite element simulation values. The loading history and strain results also indicated different mechanical behaviors between Mersen 2114 and IG-110, particularly the nonlinear behavior of the IG-110 samples. Good agreement was observed by comparing the splitting tensile strengths of two superfine-grained grades with results from other work. The specimen size effect is discussed when comparing the splitting tensile strength with corresponding uniaxial tensile strength of these two graphite grades.

Lin, Lianshan [ORNL] (ORCID:0000000203399219)↗
DOE OSTI2026

Advancing the quantification of aerosol-cloud interactions with the CALIPSO-CloudSat-Aqua/MODIS record

Aerosol-cloud-precipitation interactions are assessed over the non-polar ocean using more than 11 years of combined Aqua-MODIS, CALIPSO-CALIOP, and CloudSat products. The analysis first shows the benefit of incorporating vertically resolved aerosol extinction coefficient (σext) in aerosol-cloud interactions (ACI) assessments, demonstrating that: σext vertically collocated with the cloud layer () correlates best with cloud droplet number concentration (Nd), column-integrated aerosol optical depth (AOD) cannot explain the Nd variability in the extratropics, and the S-shape of the AOD-Nd relationship reported in previous studies is not replicated when using instead of AOD, with a Nd- linearity more consistent with in-situ studies over the ocean. ACI metric, estimated as the log-scale regression between CALIOP and MODIS Nd reveals that the eastern Pacific is the region with the strongest ACI, followed by the Southern Ocean. The susceptibility of clouds to changes in their liquid water path (LWP) and frequency of precipitation followed a 2-step calculation by combining the Nd- regression (ACI) with the regression between these macrophysical variables and Nd. LWP susceptibility is negative (LWP decreases with aerosol loading) and statistically significant over the eastern Pacific, eastern Atlantic, and extratropics. In contrast, vast areas of the tropical and subtropical ocean feature negligible changes in LWP with aerosol. Precipitation frequency susceptibility is negative, but the values are only significant over the coastal eastern Pacific and Atlantic. The findings suggest that previous modeling assessments relying on AOD may need to be revisited by taking advantage of the synergy between passive and active sensors.

Li, Zhujun↗
DOE OSTI2022

Disruption prediction with artificial intelligence techniques in tokamak plasmas

In nuclear fusion reactors, plasmas are heated to very high temperatures of more than 100 million kelvin and, in so-called tokamaks, they are confined by magnetic fields in the shape of a torus. Light nuclei, such as deuterium and tritium, undergo a fusion reaction that releases energy, making fusion a promising option for a sustainable and clean energy source. Tokamak plasmas, however, are prone to disruptions as a result of a sudden collapse of the system terminating the fusion reactions. As disruptions lead to an abrupt loss of confinement, they can cause irreversible damage to present-day fusion devices and are expected to have a more devastating effect in future devices. Disruptions expected in the next-generation tokamak, ITER, for example, could cause electromagnetic forces larger than the weight of an Airbus A380. Furthermore, the thermal loads in such an event could exceed the melting threshold of the most resistant state-of-the-art materials by more than an order of magnitude. In this work, to prevent disruptions or at least mitigate their detrimental effects, empirical models obtained with artificial intelligence methods, of which an overview is given here, are commonly employed to predict their occurrence—and ideally give enough time to introduce counteracting measures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

UNCERTAINTIES FOR IN-REACTOR TESTING OF TRITIUM LOADED GETTERS WITH TEMPERATURE CONTROL

The fourth experiment in TPBAR Materials Irradiation Separate-Effects Test (TMIST) series (TMIST-4) is planned to experimentally determine the effect of neutron irradiation on the equilibrium vapor pressure of tritium over a tritium loaded Nickel-Plated Zirconium (NPZ) getter in a simulated TPBAR service environment. This pressure will be used to estimate the tritium permeation rate. In the thermal design in the TMIST-4, two evaluations are required. 1) The cladding temperature is calculated to estimate the feasibility of detecting potential irradiation enhancement of the partial pressure for an NPZ getter. This work is required to obtain a precise cladding temperature for measurement of a low tritium partial pressure. This work has been completed and indicated that the required/target cladding temperatures based on the TMIST-4 test plan can be satisfied under the test condition range by modifying the capsule component shapes and inserting a thermally conductive metal sleeve. 2) The evaluation of uncertainties associated with the test design and operation is performed. This includes fabrication tolerances, and power and nuclear operating conditions. These factors will impact the component temperatures inside capsules and resulting the tritium partial pressure. The tritium permeation rate is a function of temperature and tritium partial pressure. The prediction of tritium permeation rates affects the feasibility of obtaining measurement resolutions required to reach reasonable conclusion for potential irradiation enhanced performance. This work will investigate the uncertainties.

TTP, TPBAR, TMIST-4, TPBAR Materials Irradiation S↗
DOE OSTI2021

Anisotropy evolution of elastomeric foams during uniaxial compression measured via in-situ X-ray computed tomography

We have characterized the three-dimensional evolution of microstructural anisotropy of a family of elastomeric foams during uniaxial compression via in-situ X-ray computed tomography. Flexible polyurethane foam specimens with densities of 136, 160 and 240 kg/m 3 were compressed in uniaxial stress tests both parallel and perpendicular to the foam rise direction, to engineering strains exceeding 70%. The uncompressed microstructures show slightly elongated ellipsoidal pores, with elongation aligned parallel to the foam rise direction. The evolution of this microstructural anisotropy during deformation is quantified based on the autocorrelation of the image intensity, and verified via the mean intercept length as well as the shape of individual pores. Trends are consistent across all three methods. In the rise direction, the material remains transversely anisotropic throughout compression. Anisotropy initially decreases with compression, reaches a minimum, then increases up to large strains, followed by a small decrease in anisotropy at the largest strains as pores collapse. Compression perpendicular to the foam rise direction induces secondary anisotropy with respect to the compression axis, in addition to primary anisotropy associated with the foam rise direction. In contrast to compression in the rise direction, primary anisotropy initially increases with compression, and shows a slight decrease at large strains. These surprising non-monotonic trends and qualitative differences in rise and transverse loading are explained based on the compression of initially ellipsoidal pores. Microstructural anisotropy trends reflect macroscopic stress-strain and lateral strain response. These findings provide novel quantitative connections between three-dimensional microstructure and anisotropy in moderate density polymer foams up to large deformation, with important implications for understanding complex three-dimensional states of deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Attributes of Bi-Directional Turbomachinery for Pumped Thermal Energy Storage

Abstract In this paper, we (i) present a methodology for determining the aerodynamic performance of bi-directional turbomachines for pumped thermal energy storage, i.e., turbomachines designed to operate as a compressor in one direction, and then as a turbine in the opposite direction, (ii) carry out performance computations for such turbomachines, and (iii) propose principles for conceptual design of these devices. Focus is placed on using the energy storage cycle not only to identify the novel requirements placed on bi-directional turbomachines but also to estimate the effect of these requirements on the efficiency of the energy storage process. In particular, the difference between aerodynamic loading in forward and in backward operation causes the blading to work at incidences leading to the performance below maximum efficiency, resulting in a lower round-trip efficiency. The description of the design principles includes determination of the number of stages, definition of nondimensional parameters for blading selection, and optimization of two-dimensional blading for bi-directional operation. The assessment of stage count shows the relationship between relative Mach number, pressure ratio, and round-trip efficiency. The nondimensional parameters are assessed through a bi-directional analogue to existing “Smith charts,” for the efficiency of single-direction turbomachines, as a function of camber and stagger. The blade shape evaluation and optimization show how the blade profile can be modified to address the requirements of a bi-directional turbomachine, enabling an increase in round-trip efficiency of two percentage points compared to a baseline double circular arc configuration.

Engineering↗
DOE OSTI2025

Adaptive PID Gain Scheduling Control for Hydropower Turbine Using Neural CDE and Stochastic Distribution Shaping

This paper introduces a gain-scheduling PID controller design strategy for hydroturbine frequency control mode. This scheme first uses real data to learn the nonlinear dynamics of the hydroturbine using neural controlled differential equations and then perturbs the obtained nonlinear system at different equilibrium points, based on which a static output feedback adaptive dynamic programming algorithm is then used to optimize the PID gains for each equilibrium point. Moreover, a continuous-time version of stochastic distribution control is proposed to further fine-tune the optimized PID gains. Finally, the controller is obtained by implementing linear interpolation between the optimized PID control gains. The simulation results show that the proposed gain-scheduling PID controller can control a larger range of operation points compared with the given fixed PID controller and the baseline method. Compared with the given fixed PID controller, the proposed gain-scheduling PID controller can regulate hydroturbine frequency against disturbances induced by power-load variation with over 50% less overshoot for some operation points.

13 HYDRO ENERGY↗
DOE OSTI2022

Computational Modeling of Multi-Pass Rolling Parameters Effect on Resulting Fuel Foil Shape

A focus of the U.S. Department of Energy is to improve production yield and reduce the cost of Low Enriched Uranium (LEU)-molybdenum alloy (U-10Mo) monolithic fuel plates that will be replacing High Enriched Uranium (HEU) oxide dispersion fuels used currently in the United States High Performance Research Reactors (USHPRR). One area of improvement is lowering the high transverse waviness and longitudinal waviness currently present within rolled foils prior to cladding to produce fuel plates. Traditional rolling manufacturing techniques for other metal foils use winders to pull and straighten the foil as it is rolled back and forth to the final thickness. This approach cannot be used to roll thin U-10Mo foils (0.008-0.025” thick) because only small castings can be rolled due to nuclear criticality safety concerns. As a result, the fuel foils are too short (1 m in length) to use traditional winders. Therefore, it is crucial to identify other rolling parameters (i.e., roller friction, axial tension load, roller diameter, and roll pass reduction percent) that might reduce transverse waviness and longitudinal waviness in the rolled fuel foil and develop a high-yield, low-cost multi-pass rolling manufacturing process. This report documents a systematic finite element modeling study to investigate the effects of numerous rolling parameters to reduce resulting transverse waviness and longitudinal waviness in the fuel foil during multi-pass rolling of U-10Mo foils. The rolling of a U-10Mo plate with initial dimensions of 1”x1”x 0.048” is modeled using Abaqus CAE. This rolling is modeled to undergo eight 20% reduction roll passes to a final fuel foil thickness of 0.01”. The elastic-plastic constitutive model of the U-10Mo alloy was input to the fuel foil rolling model. The rollers were modeled as rigid bodies. A comparison of rolling friction coefficients of 0.3 and 0.7 over a wide range of applied axial tension loads were investigated in order to evaluate the effect of using a lubricant during rolling. The effect of roller diameter on the resulting transverse waviness and longitudinal waviness of the fuel foil over a wide range of axial tension loads were also investigated by modeling rollers 7/8” and 3.75” in diameter. The results of this systematic finite element method study will aid manufacturers in producing low transverse waviness and reduced longitudinal waviness in U-10Mo fuel foils.

U-10Mo, FEA, Rolling, Residual Stress, Fuel Foil↗
DOE OSTI2025

Measurements of attenuation and scattering properties of water-based liquid scintillator

Water-based liquid scintillator (WbLS) is a hybrid detector medium which has been proposed as a fill material for future large-volume rare event searches, including ν − detectors. This family of scintillating suspensions promises waterlike attenuation and scattering behavior, while offering scintillation light yield for sub-Cherenkov-threshold events. These features may allow for improved vertex and energy resolution, and concomitant improvements in background rejection and particle identification. While subscale measurements of light yield, timing and pulse-shape have been performed in several WbLS formulations, detailed measurements of the attenuation and scattering properties of WbLS remain a critical, unresolved step along the pathway to deployment. In pursuit of a better understanding of these parameters, a “long-arm” attenuation and scattering instrument has been developed at Lawrence Livermore National Laboratory, dubbed LASE (Livermore Attenuation and Scattering Experiment). Optical property measurements have been performed using LASE in DI water, Gd-water, WbLS and Gd-WbLS. Measurements of the optical properties of WbLS (1% LAB-PPO scintillator) provided by Brookhaven National Laboratory demonstrate an attenuation minimum of 4.87 × 10 − 4 ± 1.42 × 10 − 4 cm − 1 at 450 nm, while Gd-loaded WbLS (1% LAB-PPO scintillator, 0.1% Gd) demonstrated a minimum attenuation coefficient of 4.77 × 10 − 4 ± 1.27 × 10 − 4 cm − 1 at the same wavelength. Measurements of scattering in WbLS show a scattering coefficient of 3.39 × 10 − 4 ± 1.98 × 10 − 5 cm − 1 at the 450 nm attenuation minimum, while Gd-WbLS has a scattering coefficient of 2.8 × 10 − 4 ± 1.63 × 10 − 5 cm − 1 at that wavelength. These scattering and attenuation coefficients are significantly larger than measured for DI water at similar wavelengths ( 8.8 × 10 − 5 cm − 1 ± 3.8 × 10 − 5 attenuation at 430 nm, 2.2 × 10 − 5 ± 1.28 × 10 − 6 cm − 1 scattering). Though the WbLS attenuation and scattering coefficients measured are lower than the corresponding values published for ultrapure oils and LAB-PPO liquid scintillator, the material tested attenuates significantly more than water across the blue and green portions of the visible spectrum. This attenuation is markedly stronger at wavelengths under 430 nm. Published by the American Physical Society 2025

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

EV charging site day-ahead load prediction in a synthetic environment for RL based grid-informed charging

Ensuring grid health in the face of increasing demand for power is an emerging challenge especially due to transportation electrification. A free market approach to influencing electric vehicle (EV) load through grid-informed hourly dynamic pricing is introduced in this work. The setting of charging price is done by a reinforcement learning (RL) agent that learns the complicated dynamics by interacting with a synthetic environment. This synthetic environment is a combination of distribution feeder simulation, EV charger user behavior dynamics, and EV charging simulation. A key module in this synthetic environment involves obtaining the day-ahead charging profile of EV charging stations based on real-world past data. The day-ahead prediction is also useful in other traditional optimizations related to EV charge scheduling. The proposed approach involves using EV charging data from two different past time horizons – one to determine the shape of the daily profile and the other to determine a scaling value to capture actual energy consumption. Real-world charging data over many years from the ACN charging network has been used to demonstrate the ability to predict the day-ahead profile with only charge session data. Both Python and MATLAB have been used for data cleaning, processing, analysis, and prediction.

Suryanarayana, Harish↗
DOE OSTI2020

Experimental verification of a crystal plasticity-based simulation framework for predicting microstructure and geometric shape changes: Application to bending and Taylor impact testing of Zr

This paper is concerned with experimental verification of a recently developed multi-scale simulation framework for plastic deformation of metallic materials from quasi-static to impact deformation conditions. The framework is a visco-plastic self-consistent (VPSC) polycrystalline model embedded in an implicit finite element method (FE-VPSC) to provide a microstructure-sensitive constitutive response at each material point. Each material point of the FEM model is a polycrystalline aggregate with crystallographic deformation mechanisms operating at the single crystal scale with their evolving activity based on a dislocation density-based hardening law and texture. Four beams and three cylinders machined in different orientations from a textured plate of high-purity zirconium are tested quasi-statically in 4-point bending and at speeds of 100 m/s, 170 m/s and 243 m/s during Taylor impact tests, respectively. The variation in dimensional changes resulting from different sample orientations in the plate with respect to loading directions is measured for each sample. Moreover, texture and twinning characterization is performed using electron backscattered diffraction (EBSD). The deformation processes and underlying evolution of microstructure are successfully simulated using the FE-VPSC framework. In doing so, the model parameters are optimized and validated across a broad range of strain rates and temperatures. Simulation results in terms of geometrical changes and microstructural evolution are compared with the experimental measurements. The model predicts anisotropic material flow resulting from the hard-to-deform crystallographic directions, the development of gradients in texture and twinning through the geometries, tension–compression asymmetry, as well as the extent of plasticity under impact.

42 ENGINEERING↗
DOE OSTI2023

ReaxFF reactive force field model enables accurate prediction of physiochemical and mechanical properties of crystalline and amorphous shape‐memory polyurethane

Abstract We investigate reactive force field (ReaxFF) prediction of physiochemical and mechanical properties, and stimuli response behavior of a shape‐memory polyurethane (SMPU) rigid segment. We used SMPUs as a platform to link the rigid/hard segment domain crystallinity with properties that can also be experimentally measured. Specifically, we developed models for crystalline and amorphous forms of a common rigid segment of SMPUs, for example, 4,4′‐diphenylmethane diisocyanate (MDI) with n‐butanediol (BDO), denoted as MDI‐BDO. ReaxFF molecular dynamics simulations with equilibrium dynamics are performed by controlling mass, temperature, and pressure or volume to predict the structural and physiochemical properties of crystalline and amorphous rigid segments of SMPU systems, followed by simulating their XRD and FTIR patterns, respectively. Non‐equilibrium ReaxFF simulations are implemented by using uniaxial box deformation technique to study response‐behavior of crystalline SMPU under tensile loading and during stress‐relaxation. The model and simulation scheme presented in this study demonstrates a robust theoretical pathway toward application‐dictated, on‐demand design of high‐performance SMPUs in applications such as soft robotics.

Afroz, Mohammad M.↗
DOE OSTI2024

End-Use Savings Shapes Measure Documentation: Variable Refrigerant Flow With 25% Upsizing Allowance for Heating

This documentation focuses on a single end-use savings shape upgrade - a variable refrigerant flow with heat recovery (VRF HR) heating and cooling system coupled with a dedicated outdoor air system (DOAS) for ventilation and where 25% upsizing (or 125% of the original size) is allowed for heating in colder climates (i.e., heating dominant regions). This document will primarily discuss the additional changes to the sizing algorithm and modeling approach, while a comprehensive overview of the fundamental modeling methodology and background of the VRF HR DOAS upgrade, including applicability and other key assumptions can be found in the original documentation: Variable Refrigerant Flow with Heat Recovery and Dedicated Outdoor Air System. To provide more context on the 25% upsizing algorithm, if the building is a cooling dominant (e.g., design cooling load higher than design heating load ), then the outdoor unit capacity of the VRF heat pump is sized based on the design cooling load. However, if the building is heating dominant, then this building is a candidate for the 25% upsizing allowance. Once the 25% upsizing is allowed, if the 25% upsized capacity (or 125% from the original size) represented with design condition exceeds the design heating load, then the design heating load is used for calculating the rated capacity of the outdoor unit. And if the 25% upsized capacity represented with the design condition does not exceed the design heating load, then the 25% upsized capacity represented with rated condition is used for the capacity of the outdoor unit while the remaining heating load is handled with the supplemental/backup electric resistance coil.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

The three-way junction structure of the HIV-1 PBS-segment binds host enzyme important for viral infectivity

HIV-1 reverse transcription initiates at the primer binding site (PBS) in the viral genomic RNA (gRNA). Although the structure of the PBS-segment undergoes substantial rearrangement upon tRNALys3 annealing, the proper folding of the PBS-segment during gRNA packaging is important as it ensures loading of beneficial host factors. DHX9/RNA helicase A (RHA) is recruited to gRNA to enhance the processivity of reverse transcriptase. Because the molecular details of the interactions have yet to be defined, we solved the solution structure of the PBS-segment preferentially bound by RHA. Evidence is provided that PBS-segment adopts a previously undefined adenosine-rich three-way junction structure encompassing the primer activation stem (PAS), tRNA-like element (TLE) and tRNA annealing arm. Disruption of the PBS-segment three-way junction structure diminished reverse transcription products and led to reduced viral infectivity. Because of the existence of the tRNA annealing arm, the TLE and PAS form a bent helical structure that undergoes shape-dependent recognition by RHA double-stranded RNA binding domain 1 (dsRBD1). Mutagenesis and phylogenetic analyses provide evidence for conservation of the PBS-segment three-way junction structure that is preferentially bound by RHA in support of efficient reverse transcription, the hallmark step of HIV-1 replication.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Multi‐Stage Optimization of Pore Size and Shape in Pore‐Space‐Partitioned Metal–Organic Frameworks for Highly Selective and Sensitive Benzene Capture

Abstract Compared to exploratory development of new structure types, pushing the limits of isoreticular synthesis on a high‐performance MOF platform may have higher probability of achieving targeted properties. Multi‐modular MOF platforms could offer even more opportunities by expanding the scope of isoreticular chemistry. However, navigating isoreticular chemistry towards best properties on a multi‐modular platform is challenging due to multiple interconnected pathways. Here on the multi‐modular pacs (partitioned acs) platform, we demonstrate accessibility to a new regime of pore geometry using two independently adjustable modules (framework‐forming module 1 and pore‐partitioning module 2). A series of new pacs materials have been made. Benzene/cyclohexane selectivity is tuned, progressively, from 4.5 to 15.6 to 195.4 and to 482.5 by pushing the boundary of the pacs platform towards the smallest modules known so far. The exceptional stability of these materials in retaining both porosity and single crystallinity enables single‐crystal diffraction studies of different crystal forms (as‐synthesized, activated, guest‐loaded) that help reveal the mechanistic aspects of adsorption in pacs materials.

Chen, Yichong↗
DOE OSTI2020

Geochemistry and Microbiology Predict Environmental Niches with Conditions Favoring Potential Microbial Activity in the Bakken Shale

The Bakken Shale and underlying Three Forks Formation is an important oil and gas reservoir in the United States. The hydrocarbon resources in this region are accessible using unconventional oil and gas extraction methods, including horizontal drilling and hydraulic fracturing. However, the geochemistry and microbiology of this region are not well understood, although they are known to have major implications for productivity and water management. In this study, we analyzed the produced water from 14 unconventional wells in the Bakken Shale using geochemical measurements, quantitative PCR (qPCR), and 16S rRNA gene sequencing with the overall goal of understanding the complex dynamics present in hydraulically fractured wells. Bakken Shale produced waters from this study exhibit high measurements of total dissolved solids (TDS). These conditions inhibit microbial growth, such that all samples had low microbial loads except for one sample (well 11), which had higher 16S rRNA gene copies. Our produced water samples had elevated chloride concentrations typical of other Bakken waters. However, they also contained a sulfate concentration trend that suggested higher occurrence of sulfate reduction, especially in wells 11 and 18. The unique geochemistry and microbial loads recorded for wells 11 and 18 suggest that the heterogeneous nature of the producing formation can provide environmental niches with conditions conducive for microbial growth. This was supported by strong correlations between the produced water microbial community and the associated geochemical parameters including sodium, chloride, and sulfate concentrations. The produced water microbial community was dominated by 19 bacterial families, all of which have previously been associated with hydrocarbon-reservoirs. These families include Halanaerobiaceae, Pseudomonadaceae, and Desulfohalobiaceae which are often associated with thiosulfate reduction, biofilm production, and sulfate reduction, respectively. Notably, well 11 was dominated by sulfate reducers. Our findings expand the current understanding of microbial life in oil and gas related environments, provide new insights into how the unique produced water conditions shape microbial communities, and demonstrate thatproduced water chemistry is tightly linked with microbiota in the Bakken Shale. Finally, our study shows that additional research efforts are necessary to understand this ecosystem and highlights the need for work in other unconventional reservoirs.

Tinker, Kara↗
DOE OSTI2025

Strain-Based Assessment of Shale Caprock during Cyclic Underground Hydrogen Storage

Successful large-scale underground hydrogen storage (UHS) in depleted gas reservoirs depends on the integrity of the overlying caprock to prevent hydrogen loss during cyclic injection and depletion. Prior studies on crushed Marcellus shale, a potential caprock, indicate that cyclic hydrogen injection and depletion induces microstructural changes, increasing porosity and permeability. However, the extent of these changes in intact shale remains unclear. This study presents a strain-based experimental approach to quantify volumetric strain evolution in intact Marcellus shale matrix under unconstrained stress conditions. A quadrant-shaped shale sample without visible fractures underwent eight hydrostatic pore-pressure cycles (injection to 1500 psi and depletion to 500 psi in 250 psi steps). Linear strain gauges measured strain in three orthogonal directions. Results indicate progressive plastic strain accumulation, leading to an ∼12% increase in matrix porosity after eight cycles, with an estimated 19% increase after 30 cycles. This porosity increase follows a logarithmic trend, suggesting a diminishing effect in later cycles. Additionally, permeability and diffusive mass flux are projected to rise by ∼70% over 30 cycles, enhancing hydrogen migration risk. The shale matrix also exhibited mechanical stiffening over successive cycles, limiting large-scale deformation but not preventing porosity enhancement. A new parameter, α, was introduced to characterize shale sensitivity to cyclic loading, aiding UHS caprock assessments. These findings underscore the necessity of incorporating cyclic loading effects in UHS site selection and operational strategies to ensure long-term storage integrity.

08 HYDROGEN↗
DOE OSTI2022

Multiphysics Analyses of the Bottom Components of the 3D Printed Transformational Challenge Reactor

This research presents multi-physics analyses on the bottom components of the Transformational Challenge Reactor (TCR) facility. These components include the bottom axial reflector and the steel exit cone. The bottom axial reflector is made of pure silicon carbide elements hosting helium cooling channels These elements are 3D printed and therefore can host any arbitrary shape of the helium cooling channels. The design of the bottom reflector considers the neutronics and thermo-fluid dynamics performances as well as the manufacturing process optimization. More precisely, the best design of the bottom reflector reduces neutron leakage by avoiding straight cylindrical helium channels that facilitate neutron leakage, minimizes the helium flow pressure drop, and reduces the number of 3D printed silicon carbide pieces. The exit cone steel structure collects the hot helium from the bottom fuel assemblies and channels the cold helium to the top of the fuel assemblies. The steel simultaneous contact with hot and cold helium flows sets a large thermal gradient. Different designs of the exit cone are proposed to reduce the steel equivalent stress from the helium thermal load. The multi-physics analyses have been performed using Ansys Fluent, Ansys Mechanical, STAR-CCM+, and Serpent computer programs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗