Issue 
EPJ Nuclear Sci. Technol.
Volume 2, 2016



Article Number  22  
Number of page(s)  9  
DOI  https://doi.org/10.1051/epjn/2016016  
Published online  18 April 2016 
https://doi.org/10.1051/epjn/2016016
Regular Article
2D simulation of hydride blister cracking during a RIA transient with the fuel code ALCYONE
^{1}
CEA, DEN, DEC, Bâtiment 151, 13108 SaintPaullezDurance, France
^{2}
CEA, DEN, DMN, 91191 GifsurYvette, France
^{3}
EDF, SEPTEN, 69628 Villeurbanne Cedex, France
^{⁎} email: jerome.sercombe@cea.fr
Received:
16
September
2015
Received in final form:
1
March
2016
Accepted:
8
March
2016
Published online:
18
April
2016
This paper presents 2D generalized plain strain simulations of the thermomechanical response of a pellet fragment and overlying cladding during a RIA transient. A fictitious hydride blister of increasing depth (25 to 90% of the clad thickness) is introduced at the beginning of the calculation. When a predetermined hoop stress is exceeded at the clad outer surface, radial cracking of the blister is taken into account in the simulation by a modification of the mechanical boundary conditions. The hoop stress criterion is based on Finite Element simulations of laboratory hoop tensile tests performed on highly irradiated samples with a throughwall hydride blister. The response of the remaining clad ligament (beneath the cracked blister) to the pellet thermal expansion is then studied. The simulations show that plastic strains localize in a band orientated at ∼45° to the radial direction, starting from the blister crack tip and ending at the clad inner wall. This result is in good agreement with the ductile shear failures of the clad ligaments observed postRIA transients. Based on a local plastic strain failure criterion in the shear band, ALCYONE simulations are then used to define the enthalpy at failure in function of the blister depth.
© J. Sercombe et al., published by EDP Sciences, 2016
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The behavior of high burnup fuel during a Reactivity Initiated Accident (RIA) has been studied experimentally in the NSRR [1,2] and CABRI reactors [3,4]. It is now well established that the accumulation of hydrides beneath the thick outer zirconia layer that can form in Zircaloy4 claddings during base irradiation is a key factor with respect to fuel rod failure during the Pellet Cladding Mechanical Interaction (PCMI) phase of a RIA [5]. In the extreme case of outer zirconia spalling, the local cold spot that appears triggers hydrogen diffusion in the cladding resulting in a massive hydride precipitation and eventually to a hydride blister (or lens).
Many experimental works have shown that precipitated hydrides result in a loss of ductility of zirconium alloys, especially at low temperatures. In the extreme case of a throughwall hydride blister, the failure can be brittle with no residual strains [6]. During simulated RIA transients on Zircaloy claddings, it has been reported that the rod failure proceeds in a mixed mode with a brittle fracture of the heavily hydrided periphery of the cladding and a ductile propagation in the remaining clad ligament [1–4]. Ductility is here associated to the change of direction of the throughwall crack, radially orientated in the hydride rim or blister and then bifurcating at ∼45° until the clad inner wall.
In this paper, the failure of a fuel rod containing a fictitious hydride blister of varying thickness during a simulated RIA transient is studied with the 2D generalized plain strain scheme of the fuel code ALCYONE. The relationship between the blister depth and the maximum fuel enthalpy is seeked by multiple simulations of the CABRI REPNa8 test [3,4].
2 The 2D model of the fuel code ALCYONE
ALCYONE is a multidimensional fuel code codeveloped by the CEA, EDF and AREVA within the PLEIADES environment which consists of three different schemes [7]: a 1.5D scheme to model the complete fuel rod, a 3D scheme to model the behaviour of a pellet fragment with the overlying cladding, a 2D(r,θ) scheme to model the midpellet plane of a pellet fragment, see Figure 1. The different schemes use the same Finite Element (FE) code CAST3M [8] to solve the thermomechanical problem and share the same physical material models at each node or integration points of the FE mesh. A detailed description of the main models and material parameters considered in the thermomechanical code ALCYONE can be found in references [9,10].
Postirradiation examinations performed on Pressurized Water Reactors (PWR) pellets after 2 to 5 cycles of base irradiation show that the pellets are usually broken in ∼6–10 pieces of irregular size [9]. In the 2D simulation, the behavior of an average fragment representing one eighth of the pellet is studied. Because of the geometrical symmetries, only one sixteenth of the pellet and of the overlying piece of cladding is meshed. The mechanical boundary conditions considered in the 2D calculations are shown in Figure 1. The opening and closing of the radial cracks between the pellet fragments is allowed by applying a unilateral condition (u_{y} ≥ 0) on the nodes of the (0x) line. At the pelletcladding interface, unilateral contact is assessed and a Coulomb model is introduced to simulate frictionslip or adherence.
The fictitious blister crack is introduced by modifying the boundary conditions on the axis of symmetry of the pellet fragment (0r). Initially, the tangential displacement u_{t} of the nodes is set to zero. When the hoop stress on the external clad wall reaches a predefined threshold, the boundary condition on the nodes included in the depth of the hydride blister is modified. A unilateral condition is applied to avoid nonphysical interpenetration with the symmetric part of the blister (u_{t} ≤ 0). Note that this simplified approach implies that blister cracking has an infinite length in the axial direction (that of the rod axis of symmetry).
Fig. 1
Mesh and mechanical boundary conditions in the 2D scheme of ALCYONE. 
3 Material properties for the cladding and the hydride blister
To model the behavior of fresh and irradiated Zircaloy4, the constitutive law developed in reference [11] is used in ALCYONE. It consists in a unified viscoplastic formulation with no stress threshold between the elastic and viscoplastic regimes. The textureinduced plastic anisotropy of Zircaloy4 is described by a Hill's quadratic criterion. The model includes four parameters (strain rate sensitivity exponent, strength coefficient, strain hardening coefficient, Hill's coefficients) that have been adjusted on an extensive database of laboratory test results (axial tensile tests, hoop tensile tests, closedend internal pressurization tests) essentially obtained from the PROMETRA program, dedicated to the study of zirconium alloys under RIA loading conditions [12]. The model is able to account precisely for the impact of temperature, strain rate, and irradiation damage on the ultimate stress, on the strain hardening exponent (up to uniform elongation) and on the plastic anisotropy of the material.
The explicit modeling of a hydride blister is a complex problem which would require the realistic modeling of outer zirconia formation and the partial spalling of the layer, the thermodiffusion of hydrogen and the volume expansion associated with the precipitation of δhydrides [13]. Such a work is far beyond the goal of this paper. In the simulations, we assume that a stable and nonevolving hydride blister is present at the beginning of a RIA pulse test. In this respect, it is implicitly assumed that irradiation creep of Zircaloy4 during base irradiation is sufficient to relax internal stresses generated by the precipitation of δhydrides. In the simulations, the thermal (heat capacity, thermal conductivity) and mechanical (Young modulus, Poisson ratio) properties of the cladding zone where the hydride blister is located are furthermore identical to those of the remaining cladding.
The only specific parameter required in the 2D simulations is the stress to failure of the hydride blister. An approximate value of ∼145 MPa was deduced by Desquines et al. [6] from a hoop tensile test performed on an irradiated highly corroded clad sample containing a throughwall hydride blister (PROMETRA test 2468, Zircaloy4, strain rate 5/s, temperature 480 °C). The failure of the sample actually took place outside of the gage section. Interpretation of hoop tensile tests is however complex due to structural effects that occur during the experiment (bending, friction…). A detailed Finite Element analysis where the clad section and the half cylinder inserts are considered can nevertheless provide realistic estimate of the plastic strains [11,14]. The simulation of test 2468 has therefore been undertaken and shows that the stress state is far from being homogeneous in the clad thickness and width and depends greatly on the exact position of the blister (Fig. 2). With a friction coefficient of 0.1, the hoop stress on the clad outer wall at failure and out of the gage section varies between 150 and 250 MPa.
Fig. 2
Hoop stresses calculated at failure time during the PROMETRA test 2468 (left: friction coefficient 0.1, right: friction coefficient 0.4). 
4 Simulation of the CABRI REPNa8 test
The CABRI REPNa8 test was performed on a highly corroded UO_{2}/Zircaloy4 fuel rod (maximum corrosion thickness 84–126 μm) with partial spalling detected before the test. The main characteristics of the test are recalled in Table 1 (from Ref. [3]).
The REPNa8 test led to the loss of tightness of the rod at an enthalpy of 78 cal/g. Several microphone (or acoustic) signals were however recorded before the gas ejection in the coolant. At an enthalpy level of 44 cal/g, a microphone event located near the Peak Power Node (PPN) has been correlated to a limited axial crack extension inside a hydride blister (depth ∼50% of the clad wall thickness), suggesting a possible failure initiation without loss of tightness [3,4].
A preliminary 2D simulation of the base irradiation prior to the REPNa8 pulse test was first performed with ALCYONE. Note that ALCYONE ensures a continuity in the physical and material models between base irradiation and RIA calculations. There is therefore no specific initialization of the variables prior to pulse simulations (fragment relocation, intragranular or intergranular gas bubbles, pellet cracking…). In particular, the pulse t_{0} pelletclad gap is close to 2 μm and is therefore not artificially closed as it is the case in most of the transient fuel performance codes.
The REPNa8 pulse test was then simulated with ALCYONE. The hoop stress distribution in the cladding calculated 5 ms before and at the time of the microphone event related to the blister cracking (average fuel enthalpy 44 cal/g) are shown in Figure 3 (at PPN). The stresses are maximum in front of the pellet fragment symmetry axis where the pelletclad gap was minimum at the beginning of the pulse. They reach 170–210 MPa and are therefore of the same order as the hydride blister tensile strength deduced from the PROMETRA tests. The stress level is however too small to induce significant plastic strains. The temperature of the clad external wall does not exceed 320 °C at the time of the microphone event.
The 2D simulation is then carried on assuming the complete failure of a hydride blister of half the clad wall thickness (50%). As explained in Section 2, the boundary conditions (u_{t } = 0) are partly released on the clad line situated in front of the pellet fragment symmetry plane. It results in the opening of the fictitious blister crack with a bending moment on the clad inner surface, as shown in Figure 4. This bending moment leads to the local reopening (during the pulse) of the pelletclad gap on a circumference of ∼400 μm. Thinning of the remaining clad wall is also induced by the blister cracking.
The localization of plastic strains in a band making an angle of ∼45° with the radial direction seems consistent with the reopening of the pelletclad gap. Plastic strains develop between the blister crack tip and the first location where the pellet is still in contact with the cladding. The 45° bifurcation observed after RIA pulse tests is characteristic of a ductile failure in the plane of the maximum shear stresses. The qualitative agreement between our simulation and posttest metallographic observations is illustrated in Figure 5.
Overall, the introduction of a 50% thick hydride blister in the 2D calculation has some impact on the (average) clad outer diameter variation during the test. The loss of stiffness induced by the blister cracking leads to an increase of the (average) clad outer diameter from 0.4% to 0.6%. The latter is to be compared with the 0.5% residual strains estimated posttest from the metallographic radial cut close to the PPN [3].
Main characteristics of the CABRI REPNa8 test.
Fig. 3
Hoop stresses (in MPa) calculated in the cladding 5 ms (left) and at the time of the microphone event attributed to the cracking of a hydride blister at PPN (right). 
Fig. 4
Local reopening of the pelletclad gap during the pulse transient and distribution of the clad equivalent plastic strains (clad displacements are multiplied by a factor 5). 
Fig. 5
Clad state at PPN after the REPNa8 test and calculated equivalent plastic strains at the end of the 2D simulation. 
5 Impact of hydride blister depth on clad strains
The 2D simulation of the REPNa8 test has been used to study the impact of the hydride blister depth on the clad strains. The onset of blister cracking (at the time of the microphone event) has not been modified since the blister is assumed to behave as the rest of the cladding. Only the number of nodes where the boundary conditions are released has been changed in the simulations. As illustrated in Figure 6, six configurations with blisters depths equal to 25, 50, 60, 70, 80 and 90% of the clad wall thickness, have been considered.
As can be expected, the increase of the blister depth leads to increasing plastic strains in the 45° band. To compare the strain levels in the uncracked clad ligament situated beneath the hydride blister, the average plastic hoop strain in the 45° band has been used. It was preferred to the maximum plastic strain which obviously depends greatly on the mesh refinement and to the average hoop strain in the cladding which does not account much for the pronounced strain localization at the blister crack tip.
Figure 7 depicts the clad zone where the average plastic hoop strain in the 45° band is calculated. In each concentric ring of the clad mesh, the maximum plastic hoop strain is determined. Since an equidistant mesh is used in the radial direction, the sum of the maximum plastic hoop strains is then divided by the number of elements situated beneath the blister crack tip. The hoop strain was chosen because it can be compared directly to the local hoop strain measured from wall thinning at the blister crack tip [15]. This is not the case of the deviatoric plastic strain even if the latter is more relevant to assess the extension of clad damage in hydrided cladding [16].
The calculated time evolutions of the average plastic hoop strain in the 45° shear band are plotted in Figure 8. The times of the microphone events associated to blister cracking and to the rod failure are indicated. If we assume that the 50% deep hydride blister found at PPN is representative of the initial state of the REPNa8 cladding, it appears that the allowable maximum plastic strain in the clad ligament beneath the hydride blister is of the order of 0.5%.
According to our thermomechanical simulation of the REPNa8 test, failure of the rod is reached at a very low average plastic strain level in the 45° shear band (0.5%). This value can be compared to the 1 to 5% fracturetip wall thinning estimated by Chung and Kassner from the REPNa1 posttest metallographies [17] and to a lesser extent to the 3 to 10% local plastic strains measured by Hermann et al. [15] from burst tests performed at 350 °C on irradiated Zy4 cladding samples containing large hydride lenses (40–50% of the clad wall thickness). These local strains were estimated from the local thinning of the clad ligaments situated beneath the hydride lenses and are therefore close to our analysis of the calculated plastic strains. The burst tests performed by Hermann et al. [15] were pressure driven tests which might not lead to experimental strains directly comparable to our calculated strains. It may be argued that the viscoplastic model of Le Saux et al. [11] does not account for the development of cavities and voids in the highly strained parts of the clad ligament and hence does not allow to capture the softening induced by material damage. The treatment of softening in the framework of continuum mechanics requires the development of more sophisticated damage or micromechanically based models [16,18].
In Figure 8, the time evolution of the average plastic hoop strain in the 45° shear band has been plotted for thicker (60–70–80–90%) and thinner (25%) hydride blisters. Assuming that the failure of the clad ligament beneath the blister can be related to an average plastic strain of 0.5%, the potential variation in failure time and hence in average fuel enthalpy can be deduced from the simulations. The impact of blister depth on the fuel enthalpy at failure can hence be summarized in Figure 9. In our simulations of the REPNa8 test, the first microphone event associated to blister cracking occurs at an enthalpy of 42 cal/g. In case of a throughwall blister (100%), it gives the enthalpy at failure. In the case of the assumed 50% thick hydride blister of REPNa8, the maximum enthalpy reaches 75 cal/g. This value can be compared to the 78 cal/g obtained from a SCANAIR simulation of the REPNa8 test (Tab. 1). Between 50 and 90%, the evolution of the enthalpy at failure is close to linearity, reflecting the almost constant rate of straining by the pellet and the small plastic strains at failure.
Fig. 6
Equivalent plastic strains calculated at the time of the REPNa8 fuel rod loss of tightness in function of the hydride blister depth (in % of the clad wall thickness). 
Fig. 7
Location where the average plastic hoop strain in the 45° band is calculated. 
Fig. 8
Calculated evolution of the average plastic hoop strain in the 45° shear band during the REPNa8 simulation in function of the hydride blister depth (in % of the clad thickness). 
Fig. 9
Calculated evolution of the enthalpy at failure in function of the hydride blister depth (in % of the clad wall thickness), based on the REPNa8 test conditions. 
6 2D model of the whole cladding circumference with a single hydride blister
In spite of its qualitative agreement with postRIA observations of failed rods, the 2D simulations performed with ALCYONE tend to underestimate the plastic strains in the clad ligament beneath the blister for the following reasons: the model is constrained by the fragmentation of the pellet in 8 identical pieces which implies that the simulation represents the failure of the whole circumference of the clad tube with 8 identical blisters; the blister crack is introduced on a plane of symmetry, meaning that two identical 45° shear bands develop at the same time. To improve the simulation of rod failure, a finite element model of the whole circumference of the cladding tube with a single hydride blister has been undertaken. The FE mesh is illustrated in Figure 10.
The mesh in the vicinity of the hydride blister is much more refined than in the ALCYONE calculation since only the cladding is considered. The cracking of the blister is assumed to be slightly dissymmetric (it corresponds to the blister crack angular position in the REPNa8 test, see Fig. 5), the dissymmetry being a possible input parameter via the angles θ_{1} and θ_{2}. The loading consist in the prescribed radial displacement of the fuel external surface calculated by ALCYONE. Friction and slipping between the pellet and cladding are accounted for by modeling the fuel external and clad internal surfaces by distinct elements. In this respect, the fuelclad gap reopening observed in ALCYONE simulations (Fig. 4) can be reproduced. The time evolutions of the clad external and internal temperatures are also extracted from ALCYONE simulations.
Crosscomparisons with ALCYONE simulations have shown that the time evolution of the average plastic hoop strain in the 45° shear band beneath the blister is correctly reproduced by the present calculation if only one eighth of the pellet is considered and if the blister is not dissymmetric. The mesh refinement was found to have very little impact which confirms that the chosen plastic strain criterion is numerically sound. It was also checked that the angular position of the blister crack had no impact on the results (see Fig. 11) meaning that a prescribed radial displacement based on the average of the fuel external surface is adequate. In this respect, the important stress concentration in the cladding in front of the pellet radial cracks usually considered of great importance in PCI transient [7,9] (power ramps) appears of secondary importance for RIA calculations.
The crack in the blister is meshed explicitly with distinct facing nodes. At the beginning of the calculation, the displacements of the facing nodes are equal. At the time of the microphone event related to the blister cracking, these conditions are released. Figure 10 illustrates the dissymmetric development of plastic strains in the simulation with a 50% thick blister (θ_{1} = 5°, θ_{2} = 10°, blister length on clad outer surface ∼1.3 mm). Interestingly, very high plastic strains tend to develop at the blister crack tip but also along the blisterclad interface which might be the reason for the blistercrack decohesion observed in Figure 5.
From the calculation, it appears that the plastic strain in the 45° shear band reaches in this case 2.8% at the time of failure of the REPNa8 test, to be compared to the previous estimate with ALCYONE of the plastic strain at failure (0.5%, see Fig. 8). A realistic approach to clad failure induced by hydride blister cracking during a RIA obviously requires the development of a simulation tool able to model the whole circumference of the cladding tube and the evolving PCMI. The impact of blister geometry (length on outer clad surface, dissymmetry of the crack) has not been studied yet in spite of its potential great importance on plastic strains. Moreover, the limited axial length of the blisters was not considered in this study since only 2D simulations were performed. It might as well contribute to the calculation of greater plastic strains more consistent with experimental measures.
Fig. 10
Finite Element mesh of the 2D simulation of the whole cladding circumference with a single dissymetrically cracked blister. Calculated distribution of equivalent plastic strains in the simulation. 
Fig. 11
Plastic strain distributions calculated at the end of the REPNa8 test in case the blister crack is situated in front of the plane of fracture (left) or the plane of symmetry (right) of the pellet fragment. 
7 Conclusions
In this paper, the thermomechanical response of a fuel rod containing a fictitious hydride blister has been studied by 2D plane strain simulations of RIA. The blister was assumed to behave as the rest of the cladding material till a prescribed hoop stress is reached on the outer clad surface. The stress criterion is based on the Finite Element analysis of a PROMETRA laboratory hoop tensile test performed on a highly irradiated sample with a throughwall hydride blister. The 2D simulations of the REPNa8 pulse test performed with ALCYONE led to a qualitatively good representation of the mixed failure mode encountered in RIA pulse transients on highly corroded fuel rods: the brittle failure of the prepulse 50% thick hydride blister was followed by the development of a diagonally oriented plastic shear band in the remaining clad ligament. The average plastic hoop strain in the shear band at the time of failure of the REPNa8 test (∼0.5%) was then used to quantify the enthalpy at failure in fuel rods with blister of increasing depth (25 to 90% of the clad thickness).
FE simulations of the whole cladding circumference with a single hydride blister were then performed to study the impact of the blister geometry and of the crack position on the results. In case of a dissymmetric radial crack, high plastic strains were obtained at the blister crack tip but also along the blisterclad interface, which might be the reason for the blisterclad decohesion that can be observed after RIA tests.
Acknowledgments
The authors would like to thank AREVA and EDF for the financial and technical support to this work.
References
 T. Fuketa, H. Sasajima, Y. Mori, K. Ishijima, Fuel failure and fission gas release in high burnup PWR fuels under RIA conditions, J. Nucl. Mater. 248, 249 (1997) [CrossRef] [Google Scholar]
 T. Fuketa, T. Sugiyama, Nuclear fuel behavior during RIA, in OCDE/NEA Workshop, Paris, France (2009) [Google Scholar]
 J. Papin et al., in Eurosafe Meeting, Paris, France (2003) [Google Scholar]
 J. Papin, B. Cazalis, J.M. Frizonnet et al., Summary and interpretation of the CABRI REPNa program, Nucl. Technol. 157, 230 (2007) [CrossRef] [Google Scholar]
 V. Georgenthum et al., in WRFPM Conference , Seoul, Korea (2008) [Google Scholar]
 J. Desquines et al., in Proceedings of the ASTM conference on Zirconium in the Nuclear Industry , Stockholm, Sweden (2004) [Google Scholar]
 B. Michel, C. Nonon, J. Sercombe, F. Michel, V. Marelle, Simulation of pelletcladding interaction with the pleiades fuel performance software environment, Nucl. Technol. 182, 124 (2013) [CrossRef] [Google Scholar]
 CAST3M, http://wwwcast3m.cea.fr/ [Google Scholar]
 J. Sercombe, I. Aubrun, C. Nonon, Power ramped cladding stresses and strains in 3D simulations with burnupdependent pellet–clad friction, Nucl. Eng. Des. 242, 164 (2012) [CrossRef] [Google Scholar]
 J. Sercombe et al., in TopFuel Conference , Orlando, Florida, USA (2010) [Google Scholar]
 M. Le Saux, J. Besson, S. Carassou, C. Poussard, X. Averty, A model to describe the anisotropic viscoplastic mechanical behavior of fresh and irradiated Zircaloy4 fuel claddings under RIA loading conditions, J. Nucl. Mater. 378, 60 (2008) [CrossRef] [Google Scholar]
 B. Cazalis, J. Desquines, C. Poussard et al., The PROMETRA program: Fuel cladding mechanical behavior under high strain rate, Nucl. Technol. 157, 215 (2007) [CrossRef] [Google Scholar]
 A.H. de Ménibus, Q. Auzoux, O. Dieye et al., Formation and characterization of hydride blisters in Zircaloy4 cladding tubes, J. Nucl. Mater. 449, 132 (2014) [CrossRef] [Google Scholar]
 V. Macdonald, D. Le Boulch, A.H. de Ménibus, J. Besson, Q. Auzoux, J. Crépin, T. Le Jolu, Fracture of Zircaloy4 fuel cladding tubes with hydride blisters, Procedia Mater. Sci. 3, 233 (2014) [CrossRef] [Google Scholar]
 A. Hermann et al., in Proceedings of the 15th ASTM conference on Zirconium in the Nuclear Industry , Sunriver, Oregon, USA (2007) [Google Scholar]
 M. Le Saux, J. Besson, S. Carassou, A model to describe the mechanical behavior and the ductile failure of hydrided Zircaloy4 fuel claddings between 25 °C and 480 °C, J. Nucl. Mater. 466, 43 (2015) [CrossRef] [Google Scholar]
 H.M. Chung, T.F. Kassner, Cladding metallurgy and fracture behavior during reactivityinitiated accidents at high burnup, Nucl. Eng. Des. 186, 411 (1998) [CrossRef] [Google Scholar]
 Y. Udagawa, T. Mihara, T. Sugiyama, M. Suzuki, M. Amaya, Simulation of the fracture behavior of Zircaloy4 cladding under reactivityinitiated accident conditions with a damage mechanics model combined with fuel performance codes FEMAXI7 and RANNS, J. Nucl. Sci. Technol. 51, 208 (2014) [CrossRef] [Google Scholar]
Cite this article as: Jérôme Sercombe, Thomas Helfer, Eric Federici, David Leboulch, Thomas Le Jolu, Arthur Hellouin de Ménibus, Christian Bernaudat, 2D simulation of hydride blister cracking during a RIA transient with the fuel code ALCYONE, EPJ Nuclear Sci. Technol. 2, 22 (2016)
All Tables
All Figures
Fig. 1
Mesh and mechanical boundary conditions in the 2D scheme of ALCYONE. 

In the text 
Fig. 2
Hoop stresses calculated at failure time during the PROMETRA test 2468 (left: friction coefficient 0.1, right: friction coefficient 0.4). 

In the text 
Fig. 3
Hoop stresses (in MPa) calculated in the cladding 5 ms (left) and at the time of the microphone event attributed to the cracking of a hydride blister at PPN (right). 

In the text 
Fig. 4
Local reopening of the pelletclad gap during the pulse transient and distribution of the clad equivalent plastic strains (clad displacements are multiplied by a factor 5). 

In the text 
Fig. 5
Clad state at PPN after the REPNa8 test and calculated equivalent plastic strains at the end of the 2D simulation. 

In the text 
Fig. 6
Equivalent plastic strains calculated at the time of the REPNa8 fuel rod loss of tightness in function of the hydride blister depth (in % of the clad wall thickness). 

In the text 
Fig. 7
Location where the average plastic hoop strain in the 45° band is calculated. 

In the text 
Fig. 8
Calculated evolution of the average plastic hoop strain in the 45° shear band during the REPNa8 simulation in function of the hydride blister depth (in % of the clad thickness). 

In the text 
Fig. 9
Calculated evolution of the enthalpy at failure in function of the hydride blister depth (in % of the clad wall thickness), based on the REPNa8 test conditions. 

In the text 
Fig. 10
Finite Element mesh of the 2D simulation of the whole cladding circumference with a single dissymetrically cracked blister. Calculated distribution of equivalent plastic strains in the simulation. 

In the text 
Fig. 11
Plastic strain distributions calculated at the end of the REPNa8 test in case the blister crack is situated in front of the plane of fracture (left) or the plane of symmetry (right) of the pellet fragment. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.