ML18092A999: Difference between revisions
StriderTol (talk | contribs) (Created page by program invented by StriderTol) |
StriderTol (talk | contribs) (Created page by program invented by StriderTol) |
||
| Line 41: | Line 41: | ||
==Dear Mr. Varga:== | ==Dear Mr. Varga:== | ||
FINAL RULE | FINAL RULE 10CFR 50.61 FRACTURE TOUGHNESS REQUIREMENTS FOR PROTECTION AGAINST PRESSURIZED THERMAL SHOCK EVENTS SALEM GENERATING STATION UNITS 1 AND 2 DOCKET NOS. 50-272 AND 50-311 As required by 10CFR 50.61, PSE&G submits the enclosed report, NFU-060 Revision 0 dated January 10, 1986, which gives values of reference temperature for pressurized thermal shock at present , and at the expiration of the operating licenses for Salem Units 1 and 2. The report includes the bases for the projection including assumptions regarding core loading patterns. | ||
If you need additional information, please let us know. | If you need additional information, please let us know. | ||
Sincerely, | Sincerely, | ||
Revision as of 19:10, 7 November 2019
| ML18092A999 | |
| Person / Time | |
|---|---|
| Site: | Salem  |
| Issue date: | 01/20/1986 |
| From: | Corbin McNeil Public Service Enterprise Group |
| To: | Berkow H Office of Nuclear Reactor Regulation |
| Shared Package | |
| ML18092B000 | List: |
| References | |
| NLR-N86010, NUDOCS 8601280218 | |
| Download: ML18092A999 (83) | |
Text
{{#Wiki_filter:Public Service Electric arid Gas Company Corbin A. McNeill, Jr. Public Service Electric and Gas Company P.O. Box236, Han cocks Bridge, NJ 08038 609 339-4800 Vice President - Nuclear January 20, 1986 NLR-N86010 Off ice of Nuclear Reactor Regulation U.S. Nuclear Regulatory Commission 7920 Norfolk Avenue Bethesda, MD 20014 Attention: Herbert N. Berkow Standardization and Special Projects Directorate Gentlemen: TOPICAL REPORT - "SALEM REACTOR PHYSICS METHODS" (NFU-0039) SALEM GENERATING STATION UNIT NOS. 1 AND 2 DOCKET NOS. 50-272 AND 50-311 PSE&G requests that the NRC review its licensing topical report "Salem Reactor Physics Methods" which will be used in reload analysis submittals for Salem Units 1 and 2 after September 1987. To permit us to meet this schedule your review and acceptance should be completed by June 1, 1987. The report identification symbol is NFU-0039 and the issue date is January 10, 1986. Enclosed are twenty-three copies of Volume 1: MODEL QUALIFICATION. This document describes the comparisons of physics calculations to measurements, and the determination of associated reliability factors. These factors are applied as uncertainties in the calculation of physics parameters associated with reload safety evaluations. Revision 0 of this document does not include verification of transient Power Distribution Simulation Capability. , This section, together with Volume II, will be provided in March 1986. These documents are described below.
-----**-~~~
/*
i 8601280218 860120 . I PDR ADOCK 05000272 P PDR
1-Mr. H. N. Berkow
- 1/20/86 Section 3.7 Verification of Transient Power Distribution Simulation Capability This revision describes the comparisons of physics model calculations to measurements of non-steady state power distributions.
VOLUME II: RELOAD SAFETY EVALUATION This document describes the methods and procedures for calculating physics parameters associated with reload safety evaluations for the Salem cores. These procedures utilize the models and reliability factors determined in Volume I. With regard to the application fee, it is requested that you apply the fee transmitted with LCR 85-20 dated October 16, 1985, which was withdrawn by letter dated November 19, 1985 (C. A. McNeill, Jr. to s. A. Varga). We would appreciate your acknowleging receipt of this submittal and advising if you can meet the schedule discussed above. If you require additional information please let us know. Sincerely, Enclosures C Mr. Donald c. Fischer Licensing Project Manager Mr. Thomas J. Kenny Senior Resident Inspector
Public Service Electric and Gas Company Corbin A. McNeil!, Jr. Public Service Electric and Gas Company P.O. Box236, Han cocks Bridge, NJ 08038 609 339-4800 Vice President - Nuclear January 20, 1986 NLR-N86009 Off ice of Nuclear Reactor Regulation U.S. Nuclear Regulatory Commission 7920 Norfolk Avenue Bethesda, MD 20014 Attention: Mr. Steven A. Varga, Director PWR Project Directorate #3 Division of PWR Licensing A
Dear Mr. Varga:
FINAL RULE 10CFR 50.61 FRACTURE TOUGHNESS REQUIREMENTS FOR PROTECTION AGAINST PRESSURIZED THERMAL SHOCK EVENTS SALEM GENERATING STATION UNITS 1 AND 2 DOCKET NOS. 50-272 AND 50-311 As required by 10CFR 50.61, PSE&G submits the enclosed report, NFU-060 Revision 0 dated January 10, 1986, which gives values of reference temperature for pressurized thermal shock at present , and at the expiration of the operating licenses for Salem Units 1 and 2. The report includes the bases for the projection including assumptions regarding core loading patterns. If you need additional information, please let us know. Sincerely,
.Ml - J. KNIGHT (ltr only)
EB (.BALLARD} Enclosure E:ITCSB (ROSA) PSB (GAMMILL)
,r- - - - - - - - - - - - - -
I 8601280209 860120 RSB (BERLINGER) PDR ADOCK 05000272 FOB (BENAROYA) P PDR
, ... Mr. Steven A. Varga 1-20-86 C Mr. Donald c. Fischer Licensing Project Manager Mr. Thomas J. Kenny Senior Resident Inspector
I I 0 I PS~G NFU-0039 Revision 0 I The En rgy People July 31, 1985 I I I SALEM REACTOR PHYSICS METHODS I VOLUME I I MODEL QUALIFICATION I NOTICE - I THE ATTACHED FILES ARE OFFIC I AL RECORDS OF T HE DIVISION OF DOCUMENT CONTROL. THEY HAVE BEEN CHARGED TO YOU FOR A LIMITED TIME PERIOD AND 'I MUST BE RETURNED TO THE RECORDS FACILITY BRANCH 016 . PLEASE DO NO T SEND DOCUMENTS OHARGED OUT THROUGH THE MAIL. REMOVAL OF ~ N Y PAGE(S) FROM DOCUMENT FOR REPRODUCTION M ST BE REFERRED TO FILE PERSONNEL . DEADLINE RETURN DATE I 1_ ______ RECORDS FACILITY BRANCH I I mtATORY oorm f1L( COPY Docket# ~0- DI I
-Oontrol # C:Z lo O I 'Z- "2" 0 '2.- / l I Date cf* Z.Q
- 46 ~ of Document RF.GULATORY DOCKET FILE
I NFU-0039 I Revision 0 July 31, 1985 -I I I SALEM REACTOR PHYSICS METHODS VOLUME I I MODEL QUALIFICATION I I Prepared by: R. T. Brown Date: I Senior Engineer Nuclear Department Prepared by: ~ Date: , . 12/!t/es-I R. S. Kent I Reviewed by: N ear A. Blake zrtm Senior Staff Engineer k Date: -/-1C-8fc,
----
clear Tee nology Engineer I Approved by: c-*
'
clear ~artmez:. t :;*
/1
.4 /*/L /
~ / ..!-
- t/ i;*~--
*#*
0.
-;,,..,
Date: I./ ' I .. *1~- 1
/
._ *~
(;:.
,_.
I J E. S. Rosenf ld Manager - Nuclear Fuel Nuclear Department I Copy No.~~-{>=-~~- I NFU31/l 1 I I UtATORY*DOCKET .FILE COPY I
I NFU-0039 Revision 0 July 31, 1985 ABSTRACT 'I This topical report describes the methodology used by Public Service Electric and Gas Company (PSE&G) to determine calcula-I tional uncertainties and the resultant reliability factors associated with the PSE&G ARMP reactor physics model of the Salem pressurized water reactors. I I I I I' I NFU31/l 2 i 1* I I
I NFU'-0039 Revision 0 July 31, 1985 TABLE OF CONTENTS PAGE
1.0 INTRODUCTION
1 2.0 OVERVIEW OF THE CALCULATIONAL MODEL 2 3.0 MODEL VERIFICATION AND RELIABILITY DETERMINATION 5 3.1 Rod Worth Benchmarking 7 3.2 Isothermal Temperature Coefficient Benchmarking 13 I 3.3 Doppler Coefficient Benchmarking 16 I 3.4 3.5 Isotopics Reliability Factors for Delayed Neutron 20 Parameters 25 I 3.6 Power Distribution Benchmarking 29 .I 3.7 Verification of Transient Power Distribution Simulation Capability
- 59
4.0 REFERENCES
60 APPENDIX A Statistical Methods for the Determination and Application of Uncertainties A-1 I APPENDIX B Computer Code Summary Description B-1 I ii NFU31/l 3
NFU-0039 Revision O July 31, 1985 *1.
.11 I
Table LIST OF TABLES ii I 3.0.l Reliability F~ctors and Biases for PSE&G Model Applied to Salem 6 l 3 .1.1 Dilution Mode Rod Worth Comparisons 9 3 .1. 2 Rod Exchange Rod Worth Comparisons 10 3 .1. 3 Rod Worth Reliability Factors 12 3.2.1 Measured and Calculated Isothermal Temperature Coefficients 15 3.3.1 Comparison of Measured and Calculated Doppler Test P~rameters 18 3.4.1 Comparison Between EPRI-CELL and SAXTON Experimental Data 21 3.6.1 Reactor State Points 36 3.6.2 Reactor State Points 37 3.6.3 Reactor State Points 38 3.6.4 Mean Observed Differences Axial Model 3.6.5 Bias Axial Region Definitions 48 49 I 3.6.6 Confidenc~ Limits for X(i,k,m) Distri-bution by Subgroup 52 I 3.6.7 Confidence Limits for X(i,m) Distribution by Subgroup 53 I I NFU31/l 4 ,- i i i
.I
I NFU-0039 Revision 0 July 31, 1985 LIST OF FIGURES Figure
- 2. 0. 1 Salem Physics Model 4
- 3. 3. 1 Comparison of Measured and Calculated Doppler Test Parameters 19 3.4.1 Comparison of EPRI - CELL to Yankee Pu239/Pu240 Isotopic Ratios 22
- 3. 4. 2 Comparisons of EPRI - CELL to Yankee Pu240/Pu241 Isotopic Ratios 23 3.4.3 Comparisons of EPRI - CELL to Yankee Pu241/Pu242 Isotopic Ratios 24
- 3. 6. 1 Salem Unit 1 and Salem Unit 2 Moveable Incore Detector Locations 34 3.6.2 Axial Locations of Grids and Detectors 35 3.6.3 Measured and Calculated Integrated Detector Responses Salem 1 Cycle 4 MAP 1411 39 3.6.4 Measured and Calculated Detector Responses Salem 1 Cycle 4 MAP 1411 40
- 3. 6. 5 Measured and Calculated Detector Responses Salem 1 Cycle 4 MAP 1411 41 3.6.6 Measured and Calculated Integrated Detector Responses Salem 1 Cycle 5 MAP 1522 42 3.6.7 Measured and Calculated Detector Responses Salem 1 Cycle 5 MAP 1522 43 NFU31/l 5 iv
NFU-0039 Revision 0 July 31, 1985 LIST OF FIGURES (continued) I Figure 1* 3.6.8 Measured and Calculated Detector 3.6.9 Responses Salem 1 Cycle 5 MAP 1522 Measured and Calculated Integrated 44 I 3.6.10 Detector Responses Salem 2 Cycle 1 MAP 2133 Measured and Calculated Detector 45 1* Responses Salem 2 Cycle 1 MAP 2133 46 3.6.11 Measured and Calculated Detector Responses Salem 2 Cycle 1 MAP 2133 47 I 3.6.12 Distribution of Errors X(i,k,rn) 50 I 3.6.13 Distribution of Errors for Integral 3.6.14 X(i,rn) Confidence Levels for X(i,k,rn) versus 51 I. Reactor Power % 54 3.6.15 Confidence Levels for X(i,k,m) versus Cycle Exposure 55 3.6.16 Confidence Levels for X(i,k,rn) versus Axial Height 56 3.6.17 Confidence Levels for X(i,rn) versus Reactor Power 57 3.6.18 Confidence Levels for X(i,m) versus Cycle Exposure 58 NFU31/l 6 v
I NFU-0039 Revision 0 July 31, 1985
1.0 INTRODUCTION
This report describes the Salem reactor physics model and addresses the qualification and quantification of reliability factors for application of the model to operations and reload safety evaluatio*ns of the Salem Nuclear Reactors. A summary description of the computer codes used to model the Salem reactors is given in Section 2. The qualification of the model is described in Section 3. Whenever possible, directly observable parameters (such as rod worths, and incore detector fission rates) are utilized for this qualification. The data used in this evaluation I span seven (7) reactor operating cycles. The reactor cycles included are Cycles 1 through 5 for Salem Unit 1, and Cycles 1 and 2 for Unit 2. I After the measured data to be* used in the benchmark process are defined, the model calculations are performed and are compared to measurements. These comparisons are presented in this report as part of the quantification of the PSE&G model calculational uncertainties and reliability factors. A statistical approach is used to derive the uncertainties and reliability factors. These uncertainties and reliability factors are consistent with the model application procedures and methodology. I The uncertainties and reliability factors are evaluated by direct comparison to experimental data. I In. order to provide. a continuing verification of t.he conservatism of the reliability factors determined herein, ongoing comparisons are made each cycle using statistical methods consistent with those described in this report. I* I NFU31/l 7 1 I
I NFU-0039 Revision 0 July 31, 1985 2.0 OVERVIEW OF THE CALCULATIONAL MODEL The model used to analyze the Salem Units was constructed using the Advanced Recycle Methodology Program (ARMP) system developed under EPRI sportsorship by UAI. (Reference l) A flow diagram for this model is shown in Figure 2.0.l. The spectral code, EPRI-CELL (ARMP, Part II, Chapter 5), produces initial nuclide concentrations, depletion and fission product chain data, and tables of microscopic and macroscopic cross sections varying with burn-up for input to the XY diffusion - depletion code, PDQ7/HARMONY (Reference 2 and 3). Lumped absorber data for PDQ7/HARMONY are generated by a capture fraction matching procedure between PDQ7 and either EPRI-CELL (ARMP, Part I, Chapter 6,
- 1 Section 4) for burnable poisons or CPM (ARMP, Part I, Chapter 6, Section 3) for control rods. PDQ7/HARMONY is run both in the full core (XY) geometry representation and the fuel type (color set) representation. The full core I. representation is used for nodal code normalization, local peaking factor generation, and for the establishment of assembly loading patterns.
/I, In the fuel type (color set) mode, PDQ7/HARMONY supplies
., input data for PSE&G's nodal code, TRINODE, a derivative of the EPRI-NODE-P program (ARMP, Part II, Chapter 14). The I TRINODE program contains improvements over the EPRI-NODE-P program which include input/output changes, execution options, and file management. However, the primary 1* calculational sequence and physics methodology have been preserved from the EPRI-NODE-P program.
It is recognized that the methods used ~or the construction I and application of the Salem model are as much a part of the model definition as are the codes. It is essential, therefore, that the methods used to calculate core safety margins be consistent with those used in the model benchmarking and qualifications process. This is particularly true in the calculation of core power distribution and local peaking factors in which the results are heavily dependent on the methods used to normalize the nodal model. NFU31/l 8 2 I
I f'i!FU-0039 Rev is ion .0
,Tu 1 y 31 , 1 985 I I
2.0 OVERVIEW OF THE CALCULATIONAL MODEL (continued) *1 The TRINODE model is normalized to the PDQ model. A consistent methodology is used for this normalization throughout the benchmark calculations and will be used in future safety related calculations. I In addition to the main sequence computer codes, a number of auxiliary computer codes are employed to provide a user tailored code package. These auxiliary computer codes are I not basic to the physics methodology, but are vital for automation and transformation of the large volume of calculated and measured parameters required for core analysis. The auxiliary computer codes are summarize~ in Appendix B. All comparisons to measurement data presented in this topical report are based on TRINODE calculations. NFU31/l 9 3
I I FIGURE 2.0.1 NFU-00 39 Revision 0 SALEM PHYSICS MODEL July 31, 1985
,I,
-1 I EPRI-CELL CPM
,i
\
EPRI-CELL /I NUPUNCHER LUMPED ABSORBER PROCEDURE 1- " PDQ/HARMONY 1 FUEL TYPE (COLOR SET) FULL CORE I. , 1 EPRI FIT
'.
' I NORMALIZATION I REACTION RATES
,. SUPER LINK PIN TO BOX I*
TRI NODE Ii I SIGMA 1<11-r-----' I r PLANT I MEASUREMENTS I 4
I NFU-0039 I *Revision 0 July 31, 1985 ,, I 3.0 MODEL VERIFICATION AND RELIABILITY DETERMINATION The PSE&G model is benchmarked against Salem Unit 1 I measurements made during Cycles 1 through 5, and Salem Unit 2 measurements made in Cycles 1 and 2. This benchmark ,,, serves as the basis to quantify the reliability factors to be used in safety related calculations. The term reliability factor (RF) is used to describe the /I allowances (either absolute or relative) to be used in safety related calculations to assure conservatism. The term uncertainty factor is used to describe the actual model precision and is defined as the standard deviation I (a). The reliability factor is always larger than the uncertainty factor. The term bias is used to describe the statistical I difference between an observed.or measured distribution and the calculated value. _I, Table 3.0.l summarizes the model reliability factors and ,, biases computed as a result of ihe model benchmark. remainder of Section 3 is a detailed account of the derivation of these factors. The statistical metQods employed are described in Appendix A to this report
- The
.1. I I I; I I I NFU31/l 11 5 I
NFU-0039 Revis.ion o July 31, 198S TABLE 3.0.1 RELIABILITY FACTORS .AND BIASES FOR PSE&G MODEL I APPLIED TO SALEM PARAMETER RELIABILITY FACTOR I BIAS
. Rod worth .J MEAS>600 PCM RF ROD = 1S% 0 MEAS(600 PCM TOTALS RF ROD RF ROD
=
=
lOOPCM 10% 0 0 I
. Temperature Coefficient \I Moderator (MTC) RFMTC = . 2.1 PCM/°F 0 Isothermal (ITC)
Doppler RF ITC RFDC
=
=
- 2. 1 PCM/°F 10%
0 0 I
.
Doppler Defect Delayed Neutron RFDD = 10% 0 I Parameters B eff Q.,
- RFB RFL
=
=
4% 4% 0 0 I
. Power Distribution FQ P>
P<
.so
.so RFFQ RFFQ
=
=
0.10 0.16 -(0.12*P)
**
**
,,
I. F /':, H P> P<
- .30
.30 RFF H RFF H
=
=
0.08 0.09 -(P/30) 0 0 I
** See Table 3.6.4 t
I NFU31/l 12 6 :1
!I
I NFU-0039 Revision 0 July 31, 1985 I 3.1 Rod Worth Benchmarking The purpose of this section is to benchmark the PSE&G
- 1 Sa~em model to rod worth measurements. This is accomplished by first presenting and qualifying the available measurements and second by computing model I. reliability factors.
Rod worth measurements have been performed at Salem usinq two techniques; the boron dilution method and the rod exchange technique. (Reference 9) Boron dilution rod worth measurements were performed on Unit 1 Cycles 1 through 5 and Unit 2 Cycles 1 and
- 2. The results of these measurements are summarized on Table 3.1.1 along with PSE&G model predictions.
Rod exchange measurements were performed on Unit 1 Cycles 1, 3, 4, and 5, and Unit 2 Cycle 2. These results, along with model calculations are tabulated on Tables 3.1.2 (a) and (b). For purposes of model benchmarking, some rod worth measurements are disqualified on the basis of known measurement. errors. Measurements disqualified are I the boron dilutions for Unit 1 Cycles 1 and 2, and Unit 2 Cycle 1. Additionally, rod exchange measurements for Unit 1 Cycle 1 are disqualified. The basis for this disqualification is measurement errors I discovered in dilution measurements made prior to Cycle 3. These errors are due to the effects of spatial flux redistribution caused by rod motion during the dilution (Reference 9). Test procedure changes were implemented prior to Cycle 3 measurements to reduce these effects. Since rod exchange measurements use the reference bank dilution measurement to interpret exchange worths, rod exchange measurements for Unit 1 Cycle 1 are disqualified on the same basis. I Support for the disqualification of dilution measurement made prior to Cycle 3 is available using I comparisons to calculated worths. The average difference between measured and calculated rod worths for dilution measurements performed prior to Cycle 3 and those using the improved test procedure are 11% I and 1% respectively. This difference is significant at the 99.9% confidence level, and is attributed to the known measurement errors. I I NFU31/l 13 7
I NFU-0039 Revision 0 July 31, 1985 I
.1 Rod worth reliability factors were obtained by bounding the results of the comparisons between measured and calculated rod worths. These factors are tabulated on Table 3.1.3.
I Comparisons were taken from 7 dilution measurements and 24 exchange measurements spanning 4 reactor cycles, and represent all measurements through Cycle 5 of Unit 1 and I Cycle 2 of Unit 2, except those disqualified above. Using a bounding value for the reliability factor is justi-J, _fied due to its conservatism relative to normal statis-
*tics. Calculation of reliability factors representing 95/95 confidence levels using normal statistics yields a 91 PCM reliability factor adder for rods worth less than 600 I
PCM, and 12% reliability factor for rods worth more than 600 PCM. The exception is the reliability factor for rod worth totals, which is computed to be 16%. However, this large value is due to the small sample size of only 5 values. Since the error for rod worth totals can be no larger than the largest error for the individual rod banks, the reliabilty factor for rod worth totals should be
,.
I bounded by the maximum observed error for individual banks with worth greater than 600 PCM; 10%. Thus, the factors tabulated on Table 3.1.3 conservatively bound the observed J, data, and will be used as model rod worth reliability factors. I I:
*1.*
.I I'
NFU31/l 14 8 ,, I
*1
I NFU-0039 J Revision 0
,July 31, 1985 I
1 TABLE 3. 1. 1 /I DILUTION MODE ROD WORTH COMPARISONS ,, DATE UNIT/ CYCLE BANK MEAS (PCM) CALC (PCM) DIFFERENCE** M>600
%
M<600 {:, /I 12/76 1/1*** D c
.1107 1183 1030 1005 7.5 17.7 ,, B A
SD SC TOTAL 766 1241 745 1181 6223 724 1114 681 1060 5614 5.8
- 11. 4 9.4
- 11. 4 10.8 1, 12/79 1/2*** D c
1041 938 924 846 12.7 10.9 ,,, B A TOTAL 534 1163 3676 599 973 3342 19.5 10.0
-65 8/80 2/1*** D 1391 1241 12. 1 I c B
A 1185 1359 501 1026 1262 385
- 15. 5 7.7 116 5.3 I SD SC TOTAL 750 1052 6238 712 961 5587
- 9. 5
- 11. 7 I 12/80 1/3 D c
B 834 960 565 797 900 600
- 4. 6 6.7
-35 I A TOTAL 1023 3382 1058 3355
-3.3 0.8 4/82 1/4 D 862 860 0.2 2/83 1/5 D 926 939 -1. 4 7/83 2/2 0 878 835 5.1
* % = ( ( M-C) /C)
- 100 for measurements >600 PCM.
{:, = (M-C) for measurements (600 PCM.
*** Data disqualified as discussed in text.
I NFU31/l 15 9
I NFU-0039 Revision 0 July 31, 1985 I l TABLE 3.1. 2 ROD EXCHANGE ROD WORTH COMPARISONS I I DATE UNIT/ CYCLE BANK MEAS PCM CALC PCM DIFFERENCE M>600
-
M< 600 t::.
.I
%
I 12/76 1/1*** *D c B 1107 825 522 1030 741 467 ( 7. 5 )
- 11. 3
- 55
,,
,,
A 924 858 7.7 SD 469 403 - 66 SC 351 305 - 46 TOTAL 4198 3804 10.3 12/80 1/3 *D 834 797 ( 4. 6) J,
,,
c 696 674 3.3 B 395 450 - -55 A 816 789 3.4 TOTAL 2741 2710 -1. 1 J 4/82 1/4 *D 862 860 ( 0. 2) c 596 588 - 8 370 - ~1
,,,
B 407 -37
.A 818 789 3.7 SD 265 316 - -51 SC 285 281 - 4 SB 614 649 -5.4 SA 750 733 2.3 TOTAL 4560 4623 -1. 4 I
2/83 1/5 *D c B A SD 926 613 331 784 269 939 617 361 814 292 (-1.4)
-0.6
-3.7
-
-
-30
-23
,,
SC SB SA 317 769 735 291 793 779
-3.0
-5.6
- 26 II*
TOTAL 4 744 4886 -2.9 (continued) *I' NFU31/l 16 10 i I
- - - - - - - - - - - - - - - - - ------ - --------- I NFU-0039 I Revision 0 July 31, 1985 I I TABLE 3. 1. 2 ROD EXCHANGE ROD WORTH COMPARIOSNS
'I (continued) .I. DATE UNIT/
CYCLE BANK MEAS (PCM) CALC (PCM) DIFFERENCE* M>600
%
M<600
~
7/83 2/2 *D 878 835 ( 5. 1) c 770 731 5.3 B 660 603 9.5 A 252 233 19 SD 299 287 12 SC 292 275 17 SB 787 757 4.0 SA 562 491 71 TOTAL 4500 4212 6.8
.,. ,I
- Measurement perfor!l1ed by dilution
** % = ((M-C)/C)
- 100 for measurements >600 PCM
~ = (M-X) for measurements <600 PCM
*** Data disqualified as discussed in text.
I f, I* J I" I I 11 NFU31/l 17
NFU-0039 I Revision 0 July*31, 1985 I 1 TABLE 3. l . 3 RODWORTH RELIABILITY FACTORS I Individual Rod Worth I a) Rod worth <600 pcm RF ROD = 100 pcm I b) Rod worth >600 pcm RF ROD = 15% I Total Rod worth RF ROD = 1U% I\
*I
,,,
I'
.,
I
\I
,,
NFU31/l 18 12
,,
I
*i I
t NFU-0039 I Revision 0 July 31, 1985 f, l 3.2 Isothermal Temperature Coefficient Benchmarking The objective of this section is to benchmark ~he PSE&G I model to measured isothermal temper~ture coefficients (ITC). Based on comparisons between measured and calculated coefficients, a reliability factor for both the isothermal and the moderator temperature coefficient I/ (MTC) is inferred.
- A total of 19 ITC measurements are tabulated on Table
- 3. 2. 1. These measurements span 7 reactor cycles and range from unrodded conditions to all control banks inserted.
The PSE&G model calculations for ITC are presented on Table 3.2.1 along with the corresponding measurement. Statistical tests were performed on the comparisons to evaluate normality and pooleability. Normality was demonstrated using the W-test (Reference 8), while pooleability was assured using the Bartlett test (Reference 4). The computed standard deviation of the comparisons between measured and calculated ITC's is 0.85 PCM/F. The observed standard deviation of 0.85 PCM/F ( oOBSV) is assumed to be made up of three independent components; measurement uncertainty, model calculational uncertainty on moderator te~perature coefficient, and, model calculational uncertainty on Doppler temperature coefficient. This relationship is expressed as: a 2 + a 2 + a2 a 2 = (0.85)2 I MEAS MTC DC OBSV Since each component is greater than or equal to zero,each component is bounded by the observed error, I Therefore, a conservative estimate of the model uncertainty ( a ) for both the isothermal and moderator temperature coefficients is 0.85 PCM/F. This is I* summarized as:
= 0.85
= 0.85 l
11 13 NFU31/l 19 I I
NFU-0039 I Rev is ion 0 Ju 1 y 31 , 1985 I J PSE&G model reliability factors for both ITC and MTC are computed as the product of the standard deviation and the one-sided critical factor (K c) for a 95/95 1*: confidence level using nineteen (19) samples. This product yields reliability factors for.ITC and MTC of 2 .1 PCM/F. I\ RFITC = 0.85
- 2.42 = 2.1 PCM/F RFMTC = 0.85
- 2.42 = 2.1 PCM/F J
I\ l\
,,,
J, I' NFU31/l 20 14
,,
I;* NFU-0039 I Rev is ion 0 July 31, 1985 I' TABLE 3. 2. 1 I MEASURED AND CALCULATED ISOTHERMAL TEMPERATURE COEFFICIENTS I ROD POSITION ITC PCM/°F I/ UNIT 1 CYCLE 1 BANK D (STEPS) 197 BORON 1369 MEAS
-3.51 CALC
-3.59 DIFF 0.08 c 201 1264 -4. 11 -4.34 0.23
'I B 175 1151 -6.17 -6.45 0.28
,, A SD 175 175 1085 965
-7.85
-11.25
-9.06
-11.90
- 1. 21 0.65
., 1 1 2 3 D c D c 219 214 219 206 1137 1025 1258 1157
-6.06
-5.79
-3.33
-4.85
-4.45
-5.45
-3.26
-4.10
-L61
-0.34
-0. 07
-0.75
.I. 1 1 4 5 D D 202 214 1309 1499
-3.61
-1. 52
-5.29
-2.60
- 1. 68
- 1. 0 8
- 1 2 1 D D
205 188 1334 1329
-0.65
-0.84
-0.59
-0.70
-0.06
-0. 14 D 102 1285 -2.68 -1. 89 -0.79 I c 184 1197 -4.34 -4.85 -0.51
., 2 2 B
A SD D 203 198 192 218 1083 955 910 1362
-10.53
-10.50
-13.48
-4.16
-9.09
-9.83
-13.22
-4.55
-1. 44
-0. 67
-0.26 0.39
.I Mean Standard Deviation 0.00
- 0. 8 5 I\
l I~- NFU31/l 21 15 I I
I: NFU-0039 I Revision 0 July 31, l<:J85 I 1 3.3 Doppler Coefficient Benchmarking The objective of this section is to make comparisons between measured and calculated Doppler coefficients and I
,,
establish ~odel reliability factors for Doppler reactivity calculations. Doppler reactivity coefficient measurements have been performed at the Salem Units using two measurement techniques. Both test procedures require the compensation of a reactivity imbalance induced by a reactor power change. The first test procedure balances reactivity using control rods and measures the reactivity changes using a reactimeter. This technique was used in Cycle 1 of both Salem Units. The second measurement procedure maintains reactivity balance with changes in moderator temperature. The ratio of power change to moderator temperature change is then converted I to reactivity using an isothermal temperature coefficient. This technique has been used for Cycles 2 through 5 on Unit 1, and Cycle 2 for Unit 2. l The measurements using rod banks for reactivity control are not used for purposes of model benchmarking. The basis for this disqualification is the large uncertainties associated with reactimeter interpretation for at-power measurements. I The results of all Doppler coefficient measurements performed using the moderator temperature control procedure have been tabulated on Table 3.3.1. This measurement technique requires calculated isothermal I.- temperature coefficients to infer the Doppler coef~icient temperature. Since it is the ratio of the changes of these two quantities that is actually
.I measured, this ratio is tabulated along with the inferred Doppler coefficient on Table 3.3.1. The precision associated with each measured ratio has been determined based on the standard deviation of multiple I measurements.
I NFU 31/1 22 16 I 1 I
NFU -0 0 39 Revision 0 July 31, 1985 I t Calculations of the ratio of power to moderator temperature changes have been made using the PSE&G model. Comparisons of the measured and calculated ratio are shown on Table 3.3.l and also Figure 3.3.1 in which the vertical bars represent the measurement precision. Figure 3.3.1 demonstrates that the measured and calculated ratios typically agree to within the measurement precision, and therefore confirms model capability to calculate these ratios. The scatter in the J data shown in Figure 3.3.1 is due primarily to the poor measurement precision. It is apparent from Figure 3.3.1 that the measurement precision is of the same order of magnitude as the observed differences between measurement and calculation. Thus, the model calculational uncertainty is assumed to be small. For purposes of assigning a model reliability factor for*Doppler coefficient (RFDC), a conservative value of 10% is assumed. The same reliability factor will be -assigned to the moael for Doppler only power defect (RFDD). Thus: RFDC = 10% RFDD = 10% I I I I f I t NFU 31/1 23 17 I
zt-z:j d w f-' '-.. f-' N Vl TABLE 3.3.1 COMPARISCN OF MEASURED AND CALCULATED DOPPLER TEST PARAMETERS UNIT/ PWR NUMBER an (6 p ) MEAS (6 P) CYCLE % OF MEAS PCM/% (6 T )MEAS PRECISION rrT"JCALC MEAS-CALC 1/2 39 6 -13.67 -0.95 0.04 -0.86 -0.09 93 6 -13 .15 -1.39 0.09 -1.41 0.02 f-' 1/3 44 6 -13 .31 -0.77 0.02 -0.83 0.06 ()) 94 6 -10.91 -1.40 0.17 -1.28 -0.12 1/4 43 4 -10 .11 -0.90 0.23 -0.84 -0.06 99 4 -11.49 -1.28 0.09 -1.31 0.03 1/5 46 4 -11.45 -0.66 0.04 -0.70 0.04 97 4 -12.01 -0.99 0.11 -1.10 0.11 2/2 98 2 -11.33 -1.34 0.35 -1.33 -0.01 y !:ti z i::: CD t-z:l f-' -G d
~I-'* I Ul 0 Wl-'*O f-' 0 w
.. ::J ~
f-' 0
~
()) Vl
NFU-0039 Revision 0 FIGURE 3.3.1 July 31, 198~ I COMPARISON OF MEASURED AND CALCULATED DOPPLER TEST PARAMETERS l1 I
- 1. 7 I
- 1. 6 J
- 1. 5
- 1. 4 I
- 1. 3
,I µ..
0 .......... 0\0 *(
. 1. 2 E-i <J
.......... r:i..
<J 1.1
,,,
'd I
Q) H 1. 0
- l
,.,
Ul rd Q)
- 2: 0.9 0.8 l
0.7 0.6
,,
0.5 . ~~--'~~-'-~.......JL.-~.......~~..._~_...~~..a...~--1~~--~~ 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5
*f I
Calculated 19
~P/ ~T, %/°F
.,
I
I NFU-0039 I Revision 0 July 31, 1985 I I 3.4 Isotopics Isotopic compositions calculated by EPRI-CELL have been I compared with spent fuel isotopic data obtained from Yankee Rowe fuel rods irradiated beyond 35 GwU/MTU. The reactor representation used for the EPRI-CELL benchmarking in the I calculations is described in the ARMP documentation (Part 1, Chapter 1, Section 4.0). Experimental and analytic isotopic ratios for plutonium from
- 1 the Yankee Rowe spent fuel are plotted versus accumulated fissions in Figures 3.4.1-3.4.3. The dots are experimental results and the line the EPRI-CELL results.
The agreement between calculated and experimental isotopic ratios is good. The calculated ratios Pu-239/Pu-240 and Pu-240/Pu-241 are within the scatter of the experimental results and the ratio Pu-241/Pu-242 is slightly over-predicted. Calculated and measured isotopic compositions for MU2 fuel (6.6w/o) irradiated in the Saxton Core II (pellet, rod MY, zone 6) are compared in Table 3.4.1. The agreement is good for the most important uranium and plutonium isotopes as well as for americium and curium. The measured burnup ranged from 7 to 22 GWD/MTU. The reactor representation for the EPRI-CELL benchmarking is described in the ARMP documentation (Part 1, Chapter 3, Section 4.) I 1* I .,
'-..."
..
l I NFU31/l 26 20 I I
I NFU-00 39 Revision 0 July 31, 1.985 I l TABLE 3.4.l I' COMPARISON BETWEEN EPRI-CELL AND SAXTON EXPERIMENTAL DATA 1* EXPERIMENTAL UNCERTAINTY ~PRI-CELLJ *100% NUCLIDE EXPERIMENT % Exp. I I
,,
ATOM ~ U-234 .00465 28.7 - 3.7 U-235 .574 .9 + 1. 7 U-236 .0355 5.6 - 7.6 99.386 0 0 U-238 Pu-238 Pu-239
.109 73.77 2.2 0
-32.1
+ 0.6
*I' Pu-240 19.25 .2 - 3.5 Pu-241 Pu-242 6.29
.579
*3
*9
+ 4.6
-11.7 I, I
Np-237/ U-238 1.14 x lo-4 ATOM RATIOS 15 -26.9
,,
Pu-239/ .04383 .7 + 2.0 U238 Pu-238/ Pu-239
- 1. 75 x lo-3
.4 -17.6 I Am-241/
Pu-239
.0123 4
15 10
+ 2.4
+ 0.4
,I_
Cm-242/ 1.05 x lo-Pu-239 Cm-244/ Pu-239 1.09- 4 20 - 3.2 if
*1; NFU31/l 27 21
'I I
I NFG -0039 Revision 0 July 31, 1985 FIGURE 3.4.l COMPARISON OF EPRI-CELL TO YANKEE Pu239/Pu240 ISOTOPIC RATIOS 10.0 .,
,.
'
9.0
\
s.o
\..
7.0
\ **'
.
0 H 0
~ '<l' N
6.0
\ ..
'.
\.
I ::> A<
......... ....,
O'\ N
'\
- >
~
I A< 5.0 I ~
~
I 4.0
~
~.
I 3. 0
*'
I I o. J o.o 5.0 10.0 15.0 20.0 25.0 30.0 I Accumulated Fissions (barn-cm)
-1 x 10 5
22 I
NFG-0039 Revision 0 I July 31, 1985 FIGURE 3.4.2 COMPARISON OF EPRI-CELL TO YANKEE Pu240/Pu241 ISOTOPIC RATIOS I 9.0 l I a.a
,,.
I.
\
I I 7.0 I* i
-.*
*\ l 6,0 I
0 H I I ~ 5.0 I
\ 'I
\\* ..
I II
\\ .I 1..:*
3.0
...
2.0
..
.... .
f
~
~ . .
~ .p. .....
I
.. . - .
. 1.0 I
-
o.o s.o 10.0 15.0 20.0 25.0 30.0
-1 5 Accmnulated Fissions (barn-cm) x 10 23
I NFG-0039 Revision 0 I FIGURE 3.4.3 July 31, 1985 COMPARISON OF EPRI-CELL TO YANKEE Pu241/Pu242 ISOTOPIC RATIOS I, I 10. a ..
. ...
"
I
- 9. a I *'
.\
--
'I
- 8. a
,f u:i 'I 0 H
~
N
- 7. 0
.... ~ \
,I '<l' N
- i .
I Pl
....__ ~ '<l' N
- i 6. a I\
'~
Pl
..
,, I
- 5. 0 :
""-~
...
I* 4. 0
...
I t a. a a.a 5.0 10.0 15.0 2a.a 25.0 3a.a I Accumulated Fissions (barn-cm)-l x ia 5 I I 24
I NFU-0039 I Revision 0 July 31, 1985 I I 3.5 Reliability Factors for Delayed Neutron Parameters This section deals with determining reliability factors for the effective delayed neutron fraction and I the effective neutron lifetime which are values which can be calculated but whose measurement is not practical. In these cases, an argument is made for the I general magnitude of the reliability factor without making direct comparisons between measured and predicted values. The importance of the reliability of the calculated values of the delayed neutron parameters is primarily associated with the core Beff* The .I uncertainties in the calculation of Beff composed of several components, the most of which are listed below: are important I a. b. Experimental values of S, and A, by nuclide; Calculation of the spatial nuclide inventory;
- c. Calculation of core average S as a flux weighted
.l d. average over the spatial nuclide inventory; Calculation of Beff from the core averaqe as Beff = I*S, where I = importance factor. The experimental determination of the S's and A's are assumed to be accurate to within 1%. The most important nuclide concentrgtions with respect to core I S are u23 8 , u 235 , and Pu 23 -
- Tables 3. 4. 1 and 3.4.2 indicate that the difference in the calculation of these ~oncentrations is about 1.7% for ECELL.
Therefore, components (a) and (b) above are combined as
- 2. 7 %*
I NFU31/l 31 25 I I
NFU-0039 Revision 0 Ju 1 y 31 , 1'9 8 5 The uncertainty in the calculation of a core average S depends on the relative flux weighting of the individual assemblies in the core. For demonstration purposes, consider a three region core, each with a different average burnup and average S. This is typical of advanced PWR cycles in that about a third of the core has seen two previous cycles, a third only one previous cycle and a third is the feed fuel. Typical regional B 's are given below: Region 1 (third cycle fuel) S (1) = 0.005 Region 2 (second cycle fuel) S (2) = o.006 Region 3 (feed fuel) B (3) = 0.007 The effect of errors in the calculated flux distribution can be evaluated in terms of the effect on the core average s. As a base case, flux weighting I, factors (FWF) are all set to 1.0. In this case, the core average S = 0.006. Using a maximum error in the calculation of the core average is obtained by increasing the weight of the Region 1 fuel and I decreasing the weight of the Region 3 fuel. revised S is calculated as follows: The I. (l)xl.07 = .00535 (2)xl.O (3)x0.93
=
=
.0060
.00651 I S = .00595, which yield a -0.8% error for component (C) above. I NFU31/l 32 26 I
I
I NFU-0039 I Revision 0 July 31, 1985 I The last uncertainty component, (d), concerns the I reduction of core average B to obtain Beff by using the importance factor. Since this reduction is typically about 3% to 4%, an error of 10% in this I component would lead to an error in Seff of less then 0.5%. The sum of the errors for these four factors for I ECELL are as follows: 2.7%(a+b) + 0.8%(c) + 0.5%(d) = 4.0% I So the reliability factor for delayed neutron parameters (RFB)is set at 4%. I An argument similar to the delayed neutron parameter argument is applied to the determination of the effective neutron lifetime ( i*) uncertainty. The 1 uncertainty components which go into the calcul~tion of i* are as follows: (a) Experimental values of microscopic cross I sections; (b) Calculation of the spatial nuclide inventory; I (c) and Calculation of the core average effective neutron lifetime, i*, as a flux weighted average over the l spatial nuclide inventory which includes the effects of leakages. J Uncertainties for components (a) and (b) are assumed to be the same as described for the calculation of S eff, that is, a combination of 1% uncertainty in the experimental determination of nuclear cross I sections and 1.7% uncertainty in the determination of I t I NFU31/l 33 27 I I
I NFU-0039 Revision 0 Ju 1 y "3 1 , 1 9 8 5 I I the spatial nuclide inventory of ECELL. The core average neutron lifetime depends on flux weighting of local absorption lifetimes i*. If a conservative I estimate of the error in regional power sharing (7%) is used in determining the impact on the core average lifetime ( i*), the error in lifetime is on the order I of 1.0%. Combining all of these uncertainties linearly results in a total uncertainty of 3.7%. Therefore, a 4% reliability factor (RFL) will be applied to the neutron lifetime calculation when applied to safety I related calculations. I I I I l
"
.,
I
*I NFU31/l 34 28
,.
I I
I I NFU-0039 Revision 0 July 31, 1985 I I 3.6 Power Distribution Benchmarking I It is the purpose of this section to quantify the PSE&G Salem model power distribution calculations. This is accomplished by first presenting the measurement data base, followed by a description of the calculational_ I methodology. Second, comparisons are made between the measured and calculated quantities, and lastly, model reliability factors for power distribution calculations I are computed. The primary source of power distribution measurements for Salem Units 1 and 2 is the incore detector system. I This system consists of moveable incore fission chambers which respond to neutron flux. These neutron detectors traverse through instrument guide thimbles I which are located at 58 positions throughout the core as shown in Figure 3.6.1. Measurement signals from these detectors are taken at 61 axial positions up the fuel assembly as illustrated in Figure 3.6.2, and are I corrected by the on-site process computer to account for detector sensitivity, drift, and background. The corrected signals are then -usetj to compute "measured" I power distributions using analytical data to convert the detector signals to interpreted powers in both instrumented and uninstrumented assemblies. I A total of forty-nine (49) flux maps were chosen for the purpose of benchmarking the PSE&G Salem model. These flux maps span six reactor-cycles and represent I typical steady state operation conditions. These include maps taken at powers ranging from near zero to 100 percent, and cycle exposures of zero to end of cycle, including some coastdown state points. A I description of reactor conditions for each flux map chosen is given in Tables 3.6.1 through 3.6.3. I The approach taken to benchmark and qualify the PSE&G Salem models for power distributions was to compare calculated and measured detector signals. The basis for this is twofold. First, the detector signals t I NFU31/l 35 29 I I
NFU-00 39
- Revision 0 July 31, 1985 I I
represent raw measurements and do not include interpretation, unlike "measured" power distributions. Second, the ability of the model to compute the detector signal requires the same processes as required to compute pin powers. Both calcuations require the prediction of the localized fission rate, one in a pin pellet, the other in a fission chamber. The accuracy of the two calculations is essentially the same. The only difference is that there is a small self-shielding or flux depression in the pin which is not in the I detector. Th~ impact of this difference on the power distribution reliability factor is assumed to be negligible. The simulated detector signals are calculated in a manner which is consistent with the calculation of local power peaking factors for the purpose of safety evaluation. The first step is to compute the power distribution under consideration. The resolution used is one node per fuel assembly, with 12 axial levels. The simulated detector signals are obtained by using the nodal power at each axial level to predict a signal power density for that assembly a~ that level. This power is then converted to a relative reaction rate. The conversion factors are calculated for each asembly location as a function of assembly exposure using a two-dimensional, full core PDQ7, fine mesh model. The 12 axial values in each assembly location are then synthesized using a truncated fourier sine series. I Grid flux depressions are then superimposed on the synthesized function using an empirical function designed to match the characteristics of flux depressions measured with in-core fission detectors. I The effect of the grid flux depressions is to raise the flux level in the axial region between grids while depressing the flux in the grid region. Consistency I between the above calculations of instrument signals and the calculation of local peaking factors is assured by: I t I NFU3l/l 36 30 I I
I NFU-0039 I Revision 0 July 31, 1985 I I A. B. C. Using a common Using a common Using a common full core PDQ7 model, nodal model, and procedure to account for axial J. flux gradients and grid effects. Typical comparisons of measured and calculated detector signals are shown in Figures 3o6.3 through 3.6.11. The I figures are in sets of three and are representative of various core exposures including coastdown conditions. For each statepoint the first figure of the set I presents the differences between the measured and predicted signal integrals for all instrumented locations. The instrumented core locations are indicated with circles in each of the figures. The I second and third figures of each set present axial comparisons in two specific instrumented core locations. The measurements are shown as a solid continuous line over 61 axial levels. The predicted reaction rates are represented as open circles. The two core locations were chosen as typical of regions on ., the interior of the core and on the core periphery. I In all comparisons, both the predicted and measured reaction rates have been normalized to a core average value of unity for each map. For purposes of quantifying comparisons, it is convenient to define the variable ORR (I,K,M) which I represents the difference between measured and calculated detector signals or reaction rates at location I,K, and map M.
- Thus, I DRR(I,K,M) = RRM(I,K,M) - RRC(I,K,M)
Where I I = K = M = Radial Detector Location Axial Detector Location Map Index I RRM = RRC = Measured Detector Reaction Rate Calculated Detector Reaction Rate I ,. NFU31/l 37 I 31 I
I NFU-0039 Revision 0 July 31, 1985 I J An average difference between measured and calculated reaction rates can be computed for each axial level as:
,,
I I D RR ( I , K , M) DRR( K) = I M
~~~~~~~~~~~~-
I I 1 I I M where the summation over I is performed for each available radial location, and M represents all flux map data except zero power maps. The mean observed J differences thus computed are the axial model bias and are listed in Table 3.6.4. since it is easier to describe the model uncertainties in terms of aeviations I relative to the observed bias, a second variable can be defined as
= ( RRM ( I , K, M) - RRC ( I , K , M) ) - b RR ( K )
I X(I,K,M) Where X(I,K,M) = The difference between *a, measured and calculated reaction rates adjusted for the observed bias. I All model power distribution behavior can now be characterized by quantifying the difference population X(I,K,M) or the integral of X(I,K,M); X(I,M). This latter quantity is the biased c ifference between the I measured and calculated detector signal integrals. To better evaluate the behavjor of the distributions of I X(I,K,M) and X(I,M), the difference population was divided into selected subgroups. The subgroups were chosen to parameterize the difference behavior as a function of axial height, reactor power level, and I cycle exposure. These subgroups were defined in a manner to exclude the axial points at grid locations and the upper and lower six axial points. I The axial regions are defined on Table 3.6.5. The difference population was evaluated for normality I using the chi-squared test. This test demonstrates I I NFU 31/l 38 32 I I
I I NFU-0039 Revision 0 July 31, 1985 I I that most of the subgroups cannot be considered normal. Typical comparisons of the difference population and a normal distribution is illustrated in I Figures 3.6.12 and 3.6.13. As indicated in Tables 3.6.6 and 3.6.7 and Figures I. 3.6.12 through 3.6.17, 95/95 confidence limits assuming normal statistics and 95/95 confidence limits based on non-parametric statistics are in good agreement. In I\ some cases the non-parametric limit is somewhat lower than the normal limit which simply indicates that the actual (not normal) distribution is slightly more peaked with fewer samples in the upper (higher M-C I values) tail of the distribution than is predicted by the normal distribution. To be conservative, 95/95 confidence limits were I evaluated using both normal and non-parametric statistics as described in Appendix A. limits (95/95) were computed for each subgroup. Confidence These I results are summarized on Tables 3.6.6 and 3.6.7 and Figures 3.6.14 to 3.6.18. Inspection of the figures show that the confidence limits are a function of axial height, reactor power level, and cycle exposure. I Generally, the confidence limits decrease with increased power and exposure. The approach taken to compute PSE&G model reliability factors, was to bound I the computed confidence limits. reliability factors are: Thus, the model RFFQ = 0.10 p > .so I RFF ti H = 0~16
- 0. 0 8
- (0.12*P) p P>
<
- 50
*30 P<
I 0.09 - (P/30) .30 In order to assess the impact of possihle dependence among data samples on reliability factors for local I peaking factors, the effects of reducing the sample sizes by two thirds (2/3) were evaluated. It was found that the reliability factors are relatively insensitive I to this reduction in sample 'size resulting in an increase1 rel iab il i ty factors of approximately
- 0 02 units for X(I,K,M) and .004 for X(I,M). As can be seen in Figures 3.6.14 through 3.6.18, the PSE&G reliability I factors remain bounding and are therefore not impacted significantly by possible dependence arnonq data samples.
,, I 33
I NFU-0039 Revision 0 July 31, 1"985 I FIGURE 3.6.l SALEM UNIT l AND SALEM UNIT 2 I MOVABLE INCORE DETECTOR LOCATIONS I I R p N M L K J H G F E D c B A I 28 15 l ti 4 3 51 2 1**: 10 30 39 52 3 5 36 43 4 I' 11 38 31 24 5 I 17 54 44 14 32 6 16 47 8 6 7 I
.
23 58 29 46 48 50 49 34 8 I 57 . 22 9 56 9 I 4 l 12 10 33 40 26 21 13 11 I 41 55 7 12 I 45 35 20 25 13 I 14 18 53 27 2 42 37 15
*I I
34 I
I I FIGURE 3.6.2 AXIAL LOCATIONS OF GRIDS AND DETECTORS NFU-0039 Revision 0 July 31, 1985 I I T 21!l.5<:11!l 21!l.661!l I 21!l.551!l I I ,, I-Zli i:i:i
- r:
AXIAL GRID 2111.55111 LlJ a: AXIAL DETECTOR LOCATIONS 0 u SIGNAL LOCATIONS LlJ
~
I-u
<C en I LlJ
- r:
~
....
- !;
21!l.551!l I I I 21!l.551!l I ALL ME~SUREMENTS I 24.43QI ARE IN INCHES I 1.-----v 1.243 BOTTOM OF FUEL ROD I 35
I NFU-0039 Revi~ion 0 July 31, 1985 .I I TABLE 3.6.1 I REACTOR STATE POINTS I SALEM 1 CYCLE 2 I MAP NO. CYCLE EXPOSURE ( MWD/MTU) POWER LEVEL ( %) D BANK POSITION (STEPS) 11 o.o 174 188 190 2160 3097 0 100.0 100.0 228 228 228 I:'. 194 4382 100.0 228 196 198 1201 6250 7275 7945 100.0 82.0 67.0 225 206 218 I SALEM 1 CYCLE 3 I MAP NO. CYCLE EXPOSURE ( MWD/MTU) POWER LEVEL ( %) D BANK POSITION (STEPS) I 1300 1313 1315 0 500 1040 0 100 212 228 I 100 222 1324 1330 1333 3165 4100 5670 99.5 97.0 97.0 22"3 220 228 I 1338 7060 96.8 228 1342 8800 75.0 202 I I I t NFU31/l 41 36 I
.,
I I NFU-0039 Revision 0 July 31, 1985 I I TABLE 3.6.2 I REACTOR STATE POINTS I SALEM 1 CYCLE 4 ,,, MAP NO. CYCLE EXPOSURE PmvER LEVEL 0 BANK POSITION ( WtID/MTU) ( %) (STEPS) o.o I 1400 1408 1411 180 560 0 84.0 100.0 211 228 228 1412 1580 100.0 225 I 1413 1414 2589 3715 98.6 100.0 228 215 1416 3836 100.0 228 I 1417 4998 100.0 228 ,, SALEM 1 CYCLE 5 MAP NO. CYCLE EXPOSURE POWER LEVEL 0 BANK POSITION I 1500 ( MWD/MTU) 0 ( %) o.o (STEPS) 216 1503 25 47.3 228 I 1507 1509 140 1391 99.3 100.0 228 228 1512 2531 99.9 226 I 1517 4662 99.5 228 ,, 1520 1522 1524 5444 7185 8923 100.0 99.9 100.0 218 228 228 I I I NFU 31/1 42 37 I'
I NFU-0039 Revision 0 July 31, 1985 I I TABLE 3.6.3 I REACTOR STATE POINTS I SALEM 2 CYCLE 1 I MAP NO. CYCLE EXPOSURE ( MWD/MTU) POWER LEVEL ( %) D BANK POSITION (STEPS) 'I. 2004 0 o.o 2102 2115 2435 4677 100.0 96.7 206 222 228 I\ 2120 7386 99.8 228 2122 2127 2129 9196 11755 13357 100.0 82.2 82.0 224 21*9 220 I 2131 14192 82.8 228 2133 15403 82.5 219 I CYCLE 2 SALEM 2
*I MAP NO. CYCLE EXPOSURE POWER LEVEL D BANK POSITION 2201
( MwD/MTU) 0 ( %) o.o (STEPS) 220 I 2203 21 48.6 2205 2209 47 292 72.1 98.4 180 214 228 I 2210 564 99.0 228 2213 2214 2217 1120 2106 3195 99.0 99.1 99.2 228 228 228 I I I I Nl:-,U 31/1 43 38 I I
I I NF'U-0039 Revision O July 31, 1985 I FIGURE Measured and Calculated Integrated Detector Responses 3.6.3 I SALEM 1 CYCLE 4 MAP 1411 Absolute Differences I Power = 100.0% Exposure 560 MWD/MTU I R p N M L K J H G F E D c B A .I /
,4.1/ '
Abs. Dif f = (Meas.-Calc)*lOO l v , z*9...1' I I '
\!. 8..1 -0.1'
./
2
' '
/' /
""0.6' v '\
I ,
-2.3
'- ,)
-0.5
\.. '
/ I\.. ,)
3.3 r'\. ./ 3
'
/
.t*4 .J ~1. 5/
' /
-2.l' 4 I ' /
v ... / '\ / v
- 1. 5
[\.,_
-2.7 r'\. /
- 1. 8
/
"\
~
- 1. 2 '
..I 5
/ /
'
/
'
' I
' '
/' / "\
I 4.6 [\.,_ ..I
-O.l
\..
2.9 r\. ..I
-1.1 1. 7
./
6
/ ,)
' v '
/
-1.2
'\
/
/
-1.8 ' -1. 7
\.. /
/
-1. 3
/
7
' '
,)
v "/ v
'v ' / v "\/
3.2'
,/
-2.3 ' ,) I\..
2.1 ' -2.7
'\
/
-1. 3
'\ /
-1.8 '-2. 7 -0.5
./\..
8
/ ,) /'\. ,)
' v 3.3
'\ ' ' v
-4.2 v
... 4.0
'\ ' /
2.6 ' 9 r--.. r\. /
/
' '
,) ,)
v /
' '
/
-1.4 ' -3.3 2.1 10 r'\. ,)
'- / I\.. /
v '\ v
' v /'
'
4.4 -2.3 ... 2. 6' 2.9 11 I \.. ,/ I\.. /
/
'
-2.7 v-1.3
'
r\. / \... / 12
, "- ,) I\.. ,)
I ~ 0.4 ' / v
~
"
-1. 2
/
r
-3.8
"- ,/
' ./
~.3,) ' 13
/
" v v 2.9 "
/
' "
I '
,)
v ' 0.4 I\.. / I\..
-2 .6
./ I\..
-0.l
/
14 3.5
/
1.1 ' 15 I ~ ./ ~ ./ 39 I
I Nf'U-00 39 Revision 0 July 31, 1~85 I FIGURE 3.6.4 I MEASURED AND CALCULATED DETECTOR RESPONSES SALEM 1 CYCLE 4 MAP1411 I 2.5-P--~~-~-~-~-.---~-~-------- I Legend I ll-t.8.£: Lt> 0 MEASURED PREDICTED 2 ..............;;,,;;;;m;;;;;-....,;;;;_ _ _ _ ___._ ............ ****- *************** * * * * * * * * * * - - - - 1-1 LU (/) z I I I I I 0 c.. (/) LU a:: a:: 1.5 ........ *.......... ********:-********. ********:********* *******-:-********:********:********* *******-***************** I
I I I I
I
I I I I
I I I I
I I I I
I I
,,
g ~l u LU
~
LU I c . I
.. ..
,,
LU 1 ......**. **** .* * **...*~ ***.*... ~ .**..*.* ! ........ ~-*-*****.:. ........* ......... ........* ..... -:-*****************
>
~
. ..
.....J LU a::
0.5 .. ...... ........................... ........ ........ *******-******** ................... *******-******* I 0 I I o_....,__..__,..__...,..__.,.__,....._,.._...,..__...---..--...--...... 0 5 10 15 20 25 30 AXIAL POINTS 35 40 45 50 55
--..-~
60 65 I, 40
. l_
,,
I
---------*------
- I FIGURE 3.6.5 NFU-0039 Revision 0 July 31, 1985
~ MEASURED AND CALCULATED DETECTOR RESPONSES SALEM 1 CYCLE 4 MAP1411 ~ 2.5-----------~---~---------- ~I Legend MEASURED 2 ***** *** 11-taE =N13 ......... *******-******* ......... ********----
0 PREDICTED w
~
Vl z 0 0-
~ 1.5 ....*.*.........*................**......**.**....* *******-******* *................* *******-******* *********
a:::
- a:::
0 t>
~ ~
w Cl w 1 ........ **************
>
~
_J w a:::
' .
0 *5 . ...... .********* ....... **************************~********-******* ******************41 *******-***
- o 0 .......--.,....._..,--~___,,___;..---;.--......--.;-~~--;.--..;---.;---t 0 5 10 15 20 25 30 35 40 45 50 55 60 65 AXIAL POINTS 4 J_
NFU-0039 I Revision 0 July 31, , 1985 I FIGURE 3.6.6 Measured and Calculated Integrated Detector Responses I SALEM 1 CYCLE 5 MAP 1522 I Absolute Differences Exposure Power = 99. 7% 7185 MWD/MTU
,,
R p *N M L K J H G F E D c B A I
/'
0.3
\,,
',)
/
-0.0
\,,
'
)
Abs. Di ff = (Meas.-Calc)*lOO 1 I v ' I' 0.8
'- ..I
'-0.~ ' 2 I
v I'
-1. 7' /'
-1.8
"\
0.6 " 3 v ,,,, ...,_ v
,)
'- .J ...,_ ,)
I> 0.1'
'-
- 1. 7
.J I\.. ..I 0.6
'- /
'v 4
/'
'-
0.9
,) ~-6 ..I
' /
'9*1/ ' 5 I
/' I\..
'
- 1. 2
/
/'
-0.0
'-
'
.J 0.9 ')
/'
2.2
\,,
'
)
6 I v I' '
'
/
-0.4 0.6
'1.~
" /'
0.3 ' 7
.J ...I I
\,, \,, \,, ,)
/ '- 1.3'
/
/'
'
-1. 7,,
\,,
/'
\,,
-0.1 ',)
/'
'
'1. 7,;
... 3.0' 0.7 '
/'
f\...
/'
,) '- .J 8
I
/'
' '
/'
- 1. 2 9
~0*5;
'- ..I v v ....
/'
1.5' -1.5 " 0.8 10 v I\. ,) I\.. ,J
/'
'- ..I
/'
'
I '- 1.1
"\
(2*~ ' /' 0.9 '
.J
\,, 0. 5.J 11 I
,) \,,
/' v v v 0.7
\,,
'.J
-0.9 I\.. '
./
. 12
-1. 9
"\
0.6 I\..
' /'
-2.3'/
/
~-1)
'\
13 I
'- ./ ,,)
'-
v v
'
"\ /' /*,'
' 14 0.7
'- / '-
-1.5
,)
/'
'
-1.9 I\.. / !\,.
0.3
/
15 I
~-8./
I 42 I
I I NFU-0039 Revision O July 31, 1985 I FIGURE 3.6.7 I i MEASURED AND CALCULATED DETECTOR RESPONSES SALEM 1 CYCLE 5 MAP1522
~ 2.5--~--~-~-~-~---'!"--"!"'"-~-~~----'!"--.
I I POWER= 99.73 Legend EXPOSlff: = 7185 MWD;MTU MEASURED
.. . . ... . lHM3L£ =LD .
0 PREDICTED
..... ....
2 ********~******** :-********~*****************-----
... ....
... ...
.. ...
w . (/) z
..... .....
.0 a_
..... ...
...
. .
~ . 1.5 ........ **************** ....... ................
-~
0
~-.......
o
....... *: ................. -:*. ............... ........ .
.
I o
-~
o 0::: . o I
. o o
. I o
.
o o
. o I
0::: .. 0 0
..
I 0
..
o f
.. I 0
..
o I 0 ... .. ... ... ...
...
o o 0 0 I t5 .. .. .. .. w .. 0 I
. o
..I 1-w Cl w
>
~
......J w
0::: 0 5 ********* ........ ; .................. ~ .....................................................................................
. ..:.. : ...:. .:... ..:. :o
....
.... ... .
... ..... ..
'
.. .. .
.. ..
...
~
.. o
..o
.. o
.. I o
.... .... .... .. ..
0 I 0 I 0
.. .
... ... ... . ....
o..,.-.___.____.,,___.___..____. __,..____ ... o o I o I ___ .
. _ _. . ____,..___.
I'
~*---i-0 5 10 1S 20 2S 30 35 40 4S so SS 60 6S AXIAL POINTS I 43 I
*------*I I
i I NFU-003~ Revision 0 July 31, 1985 I I FIGURE 3.6.8 I I
,,
MEASURED AND CALCULATED DETECTOR RESPONSES SALEM 1 CYCLE 5 MAP1522 2.5-...-~~-~-~-~-~--~------~--. POWER= 99.73 I Legend EXPOSlR: = 7185 MWD/MlU 2 ........ 11-M31...E: N13 0
***************** ......................... ~--........
MEASURED PREDICTED I I Lt.I I (/) z 0 a_ (/) Lt.I 0:: 1.5 ................ -:* ................. *: ........
- I o
- I o
........ ................... ................ -:* .........
~- I o
-~-
- I o
- o
;. ....... ........ .
-~
- I o
I
.:... .:.. :.. :.. .:... :..
I 0
.
0:: .. .. 0 . . '
~ ~ ~ I I-(.)
Lt.I 1 1 1 I- .. I
.
0
..
o
..
0
..I
..
o I o I o o I O Lt.I I I o I o 0 6 . .
*
- I
- Lt.I
>
~
....J Lt.I 0::
.1
* ............... T........ ........ r................ T........ ........ r................. r....... T..... *r ...... .
0 5 I o
~
I
~
0 I O I I I I o o I o o I o I o o I I
.
o o
.
o 0
. o o
.
I I
.
o o I I
'
- o I O O O O I
.... .... ... ... .....
o o o o I I o
.. o o o o I
.. . .
I o o I o o I
.0 o
o o
. o o
I 0
.0 o
o o
.
o o o o
.I o
o I
. I I
I 0
..
o o
..I o
..
o o
..
o o
..
o o
..
I I o--.__..._..,_..___...,__ I O __..._. O O O I ___, 0 5 10 15
~*--__..
20 25 30
--..;.---~*--~--~*--~*i--__,.
35 40 45 50 55 60 65 I AXIAL POINTS I 44 I
I NFU-0039 I Revision 0 July 31, 1985 FIGURE 3.6.9 I Measured and Calculated Integrated Detector Responses I SALEM 2 CYCLE 1 MAP 2133 Absolute Differences I Power 82.5% Exposure = 15403 MWD/MTU I R p N M L K J H G F E D c B A I Abs. Diff = (Meas.-Calc)*lOO 1 I ~
\....
0.5
'
/
/
0.1
\....
'
/
2
/
-5.5
' v-1. 3' /'
0.8
' 3 I '- /
/ '
~ ,) i\.. /
4 1.4 I /
"-
2.2
/
'
'- .I I'
- 1. 8
\.... /
'
5 v ' I '\ ' I -1.4
\.... /
~L9)
/
"-
-0.8
./
v
'-
'
0.7
./
6 v-0.0' /'
-1.5' I' 2.2 ' 7
/
-1. 5
' /
-3.0
'
i\.. ,)
/
'-1.5
'
\..
/'
\..
0.8
'/
,/
/
~
0.1
,,-2.2'
'-
/'\.
.J
./
8
'- .I '- / i\.. /
I 9
/ ' v '
10 I -0.7
'- /
/ ' /'2.1 '
",
3.2
/
,r ' 11
-1.1 1. 2 I '- .I
/
-0 .1)
'
\.... /
/ '
.... 0 5/
['-,_.
"' .I 12
'*...... _/
I / ' 13
~.7 /
v0.1' / ' v 2.4' I I'.. /
/ '
- 1. 9
~ /
14 15 0.2 I \,,. / 45 I
I
~
i NFU-0039 Revision Q. July 31, 1985 I FIGURE 3.6.10 I MEASURED AND CALCULATED DETECTOR RESPONSES SALEM 2 CYCLE 1 MAP2133 I I I 2 ........
- -"5403 1}"'8.E =lJ)
~]
UI
. .......* *******-******** ........ ********-----
0 Legend MEASURED PREDICTED
,,
I w (/) z 0 a.. I (/) w 1.5
.... ..
ee * * * * * * * * * * * * * * - : * * * * * *
- e ~-* D** ******* e
.
*<>-:-Cl** *** e *fl>********
I
********* e 0 e . . . . ** * * * * - * * ** e ** ********
- e ** * * * * ** e * ***
a:: ... .. .... a:: .. .. ... 8 f-w
. I 0
.w
> I ~ -'
w a:: I 0 I
.....
0.5 ......... ...
*******-:*********~***************** ......... *******-***************** ********* *******-********:********
. :o
....
... ..... I
..... ...
...
...
...
..
I
.
0-'---;...-....;.*--~*=--*~*--.;..---;-.~-- .....--'9i---ti---i---..;.---t 0 5 10 15 20 25 30 AXIAL POINTS 35 40 45 50 55 60 65 I I 46 I
'I I
NFU-0039 Revision O July 31, 1985 FIGURE 3.6.11 MEASURED AND CALCULATED DETECTOR RESPONSES SALEM 2 CYCLE 1 MAP2133
2.5-------~-~-~----i----.---"--i'--. POWER= 82.5% Legend
' 2 -------- 11-Ml.E =N13 --------* ---------------* 0
**--**-- ________ ...______
MEASURED PREDICTED I' w (/) z 0 a.. (/) w 1.5 ........................... **************************- ............................ ********~********-******* ........ . a:: a:: 0 t> w I-
*I w
Cl w
>
~
_. w
~
a::
. o .......................................................
05
. r.............................................
. . ,.. ..
' :o
- 0-i---.,......_,~_,.--_,_
0 5 10 15 20 __'"'-"__,__~-.......---;---.---P....~--.... 25 30 35 40 45 50 55 60 65
- ~L-.--+-------47 AXIAL POINTS
_ _ _ _ __
I
.NFU -0 l) 3 9 Revision 0 July 31, 1985 I I
TABLE 3.6.4 MEAN OBSERVED DIFFERENCES AXIAL MODEL BIAS I AXIAL LEVEL MEAN DIFF AXIAL LEVEL MEAN DIFF I I X( I) I X (I) (TOP) 61 60
.119
-.024 31 30
-.012
-.014 I
59 -.023 29 58 57
.005
.018 28 27
-.016
-.062
-.039 I
56 .015 26 -.015 55 54 53
.018
-.031
-.020 25 24 23
.007
.010
.009 I
52 -.021 22 .010 51 50
-.009
.003 21 20
.031
-.041 I
49 -.001 l~ -.033 48 47 46
-.005*
.011
-.016 18 17 16
-.018
-.009
-.003 I
45 -.040 15 -.001 44 43
-.024
-
- 011.J 14 13
-.001
- l) 12 I
42 -.017 12 .029 41 40 39
-.023
-.033
-.033 11 10 9
-.024
.015
.035 I
38 -.004 8 .052 37 36
-.069
-.036 7
6
.066
.068 I
35 -.009 5 .064 34 33 32
-.001
.005
.007 4
3 2
.056
.038
.024 I
1 .016 I I I NFU31/l 53 48 I I
I I NFU-0039 Revision 0 July 31, 1~85 I I TABLE 3.6.5 I AXIAL REGION DEFINITIONS I REGION AXIAL POINTS l I 2 7 - 10 14 - 18 I 3 4 22 - 31 - 27 35 I 5 39 - 44 6 48 - 52 I I I ,I I I I I I NFU 31/1 54 49 I
-------~~-- .. -..o .. ,- . .
l I FIGURE 3.6.12 I DISTRIBUTION OF ERRORS X(i,k,m) 0.5 NON PARAMETRIC STATISTICS 95/95 CONFIDENCE LIMIT I I OBSERVED 0.4 DISTRIBUTION
/~
I NORMAL STATISTICS I I 95/95 CONFIDENCE LIMIT I
>-
() I z 0.3 w
- > PSEG RF FQ Ul 0 0.
w a:::
\
LL. w
>
~ 0.2
\
~
w a:::
\
~ :Al z 0.1 c (1) '"IJ NORMAL f-' < c DISTRIBUTION '< ~- I en o w ~-0 f-. 0 w
.. ::i \.!)
o.oJ--=.....,~~===-~~--~~--~-----~-_t_i______c-=~~~---=~
-4 -3 -2 -1 0 1 2 3 4 STANDARD ERROR UNITS (Z)
---------
~-----*I FIGURE 3.6.13 DISTRIBUTION OF ERRORS FOR INTEGRAL X(i,m) 0.5 NON PARAMETRIC STATISTICS 95/95 CONFIDENCE LIMIT 0.4 OBSERVED DISTRIBUTION I
,,--"' . NORMAL STATISTICS 95/95 CONFIDENCE LIMIT
>-
() I zl..&J 0.3
- ..11 :::>
\
I-' 0 PSEG RF FAH l..&J e:::: La... I w
>
~
_J 0.2 I \ l..&J e:::: I \ 0.1 NORMAL I c..i ::0 c: ro I-' <: c z
"1 DISTRIBUTION I ~I-'*
Wl-'-0 rn o I I I-'
~
I-' 0 0 w
- J l.C' l.O
---
O'.l U1 0.0
-4 -3 -2 -1 0 1 *2 3 4 STANDARD ERROR UNITS (Z)
I NFU-00 39 Revision 0 July 31, 1985 I I TABLE 3.6.6 CONFIDENCE LIMITS FOR X(I ,K,M) DISTRIB.JTION I BY SJBGRaJPS I REACTOR CYCLE AXIAL i'UMBER ST. DEV 95/95 CONFIDENCE LIMITS POWER (%) 0 EXP0SJ RE ( G/T) ALL R.EX;ICNS 1-6 SAMPLES 10075 .075 NORMAL NON - PARAMETERIC I
.125 .139
,
50< p <70 I ALL 1-6 8059 .045 .076 .063 100 ALL 1-6 49573 .036 .059 .063 100 E<2.5 1-6 22966 .041 .067 .072 .. , 100 2.5<E<6 1-6 19105 .031 .051 .052 100 6(E 1-6 7502 .028 .047 .040 100 100 100 ALL ALL ALL 1 2 3 6396 7991 9587
.036
.033
.034
.061
.055
.057
.069
.065
.069 I!
100 ALL 4 7995 .035 .059 .067 100 100 ALL ALL 5 6 9595 8009
.037
.039
.062
.066
.079
.073 I
I I I I
,,
I NFU 31/l 57 52 I. I I
I I NFU-0039 Revision 0 July 31, 1985 I I TABLE 3.6.7 CONFIDENCE LIMITS FOR X(I,M) DISTRIBUTION I BY SUBGROJP I REACTOR POWER (%) CYCLE EXPOSURE (G/T) NUMBER SAMPLES ST. DEV 95/95 CONFIDENCE LIMITS NORMAL NON - PARAMETERI 1 0 I 50< p <70 100 ALL ALL ALL 322 258 1593
.045
.034
.028
.081
.062
.048
.075
.056
.055 I 100 100 100 E<2.5 2.5<E<6 6(E
. 739 614
.033
.CJ24
.057
.042
- 066
.042 240 .021 .038 .045 I
I ! NFU 31/1 59 I I I ,I 53 NFU 31/1 58 I I
FIGURE 3.6.14 CONFIDENCE LIMITS FOR X(i,k,m) VS REACTOR POWER % 0.20 0.18 0.16 0.14
~
(Jl ~
~
_J 0.12
~' '
~ ~,_
PSEG RffQ w (.) zw 0.10 NORMAL STATISTICS ~--- 95/95 CONrlDENCE LEVEL '~ ::---.__ 0 G: z 0 (.) 0.08
----- -------
0.06 ~ ' ' - - - - ----- -==-==--- L NON-PARAMETRIC STATISTICS 95/95 CONrlDENCE LEVEL
~ :;o :z
. 0.04 C: .CD t'IJ I-' <: c:
. '< I-'* I Ul 0 Wl-'*0
'1-'QW
~ ::l l.D 0.02 0.004----------~-------~--
-.
o 25 50 75 100 .. REACTOR POWER %
FIGURE 3.6.15 CONFIDENCE LIMITS FOR X(i,k,m) VS CYCLE EXPOSURE 0.14 0.12 PSEG RFFQ 0.10-r----------------------------------- (/) I-Ul
~ NON-PARAMETRIC STATISTICS
-----r-Ul
__..J 0.08 95/95 CONFIDENCE LEVEL -=- w u z w 0 LL z 0.06 lL NORMAL STATISTICS 0 u - - - - - - __ *l ---95/95 CONFIDENCE LEVEL---. --- 1 0.04
~ :;cl z cf-' (])<: c'""
'< I-'* I Cfl 0 Wt-'*C 0.02 f-' 0 w
.. :J \!)
.._. 0 l.O cc lJ1 0.00-+------~- --,
o 2 4 6 8 10 12 CYCLE EXPOSURE (GWD/MTU)
I,. FIGURE 3.6.16 CONFIDENCE LIMITS FOR X(i,k,m) VS AXIAL HEIGHT I I 0.14 j I I 0.12
!
PSEG Rf"FQ 0.10-1------------------------------------
~
Ul
~
O"\ _J w 0.08 NON-PARAMETRIC STATISTICS 95/95 CONFIDENCE LEVEL ~ 7 -- - ...... ......
-
-----
/
---
"
--- ----
(.)
- ------
/
zw 0 G: 0.06
- --- --- - --- ......
z 0 (.)
..........___ --N-ORMAL STAT-15-TIC_S_ - + .
95/95 CONF"IDENCE LEVEL_j 0.04
~ ::0 z i:: <D t"JJ I-' <: c
'<.I-'* I en o Wl-'*0 I-' 0 VJ 0.02 ' ::J l.D f-'0 l.C CXl Ul o.oo---~~.--~-.--~~-.-~~..--~--..~~-.-~~-.--~___,~~-.-~~---~~~~~
0 1 2 3 4 5 6 7 8 9 10 11 12 AXIAL HEIGHT (FEET)
- .. - ----
FIGURE 3.6.17 CONFIDENCE LIMITS FOR X(i,m) VS REACTOR POWER % 0.12 0.10 PSEG RF F"AH
~
0.08
- - ---- ---- NORMAL STATISTICS
\_ ---- 95/95 CONFIDENCE ~
.
~ ::::E -...J ::J - - _ - - ---- ---- .
---=-=----~
w (.) zw 0.06 0 u.... NON-PARAMETRIC STATISTICS - - -:- - .- ~ -----* - 95/95 CONFIDENCE LEVEL - -- z 0 (.) 0.04 c:... ::0 z c:: ro r-r; I-' <: c: 0.02 '< I-'- I en c Wl-'*0 I-' 0 w
' ::i ~
o.oo--~~~~~~~~~~~~~~~~.--~~~~~~~-T-~~~~~~~---, 0 25 50 75 100 PERCENT REACTOR POWER%
FIGURE 3.6.18 CONFIDENCE LIMITS FOR X(i,m) VS CYCLE EXPOSURE 0.12 0.10 PSEG RF" F"AH 0.08
~
Ul ~ 00 - _J w - - - - u 0.06 z w 0 G: z 0 1I NON-PARAMETRIC STATISTICS\ 95/95 CONF"IDENCE LEVEL
- u 0.04 NORMAL STATISTICS "
95/95 CONF"IDENCE LEVEL_/
~ :;c z c: C1l 'Tl f-'<:C 0.02 '< I-'* I (fl 0 Wt-'-0 t-' 0 w
~ :J \.0 t-' 0
\D cc Ul o.oo---~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
0 2 4 6 8 10 12 CYCLE EXPOSURE (G/T)
I NFU-0039 I Revision 0 July 31, 1985 I I 3.7 Verifiication of Transient Power Distribution Simulation Capability I I (To be completed later) I I I I I 1* I I I I I NFU31/l 61 59 I I
I I NFU-0039 Revision 0 July 31, 1985 I I ,, 4.0 References
- 1. Advanced Recycle Methodology Program (ARMP)
System Documentation CCM-3 Research Project 118-1, September 1977. I 2. Pfeifer, C. J., "PDQ-7 Reference Manual II", WAPD-TM-947(L), Westinghouse Electric Corporation, February 1971. I 3. Breen, R. J., o. J. Marlowe, and c. J. Pfeifer, "HARMONY: System for Nuclear Reactor Depletion Computation," WPAD-TM-478, Westinghouse Electric I Corporation, January 1965.
- 4. Walpole, R. E., Myers, R. H., "Probability and I Statistics for Engineers and Scientists, MacMillan Publishing Company, New York, 1978.
- 5. Owen, D. B., "Factors for One-Sided Tolerance Limits and I for Variables Sampling Plans", SCR-607, Sandia Corporation, March 1963. (Available from office of Technical services, Department of Commerce, Washington I 6.
D.C.) USNRC Regulatory Guide 1.126, "An Acceptable Model and Related Statistical Methods for the Analysis of Fuel I Densification.", March 1978.
- 7. Somerville, P. N., "Tables for Obtaining Non-Parametric I Tolerance Limits", Annals of Mathematical Statistics 29, 599 (1958).
- 8. Assessment of the Assumption of Normality (Employing I Individual Observed Values), ANSI Nl5.15-1974.
- 9. Safety Evaluation of the PSE&G Rod Exchange Methodology, I NFU-004, Revision 2, August 22, 1984.
I I I 60 I NFU 31/l 61
I I NFU-0039 Revision 0 July 31, 1985 I I APPEND IX A I STATISTICAL METHODS FOR THE DETERMINATION AND APPLICATION OF UNCERTAINTIES I I I I I I I I I I I I I NFU 31/1 62 A I
I NFU-0039 Revision O July 31, 1985 I I APPENDIX A I STATISTICAL METHODS FOR THE DETERMINATION AND APPLICATION OF UNCERTAINTIES I The purpose of using statistical methods is to compute the value X such that there is a 95% probability at I the 95% confidence level that XR will be conservative with respect to X (true value) when applying the calculational met~ods to safety related reactor analyses. I The first step is to determine whether or not a distribution is normal. in Section A.l are used. If it is, the methods described If the distribution cannot be I treated as normal, but the distributions are known, then the methods described in Section A.2 are used. I I I I I I I I I NFU 31/1 63 Al I
I I NFU-0039 Revision 0 July 31, 1985 I I A.l Application of Normal Distribution Statistics Treatment of Measurement and Calculational I Uncertainties Comparison of measured and calculated reactor I parameters include the effects of both the measurement and calculational uncertainties. Methods used in this report to isolate the calculational uncertainties are described below in terms of the following definitions: I XT = true reactor parameter I XM = measured reactor parameter XC = calculated reactor parameter I eM = XM XT = measurement error ec = XC XT = calculation error I eMC = XM - Xe = observed differences I µi = ei =mean error (i = M, C, or MC) I I If eM and ec are independent, then the I following relationships exist. (Note that these relationships apply for non-normal distributions as I well)
- I I
I I
I NFU-0039 Revision .0 July 31, 1985 I I 2 2 2
µ aMC =
µ a
c
+
µ aM (A-1)
I MC = M c (A-2) These equations can be solved for a~ andµc. ar and µC are calculated from hitor1cal data Once I tney can be used to apply conservatism to future calculations of reactor parameters, XR, as follows: I xR = xc - Mc +/- Kc ac The factor Kc is defined as described in Table I A.l to proviae a 95% probability at the 95% confidence level that XR is conservative with respect to the true value, XT. I Alternatively, as done in most instances this report, it can be noted that since each term in equation (A-1) is greater than or equal to zero each I term is bounded by the variance in the observed differences. Thus the calculational uncertainty can be conservatively estimated as the uncertainty in I the observed differences between measured and calculated values: a 2
=a MC 2
(A-3) I or c ac - aMC (A-4) I In the later alternative, once crMc and µMC are calculated from historical data they can be used I to apply conervatism to future calculations of reactor parameters, XR as follows: I where RF = Kc aMC (A-5) I The quantity µMC is used as a bias on the calculated parameter and K is as defined I above. The term RF is called the reliability factor as described below. I I NFU31/l 65 A3 I
I I NFU-0039 Revision 0 July 31, 1985 I I Reliability Factors It is the objective to define reliability factors I which are to be used to increase or decrease calculated results to the point where there is a 95% probability at the 95% confidence level that they I are conservative with respect to actual parameters. For any given application, we are only concerned with one side of the component; that is, if the I calculated value is too large or too small. we may therefore use one-sided tolerance limits based on normal distributions to find a Kc which will give a I 95% probability at the 95% confidence level to the reliability factor defined by RF = Kc* 0c I Numerical values of Kc for various sample sizes used to calculated crc are provided on Table A.l I (Reference 5). I I I I I I I I NFU31/l 66 A4 I
I NFU-0039 Revision 0 July 31, 1985 I I TABLE A. l I SINGLE-SID ED TOLERANCE FACTORS N Kc I
- -
2 3 26.26 7.66 I 4 5.14 5 6 7 4.20 3.71 3.40 I 8 3.19 10 9 3.03 2.91 I 11 2.82 12 15 20 2.74 2.57 2.40 I 25 30 40 2.29 2.22 2.13 I 60 2.02 100 200 500
- 1. 93 1.84
- 1. 76 I
1.645 I OJ N = Number of samples used to calculate oc I I I I I NFU 31/l 67 AS I I
I NFU-0039 Revision 0 July 31, 1985 I A.2 Application of Non-Normal Distribution Statistics This section documents the procedure used to I determine the value XR such that there is a 95% confidence level that XR will be conservative relative to the actual value (X ) when the distribution of X is not assume~ to be a normal I distribution. The approach taken is consistent with non-parametric methods given in reference 6 and 7. I In general, the procedure requires the ordering of N samples taken from a continuous but unknown distri~ution function. The statistic "rn" is determined such that, at the 95% confidence level, I 95% of the population lies between the rth smallest and the sth largest value in the ordered N samples, where m = r + s. The statistic m can be determined I from Table A.2 (Reference 7). Since, for any given application we are usually only concerned with one side of the component, one-sided tolerance limits are required. Therefore, for upper one-sided I tolerance limits, r is set to zero, and m = s. This procedure has been implemented to obtain I reliability factors using th~ following steps. First, the mean error µMC = eMc was determined, where eMc = XM - Xe. (See Section A.l for definitions). Next, the population I of N errors eM - µ was computed, and the resulting dis~ribu~~on ordered. Using Table A.2, the mth value of the error distribution defines I error eR,which, at the 95% confidence level, 95% of the error distribution will be less than eR. I I I I I NFU31/l 68 A6 I
NFU-0039 Revision o July 31, 1985 I I once e and µMC are calculated from histor~cal data they can be used to apply I conservatism to future calculations of the rector parameter, XR, as follows: (A-6) I XR = XC + µ MC .:!:. RF where RF = eR. I The term RF is the reliability factor which provides the desired 95% probability at the 95% confidence level for the com'puted parameter X. I I I I I I I I I I I I NFU 31/1 69 A7 I I
I NFU-0039 I Revision 0 July 31, 1985 I I TABLE A.2 Values of m for 95% Confidence and 95% I Probability Tolerance Limits Number of I Observations (n) m I 50 55 60 1 65 1 I 70 75 1 1 80 1 I 85 90 95 1 1 2 100 2 I 110 120 2 2 130 3 I 140 150 170 3 3 4 200 5 I 300 400 9 13 500 17 I 600 700 800 21 26 30 900 35 I 1000 39 For n> 1000 I Increase m by 4 for each additional 100 observations I I I NFU31/l 70 AB I
NfU -00 39 Revision U July 31, 1985
I NFU-0039 Revision 0 July 31, 1985 I I APPEND IX B COMPUTER CODE
SUMMARY
DESCRIPTION I COMPUTER I CODE DESCRIPTION CPM CPM is a multigroup two-dimensional collision probability code for I depletion and branch calculations for a single assembly. I Reference 1 EPRI-CELL EPRI-CELL computes the space, energy and burnup dependence of the neutron I spectrum within cylindrical cells of light water reactor fuel rods, used to generate cross sections for It is I PDQ on a ECD AT A f i le
- Reference 1 I INTEGRAL INTEGRAL edits PDQ files to obtain pin and assembly powers. Pin to box ratios are then input to TRINODE.
NU PUNCHER NUPUNCHER prepares HARMONY cross section tables from cross section data on an ECTIATA file. I Reference 1 I PDQ7/ PDQ7/HARMONY is a nuclear reactor HARMONY analysis program which solves the neutron diffusion equations and performs depletion calculations. I Reference 2, 3 I I I NFU 31/l 72 Bl I I
1-I NFU-0039 Revision 0 July 31, 1985 I I COMPUTER CODE DESCRIPTION I SHUFFLE SHUFFLE is the same as EPRI-SHUFFLE and will read a PDQ7 concentration file and write a new updated I- concentration file. It is used to simulate assembly movement between cycles. I Reference 1 SIGMA SIGMA calculates the predicted I detector reaction rates using nodal power distribution and PD07 detector reaction rate to assembly power factors. The predicted detector I reaction rates are then compared to measured detector reaction rates. I TRINODE TRINODE is a modified version of the EPRI-NODE-P computer code program. Modifications are summarized as follows: I a) Automated file management b) User friendly input c) Rod search for constant axial I d) offset control Separate BP reactivity insertion equations e) Flexible edit options I TAU TAU is a computer code used to compute statistics from residual I reaction rates generated by SIGMA. I I I I . NFU 31/l 73 I B2}}