Preliminary modelling of crack nucleation and propagation in SiC/SiC accident-tolerant fuel during routine operational transients using peridynamics

Abstract

Silicon carbide fibre in silicon carbide matrix composites (SiC/SiC) are a promising cladding for use in accident tolerant fuels (ATF) in current light water reactor (LWR) designs. However, as they are a radically different material from current metal clads, current thermomechanical simulation methods struggle to accurately predict their behaviour, especially regarding the potential development of cracks. Thus, a new peridynamic model for SiC/SiC cladding has been developed in the Abaqus finite element code. The material model was isotropic and considers matrix cracking and fibre pull-out. The thermal expansion, swelling and the degradation of the thermal conductivity are modelled under typical LWR irradiation conditions. The swelling on the outer surface is predicted to be greater than the inner surface due to the lower irradiation temperature, causing a tensile stress on the inside of the cladding; tension being more challenging for a ceramic than a metal. This stress increases during the decrease in power at the start of a typical pressurised water reactor refuelling outage and causes microcracking of the matrix on the cladding inner surface. In models without fibres, cracks would propagate through the cladding. If fibres are modelled, matrix cracking will extend to a depth of around 20% through the cladding from the inner surface, which is unlikely to be an acceptable design. If an inner monolith of SiC is additionally modelled, cracking propagates through the monolith and acts as a stress raiser for matrix cracking in the composite, and therefore does not constitute a design improvement. If an outer SiC monolith is modelled, fibre pull-out strain on the inner surface of the cladding was increased by just under 70%. No cracks are predicted in an outer monolith which may therefore remain gas-tight and thus a more suitable design. These predictions are consistent with experimental findings. Highlights• Irradiation of accident tolerant SiC/SiC cladding is modelled using peridynamics.• Due to swelling, a tensile stress is predicted on the inside of the cladding. • In the absence of a monolith, micro-cracking extends to 20% from the inner surface.• Cracks in an inner monolith are a stress raiser in the composite.• No cracks are predicted in an outer monolith, hence it's the most promising design.

Application of Weibull Fracture Strength Distributions to Modelling Crack Initiation Behaviour in Nuclear Fuel Pellets Using Peridynamics

“…The two-dimensional bond-based implementation of peridynamics in the commercial finite element software Abaqus, was initially presented by Macek and Silling [21]. This implementation has subsequently developed by the authors [19,[22][23][24] and has recently been extended to model solids with failure described by Weibull statistics [19]. The approach is described in detail elsewhere [19,23], and only a brief overview is presented here.…”

“…This implementation has subsequently developed by the authors [19,[22][23][24] and has recently been extended to model solids with failure described by Weibull statistics [19]. The approach is described in detail elsewhere [19,23], and only a brief overview is presented here. Figure 2 shows a 2D section of material decomposed into a series of points, termed 'material points' in peridynamics.…”

“…Also visible in Figure 14 is a tendency for some, but not all, cracks to initiate on geometric features of the mesh. This is in part a feature of modelling a curved surface using a rectangular mesh and is in-part offset by the use of the edge correction term detailed in [23].…”

Peridynamic Modelling of Cracking in TRISO Particles for High Temperature Reactors

“…In peridynamics, objects are decomposed into a number of 'material points'. The behaviour of each material point is dependant upon those in a horizon surrounding the material point, not just its nearest neighbours [37] . An example of this is shown Fig.…”

“…In the Abaqus implementation presented here, each 'material point' is represented by a node, and truss elements represent the 'bonds' linking the material points. This approach is discussed in more detail in [37] .…”