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.     

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.     

BEAF - a computer program for analysis of light water reactor fuel rod behavior

BEAF - a computer program for analysis of light water reactor fuel rod behavior was developed. The BEAF code, which is appropriate for on-line prediction of fuel rod behavior, can analyze fuel rod thermal and mechanical behaviors using the axisymmetric, plane strain approximation and finite difference method to realize a fast running time.In the mechanical analysis, a new cracked pellet compliance model is introduced, in which pellet cracking and crack healing, pellet initial relocation, modified elastic moduli of a cracked fuel pellet, and stress dependent hot pressing are considered. Adding to those capabilities, fission gas flow and diffusion in the fuel-clad gap are analyzed to take into account the slow fission gas dilution effect on the gap conductance during power ramp.     

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.     

FEMAXI-IV: A computer code for the analysis of fuel rod behavior under transient conditions

A fuel rod behavior code FEMAXI-IV, presently under development, is an improved version of the FEMAXI-III code for the analysis of fuel rod behavior under transient conditions. To apply the FEMAXI-III code to transient conditions, the following additional models have been incorporated into the FEMAXI-III code: transient heat transfer model: axial gas mixing model; diffusion-type fission gas release model. This paper summarizes the above additional models, and the comparison of the FEMAXI-IV calculations with the experimental data.     

Numerical treatment of the cladding mechanics in the fuel rod computer code IAMBUS

IAMBUS (INTERATOM Model for Burn-Up Studies on Fuel Rods) is a digital computer code for the thermal and mechanical design, in-pile performance prediction and post-irradiation analysis of fuel rods. The mechanical analysis of the cladding is approximated by a state of generalized plane strain. The analysis includes routines for plasticity, creep and swelling due to void nucleation and growth. The numerical integration of these equations is described in detail and the accuracy of the results is discussed.     

The SABRE code for fuel rod cluster thermohydraulics

This paper describes the capabilities of the SABRE code for the calculation of single phase and two phase fluid flow and temperature in fuel pin bundles, discusses the methods used in the modelling and solution of the problem, and presents some results including comparison with experiments.The SABRE code permits calculation of steady-state or transient, single or two phase, flows and the geometrical options include general representation of grids, wire wraps, multiple blockages, bowed pins, etc. Transient flows may be calculated using either semi-implicit or fully implicit time solution methods and the temperature distributions within the fuel pins are determined as well as the velocity and temperature of the coolant. Two phase flows are calculated using a homogeneous boiling model, with the possibility of a specified slip between the two phases. General inlet boundary conditions are available (including pressure, velocity, total mass flow) and these may vary linearly with time; the outlet boundary condition is taken as constant pressure. The treatment of grids allows for irreversible head losses at entry and exit. The wire wrap model introduces a grid resistance tensor with its principal axes along and perpendicular to the wire, resulting in a very satisfactory modelling of inducement of swirl.The derivation and solution of the difference equations is discussed. Emphasis is given to the derivation of the spatial differences in triangular subchannel geometry, and the use of central, upwind or vector upwind schemes. The method of solution of the difference equations is described for both steady state and transient problems. Together with these topics we consider the problems involved in turbulence modelling and how it is implemented in SABRE. This includes supporting work with a fine scale curvilinear coordinate programme to provide turbulence source data. The problem of modelling boiling flows is discussed, with particular reference to the numerical problems caused by the rapid density change on boiling.The final part of the paper presents applications of the code to the analysis of blockage situations, the study of flow and power transients and analysis of natural circulation within clusters to demonstrate the scope of the code and compare with available experimental results. The comparisons include the calculation of a flow pressure drop characteristic of a boiling channel showing the Ledinegg instability, examples of overpower and flow rundown transients which lead to coolant boiling, and calculation of natural circulation within a rod cluster.     

Application of adjoint sensitivity analysis to nuclear reactor fuel rod performance

Adjoint sensitivity analysis in nuclear fuel behavior modeling is extended to operate on the entire power history for both Zircaloy and stainless steel cladding via the computer codes FCODE-ALPHA/SS and SCODE/SS. The sensitivities of key variables to input parameters are found to be highly non-intuitive and strongly dependent on the fuel-clad gap status and the history of the fuel during the cycle. The sensitivities of five key variables, clad circumferential stress and strain, fission gas release, fuel centerline temperature and fuel-clad gap, to eleven input parameters are studied. The most important input parameters (yielding significances between 1 and 100) are fabricated clad inner and outer radii and fuel radius. The least important significances (less than 0.01) are the time since reactor start-up and fuel burnup densification rate. Intermediate to these are fabricated fuel porosity, linear heat generation rate, the power history scale factor, clad outer temperature, fill gas pressure and coolant pressure. Stainless steel and Zircaloy have similar sensitivities at start-up but these diverge as burnup proceeds due to the effect of the higher creep rate of Zircaloy which causes the system to be more responsive to changes in input parameters. The value of adjoint sensitivity analysis lies in its capability of uncovering dependencies of fuel variables on input parameters that cannot be determined by a sequential thought process.     

Analysis of MOX fuel behavior in reduced-moderation water reactor by fuel performance code FEMAXI-RM

To assess the feasibility of the 31% Pu-MOX fuel rod design of reduced-moderation water reactor (RMWR) in terms of thermal and mechanical behaviors, a single rod assumed to be irradiated in the core of RMWR up to 106 GWd/tHM has been analyzed by a fuel performance code FEMAXI-RM which is an extended version of FEMAXI-6 code. In the analysis, design specifications of fuel rod and irradiation conditions have been input, and available models of both MOX fuel and UO₂ fuel have been used as appropriate. The results are: fission gas release is several tens of percent, rod internal pressure does not exceed the coolant pressure, and the highest fuel center temperature is 2400 K, while cladding diameter increase caused by pellet swelling is within 1% strain. These predictions suggest that the MOX fuel rod integrity will be held during irradiation in RMWR, though actual behavior of MOX pellet swelling and cladding oxidation require to be investigated in detail.     

Thermo-mechanical analysis of supercritical pressure light water-cooled fast reactor fuel rod by FEMAXI-6 code

Thermo-mechanical behaviors of supercritical pressure light water cooled fast reactor (SWFR) fuel rod and cladding have been investigated by FEMAXI-6 (Ver.1) code with high enriched MOX fuel at elevated operating condition of high coolant system pressure (25 MPa) and high temperature (500 °C in core average outlet temperature). Fuel rod failure modes and associated fuel rod design criteria that is expected to be limiting in SWFR operating condition have been investigated in this fuel rod design study. Fuel centerline temperature is evaluated to be 1853 °C and fission gas release fraction is about 45% including helium production. Cumulative damage fraction is evaluated by linear life fraction rule with time-to-rupture correlation of advanced austenitic stainless steel. In a viewpoint of mechanical strength of fuel cladding against creep rupture and cladding collapse at high operation temperature, currently available stainless steels or being developed has a potential for application to SWFR. Admissible design range in terms of initial gas plenum pressure and its volume ratio are suggested for fuel rod design The stress ranges suggested by this study could be used as a preliminary target value of cladding material development for SWFR application.     

Improvements to COMETHE III-J fuel rod modeling code

Current improvements to the COMETHE fuel performance code focus on pellet-clad axial interaction and Zircaloy cladding failure predictions. Slipping and sticking between pellets and clad as well as trapped stack are evaluated. The main conclusions are that slipping with friction concerns only local effects while axial PCMI is primarily dependent on pellet expansion with a strong ‘strain biaxiality’ effect dictated by the dishing. The notion of locking prior to radial PCMI is also introduced, which explains experimental features not previously understood. Benckmarking of the version of COMETHE against ramp tests has been initiated and will enable assessment of the code capability in Zircaloy clad failure predictions.     

Plate vs rod fuel elements in a tight-lattice water reactor

In this paper, the possibility of using plate fuel elements in tight-lattice water reactors instead of rod fuel elements is explored. For this purpose we compared the thermal-hydraulic and neutronic aspects of plates with those of rod fuel elements. Because of the fact that plates have a higher mechanical strength than rods with equal dimensions, it is possible to reduce the thickness of the plates and thus reduce the centerline temperature. This fact and other thermal-hydraulic considerations demonstrate that plate fuel elements are safer than rods with respect to their resistance to LOCA.Preliminary neutronic calculations indicate that plates and rods have about the same neutronic characteristics.     

A new mechanistic code SFPR for modeling of single fuel rod performance under various regimes of LWR operation

A new mechanistic code SFPR for modeling of single fuel rod behavior under various regimes of LWR reactor operation (normal and off-normal, including severe accidents) is under development at IBRAE. The code is designed by coupling of two stand-alone mechanistic codes MFPR (for modeling of irradiated UO₂ fuel behavior and fission product release) and SVECHA/QUENCH, or S/Q (for modeling of Zr cladding thermo-mechanical and physico-chemical behavior). Both codes were initially designed for accident conditions (and for this reason, are rather mechanistic) and later extended to various normal operation conditions. On the base of thorough validation against various out-of-pile and in-pile experiments, development of an advanced fuel performance code for best estimate code calculations for both normal and off-normal LWR operation regimes is foreseen.     

FARST — A computer code for the evaluation of FBR fuel rod behavior under steady-state/transient conditions

FARST, a computer code for the evaluation of fuel rod thermal and mechanical behavior under steady-state/transient conditions has been developed. The code characteristics are summarized as follows:

1. (i) FARST evaluates the fuel rod behavior under the transient conditions. The code analyzes thermal and mechanical phenomena within a fuel rod, taking into account the temperature change in coolant surrounding the fuel rod.

2. (ii) Permanent strains such as plastic, creep and swelling strains as well as thermoelastic deformations can be analyzed by using the strain increment method.

3. (iii) Axial force and contact pressure which act on the fuel stack and cladding are analyzed based on the stick/slip conditions.

4. (iv) FARST used a pellet swelling model which depends on the contact pressure between pellet and cladding, and an empirical pellet relocation model, designated as “jump relocation model”.

The code was successfully applied to analyses of the fuel rod irradiation data from pulse reactor for nuclear safety research in Cadarache (CABRI) and pulse reactor for nuclear safety research in Japan Atomic Energy Research Institute (NSRR).The code was further applied to stress analysis of a 1000 MW class large FBR plant fuel rod during transient conditions. The steady-state model which was used so far gave the conservative results for cladding stress during overpower transient, but underestimated the results for cladding stress during a rapid temperature decrease of coolant.     

Cold-neutron tomography of annular flow and functional spacer performance in a model of a boiling water reactor fuel rod bundle

Dryout of the coolant liquid film at the upper part of the fuel pins of a boiling water reactor (BWR) core constitutes the type of heat transfer crisis relevant for the conditions of high void fractions. It is both a safety concern and a limiting factor in the thermal power and thus for the economy of BWRs. We have investigated adiabatic, air-water annular flows in a scaled-up model of two neighboring subchannels as found in BWR fuel assemblies using cold-neutron tomography. The imaging of the double suchannel has been performed at the ICON beamline at the neutron spallation source SINQ at the Paul Scherrer Institute, Switzerland. Cold-neutron tomography is shown here to be an excellent tool for investigating air-water annular flows and the influence of functional spacers of different geometries on such flows. The high-resolution, high-contrast measurements provide the spatial distributions of the coolant liquid film thickness on the fuel pin surfaces as well as on the surfaces of the spacer vanes. The axial variations of the cross-section averaged liquid hold-up and its fraction in the gas core shows the effect of the spacers on the redistribution of the two phases.