Analysis of Halden overpressure tests using the FALCON code

Transcription

1 WIR SCHAFFEN WISSEN HEUTE FÜR MORGEN G. Khvostov, W. Wiesenack Analysis of Halden overpressure tests using the FALCON code Enlarged Halden Programme Group Meeting 216

2 Contents Scope and goals of the study Main data Models and methods Results Thermal response to overpressure in IFA-61.1 (BWR) Thermal response to overpressure in IFA-61.1 (/3/5) (PWR) Response of cladding elongation to overpressure Rod FGR and cladding failure Conclusions /STARS/KG41 - ( 2 / 13 )

3 Scope and goals of the study Four tests of the IFA-61 series, using uranium dioxide as fuel material, were analyzed with the FALCON code Sample origin Local parameters after BI Test Vendor/NPP Fuel type Segment type Burn-up, MWd/kgU IFA-61.1 WSE / KKL 1X1 SVEA-96+ BWR IFA-61.1 Siemens- KWU GE / KKG 14X14 PWR A segment cut from integral rod A segment from specially designed segmented rod A segment from specially designed segmented rod Fast fluence, 1 21 n/cm2 Oxide thickness, microns Overpressure sequence, bar (5-1-15) IFA IFA-61.5 Framatom / EDF Gravelines 5 17X17 PWR General goals: verify the code s ability to treat lift-off find out governing processes and specific features of high burnup fuel behaviour explore critical conditions for rod failure caused by very high internal gas pressure /STARS/KG41 - ( 3 / 13 )

4 Cladding elongation, mm Temperature, o C Main data CL temperature in TC-section Normalized temperature Data for fuel temperature in TC-section of IFA-61.1: Measured temperature. Normalized temperature at LHRTC=11 kw/m. T fn = T f T c LHR LHR + T c 65 Normalized temperature against calculated temperature-variation across the pellet radius 6 55 T fn = T f LHR LHR Time, eph SST cladding elongation shows significant difference, likely depending on degree of bonding PWR: increasing BWR: decreasing Note: The data taken into account, but remaining beyond the scope of the paper, are: Cladding elongation during power ramps; Noise analysis for transient elongation; measurement of TH diameter Raw data for cladding elongation in: IFA IFA IFA IFA Time, eph /STARS/KG41 - ( 4 / 13 )

5 Models and methods Use of a research version of EPRI s FALCON code coupled with the mechanistic model, GRSW-A, for fission gas release and gaseous-bubble swelling RIP is treated as code input, as function of time Use of code restart to simulate re-fabrication after BI Rod overpressure, bar 35 Raw data for rod overpressure in: IFA IFA IFA IFA Time, eph FE model is designed to simulate thermal response in TC section, in consideration of its details /STARS/KG41 - ( 5 / 13 )

6 Cladding elongation, mm Models and methods TM calculation for two limit cases: without bonding and with complete bonding Cladding elongation: Measurement. Calculation without bonding. Calculation with bonding. AAFB=33% Burn-up, MWd/kgU Use of the two calculation for decomposition of measured elongation dl = AFB dl db exp db w bonding + (1 AFB) dl db w/o bonding Proposed estimation of Apparent Fractional Bonding (AFB) in the sample /STARS/KG41 - ( 6 / 13 )

7 Results 5 4 No thermal feedback up to overpressure of 15 bar. Thermal feedback starts at overpressure of 175 bar. Gap opening at overpressure of 25 bar. Rod pressure, MPa Thermal response to overpressure in IFA61.1 (BWR) Rod pressure 3 2 Coolant pressure Burnup, MWd/kgU Normalized temperature in TC-section: Measurement. Running average for measured temperature. Calculation for normalized centreline temperature. Calculation for normalized temperature-variation accross pellet radius (Tc-Ts) Burnup, MWd/kgU 6 Calculated diametral gap size. Pellet-cladding contact pressure Radial gap size, microns Contact pressure, MPa Normalized temperature, oc Normalized temperature variation, oc Burnup, MWd/kgU /STARS/KG41 - ( 7 / 13 )

8 Normalized temperature in TC-section: Measurement. Calculation w/o bonding. Calculation assuming bonding. 12 Thermal response to overpressure in IFA-61.1 (/3/5) (PWR) The model, assuming pellet-cladding bonding, inhibits gap opening. For the periods with overpressure maintained at 2 bar, 25 bar and 3 bar, a perfect agreement is obtained when bonding is assumed. Normalized temperature, oc Over-pressure, bar Results Overpressure Time, hours Splashes in temperature at power fluctuations A hysteresis: rearrangement of the fuel fragments resulting in circumferential crack opening LHGR in TC-section: 74 Raw data. FALCON input (condesed). 73 Temperature, oc LHRTC, kw/m P = 2 bar Normalized temperature, oc Normalized temperature in TC-section: Measurement. Linear fit to measurement. Running average for selected data-set. Calculation for normalized temperature variation accross pellet radius (Tc-Ts) Time, h LHRTC, kw/m Burnup, MWd/kgU /STARS/KG41 - ( 8 / 13 ) Normalized temperature variation, oc

9 Cladding elongation, mm Cladding elongation, mm Results Response of cladding elongation to overpressure Assumption of bonding in calculation maximizes impact of fuel swelling on cladding elongation under over-pressure The degree of bonding is quantified with the AFB number (current and test average), viz.: an effective fraction (%) of pellet surface bonded with cladding Test Bu dl_db-exp dl_db-nb dl_db-b AAFB CAFB STD,% The degree of pellet-cladding bonding in the BWR fuel segment is lower than in the PWR fuel segments (ca. 3 % vs 7 %) Bonding was not broken in the PWR tests throughout experiments, even though the applied maximum overpressure was very high (3 bar in IFA-61.1) Burn-up, MWd/kgU /STARS/KG41 - ( 9 / 13 ) Cladding elongation: Measurement. Calculation without bonding. Calculation with bonding. Cladding elongation: Measurement. Calculation without bonding. Calculation with bonding. T c =36 o C; AAFB=83 16 % P=175 bar AAFB=33% IFA-61.1 T c =316 o C; AAFB=61 13 % IFA Burnup, MWd/kgU

10 Temperature, K FGR, % Fission gas release, % Results Rod FGR and cladding failure The calculated FGR is a smooth continuation of the FGR at the end of BI No effect of thermally induced feedback is predicted Base Irradiation Overpressure tests Centreline temperature during BI in segments used for: IFA IFA IFA Empirical FGR threshold. Relative FGR during BI in a segment used for: IFA IFA IFA Calculated relative FGR during: IFA IFA IFA IFA Burn-up, MWd/kgU Pellet burnup, MWd/kgU According to the modelling, no considerable deterioration of rod reliability resulting from applied test conditions of overpressure can be expected, in spite of rather severe overpressure utilized in the tests. none of the test rods failed. This is consistent with the results of the calculations. Need of additional single effect experiments resulting in cladding failure due to the predicated lift-off phenomenon /STARS/KG41 - ( 1 / 13 )

11 Conclusions The FALCON code has been shown to be applicable to best-estimate predictive analysis of the overpressure tests using fuel without, or weak pellet-cladding bonding, and scoping analysis of the tests with fuels where stronger pellet-cladding bonding occurs. Significant effects of bonding and fuel cracking/relocation on the thermal and mechanical behaviour of highly over-pressurized rods, with a fuel burnup of ca. 6 MWd/kgU, have been identified. The effect of bonding is particularly pronounced in the tests with the fuel segments preirradiated in PWRs. The identified effects are beyond the current traditional fuel-rod licensing analysis methods. Derivation of appropriate acceptance criteria and model parameterization, which is important for the corresponding licensing analysis, would need additional single effect experiments resulting in cladding failure due to lift-off phenomenon /STARS/KG41 - ( 11 / 13 )

12 Wir schaffen Wissen heute für morgen The main goals of the study are basically reached: - All the tests, as planned, successfully analyzed. - Considerable effects of bonding and fuel cracking/relocation identified. - The identified effects are beyond current traditional licensing analysis. The paper and presentation outline the main results, while the whole activity is described in PSI Report TM

13 Wir schaffen Wissen heute für morgen My thanks go to ESB OECD HRP

14 Back-up slides /STARS/KG41 - ( 14 / 13 )

15 Contact pressure, MPa Radial gap, m Normalized temperature, o C IFA-61.1 (BWR) 7 Thermal response to overpressure in IFA-61.1 (BWR) Calculated and measured normalized temperature (top) in TC-section against contact pressure and gap size (bottom) at the BOC4 in IFA-61.1 Complex interplay of effects of temporary irradiation with lowpower and overpressure and use of He at EOC3, followed by gap opening after overpressure-ramp to 25 bar Signs of secondary gap closure after power/overpressure dips and rises, presumably due to fuel fragment relocation into the open gap space and reduced cladding creep-out rate 35 3 Cycle 3 Cycle E-6 4.5E-6 4E-6 Cycle 3 Ar He Ar Cycle 4 Calculation assuming Ar as filler gas. Calculation assuming He as filler gas. Measurement. Running average for measurement Cycle 3 Cycle E-6 3E E E E-6 1E-6 5 5E Time, h Time, h /STARS/KG41 - ( 15 / 13 )

16 Normalized temperature, o C Rod pressure, MPa Gap size, micron Rod pressure, bar LHR, W/cm IFA-61.1 (BWR) Thermal response to overpressure in IFA-61.1 (BWR) 4 3 Assumed idealized BCs for simplified analysis of IFA-61.1 test Time, hours Calculated normalized temperature in IFA assuming hypothesized BCs and initial open-gap condition Calculation with idealized BC.: Normalized temperature. Daimetral gap size. LHR = 11 kw/m +5 bar Effects of primary- and secondary creep. Effect of irradiation at zero-overpressure that results in temperature decrease due to gap reduction P= 15 bar +25 bar +25 bar Rod peressure Coolant peressure Time, hours /STARS/KG41 - ( 16 / 13 )

17 Hoop stress, MPa IFA-61.1 (BWR) Mechanical behaviour of IFA-61.1 (BWR) Calculated cladding hoop stress as function of LHR in TCsection for three power-ramps in IFA-61.1 test Mitigation of the contact and gap opening (lift-off) would start, according to the calculation, at LHR of 8.6 kw/m and 1.4 kw/m at overpressure of 15 bar and 2 bar, respectively Calculated hoop stress at different overpressure levels: at 15 bar (2726 eph). at 2 bar (463 eph). at 25 bar (6577 eph) LHR, kw/m Calculated and measured cladding diameter during IFA-61.1 test Calculation shows good agreement with post-test measurement /STARS/KG41 - ( 17 / 13 )

18 Effective bonding fraction B, % Effective bonding fraction B, % Relative linear elongation rate, % per 1 MWd/kgU IFA-61.1 (BWR) Mechanical behaviour of IFA-61.1 (BWR) Measure and calculated elongation rate in IFA-61.1 test using two FALCON models Estimated SST elongation rate at P=2 bar: based on measurement. based on the model with bonding. based on the model without bonding Burn-up, MWd/kgU Estimated current values of apparent fractional bonding (CAFB) in IFA P=175 bar 8 P=2bar Burn-up, MWd/kgU Burn-up, MWd/kgU /STARS/KG41 - ( 18 / 13 )

19 Diametral gap size, microns Over-pressure, bar Temperature, o C IFA-61.1 (PWR) Thermal response to overpressure in IFA-61.1 (PWR) Centreline temperature in TC-section: Measured. Calculated assuming bonding. Calculated after assumed de-bonding. 2 bar 25 bar 3 bar 76 Measured response of temperature in TCsection to overpressure ramps in IFA-61.1 test vs. calculation using two FALCON models Calculation assuming breach of bonding predicts the geometric gap opening immediately after the overpressure ramp to 25 bar, followed by the fast increase due to primary- and secondary- creep of the cladding Consequently, the calculated normalized temperature grows up drastically, significantly exceeding the measurement Calculation Restart Time, hours Rod overpressure. Calculated gap-size assuming bonding. Calculated gap-size after assumed de-bonding. Calculation Restart Overpressure Gap size Time, hours /STARS/KG41 - ( 19 / 13 )

20 Normalized temperature, o C Over-pressure, bar IFA-61.5 (PWR) Thermal response to overpressure in IFA-61.5 (PWR) Calculated normalized temperature in TCsection against rod overpressure during IFA test, using two FALCON models Normalized temperature in TC-section: Measurement. Calculation w/o bonding. Calculation assuming bonding The model assuming pellet-cladding bonding shows good agreement with the measurement, whereas the model with the bonding switched off overpredicts temperature T c = 36 o C Overpressure T c = 316 o C Time, hours /STARS/KG41 - ( 2 / 13 )

21 Effective bonding fraction B, % Effective bonding fraction B, % Relative linear elongation rate, % per 1 MWd/kgU IFA-61.1 (PWR) Mechanical behaviour of IFA-61.1 (PWR) Measure and calculated elongation rate in IFA-61.1 test using two FALCON models P = 2 bar Estimated SST elongation rate: based on measurement. based on the model with bonding. based on the model without bonding Burn-up, MWd/kgU Estimated current values of apparent fractional bonding (CAFB) in IFA % P = 2 bar P = 15 bar ca. 88 % Burn-up, MWd/kgU Burn-up, MWd/kgU /STARS/KG41 - ( 21 / 13 ) 1