SEARCH · Engineering Papers

Results for “variational fracture modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Mesoscale simulations of pressure-shear loading of granular tungsten carbide

Numerical simulations of pressure-shear loading of a granular material are performed using the shock physics code CTH. A simple mesoscale model for the granular material is used that consists of a randomly packed arrangement of solid circular or spherical grains of uniform size separated by vacuum. The grain material is described by a simple shock equation of state, elastic perfectly plastic strength model, and fracture model with baseline parameters for WC taken from previous mesoscale modeling work. Simulations using the baseline material parameters are performed at the same initial conditions of pressure-shear experiments on dry WC powders. Except for some localized flow regions appearing in simulations with an approximate treatment of sliding interfaces among grains, the samples respond elastically during shear, which is in contrast to experimental observations. By extending the simulations to higher shear wave amplitudes, macroscopic shear failure of the simulated samples is observed with the shear strength increasing with increasing stress confinement. The shear strength is also found to be strongly dependent on the grain interface treatment and on the fracture stress of the grains, though the variation in shear strength due to fracture stress decreases with increasing stress confinement. At partial compactions, the transverse velocity histories show strain-hardening behavior followed by formation of a shear interface that extends through the transverse dimensions of the sample. Near full compaction, no strain hardening is observed and, instead, the sample transitions sharply from an elastic response to formation of an internal shear interface. Agreement with experiment is shown to worsen with increasing confinement stress with simulations overpredicting the shear strengths measured in experiment. The source of the disagreement can be ultimately attributed to the Eulerian nature of the simulations, which do not treat contact and fracture realistically.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Completion Effects on Diagnosing Multistage Facture Treatments with Distributed Temperature Sensing

Summary Distributed temperature sensing (DTS) is a valuable tool to diagnose multistage hydraulic fracture treatments. When a stage interval is shut in, the clusters that take more fluid during pumping warm up more slowly. Therefore, the fluid volume injected into each cluster can be quantitatively interpreted by numerical inversion of the warm-back temperature behavior. This general concept assumes that the different warm-back behavior is controlled by only the injected fluid volume; however, recent observations of DTS data indicate that completion configurations significantly influence the warm-back behavior. This paper investigates the completion effects on the DTS interpretation. In ideal conditions, when a stage is fractured, the upstream stage intervals should show an almost uniform temperature that is close to the injected fluid temperature. This is due to the high fluid velocity of injected fluid in the wellbore, and the upstream intervals have not been perforated (noncommunicating intervals), so the only heat transfer is heat conduction between the wellbore fluid and the surrounding reservoir. But the field DTS data show considerably irregular variations in temperature along the upstream stage intervals. These variations are caused by the completion effects. The nonuniform temperature profile is caused by different heat transfer behavior induced by completion hardware along the production casing string, such as joints, clamps, and blast protectors, and by the sensing cable location in the cement, as well as the cement quality. Because the varying heat transfer behavior impacts the warm-back behavior as well as the temperature profile, the completion effects need to be considered in DTS interpretation. A method of DTS interpretation considering the completion effects to diagnose multistage fracture treatments was developed. Because the heat transfer between a wellbore and a reservoir depends on the overall heat transfer coefficient describing heat conduction through the completion in a forward model, this parameter needs to be tuned all along the wellbore. To calibrate the completion effect, the temperature inversion is conducted using the temperature measured at a stage interval that is upstream of a stage interval currently being treated. Because the interpreted stage interval is not perforated at that time, the thermal behavior at the noncommunicating interval is governed by only the heat conduction through the completion environment. Once the effective values of the overall heat transfer coefficient are estimated along the interpreted stage interval, they can be assumed to be constant physical parameters. Then, the fluid volume distribution is interpreted by using the effective overall heat transfer coefficient profile along each interval. This study provides a field application of the developed interpretation method. The new interpretation method provides more accurate diagnosis of fracture treatments by DTS interpretation.

Engineering↗

DOE OSTI2024

Variations in Bedrock and Vegetation Cover Modulate Subsurface Water Flow Dynamics of a Mountainous Hillslope

Abstract Predicting the hydrological response of watersheds to climate disturbances requires a detailed understanding of the processes connecting hillslopes and streams. Using a network of soil moisture and temperature sensors, electrical resistivity tomography monitoring, and a weather station we assess the above and below‐ground processes driving the hydrological response of a hillslope during snowmelt and summer monsoon. The transect covers bedrock and vegetation gradients, with a steep upper part characterized by shallow bedrock, and gentle lower part underlain by colluvium. The main vegetation cover is conifers on the upper, and grass and veratrum on the lower part. Combined with a simplified hydrological model, we show that the thin soil layer of the steep slope acts as a preferential flow path, leading to mostly shallow lateral flow, interrupted by vertical flow, mostly at tree locations, and likely facilitated by flow along fractures and roots. Vertical flow and upstream‐driven groundwater dynamics are prevailing at the colluvium, presenting a very different hydrological behavior compared to the upper part. These results show that subsurface structure and features have a strong control on the hydrological response of a hillslope and that those can create considerably varying hydrological dynamics across small spatial scales.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Mesoscale Modeling for Restructuring and Fragmentation in High Burnup UO 2

This report summarizes the mesoscale modeling conducted in fiscal year 2025 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, focusing on the microstructural evolution and restructuring in high burnup UO 2 nuclear fuel and its impact on fuel fragmentation. We developed a pioneering phase-field model to simulate restructuring behavior across different regions of high burnup fuel, including the dark zone and rim region. A grand-potential-based phase-field model is employed to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of temperature and burnup rate were studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to initiate around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. The restructuring occurs faster with higher burnup rate. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. The predicted microstructures are consistent with experimental observations of the restructured regions. Based on the mesoscale simulations, a mechanistic model for restructuring and grain size evolution was implemented in BISON. Thus, this work provides a first-of-its-kind restructuring model for different regions of high-burnup fuel to BISON. It was shown that the model predicts appropriate grain size evolution along fuel radius as those observed in experiments. This model enables BISON to predict the effect of restructuring on the fission gas release. Finally, the phase-field fracture simulation with dark zone specific microstructures were presented to provide the fragmentation criteria for the dark zone. This work introduces a first-of-its-kind restructuring model for high burnup fuel, significantly enhancing BISON's predictive capabilities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

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

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

58 GEOSCIENCES↗

DOE OSTI2025

Orientation-Dependent Nanomechanical Behavior of Pentaerythritol Tetranitrate as Probed by Multiple Nanoindentation Tip Geometries

Nanoindentation can be leveraged to aid in the high fidelity modeling of dislocation mediated plasticity in pentaerythritol tetranitrate (PETN), an anisotropic energetic molecular crystal. Moreover, nanoindentation tip parameters such as tip geometry, size, and degree of acuity can be utilized to target anisotropic behavior. In this work, nanoindentation was conducted across a range of orientations on the (110) face of PETN to characterize resultant yield behavior, mechanical property measurements, and resultant slip behavior and fracture initiation. Three different indentation tips were utilized: a 3-sided pyramidal Berkovich tip, a 4-sided high aspect ratio Knoop tip, and a 90° conical tip. Ultimately, indenter tip radius was documented to impact yield behavior, whereas tip geometry affected larger scale processes such as slip, and tip acuity was the dominating factor that led to fracture. The axisymmetric conical tip, serving as a baseline, showed the least amount of variation in mechanical property measurements but also the largest distribution of maximum shear stress at which initial yielding occurred. Its high degree of acuity, however, was more prone to induce fracture at higher loads. The Knoop tip was shown to be suitable for average measurements, but also for elucidation of certain anisotropic features. A distinctly higher perceived hardness at 45° was measured with the Knoop tip, indicating less dislocation motion in that direction also observed in this work via scanning probe microscopy. Lastly, the commonly used Berkovich tip was a good compromise whereby it provided a representative volume element describing the average behavior of the material. These results can be utilized to target desired anisotropic behavior in a wider range of molecular crystals, as well as to inform theoretical considerations for dislocation mediated plasticity in PETN.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Experimental Investigation of Low-Frequency Distributed Acoustic Sensor Responses to Two Parallel Propagating Fractures

Low-frequency distributed acoustic sensing (LF-DAS) is a diagnostic tool for hydraulic fracture propagation with far-field monitoring using fiber optic sensors. LF-DAS senses strain rate variation caused by stress field change due to fracture propagation. Fiber optic sensors are installed in the monitoring wells in the vicinity of a fractured well. From the strain responses, fracture propagation can be evaluated. To understand subsurface conditions with multiple propagating fractures, a laboratory-scale hydraulic fracture experiment was performed simulating the LF-DAS response to fracture propagation with embedded distributed optical fiber strain sensors under these conditions. The experiment was performed using a transparent cube of epoxy with two parallel radial initial flaws centered in the cube. Fluid was injected into the sample to generate fractures along the initial flaws. The experiment used distributed high-definition fiber optic strain sensors with tight spatial resolutions. The sensors were embedded at two different locations on opposite sides of the initial flaws, serving as observation/monitoring locations. We also employed finite element modeling to numerically solve the linear elastic equations of equilibrium continuity and stress–strain relationships. The measured strains from the experiment were compared to simulation results from the finite element model. The experimentally derived strain and strain-rate waterfall plots from this study show the responses to both fractures propagating, while the fracture at the lower position took most of the fluid during the experiment. Interestingly, a fracture first began propagating from the upper flaw of the two flaws, but once the lower fracture was initiated, it grew much faster than the upper fracture. Both fibers were intercepted by the lower fracture, further verifying the strain signature as a fracture is approaching and intersecting an offset fiber.

Chemistry↗

DOE OSTI2022

Polymorphism of garnet solid electrolytes and its implications for grain-level chemo-mechanics

Understanding and mitigating filament formation, short-circuit and solid electrolyte fracture is necessary for advanced all-solid-state batteries. Here, we employ a coupled far-field high-energy diffraction microscopy and tomography approach for assessing the chemo-mechanical behaviour for dense, polycrystalline garnet (Li 7 La 3 Zr 2 O 12 ) solid electrolytes with grain-level resolution. In situ monitoring of grain-level stress responses reveals that the failure mechanism is stochastic and affected by local microstructural heterogeneity. In this work, coupling high-energy X-ray diffraction and far-field high-energy diffraction microscopy measurements reveals the presence of phase heterogeneity that can alter local chemo-mechanics within the bulk solid electrolyte. These local regions are proposed to be regions with the presence of a cubic polymorph of LLZO, potentially arising from local dopant concentration variation. The coupled tomography and FF-HEDM experiments are combined with transport and mechanics modelling to illustrate the degradation of polycrystalline garnet solid electrolytes. The results showcase the pathways for processing high-performing solid-state batteries.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Stratigraphically Controlled Stress Variations at the Hydraulic Fracture Test Site-1 in the Midland Basin, TX

We investigated the relationship between stratigraphy, stress, and microseismicity at the Hydraulic Fracture Test Site-1. The site comprises two sets of horizontal wells in the Wolfcamp shale and a deviated well drilled after hydraulic fracturing. Regional stresses indicate normal/strike-slip faulting with E-W compression. Stress measurements in vertical and horizontal wells show that the minimum principal stress varies with depth. Strata with high clay and organic content show high values of the least compressive stress, consistent with the theory of viscous stress relaxation. By integrating data from core, logs, and the hydraulic fracturing stages, we constructed a stress profile for the Wolfcamp sequence, which predicts how much pressure is required for hydraulic fracture growth. We applied the results to fracture orientation data from image logs to determine the population of pre-existing faults that are expected to slip during stimulation. We also determined microseismic focal plane mechanisms and found slip on steeply dipping planes striking NW, consistent with the orientations of potentially active faults predicted by the stress model. This case study represents a general approach for integrating stress measurements and rock properties to predict hydraulic fracture growth and the characteristics of injection-induced microseismicity.

microseismicity↗

DOE OSTI2023

Utah FORGE: Optimization of a Plug-and-Perf Stimulation (Fervo Energy)

Information around the plug-and-perf treatment design at Utah FORGE by Fervo Energy. Objective and Purpose: - Develop a multistage hydraulic stimulation approach designed specifically to target the top three factors that control the technical and commercial viability of an EGS system: i) Achieving sufficient injectivity to support high cross-well flow rates ii) Distributing flow evenly across the wellbore and reservoir to maximize heat mining efficiency, ensure sustained heat transfer, and mitigate thermal breakthrough iii) Overcoming the effects of stress heterogeneity, stress shadowing, and variations in natural fracture properties during the stimulation treatment, leading to a more predictable stimulated reservoir volume and offset well placement - The following activities will be performed: i) Design, plan, and execute a multistage plug-and-perf stimulation treatment at a Fervo site with data acquisition and well testing activities aimed at addressing key technical aspects of the issues above ii) Perform data processing and interpretation of field results to translate the results form the Fervo site to a site-specific design at the Utah FORGE site iii) Design, plan, and execute a multistage plug-and-perf stimulation treatment design at the Utah FORGE site Methods and Approach: - Design a detailed data acquisition plan to maximize learning around: i) DFIT testing ii) Petrophysical logging, image logging iii) Permanent DAS/DTS fiber optic monitoring iv) Deep borehole microseismic monitoring v) Shallow borehole induced seismicity monitoring vi) Injection/production testing (RTA analysis, tracer testing) vii) Integrated numerical modeling and production forecasting

15 GEOTHERMAL ENERGY↗

DOE OSTI2022

Effect of lot microstructure variations on detonation performance of the triaminotrinitrobenzene (TATB)-Based insensitive high explosive PBX 9502

PBX 9502 is an important insensitive high explosive due to its combination of safety properties and detonation performance. It is a polymer-bonded formulation consisting of 95 wt.% 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as the high explosive crystal, bound with Kel F-800 (FK-800), a co-polymer of chlorotrifluoroethylene and vinylidene-fluoride. Two different types are used, one known as virgin that uses only pristine manufactured TATB, and a second known as recycled that has 50 wt.% of its TATB reclaimed from machining scraps of previously pressed virgin PBX 9502. Recycled lots have a higher percentage of fine particles compared to virgin lots, due to the fracturing and damage sustained by TATB crystals during pressing. We examine the influence that this TATB microstructure difference between virgin and recycled PBX 9502 has on detonation performance properties. New rate-stick geometry diameter effect, detonation front shape and cylinder expansion test data are obtained for two previously uncharacterized virgin lots and one recycled lot of PBX 9502. This is combined with previously published data for one virgin and one recycled lot for evaluation purposes. Detonation shock dynamics model calibrations are conducted on each of the five lots to provide an assessment of the detonation timing characteristics of virgin versus recycled PBX 9502 lots. For two of the virgin lots, detonations propagate slower in the rate-stick geometry than those in the two recycled lots. However, the other virgin lot tested has propagation rates comparable to that of the recycled lots for larger diameter rate-sticks. The latter result, though, is shown to be geometry dependent and depends on the range of detonation curvatures accessed in different geometries. New copper-confined, cylinder expansion tests are conducted on each of the three virgin and two recycled lots to obtain detonation product Jones-Wilkins-Lee equations of state, enabling an assessment of the metal push capabilities for each of the five lots. We find that the metal push capabilities, characterized by the evolution of the heat of detonation with volume, are similar between the virgin and recycled lots. Thus, changes in microstructure between different PBX 9502 lots seemingly affect the rate of reaction, but not the overall energy content.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

University of Hawai‘i, Shallow Geothermal Resources: Energy Technology Innovation Partnership Project (Final Report)

Scientists at Lawrence Berkeley National Laboratory (Berkeley Lab) have teamed up with the University of Hawai‘i at Manoa (UH Manoa) through the U.S. Department of Energy’s Energy Technology Innovation Partnership Project to evaluate the technological and market feasibility of shallow geothermal heat exchanger (GHE) technology. UH requested this analysis to evaluate opportunities in building cooling, energy efficiency, and emissions reduction applications in Hawai‘i. UH has an abundance of geologic and geothermal data and is looking to the national labs’ expertise to execute this analysis. UH is also interested in investigating policy, regulatory, and business conditions advantageous for implementation of a pilot project and more broad deployment of this technology in Hawai‘i. In many locations around the world, the demands for heating and cooling are roughly balanced over the course of the year, so GHEs do not cause significant long-term changes in subsurface temperature. This is not the case in Hawai’i, where the demand for heating is very small, meaning that, over time, GHEs will add heat to the subsurface. If temperatures increase significantly, GHE systems will not work as designed. Regional groundwater flow has the potential to sweep heated water away from boreholes, thereby maintaining the functionality of the GHE system. Significant regional groundwater flow requires two things: a sufficiently large driving hydraulic head gradient (usually closely related to surface topography), and sufficient porosity and permeability to enable groundwater to flow in large enough quantities to enable near-borehole temperatures to be maintained at ambient values. Hawai‘i’s volcanic terrain offers ample surface topographic variation. The lava itself shows an extremely large range of porosity and permeability, so sites with large enough values of these properties must be selected. Numerical modeling of coupled groundwater and heat flow can be used to determine how large is large enough. Primarily, closed-loop systems have been investigated. Other options considered are open-loop systems and using cool seawater as the chilling source. Project work investigated the feasibility of GHE technology at two scales. At the island scale, GIS layers of various attributes relevant for GHE were combined to develop an overall favorability map for employing GHE in Hawai‘i. At the local scale, a hydrogeologic model for the subsurface component of a closed-loop system was developed for the Stan Sheriff Center at the UH Manoa campus. This site is considered promising because the rock below and immediately downgradient of the borefield is highly permeable, consisting of a subsurface karst system (limestone containing high-permeability open channels), which is underlain by a thick, high-permeability fractured basalt. Moreover, the site is near the base of the Ko‘olau Range, providing a large hydraulic head gradient. Thus, groundwater flow through the site is expected to be large, enabling efficient removal of heated groundwater. A full-GHE-system model of the site was also developed, with a simplified representation of the subsurface, in which groundwater flow is not considered and heat transfer is purely by conduction. Using the building cooling load data provided by UH, simulation results show that with groundwater flow present, a GHE can operate successfully for at least 10 years, but with no groundwater flow, the subsurface begins to heat up after only one year of operation, making the GHE unviable within 2-6 years. The team also developed a techno-economic model for this site to compare the cost of cooling using a GHE system with the costs of operating the current air-conditioning system. The GHE system is advantageous economically if favorable tax incentives and interest rates can be obtained.

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Microstructural Engineering and Accelerated Test Method Development to Achieve Low Cost, High Performance Solutions for Hydrogen Storage and Delivery

This project made advancements in developing lower cost steel alloys with novel microstructural design for use in hydrogen refueling infrastructure such as storage, compressors, and dispensing components, and utilizing accelerated test methods to efficiently evaluate variations in alloy and microstructure design. The project specifically sought to design alloys with lower nickel contents to reduce alloy cost, which was accomplished through substituting manganese for nickel along with other alloy additions to control deformation characteristics known to be important for hydrogen embrittlement resistance. Through Mn substitutions for nickel, austenitic and duplex austenite-ferrite steels were successfully developed with lower cost than currently available commercial stainless steel products that are employed for hydrogen refueling infrastructure. The steels were processed to achieve comparable strength and toughness in hydrogen environments as the commercially available steels containing high Ni contents. To evaluate mechanical performance in hydrogen, a testing methodology was employed to compare ubiquitous laboratory testing using electrochemical hydrogen charging in a liquid electrolyte to less accessible high pressure gaseous testing. While the application of these steels is in high pressure gaseous environments, the electrochemical hydrogen charging tests produced comparable results. Additionally, the Los Alamos Neutron Scatting Center enabled characterization of deformation mechanisms of the steel alloys in the presence of hydrogen, which has been associated with steel alloy characteristics associated with hydrogen embrittlement. Finally, a fracture mechanics based test bed model was developed to predict the influence of hydrogen gas pressure and fatigue conditions on fatigue lifetimes of pressure vessel steels. Together, these developments can be employed to enable lower cost hydrogen fueling infrastructure and more reliable prediction of steel alloy components in hydrogen service conditions. In particular, the newly alloys have the potential to replace stainless steels, having demonstrated comparable performance at substantially reduced cost.

08 HYDROGEN↗

DOE OSTI2022

Irradiation damage concurrent challenges with RAFM and ODS steels for fusion reactor first-wall/blanket: a review

Reduced activation ferritic martensitic (RAFM) and oxide dispersion strengthened (ODS) steels are the most promising candidates for fusion first-wall/blanket (FW/B) structures. The performance of these steels will deteriorate during service due to neutron damage and transmutation-induced gases, such as helium/hydrogen, at elevated operating temperatures. Here, after highlighting the operating conditions of fusion reactor concepts and a brief overview, the main irradiation-induced degradation challenges associated with RAFM/ODS steels are discussed. Their long-term degradation scenarios such as (a) low-temperature hardening embrittlement (LTHE)—including dose-temperature dependent yield stress, tensile elongations, necking ductility, test temperature effect on hardening, Charpy impact ductile-to-brittle transition temperature and fracture toughness, (b) intermediate temperature cavity swelling, (c) the effect of helium on LTHE and cavity swelling, (d) irradiation creep and (e) tritium management issues are reviewed. The potential causes of LTHE are discussed, which highlights the need for advanced characterisation techniques. The mechanical properties, including the tensile/Charpy impact of RAFM and ODS steels, are compared to show that the current generation of ODS steels also suffers from LTHE, and shows irradiation hardening up to high temperatures of ~400 °C–500 °C. To minimise this, future ODS steel development for FW/B-specific application should target materials with a lower Cr concentration (to minimise α'), and minimise other elements that could form embrittling phases under irradiation. RAFM steel-designing activities targeting improvements in creep and LTHE are reviewed. The need to better understand the synergistic effects of helium on the thermo-mechanical properties in the entire temperature range of FW/B is highlighted. Because fusion operating conditions will be complex, including stresses due to the magnetic field, primary loads like coolant pressure, secondary loads from thermal gradients, and due to spatial variation in damage levels and gas production rates, an experimentally validated multiscale modelling approach is suggested as a pathway to future reactor component designing such as for the fusion neutron science facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Micro-tensile characteristics of As-fabricated and irradiated AGR-2 TRISO fuel particle buffer, IPyC, and buffer-IPyC interlayer regions

A recently developed micro-tensile sample preparation technique was implemented to evaluate the tensile strengths of the buffer, IPyC, and buffer-IPyC interlayer regions of the unirradiated and irradiated AGR-2 TRISO fuel particles. Understanding the mechanical properties of the buffer-IPyC interlayer is essential for developing thermomechanical models of buffer-IPyC separation, yet there is a lack of experimental data on its micro-tensile properties. TEM analysis was conducted on these regions to determine the microstructural changes relevant to the samples' tensile properties. In the unirradiated TRISO particle samples, the buffer layer demonstrated the weakest tensile strength, while the IPyC layer exhibited the highest. Conversely, in the irradiated TRISO particle samples, the buffer-IPyC interlayer region showed the lowest tensile strength, with the IPyC layer being the strongest. Fractures in the samples from the buffer-IPyC region predominantly occurred either in the buffer layer or at the buffer-IPyC interface. However, some buffer-IPyC interlayer samples displayed stress-strain and fracture behaviors more akin to the IPyC layer than the buffer layer. Analysis of diffraction patterns suggests that irradiation may have increased anisotropy in the three regions tested. Despite this suggested increase in anisotropy, there was no evidence that it affected the measured strengths. The irradiated TRISO particles demonstrated a considerable increase in void space and a decrease in ultimate tensile strength within the buffer-IPyC interlayer region due to the densification and contraction of the buffer layer. Minor variations in diffraction ring patterns were also observed. These changes, coupled with a significant reduction in the Weibull modulus/shape parameter, imply that irradiation-induced densification leads to tearing between the buffer and IPyC layers at locations of elevated porosity in the buffer-IPyC interlayer region.

Tristructural isotropic (TRISO)↗

DOE OSTI2022

HFTS-1 Natural Joint Data and Engineering Summary

This report provides an analysis of engineering and geologic data collected from the Hydraulic Fracturing Test Site #1 (HTSF-1) project in the southern Midland Basin, Reagan County, Texas. The site is being studied as part of the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative at the National Energy Technology Laboratory (NETL). The data collected is intended to provide a basis of understanding of the site and to construct simulation models. This report provides a summary analysis of the engineering and geologic data collected from the project to provide a basis of understanding of the site and to construct simulation models. The report examines various possible correlations in engineering properties based on natural fracture data from the program, which consisted of 11 horizontal wells, one vertical well, and one slant well. The horizontal wells are in two horizons: 1) Upper and 2) Middle Wolfcamp formations. Program data include: fracture frequency and fracture orientation data from four core runs in a slant well; fracture spacing and orientation from a vertical pilot well; laboratory triaxial testing and mineralogical determinations; and porosity results from magnetic resonance analyses from various wells, together with observations based on the data collection. In addition, available references were reviewed on the site for additional insights. As the focus of the report is on the natural system, hydraulic fracture data from the site were not examined in detail in this report. Data variability is the chief observation in examination of the database. Fracture frequency in the slant well can range from sections with values as high as five fractures per ft to sections up to 100+ ft in length with no natural observed fractures. Fracture spacing across is typically less than 10 ft, but can range up to hundreds of feet. Apparent fracturing shows the trends in two predominate orientations, E-SW and WNW-ESE, but minor variations exist. The rock units vary across the site from siliceous mudstones to calcareous mudstones, showing a general layering with depth. The laboratory properties such as strength and modulus show no apparent trend with depth, but appear to correlate with rock mineralogy with high strength and modulus values where calcium content is high. In addition, an attempt to examine variability and mineralogy on a larger scale was made using a color-coded system based on gamma ray measurements. A staged colored approach was adopted, presuming that lower gamma ray values indicate higher value of calcium content (blue scale) and that higher gamma ray values indicate higher clay mineral content (orange scale). As provided in report appendices, the system correlated well with visual examination of the slant core and the petrofabric analyses of the vertical pilot well. The results showed large variability in mineral content along the horizontal plane across the site. The change in mineral content was also rapid, on a scale less than that of the average hydraulic fracture stage length of about 180 ft.

42 ENGINEERING↗

DOE OSTI2025

LATTE: open-source, high-performance traveltime computation, tomography and source location in acoustic and elastic media

Traveltime-based tomography and source location are fundamental approaches for imaging subsurface structures and understanding the spatiotemporal distribution of seismicity from local to global scales. We present an open-source, high-performance framework integrating eikonal equation solvers and adjoint-state theory for traveltime computation, velocity tomography, source location and joint tomography-location in 2-D/3-D acoustic and elastic media. We introduce novel regularization schemes based on total generalized p-variation, structural similarity and multitask machine learning to enhance the fidelity and interpretability of inverted models and source locations. Key features of our implementation also include the ability to leverage both absolute-difference and double-difference traveltime misfits for high-fidelity velocity tomography and source parameter estimation; support for traveltime computation and inversion in diverse 2-D/3-D scenarios with arbitrary source and receiver distributions; and a perturbation-based optimal step-size estimation method to reduce computational costs. In addition, our implementation employs shared-memory and distributed-memory parallelization to provide an efficient solution for traveltime computation, tomography, and source location. In conclusion, we validate the efficacy and accuracy of our approach through multiple synthetic data examples.

58 GEOSCIENCES↗

DOE OSTI2020

Molecular scale insight of pore morphology relation with mechanical properties of amorphous silica using ReaxFF

Porous materials are typically heterogeneous and they contain large variations of micro-/nano-pore structures, causing complicated behaviors. In continuum models, most mechanical properties of porous materials are estimated based on porosity, while the variations of micro/nano structures are ignored. That could be problematic as the microscopic heterogeneity may affect the mechanical response of porous materials. Thus, understanding micro/nano heterogeneity impact has been the focus in many scientific and engineering subjects. In the present study, we investigated the effect of nanopore structure (including pore shape and orientation) as well as porosity on mechanical properties of amorphous silica (a-SiO 2 ). The pore sizes in our simulations are comparable to the corresponding ones observed in a-SiO 2 based materials. We found that the existing of nanopores strongly influences Young’s modulus (E) and critical energy release rate (G IC ). These properties decrease with increasing porosity. Importantly, the impact of nanopores was characterized by structural parameters of porous materials. In addition to dependency on porosity, Young’s modulus also was found to vary as a function of potential energy per atom, which highly depends on nanopore shape. Furthermore, critical energy release rate was found to increase with increasing ligament length (also known as pore wall thickness). Finally, the results highlighted the importance of nanopore structures, which must be taken into account when studying fracture mechanisms in porous materials. Based on our findings, it was proposed that mechanical properties of porous materials can be controlled by nano-engineering pore structures.

36 MATERIALS SCIENCE↗