Experimental evaluation and modeling of fatigue of a 316L austenitic stainless steel in high-temperature water and air environments

Austenitic stainless steels is used in many components of nuclear power plants, particularly in the pipes of cooling systems. Owing to power transients and to start-ups and shutdowns, these components are subjected to thermo-mechanical loadings (low-cycle fatigue LCF & high cycle fatigue HCF) and flow-induced vibration (HCF). The corresponding fatigue design curves were established with solid smooth specimens tested in air at room temperature under strain-control. However, these curves do not consider the influence of LWR water environments that were shown to reduce fatigue life. US NRC Regulatory Guide 1.207 or other national equivalents (JSME Code in Japan) then were established taking strain rate, dissolved oxygen and temperature into consideration. Besides, a significant number of issues have been recently identified, which may have an impact on fatigue life but have not been sufficiently investigated, e.g., the potential negative effects of mean stress, mean strain, surface finish, long-term static hold, specimen geometry, multiaxial stress state were sensitivities not explicitly addressed. This study was launched to assess the effect of mean stress on fatigue behavior of a 316L austenitic steel in boiling water reactor environment with hydrogen water chemistry (BWR/HWC) at 288°C.

Load-controlled fatigue tests were selected as the easiest experimental technique to impose a pre-defined mean stress. The tests were carried out on hollow specimens at different stress amplitudes and mean stresses in BWR/HWC environment and in air at 288°C, the latter serving as reference environment to evaluate the life reduction induced by the BWR/HWC medium. Few tests in strain-controlled were performed to derive a consistent way to represent the data obtained with the two control modes together.

The test results reported in the form of stress-life curves showed that, in LCF regime (< 1e5 cycles), positive and negative mean stresses increase and BWR/HWC environment decreases fatigue life. In the HCF regime, negative mean stress is always beneficial for fatigue life but positive mean stress decreases fatigue life in BWR/HWC. The beneficial effect of mean stress is attributed to its enhanced cyclic hardening on material, which leads to smaller strain amplitude at given stress amplitude. The load-controlled data were converted into strain-life presentation by considering the average strain amplitude over the whole loading cycles. Additionally, a modified Smith-Watson-Topper mean stress correction method was successfully used to correlate the results with and without mean stress. Finally we showed that strain energy density criteria are good alternatives to correlate all data.

Crack growth rates were determined from the striation spacing on the fracture surfaces with high resolution scanning electron microscopy (HRSEM); crack initiation sites were also studied with HRSEM, and the microstructures were observed with transmission electron microscopy and electron channeling contrast imaging. From the striation spacing measurement, we concluded that BWR/HWC environment essentially reduces the number cycles needed to initiate a physical crack (crack depth < 50 microns) but modifies slightly the crack growth rate, which correlates well with a strain intensity factor and J-integral method.

Life predictions were also done by a well-trained artificial neuron network that considers mechanical, environmental and material properties factors.

Advanced analytical electron microscopy applied to Solid Oxide Cell materials and their degradation

Stéphane François Jean Poitel

Electron microscopy offers a wide range of techniques to study materials. Environmental Scanning Electron Microscopy (ESEM), in particular, allows the characterisation of samples at high-temperature using a laser heating stage, and under various gaseous atmospheres leading to direct observation of the sample surface morphology dynamics in conditions that are close to those in real operating devices (i.e. at high temperature in various gases). Solid Oxide Fuel Cells (SOFC) operating around 800°C in both oxidising and reducing gases, with proven co-generation efficiencies of 90%, offer very promising solutions in the ongoing energy transition. However, production costs and materials degradation in those demanding operating conditions remain significant barriers to their large-scale deployment. Studying SOFC materials' degradation through a combination of ESEM with other advanced microscopy techniques fulfilled the double aim of demonstrating the power of an ESEM and improving the understanding of SOFC materials' degradation during operation. Two materials from the SOFC stack were selected for degradation case studies at around 800°C in oxidising/reducing gases: coated ferritic steel interconnect and Ni-YSZ anode. The objectives were to explore the parameters' influence (temperature, heating ramp, pressure, gas exposure...), set the boundaries and define the guidelines of the ESEM technique while demonstrating its capabilities. The two different materials sets were selected for their relevance both for the SOFC field and for their different specific characteristics allowing a wide range of experiments to be explored. A specific method for the study of the oxidation of coated steel was developed combining ESEM, Focused Ion Beam (FIB) and Scanning Transmission Electron Microscopy (STEM). Post-test segmentation of the in-situ observed images demonstrated the possibility to obtain quantitative data. FIB sample preparation from the exact site of ESEM observation allowed to correlate the underlying microstructure with the in-situ surface observation. Applying this method, oxidation mechanisms of several coatings were characterised and specific behaviours observed and explained such as the appearance of oxide ridges and the behaviour of a cerium barrier layer in case of cobalt-cerium coated steel. The heating ramp speed and the oxidation temperature were demonstrated to be crucial parameters to obtain the desired microstructure. The nitriding pretreatment of interconnect was investigated too and shows promising results for improving the steel oxidation resistance. Another combination of techniques was developed for the study of the Ni-YSZ anode by complementing the ESEM observation with continuous Mass Spectrometer (MS) acquisition. The first reduction of the NiO-YsZ anode was analysed under gas flows of hydrogen and hydrocarbons. Redox cycling effect on the anode microstructure was investigated and the resulting observations were in line with models from the literature for the reduction and re-oxidation kinetics. Finally, the demonstration of an active solid oxide cell in the microscope operating at more than 700°C under mixed gas paves the way to further possible studies approaching the device's observation in real operating conditions. Limitation of the ESEM techniques are discussed and recommendations for further improvements and perspectives given.

EPFL2020

Finite element modeling and experimental study of brittle fracture in tempered martensitic steels for thermonuclear fusion applications

Pablo Federico Mueller

In this work we have studied brittle fracture in high-chromium reduced activation tempered martensitic steels foreseen as structural materials for thermonuclear fusion reactors. Developing the adequate materials that can withstand the severe irradiation conditions of the burning plasma in a future fusion reactor is one of the major challenges to be solved in order to make profit from the great advantages of thermonuclear fusion as an energy source. High-chromium tempered martensitic steels such as F82H and the most advanced version Eurofer97 are among the main candidate materials for structural applications in future fusion power plants due to low irradiation-induced swelling, good mechanical and thermal properties, and reasonably fast radioactive decay. The most concerning drawback of these kind of steels is irradiation embrittlement, which is manifested by a ductile-to-brittle transition temperature shift to higher temperatures after irradiation whose amplitude depends on the irradiation conditions (temperature, neutron flux, neutron fluence, etc). The aim of this work was to study and model brittle fracture in the ductile-brittle transition region of this kind of steels in the as-received unirradiated conditions. It is necessary to be able to transfer laboratory specimen fracture data to real components and structures in order to assess the performance of these steels in the different operating and transient conditions they could find during the operation life of a fusion reactor. In order to do so, the specimen geometry effects and specimen size effects on measured fracture toughness need to be properly understood, taken into account and predicted with an appropriate model. In particular, specimen size effect on measured toughness is a major concern for the nuclear materials research community owing to the limited irradiation volume in current and planed materials irradiation facilities. The main results of this PhD work are summarized below. The microstructure of Eurofer97 and F82H has been characterized and compared by means of optical microscopy, scanning electron microscopy, transmission electron microscopy and energy-dispersive X-ray spectroscopy in order to identify microstructural features that could play a role in the measured fracture toughness. Both steels have similar but slightly different chemical composition and final heat-treatments but the prior austenitic grain size measured in F82H is approximately 8 times larger than in Eurofer97. It was shown that the alloying element Tantalum, added to stabilize the austenite grain size, played a different role in both steels. After a careful analysis of the particles present in both steels, it was found that Tantalum in Eurofer97 formed carbides of an average size around 100 nanometers. In contrast in F82H it did not form small carbides but formed big oxide inclusions with a size up to 30 µm. These large particles do not effectively pin the grain boundaries. The different behavior of Tantalum in these steels is believed to be mainly a consequence of the larger content of Oxygen present and the smaller amount of Aluminum in F82H compared to Eurofer97. The Master-Curve ASTM-E1921 standard is a method initially developed to determine the ductile-to-brittle transition reference temperature, T0, in fission reactor ferritic steels from a small number of experiments. In this work the applicability of the Master-Curve method to reduced activation tempered martensitic steels such as Eurofer97 and F82H was studied in detail. Fracture tests with pre-cracked sub-sized compact tension specimens (three different sizes, 0.18T, 0.35T and 0.87T of Eurofer97) were carried out in the temperature range [-196 °C, -40 °C]. The toughness-temperature behavior and scatter were shown to deviate from the ASTM-E1921 standard predictions near the lower shelf. Using the method of maximum likelihood, the athermal component of the Master-Curve was calculated to better fit our fracture toughness data from the lower to the middle transition region. We showed that these Master-Curve adjustments are necessary to make the T0 values obtained near the lower shelf with 0.35T size C(T) specimens consistent with those obtained in the middle transition region with 0.87T C(T) specimens. The ASTM-E1921 specimen size limitations, setting the maximum toughness measurable with a given specimen size, were found to be too lenient for this kind of steels. This problem was especially evident in fracture toughness data of Eurofer97 and F82H obtained in the upper transition temperature range with two different specimen sizes. Thus a more stringent specimen size requirement was proposed to avoid inconsistent transition temperature determinations. A promising local fracture model with the potential of predicting cleavage fracture toughness was studied. Finite element simulations were undertaken for compact specimens, notched specimens and tensile specimens of Eurofer97 steel tested from 20 °C down to -197 °C. Three and two dimensional as well as axisymmetric simulations were run in order to calculate the stress and strain fields at the onset of brittle fracture. For each tested temperature, the calculated load-displacement curves were found to reproduce very well those of the experiments. A local approach fracture criterion was studied. This criterion states that when the maximum principal stress is larger than a critical stress within a critical volume ahead of the crack tip, notch or neck, cleavage is triggered leading to macroscopic fracture of the specimen. It was shown that this model is able to predict the minimum fracture load of the notched specimens with the same values of critical stress and critical volume that were calibrated to predict the lower bound of fracture in compact specimens. This local approach model was also successfully used to predict the strong size effect observed experimentally in pre-cracked compact tension specimens in the upper transition region. The critical fracture stress determined in the standard tensile specimens was found higher than that of the fracture specimens. It was suggested that this difference stems from a significant difference in the stress state between the different specimens (triaxiality). Finally, the comparison between the overall fracture behavior in the transition between F82H and Eurofer97 steels indicated that these two materials are quite similar. A difference in the reference temperature T0 of about 20 °C was found, with a nominal value of -100 °C and -80 °C for F82H and Eurofer97 respectively. However, some excessive scatter in the toughness data was found in F82H with more data points than expected lying below the 1% tolerance bound; this was not observed in Eurofer97. From an engineering point of view, the entire fracture databases of these two steels are encompassed practically by the same lower tolerance bound.

EPFL2009

Constitutive behavior and fracture properties of tempered martensitic steels for nuclear applications

Raul Alejandro Bonadé

In this study, the tensile and fracture properties and the microstructure of the reduced-activation tempered martensitic steel Eurofer97 have been investigated. This technical alloy is a 9%Cr steel developed within the European fusion material research program. To a lesser extend, the plastic flow properties of a equiaxed ferritic Fe9%Cr model alloy were also studied for comparison with those of the tempered martensitic structure. The main objectives of this work are described hereafter: to correlate the microstructural features with the plastic flow properties measured by the tensile tests for both the Eurofer97 steel and the model alloy. The correlation established should be reflected in a physically-based model of plastic flow. to study the fracture properties of the Eurofer97 steel in details in the lower ductile-to-brittle transition region. to calculate with finite element modeling the stress fields at the crack tip. This information is further used in conjunction with a local approach model for quasi-cleavage to reconstruct the fracture toughness-temperature curve. Tensile tests were carried out at different imposed nominal strain rate at several temperatures from 77K up to 473K, on both the Eurofer97 steel and the Fe9%Cr model alloy. The temperature dependence of the yield stress was precisely determined. As expected for body-centered cubic (BCC) materials, a strong increase of the yield stress by decreasing temperature, below 200 K, was observed. At higher temperatures, the temperature dependence of the yield stress was found much weaker, being associated with the temperature dependence of the shear modulus. Efforts were made to analyze in details the post-yield behavior (strain-hardening) as a function of temperature. The post-yield behavior was modeled using the Kocks phenomenological model based on the competition of storage and annihilation of dislocations. While this model was originally developed for face-centered cubic (FCC) metals where the rate-controlling mechanism of dislocation motion is the dislocation-dislocation interaction, we used is model in the high temperature domain (T>200K) of BCC materials, to model the strain-hardening evolution of the Eurofer97 steel and the model alloy. The values obtained for the key parameters of the model, namely the dislocation mean free path and annihilation coefficient, were found consistent with the microstructural features. The parameters temperature-dependence observed was also consistent with the physics of two basic mechanisms of dislocation storage and annihilation (dynamical recovery). For the low temperature domain, the strain-hardening model was modified to account for the strong Peierls lattice friction, based on an original idea of Rauch. Transmission electron microscopy observations were done to characterize the undeformed microstructures and their evolution with strain. A clear correlation was established between the stress-dependence of the strain-hardening and the microstructures. Fracture toughness tests on the Eurofer97 steel were performed with 0.2T C(T) and 0.4T C(T) specimens in the lower transition. Five temperatures were selected at which many tests were repeated to determine the amplitude of the inherent scatter of quasi-cleavage. The temperature were chosen in such a way that the measure fracture toughness remain below 150 MPa m1/2, which correspond for the 0.4T C(T) at a M value larger than 70 at the highest temperature. Such a M value for 0.4T C(T) specimens is known to ensure enough constraint. An attempt to analyze the experimental data in the framework of the master-curve approach by following the ASTM E-1921 standard was done. It was clearly demonstrated that the assumed shape of the toughness-temperature curve as described in the ASTM E-1921 standard for the reactor pressure vessel steels is not adequate for the tempered martensitic Eurofer97 steel, which present a particularly steep transition. Our Eurofer97 fracture database was then compared to the existing one on another similar steel, the F82H. Differences in the amplitude of the scatter of both steels was found while the lower bound of the toughness-temperature curve, describing the 1% failure probability was shown to be the same. 2D finite element simulations of the compact tension specimens were performed at various temperatures using the constitutive laws determined previously. The stress field around the crack tip were calculated and used to determine a local criterion of quasi-cleavage. The criterion is defined by the attainment of a critical stress encompassing a critical area. The lower bound of the toughness-temperature curve was then successfully reconstructed by using this local criterion. Finally, the relationship between the critical area and the applied stress intensity factor for the C(T) specimen was shown to follow a power law whose coefficients are dependent on the real dimensions of the specimen. Such a relationship allows scaling the C(T) toughness data from one size to another in case of in-plane constraint loss.

EPFL2006