Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS MEASUREMENT”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Versatile Liquid Crystal Elastomer Formulations Using Amine-Acrylate Chemistry and Processing for Advanced Manufacturing

This study aims to tune rheological properties of liquid crystal elastomers (LCEs) through a strategic approach, leveraging the differing reaction rates of thiols and primary and secondary amines with acrylate monomers or the inclusion of processing additives to enable amenability to advanced manufacturing. We found that varying the time of oligomerization at elevated temperatures of amine-functional monomers with acrylate-functional monomers enabled simple control of rheological properties. Additionally, we could also control rheological properties by introducing reinforcing fillers or organic solvents. By varying the oligomerization conditions, we achieved a broad range of viscosities, from 1.1 Pa s to 560 Pa s as measured at 60 °C. When the rheological properties were varied through the incorporation of various solid and liquid inclusions, viscosities could range from less than 1 Pa s to nearly 1000 Pa s when measured at 60 °C. The shape morphing capabilities that are intrinsic to LCEs and the mechanical properties were also characterized. Here, we find that the oligomerization conditions and the addition of fillers influence shape morphing and energy absorption, as characterized by the stress–strain characteristics, further suggesting their potential to be used in energy dissipation applications where performance can be turned. Last, we demonstrate how the range of accessible rheological properties grants formulations amenable to various advanced manufacturing techniques.

Liquid chromatography↗
DOE OSTI2025

Isomeric yield ratios of fission products: A missing piece in reactor antineutrino summation calculations

The calculation of the spectrum of antineutrinos ($\bar{v}_e$) from a reactor is a complicated problem requiring several nuclear data and physics inputs. Many of these have been investigated thoroughly to improve calculations and properly account for uncertainties. One input which has heretofore escaped consideration is the fission-yield distribution between ground and isomeric states. Here, in this work, we explore the effect of incorporating newly evaluated isomeric yield ratios (IYR) for 43 fission products into summation calculations and identify the disproportionate importance of certain isotopes, particularly at higher energies. Our analysis indicates that updated IYRs contribute to a significant increase in the $\bar{v}_e$ spectrum around and above 7 MeV, with increases of more than 50% at higher energies. Through a detailed sensitivity study, we highlight a number of isotopes for which the IYR has a substantial effect on the $\bar{v}_e$ spectrum. This work stresses the critical role of isomeric yields in calculations of reactor $\bar{v}_e$ spectra and points to the necessity for their accurate experimental determination, especially for key fission products, in order to refine our understanding and address the observed anomalies between measured and calculated $\bar{v}_e$ spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Mechanisms for Microseismicity Occurrence Due to CO 2 Injection at Decatur, Illinois: A Coupled Multiphase Flow and Geomechanics Perspective

Here, we numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO 2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologically consistent model that honors key stratigraphic horizons and 3D fault surfaces interpreted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO 2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO 2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection interval to the basement is essential to explain the destabilization of the regions where microseismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO 2 sequestration operations.

58 GEOSCIENCES↗
DOE OSTI2025

Isomeric Yield Ratios of fission products: a missing piece in reactor antineutrino summation calculations

The calculation of the spectrum of antineutrinos ( $\overline{V}$ e ) from a reactor is a complicated problem requiring several nuclear data and physics inputs. Many of these have been investigated thoroughly to improve calculations and properly account for uncertainties. One input which has heretofore escaped consideration is the fission yield distribution between ground and isomeric states. In this work, we explore the effect of incorporating newly evaluated isomeric yield ratios (IYR) for 43 fission products into summation calculations and identify the disproportionate importance of certain isotopes, particularly at higher energies. Our analysis indicates that updated IYRs contribute to a significant increase in the $\overline{V}$ e spectrum around and above 7 MeV, with increases of more than 50% at higher energies. Through a detailed sensitivity study, we highlight a number of isotopes for which the IYR has a substantial effect on the $\overline{V}$ e spectrum. This work stresses the critical role of isomeric yields in calculations of reactor $\overline{V}$ e spectra and points to the necessity for their accurate experimental determination, especially for key fission products, in order to refine our understanding and address the observed anomalies between measured and calculated $\overline{V}$ e spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Materials Characterization: A Primer for Solid Phase Processing Applications

The Pacific Northwest National Laboratory (PNNL) undertook the Materials Characterization, Prediction, and Control (MCPC) Laboratory Directed Research and Development (LDRD) Project to advance understanding of nuclear material processing and enable multifold acceleration in the development and qualification of new material systems produced via advanced manufacturing methods, such as solid phase processing, for use in national security and advanced energy applications (Smith 2021). As a two-year LDRD investment requiring focused research, the MCPC project applied only a subset of the wide range of available destructive and nondestructive characterization methods to provide data to the predictive modeling and data analytics tasks. The purpose of this report is to review a wide range of destructive and nondestructive characterization methods that are relevant in solid-phase processing (SPP) applications, but not necessarily applied in the MCPC Project as a guide to the planning of characterization activities in future research. Particular attention is given to measured characteristics that can correlate to other material characteristics, with a particular interest in nondestructive evaluation (NDE) that can be applied to samples obtained in the MCPC Project. Destructive examinations include tensile tests, optical and electron microscopy, micro-hardness, and residual stress tests. NDE tests include surface visual inspection, eddy current examination for cracks, 4-point potential drop, ultrasound, x-ray, and computed tomography.

36 MATERIALS SCIENCE↗
DOE OSTI2025

PET-FBA: A lightweight enzyme allocation and thermodynamics-constrained flux analysis approach to explore Escherichia coli metabolic adaptation to intracellular acidification

Escherichia coli employs diverse strategies to adapt to acidic environments that disrupt enzyme activity and the thermodynamic feasibility of essential reactions. To understand the impact of pH stress on cell metabolism, we present the PET-FBA (pH-, Enzyme protein allocation-, and Thermodynamics-constrained Flux Balance Analysis) framework. PET-FBA extends genome-scale modeling by integrating enzyme protein costs and reaction Gibbs free energy changes. Additionally, by incorporating pH-dependent enzyme kinetics in response to intracellular acidification, this framework enables the simulation of E. coli's metabolic adjustments across varying external pH levels. The model's accuracy is validated by comparing in silico growth simulations with experimental measurements under both anaerobic and aerobic conditions, as well as in silico gene knockouts of essential genes. By explicitly incorporating pH effects, our model accurately replicates the metabolic shift towards lactate production as the primary fermentation product at low pH in anaerobic conditions. This shift is only predicted when enzyme kinetics are dynamically adjusted as a function of pH. Further analysis revealed that this shift can be attributed to the reduced protein efficiency of the acetyl-CoA branch compared to lactate dehydrogenase under acidic stress, which then becomes crucial for maintaining NAD regeneration and cell growth at low pH. Furthermore, we identified strategies for enhancing cell growth under acidic anaerobic conditions by improving the enzyme activity of lactate dehydrogenase and pyruvate formate lyase, which increases NAD production efficiency and reduces enzyme protein allocation costs. Designed as a lightweight yet versatile framework, PET-FBA enables efficient genome-scale metabolic analysis. Using E. coli as a model system, our framework provides a systematic approach to understanding metabolic responses to environmental stress, pinpointing key metabolic bottlenecks, and identifying potential targets for strain optimization.

42 ENGINEERING↗
DOE OSTI2024

Concurrent measurement of strain and chemical reaction rates in a calcite grain pack undergoing pressure solution: Evidence for surface-reaction controlled dissolution

Pressure solution is inferred to be a significant contributor to sediment compaction and lithification, especially in carbonate sediments. For a sediment deforming primarily by pressure solution, the compaction rate should be directly related to the rate of calcite dissolution, transport along grain contacts, and calcite reprecipitation. Previous experimental work has shown that there is evidence that deformation in wet calcite grain packs is consistent with control by pressure solution, but considerable ambiguity remains regarding the rate limiting mechanism. We present the results of laboratory compaction experiments designed to directly measure calcite dissolution and precipitation rates (recrystallization rates) concurrently with strain rate to test whether measured rates are consistent with predicted rates both in absolute magnitude and time evolution. Recrystallization rates are measured using trace element chemistry (Sr/Ca, Mg/Ca) and isotopes (87Sr/86Sr) of fluids flowing slowly through a compacting grain pack as it is being triaxially compressed. Imaging techniques are used to characterize the grain contacts and strain effects in the post-experiment grain pack. Our data show that calcite recrystallization rates calculated from all three geochemical parameters are in approximate agreement and that the rates closely track strain rate. The geochemically inferred rates are close to predicted rates in absolute magnitude. Uncertainty in grain contact dimensions makes distinguishing between surface reaction control and diffusion control difficult. Measured reaction rates decrease faster than predicted from standard pressure solution creep flow laws. This inconsistency may indicate that calcite dissolution rates at grain contacts are more complex, and more time-dependent, than suggested by geometric models designed to predict grain contact stresses.

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

Modeling the contributions to acoustic nonlinearity from complex dislocation networks using 3D dislocation dynamics

Nonlinear ultrasonic parameters are highly sensitive to microstructural features that affect macroscale material behavior, providing a nondestructive means to characterize their evolution. Although dislocations are known to be a strong source of acoustic nonlinearity, establishing quantitative links between the acoustic nonlinearity parameter (β), measured via Second Harmonic Generation, and dislocation morphology—such as dislocation length and density—remains an open challenge. This work advances the numerical modeling of dislocation–β relationships using 3D dislocation dynamics (DD) simulations in two approaches: a “static” method computing strain and stress fields from dislocation configurations in the absence of external loading, and a “quasi-static” method to estimate β from the curvature of dislocation lines under applied load. First, the static method is combined with finite element analysis to investigate a recent assertion that heterogeneous initial strain fields can induce higher harmonic generation in a linear elastic medium; the present results do not corroborate this outcome. Then, the quasi-static method is applied to multiple-dislocation scenarios through parametric studies, revealing behaviors not predicted by analytical models, such as the competing interactions of edge and screw dislocations and the significant influence of applied stress on β. Finally, the simulations are used to model SHG experimental results and validate the hypothesis that β can decrease during plastic deformation, despite increasing dislocation density. As the DD code used here is open-source, it provides a practical platform for future investigation into microstructure–β relationships important to the interpretation of SHG results.

Materials science↗
DOE OSTI2024

Understanding the L-H isotope effect at the DIII-D tokamak and advancements in synthetic turbulence diagnostics

Abstract It is determined that while heat flux differences between hydrogen and deuterium isotope experiments result from natural differences in carbon impurity content at DIII-D, it is not the origin of the low to high confinement mode (L-H) transition isotope effect. More specifically, a two times larger edge radial electric field in hydrogen compared to deuterium is uncovered and believed to play an important role. The origin of this radial electric field difference is determined to have two possible origins: differences in poloidal rotation and turbulent Reynolds stress in the closed field line region, and increased outer strike point temperatures and space potentials on open field lines. Experimental observations from both profile and turbulence diagnostics are supported by nonlinear gyrokinetic simulations using the code CGYRO. Simulations illustrated heat transport isotope effects in the plasma edge and shear layer resulting from differences in impurity content, electron non-adiabaticity, and main ion mass dependent E × B shear stabilization. Turbulence prediction comparisons from flux-matched CGYRO simulations to experimental measurements including electron temperature, density and velocity fluctuations are found to be in good agreement with available data. A dedicated DIII-D experiment in hydrogen was performed to seed more carbon than naturally occurring, to match deuterium experiments, and possibly reduce the L-H power threshold based on gyro-kinetic predictions. To our surprise, while ion temperature gradient (ITG) turbulence was stabilized, nodiscernible change in L-H power threshold were observed in these special hydrogen experiments. In particular, it is noticed that the edge radial electric field and Reynolds stress were observed as nearly unchanging in the presence of ITG stabilization. These experimental data have enabled a more comprehensive picture of the multitude of isotope effects at play in fusion experiments, and the important potential connection between the confined and unconfined plasma regions in regulating L-H transition dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Refractive index of lithium fluoride at high dynamic stresses

Alkali halides are prototypical ionic solids whose refractive index is a fundamental property related to their lattice structure and ionic polarizability. In particular, lithium fluoride (LiF) has the largest band gap of any known transparent material and maintains transparency at high pressures, making it well suited for use as an optical window in dynamic compression experiments. While empirical fits to the measured density dependence of the refractive index have been provided, a model that is valid over the entire pressure range of experiments and based on a polarization-based theoretical description is lacking. We present a refractive index model for dynamically compressed LiF based on the Lorentz-Lorenz equation, where the molecular polarizability is determined using a single oscillator model with a strain polarizability parameter Λ=0.73. Here, we show that our modeling approach provides an excellent match to the LiF refractive index data for both shock and ramp compression experiments to 900 GPa (density compression greater than threefold). Additionally, updated fits are provided to determine the refractive index correction for LiF windows used in dynamic compression experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Delamination-informed lifecycle decisions: A dielectric and machine learning framework for composite sorting and recycling

Composite materials are widely used in aerospace, marine, and automotive sectors due to their high strength-to-weight ratio and durability. However, their long-term reliability can be compromised by damage accumulation. Specifically, delamination initiation serves as a precursor to structural failure, which is often difficult to detect during damage inspection. Identifying and sorting delamination initiation in samples not only increases operational safety while providing critical information for end-of-life decisions, which influences both the service life extension value and the efficiency of fiber extraction during recycling. This research addresses two challenges: (1) developing a nondestructive, ex-situ framework to sort composite materials based on damage severity, particularly delamination, and (2) understanding how damage in composites influences resin removal during pyrolysis. Both experimental work and finite element analysis were performed to predict critical stress levels that are associated with delamination onset. Based on these results, three loading levels 50 %, 75 %, and 90 % of maximum stress, were selected for controlled experiments, generating composite samples with varying extents of damage for machine learning model training. Microscopic imaging of these samples confirmed the damage progression from matrix cracking to delamination, validating the computational predictions. We explored supervised machine learning using dielectric measurements to classify damage states. Preliminary results show an artificial neural network can identify early delamination which is a potential precursor to failure, with 94.44 % accuracy on our dataset. A parallel investigation into the effect of damage severity on pyrolysis recycling showed that heavily delaminated samples required significantly less energy for comparable matrix removal than undamaged samples.

dielectric variables↗
DOE OSTI2025

Effect of 3D-printed surface textures on wear mechanism in 3-body abrasion of soil

This study systematically investigates the enhancement of wear resistance in 3D printed surface textures through both experimental and theoretical approaches. Three distinct surface morphologies (Smooth Surface, Surface with uniformly distributed Pits, and Surface with uniformly distributed Bumps) were fabricated using High-Impact Polystyrene, where the meso-scale textures were precisely controlled through the 3D printing process. Wear behavior was evaluated using a 3-body wear tester in an abrasive particle environment, analyzing the influence of surface textures under various operating conditions. Systematic wear tests revealed that optimally designed surface textures achieved a remarkable 77 % reduction in wear compared to the worst-performing sample. The wear mechanisms were comprehensively characterized through weight loss measurements, Scanning Electron Microscopy (SEM), and Energy Dispersive Spectroscopy (EDS) analyses, elucidating the surface morphology changes and their interaction with wear particles. Notably, the study identified how the geometric characteristics of surface textures influence the movement of wear particles and the distribution of contact stresses. Discrete element method simulations corroborated the experimental findings, providing theoretical validation for the enhanced wear resistance of the optimal structure. The high correlation between simulated wear patterns and experimental results validates the reliability of the proposed design methodology. In conclusion, these results demonstrate that 3D printed surface texturing offers a cost-effective and scalable approach to significantly improve wear resistance in engineering applications, presenting a practical alternative to conventional, high-cost surface engineering methods.

3-body abrasion↗
DOE OSTI2026

Degradation Dynamics of Perovskite Solar Cells Under Fixed Reverse-Current Injection

Previous studies of reverse-bias stability in perovskite solar cells have focused primarily on voltage-controlled tests. Here we instead investigate perovskite solar cell degradation under well-defined, constant reverse-current stress. We show that the choice of hole-transport layer dictates the dominant degradation pathway: cells using thick poly(triphenylamine) (PTAA) layers can tolerate high reverse bias but quickly undergo catastrophic breakdown under fixed reverse current near their one-sun maximum power-point. In contrast, cells modified with phosphonic-acid interface layer MeO-2PACz, with poorer indium-doped tin oxide (ITO) coverage compared to PTAA, exhibit gradual and recoverable degradation under certain stress and recovery conditions. For MeO-2PACz devices, degradation increases with both current magnitude and duration. Importantly, when normalized by injected charge (current times duration), lower currents applied over longer times cause more degradation than higher currents over shorter periods. Here, combining electrical measurements with spatially resolved photoluminescence imaging, our observations support an ion- and charge-mediated interfacial electrochemical degradation mode.

Degradation↗
DOE OSTI2025

The role of specimen size and grain boundary characteristics in the yield strength of tungsten in microtensile tests

To effectively use the measured properties from small-scale tensile tests for bulk material performance predictions, it is essential to understand the threshold of specimen size-effect strengthening and the interaction between dislocations and microstructures within miniaturized specimens. This study uses pure tungsten to investigate the size effect in terms of specimen size, grain size, and grain boundary characteristics relative to the yield strength of tungsten at room temperature. We evaluate the transition from miniaturized specimen properties to bulk properties and the deformation behavior through small-scale tensile tests of three specimen sizes (large: 80 × 100 × 233 µm³; medium: 7 × 7 × 18 µm³; and small: 2 × 2 × 5 µm³). The testing results reveal that the small and medium specimens exhibit high yield strength with ductile behavior, while the large specimens exhibit brittle failure, consistent with the room temperature strength of tungsten, indicating bulk behavior. We further explore the specimen size-effect sensitivity to yield stress and the scaling relationship between yield strength and the number of grains involved in the deformation. A power-law relationship with the exponent value of approximately -0.5 was found in the yield strength–grain number scaling, implying the Hall-Petch like behavior. A minimum of 7–17 effective grain boundaries across the tensile gauge dimension is required to accurately measure bulk properties.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

Geophysical Signatures of Crack Network Coalescence in Rocks at Multiple Length Scales

The main goal of the research project was to identify the geophysical signatures of fracture growth in natural rocks by utilizing novel geophysical techniques. The research objectives were to (a) investigate the potential for geophysical methods to determine when cracks initiate, the types and locations of propagated cracks, and the coalescence of networks of cracks in natural rocks at multiple scales, (b) determine how damage at the microscale evolved into damage at the macroscale and then link the microscopic and macroscopic observations, (c) quantify crack coalescence in rocks under realistic stress conditions using coupled mechanical-geophysical-optical visualization, and (d) identify the precursors in geophysical signals to crack coalescence. The following research thrusts were explored to achieve the research objectives: (1) uniaxial compression testing of rock specimens with and without a set of pre‐existing flaws and (2) triaxial compression testing of natural rock specimens. These thrusts allowed for exploring fracturing in rocks under realistic in situ environments and at multiple scales. This project provided educational opportunities for nine graduate and undergraduate students and resulted in 27 peer-reviewed publications. This first research thrust focused on investigating the micromechanics of fractures in rocks through uniaxial compression testing combined with advanced geophysical and imaging techniques, specifically acoustic emission (AE) monitoring, ultrasonic imaging, and 2-dimensional Digital Image Correlation (2D-DIC). By examining damage processes under time-independent and time-dependent loading conditions, insights into damage localization, crack initiation, and fracturing mechanisms were gained. It was observed that the AE signals and the strain-based measurements directly reflect the state of damage in the rock specimen and could be used to identify the cracking levels, such as the crack initiation (CI) and crack damage (CD), and the mode of deformation. A novel calibration apparatus was developed to enhance the accuracy of AE sensors, allowing for the estimation of key parameters such as magnitude, source dimension, stress drop, and radiated seismic energy associated with the fractures. The findings highlighted significant variations in the temporal evolution of AE source parameters during the primary, secondary, and tertiary stages of creep, identifying tensile cracking as the primary deformation mode. The second research thrust focused on enhancing the understanding of fracturing processes in natural rocks through triaxial compression testing, real-time AE monitoring, and ultrasonic monitoring. We investigated the impact of various factors such as fracture propagation regimes, injection parameters, rock types, and pre-existing conditions on the hydraulic fracture (HF) behavior using scaled true-triaxially loaded specimens of Barre granite and Lyons sandstone. Custom sensor housing facilitated concurrent active and passive monitoring to analyze hydro-mechanical responses and microseismicity associated with different HF scenarios. A coupled investigation of passive microseismicity and active signal attributes permitted a detailed comprehension of the various HF processes (aseismic deformation, fracture initiation and propagation, fluid permeation, and leak-off) and their dependence on the specific rock type. The findings of this research demonstrated the effectiveness of AE monitoring techniques in providing valuable insights into the impact of various factors on the behavior and dynamics of HF processes. The advancements in monitoring techniques, offering a more thorough and precise approach, represent a significant step towards optimizing HF practices and ensuring sustainable resource extraction.

58 GEOSCIENCES↗
DOE OSTI2024

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Atomic force microscopy monitoring of subcritical slow crack growth in soda lime silicate glass

Prediction of brittle fracture in silicate glasses requires a fundamental understanding of slow crack growth rates, far below the critical stress intensity (K IC ) of the material. Here, in this study, atomic force microscopy (AFM) was used to characterize slow crack propagation in soda lime silicate glasses. Using large-sample AFM systems with specialized tips and scanners, fractures less than 20 nm wide were measured as a function of relative humidity (RH). The results were found to be sensitive to the type of AFM tip used and to the dryness of the surface. Continuous scanning and image capture over several days allowed the determination of the crack tip velocity. The dependance of crack tip velocity on RH was measured by enclosing the sample and AFM head in a humidity-controlled chamber. The crack tip velocity was measured as low 4.0×10 -12 m/s at 15 % RH and as high as 3.7×10 -10 m/s at 40 % RH. Overall crack growth was smooth during the AFM scan. Additionally, recession of the crack tip was noted as the sample dried, indicating crack tip closure may be occurring. Further studies on alternative compositions will provide details on subcritical fracture processes in a range of different environmental conditions.

Atomic force microscopy↗
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)↗