# Coupled Experimental / Numerical Approach to Determine the Creep Behavior of Zr-4 Cladding Under LOCA Condition

Damien Campello, Nicolas Tardif, Marie-Christine Baietto, Michel Coret, Jean Desquines

► To cite this version:

Damien Campello, Nicolas Tardif, Marie-Christine Baietto, Michel Coret, Jean Desquines. Coupled Experimental / Numerical Approach to Determine the Creep Behavior of Zr-4 Cladding Under LOCA Condition. iDICs 2016, Nov 2016, Philadelphia, United States. 10.1007/978-3-319-51439-0_54 . hal-01863191

HAL Id: hal-01863191

https://hal.science/hal-01863191

Submitted on 28 Aug 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Coupled Experimental / Numerical Approach to Determine the Creep Behavior of Zr-4 Cladding Under LOCA Condition.

D. Campello<fnref n="1" />, N. Tardif<fnref n="1" />, M.-C. Baietto<fnref n="1" />, M. Coret<fnref n="2" />, J. Desquines<fnref n="3" />

## ABSTRACT

The thermo-mechanical behavior of Zircaloy-4 fuel rods under Loss-Of-Coolant Accident (LOCA) conditions is investigated. A custom experimental setup is dedicated to the high-temperature creep ballooning study of 90mm long cladding samples. Creep tests were performed under an inert environment (argon), for temperatures from 750 to 850 °C and internal pressures ranging from 1 to 5MPa. As the high-temperature creep of metals is strongly influenced by the temperature, the setup allow for a heterogeneous thermal distribution along the specimen. A unique test provides a rich database about the steady-state creep of the alloy. A first campaign is dedicated to bare Stress Relieved Annealed Zr-4.

Creep-rates are computed along a generatrix of the tube using 2-Dimensional Digital Image Correlation (2D-DIC) and are correlated to thermal distribution measurements allowed by Near Infra-Red Thermography (NIRT). As-received Zr-4 behavior laws are determined using Finite Element Model Updating (FEMU). Norton exponents and activation energies are determined for thermal-mechanical conditions of the performed tests.

**Keywords:** 2D-DIC, NIRT, FEMU, creep power-law, Zircaloy-4.

## INTRODUCTION

Loss-Of-Coolant of the primary loop is a complex multi-physical design basis accident. The study of the thermal mechanical and thermal hydraulics behavior of a fuel rod assembly is the context of PERFROI project [1]. It relies on both experimental and modeling activities about the flow blockage caused by fuel rods deformation, fuel fragment relocation and the potential loss of cladding integrity. In this framework, the ballooning behavior of fuel rod claddings submitted to high temperature and internal pressure is investigated.

Creep-tests and thermal ramps are usually performed in order to investigate the cladding behavior under LOCA conditions [2-5]. The setups proposed in this literature provide on a test a single temperature condition assuming thermal field is homogeneous when testing. A unique stress condition is often investigated. The number of the tests to be performed is consequent in order to accurately determine material parameters.

The study focuses on creep ballooning tests. Creep power-law in equation (1) is assumed to well model the secondary-creep behavior $\dot{\epsilon}_{eq}$ of zirconium alloys [6] depending on the temperature $T$ and the von Mises stress $\sigma_{vM}$. The multiplicative coefficient $A$, the Norton exponent $n$ and the activation energy $Q$ are widely investigated with experimental studies.

$$ \dot{\epsilon}_{eq}(\sigma, T) = A \cdot \sigma_{vM}^n \cdot e^{\frac{-Q}{R \cdot T}} \qquad (1) $$

These creep parameters depend on the temperature and applied stress conditions.

## SETUP

The experimental setup initially designed by Tardif et al. [1] is upgraded to test Zircaloy-4 claddings under inert atmosphere using argon gas [8]. The thermo-mechanical loading applied to the specimen is represented in Fig. 1.

The tube is first heated using an induction system, and then 3 successive internal pressure levels are applied using argon gas. The bottom end effect is compensated. The bench is fitted to perform 2D-DIC [2] and NIR measurements [3].

Axial and hoop displacements fields are calculated in the ROI using a 2D-DIC algorithm. It is performed using Ufreckles software developed at LaMCoS by J. Réthoré [15]. The conservation of the optical flow is solved using a non-linear least-square method relying on a finite element basis. Performing 2D-DIC on tubular specimen is not common considering the 3D geometry. Uncertainties related to out-of-plane displacement (ballooning) have been estimated and are negligible with respect to the cladding deformation when steady-state creep rates are calculated. The equivalent steady-state creep-rates are calculated assuming material incompressibility and using the von Mises criterion.

<footnote n="1">¹ Université de Lyon, CNRS, INSA-LYON, LaMCoS UMR 5259, 20 Avenue Albert Einstein, F69621 Villeurbanne Cedex, France</footnote>

<footnote n="2">²GeM (UMR 6183), Ecole Centrale de Nantes, 1 rue de la Noë, 44321 Nantes, France</footnote>

<footnote n="3">³IRSN/PSN-REX/SEREX/LE2M, Saint Paul lez Durance, France</footnote>
<figure><img src="image_2.png" /><figcaption>Fig. 1 Schematic thermal mechanical loadings applied to the specimen</figcaption></figure>

The second digital image processing is dedicated to near infra-red measurements [9]. Standard cameras (with silicon detectors) are preferred to infra-red cameras in order to higher the image definition. The K₁ and K₂ parameters of the radiometric model detailed in equation (2) are calibrated using the data measured during the heating part of the test on the specimen itself.

$$T = \frac{K_1}{\ln(\frac{K_2}{I} + 1)} \qquad (2)$$

## RESULTS

The steady-state creep rate and temperature axial distributions measured in the gauge length (z₀=20 mm) on a specimen heated up to 850 °C are plotted in Fig. 3.

<figure><img src="image_3.png" /><figcaption>Fig. 2 Kinematic and thermal axial distributions</figcaption></figure>

The thermal distribution clearly induces a variation of the creep-rates along the height of the specimen. The heterogeneous thermo-mechanical conditions are fully capitalized with the modeling of the tests. An iterative Gauss-Newton algorithm is then performed in order to lower the error between experimental creep-rate logarithms and calculated ones. A constitutive behavior law is determined using the data of a unique test. The measurements of nine experiments were processed in order
to map the Norton exponent and the activation energy as represented in Fig. 3.

<figure><img src="image_4.png" /><figcaption>Fig. 3 Norton exponent and activation energy determined in the domain of the study using data of 9 experiments</figcaption></figure>

The results are in a good agreement between each other and are consistent with literature data. A stress of 25 MPa is determined in the α-phase domain of Zr-4 as a boundary between diffusion and dislocation creep mechanisms. The microstructure of the alloy is changing from α to α+β phase at a temperature of 820 °C. As an effect, the determined activation energy is changing from around 180 kJ/mol to 300-400 kJ/mol in this domain.

## CONCLUSION

A FEMU methodology was performed in order fully capitalize high-temperature creep ballooning tests data and determine the parameters a steady-state creep power-law. The ballooning structural effect are accounted for using FEM granting access to the true stress level and state over a 20 mm gauge length. Results of seven tests are detailed in this paper addressing to thermal conditions ranging from 750 to 850°C and von Mises stresses in a range of 7 to 48MPa.

The effect of applied stress on the Norton exponent is evidenced in the study for both α-phase and (α+β)-phase domains. The determined stress exponent magnitudes are very consistent with literature data depending on the thermo-mechanical test conditions.

The activation energy parameter is very sensitive to the thermo-mechanical conditions and remains valid within a limited temperature and stress ranges to the authors' knowledge. The beginning of the (α+β)-phase domain is clearly identified with the activation energy increase.

## ACKNOWLEDGMENTS

The study was performed in the framework of the PERFROI ANR project (n° ANR-11-RSNR-0017-01). The authors would like to acknowledge all the participants of the project: EDF, IRSN-SEREX, LEMTA, Ecole Centrale of Paris.

## REFERENCES

[1] Repetto, G., Dominguez, C., Durville, B., Carnemolla, S., Campello, D., Tardif, N., and Gradeck, M., "The R&D PERFROI project on the thermal mechanical and thermal hydraulics behaviors of a fuel rod assembly during Loss-Of-Coolant Accident," *Proceedings of the Nureth16*, Chicago, USA, August 30-September 4, (2015).

[2] Busby, C., and Marsh, K., "High temperature, time-dependent deformation in internally pressurized Zircaloy-4 tubing (LWBR Development Program)," in *Tech. report, Colorado Cooperative Wildlife Research Unit*, Fort Collins, USA, (1974).

[3] Kaddour, D., Frechinet, S., Gourgues, A., Brachet, J., Portier, L., and Pineau, A., "Experimental determination of creep properties of zirconium alloys together with phase transformation," *Scripta Materialia*, **51**, 6, (2004), 515-519.

[4] Rosinger, H. E., Bera, P. C., and Clendening, W. R., "The steady-state creep of Zircaloy-4 fuel cladding from 940 to 1873K", in *Tech. report, Atomic Energy of Canada Limited*, (1978).

[5] Réocreux, M., and de Martinville, E. S., "A study of fuel behavior in PWR design basis accident: and analysis of results from the PHEBUS and EDGAR experiments," *Nuclear Engineering and Design*, **124**, 3, (1990), 363-378.

[6] Hayes, T. A., and Kassner, M., "Creep of zirconium and zirconium alloys," *Metallurgical and Materials Transactions A*, **37A**, (August 2006), 2389-2396.

[7] Tardif, N., Coret, M., and Combescure, A., "Experimental study of the fracture kinetics of a tubular 16MnNiMo5 steel specimen under biaxial loading at 900 and 1000°C. Application to the rupture of a vessel bottom head during a core meltdown accident in a pressurized water reactor," *Nuclear Engineering and Design*, **241**, 3, (2011), 755-766.

[8] Campello, D., Tardif, N., Baietto, M.-C., Coret, M., Desquines, J., "Secondary creep behavior of Zr-4 claddings under LOCA conditions," *Proceedings of the TOPFUE 2016*, Boise, USA, September 11-16, (2015).