ML042260050

From kanterella
Jump to navigation Jump to search
Hearing 07/14/04 - Intervenor Exhibit F (Official Exhibit 30), Technical Paper, Use of CATHARE2 Reactor Calculations to Anticipate Research Needs - Recd 07/14/04
ML042260050
Person / Time
Site: Catawba  Duke Energy icon.png
Issue date: 05/25/2004
From: Bourdon S, Chatelard P, Grandjean C, Guillard V
Govt of France, Institute for Radiological Protection & Nuclear Safety (IRSN)
To:
Office of Nuclear Reactor Regulation
Byrdsong A T
References
50-413-0LA, 50-414-OLA, ASLBP 03-815-03-OLA, Catawba-Intervenor-F, RAS 8292
Download: ML042260050 (13)


Text

-OA e

~c~27 DOCKETED USNRC August 9, 2004 (11:45AM)

OFFICE OF SECRETARY RULEMAKINGS AND ADJUDICATIONS STAFF EXHIBIT F NUCLEAR REGULATORY COMMISSION Docket No. 5pLaI3/P4

-0Lk OffaIl Exh.No.

3° In the matter of Puke Grfawiu Staff IDENTIED 7/i i

Applicant RECEIVED _

Intervenor V

REJECTED_

Conf'g Of V' Contractor DATE Other_

Witne__

Reporter F4&,

i.A d

/Iz-ec,' o se-6CY-O- 0

SEGFSM Topical Meeting on LOCA Issues Argonne National Laboratory, May 25-26, 2004.

USE OF CATHARE2 REACTOR CALCULATIONS TO ANTICIPATE RESEARCH NEEDS V. Guillard', C. Grandjean, S. Bourdon, P. Chatelard Institut de Radioprotection et de Sfret6 Nuclaire IRSN/DPAM: BP3 - 13115 St Paul-Lez-Durance, France valia.guillardO~irsn.fr ABSTRACT To analyze the consequences of the introduction, in Nuclear Power Plants, of advanced fuels at high bum-up, decided by most of the utilities in western countries in order to reduce the fuel cycle costs, IRSN has initiated a research program focused on the study of such PWR fuel behavior in LOCA conditions.

A first step of this program, comprising analytical and experimental parts, has been to identify the main physical phenomena, linked with thermomechanical behavior of irradiated rods in bundle geometry, to be taken into account in reactor safety analysis: cladding deformation and flow section restriction in bundle geometry, mechanical interaction between neighbor rods or structures, axial extension of balloons ; cladding oxidation and secondary hydriding ; fuel fragmentation and relocation, balloons filling rate and FP release ; fuel rods thermal behavior in bundle geometry during reflooding conditions, rewetting of the claddings around ballooned regions with fuel relocation mechanical resistance of irradiated claddings in post-quench conditions.

This paper summarizes an analysis of sensitivity calculations performed with CATHARE2 "Best-Estimate" code, used in France in the frame of realistic methodology to evaluate safety margins.

The objective of these calculations is to point out, among parameters affecting last-mentioned phenomena, those for which taking into account basic uncertainties lead to important uncertainty on global code response (Peak Cladding Temperature, oxidation rate...). That is the case of fuel relocation phenomena, whose impact is highly dependent on parameters such as, in the example of LB LOCA transient, cladding radial and axial deformations in bundle geometry, burst criteria, balloon filling rate, thermalhydraulics around balloons. A lack of knowledge on theses parameters for irradiated U02 and particularly MOX fuel may lead to reduce safety margins.

This study may provide some elements to identify future research needs to complement present experimental data base, reduce uncertainties and develop more realistic calculation models, which may better fit the thermomechanical behavior of advanced irradiated fuels.

' Corresponding author

?SEGS TPCLETNONLA ISSUE MEET 4"'..'

AroneNtonlLboaoyMy25-27,2004 USE OF CATHARE2 REACTOR CALCULATIONS TO ANTICIPATE RESEARCH NEEDS V. GUILLARD, C. GRANDJEAN, S. BOURDON, P. CHATELARD Institute for Radiological Protection and Nuclear Safety (IRSN)

Major Accident Prevention Division I Fuel in Accident Situations Department BP3 - 13115 St-Paul-Lez-Durance Cedex - France Content of the presentation

> Introduction

> Main physical phenomena to be modeled

> Main hypotheses of the calculations

> Main results of the LB LOCA calculations

> Conclusion and perspectives fia0ITTHMP ANLOn May 2527, 2004

^

vala~qullar@.irn~frANL, May 25-27, 2004

111 U,;,

tQ General Background

~ oemanagemnjotJ.-1 StudyV'f-PWR Jqs eatr stuie Kfue Thterals th~m omechnicI bh'avior withCATHAE M CATHARE2

> French thermalhydraulics system code with CATHACOMB fuel module 0

Objectives

,s..

  • ..j..

L..f I_

r_

J iuerXITy suture researcn needs Tor new generation of tuels

> Improve knowledge, models and calculation methodologies 1i,

OApproach used to reach this objective

> Identify physical phenomena involved in thermo-mechanical behavior of advanced irradiated rods in bundle geometry (State Of the Art by C. Grandjean & G. Hache)

> Take them into account in the modeling

> Quantify basic uncertainties (CIRCE tool and sensitivity calculations)

> Evaluate global uncertainty on CATHARE2 code response (SUNSET tool)

E usdiboundary conditionsfor further more precise lu ations and analysis: under LOCA conditions V

yal.a auillard(cirsn fr ANL, May 25-27, 2004 I.

MAI PHY;SIC-AL PHENOM ENA TOEBE MODELED;`

~; liked-ith' the use ofavanted fesa ihbri-p Reference C. GRANDJEAN &GIACHE "LOCA Issues Related to Ballooning, Relocation, Flow Blkagte and Coblabiity.:

Main Findings from a Review of Past p

tlgrn Temperatures ('C) ifit-Fuel i

'-Clad TYPICAL LB LOCA TRANSIENT EASY NOTSO EASY Ox datin kiec Fuel f ragmentation and....

-seonar hydrid' 4ng<

(debris size,10 granulometry, porosity) 800 Clad ballooning Mechanical

: interactions Ar0 between neighboT"'

rods or structures Axial extension of balloon

.,Fi Flowrestrion dL (instant of -occurrence, balloons fillin ti) 50 100 150 Time (s)

Clad coolability around ballooned regions with fuel relocation jel rods thermal behavior iring ref looding conditions Post-quench

  • .V residual ductility Fuel Clad Bundle Haua rrES Meet. oi LOcA issues Vcatiaqau a dis~rANL, May 25-27, 2004 o.-

I,.CAThAE2CLCLAINS Man ypothess X l Calculation of several transients including Large & Intermediate breaks LOCA i l Use of standard CATHARE2 versions (V1.3L, V2.5) o Basic uncertainties on CATHARE2 models taken into account

> Consistent with PIRT implications El Zry-4 cladding without hydrogen uptake effect on mechanical properties Q Hypothesis linked with the use of irradiated fuel at high burn-up

> Modification of thermal properties conductivity and thermal capacity laws from SCANAIR code

> Deletion of protective effect of initial oxide layer on transient oxidation
outcome from experimental ANL and TAGCIR program analysis on irradiated Zry

> Introduction of Baker-Just correlation to calculate oxidation rate O Initial state of irradiated rods given by METEOR code

> Gap width and pressure, radial and axial fuel power profile, external cladding oxidation profile, cladding thickness...

.kI2-Exam)le:

I ofthe ct of irradiated fuel relocation ditcndiions on a PWR9OOMwe U.e Of calcuiations GFSM han.e.s va~i~pui~ar~ir~~frANL, May 25-27, 2004 to

CATHR2LL

~AcLT ON S Siml b

4afeigcatitn!O ph enomena U Main hypothesis Balloon filling rate Rel o mh by relocalized fragments

  • Relocation mesh g

m90

= burst location mesh filled by fuel debris 80 U

from upper meshes

  • Filling ratio of the balloon

= function (balloon size) 60 A

A Balloon size derived from 5 0 CATHACOMB calculation of o

mean elongation 040 0 PBFILOC-gammascanning U PBFILOC-micrographies Basic calculation value = 61.5 %

  • Modification of the power factor E 30 A F12 and the fuel mass in the X

-Upper bound relocation mesh 20 0 EdF calculation K IRSN reference calculation

  • Considering relocated fuel as 10

+ IRSN sensitivity calculation a porous medium (Imura-Yagi model for the conductivity calculation) 0 20 30 40 50 60 70 80

  • Flow blockages: impact on T/H Ballooning (%)

and fuel cooling not taken into account t

SEGES

pMt, Asues

-uillardO<irsnfr ANL, May 25-27, 2004 r.

CATHARE LB A C ALCLTION L I High sensitivity of the relocation impact on PCT

> Burst at 24 s on stress criterion deduced from EDGAR experiments (£ at rupture = 57.5%)

LB LOCA - Irradiated U02 Hot rodcklddihg feinpEreturC orburtjpnecsh 1200.0 1100.0 1000.0 5'

a.

E

,O-.

900.0 800.0 700.0 600.0 Yl rel~o~ctation impact depends on:

bUrstttime, which determines the energy

~ re-di6tributionJin the burst mesh o

ifilling

ratio, which d~terrines relocated fuel mass 3:O to 150 C on PCT
z dbaoonfilling. ratio 61.5 to 70%)

500.0 400.0 300.0 200.0 100.0 flme Ls]

WaiaaIlardc~irsn. r ANL, May 25-27, 2004 I.

I,

4 t,

XO!'S'"S gs P,4v'29' itX n

M;'^.a

,1E X_

{Ss, ;,

At '." t...

s _

<7.,/",

1,34 2

1 i !,,,,;

i V i

!' t:. S.;:,,:

St.;

':"'::t:-''^"'"'.'S 1

!.:,,:.:.;: X 1

f".'

1 Fuel relocation impact on ECR and residual P-layer

> Burst at 24 s on stress criterion deduced from EDGAR experiments (£ at rupture = 57.5%)

LB LOCA - Irradiated U02 LB LOCA - Irradiated U02 Hfot rod oxildatilon rate oab ursrron esh Hor rod betia-&^er tickhuer at burrimei 20.0 8.0 5.0 4.0 12.0 600.0 550.0 500.0

  • r 10.0 0a UJ.

5.0 6.0 4.0 2.

0.0 o.o U,2 PE p

I I

50.0 50.0 100.0 1200 Time [s]

rime [as

+

t 7% on ECR (depending on balloon filling ratio: 61.5 to 70%)

5o <o.toon residual P-layer thickness j<.

(dpenin i ba~loon filling ratio: 61.5 to 70%)

R aivaha.qui1ardcrsnr:

ANL, May 25-27, 2004

'<CAHARE2 BQA$CLTQSlt EC Fuel relocation impact on PCT Burst at 18 s on stress criterion deduced from EDGAR experiments (s at rupture = 58.7 %)

3t,-i, 3.

'3

.K

-:E:!

'p 3

3e,+

\\: '

. ';(2

,4, LB LOCA - Irradiated MOX Hot rod aedding rempenietu e tit butesh n

I I

I v

T I

I I

l 43 CLAD BURST

~>

FUEL RELOCATION Withtt rebcatan l

With rebcation n 4

  • e9Enh'ance of fuel relocation impact k ~ii i~.
M4 Ms3 W

4^

3 i

I.

.. I.I

°o

+I6O 0,.C on PCT 1

I Time [a]

vat~~auitar~irs~frANL,

.M ay 25-27, 2004

7 CATHRE2 LBLOCA CALCU LATIONS

s ~~~

~

esiivt to fXe tyes MO (2/2)*

$S X 0 ':^

El Fuel relocation impact on ECR and residual P-layer

> Burst at 18 s on stress criterion deduced from EDGAR experiments (Eat rupture 58.7 %

LB LOCA - Irradiated MOX LB LOCA - Irradiated MOX Hot rod oxiJatioh n-te atburrtrnEs/

Hot rod beta-layer residual thickness at burst inesh

4.

I

,.,.,.,.,.I D) t

wCZ FE i,-,"BUT;'

Time [a]

Hige riimpact tonEM-AfidFresidual,B-layer Time Is]

iL ia;rg uncertainty ion: hydrgen uptake effect nd mecanical properties depending on clad alloy

-'.Need to reduce crtain ties u s e ~ o, 6 6e d

E F M T p M e, a L c i ~ e yalta.guiI ard~irsin.fr ANL, May 25-2 7, 2004 I.

4

-~ r

. I bnt iL.

1,,.

U El, Use of CATHARE2 code for high burn-up fuel analysis under LOCA conditions

> Example of Large Break transient

> Emphasis on fuel relocation phenomena impact

> Main results: PCT and ECR increased High uncertainty on global code response due to identified lack of knowledge

> Instant of fuel movement ^-

HALD

?

  • imposed as burst time in the simulation
  • depends on clad ballooning/deformation and burst criterion

> Balloon size HALDE:'

i s..~-,s,

  • which is also linked with ballooning/deformation model

> Filling ratio FR2, PBF, HALDEN?

> Relocated fuel properties

  • fragments size, granulometry, porosity, conductivity, ANL, HALING

> Bundle effects Need of inte tes Axial extension of balloon Flow blockages Clad coolability around ballooned regions with fuel relocation e

CUeo ir'R t J

n neec S

LAi' U' ~

J varia.guill rd~isa f

ANL, -May 25-27, 2004

i1

'KG~~ONPLU IONANDPRPCIE

'A.

t U Perspectives S Modification of clad mechanical properties to take into account hydrogen Gi

.; N uptake effect

> Study to be complemented by calculations using CATHARE / FRETA E rods 3D thermo-mechanics

. rod-to-rod interactions models

  • cooling and reflooding models for overall bundle

> Use of NEPTUNE 3D local module for flow blockage cooling calculations

> IRSN plans to develop a new code for fuel LOCA calculations Stn.

lache) + Analytical studies

~

lentfiatin f'kowedg> mprovement on iradited'rod ehavor uderLOCA cond~itions on, S

> evlp ewgneration offel a 'nd blundle moddels Reduce ucin new, margins

-j-p

---Mrraa'bb '

ag

,~ se~fCAHAF~2 eacor a~c~at~nso~n~ct at esarc neds~

SEG FSMUW Top.. Meet on J'LQ "s~e

~

I~J Vaia.~it~rd(~irs~frANL, May 25-27, 2004