Search NASA⌕ Search

SEARCH · Search NASA

Results for “TENSILE STRESS”

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 73 records · Page 4

Influence of Shear Strength Assumptions on BISON Debonding Simulations

Accurately predicting the thermomechanical response of buffer–IPyC debonding in TRISO fuel particles requires reliable mechanical property inputs for each coating layer, particularly the normal and shear strengths that influence interlayer delamination and stress concentrations. Micro tensile testing of AGR-2 fuel particles provided experimentally measured normal strengths for the buffer, IPyC, and buffer–IPyC interface; however, shear strength was not measured. As a result, BISON simulations of interface debonding must rely on assumed shear strength values, typically estimated as 20–40% of the measured ultimate tensile strength. This study evaluates how these assumed shear strength values influence cohesive zone model (CZM) predictions of buffer–IPyC separation in AGR 2 TRISO particles. Using micro tensile data from three AGR 2 compacts (2 1 3, 5 1 3, and 6 3 3), BISON simulations were performed with multiple shear strength assumptions to quantify their effect on radial and tangential stress evolution, debonding, and gap propagation. The results show that shear strength is a high sensitivity parameter: increasing the assumed shear strength significantly alters the stress distribution at the buffer–IPyC junction, shifts the predicted debonding location, and changes the extent of partial gap formation. While normal strength controls the initiation of interface separation, shear strength strongly influences the mode mixity of the failure process and the resulting stress concentrations transmitted to the IPyC and SiC layers. These findings highlight a critical gap in current TRISO mechanical characterization. Without experimentally measured shear strength, BISON simulations must rely on approximations that introduce uncertainty into predictions of coating layer integrity and fission product barrier performance. Future fuel qualification campaigns should therefore consider measurement of shear strength at the interlayer interfaces to reduce model uncertainty and improve the fidelity of TRISO fuel performance simulations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Stress and Strain Heterogeneity and Persistence in Uniaxially‐ and Triaxially‐Loaded Sandstone

Two critical questions in brittle rock mechanics are how rocks developed localized strains and to what extent internal stress heterogeneity controls this localization and subsequent macroscopic failure. Definitive answers have not yet been found, but would provide insight into rock fracture mechanics as relevant to hydrocarbon extraction and sequestration. Here, we use synchrotron X‐ray tomography (XRT) and 3D X‐ray diffraction (3DXRD) during uniaxial and triaxial tests on Nugget and Bentheimer sandstones to examine strain and stress localization prior to mechanical failure. 3DXRD was used to measure intra‐granular lattice strains which were used to compute elastic stress tensors of each grain. Digital volume correlation (DVC) was applied to XRT images to determine the strain field in the sample. Both samples featured marked spatial heterogeneity, localization, and temporal persistence of elevated stresses and strains during their mechanical deformation toward failure. Both samples featured a majority of grains with at least one principal stress component that was tensile, a signature of the influence of heterogeneity on stress transmission. Measurements further revealed that compressive stress orientations and statistics evolved in a similar manner to those of inter‐particle forces in loose granular materials, with triaxially‐compressed rock exhibiting enhanced grain stress heterogeneity compared to uniaxially‐compressed rock. Our results complement recent work by others who employed XRT and scanning 3DXRD to study triaxially‐compressed sandstone, but extend those results to uniaxial compression, sandstones of varied porosity, and grain stress measurements throughout the 3D full extent of the samples rather than in a single layer examined with scanning 3DXRD.

58 GEOSCIENCES↗

Imaging Bragg Edge Analysis TooLs for Engineering Structures (iBeatles)

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) provides pulsed neutrons with energies varying from epithermal to cold. In preparation for VENUS, the neutron imaging beamline to be located at beam port 10, we have performed a series of experiments focused on wavelength-dependent radiography and computed tomography for a broad range of applications, from materials science to biological tissues.One of the time-of-flight (TOF) techniques that is of interest to the scientific community is the 2-dimensional mapping of phases and average crystalline plane orientation in samples both ex-situ and during applied stresses such as tensile loading and heating. This technique is known as Bragg edgeimaging and relies on the identification of changes of transmission values, fitting of the edge to measure its displacement, and thus identify the shift in lattice parameter due to stresses. One of the challenges of TOF imaging measurements is the amount of data and the inability to observe Bragg edge shifts in real time during an experiment. Thus, we have been focusing on creating a Python-based interface that allows fast data processing and instantaneous mapping and fitting of the Bragg edges, and their evolution through time. Python libraries and Jupyter notebooks have been implemented to facilitate decision making during an experiment. The advantage of the notebooks is the possibility to guide an experiment as they can quickly process and display Bragg edge data. These notebooks can be used independently, or can be combined in a Python Graphical User Interface (GUI) tool called iBeatles. This interface permits visualization and fitting of the Bragg edges, and ultimately back-projects the fitting results onto the radiographs to display a strain map. Assuming data collection has sufficient statistics, the strain mapping analysis can be performed on a pixel-by-pixel basis. This development is a step forward toward a better user experience at the future VENUS beamline in terms of live feedback and productivity. Analysis that used to take days of switching between different applications can now be done in minutes within the

Bilheux, JeanChristophe [Oak Ridge National Labora↗

Modeling low cycle fatigue (LCF) of additively manufactured Hastelloy X using An accelerated crystal plasticity fatigue damage model

This paper presents a microstructure-based model for low cycle fatigue (LCF) behavior and life of Nickel-based alloy Hastelloy X manufactured using laser-powder bed fusion (L-PBF) additive manufacturing (AM). AM Hastelloy X, a solution-strengthened alloy, is tested at elevated temperature under fully reversed LCF conditions at different strain levels. A generalized plane strain finite element model is generated from electron backscatter diffraction (EBSD) characterization. The constitutive behavior of the material under fatigue is modeled using crystal plasticity and calibrated with both monotonic tensile and cyclic stress–strain data. The fatigue micro-crack initiation and propagation in the microstructure is modeled using a modified Chaboche fatigue damage model. An embedded boundary condition with a homogenous medium is used to apply the cyclic deformation and prevent numerically introduced over-constraints during fatigue simulation. A ‘cycle-jump’ method is used to accelerate the fatigue simulation and reduce the computational cost. The simulation results are compared to LCF experiments, showing satisfactory matches in cyclic stress behavior and number of cycles to macro-crack initiation for all applied strain ranges. In addition, the model illustrates the potential for quantifying microscale fatigue life impacting factors such as microstructure and surface roughness, which is needed to accurately quantify the reliability of AM components in service.

36 MATERIALS SCIENCE↗

Impact of post-irradiation annealing on mechanical performance of irradiated 718 alloy

Here, the effect of post-irradiation annealing on solution-annealed 718 alloy was investigated using advanced mechanical testing, fractography, scanning electron microscopy, and transmission electron microscopy. Specimens were extracted from a proton beam window operated at the Spallation Neutron Source, irradiated with 940 MeV protons to a maximum dose of approximately 9.7 displacements per atom (dpa) at a calculated temperature not exceeding 110 °C while in service. Helium and hydrogen concentrations reached about 1700 and 6900 atomic parts per million (appm), respectively. Despite irradiation and high tensile strength (yield stress over 1 GPa), the material exhibited significant ductility. Annealing at 500 °C, 700 °C, and 900 °C for 30 min resulted in an appreciable decrease in yield strength and an increase in ductility for annealing treatments at 500 °C and 900 °C relative to the strength and ductility of the as-irradiated material. The presence of helium and hydrogen led to cavity formation and cleavage-like brittle features on fractured surfaces; however, high-magnification imaging revealed the presence of small-scale ductile dimples, indicating that the fracture mechanism remained mixed. The annealed specimens retained total elongation levels of 14–33 %, and there was no sudden drop in ductility after heat treatments. The ductility level in the irradiated and annealed material is notable despite the presence of helium and hydrogen.

36 MATERIALS SCIENCE↗

Laser damage of crazed electron-beam high-reflectors following infrared and ultraviolet irradiation in the nanosecond pulse regime

Laser damage of optical components can be a limiting factor in scaling the energetics of high-peak and average power laser systems. Specifically for optical coatings, damage under nanosecond pulsed irradiation is initiated by pre-existing defects in the coating layers, including those that cause discontinuities in the structure, like craze lines. Crazing or cracking in a multilayer dielectric optical coating is induced when the overall coating stress is sufficiently tensile, and is an occasionally observed issue when employing more porous deposition techniques like electron-beam evaporation. Here, in this study, electron-beam high-reflectors were fabricated utilizing process parameters that are known to induce crazing based on prior processing history to systematically evaluate the impact of crazing on reflector damage performance for 1064 and 355 nm lasers. The crazing that was observed was apparently nucleated at nodular defects. When the cross-section of these nodules was investigated, it was observed that there were cracks into the fused silica substrate of approximately 5 µm in depth. The craze lines were irradiated with 1064 and 355 nm light at fluences slightly above the onset of damage initiation fluence of the coating. The 1064 nm irradiated sub-apertures exhibit laser damage but with no spatial correlation with the craze line, whereas the 355 nm irradiated area exhibited many damage sites along the craze line. Finite-difference time-domain electric-field simulations were conducted, and ∼2× field amplification in hafnia was observed for the 355 nm wavelength case. The laser damage can be attributed to a slight electric-field intensification coincidental with an area where UV damage-prone precursors are known to occur. The 355 nm laser damage in uncoated fused silica substrates has been previously correlated to initiate through localized UV absorption at the broken silica bonds in the tips of fractures.

Harthcock, Colin [Lawrence Livermore National Labo↗

Stress due to electric charge density distribution in a dielectric slab

The spatial distribution of electric field due to an imposed electric charge density profile in an infinite slab of dielectric material is derived analytically by integrating Gauss’s law. Various charge density distributions are considered, including exponential and power-law forms. Here, the Maxwell stress tensor is used to compute a notional static stress in the material due to the charge density and its electric field. Characteristics of the electric field and stress distributions are computed for example cases in polyethylene, showing that field magnitudes exceeding the dielectric strength would be required in order to achieve a stress exceeding the ultimate tensile strength.

42 ENGINEERING↗

Energy absorption of architectured PrintCast interpenetrating composites in tension

Additively manufactured (AM) metal-metal composites consisting of PrintCasted 316 L austenitic stainless-steel lattice structures infiltrated with A356 casting alloy, have recently been developed for use in high energy absorption systems with potential applications ranging from static load bearing to dynamic blast containment structures. This system has a unique mechanical behavior as the volume fraction of lattice increases showing a transition from localized to de-localized failure and dramatic increase in energy absorption capability. In this work, PrintCast A356/316 L composite tensile specimens were produced with lattice volume fractions ranging from 20 % to 50 % to capture the range of this behavior. Finite element simulations support neutron diffraction measurements of stress state. Results illustrate that in tension, the reinforcement material is in tension while the matrix support material is in compression, information offering significant insight into the transition to de-localized failure. Moreover, the simulation results provide further insight into how interfacial bonding (or lack of bonding) affects the energy absorption capabilities of the PrintCast composites.

36 MATERIALS SCIENCE↗

Structure-performance relationships in lignin-based transesterification vitrimers: The role of lignin structural features

Lignin has been hailed as an ideal renewable alternative for petrochemical-based prepolymers in material synthesis for a sustainable and circular economy, due to its abundant aromatic network and high carbon content. However, the properties and performance of lignin-derived macromolecules are strongly influenced by the lignin itself. While numerous studies have explored the impact of lignin content on the thermomechanical performance of lignin-based vitrimers, literature on how the inherent structural features of lignin affect these properties is scanty. In this study, hardwood organosolv lignin was fractionated in ethyl acetate, ethanol, and acetone to obtain lignin fractions with varying structural characteristics. These fractions were then modified through carboxylation and crosslinked with epoxidized soybean oil (ESO) at a hydroxyl to epoxy group ratio of 1:1 to produce lignin-based transesterification vitrimers (LVs). The thermal properties (i.e. glass transition temperature and thermal stability), tensile strength, storage modulus, and stress relaxation behavior of the LVs were studied and carefully related to the structural features of lignin. The results revealed a positive relationship between strong hydroxyl content in modified lignin and the tensile strength (5.10–9.71 MPa), storage modulus (1099.4 – 1372.8 MPa), crosslinking density, and stress relaxation of the LVs. Additionally, both the storage modulus and tensile strength exhibited a positive relationship with the ratio of rigid linkages in modified lignin, while lignin molecular weight was found to significantly impact the thermal properties of LVs (i.e Tg and thermal stability). This study not only highlights the valorization of lignin in vitrimer synthesis but also provide insights for designing lignin-based materials with tailored properties for specific applications.

Bio-based polymer↗

Femto-second Laser’s Enabling New Length Scale Fabrications for Rapid Post Irradiation Examination of Materials: Concluding LDRD Project Poster

Mechanical testing campaigns are required to qualify materials for advanced reactor conditions, yet economical and safety limitations restrict the number of standardized mechanical tests that can be performed. Reducing the size of the sample is one approach to addressing these challenges and to accelerating testing. Previous research has shown that smaller mechanical test samples produce higher yield and ultimate stress values compared to values measured from standard sample sizes: the “smaller is stronger” effect. Specimens used in accelerated material testing campaigns must reflect bulk material performance to enable engineering scale material property measurement. The objective of this research project was to determine if engineering scale mechanical behavior—the yield stress—could be measured with micro-tensile test samples smaller than traditional standard testing geometries. The relationship between yield stress and sample size was explored with two different nuclear-relevant structural materials: Zircaloy-4 and tungsten. Mechanical testing of both metals demonstrated decreasing yield stress values with increasing sample gauge size across three different sizes. Yield stress values from the largest gauge size, fabricated with a femto-second laser ablation system, approach bulk material yield stress values reported in published literature. Preliminary analysis of the tungsten samples indicates the yield stress value depends on the grain characteristics within the gauge section, in addition to the gauge size. Accompanying modeling efforts, including response surface generation and crystal plasticity approaches, further demonstrated that the size of the sample gauge section alone cannot explain the change in yield stress values.

36 MATERIALS SCIENCE↗

Evaluation of the Simplified Assessment Peak Equivalent Stress Design Limit Probability of Failure

The ASME Boiler and Pressure Vessel Code Section III Division 5 Subsection HH Subpart A (HHA) outlines requirements and qualification methods for graphite core components in high temperature reactor environments. One such structural qualification method is the Simplified Assessment outlined in HHA-3220, which represents a comparison of service level specific equivalent stresses to allowable stress values derived from target probability of failures (POFs). The target POF is based on the structural reliability class (SRC) for the components and service condition of interest, and is used to establish design limits for combined membrane and peak equivalent stresses. Combined membrane stress limits are derived from two parameter Weibull distribution fits of the graphite grade specific tensile strength data and the target POF. The combined membrane stress limits are then scaled based on the ratio of flexural strength to tensile strength (referred hereinafter as Rtf) to establish the corresponding peak equivalent stress limits. It is desired that the target POF is at least maintained for both combined membrane and peak equivalent design limits, especially with the use of Rtf. The assessment presented herein examines tensile and flexural data for medium and fine grade graphites relative to the calculated combined membrane and peak equivalent stress limits per HHA-3220 and illustrates the simplified assessment method results in anticipated, or at least conservative, peak equivalent design stress limits.

36 - MATERIALS SCIENCE↗

Role of mechanical stress localizations on the radiation hardness of AlGaN/GaN high electron mobility transistors

Multi-material, multi-layered systems such as AlGaN/GaN high electron mobility transistors (HEMTs) contain residual mechanical stresses that arise from sharp contrasts in device geometry and materials parameters. These stresses, which can be either tensile or compressive, are difficult to detect and eliminate because of their highly localized nature. We propose that their high-stored internal energy makes potential sites for defect nucleation sites under radiation, particularly if their locations coincide with the electrically sensitive regions of a transistor. In this study, we validate this hypothesis with molecular dynamic simulation and experiments exposing both pristine and annealed HEMTS to 2.8 MeV Au +3 irradiation. Our unique annealing process uses mechanical momentum of electrons, also known as the electron wind force (EWF) to mitigate the residual stress at room temperature. High-resolution transmission electron microscopy and cathodoluminescence spectra reveal the reduction of point defects and dislocations near the two-dimensional electron gas region of EWF-treated devices compared to pristine devices. The EWF-treated HEMTs showed relatively higher resilience with approximately 10% less degradation of drain saturation current and ON-resistance and 5% less degradation of peak transconductance. Both mobility and carrier concentration of the EWF-treated devices were less impacted compared to the pristine devices. Our results suggest that the lower density of nanoscale stress localization contributed to the improved radiation tolerance of the EWF-treated devices. Intriguingly, the EWF is found to modulate the defect distribution by moving the defects to electrically less sensitive regions in the form of dislocation networks, which act as sinks for the radiation induced defects and this assisted faster dynamic annealing.

AlGaN/GaN HEMTs↗

Wrought FeCrAl alloy (C26M) cladding behavior and burst under simulated loss-of-coolant accident conditions

Cladding burst experiments for FeCrAl cladding were performed in the Severe Accident Test Station facility at Oak Ridge National Laboratory. These experiments were simulated using the BISON fuel performance code to better understand the cladding plastic behavior and failure under simulated loss-of-coolant accident conditions. 3D cladding surface boundary conditions were generated using composite axial and azimuthal profiles from experiment thermocouple data. To improve the simulation analysis capabilities in BISON for cladding burst behavior, new thermal creep, plasticity, and failure stress models specific to C26M, a wrought FeCrAl alloy, were developed and implemented. Initial cladding burst results indicated a general underprediction in the failure temperature of the six cladding burst simulations versus the observed failure temperatures. Close investigation of the experiment timing versus the underlying tensile test data revealed that, compared with the tensile specimens, the cladding tubes did not experience the same long holding time at high temperatures. New tensile tests were performed at high temperatures using a temperature ramp similar to the simulated loss-of-coolant accident experiments. These new tensile curves showed an approximately 80% increase in the ultimate tensile strength of the C26M alloy, indicating that a holding time of 10 min at 700 °C and 800 °C allows annealing to change the material microstructure. Using the updated tensile properties, the burst temperatures and stresses from the simulations showed remarkable agreement with the experimental results. This study was then extended by varying the initial pressure to highlight the burst temperature difference between standard Zircaloy-4 and C26M cladding under equivalent conditions. The results show that C26M has a burst temperature that is approximately 70–130 K greater than that of Zircaloy-4. In conclusion, these modeling predictions can be further improved by collecting high-temperature tensile data for C26M beyond the temperature ranges used in this work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Creep and failure at metal-oxide interfaces

Small-scale bicrystal creep experiments were performed on contacts formed via in situ high-temperature diffusion bonding of metal-oxide interfaces including Ag-ZrO 2 , Pd-ZrO 2 , Pt-ZrO 2 , and Ag-high entropy oxide. This work characterizes deformation and failure at metal-oxide interfaces during mechanical loading. Interfacial sliding can be activated easily, while tensile interfacial creep was not observed at any condition of stress or temperature measured. Plastic strain, instead, localizes within the metal under tensile loading. A variety of mechanisms for plastic strain occur in the metal including lattice dislocation-mediated plasticity, twinning, low-angle grain boundary formation, and low-angle grain boundary creep. Surface and low-angle grain boundary diffusion occur under conditions where no metal-oxide tensile creep is observed, highlighting the significant differences in their interfacial mechanical response. High-temperature interfacial failure occurs when the mean curvature at the contact neck is approximately zero and the applied stresses comparable to brittle fracture stresses. The brittle fracture stresses were measured to be σ ƒ = 180 ± 90 MPa at the Ag-ZrO 2 interface at 225 °C, σ ƒ = 460 ± 160 MPa at the Pd-ZrO 2 interface at 680 °C, and σ ƒ = 640 ± 440 MPa at the Pt-ZrO 2 interface at 1010 °C.

Creep↗

Evaluating design safety margins in the American Society of Mechanical Engineers graphite core components design-by-analysis assessments

Graphite is an important material being used for core components in next-generation high-temperature gas-cooled nuclear reactors. The selection of graphite grade for a specific Designer is a complex task, dependent on reactor conditions, component functionality, and required reliability. The American Society of Mechanical Engineers (ASME) provides two semi-probabilistic design-by-analysis assessments to evaluate graphite core components against design reliability targets. The simplified assessment uses a 2-parameter Weibull distribution to describe the graphite grade’s tensile-strength distribution to establish component stress limits. The full assessment uses the 3-parameter Weibull distribution and a modified Weakest-Link Theory approach to calculate a component design probability of failure. The paper defines recommended assessment rules, which are the as-written simplified assessment and the full assessment with parameter lower bounds, the modulus update with threshold reduction, and the 2027 grouping rules. Code rules are applied to three grades: 2114, IG-110, and NBG-18. The baseline margin calculation is developed using the experimental tensile dogbone specimen. Percent margin is defined as the percent reduction in the median experimental load to obtain the allowable load per ASME assessments. Under the recommended rules, the SRC–1 margin in the simplified assessment ranged from 40.2 % to 52.7 % among the grades in this study and from 36.1 % to 49.8 % in the full assessment. The full assessment only decreases the margins by 2.5–4.5 % for the SRC-1 components and 0–1.5 % for the SRC-2 components for this baseline case. Margin is inversely related to material median strength (i.e., the strongest grade, 2114, has the lowest margin).

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Recrystallization driven softening and heating rate dependencies of FeCrAl nuclear fuel cladding during accident transients

A refined understanding of FeCrAl cladding behavior during rapid transients is critical for its potential deployment in light-water reactors. Current assessments focus on transient burst testing metrics such as balloon geometry, burst temperature, and hoop stress, often used as proxies for simpler conventional tensile properties. However, directly correlating isothermal tensile and creep data with accident transient scenarios remains a challenge, although it is essential for high fidelity model development. Recent modeling based on tensile tests up to 800 °C, conducted with both immediate loading and a 10-minute soak, showed that immediate loading better predicts experimental burst temperatures, indicating a thermal softening effect. Building upon this observation, the current study connects transient performance, microstructural evolution, and high-temperature tensile properties by leveraging results from C26M claddings burst tests performed at heating rates of 1–50 °C/s and hoop stresses from 25 to 100 MPa. At 25 MPa, rupture temperatures varied by only 6 °C, but at 100 MPa, the difference reached 116 °C, with faster heating yielding higher burst temperatures. Microstructural analysis identified recrystallization as the primary cause of heating rate-dependent softening, eliminating prior cold-working. In-situ thermomechanical data linked ballooning onset to localized instabilities, similar to ultimate tensile strength behavior in conventional tensile tests. High heating rates correlated with immediate loading tensile data, while lower rates matched soaked data. Furthermore, by linking burst performance to microstructural evolution and tensile properties, this work provides a foundation for more accurate modeling of FeCrAl claddings and potentially other Fe-based materials under accident conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Critical mineral substitutions in IN617: A combined computational and experimental approach to performance evaluation and feasibility

Addressing the escalating demand for critical minerals (CMs) driven by global climate change initiatives, this study explores compositional modifications to Inconel 617 (IN617) by substituting cobalt (Co) with manganese (Mn) across various atomic percentages. Here, we conducted a computational feasibility study employing Molecular Dynamics (MD) simulations to provide strategic guidance for experimental validation. The simulations analyzed tensile strength and corrosion resistance for five modified compositions (M1 to M5) to identify optimal properties. Tensile tests on cubic simulation cells were performed to generate stress-strain curves, revealing the impact of Co replacement with Mn on tensile strength—a metric correlated with hardness. Oxygen penetration simulations were conducted to evaluate corrosion resistance, indicating that reduced oxygen penetration depth corresponds to enhanced resistance. Promising compositions underwent phase diagram calculations for assessing phase stability. The optimal composition (M1), characterized by high tensile strength and minimal oxygen penetration, was chosen for experimental validation using induction melting and friction stir consolidation techniques. The materials were further characterized using SEM-EDS, XRD, and Vickers hardness testing. Our findings suggest that Mn substitution in IN617 can yield mechanical performance at par with high-Co alloys in energy-critical applications.

36 MATERIALS SCIENCE↗

DECOVALEX-2023: Task A Final Report

Task A, also known as HGFrac, examines the fracturing processes that may occur in the Callovo-Oxfordian claystone (COx) in the context of the high-level (HLW) and intermediate-level long-lived (ILW-LL) radioactive waste repository in France. Understanding these processes and improving numerical models to reproduce them will aid in the design, optimization, and safety of the repository. Heat and gas fracturing are studied in two independent subtasks following a stepwise approach: laboratory tests/benchmark exercises, in-situ experiment, and finally, an application case. The in-situ heater experiment aimed to thermally induce a hydraulic fracture; temperature and pore pressure were monitored to detect any evidence of fracturing. The in-situ gas injection experiment aimed to study the effect of the stress orientation and gas injection kinetics on the gas fracturing process. The occurrence of fracturing was monitored by gas pressure measured in the injection interval. In both experiments, the excavation-induced fracture network around the heater/injection boreholes played an important role in the reduction of the compressive stress state, leading to both a tensile and shear failure response of the COx. During the first two years of the project, the research teams working on each task developed and/or proposed numerical approaches for reproducing the occurrence of fracturing in the in-situ experiments. The failure criteria were defined by reproducing the measurements from laboratory extension tests for the heat fracturing subtask. In the gas fracturing subtask, their approaches were used for simulating several benchmark exercises, and an inter-comparison between models was carried out. In both tasks, the developed approaches were compared with a simplified approach considering poro-elasticity for the mechanical behaviour of the COx. Most of the developed approaches are based on a continuous medium that takes into account variations in hydraulic properties due to mechanical degradation, such as plastic deformation or damage. Other approaches implicitly modelled weak planes or embedded discontinuities to reproduce fracture propagation. The potential for fracture initiation was also studied through of a discrete approach. In the second half of the project, the research teams mainly focused on interpretative modelling of two in-situ experiments and a blind prediction exercise to test their respective approaches. The models developed by the research teams were also applied at the repository scale to evaluate fracture initiation in a case study under vi unfavourable conditions, particularly in terms of spacing between High-Level Waste cells. The results showed that the poro-elasticity approach could be an efficient tool for understanding the main processes occurring in the COx. One example is the explicit representation of the excavation-induced fracture network around the boreholes, which yielded acceptable results compared to the measurement data. However, advanced approaches were needed to evaluate the potential increase of the excavation-induced fracture network extend and better understand fracture initiation. The stress analyses carried out by the teams revealed that hydraulic boundary conditions had a strong impact on fracture initiation in the heater experiment. Furthermore, in most cases, the results required higher pore pressure increments to reach fracturing than those measured in the experiment. This implies that the measurements may have been biased by the packer’s capacity to fully isolate the piezometric chambers, leading to lower pressures. On the contrary, there was no agreement on the fracturing mode, as some reported either shear or tensile fracturing, while others reported a combination of the two modes. In the case of gas fracturing, the research teams were limited to the comparison of a single point, which complicated their task. Nonetheless, the numerical results were able to reproduce the measurements and capture processes such as longitudinal gas flow through the excavation-induced fracture network, as suggested by some evidence in the observation piezometric chambers. The numerical models also agreed with the measurements in the sense of higher probability of developing along the injection borehole than radially towards the sound rock. The approaches developed by the research teams showed that they are capable of analysing and reproducing fracture initiation in the COx. However, areas of future work should focus on the fracture propagation and fracture aperture, which were out of the scope of this task. To this end, additional data must be gathered for the parameter characterisation and validation of the numerical models. Nonetheless, various approaches showed promising results as they were able to reproduce fracture development under certain conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗