Mathematical description of WAFER-1, a three-dimensional code for LWR fuel performance analysis

This article describes in detail the mathematical formulation used in the WAFER-1 code, which is presently used for three-dimensional analysis of LWR fuel pin performance. The code aims at a prediction of the local stress-strain history in the cladding, especially with regard to the ridging phenomenon. To achieve this, a clad model based on shell theory has been developed. This model interacts with a detailed finite difference pellet model which treats radial and transversal cracking in the pellet in a deterministic way, based on certain assumptions with respect to the cracking pattern. Pellet and clad creep are taken into account. The inner core of the pellet, bounded by a specified isotherm, may be treated as a viscous material. Axial force exchange between pellet and clad is also included. The axial loading is distributed on the pellet end face with due regard to any pellet dishing. An arbitrary power history may be used as input to the model.     

Modeling of fission gas swelling in LWR UO2 fuel under irradiation

An understanding of the behavior of fission gas in uranium dioxide (UO₂) fuel is necessary for the prediction of the performance of fuel rods under irradiation. A mechanistic model for matrix swelling by the fission gas in LWR UO₂ fuel is presented. The model takes into account intragranular and intergranular fission gas bubbles behavior as a function of irradiation time, temperature, fission rate and burn-up. The intragranular bubbles are assumed to be nucleated along the track of fission fragments, which play the dual role of creator and destroyer of intragranular bubbles. The intergranular bubble nuclei is produced until such time that a gas atom is more likely to be captured by an existing nucleus than to meet another gas atom and form a new nucleus. The capability of this model was validated by a comparison with the measured data of fission gas behavior such as intragranular bubble size, bubble density and total fuel swelling. It was found that the calculated intragranular bubble size and density are in reasonable agreement with the measured results in a broad range of average fuel burn-ups 6–83 GW d/tU. Especially, the model correctly predicts the fuel swelling up to a burn-up of about 70 GW d/tU.     

FEMAXI-III,a computer code for fuel rod performance analysis

This paper presents a method of fuel rod thermal-mechanical performance analysis used in the FEMAXI-III code. The code incorporates the models describing thermal-mechanical processes such as pellet-cladding thermal expansion, pellet irradiation swelling, densification, relocation and fission gas release as they affect pellet-cladding gap thermal conductance. The code performs the thermal behavior analysis of a full-length fuel rod within the framework of one-dimensional multi-zone modeling. The mechanical effects including ridge deformation is rigorously analyzed by applying the axisymmetric finite element method. The finite element geometrical model is confined to a half-pellet-height region with the assumption that pellet-pellet interaction is symmetrical. The 8-node quadratic isoparametric ring elements are adopted for obtaining accurate finite element solutions. The Newton-Raphson iteration with an implicit algorithm is applied to perform the analysis of non-linear material behaviors accurately and stably. The pellet-cladding interaction mechanism is exactly treated using the nodal continuity conditions. The code is applicable to the thermal-mechanical analysis of water reactor fuel rods experiencing variable power histories.     

Calculational procedure for determining creep collapse of LWR fuel rods

Observed collapses in pressurized water reactor fuel rods have been attributed to the radiation enhanced creep of Zircaloy cladding into regions where separations in the fuel pellet stack have occurred. A computer code, COLAPX, has been written to determine the growth of ovality and the ultimate collapse of fuel rod cladding under reactor operating conditions. This paper describes the theoretical bases of this code, the finite element formulation used, the constitutive relations between the displacement fields and the element forces, and the radiation, temperature and stress dependent material model for creep of Zircaloy tubing. Comparisons of the creep rate predictions and of the ovality predictions with data from irradiated tubes and fuel cladding are presented.     

Irradiation creep and interrelation with swelling in austenitic stainless steels

The available experimental data on irradiation-induced creep in austenitic stainless steels are summarized and the existing theories reviewed. Attention is paid to the influence of material composition and pretreatments on irradiation creep. In particular the stress, flux, fluence and temperature dependencies are reported and possible correlations of irradiation creep with the microstructural evolution, the swelling behaviour and the precipitation kinetics of the materials are outlined. The consequences of stress effects connected with swelling for the irradiation-creep behaviour, especially the stress-dependence, are discussed.     

MIPAC, a computer code for fuel performance analysis by the finite element method

An axisymmetric finite element computer code named MIPAC has been developed for analysis of the mechanical interaction behaviour between a fuel pellet and cladding. This computer code can deal with elastoplasticity of the pellet and cladding materials, creep effects for the both materials, pellet-cladding and pellet-pellet contact problems, hot pressing effect of the fuel pellet, fuel pellet cracking, and the cracked pellet's stiffness. A cyclical boundary condition is introduced to deal with one pellet length instead of the full-size fuel rod. The contact problems are solved without a fictitious contact element. In the fuel pellet cracking model the crack opening and closing behaviour under arbitrary power changes can be treated by introducing five kinds of crack modes. Mismatch of irregular crack surfaces is taken into account in the evaluation of the cracked pellet's stiffness. Finally, calculated results are compared with experimental data to show validity of the computer code.     

Irradiation devices for the study of creep and swelling in ceramic fuels

In fuel element design for advanced nuclear reactors perfect knowledge of fuel behaviour under irradiation plays a decisive role, above all for long service lives and high burnups. Therefore, the development of fast breeder fuel elements within the framework of the Karlsruhe Fast Breeder Project included various irradiation rigs which allow continuous measurement during irradiation of fuel specimen creep and swelling. A survey is presented of some of these irradiation rigs.     

Coupled thermal structural analysis of LWR vessel creep failure experiments

Considering the hypothetical core melt down scenario for a light water reactor (LWR) the failure mode of the reactor pressure vessel (RPV) has to be investigated to determine the loadings on the containment. The failure of reactor vessel retention (FOREVER)-experiments, currently underway, are simulating the thermal and pressure loadings on the lower head for a melt pool with internal heat sources. Due to the multi-axial creep deformation of the vessel with a non-uniform temperature field these experiments are an excellent source of data for validation of numerical creep models. Therefore, a finite element (FE) model has been developed based on a commercial multi-purpose code. Using the computational fluid dynamics (CFD) module the temperature field within the vessel wall is evaluated. The transient structural mechanical calculations are performed using a new numerical approach, which avoids the use of a single creep law employing constants derived from the data for a limited stress and temperature range. Instead of this a three-dimensional array is developed where the creep strain rate is evaluated according to the values of the actual total strain, temperature and equivalent stress. Care has to be exercised performing post-test calculations particularly in the comparisons of the measured data and the numerical results. Considering the experiment FOREVER-C2, for example, the recorded creep process appears to be tertiary, if a constant temperature field is assumed. But, small temperature increase during the creep deformation stage could also explain the observed creep behavior. Such considerations provide insight and better predictive capability for the vessel creep behavior during prototypic severe accident scenarios.     

Modeling of gap conductance for LWR fuel rods applied in the BISON code

ABSTRACT

A new gap conductance model is proposed in this study as a combination of Toptan’s model and the Ross-Stoute model. A variance-based sensitivity analysis is performed to understand how simulation results depend on all input parameters of the proposed model. Additionally, new modeling options (e.g. fill gas thermal conductivity, temperature jump distance, thermal accommodation coefficient, etc.) are added into the nuclear fuel performance code, BISON. The need for further investigation of the gap heat transfer between fuel and cladding in BISON motivated this study to evaluate its impact on the code’s predictions. New gap conductance modeling is proposed. A series of integral-effects validation tests is performed: (1) to demonstrate the impact of the proposed model on the code’s fuel temperature predictions at the beginning of life and through the reactor’s life; (2) to ensure that the proposed model is capable of accurately modeling gap heat transfer characteristics in real-world problems; and (3) to investigate the impact of the estimation of fission gas release on the fuel temperature predictions with the proposed model. The results indicate that the proposed gap conductance model improves BISON’s predictions.     

Enhancing the ABAQUS thermomechanics code to simulate multipellet steady and transient LWR fuel rod behavior

A powerful multidimensional fuels performance analysis capability, applicable to both steady and transient fuel behavior, is developed based on enhancements to the commercially available ABAQUS general-purpose thermomechanics code. Enhanced capabilities are described, including: UO₂ temperature and burnup dependent thermal properties, solid and gaseous fission product swelling, fuel densification, fission gas release, cladding thermal and irradiation creep, cladding irradiation growth, gap heat transfer, and gap/plenum gas behavior during irradiation. This new capability is demonstrated using a 2D axisymmetric analysis of the upper section of a simplified multipellet fuel rod, during both steady and transient operation. Comparisons are made between discrete and smeared-pellet simulations. Computational results demonstrate the importance of a multidimensional, multipellet, fully-coupled thermomechanical approach. Interestingly, many of the inherent deficiencies in existing fuel performance codes (e.g., 1D thermomechanics, loose thermomechanical coupling, separate steady and transient analysis, cumbersome pre- and post-processing) are, in fact, ABAQUS strengths.     

FRAP-T6: A computer code for the transient analysis of oxide fuel rods

The FRAP-T6 computer code was developed to model the transient performance of light water reactor fuel rods during reactor transients ranging from mild operational transients to large break loss-of-coolant accidents. The code models all of the thermal, structural, and chemical phenomena needed for the complete evaluation of light water reactor fuel rod performance. The code was developed using rigorous quality assurance procedures and a large assessment data base. The results of assessment show that the code accurately models the response of light water reactor fuel rods.     

SAMURA—A fast computer code for reactor fuel pin performance modelling and analysis

A very fast integral numerical computer code for the modelling of transient and steady-state thermal and mechanical behaviour of Zircaloy-clad UO₂ fuel pins in water reactors has been developed. The computational technique which determines the stress and deformation state of the fuel pin is based upon an extremely efficient finite difference scheme, i.e. the non-linear terms in the constitutive equations which produce a non-linear system of equations have been linearised using a Taylor expansion technique coupled with a very sophisticated error minimization algorithm and then solved with great accuracy. An improved numerical method has also been developed for the fast and efficient solution of the transient heat conduction equation. In this way a very stable and economical one-dimensional code (with appropriate provisions made for its conversion to a quasi two-dimensional code) has been obtained. The physical processes included are thermo-elastic deformation, thermal and irradiation creep, plasticity, fission gas swelling and release, formation of cracks in the fuel, hot pressing, densification, pore migration and dish or central void filling. Here the mathematical basis of SAMURA is presented along with some preliminary calculations and benchmarkings. It is concluded that SAMURA is quite fast indeed, converges to accurate results and within the margins of the error criterion chosen has very reasonable computer demands. It is also stable under all conditions tested.     

Irradiation swelling in UO2-SiO2 fuel under low axial stress

Samples of UO₂-SiO₂ nuclear fuel were irradiated in special capsules, developed for the recording of length changes by very low axial stress (1 kg f/cm²) without any radial restriction. The results of irradiations at temperatures between 200 and 1300°C and burnups up to 20 000 MWd/t are presented. The experimental results emphasize a long-term compaction which is an irradiation-induced creep, and an important swelling occurring at higher burnups for lower temperatures. Between 1100 and 1200°C and at 1 atm this swelling appears at about 10 000 MWd/t.     

高温气冷堆包覆燃料颗粒破损机制及失效模型

高温气冷堆的燃料元件的基本构成单元是全陶瓷型的包覆燃料颗粒,其性能决定了高温气冷堆的安全性。除了传统的辐照实验检测外,建立理论模型对其研究具有重要的意义。本文主要介绍了TRI-SO型包覆燃料颗粒的结构及破损机制,以及国外现有的几个主要模型的基本假设,计算原理和特点,通过对比几个模型的优缺点,提出今后研究的方向。     

Low temperature swelling and densification properties of LWR fuels

Analysis of recent data reported in the literature on low temperature swelling and densification behavior of light water reactor (LWR) fuel suggests that, at low irradiation temperatures, the extent of irradiation induced densification shows a simple exponential dependence on burnup with a rate constant that is insensitive to temperature and flux level. These data are also consistent with the view that the total porosity controls the kinetics of irradiation-induced densification process. In the same temperature range the rate of swelling was found to be constant to a burnup near 50 GWD/MTU, with a value of approximately 1% per 10 GWD/MTU. The significance of these results is discussed in terms of current theories.