ML19212A526

From kanterella
Jump to navigation Jump to search
Research Info Ltr 25:transmits FRAP-S3 Code Description & Verification Documentation.Forwards Id Natl Engineering Lab Repts TFBP-TR-164,Revision 2 & TFBP-TR-228
ML19212A526
Person / Time
Issue date: 03/24/1978
From: Levine S
NRC OFFICE OF NUCLEAR REGULATORY RESEARCH (RES)
To: Case E, Minogue R
Office of Nuclear Reactor Regulation, NRC OFFICE OF STANDARDS DEVELOPMENT
Shared Package
ML19211A032 List:
References
RIL-025, RIL-25, NUDOCS 7912140454
Download: ML19212A526 (17)


Text

Nw v L-

.y y

'o UNITED STATES 8

NUCLEAR REGULATORY COMMISSION 7~ p r l'37

~~

E WASHINGTON, D. C. 20555 k..... p n.

.,i c - ;

Mg 21 197g Jin &d;

. M..

B" MEMORANDUM FOR:

Edson G. Case, Acting Director Office of Nuclear Reactor Regulation Robert B. Minogue, Director Office of Standards Development FROM:

Saul Levine, Director Office of Nuclear Regulatory Research

SUBJECT:

RESEARCH INFORMATION LETTER - #25 - FRAP-S3

REFERENCES:

1.

P. E. MacDonald, et.al., "MATPRO-Version 09: A Handbook of Materials Properties for Use in the Analysis of Light Water Reactor Fuel Rod Behavior,"

TREE-NUREG-1005, December 1976 2.

J. Rest, " GRASS-SST:

A Comprehensive, Mechanistic Model for the Prediction of Fission-Gas Behavior in UO2 Base Fuels During Steady-State and Transient Conditions," to be published This Research Information Letter transmits the FRAP-S3 code description and verification documentation.

Introduction FRAP-53 is a best-estimate computer code that calculates the thermal and mechanical response characteristics of a nuclear fuel rod operating under steady-state power conditions.

It is the third version of a code developed to provide accurate initial values of fuel-rod parameters for input into transient analysis codes such as FRAP-T and RELAP.

It is capable of supplying the hot-state values of such quantities as:

1.

stored energy 2.

radial temperature distributions at given axial locations 3.

total fission gas release 4.

rod internal gas pressure and composition 5.

clad deformation 6.

amount of pellet-clad interaction (PCI) 7.

fuel deformation (swelling, densification, relocation, and thermal expansion) 8.

fuel-clad gap size and gap conductance 9.

clad corrosion and hydriding.

7821903o4 H

1 4'5//

1568 294 79 m o u

., 9 ' '.

n c-Edson G. Case 2

Robert B. Minogue All of these quantities are strongly dependent upon the operating history of the rod, and each will have a large effect on the predicted and measured response of a fuel rod during a transient. The code, therefore, has been designed to provide these and other quantities for any given power history as initial conditions to the transient codes.

The required material properties such as thermal conductivity, thermal expansion, etc., are obtained from the MATPRO package (reference 1).

The verification of the FRAP-S3 code had two major objectives:

(1) to determine the code performance in predicting the available, qualified, experimental data, and (2) to identify those areas that require more sophisticated modeling or more experimental data.

For the first time, the code performance and data were analyzed using statistical methods, since sufficient data are now available for significant results. Thus, all of the major response variables are presented alocg with their corresponding standard error bounds.

A stat %tical spread can now be put on such input as stored energy when coapting the behavior of a rod during a transient. The verification procer.ure used information from over 700 fuel rods containing a wide rangr. of operating and design parameters.

Discussion Code Description.

A typical PWR or BWR fuel rod is divided into a maximum of 15 axial segments, each operating at a power level averaged over its length.

It is also divided into a maximum of ten radial rings.

Each ring-segment volume element is assumed to possess averaged properties such as temperature and power. The fuel rod power history is approxi-mated by a series of steady-state power levels with instantaneous jumps from one level to another.

Five major calculational models are used in the code:

thermal, gas pressure, fuel deformation, clad deformation, and fuel-clad interaction. Although these models are interdependent, the calculational techniques used in them are distinct, and can be described separately.

The fuel rod temperature model is based on the following assumptions:

(1) there is no axial or circumferential heat conduction, (2) steady-state boundary conditions exist during each power step, and (3) gamma heating effects are negligible.

The calculation begins with the deter-mination of the bulk coolant temperature at each axial segment using the standard thermal-hydraulic equation for channel flow, and specifying the needed thermal-hydraulic parameters.

The temperature drops between the coolant and rod surface, across the clad oxide layer, and across the cladding are then computed in sequence using the Fourier heat flow law as applied to thin layers.

Next, the temperature drop across the 1568 295

Edson G. Case 3

Robert B. Minogue fuel-clad gap is computed. This calculation requires the use of iterative methods since the gap conductance depends strongly on tem-perature through its two major parameters, gap size and gas thermal conductivity. Each of these parameters depends, in turn, on the other major models of the code through such processes as fission gas release, fuel swelling, fuel densification, fuel relocation, and fuel and clad thermal expansion.

As a result, the gap temperature calculation and the subsequent fuel temperature calculation is done within an iteration loop which encloses the calculations of the above phenomena.

The fuel pellet temperature distribution is computed using the /KdT method (see ) in which account is taken of the effects of neutron flux depression on the volumetric heat generation rate.

The fuel rod internal gas pressure model assumes:

(1) the perfect gas law inside the rod, (2) gas pressure is constant throughout the rod, and (3) the gas in fuel pellet cracks is at the fuel average temperature.

The rod is considered to consist of regions along its axis, each having an average temperature and gas volume determined by the average fuel and clad dimensions over the axial segment. The plenum gas temperature is calculated via convective heat transfer between the plenum gas and its contacting surfaces; i.e., the top of the pellet stack, the holddown spring, and the cladding. Ganma heating is accounted for in the plenum spring.

Fission gas production and release are computed using empirical correlation models. Gas production is burnup-dependent only, whereas its release utilizes a correlation equation containing time, temperature, and fuel density.

The fuel rod deformation model contains three submodels:

thermal expan-sion, irradiation-induced swelling, and densification.

It, therefore, assumes that no mechanical deformation of the fuel occurs via fuel-clad contact or thermal stresses, and that no fuel creep occurs under applied stress. The thermal expansion model assumes that each fuel ring expands fully in both the axial and radial directions with no constraining influence arising from the circumferential expansion.

The fuel swelling model is obtained from MATPRO via a correlation equation dependent on burnup and temperature.

The fuel densification calculation is also obtained from MATPR0 correlation equations, which relate the density change to the sintering temperature, initial fuel density, and burnup.

The cladding deformation model consists of many submodels which account for time-independent elastic and plastic deformation, and time-dependent plastic deformation (creep).

It uses the incremental theory of plasticity and the Prandtl-Reuss flow rule (see Enclosure 1) to determine the plastic strain increments for each set of load conditions; it obtains 1568 296

t 4

Edson G. Case Robert B. Minogue The model the required mechanical properties from MATPRO correlations.

(1) work-hardening is isotropic; contains the following assumptions:

(2) stress, strain, and temperature are uniform through the cladding thickness; (3) no slippage occurs at the fuel-clad interface; (4) bend-ing stresses and strains in the clad are negligible; and (5) the loading In general, the methods used to solve and deformation are axisymetric.

for the stress and strain in the clad depend upon the structural relation-If the ships between the fuel and clad at the time of the load step.

fuel-clad gap is open, the stresses are calculated directly from the stress theory of a thin cylindrical shell with specified internal and If the gap is closed, the model used is a thin external pressures.

cylindrical shell with a prescribed external pressure and a prescribed radial displacement of its inside surface. This radial displacement is obtained from the fuel deformation model. Moreover, since no slip is allowed, the axial fuel strain is transmitted directly to the cladding, thereby prescribing the axial strain in the clad.

If the gap is par-tially closed (i.e., pinched off) along the rod axis, a " trapped stack" regime results. The model used in this situation is a thin cylindrical shell with specified internal and external pressures and a prescribed total change in length of the cylinder. Whereas, the first two models solve for the stresses and strains at one axial segment at a time, the trapped stack model simultaneously solves for all the axial segments Details of the iteration loops and computational within the regime.

methods used are given in Enclosure 1.

The fuel clad interaction model was developed to predict cladding failure as a result of any of the following phenomena:

1.

clad melting 2.

oxide layer wall thinning 3.

ballooning (without rupture) 4.

eutectic melting 5.

clad collapse 6.

overstress or overstrain failure 7.

cumulative stress damage.

Actually, the model consists of several subroutines, each of which com-The putes the probability of failure via one of the above mechanisms.

model then combines these probabilities into a single overall failure probability of the cladding.

Results - Verification of FRAP-S3 As stated earlier, the verification of the FRAP-53 code utilized data from over 700 fuel rods, enconpassing a wide range of operating and

,ren OGJ

}.j o O

Edson G. Case 5

Robert B. Minogue design parameters.

The verification effort involved over 700 individual computer runs and extensive statistical analysis of the results.

As shown below, the overall performance of FRAP-S3 in predicting the important parameters of the models described above is excellent.

Although improvement may be desired in some areas, in other areas the code approaches the limit of accuracy of the experimental data. The code performance in predicting rod temperature, rod pressure, and rod defor-mation is summarized below.

The thermal performance was analyzed using fuel centerline temperature results of over 100 rods representing over 800 data points.

Figures 1 and 2 sumarize the results for unpressurized and pressurized rods, respectively. The standard errors between predicted and measured values based on these figures were 198 K and 254 K, respectively. On a per-centage basis for a centerline temperature of 1500 K, for example, these results translate to 14% and 17%, respectively. This indicates an accuracy within, or close to, data measurement uncertainty.

The fuel rod pressure performance is sumarized in Figure 3 for 48 rods (including both pressurized and unpressurized) and 658 data points.

The respective standard errors in the calculated pressure for pressurized and unpressurized rods was 1.34 and 0.66 MPa. The group of underpre-dictions between pressures of 7.58 and 11.72 MPa corresponds to BOL measurements on two rods which exhibited significant transducer drift.

Burnup values ranged from 3000 to 22,000 mwd /MTM.

The errors in terms of percentage correspond to -14% for a pressurized rod operating at 10 MPa, and ~34% for an unpressurized rod operating at 2 MPa.

The latter error is reduced to less than 10% if only BOL conditions are used. This reflects the strong influence of gas release on the pressure uncertainty in unpressurized rods.

The fission gas release comparison is shown in Figure 4 for 176 fuel rods.

The standard error in the gas release fraction is 0.188 which corresponds to about a 40% error in a rod releasing 50% of its gas.

As can be seen from Figure 4, the large overall error is caused by considerable averprediction at the low end (less than 10% release) of the data spectrum.

The code performance is much better for high release data.

The poor performance at the low end is believed to be caused by inadequate modeling of the diffusion processes occurring at the start of release.

It is hoped that the GRASS computer model (reference 2) for gas release, when coupled to the code (either as a fast-running approximate version or as a GRASS-derived correlation equation), will reduce the error considerably.

1568 298

Edson G. Case 6

Robert B. Minogue One of the most important (and most difficult to predict) indicators of the rod deformation model, is the prediction of the heat rating at which gap closure occurs. The performance of FRAP-S3 in this respect is shown in Figure 5 for 77 rods. The horizontal lines indicate the range of data uncertainty.

The standard error is 13.4 KW/M, which corresponds to an error of ~50% for a typical PWR rod operating at 26 KW/M (8 KW/ft).

Although this uncertainty is not small, it is a considerable improvement over previous versions of the code.

The improvement noted here (and also reflected in the thermal performance) is due to a new fuel relocation model coupled with an improved accounting of the effects of relocation-induced cracks on the pellet thermal conductivity.

Another measure of the performance of the rod deformation model is the total permanent axial expansion of the fuel stack caused by swelling and densification.

Figure 6 shows the code performance in this area for 100 rods and 368 data points. The burnup range in the data sample has been significantly extended over previous versions by inclusion of power reactor post-irradiation examination results. The calculated standard error is 0.44%

of the active length. The amount of deformation, though small, is gen-erally underpredicted, since compression effects on the fuel have not been modeled as yet.

Finally, the code performance in predicting the permanent cladding hoop and axial strains is illustrated in Figures 7 and 8.

A total of 170 rods containing 393 data points were used in the hoop strain comparison, whereas 115 rods with 161 data points were used for the axial comparison.

The computed stadard errors were 0.59% of cladding diameter and 0.47%

of the active length, respectively.

The hoop strain calculations are dominated by the creep collapse model in the code and tend to overpre-dict the data.

The reasons for this may be twofold:

the add-on fast-flux term for creep may be too high, or the free-standing clad model used neglects the possible support that the clad can receive from the fuel.

The points showing large overpredictions of the axial strain in Figure 8 correspond to small gap or high temperature conditions, under which structural gap closure is calculated.

Since no stress-induced fuel deformation is modeled, the cladding strain consequences of gap closure are overestimated.

In any event, the measured and predicted strains are quite small and are not expected to have much influence on the effective gap size and corresponding thermal calculations.

It has been shown above that the thermal model gives excellent agreement with the data, and that the onset of gap closure is reasonably well predicted; both of these predictions take into account the computed cladding strains.

Table I sumarizes the standard errors for the above response character-istics, plus the cladding corrosion and hydriding behavior and axial fuel thermal expansion results.

1568 299

9 Edson G. Case 7

Robert 8. Minogue Recommendations The above results are offered for user office consideration for applica-tion to the identified regulatory need.

The information presented herein should be especially useful in the current RSR/NRR Code Coor-dination effort. The statistical analysis used will aid in the licensing evaluation of the code output, and serve to direct the effort in the development of the evaluation models to be used in the new steady-state FRAPCON code.

For information on further evaluation of the results, G. P. Marino of the Fuel Behavior Research Brarah of Reactor Safety Research may be contacted.

')

Sau Levine, or Office of Nuclear Regulatory Research

Enclosures:

- J. A.

Dearien,

G. A. Berna, M. P. Bohn, J. D. Kerrigan and D. R. Coleman, "FRAP-S3: A Computer Code for the Steady-State Analysis of 0xide Fuel Rods, Volume 1 - FRAP-S3 Analytical Models and Input Manual," October 1977, TFBP-TR-164 - D. R. Coleman, E. T. Laats and N. R. Scofield, "FRAP A Computer Code for Steady-State Analysis of 0xide Fuel Rods, Volume 2 -

Model Verification Report," October 1977, TFBP-TR-228 See cc next page 1568 300 i

cc w/o

Enclosures:

L. Rubenstein, NRR V. Stello, 00R D. Ross, DSS R. Baer, DDR P. Check, DSS Z. Rosztoczy, DSS C. Fitzgerald, MIPC W. Reinmuth, IE R. Mattson, DSS M. Kenemuyi, A/DGES G. H. Smith, IR C. MacDonald, FCTR J. Volglewede, DSS R. Lobel, RS 1568 301

SYM80L RUN REFERENCE O

14 HPR40W O

33-36 IF A.226III 6

4345 HPR 107(8)

X 53,54 RISO-269,RISO-415(8

  • O 57.59,60,612,613 IF A 223('i 88.97)

M 6244 WC AP-2923(12.99) 6549 IF A 130,13103I 4000.0-U 219 221 AE 318(22) 319,320 IF A 116(36,39)

H 322-324 IF A 117(37 39) 370,371 IF A 116(36) 424 426.428,429 IFA 431(56-591 b55 EPRIM 3500.0 -

IF A 208(85) 587 C>

588,589

! F A-211(85,95)

ED 590A,5908 IF A 224(85) 591 IF A-230(85)

Q 592-594 IF A-410(85)

_ 3000.0 - 3CiC 595 I F A-411(85.86)

M 25 611 IF A 207(96)

/

E

+

614A-6158 THECO 1,11(89) d, $

/P=M 8

C 623.624 IFA 106(92)

M M

7 629,630 I F A-107(92)

/

$ 2500.0 -

0 586 iF A-206(85)

/ A 7

b

+ EP8 2000.0 -

6

@a-

- fk O

8 m? {

w 1500.0-o 2

O W

w h

=

r-9 de

@ 1000.0 -

E

=

500.0-4 o

/

/

/

/

0,0 0.0 500.0 1000,0 1500.0 2000.0 2500.0 3000.0 3500.0 4000.0 MEASURED FUEL CENTERLINE TEMPERATURE (K)

Figure 1. FR AP-S3 Predicted Versus Measured Center Tempaatures for Unpressurized Rods.

1568 302

4000.C q SYMBOL RUN REFERENCE 273.274,271,278.433,434 PBF(25 27.29,30.61.62.64) 437,438,446 448 276,2M,4R439,454,455, P8 F ( 28.31.62.66.68) 3500.0 G

460A-463A E

415.416 IF A-429(55) 441-445.553 IF A 418,419(63.83) 4 450,453 P8F(65)

A 596,597 IF A-418(85)

E 3000.0 -

E

/

5 7

P=M

/

/

c

$ 2500.0 -

/

d:

4 N

l

./e<

W 5

E 2000.0 -

g j

E M f *.f.

  • 1 1500.0 -

@ 1000.0 -

/s v$

500.0 -

/

/

/

/

/

0.0 4' 0.0 500.0 1000.0 1500.0 2000.0 2500.0 3000.0 3500.0 4000.0 MEASURED FUEL CENTERLINE TEMPERATURE (K)

Figure 2. FRAP-S3 Predicted Versus Measured Center Temperatures for Pressurized Rods.

1568 303

SYMBOL g

REFERENCE O

33.34.50 52 IF A 22687I

+

70 IF A.132813) 4 273.274.27s.433.437-4 3.

18.0-PBFt 2s 21.29.30.61.62.64) 446448 e

279.280.463A PBF128.31.68)

Q 314.315 MP(pie)(33 34) 319.320 I F A-116(36.39) 16.0-(

321 324 IF A-117(3 7.39)

C 331-334 IF A 178(40) o 336-339.386 I F A-181(41)

/

[/

P=M 7

U 372.373 IF A-142(40)

O 387-389 I F A-225( 48)

$ 14.0 E

41i-423 IF A-429(6D)

[

]

/

8

/

s

/

w 12.0-

/

  • * *e s C
g. # +/

I

  • /,ff**,

5 E

10.0-

  • /

4

/

z

/,*

a O

8.0 -

A,p /

-t O

e

.v. jf

+

le.0

  • 8.*

s g

Y

/.s 40-4

/

as

$e

/

8

,46 /

o v4 g 9 40/

2.0 - +

%/

/

0.0,

0.0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 16.0 18.0 20.0 MEASURED ROD INTERNAL PRESSURE (MPa)

Figure 3. FR AP-S3 Predicted Versus Measured Rod internal Pressure.

1568 304

102 f

5 OO O

O 0

e

  1. g o*/ u

= = 8 W'

+

+

w ^ ^ t.* 5 9 ',,

+ e++e, e49 6

e

+

e

/

0

+ +

o 1

+

+

+

/

10 -

+*

/*+ +Eg a

e ao

/

1m 8 e

.mB ge +

h O

m

/

w M

e a

+

E SYMBOL RUN W

g REFERENCE

_J

/

O 73 77.79 82.84.85.87 e & wi

m E

E

/

O 95A 968.99100.103107 SA XTON il '"

l C

10 a

f a

iO8.109.iiu t2 ii4.iis sA xrON uv

M 120.122 4

124.126.1245 130.133 136, i i,,

C 4

g,g 2

/

138 146 1201 X

205 208.211 AECL 1676 O

20 M

/ g D

213 21*

AECL 1685 7

233 239 AECL 2230<231

/

12 as T

251 253.256 258.264 265 GE AP 4264 1

O

/

206 295.297.298 PRTR132)

W

/

D 299.302 313.315 HPtPIE d33 34' 4

316 318 IF A 208' J

g y

319.320.321 324.325 328 IF A 116,117.HP(Pp36 3 7 34 39) 1401

/

331 334 IF A 178 w

N 336 339 IF A 18118 "

g 10-g -

3CZ 4i4 err < s3 548 o.

/

P=M s

496 499 502 507.509 513 519 MAINE YANKE E '"

I O

$43 545 BIG ROCK POINTas t a2i 3

650 EPRit e3) j

(>

555 EPRit 8 3) 557.559 561 563.566 568.

t a s)

YBWR. ORESDE N

/

570.578 588.589 IF A 211 s$ 95) t 1921 625.627 IF A 106 t921

/

9 631.632 IF A 107 A

633639 IF A 138I93 3933 a

640644 lF A 142

/

T 645.646 IF A 206'93I 10-2 10-2 10-1 100 101 102 MEASURED FISSION GAS RELEASE (%)

Figure 4. FR AP-S3 Predicted Versus Measured Fission Gas Release Fraction.

1568 305

SYM8OL Rg REFERENCE E

29,30 AECL4072(5)

O 31,32 IF A-225(6) 70.0,

G 402404 IFA4041(50)

O 405 IF A414(51.52)

+

408 4 10 IFA404lII50) 65.0 -

Q 616A,6168,6178,618A,6188, IF 227LllJil(90.91) 619 A,6198 0

3549 IF A.227(7) 60.0-

.d 40 HPR-107(8) a 4345 HPR-107 s) t X

58,60,61 I F A-223('3 E

71 IF A-132( 13)

^

E 55.0-p 148 HPR-13209) h g

273,274,277,437,433,434, P8 F(25-2 7,29,30.61,0 4 446449 0

50.0-276,460 A,279,281 P8 F (28,31.68)

E 317 IF A-208(35) k 353 IF A-230(35.42)

Q~

45.0-D 386 IF A.181(41)

Q 405,406 IF A414(51,52)

W G

676486 IF A402(96) 1 40.0-A 687491 IF A-215( 100) 692493 IF A-216( 101) 3 T

694497 IF A413(38)

XX 345,347 IFA 150(ddl o

35.0-u

/

m c

s 7

30.0 -

7

_e

/.

T-5 I I 15 d t~

25.0

. ~_

c-

=

- = _

Q f_L " ^

o h

20.0-gv 4 Q

___ _ _l,_

d 5

g s

2 w

/% ___. o

_ f

=

g 15.0 -

g 3 -

sg

.s g _ s__.

e tj m

10.0 -

/ -

/

5.0 -

/

/ P=M

'b 0.0 O.0 5O 15.0 1$.0 2d.0 25.0 30.0 35.0 40.0 4$.0 5d.0 55.0 60.'O 65.0 MEASURED INITIAL GAP CLOSURE HEAT RATING (KW/m)

Figure 5. FR AP-S3 Predicted Versus Measured Gap Closure Heat Rating 1568 306

SYM8OL RuN REFERENCE A

88,89 KWO' 15) idSI 0

359-362 WEST O

363 366 B & W"

391401 IF A4011 tid 9) 7.0 -

405 406 IFA414f 51.52) 9 595,594,599 IF A 419(85.86)

B 600 602 IF A40911(86)

Y 603405

  1. # #'6

Y 606 4'O IF A 4 f 9,lF A 42 7f 87) 0.5 - S 340 345

" %2

/

494499.501503,50s. Soy' MAINE YANKEEfif) 9 9 5 f.5 3,516,517.519

/

X 5 h M 8. RO81NSON(72 ao; l

G 543-548 BtG ROCK POINTl81.82)

/

0.0 - O 621626,628 IF A 106f92)

/

u ~

fa V

647456 8RP(MO )'848 p

2 g

g y

-0.5 -

  1. c
  1. o U

aUS

.a h

cp

+

j

+ $

/

+

  • - 1.0 -

/,

3

+

E

/

y

./

0 E

883 p

4 a,

/c g

-1.5 -

0

/

w

/

e 8

al/

+

su 5

-2.0 -

Yo

[ /g o

/

/

- 2.5 -

/

P=M p

/

/

s'

-3.0

-3.5

_j0

-2'S

- 2'.0

-3.5

-1.0

_o'.5 s

s s

0.0 0.5 3 '0 MEASURED PERMANENT FUEL AXI AL DEFORMATION (%)

Figure 6. FR AP-S3 Predicted Versus Measured Permanent Fuel Axial Deformation 1568 307

r 2.0 -

1.5 -

6 E<

/

x 1.0 -

/

P=M e

o D

/

/

oo

/

/

O o

o 9

++

+

=

H g

+

0.0 -

v 4

ES I 2

f v vv T

/

C)

W A

g

$vMPot As iN Rf F E HE N( E E I,E E

7 _i u2,.,

sa xiou n'".

g o g o

aE, n

ir ii4iisi2c122 sa i oN i. ' '*

  • p

.gg g g

+

1'4 125 127,128 137.143.52 MP3 132' '*

O 82t+

243 217 AEll 1685 f

Q f

23 X

2'3 239 A E CL 2230 W

C AE O

2.i 253 2ss 20a 2c4 c.s av 4264 3'

- 1.0 -

/

8 2

  • 2' '

6 317 348 IF A 208' 3

/

E A t 1414 H T R '53 'ed '

g

0 466 4t>8 4d'J SA A I UN' 'O' g

g, 4

/

494 4% 4'19 502 bul.509 bli.

M A l' y

6t 213 515.516.5i1.519

/

8 ' d *'

M G

3 A3 Sea big HOLK Pf>iN T

e um2 > cam >( 1 >e9

-1.5 -

/

EB e

v ev.H. ostsot s **

513.513.514

~

64i cs6 sa emo i

/

a z

a sa ca awm AA

+

% 1 552 555 E phi'8 3' e

s2a. 4-M e HoeissoN.iiam a

y

-2.0 /

-2.0

-1.5

- 1.0

-0.5 0.0 0.5 10 1.5 2.0 ME ASURED PERMANENT CLADDING HOOP STR AIN %)

Figure 7. F R AP-S3 Predicted Versus Measured Permanent Clat iing Hoop Strain 1568 308

SYM8OL EN REFERENCE O

96A,968,99-107 SAXTON 81(173 0

108-114.118.120,122 SAXTON lil(18I l 8-g 124,126.127,133,135,136 HPR-132(183 143.145-147,149-152

-t-205 207.211 AECL 1676(20) 7 1.6 -

X 213217 AE CL-1685(2 H Q

Q 233 239 AECL 2230t 23) 9 251-253.256 258.264 GE AP-4264(24)

E 200 293 PRTRt32)

~

316,317,348 IF A-208(35)

G G

411414 M T R( 53.54) 3 b43 548 BIG ROCK POINTt 81.82)

~U 1.2 -

g e

566-577 YBWR, DR ESDEN'84)

E

.3 V

O 620 626 IF A 106t 92) y A

647457 BRP(MO )I'#'

2 g

1.0 e

346,347 I F A-150( ddl 8

551,552 EPRl(833 h_

XX 595,599 I F A-411185.86)

U 0.8 -

00 5

/

U n

/

H 0.6 -

5 l

g O'

p P=M U

2 0.4 -

s

/

c E

d595 F 0.2 -

GG fo S

@ ~\\"

0.0 a.

,gm, p

-0.2 -

4a

/

-0.4 -

p

/

-0.6 /

-0.6

-0.4

-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 MEASURED PERMANENT CLADDING STRAIN (%)

Figure 8. FRAP-S3 Predicted Versus Measured Permanent Cladding Axial Strain.

1568 309

TABLE I FUEL BEHAVIOR VERIFICATION:

FRAP-S3 STANDARD ERROR

SUMMARY

Sample Size Output Parameter' (Rod / Pts)

Standard Error Fuel Center Temperature 33/290 254 K (Pressurized Rods) 64/511 198 K (Unpressurized Rods)

Released Fission Gas 176/176 18.8% Generated Gas Rod Internal Pressure 28/309 0.66 MPa (Unpressurized) 20/349 1.34 MPa (Pressurized)

Gap Closure Heat Rating 77/77 13.4 KW/m Axial Fuel Thermal Expansion 19/173 0.37% Active Length Permanent Fuel Axial Deformation 100/368 0.44% Active Length Permanent Cladding Hoop Strain 170/393 0.59% Cladding Diameter Permanent Cladding Axial Strain 115/161 0.47% Active Length Cladding Surface Corrosion Layer 48/84 6.6p Ln c7s Cladding Hydrogen Concentration 38/53 39 PPM CD u

CZ) f b)