Regulatory Guide 1.190: Difference between revisions
StriderTol (talk | contribs) (Created page by program invented by StriderTol) |
StriderTol (talk | contribs) (Created page by program invented by StriderTol) |
||
(One intermediate revision by the same user not shown) | |||
Line 9: | Line 9: | ||
| docket = | | docket = | ||
| license number = | | license number = | ||
| contact person = Jones W | | contact person = Jones W (301)415-7558 | ||
| case reference number = -nr, DG-1025, DG-1053 | | case reference number = -nr, DG-1025, DG-1053 | ||
| document report number = RG-1.190 | | document report number = RG-1.190 | ||
Line 15: | Line 15: | ||
| page count = 56 | | page count = 56 | ||
}} | }} | ||
{{#Wiki_filter:Regulatory guides are issued to describe and make available to the public such information as methods acceptable to the NRC staff for implementing specific | {{#Wiki_filter:U.S. NUCLEAR REGULATORY COMMISSION March 2001 REGULATORY | ||
WWW.NRC.GOV> in the Reference Library under Regulatory Guides. | GUIDE | ||
OFFICE OF NUCLEAR REGULATORY RESEARCH | |||
REGULATORY GUIDE 1.190 | |||
(Previous drafts were DG-1053 and DG-1025) | |||
CALCULATIONAL AND DOSIMETRY METHODS | |||
FOR DETERMINING | |||
PRESSURE VESSEL NEUTRON FLUENCE | |||
Regulatory guides are issued to describe and make available to the public such information as methods acceptable to the NRC staff for implementing specific parts of the NRCs regulations, techniques used by the staff in evaluating specific problems or postulated accidents, and data needed by the NRC staff in its review of applications for permits and licenses. Regulatory guides are not substitutes for regulations, and compliance with them is not required. Methods and solutions different from those set out in the guides will be acceptable if they provide a basis for the findings requisite to the issuance or continuance of a permit or license by the Commission. | |||
This guide was issued after consideration of comments received from the public. Comments and suggestions for improvements in these guides are encouraged at all times, and guides will be revised, as appropriate, to accommodate comments and to reflect new information or experience. Written comments may be submitted to the Rules and Directives Branch, ADM, U.S. Nuclear Regulatory Commission, Washington, DC 20555-0001. | |||
Regulatory guides are issued in ten broad divisions: 1, Power Reactors; 2, Research and Test Reactors; 3, Fuels and Materials Facilities; 4, Environmental and Siting; | |||
5, Materials and Plant Protection; 6, Products; 7, Transportation; 8, Occupational Health; 9, Antitrust and Financial Review; and 10, General. | |||
Single copies of regulatory guides (which may be reproduced) may be obtained free of charge by writing the Distribution Services Section, U.S. Nuclear Regulatory Commission, Washington, DC 20555-0001, or by fax to (301)415-2289, or by email to DISTRIBUTION@NRC.GOV. Electronic copies of this guide are available on the internet at NRCs home page at <WWW.NRC.GOV> in the Reference Library under Regulatory Guides. This guide is also in the Electronic Reading Room at NRCs home page, along with other recently issued guides, Accession Number ML010890301. | |||
TABLE OF CONTENTS | |||
Page | |||
==A. INTRODUCTION== | ==A. INTRODUCTION== | ||
............................................................. | . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 | ||
==B. DISCUSSION== | ==B. DISCUSSION== | ||
................................................................ | . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 | ||
==C. REGULATORY POSITION== | ==C. REGULATORY POSITION== | ||
.................................................... | . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 | ||
1.NEUTRON FLUENCE CALCULATIONAL METHODS ............................. | 1. NEUTRON FLUENCE CALCULATIONAL METHODS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 | ||
1. | 1.1 Input Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 | ||
1. | 1.2 Core Neutron Source . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6 | ||
2. | 1.3 Fluence Calculation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7 | ||
2. | 1.4 Methodology Qualification and Uncertainty Estimates . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 | ||
3. | 2. NEUTRON FLUENCE MEASUREMENT METHODS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 | ||
3. | 2.1 Measurement Procedures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21 | ||
2.2 Validation in Standard and Reference Neutron Fields . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26 | |||
2.3 Fluence Determination from Detector Measurements . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 | |||
2.4 Ex-vessel Dosimetry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 | |||
3. REPORTING . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28 | |||
3.1 Fluence Methods . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28 | |||
3.2 Multigroup Fluences . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28 | |||
3.3 Integral Fluences . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28 | |||
3.4 Comparisons of Calculation and Measurement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29 | |||
3.5 Specific Activities and Average Reaction Rates . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29 | |||
==D. IMPLEMENTATION== | ==D. IMPLEMENTATION== | ||
.......................................................... | . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29 References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40 | ||
LIST OF TABLES | |||
1. Summary of Regulatory Positions on Calculation and Dosimetry . . . . . . . . . . . . . . . . . . . . . . . . 30 | |||
2. Threshold Detectors Recommended for Pressure Vessel Dosimetry . . . . . . . . . . . . . . . . . . . . . . 34 LIST OF FIGURES | |||
1. Discrete Ordinates Calculation Methodology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35 | |||
2. Monte Carlo Calculation Methodology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36 | |||
3. Calculation Methodology Qualification Procedure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37 | |||
4. Fluence Determination . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38 | |||
5 Measurement Qualification Procedure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39 iii | |||
==A. INTRODUCTION== | ==A. INTRODUCTION== | ||
This regulatory guide has been developed to provide state-of-the-art calculations | This regulatory guide has been developed to provide state-of-the-art calculations and measurement procedures that are acceptable to the NRC staff for determining pressure vessel fluence. | ||
The U.S. Nuclear Regulatory Commission (NRC) has promulgated regulations that ensure the structural integrity of the reactor pressure vessel for light-water-cooled power reactors. | |||
Chapter 15, | Chapter 15, Accident Analysis, of the Standard Review Plan (NUREG-0800, Standard Review Plan for the Review of Safety Analysis Reports for Nuclear Power Plants) assumes the pressure vessel does not fail. Specific fracture toughness requirements for normal operation and for anticipated operational occurrences for power reactors are set forth in Appendix G, "Fracture Toughness Requirements," of 10 CFR Part 50, "Domestic Licensing of Production and Utilization Facilities." The requirements of Appendix G are imposed by 10 CFR 50.60. Additionally, in response to concerns over potential pressurized thermal shock (PTS) events in pressurized water reactors (PWRs), the NRC issued 10 CFR 50.61, "Fracture Toughness Requirements for Protection Against Pressurized Thermal Shock Events." | ||
To satisfy the requirements of both Appendix G and 10 CFR 50.61, methods for determining the fast neutron fluence (E > 1 MeV) are necessary to estimate the fracture toughness of the pressure vessel materials. Appendix H, "Reactor Vessel Material Surveillance Program Requirements," to 10 CFR Part 50 requires the installation of surveillance capsules, including material test specimens and flux dosimeters, to provide data for material damage correlations as a function of fluence. | |||
The fracture toughness of pressure vessel materials is related to a parameter called the material's "reference temperature for nil-ductility transition," or simply reference temperature, and is denoted as RTNDT. The RTNDT is defined in Revision 2 of Regulatory Guide 1.99, "Radiation Embrittlement of Reactor Vessel Materials," by a correlation of the fluence (E > 1 MeV), material chemistry (concentrations of Cu and Ni), initial reference temperature, and margin to account for uncertainties in the correlation and input values. In 10 CFR 50.61, evaluation of the reference temperature based on the best estimate of the fast neutron fluence at the end of the license period is required, and the corresponding reference temperature is termed RTPTS. | |||
This guide is intended to ensure the accuracy and reliability of the fluence determination required by General Design Criteria (GDC) 14, 30, and 31 of Appendix A, General Design Criteria for Nuclear Power Plants, to 10 CFR Part 50. The guide describes methods and assumptions acceptable to the NRC staff for determining the pressure vessel neutron fluence. | |||
These methods are directly applicable to the determination of | These methods are directly applicable to the determination of RTNDT and RTPTS. The licensee may propose alternative methods. Alternative methods will be considered on a plant-specific basis, especially in cases involving unusual plant characteristics or factors that require different methods and assumptions. It is recognized that, when the embrittlement of a reactor vessel material is not significant and there is a large margin to the RTNDT limits, more approximate methods for determining the fluence may be appropriate (Reference 1). | ||
1 | |||
1 | Compliance with this guide is not a regulatory requirement of the USNRC. However, if a licensee elects to use this guide to determine pressure vessel neutron fluence, implementation of the guide would not be satisfied unless the licensee complies with certain specific provisions identified in the Regulatory Position of the guide. The use of the following terms is explained to clarify compliance with these regulatory positions. | ||
Must - Necessary provisions if implementation is to be satisfied. | |||
position are not available).May -Provisions that are acceptable and recommended, but are to be applied at the option of the licensee. | Should - Provisions that are expected to be complied with unless it is not possible because of specific circumstances (for example, data needed to meet the position are not available). | ||
May - Provisions that are acceptable and recommended, but are to be applied at the option of the licensee. | |||
is not required to respond to, the information collection. | The information collections contained in this regulatory guide are covered by the requirements of 10 CFR Part 50, which were approved by the Office of Management and Budget, approval number 3150-0011. If a means used to impose an information collection does not display a currently valid OMB control number, the NRC may not conduct or sponsor, and a person is not required to respond to, the information collection. | ||
==B. DISCUSSION== | ==B. DISCUSSION== | ||
The methods and assumptions described in this guide are for the calculation | The methods and assumptions described in this guide are for the calculation and measurement of vessel fluence for core and vessel geometrical and material configurations that are typical of current PWR and boiling water reactor (BWR) designs. This guide does not address the determination of surveillance specimen material properties or the correlation between material properties and neutron fluence. The methodology presented is intended as a best estimate, rather than a bounding or conservative fluence determination. For example, in the RTPTS correlation called for in 10 CFR 50.61, the best-estimate fluence is used to calculate the shift in RTPTS. | ||
fluence has been selected as the exposure parameter for the | |||
Uncertainty in the shift prediction (e.g., from uncertainty in the fluence, chemistry factor, or shift correlation) has been included separately in an explicit margin term. While the E > 1 MeV | |||
fluence has been selected as the exposure parameter for the RTNDT and RTPTS correlations, the procedures described in this guide determine the damage fluence spectrum (from 0.1 to 15 MeV) | |||
and are applicable to other exposure units, such as iron displacements per atom - dpa (Ref. 2). | and are applicable to other exposure units, such as iron displacements per atom - dpa (Ref. 2). | ||
The determination of the pressure vessel fluence is based on both calculations and measurements; the fluence prediction is made with a calculation, and the measurements are used to qualify the calculational methodology. | The determination of the pressure vessel fluence is based on both calculations and measurements; the fluence prediction is made with a calculation, and the measurements are used to qualify the calculational methodology. Because of the importance and the difficulty of these calculations, the methods must be qualified by comparison to measurements to ensure a reliable and accurate vessel fluence determination. In this qualification, measurement-to-calculation comparisons are used to identify biases (i.e., systematic errors) in the calculations and to provide reliable estimates of the fluence uncertainties.1 When the measurement data are of sufficient quality and quantity that they allow a reliable estimate of the calculational bias (i.e., they represent a statistically significant measurement data base), the comparisons to measurement may be used to | ||
1 For a discussion of the terms bias and uncertainty as used in this guide, see Sections 3.2 and 3.3 of Reference 3. | |||
2 | |||
(1) determine the effect of the various modeling approximations and any calculational bias and, if appropriate, (2) modify the calculations by applying a correction to account for bias or by model adjustment or both. Typically, plant-specific data alone are not sufficient to determine a calculational bias for an individual plant. As an additional qualification, the sensitivity of the calculation to the important input and modeling parameters must be determined and combined with the uncertainties of the input and modeling parameters to provide an independent estimate of the overall calculational uncertainty. | |||
The prediction of the vessel fluence must be made by an "absolute" fluence calculation in which the transport of the neutron flux is calculated from the core out to the vessel and cavity, rather than a simple spatial extrapolation of the fluence measurements. | |||
The calculations of the pressure vessel fluence consist of the following steps: (1) | |||
determination of the geometrical and material input data, (2) determination of the core neutron source, (3) propagation of the neutron fluence from the core to the vessel and into the cavity, and | |||
(4) qualification of the calculational procedure. These steps are discussed in detail in Regulatory Positions 1.1 through 1.4. In Regulatory Position 2, the use of surveillance dosimetry as an in situ verification for the calculations is described. Reporting is discussed in Regulatory Position 3. | |||
The major regulatory positions are summarized in Table 1, Summary of Regulatory Positions on Calculation and Dosimetry. | |||
As an indication of current practice, selected codes and cross-section libraries are listed in the references; however, it is the responsibility of the licensee to demonstrate their acceptability in a specific application. | |||
As an indication of current practice, selected codes and cross-section libraries are listed | |||
==C. REGULATORY POSITION== | ==C. REGULATORY POSITION== | ||
1.NEUTRON FLUENCE CALCULATIONAL | 1. NEUTRON FLUENCE CALCULATIONAL METHODS | ||
This guide describes the application and qualification of a methodology acceptable to the NRC staff for determining the best-estimate neutron fluence experienced by materials in the beltline region of light water reactor (LWR) pressure vessels, as well as for determining the overall uncertainty associated with those best-estimate values. Both the deterministic discrete ordinates method and the statistical Monte Carlo transport method are considered acceptable for vessel fluence determination and are described in Regulatory Position 1.3. | |||
Appendix B, "Quality Assurance Criteria for Nuclear Power Plants and Fuel Reprocessing Plants," to 10 CFR Part 50 requires that this methodology be properly qualified. Qualification includes determination of the uncertainty in the reactor vessel fluence as described in Regulatory Position 1.4. The uncertainty of the fluence must be 20% (1 ) or less when the fluence is used to determine RTPTS and RTNDT for complying with 10 CFR 50.61 and Revision 2 of Regulatory Guide 1.99, "Radiation Embrittlement of Reactor Vessel Materials," respectively. It should be recognized that this 20% uncertainty value has been included in the margin term for the RTPTS. | |||
Uncertainty of fluence for other applications should be determined using Regulatory Position 1.4 and included as an uncertainty allowance in the fluence estimate, as appropriate for the specific | |||
3 | |||
application. For example, when performing probabilistic risk assessment (PRA) evaluations such as those described in Regulatory Guide 1.154, "Format and Content of Plant-Specific Pressurized Thermal Shock Safety Analysis Reports for Pressurized Water Reactors," an explicit uncertainty term is required in addition to the best-estimate fluence values. In this case, it is acceptable to use the fluence uncertainty determined per the Regulatory Position 1.4 of this guide. | |||
An overview of the calculational methods described in Regulatory Positions 1.1 through | |||
1.3 is provided in Figures 1 and 2. The procedures for qualifying these methods and determining the pressure vessel fluence are summarized in Figures 3 and 4. | |||
1.1 Input Data | |||
1.1.1 Materials and Geometry Detailed material and geometrical input data should be used to define the physical characteristics that determine the attenuation of the neutron flux from the core to the locations of interest on the pressure vessel. These data include material compositions, regional temperatures, and geometry of the pressure vessel, core, and internals. The geometrical input data include the dimensions and locations of the fuel assemblies, reactor internals (baffle, shroud, former plates, core support barrel, thermal shield, and neutron pads), the pressure vessel (including identification and location of all welds and plates) and cladding, and surveillance capsules. For cavity dosimetry, input data should also include the width of the reactor cavity and the material compositions of the support structure and concrete (biological) shielding, including water content, rebar, and steel. The data should, to the extent possible, be based on documented and verified as- built dimensions and plant-specific materials. The isotopic compositions of important constituent nuclides within each region should be based on as-built materials data. In the absence of plant- specific information, nominal compositions and design dimensions may be used; however, in this case conservative estimates of the variations in the compositions and dimensions should be made and accounted for in the determination of the fluence uncertainty (Regulatory Position 1.4.1). The determination of the concentrations of the major isotopes responsible for the fluence attenuation (e.g., iron and water) should be emphasized. The water number densities should be based on plant full-power operating temperatures and pressures, as well as standard steam tables. The data should account for axial and radial variations in water density caused by temperature differences in the core and inside the core barrel, as well as the presence of in-channel and downcomer voids in the case of BWRs. The in-channel and bypass water and fuel channel may be combined to determine a homogenized material region. This approximation has been evaluated in Reference 4 and results in less than a ~3% additional uncertainty in the vessel fluence. | |||
1.1.2 Cross-Sections The calculational method to estimate vessel fluence should use the neutron cross-sections over the energy range from ~0.1 MeV to ~15 MeV and should apply the latest version of the Evaluated Nuclear Data File (currently ENDF/B-VI). These data have been thoroughly reviewed and tested relative to experimental benchmarks.2 Cross-section sets based on earlier or equivalent | |||
2 It should be noted that in many applications the ENDF/B-IV and the first three MODs of the ENDF/B-V iron cross-sections result in as much as ~20% underprediction of the vessel inner-wall fluence and ~35% underprediction of the cavity fluence (Refs. | |||
5-7). Updated ENDF/B-VI iron cross-section data (Ref. 8) have been demonstrated to provide a more accurate determination of the flux attenuation through iron (Refs. 5, 6) and are strongly recommended. These new iron data are included in ENDF/B-VI. | |||
4 | |||
nuclear-data sets that have been thoroughly benchmarked for a specific application may be used for that specific application.3 However, when the evaluated cross-section data change, the effect of these changes on the licensee-specific methodology must be evaluated and the fluence estimates updated when the effects are significant (see reporting requirements in Regulatory Position 3). | |||
1.1.2.1 Multigroup Libraries. Transport theory codes that are used to determine the neutron fluence, which employ a multigroup approximation, require cross-section libraries in which the basic data contained within the ENDF files has been preprocessed into a multigroup structure. The development of a multigroup library should consider the adequacy of the group structure, the energy dependence of the flux used to average the cross-sections over the individual groups, and the order of the Legendre expansion of the scattering cross-section. Sufficient details of the energy and angular dependence of the differential cross-sections (e.g., the minima in the iron total cross-section) should be included to preserve the accuracy in attenuation characteristics. | |||
1.1.2 | 1.1.2.2 Constructing a Multigroup Library. The construction of the multigroup library involves the selection of a problem-independent, fine-multigroup, master library containing data for all required isotopes. This master library should include a sufficiently large number of groups | ||
(/100) that differences between the shape of the assumed flux spectrum and the true flux have a negligible effect on the multigroup data. This library typically includes ~50-100 energy groups above ~0.1 MeV. In addition, a minimum of a P-3 Legendre expansion of the scattering cross- section must be used for typical LWR configurations (Ref. 10). Several libraries satisfying these provisions are available from RSICC, the Oak Ridge National Laboratory Radiation Safety Information Computational Center (Refs. 11-13). | |||
The number of groups may be reduced, with little loss in accuracy, by collapsing the data in the master library over spectra that more closely approximate the true spectra. The vessel fluence calculations should be performed with a job library that has been determined by collapsing the master library. This reduction may be accomplished with a one-dimensional calculation that includes the discrete regions of the core, vessel internals, by-pass and downcomer water, pressure vessel, reactor cavity, shield, and support structures. The resulting job library should consist of multigroup cross-sections based on the region-specific isotopic compositions. This job library should include ~20 energy groups above ~0.1 MeV. The adequacy of the job library must be demonstrated by comparing calculations for a representative configuration performed with both the master and job libraries. In these comparisons, the threshold-detector cross-sections and reaction rates obtained with the fine-group structure should be preserved in the multigroup calculations. In addition to the validation of the collapsing of the master library, additional qualification of the job library is provided by the benchmarking comparisons of Regulatory Position 1.4.2. | |||
There are several ~50-group libraries available from RSICC that were generated using light-water reactor (LWR) spectra for the group collapsing from master libraries (Refs. 14-16) and may be used for LWR application. These libraries contain microscopic cross-sections as well as some premixed macroscopic cross-sections for relevant mixtures. Because these prepackaged | |||
3 The ASTM E1018 standard (Ref. 9) provides additional guidance on the selection of dosimeter reaction rate cross-sections. | 3 The ASTM E1018 standard (Ref. 9) provides additional guidance on the selection of dosimeter reaction rate cross-sections. | ||
5 | |||
libraries were designed specifically for LWR pressure vessel fluence calculations, their applicability in any atypical application (for example, when an accurate thermal flux is required) | |||
should be verified prior to use. | |||
1.2 Core Neutron Source The determination of the neutron source for the pressure vessel fluence calculations should include the temporal, spatial, and energy dependence together with the absolute source normalization. | |||
The spatial dependence of the source should be based on depletion calculations that simulate core operation or measured data. The depletion calculations can be performed in either two or three dimensions. Three-dimensional calculations will provide the source in both the radial and axial directions. If two-dimensional calculations are used, the axial effects should be based on measured data. | |||
The core neutron source should be determined by the power distribution (which varies significantly with fuel burnup), the power level, and the fuel management scheme. The detailed state-point dependence should be accounted for (Refs. 17, 18); however, if this is not feasible, averaging over the operating power distribution is acceptable and may be obtained by either (1) | |||
averaging representative power distributions within the cycle or (2) assuming the cycle-average assembly power distribution is well approximated by the accumulated-exposure distribution at the end of the cycle. | |||
A best-estimate power distribution may be used for reactor vessel neutron fluence calculations. However, this best-estimate power distribution must be updated if changes in core loadings, surveillance measurements, or other information indicate a significant change in projected fluence values. This updating may be avoided by using a conservative generic average power distribution, provided no measured distribution yields higher power levels for the important peripheral assemblies. | |||
The peripheral assemblies, which contribute the most to the vessel fluence, have strong radial power gradients, and these gradients should not be neglected (Ref. 19). In the case of PWRs, the fuel pins closest to the core periphery tend to have reduced relative power while for BWRs the peripheral fuel pins may have increased relative power. The pin-wise source distribution should be used for best-estimate calculations and the peripheral-assembly pin-wise source data should be obtained from core depletion calculations. The pin-wise source distribution should represent the absolute source distribution in the assembly. | |||
The local source should be determined as the product of the fission rate and the neutron yield. The energy dependence of the source (i.e., the spectrum) and the normalization of the source to the number of neutrons per megawatt must account for the fact that changes in the isotopic fission fractions with fuel exposure (caused by Pu build-up) result in variations in the fission spectra, the number of neutrons produced per fission, and the energy released per fission. | |||
These effects tend to increase the fast neutron source per megawatt of power for high-burnup assemblies. The variations in these physics parameters with fuel exposure may be obtained from | |||
6 | |||
standard lattice physics depletion calculations (Refs. 20, 21). This effect is particularly important for plants that have adopted the PWR low-leakage refueling schemes (vs. out-in-in three-batch fuel management) in which once-, twice-, or thrice-burned fuel is located in the high-importance peripheral locations (Refs. 22-25). The harder spectrum in the BWR fuel regions having a high void fraction will have a similar effect on the isotopic fission fractions and on the neutron-source normalization and spectrum. | |||
A planar-octant representation is generally acceptable for the octant-symmetric fuel- loading patterns typically employed in LWRs. However, in the case of BWRs in which the jet- pump positioning is quadrant-symmetric, a quadrant-symmetric model may be required. For determining the peak fluence, fuel-loading patterns that are not octant symmetric (e.g., as in certain low-leakage patterns) may be represented in octant geometry using the octant having the highest fluence. For evaluating dosimetry, the octant in which the dosimetry is located may be used. | |||
When the actual planar core rectangular geometry cannot be modeled (e.g., in the case of (r,) discrete ordinates calculations), the core geometry may be described using an (r,) planar representation. To accurately represent the important peripheral assembly geometry, a -mesh of at least 40 angular intervals over an octant must be applied. The (r,) representation should reproduce the true physical assembly area to within ~0.5% and the pin-wise source gradients to within ~10%. The assignment of the (x,y) pin-wise powers to the individual (r,) mesh intervals should be made on a fractional area or equivalent basis (Ref. 26). Reference 26 is particularly useful if the radial mesh is a function of (). The overall source normalization should be performed with respect to the (r,) source so that differences between the core area in the (r,) | |||
representation and the true core area do not bias the fluence predictions. A discussion of the uncertainties in determining the core neutron source is provided in Reference 27. | |||
3 | Determination of the three-dimensional (3-D) fluence at the vessel using (r,z)- and (r)- | ||
geometry calculations may also be appropriate (see Regulatory Position 1.3.4). If these calculations are used to provide an axial correction factor, since the (r,z)- and the (r)- dependent fluences enter as a ratio (see Equation 4), the source specification may be less stringent if consistent sources are used in both the (r,z) and (r) calculations. | |||
1.3 Fluence Calculation | |||
1. | 1.3.1 Discrete Ordinates Transport Calculation The transport of neutrons from the core to locations of interest in the pressure vessel may be determined with a two-dimensional discrete ordinates transport program (Refs. 28-30) in (r,) | ||
and, when appropriate, in (r,z)- and (r)- geometries.4,5 When calculating a horizontal plane of the core/vessel geometry in which the rectangular (x,y) geometry of the core boundary and the | |||
4 Additional information concerning the application of transport methods to reactor vessel surveillance is provided in ASTM | |||
Standard E 482-89 (Ref. 31) | |||
5 If DOT (Ref. 29) is used, the "-weighted" option (MODE-5 in DOT 4.3) is considered to be more accurate than the "weighted" | |||
option (MODE-3 in DOT 3.5 or 4.3) for flux extrapolation and is recommended. Also, a value of = 0.3 is generally adequate (Refs. 4, 32, 33). Reference 34 includes a description of the directional -weighted scheme which may also be used. | |||
7 | |||
( | cylindrical (r) geometry of the vessel are mixed, a more accurate description is provided by the variable (r,) mesh option (Ref. 29) and may be applied. | ||
Position 1 | An azimuthal () mesh using at least 40 intervals over an octant in (r,) geometry in the horizontal plane should provide an accurate representation of the spatial distribution of the material compositions and source described in Regulatory Position 1.2. The radial mesh in the core region should be ~2 intervals per inch for peripheral assemblies and may be much more coarse for assemblies more than approximately two assembly pitches removed from the core- reflector interface. In excore regions, a spatial mesh that ensures the flux in any group changes by less than a factor of ~2 between adjacent intervals should be applied, and a radial mesh of at least | ||
~3 intervals per inch in water and ~1.5 intervals per inch in steel should be used. Because of the relatively weak axial variation of the fluence, a coarse axial mesh of ~0.5 interval per inch may be used except near material and source interfaces, where flux gradients can be large. An S8 fully symmetric angular quadrature must be used as a minimum for determining the fluence at the vessel. However, in reactor cavity fluence calculations a higher order quadrature may be needed, depending on the width of the cavity and the axial location at which the fluence is being calculated. | |||
should be | Where computer-storage limitations prevent the implementation of these mesh densities in a single-model representation, the calculation should be performed in two or more "bootstrap" | ||
steps rather than compromising the spatial mesh or quadrature (the number of groups used usually does not affect the storage limitations, only the execution time). In this approach, the problem volume is divided into overlapping regions. In a two-step bootstrap calculation, for example, a transport calculation is performed for the cylinder defined by 0 < r < RN with a fictitious vacuum- boundary condition applied at RN. From this initial calculation a boundary source is determined at the radius RO = RN - and is subsequently applied as the internal-boundary condition for a second transport calculation from RO to R (the true outer boundary of the problem). The adequacy of the overlap region must be tested (e.g., by decreasing the inner radius of the outer region) to ensure that the use of the fictitious boundary condition at RN has not unduly affected the boundary source at RO or the results at the vessel. | |||
A point-wise flux convergence criterion of .0.001 should be used, and a sufficient number of iterations should be allowed within each group to ensure convergence. To avoid negative fluxes and improve convergence, a weighted-difference model should be used.5 The adequacy of the spatial mesh and angular quadrature, as well as the convergence criterion, must be demonstrated by tightening the numerics until the resulting changes are negligible (Ref. 32). In discrete ordinates codes, the spatial mesh and the angular quadrature should be refined simultaneously. In some cases, these evaluations can be adequately performed with a one- dimensional model. | |||
The application of the discrete ordinates methods described in this section is illustrated in detail in Reference 4. | |||
1.3.2 Monte Carlo Transport Calculation The transport of neutrons from the core to locations of interest in the pressure vessel may be determined by using Monte Carlo Transport (Refs. 4 and 35-46). The Monte Carlo method has | |||
8 | |||
the advantage of allowing an exact representation of the problem geometry (Refs. 43-44); i.e., the rectangular fuel assemblies, the barrel, thermal shield and vessel cylindrical regions, and the axial geometry including upper and lower core reflectors, axial dosimeter locations, and vessel circumferential weld locations. In addition, the Monte Carlo method allows a continuous (as well as a multigroup) energy description of the nuclear cross-sections and flux solution. It is noteworthy, however, that in order to obtain comparable accuracy this method typically requires substantially longer computing times than the discrete ordinates method. | |||
The | The cross-section data used in the Monte Carlo calculations should satisfy the requirements of Regulatory Position 1.1.2. However, if a continuous energy cross-section library is used, the requirements concerning the Legendre expansion of the scattering cross-sections and the number of energy groups are not applicable. In addition, to ensure an accurate integration over the dosimeter cross-section, the energy bin structure used to determine dosimeter response scoring should satisfy the requirements of Regulatory Position 1.1.2 concerning the selection of the multigroup library group structure. | ||
source | While the Monte Carlo programs allow a great deal of flexibility in defining problem geometry, the input is generally relatively complex, requiring the definition of a large number of surfaces and cells, and is a potential source of error. Consequently, the geometry input should be checked thoroughly (using plots, volume checks, etc.) against the model design data to validate the definition of the geometry. | ||
The | In pressure vessel evaluations, the flux is required at specific dosimeter and vessel surface and internal locations. The calculated flux is edited in area and volume tally regions at these locations. Typically, the size of these tally regions is increased to increase the number of particle histories and reduce the statistical error in the flux estimate. The sensitivity of the flux estimate to the volume or area of the tally region should be determined to ensure that the spatial dependence of the flux does not introduce a bias in the flux estimate. If the size of the tally region introduces a bias into the flux estimate, the calculation should be rerun with a smaller tally region and an increased number of particle histories or the Monte Carlo prediction should be adjusted to eliminate the calculational bias. This adjustment may be made, for example, using a correction based on a deterministic calculation. | ||
In addition to the area and volume tally regions, point detectors may be used to estimate the flux at the locations of interest. However, the use of point detectors is complicated and the following cautions should be observed. To minimize the variance in point detector edits, Monte Carlo codes may define an average-flux region surrounding the detector within which the collision estimation is determined. The neutron flux in this region is assumed to be uniform. In order to ensure that the point estimate is not biased, the flux estimate should be shown to be insensitive to the selection of the average-flux region. If the size of the average-flux region introduces a bias into the flux estimate, the calculation should be rerun with a smaller average-flux region and an increased number of particle histories. No material boundaries should be included in the average- flux region surrounding the point detector. Since the next-event particle tracking used in point detectors does not account for reflecting, periodic, or white boundary conditions, the point detector flux estimate is underpredicted when particles from these boundaries can contribute to the flux estimate. Therefore, point detectors should not be located where particles from these | |||
9 | |||
boundaries can contribute to the flux estimate. In order to ensure that a sufficient number of particle histories are included and that the region close to the detector is adequately sampled, the statistical edits should be reviewed thoroughly, especially if the relative error is greater than 5%. | |||
In the Monte Carlo analysis, a calculational uncertainty is introduced as a result of the finite number N of particle histories sampled. The Monte Carlo analysis provides estimates of the mean flux, relative error (standard deviation/mean) and terms related to the higher moments, such as the variance of the variance (VOV), based on the calculated particle histories. In order to ensure that the calculation has converged and the estimated (mean) flux is valid, a reliable estimate of the sampling uncertainty is required. However, because of the finite sampling of the physical phase space, important paths in the geometry and peaks in the cross-sections may not be adequately sampled, resulting in significant errors in the estimated mean and standard deviation. | |||
Monte Carlo programs include statistical criteria to provide assurance that the calculations have converged and that the estimated mean and relative error are valid. When performing Monte Carlo pressure vessel fluence calculations, all statistical tests provided by the Monte Carlo code should be satisfied or justification for accepting the results should be reported. The Monte Carlo code should provide sufficient statistical testing (in addition to the usual tally standard deviation) | |||
to ensure that the scoring phase space has been adequately sampled. Representative statistical tests that are used to ensure tally convergence include (e.g., Reference 43): | |||
(1) The flux mean should not have a significant monotonic dependence (either increasing or decreasing) on N for the last half of the problem, | |||
(2) The percent relative error of the mean flux should be less than ~10%, except for point detectors for which it should be less than ~5%, | |||
(3) The relative error should have a ~ 1/ N dependence for the last half of the problem, | |||
(4) The figure of merit defined as | |||
1 FOM ' (Equation 1) | |||
R 2T | |||
where R is the relative error and T is the problem computing time, should not have a significant monotonic dependence on N for the last half of the problem, and | |||
(5) The VOV should be less than ~0.10 for all tallies and decrease as ~1/N for the last half of the problem. | |||
As a result of the strong flux attenuation between the core and the pressure vessel, very few neutrons (~less than one in a thousand) actually reach the vessel. Since the relative error (and variance) in the Monte Carlo calculation only decreases as 1/ N , the number of neutrons tracked and resulting computer times required to obtain acceptable accuracy are very large. In practice, the computing time may be reduced by using variance reduction techniques that increase the number of neutrons reaching the locations of interest (such as the dosimetry capsules and vessel inner-wall) or improve the particle tracking efficiency, or both. Typically, the number of particles that contribute to the tally is increased by increasing the number of particles in the regions of the particle phase space that are judged (based on prior knowledge) to contribute significantly to the | |||
10 | |||
calculation. However, if this judgment is incorrect and the neutrons believed to be unimportant make a significant contribution to the tally, the Monte Carlo estimated fluence and relative error (and variance) will be erroneous. Consequently, the variance reduction methods used in performing pressure vessel fluence analyses should be qualified by comparing the Monte Carlo variance reduction predictions with estimates made without the application of the variance reduction technique. Since the calculation without variance reduction will typically involve a large number of particle histories and long computing times, this comparison can be made using a simpler representative configuration (e.g., one-dimensional geometry). | |||
The following Monte Carlo variance reduction methods may be used in performing pressure vessel fluence analyses: (1) neutron energy cutoff, (2) source biasing, (3) geometry splitting with Russian Roulette, and (4) weight windows. Specific concerns associated with the application of these methods are described in the following. | |||
In the application of the neutron energy cutoff technique, neutron tracking efficiency is improved by terminating tracking when the particle weight falls below a selected input energy. | |||
This method assumes that particles with energy below the cutoff do not contribute to the tally and, consequently, may underestimate or bias the tally edit. While a cutoff energy of 1-MeV can be used to calculate the E > 1-Mev fluence without introducing any error, the cutoff energy for the dosimeter response calculation depends on the specific dosimeter reaction cross-section and should be selected low enough to ensure that the calculation provides acceptable accuracy. If this method is employed, the bias introduced by the energy cutoff should be estimated by comparison with an unbiased calculation and the Monte Carlo calculation adjusted to account for the bias. | |||
The source biasing method increases the number of source particles in the regions of phase space that are believed (based on prior knowledge) to provide the most significant contribution to the Monte Carlo tally. For example, the neutrons born in the peripheral fuel assemblies contribute most to the vessel fluence and the number of these particles can be increased (and their weight decreased proportionally) to improve the particle tracking efficiency. The space- and energy- dependent importance of the core neutrons to a specific location on the vessel can be estimated using a precalculated adjoint flux and used to determine the source biasing (e.g., Refs. 41, 42, and | |||
44). To ensure that the biased sampling of the source is producing reliable Monte Carlo estimates, the source biasing method should be qualified by comparison with an unbiased calculation for a representative problem. | |||
Geometrical splitting with Russian Roulette reduces the variance of the tally edit by increasing the number of neutrons that reach the vessel. This is accomplished by assigning numerical importance values to each spatial cell such that the cell importance values increase as the neutrons approach the vessel (Ref. 46). When a neutron travels toward the vessel passing from a region of lower importance I to a region of higher importance IN, it is split into IN/I | |||
neutrons. Conversely, when a neutron travels away from the vessel passing from a region of higher importance I to a region of lower importance IN, it undergoes Russian Roulette with a survival probability of IN/I. The reliability of this technique depends on the selection of the cell importances I. It is recommended that the importance in adjacent cells should not change by more than a factor of ~3, and the cell radial optical thickness should be less than ~2 mean free paths (Ref. 43). To ensure that the selection of the cell importance distribution provides reliable Monte | |||
11 | |||
The | Carlo estimates, the geometrical splitting with Russian Roulette method should be qualified by comparison with an unbiased calculation for a representative problem. | ||
The weight windows method may also be used to improve the efficiency of the neutron tracking in Monte Carlo pressure vessel analyses (e.g., Refs. 41, 42, and 45). In the weight windows method, the neutron weight in each space-energy cell is limited to a specific range or weight window by splitting and Russian Roulette. The window is defined so that the neutron weight is inversely proportional to the importance of the space-energy cell so that the scoring variance is minimized (Ref. 43). In practice, the specification of the weight windows is difficult and requires a detailed knowledge of the space-energy-dependent cell importance. Consequently, in order to ensure that the definition of the weight windows is providing reliable Monte Carlo estimates, it is important that the weight windows method be qualified by comparison with a calculation performed without variance reduction. | |||
The application of Monte Carlo methods described in this section is illustrated in detail in Reference 4. | |||
1.3.3 Fluence Determination The transport calculations (both discrete ordinate and Monte Carlo) may be performed in either the forward or adjoint modes. When several transport calculations are needed for a specific geometry, assembly importance factors may be precalculated by either performing forward calculations with a unit source (with the desired pin-wise source distribution) specified in the fuel assembly of interest or by performing adjoint calculations. The adjoint fluxes are used to determine the fluence at a specific (field) location, while the forward fluxes from the unit-source calculations determine the fluence at all locations in the problem. Once calculated, these factors contain the required information from the transport solution, and by weighting the assembly importance factors with the source distribution of interest, the vessel (or capsule) fluence may be determined without additional transport calculations, assuming the in-vessel geometry, material, and source distribution within the assembly remain the same. | |||
In performing calculations of surveillance capsule fluence (Regulatory Position 1.4), it should be noted that the capsule fluence is extremely sensitive to the representation of the capsule geometry and internal water region (if present), and the adequacy of the capsule representation and mesh must be demonstrated using sensitivity calculations (as described in Regulatory Position | |||
1.4.1). The capsule fluence and spectra are sensitive to the radial location of the capsule and its proximity to material interfaces (e.g., at the vessel, thermal shield, and concrete shield in the cavity), and these should be represented accurately. The core shroud and baffle former plates can result in a 5-10% reduction in the surveillance capsule dosimeter response and should be included in the model. | |||
When fluence reduction schemes have introduced strong axial or azimuthal heterogeneities into the source (e.g., an axially zoned replacement of fuel with stainless steel for fluence reduction), these should be modeled in detail. When the transport calculation is performed using the discrete ordinates method, a finer spatial mesh and tighter convergence criteria may be appropriate to ensure an accurate solution. These schemes may also entail a 3-D flux calculation. | |||
12 | |||
To account for the neutron spectrum dependence of RTNDT when the E > 1 -MeV fluence is extrapolated from the inside of the pressure vessel to the T/4- and 3T/4-vessel locations, a spectral lead factor (which accounts for the change in neutron spectrum between downcomer and vessel internal locations) must be applied to the fluence for the calculation of RTNDT (Ref. 47). This spectral lead factor has been included in the Equation 3 attenuation formula of Revision 2 of Regulatory Guide 1.99, and therefore is not required when this formula is used. However, when this formula is not used, the spectral lead factor must be applied to the fluence at the vessel internal locations. Displacements per atom is used as the extrapolation parameter in Regulatory Guide 1.99 and is an acceptable parameter for extrapolating the E > 1-MeV fluence to vessel internal locations and determining this lead factor. | |||
1.3.4 Synthesis of the 3-D Fluence When 3-D effects are important and three-dimensional transport calculations are not practical, a 3-D fluence representation may be constructed by synthesizing calculations of lower dimensions using the expression Ng(r,2,z) = Ng(r,2)(Lg(r,z) (Equation 2) | |||
where g(r,) is the group-g transport solution in (r,) geometry for a representative plane and Lg(r,z) is a group-dependent axial shape factor. Two simple methods available for determining Lg(r,z) are defined by the expressions Lg(r,z) ' P(z) (Equation 3) | |||
where P(z) is the peripheral-assembly axial power distribution, and L g(r,z) = Ng(r,z)/Ng(r) (Equation 4) | |||
where g(r) and g(r,z) are one- and two-dimensional group-g flux solutions, respectively, for a cylindrical representation of the geometry that preserves the important axial source and attenuation characteristics (Refs. 4 and 48). The (r,z) plane should correspond to the azimuthal location of interest (e.g., peak vessel fluence or dosimetry locations) or a conservatively - | |||
averaged (r,z) plane. The source per unit height for both the (r,)- and (r)- models should be identical, and the true axial source density should be used in the (r,z) model. | |||
Equation 3 is only applicable when (1) the axial source distribution for all important peripheral assemblies is approximately the same or is bounded by a conservative axial power shape and (2) the attenuation characteristics do not vary axially over the region of interest. In addition, since the axial flux distribution tends to flatten as it propagates from the core to the pressure vessel, for typical axial power shapes, use of Equation 3 will tend to overpredict axial flux maxima and underpredict minima. This underprediction is nonconservative and can be large near the top and bottom reflectors, as well as when minima are strongly localized as occurs in some fluence-reduction schemes. Because of the neglect of the former plates, this method may result in an additional 5-10% underprediction of the fluence. | |||
13 | |||
Equation 4 is applicable when the axial source distribution and attenuation characteristics vary radially but do not vary significantly in the azimuthal () direction within a given annulus. | |||
For example, this approximation is not appropriate when strong axial fuel-enrichment variations are present only in selected peripheral assemblies. | |||
In summary, an (r,)-geometry fluence calculation and a knowledge of the peripheral assembly axial power distribution are needed when using Equation 3. Use of this equation may result in fluence overpredictions near the midplane at relatively large distances from the core, e.g., | |||
in the cavity, and underpredictions at axial locations beyond the beltline at relatively large radial distances from the core. Conservatism may be included in the latter case by using the peak axial power for all elevations. | |||
Both radial and axial fluence calculations are needed when using Equation 4; thus, it is generally more accurate in preserving the integral properties of the three-dimensional fluence. | |||
Both Equation 3 and Equation 4 assume separability between the axial and azimuthal fluence calculations, which is only approximately true. | |||
When these simple synthesis techniques are not applicable, multichannel synthesis methods may be used. In the multichannel synthesis calculation, the fluence is represented as the sum N | |||
Ng(r,2,z) = j aiNgi(r,2)Ngi(r,z)/Ngi(r) (Equation 5) | |||
i'1 where the gi are basis flux solutions, typically representing specific regions of the core and vessel geometries, and the weighting coefficients ai are determined to provide an optimum prediction of the vessel fluence. It should be emphasized, however, that the accuracy of this method is sensitive to the selection of the basis functions, especially at region interfaces, and three-dimensional calculations should be considered where strong axial or azimuthal heterogeneities exist. This synthesis technique has been applied to a calculational benchmark in Reference 28 and to an experimental benchmark in Reference 49. | |||
1.3.5 Cavity Fluence Calculations Accurate cavity fluence calculations are relied on to analyze dosimeters located in the reactor cavity (Ref. 50). The calculation of the neutron transport in the cavity is made difficult by (a) the strong attenuation of the E > 1 MeV fluence through the (5.5- to 10.0-inch thick) vessel and the resulting increased sensitivity to the iron inelastic-scattering cross-section and (b) the possibility of neutron streaming (i.e., strong directionally dependent) effects in the low-density materials (air and vessel insulation) in the cavity. Because of the increased sensitivity to the iron cross-sections, ENDF/B-VI cross-section data should be used for cavity fluence calculations.2 As indicated in Regulatory Position 1.1.2, properly benchmarked alternative cross-sections may also be used. However, for cavity applications, the benchmarking must include comparisons for operating reactor cavities or simulated cavity environments. | |||
14 | |||
Typically, the width of the cavity together with the close-to-midplane locations of the dosimetry capsules result in minimal cavity streaming effects. Consequently, in discrete ordinates calculations, an S8 angular quadrature is acceptable. However, when off-midplane locations are analyzed, the adequacy of the S8 quadrature to determine the streaming component must be demonstrated with higher-order Sn calculations. In addition, since the radial mesh in the (r,z) | |||
in | calculation is generally finer than the z-mesh in the cavity resulting in narrow spatial-mesh intervals, a -weighted difference model should be used.5 The cavity fluence is sensitive to both the material compositions (Regulatory Position | ||
1.1.1) and the local geometry (e.g., the presence of detector wells) of the concrete shield, and these should be represented as accurately as possible. Benchmark measurements involving simulated reactor cavities that are recommended for methods evaluation are given in Regulatory Position 1.4.2.2. The measured energy spectrum for a typical PWR cavity is described in Reference 51. When both in-vessel and cavity dosimetry measurements are available, an additional verification of the measurements and calculations may be made by comparing the vessel inner-wall fluence determined from (1) the absolute fluence calculation, (2) the extrapolation of the in-vessel measurements, and (3) the extrapolation of the cavity measurements. | |||
Measurements performed in reactor cavities and the associated calculations are described in References 5 and 6. | |||
1.4 Methodology Qualification and Uncertainty Estimates The neutron transport calculational methodology must be qualified, and flux uncertainty estimates must be determined. The neutron flux undergoes several decades of attenuation before reaching the vessel, and the calculation is sensitive to the material and geometrical representation of the core and vessel internals, the neutron source, the nuclear cross-section data, and the numerical schemes used in its determination. The uncertainty estimates are used to determine the appropriate uncertainty allowance to be included in the application of the fluence estimate. While adherence to the guidelines described here will generally result in accurate fluence estimates, the overall methodology must be qualified in order to quantify uncertainties, identify any potential biases in the calculations, and provide confidence in the fluence calculations. In addition, while the methodology, including computer codes and data libraries used in the calculations, may have been found to be acceptable in previous applications, the qualification ensures that the licensee's implementation of the methodology is valid. | |||
2 | The methods qualification consists of three parts: (1) the analytic uncertainty analysis (Regulatory Position 1.4.1), (2) the comparison with benchmarks and operating reactor measurements (Regulatory Position 1.4.2), and (3) the estimate of uncertainty in the calculated fluence (Regulatory Position 1.4.3). These three steps in the overall qualification procedure are discussed below and are outlined in Figure 3. | ||
1.4.1 Analytic Uncertainty Analysis To demonstrate the accuracy of the methodology, an analytic uncertainty analysis must be performed. This analysis includes identification of the important sources of uncertainty. For typical fluence calculations, these sources include: | |||
15 | |||
C Nuclear data (transport and dosimeter reaction cross-sections and fission spectra), | |||
C Geometry (location of internals and deviations from the nominal dimensions), | |||
C Isotopic composition of material (density and composition of coolant water, core, core barrel, former plates, thermal shield, pressure vessel with cladding, and concrete shield), | |||
C Neutron sources (space and energy distribution, burnup dependence), | |||
C Methods error (mesh density, angular expansion, convergence criteria, macroscopic group cross-sections, fluence perturbation by surveillance capsules, spatial synthesis, and cavity streaming). | |||
This list of uncertainty components is not necessarily exhaustive, and other uncertainties that are specific to a particular reactor or a particular calculational method should be considered. | |||
In typical applications, the fluence calculational uncertainty is dominated by a few easily identified uncertainty components, such as the geometry and core neutron source, which results in a substantial simplification of the uncertainty analysis. | |||
The sensitivity of the flux to the significant component uncertainties should be determined by a series of transport sensitivity calculations in which the calculational model input data and modeling assumptions are varied and the effect on the calculated flux is determined. (A typical sensitivity would be -10-15% decrease in vessel E > 1 MeV fluence per centimeter increase in vessel inside radius.) Estimates of the expected uncertainties in these input parameters must be made and combined with the corresponding fluence sensitivities to determine the total analytically determined calculation uncertainty ca (i.e., standard deviation) as indicated in Figure 3. The uncertainties should be combined in a statistical (root-sum-of-the-squares) fashion to determine the total fluence uncertainty, accounting for the correlation between uncertainties when necessary (see, e.g., Section 5.2 of Ref. 3). The known systematic errors (or biases) should be combined algebraically, recognizing the sign of each contribution, to determine the overall calculational bias Bca (Ref. 52).6 The component uncertainties should be based on measurement or on the acceptable deviations included in the design specifications. The sensitivity calculations may be performed in one dimension when the model sensitivity does not require a two-dimensional representation. | |||
If the sensitivities are determined using Monte Carlo transport calculations, the magnitude of the model variations and the number of particle histories may have to be increased to provide reliable estimates of the calculation sensitivities. If the increased model variations are large, several calculations may be required to determine the non-linear dependence of the sensitivity. | |||
Perturbation calculations using correlated sampling may also be used to determine the sensitivity to model geometry, material compositions, and cross-sections. Adjoint Monte Carlo calculations may be used to determine the sensitivity to input such as fission spectra. Alternatively, if the sensitivities can be determined with sufficient accuracy in one-dimensional geometry, the sensitivities can be determined using a one-dimensional model. | |||
6 The systematic errors affect the fluence prediction in a specific direction (i.e., increase or decrease the estimate). Consequently, when these errors (or biases) are applied, the signs (or directions) must be recognized and they must be combined algebraically. | |||
The independent random uncertainties have no specific direction (i.e., an increase or decrease) associated with them, and they are combined without recognizing the sign of the error. | |||
16 | |||
A sensitivity analysis, in which the influence of each of these uncertainties on the calculated group fluences has been considered, is included in References 5 and 6. Since the uncertainties used in these analyses are common to many pressurized water reactors, the sensitivities (including correlations) may be used to determine initial uncertainty estimates. These variance estimates should be modified as the methods and cross-section data change or if the reactor of interest differs substantially from the reactors analyzed in these references. The referenced sensitivity analysis provides guidance for such modifications. | |||
1.4.2 Comparisons with Benchmark Measurements and Calculations Calculational methods must be validated by comparison with measurement and calculational benchmarks. The fluence calculation methods must be validated against (1) | |||
operating reactor measurements that provide in-vessel surveillance capsule dosimetry or ex-vessel cavity measurements or both, (2) a pressure vessel simulator benchmark that provides | operating reactor measurements that provide in-vessel surveillance capsule dosimetry or ex-vessel cavity measurements or both, (2) a pressure vessel simulator benchmark that provides measurements at the inner surface and at the T/4 and 3T/4 positions within the vessel (see Regulatory Position 1.4.2.2 for a discussion of such benchmarks), and (3) the fluence calculation benchmark. The results of the validation should include comparisons of reaction rates, fluences, and group fluxes for the locations of interest (Refs. 53, 54). | ||
The methods used to determine the plant-specific data and calculate the benchmarks must be consistent (to the extent possible) with those used to calculate the vessel fluence. That is, the same cross-sections, transport techniques, and transport code parameters that are to be used in the reactor licensing application must be employed in the calculation of the benchmark measurements and calculations. | |||
measurements | Differences between measurements and calculations should be consistent with the combined uncertainty estimates for the measurements and calculations. (Note that the uncertainties in both the calculations and measurements will contribute to the observed measurement-to-calculation differences.) The calculated reaction rates (using the methods described in Regulatory Positions 1.1 through 1.3) typically agree with the measurements to within about 20% for in-vessel surveillance capsules and 30% for cavity dosimetry. Deviations greater than these values must be investigated and, when the cause of the deviation is determined to be an error in the calculation, the calculations must be modified. | ||
If the calculation is modified using a least-squares statistical approach (Refs. 55-62), the measurement-to-calculation (M/C) ratios should be reviewed, prior to adjustment, in order to identify any systematic trends in the M/C data. Specific examples include systematic overprediction (or underprediction) of the dosimeter response for (1) dosimeters with low-energy thresholds, (2) dosimeters in the cavity (versus inner-wall and accelerated capsule locations), (3) | |||
dosimeters with long half-life reaction products, and (4) dosimeters monitoring specific cycles. If systematic trends, interdependence, or correlations appear in the M/C data, these should be explained and removed (prior to adjustment) or accounted for in the statistical adjustment | |||
17 | |||
procedure.7 Also, since the fluence adjustment is sensitive to the input uncertainties, realistic uncertainty values should be used and documented. The dosimeter measurement uncertainty should include the response uncertainty caused by dosimeter mislocation. Additional guidance on least-squares adjustment methods is provided in Ref. 63. It should be noted, however, that the available least-squares statistical adjustment codes differ substantially in their approach and their capabilities. Consequently, the selection of the least-squares code and the details of the analysis are highly application-specific and are outside the scope of this guide. | |||
and calculations. | The simulator benchmarks provide accurate measurements but typically do not provide an accurate representation of the actual plant configuration. The operating reactor measurements, on the other hand, represent the actual as-built plant configuration but, typically, include substantial measurement uncertainties. The simulator benchmarks together with the operating reactor measurements generally provide an acceptable measurement data base. The comparisons of the calculations to the measurement benchmarks should be used to estimate the calculational bias Bcb and uncertainty cb. If the bias estimate is considered reliable, a correction equal to the bias may be applied to the calculation to determine a best-estimate fluence.8 If the calculations or measurements are adjusted to improve the M/C agreement, using either a least-squares approach or a direct multiplicative bias, the adjustment and basis must be reported. | ||
1.4.2.1 Operating Reactor Measurements. Well documented fluence dosimetry measurements for operating power reactors may be used for methods and data qualification.9 Descriptions of these configurations include three-dimensional geometry, reactor operating conditions and, in some cases, both in-vessel and ex-vessel measurement data. | |||
Comparisons of measurements and calculations must be performed for the specific reactor being analyzed or for reactors of similar design. Plant-specific measurements have the advantage of including the as-built materials and geometry and the actual reactor operating conditions. An especially accurate determination of the fluence can be obtained when both in-vessel and cavity dosimetry are available. | |||
These measurements should not be used to bias or adjust the fluence calculations unless a statistically significant number of measurements is available, the various dosimeter measurements are self-consistent, and a reliable estimate of the calculational bias can be determined.8,9 Similarly, plant-specific biases should not be used unless sufficient reliable measurement data are available. As capsule and cavity measurements become available, they should be incorporated | |||
7 For example, if the M/C data for the Np-237 dosimeters tends to be high, indicating a systematically low calculation or high measurement, the errors in the Np-237 calculated or measured responses are correlated. If this error is not removed prior to the adjustment, the important sources of correlated error should be identified (e.g., the Np-237 dosimeter reaction cross-section, reaction product half-life and Np-237 photo-fission correction) and included in the adjustment procedure to allow the proper treatment of this error. | |||
8 The adequacy of the measurement data base for determining a bias in the calculations depends on the magnitude of the bias and is therefore problem dependent. For example, if the bias is small an accurate estimate of the bias will require either (1) a large number of reasonably accurate measurements or (2) a small number of very accurate measurements. Specifically, the uncertainty in the bias should be substantially less than the bias itself. For example, if the calculated fluence is to be increased by 10% , the uncertainty in this increase should be substantially less than 10%. | |||
9 Additional guidance on benchmark testing is provided in the ASTM standard E2005-99 (Ref. 64) and examples of operating reactor benchmark measurements are provided in Reference 65. | |||
18 | |||
into the operating reactor measurements data base, and the calculational biases and uncertainties should be updated as necessary. | |||
1.4.2.2 Pressure Vessel Simulator Measurements. Several pressure vessel simulator benchmarks are available (Refs. 7, 49, 65-82) and may be used for methods qualification. These benchmark experiments were carried out by several laboratories, and dosimetry measurements using different techniques were compared to provide experimental results with well known and documented uncertainties. Three configurations were used in these experiments. | |||
C Pressure Vessel Simulator Experiments. The Oak Ridge National Laboratory (ORNL) | |||
Pool Critical Assembly (PCA) is a full-scale-section mockup of a pressure vessel with dosimetry measurements at the inner surface of the vessel and at the T/4, T/2 and 3T/4 locations within the vessel wall (Refs. 72-75). The PCA benchmark experiment includes a simulated surveillance capsule that allows the determination of the effect of the capsule on the dosimetry measurements. This benchmark is characterized by a relatively simple geometry with generally less uncertainty in region compositions, temperatures, and source distributions than in operating power reactor benchmarks. | |||
C Experiments Simulating Operating Reactor Conditions. The CEN/SK Laboratory VENUS benchmark experiment (Refs. 49, 76-79) and the ORNL Pool Side Facility (PSF) | |||
benchmark experiment (Refs. 80-82) simulate the operating reactor: (1) pin-wise power distributions in peripheral fuel assemblies, (2) three-dimensional effects caused by partial- length shield assemblies, and (3) heterogeneities caused by neutron pads attached to the core barrel. | |||
C Pressure Vessel Cavity Measurements. The H.B. Robinson-2 cavity measurement (Refs. 65-67) and the NESDIP cavity simulator experiment (Refs. 7, 68-71) may be used as benchmarks to qualify cavity fluence calculational methods. Calculations of dosimetry measurements performed in an operating reactor cavity are described in References 5 and 6. | |||
1.4.2.3 Calculational Benchmark. The vessel fluence benchmark problems provided by the NRC in NUREG/CR-6115 (Ref. 4) should be used for methods qualification. This report provides the detailed information defining the problems and corresponding numerical solutions for a set of PWR and BWR pressure vessel fluence benchmark problems. The geometry, materials, and space- and energy-dependent source are fixed by the problem specification. PWR | |||
benchmark problems have been specified for (1) a standard core loading pattern, (2) a low-leakage core loading pattern, and (3) a partial-length shield assembly core design. Since BWR fuel loading patterns are presently not being designed for vessel fluence reduction, only a single BWR | |||
problem is specified. In addition, MCNP Monte Carlo calculations have been performed for both the BWR problem and the PWR problems for the standard core and partial-length shield assembly core loadings. Comparisons of the vessel fluence determined by using both the MCNP and DORT computer codes are included. | |||
The calculation of the benchmark problems allows a detailed assessment and verification of the numerical procedures, code implementation, and the various modeling approximations | |||
benchmark problems | 19 | ||
relative to state-of-the-art solutions for representative operating configurations.10 If the differences between the benchmark problem calculation and the reference solution are substantially larger than what would be expected based on the differences in the methods approximations and nuclear data used in the two calculations, the agreement is considered unacceptable. In this case, the calculation should be reviewed and the differences between the two solutions explained. (Note that NUREG/CR-6115 includes sensitivity calculations and comparisons with Monte Carlo predictions which provide an indication of the accuracy of the methods approximations used in the NUREG/CR-6115 reference calculations.) When the cause of the deviation is determined to be an error in the calculation, the calculational method must be revised. | |||
The calculation of the benchmark | 1.4.3 Estimate of Fluence Calculational Bias and Uncertainty The overall fluence calculation bias Bc and uncertainty c must be determined by an appropriate combination of (1) the analytic uncertainty analysis results of Regulatory Position | ||
1.4.1 (Bca and ca) and (2) the results of the uncertainty analysis based on the comparisons to the operating reactor and simulator benchmark measurements of Regulatory Position 1.4.2 (Bcb and cb). | |||
This combination may be a weighted average that accounts for the reliability of the individual estimates.11 The bias Bc may be applied as a multiplicative correction to the calculated fluence to determine the best-estimate value, and c should be used (when required) as the (1-sigma) | |||
uncertainty in the calculated best-estimate fluence. | |||
The fluence accuracy requirements are generally application specific; however, a vessel fluence uncertainty of 20% (1 sigma) is acceptable for RTPTS and RTNDT determination. In the determination of RTNDT and RTPTS, if the overall calculation uncertainty c is greater than 20%, the calculational model must be adjusted or a correction must be applied to reduce the difference between the fluence prediction and the upper 1-sigma limit to within 20%. For example, if c is greater than 20% but less than 30%, the fluence should be determined from the calculated fluence c by the relation N 6 1 B c [Fc(%) 20]/100 > Nc (Equation 6) | |||
If c is greater than 30%, the methodology of this regulatory guide is not applicable and the application will be reviewed on an individual basis. The procedure used to determine the fluence is illustrated in Figure 4. | |||
10 | |||
The calculation being used as the benchmark must be the NUREG/CR-6115 original referenced benchmark calculation, and not just a second independent calculation of the benchmark. | |||
11 Because the weighting of the analytic and benchmark uncertainty estimates depends on the details of the specific application that can vary widely, it is not possible to specify a practical and generically valid prescription for determining the weights. | |||
of the | However, the following example illustrates factors that should be considered. When as-built measurements of the vessel diameter are available and reasonable estimates of the core neutron source and other input uncertainties can be determined, the analytic uncertainty estimate should be reliable and have an uncertainty of -15% (1 sigma). Assume that there are a statistically significant number of accurate ( <~ 5%) operating reactor measurements and the uncertainty estimate based on this data has an uncertainty of -20%. The uncertainty estimate based on the vessel simulator measurements is assumed to be less certain and has an uncertainty of -25%. Using a weighted mean in which the weight is inversely proportional to the square of the standard deviation of the estimate (i.e., -2), the weights are wA = 0.5, wo = 0.3 and wS = 0.2. | ||
20 | |||
1 | The fluence calculation methods of this section are summarized in Table 1, Summary of Regulatory Positions on Calculation and Dosimetry. | ||
2. NEUTRON FLUENCE MEASUREMENT METHODS | |||
Dosimetry measurements provide independent estimates of specific activities and isotopic production rates that are used for validating the neutron transport calculations. Fluence is obtained from the response of passive integral detectors placed in surveillance capsules and, more recently, in the ex-vessel cavity. Procedures for performing the measurements and assessing uncertainty, including the validation of measurement methods and detector response in standard neutron fields, is described in this section. | |||
The fluence measurement methods of this section are summarized in Table 1, Summary of Regulatory Positions on Calculation and Dosimetry. The procedure used to qualify the measurement methods is presented in Figure 5. The guidance provided in this Regulatory Position is intended for new dosimetry and, to the extent practical, for existing dosimetry. | |||
2.1 Measurement Procedures Power reactor dosimetry measurement methods use passive integral detectors, which are typically activation detectors and solid-state track recorders. Most of the detectors respond to neutrons with energies above a characteristic reaction threshold. These detectors should be selected with substantial nonoverlapping energy regions (i.e., with well-separated thresholds) to provide coarse spectrum information as well as an estimate of the neutron fluence. | |||
2.1.1 Specification and Application of Dosimeters Neutron dosimetry for pressure vessel surveillance may consist of as-built packages of threshold dosimeters placed in surveillance capsules during reactor construction. The selected dosimeter set should provide adequate spectrum coverage. A common set of fast neutron integral detectors that may be employed in these packages is listed in Table 2. These neutron detectors are discussed in References 83 through 89. The 54Fe(n,p)54Mn and 58Ni(n,p)58Co activation reactions are useful for monitoring neutrons with energies above ~2 MeV. The 46Ti(n,p)46Sc reaction is useful for monitoring neutrons with energies above approximately 4 MeV, and the 63Cu(n,)60Co reaction is useful for measuring neutrons with energies above approximately 5 MeV. Fission monitors are important since they typically have reaction thresholds below 1 MeV. Activities determined by assaying one or more fission products from the 238U(n,f)FP fast-neutron fission reaction are useful for monitoring neutrons with energies from approximately 1.5 to 7 MeV. | |||
54Fe(n,p) | |||
63Cu(n, | |||
238U(n,f)FP fast-neutron | |||
Activities determined by assaying fission products from the | Activities determined by assaying fission products from the 237Np(n,f)FP reaction are useful for measuring neutrons with energies from approximately 0.7 to 6 MeV; however, there is a significant gamma-ray background associated with neptunium dosimetry that is due to the large quantity of naturally occurring 233Pa existing in secular equilibrium with 237Np. Compton scattering of the protactinium decay photons generates a large background continuum, which reduces the effective signal-to-noise ratio and lowers the accuracy of the radioassay. The | ||
237Np(n,f)FP reaction are useful | 93 Nb(n,n)93mNb activation reaction is useful for monitoring neutrons with energies above approximately 1 MeV, but involves the counting of a very soft photon, which can complicate the | ||
233Pa existing in secular equilibrium with | 21 | ||
corrections made to the measured activity to compensate for photon attenuation within the dosimeter. | |||
In addition to the fast-neutron threshold monitors, a thermal monitor such as 59Co should be included among the selected dosimeters to allow the determination of the thermal neutron fluence. The thermal fluence is needed to assess the effects of interfering activities in threshold monitors produced by low-energy activation of impurity elements. In instances where the thermal flux is relatively high, it can also be used for assessing the effects of burnout of the activated species. | |||
criteria and related recommendations in ASTM E 844 (Ref. 90) and E 1005 (Ref. 91) should be | Taken together with a low-energy detector such as cobalt (to estimate the thermal neutron fluence for determining interference from low-energy activations), the Table 2 dosimeter set provides satisfactory neutron energy spectrum coverage for pressure vessel dosimetry. Alternative detector sets that are used should provide equivalent spectrum coverage. Detector selection criteria and related recommendations in ASTM E 844 (Ref. 90) and E 1005 (Ref. 91) should be followed (see Table 2). | ||
The application of the dosimeter should be evaluated relative to the expected neutron flux and the total irradiation time in order to assess the effects of dosimeter burnout, and the burnout and burn-in of reaction products. For example, when using 58Ni monitors in the presence of a high thermal flux, the burnout of the nickel dosimeter as well as the 58Co and 58mCo reaction products should be considered. The resonant absorption of epithermal neutrons in 238U and the breeding of 239Pu, which can undergo thermal fission and contribute to the fission product activity, should also be considered. | |||
Application of activation detectors involves aspects of the measurement process that must be carefully controlled and documented to obtain accurate results and establish reasonable uncertainty estimates. Where applicable, procedures in ASTM Standards E 181 (Ref. 92), E 844 (Ref. 90) and E 1005 (Ref. 91) and methods devoted to individual radiometric sensors must be used as indicated in Table 2. Specific regulatory positions associated with the dosimetry measurements are indicated in the following. | |||
2.1.1.1 Dosimeter Nuclear and Material Proprieties. The dosimeter nuclear and material properties are important considerations in the selection and composition of the surveillance dosimetry package. The physiochemical properties of the dosimeter materials should be compatible with the prevailing service conditions. Dosimeter melting, for example, can result in the loss of monitor activity and unreliable measurements. Consequently, the melting point of the dosimeters and any associated materials should be carefully considered relative to the service environment ambient temperature and the radiation heating within the dosimeter. Although the melting point of typical radiometric monitors is sufficiently high to preclude melting, even mild heating can be problematical for certain dosimeters. For instance, diffusion losses of 54Mn activity from iron dosimeters can occur when heated to temperatures greater than 700-C. In addition, the dosimeter should be chemically stable and corrosion resistant. Fissionable materials are especially prone to oxidation and, consequently, fission monitors should be hermetically sealed to prevent oxidation and subsequent loss of material. | |||
22 | |||
The nuclear properties to be considered in dosimeter selection include the reaction cross- section, decay half-life, gamma-ray yield, and the fission yield when employing fission dosimeters. The reaction cross-section should be of sufficient magnitude to allow the dosimeter to be used as a fluence monitor, consistent with other requirements on the dosimeter mass and geometry, as discussed in following sections. | |||
The photon yield associated with the nuclear decay should be large enough to allow an accurate radio assay, and the energy of the particular photons counted should minimize interference from neighboring spectral lines. For example, the 24Na reaction product, created from the 27Al(n,) reaction, emits photons with energies of 1.369 and 2.754 MeV. Since there is less interference from the contaminant gamma rays associated with the higher-energy photon, it is preferable to count the high-energy line during the radioassay. In the case of fission product radioassay involving the 95Zr-95Nb fission product, three dominant spectral lines are observed at | |||
724, 756, and 765 keV. However, the two high-energy lines are in close proximity and complicate the determination of the background subtraction. Consequently, counting the 724 keV | |||
95Zr- | line provides the more accurate result. | ||
line provides the more accurate result. | |||
The half-life of the activated species should be long enough, relative to the elapsed time from the end of the irradiation to the time at which the dosimeters are actually counted, to permit nuclear counting to proceed with reasonable counting statistics. This time interval is often a significant consideration. The half-life should also be long enough to provide an accurate indication of the fluence during the irradiation period. | |||
monitors in power reactor applications. | The fission product yields associated with fission dosimeters should be of sufficient magnitude to permit an accurate radioassay. The yields (and decay constants) of the fission products 99Mo and 140Ba-140La are known with greater accuracy than most fission products; | ||
95Zr- | however, their half-life is generally too short for these isotopes to be useful as neutron fluence monitors in power reactor applications. Useful fission products include 95Zr-95Nb, 103Ru, 144Ce, and especially 137Cs-137mBa, because of its 30-year half-life. | ||
2.1.1.2 Material Composition. The dosimeter materials should be pure enough to ensure there is no significant error in the response of the dosimeter from extraneous activities. | |||
residuals in | Cursory specifications of materials regarding impurities are often unreliable. Specifically, fissile residuals in 237Np and 238U and minute amounts of cobalt in copper and nickel dosimeters (fractional parts per million) should be determined by mass spectrography or radioactivation analysis. Tracking the principal decay half-life is also a useful technique for determining the presence of extraneous activities from impurity elements. | ||
Dosimetric materials should be handled with care to prevent cross contamination of | Dosimetric materials should be handled with care to prevent cross contamination of the dosimeters and, when necessary, the dosimeters should be cleaned of surface contaminants. For example, the handling of aluminum monitors can introduce sodium (24Na) as an external surface contaminant. Neutron activation of the sodium will produce 24Na activity that will interfere with the radioassay of the 24Na activity produced from the desired dosimetric reaction. | ||
24Na activity that will interfere | |||
24Na activity produced from the desired dosimetric reaction | |||
Post-irradiation quality assurance may be accomplished at the National Institute for Standards and Technology (NIST) for dosimeters already in surveillance capsules by additional | |||
23 | |||
2.1.1. | irradiation in a thermal neutron field or a standard fission spectrum. These additional exposures and analyses can provide data on the dosimeter masses and impurities (for example, the original cobalt-59 impurity content in copper used for the copper-63 (n, ) reaction). | ||
2.1.1.3 Encapsulation. The detector capsule design must take into account possible activation interference and neutron spectrum perturbation. Thermal neutron shields that eliminate interference from thermal neutron reactions in some detectors must be designed to accommodate radiation heating and should be placed apart from low-energy detectors (see ASTM Standard E | |||
844, Ref. 90). The thickness of the thermal shield should be selected based on the expected thermal flux and duty cycle. (Note that the effective cadmium cutoff energy is a function of shield thickness.) The design of the shields and their placement relative to other dosimeters should be carefully considered in order to minimize their impact on neighboring foils. For example, the use of a cadmium shield will generate a local gamma field that could affect adjacent fission foils by producing a photofission component in the measured activity. | |||
the | 2.1.1.4. Isotopic Mass. The mass of the dosimeters should be selected to permit the production of sufficient activity to allow an accurate radioassay. However, a dosimeter mass that is too large can lead to high activities that result in excessive dead-time losses during nuclear counting and should be avoided. Metal alloys can be used for diluting materials that have high nuclear cross-sections (e.g., dilution of cobalt in an aluminum alloy). However, stoichiometry and isotopic analysis should be well documented for dosimeters that are not of pure natural elements. | ||
2.1.1.5. Geometry and Location. The geometric configuration of the dosimeters can have a significant effect on the dosimeter response. Most radiometric monitors are in the shape of thin circular activation foils, although other shapes are available. Dosimeters should be small enough to allow placement at the desired location and large enough to provide sufficient activity to allow an accurate measurement. The dosimeter thickness can have a significant effect on the dosimeter self-shielding (of the neutron flux) and the dosimeter response for strongly absorbing foils. In addition, increased dosimeter thickness increases the effects of scattering and photon attenuation within the dosimeter and can lead to large response corrections. When selecting fissionable dosimeters, it should be recognized that reduced foil thickness increases fission product loss from nuclear recoil out of the foil, which reduces the measured dosimeter activity. | |||
The placement of the dosimeters can also have a significant effect on the dosimeter measurement. The surroundings of a dosimeter (e.g., adjacent dosimeters or material interface) | |||
can influence detector response. Strongly and weakly absorbing foils that respond over the same energy range should be sufficiently separated within the dosimetry package to minimize interference effects. The arrangement of the dosimeters within the package should minimize scattering effects (viz., thick dosimeters should not be stacked so that large scattering corrections are required for adjacent foils). | |||
The location of individual dosimeters must be determined accurately and recorded, because fluence gradients in out-of-core positions are generally severe. In the pressure vessel cavity, establishing azimuthal position can be as important as the radial location. Specially designed mounting arrangements that ensure accurate radial, azimuthal, and vertical positions | |||
24 | |||
should be used for cavity dosimetry. Vertical wires may be used. Comprehensive and accurate detector location information should be maintained. | |||
2.1.1.6 Solid-State Track Recorders. Solid state track recorders (SSTRs) are integral detectors that employ fission reactions. These sensors directly record the tracks of fission fragments from a thin fissionable deposit (Refs. 93 and 94). The principal advantages of these detectors are wide sensitivity ranges and a permanent measurement record. Because the application of SSTRs employs fissionable deposits in the nanogram to picogram range, details of the measurements should be well documented, and standard neutron field calibration should be performed prior to the application. ASTM Standard E 854 (Ref. 95) provides additional information concerning the use of SSTRs. | |||
2.1.1.7 Helium Accumulation Fluence Monitors. The helium accumulation fluence monitor (HAFM) is another type of stable-product neutron monitor. These monitors rely on the generation of helium gas via (n,) reactions that occur in a variety of materials when exposed to a neutron field. The amount of helium generated in a specific material is a function of the neutron fluence and is determined by vaporizing the sample and measuring the helium content by mass spectrometry. HAFMs are routinely used in conjunction with radiometric monitors and, in fact, several commonly used radiometric monitors also undergo (n,) reactions and can be analyzed for helium accumulation as well as for radiometric activity. ASTM Standard E910 (Ref. 96) provides additional guidance on the use of HAFMs. | |||
2.1.2 Detector Response Measurements A measured response must be provided to allow comparison of the calculations and measurements. Typical measured responses include the specific activity at end-of-irradiation (given in disintegrations per second per nucleus), the measured isotopic production (for example, helium atoms per initial atom of material), and the total reactions (for example, fissions per initial atom of material). When comparing calculations with measurements, corrections should be included for detector response perturbations, interfering reactions, and, when applicable, burnup and photofission. In the calculation of end of irradiation activities (or in the conversion of end of irradiation activities into reaction rates), the power-history and the half-life of the dosimeter activation products must be included (Ref. 91). Photofission corrections can vary considerably (from 2-15%) depending upon the location of the dosimetry and the type of reactor. Fission yields should be those specified in the relevant ASTM Standards, the ENDF library, or the validated job library. In situ neutron field perturbations (e.g., by the surveillance capsule and detector encapsulation) must be accounted for if they are not an integral part of the neutron transport calculation. | |||
a | When reporting measurements, these corrections should be described and quantified along with any other effects that have a significant impact on the measurements. This is especially important because pressure vessel surveillance dosimetry often involves comparison of measurements carried out by different organizations and over long periods of time. | ||
2.1.3 Measurement Uncertainties Regulatory Position 1.4 states that the calculations must be validated by comparison with measured benchmarks. In order to perform this validation, the uncertainty associated with the | |||
25 | |||
measured response must be determined for each dosimeter type. The uncertainty must be included in the documentation of the measured results. | |||
The important uncertainties must be quantified and included in an uncertainty table that summarizes the specific components of uncertainty contributing to each detector response. | |||
Typical sources of uncertainty include (1) background counting, counting efficiency, counting statistics, and fraction of the sample counted, (2) decay constant and fission yield uncertainty, (3) | |||
uncertainty in detector weight, geometry, composition and isotopic purity, (4) power history and time of irradiation, (5) dosimeter location uncertainty, and (6) cross-section uncertainty when the fluence is determined. Additional sources are identified in Regulatory Positions 2.1.1 and 2.1.2 and in the dosimeter-specific ASTM Standards of Table 2. The entries in the table should be identified as standard deviations, upper bounds, or appropriate fractions of the correction. Each entry should be described, along with the method used to combine the entries and determine the total response uncertainty. The evaluation of the measurement uncertainty for each dosimeter is important for ensuring a meaningful comparison with the calculations. | |||
2.2 Validation in Standard and Reference Neutron Fields The dosimeter measurements used to benchmark the calculations must meet the requirements in Appendix B to 10 CFR Part 50 for a quality assurance program (for surveillance measurements). To ensure long-term measurement consistency and confirm measurement uncertainties, dosimetry measurements must be performed every few years12 in well characterized neutron fields. If there are changes in dosimeter evaluation procedures, equipment, or personnel that could have a significant effect on the dosimetry, the validation must be updated. This validation may be performed using the Materials Dosimetry Reference Facility (MDRF) or using reference fission neutron sources (Refs. 97, 98). Neutron field referencing may be used as a detector response calibration (Ref. 99). | |||
The validation is accomplished by exposing each type of detector to a certified neutron fluence in the reference neutron field and by determining the fluence using the measurement method to be validated. A calculated spectrum-averaged cross-section, generally specified along with the certified neutron fluence, must be used to derive the measurement value. If measured and certified neutron fluence agree (to within the combined uncertainty of the measured and certified fluence), the detector measurement method, including the detector cross-section, is validated. If the fluences disagree, this calculation-to-experiment (C/E) ratio represents a bias associated with the detector response measurement or the detector cross-section or both. In this case, the detector measurement methods and input parameters should be re-examined in order to eliminate the bias. If after re-examination the bias is still present, the bias may be used directly as a detector response calibration factor. | |||
The results of the neutron field validation procedure should be reported in terms of C/E | |||
ratios for the individual detectors and should include an uncertainty table. The standard neutron field validation may be used, as appropriate, to simplify the uncertainty table called for in | |||
12 This measurement validation procedure does not require transport calculations. | |||
26 | |||
Regulatory Position 2.1.3 by reducing or eliminating many uncertainties in the activity measurement, nuclear decay parameters, and detector cross-sections. | |||
Aside from validating the measurement method, standard and reference neutron fields may be used for quality assurance of critical features of the detector response. Examples are activation interference by impurities, proper determination of fission product activities, and mass assay of fissionable deposits for track recorders. | |||
Interlaboratory comparisons, such as those carried out under the NRC-sponsored LWR | |||
Surveillance Dosimetry Improvement Program (Refs. 7, 49, 72, 80, 81), have been useful in validating the quality assurance programs of various laboratories. These interlaboratory comparisons have been used to identify problems in the measurement procedures and establish consistency between the various participating laboratories. | |||
calculation, the | 2.3 Fluence Determination from Detector Measurements As stated in Regulatory Position 1.4, the calculated fast-neutron fluence values must be validated by comparison with measurement benchmarks. These comparisons validate the calculational methodology by determining the bias and uncertainty in the calculated fluences, E > | ||
1 MeV. A fast-neutron fluence, E > 1 MeV, must be obtained for each detector as the quotient of the measured reaction probability and the effective E > 1-MeV spectrum averaged cross-section (Ref. 100) based on a neutron transport calculation. These fluences and a suitably weighted average fluence must be reported together with the associated measurement uncertainty determined using the methods of Regulatory Positions 2.1.3 and 2.2 and accounting for the cross- section uncertainty. An alternative to deriving a neutron fluence from the detector responses is to directly compare specific end-of-irradiation activities, measured reaction probabilities, isotopic production, total reactions, average reaction, or measured reaction rates with results from the neutron transport calculation. The response selected should provide a reliable indication of the uncertainty and bias in the E > 1 MeV fluence calculation. | |||
should be | The measurement biases and uncertainties must be documented for the selected detector response. Based on the measurement-to-calculation (M/C) comparisons, an average M/C ratio should be determined as a suitably weighted average of the individual M/C values. In determining the weighting for this average, as a minimum the measurement and spectral average cross-section uncertainties should be considered. An estimate of the uncertainty in the neutron spectrum may be included as an additional contributor to the dosimeter weighting. | ||
When justification can be provided, an individual detector may be declared suspect and given reduced weight or discarded. However, in this case the justification and procedure used must be documented. | |||
2.4 Ex-vessel Dosimetry As-built in-vessel surveillance capsule dosimetry cannot easily be updated. Furthermore, the dosimetry irradiated with metallurgical specimens is only available at infrequent intervals. | |||
However, additional and upgraded dosimetry is important for understanding and following vessel | |||
27 | |||
exposures, especially for low-leakage core modifications. The ex-vessel cavity may be used as an alternative site for installing additional improved dosimetry. Recent pressure vessel benchmark experiments (Refs. 7, 67, 94) have demonstrated that the ex-vessel dosimetry can provide useful exposure information within the pressure vessel wall (Refs. 72, 81). When placed at appropriate circumferential locations, this dosimetry is a good monitor of the effectiveness of low-leakage core strategies. | |||
The | |||
===3. REPORTING=== | |||
When fluence determinations are required by the regulations, the licensee's documentation describing the determination of pressure vessel fluence must provide a complete description of the methods used to calculate and measure the neutron fluences. In applying the methodology of this guide, the details of the application and the results should be reported as described in this section. | |||
A topical report providing the details of the fluence methodology, as outlined in Regulatory Positions 1 and 2, should be satisfactory for compliance with this guide. Plant-specific measurements and data should be reported in the pressure vessel surveillance reports. | |||
The specific regulatory positions on reporting are given in this section and are summarized in Table 1, Summary of Regulatory Positions on Calculation and Dosimetry. | |||
1. | 3.1 Fluence Methods The methods used to calculate the integral and multigroup fluences and fluence rates and associated methods qualification should be reported. The calculational uncertainty analysis, benchmark comparisons, and the determination of the overall fluence uncertainty should be described in detail. If the calculated fluence is adjusted, the fluence adjustment and associated uncertainty together with justification must be reported. A discussion of any deviations from the procedures provided in this regulatory guide should be included. The source of the cross-section data, the numerical methods (e.g., quadrature, mesh, and convergence criteria), and the treatment of special effects (e.g., fuel burnup, axial effects, and pin-power distributions) should be described in detail. | ||
3.2 Multigroup Fluences The calculated absolute multigroup neutron fluences and fluence rates at the peak wall (or other limiting) location, surveillance locations, and T/4 and 3T/4 positions within the pressure vessel should be reported. The multigroup energy boundaries should be included. | |||
1 | 3.3 Integral Fluences The calculated E > 1 MeV integral fluences and fluence rates at the vessel inner wall locations (determined as described in Regulatory Position 1.4.3), together with the uncertainties, should be reported. | ||
28 | |||
1. | 3.4 Comparisons of Calculation and Measurement If the qualification of the calculation methods is performed using fluence comparisons, the measured and calculated E > 1 MeV integral neutron fluences should be reported. If the methods qualification is performed using reaction rate comparisons (or other responses as described in Regulatory Position 2.3), the calculated and measured reaction rates (or reaction probabilities) | ||
should be reported. The spectrum-averaged reaction cross-section used to relate the fluence and reaction rate and the method for its determination should also be provided. In either case, the M/C ratios and the measurement uncertainty should be reported for the average fluence and for each detector at each measurement location. If these measurements are used to adjust the calculation, the adjustment should be described and the M/C ratios before and after adjustment should be reported. | |||
The | 3.5 Specific Activities and Average Reaction Rates The specific activities at the end of irradiation and the measured average reaction rates should be reported, together with the associated uncertainty tables and the power-time history. | ||
For each dosimeter, provide the reaction type, dosimeter material and form (wire, foil, etc.), | |||
weight-percent and isotopic-percent of the target material, fission yields, and half-lives. The corrections for the detector response perturbations, interfering reactions, and photofission should also be described. The results of the measurement validation should also be described. If thermal-neutron fluence rate measurements have been performed, these should be reported together with the uncertainty and the dosimeter thermal-neutron shield material. | |||
==D. IMPLEMENTATION== | |||
The purpose of this section is to provide information to applicants and licensees regarding the NRC staffs plans for using this regulatory guide. | |||
Except in those cases in which an applicant or licensee proposes an acceptable alternative method for complying with the specified portions of the NRCs regulations, the methods described in this guide will be used in the evaluation of applications for new licenses and for evaluating compliance with 10 CFR 50.61 and Appendix G to 10 CFR Part 50, which is imposed by 10 CFR 50.60. | |||
29 | |||
TABLE | |||
===1. SUMMARY O=== | |||
==F. REGULATORY POSITION== | |||
S ON | |||
CALCULATION AND DOSIMETRY | |||
FLUENCE CALCULATION METHODS Regulatory Position Fluence Determination. Absolute fluence calculations, rather than extrapolated 1.3 fluence measurements, must be used for the fluence determination. | |||
Modeling Data. The calculation modeling (geometry, materials, etc.) should be 1.1.1 based on documented and verified plant-specific data. | |||
1. | Nuclear Data. The latest version of the Evaluated Nuclear Data File (ENDF/B) 1.1.2 should be used for determining nuclear cross-sections. Cross-section sets based on earlier or equivalent nuclear-data sets that have been thoroughly benchmarked are also acceptable. When the recommended cross-section data change, the effect of these changes on the licensee-specific methodology must be evaluated and the fluence estimates updated when the effects are significant. | ||
1. | Cross-Section Angular Representation. In discrete ordinates transport 1.1.2 calculations, a P3 angular decomposition of the scattering cross-sections (at a minimum) must be employed. | ||
The | Cross-Section Group Collapsing. The adequacy of the collapsed job library must 1.1.2 be demonstrated by comparing calculations for a representative configuration performed with both the master library and the job library. | ||
1. | Neutron Source. The core neutron source should account for local fuel isotopics 1.2 and, where appropriate, the effects of moderator density. The neutron source normalization and energy dependence must account for the fuel exposure dependence of the fission spectra, the number of neutrons produced per fission, and the energy released per fission. | ||
End-of-Life Predictions. Predictions of the vessel end-of-life fluence should be 1.2 made with a best-estimate or conservative generic power distribution. If a best estimate is used, the power distribution must be updated if changes in core loadings, surveillance measurements, or other information indicate a significant change in projected fluence values. | |||
1.3. | Spatial Representation. Discrete ordinates neutron transport calculations should 1.3.1 incorporate a detailed radial- and azimuthal-spatial mesh of ~2 intervals per inch radially. The discrete ordinates calculations must employ (at a minimum) an S8 quadrature and (at least) 40 intervals per octant. | ||
30 | |||
1.3. | Regulatory Position Multiple Transport Calculations. If the calculation is performed using two or 1.3.1 more "bootstrap" calculations, the adequacy of the overlap regions must be demonstrated. | ||
Point Estimates. If the dimensions of the tally region or the definition of the 1.3.2 average-flux region introduce a bias in the tally edit, the Monte Carlo prediction should be adjusted to eliminate the calculational bias. The average-flux region surrounding the point location should not include material boundaries or be located near reflecting, periodic, or white boundaries. | |||
Statistical Tests. The Monte Carlo estimated mean and relative error should be 1.3.2 tested and satisfy all statistical criteria. | |||
Variance Reduction. All variance reduction methods should be qualified by 1.3.2 comparison with calculations performed without variance reduction. | |||
Capsule Modeling. The capsule fluence is extremely sensitive to the geometrical 1.3.3 representation of the capsule geometry and internal water region, and the adequacy of the capsule representation and mesh must be demonstrated. | |||
Spectral Effects on RTNDT. In order to account for the neutron spectrum 1.3.3 dependence of RTNDT, when it is extrapolated from the inside surface of the pressure vessel to the T/4 and 3T/4 vessel locations using the E > 1-MeV fluence, a spectral lead factor must be applied to the fluence for the calculation of RTNDT. | |||
Cavity Calculations. In discrete ordinates transport calculations, the adequacy of 1.3.5 the S8 angular quadrature used in cavity transport calculations must be demonstrated. | |||
Methods Qualification. The calculational methodology must be qualified by both 1.4.1, | |||
(1) comparisons to measurement and calculational benchmarks and (2) an analytic 1.4.2, uncertainty analysis. The methods used to calculate the benchmarks must be 1.4.3 consistent (to the extent possible) with the methods used to calculate the vessel fluence. The overall calculational bias and uncertainty must be determined by an appropriate combination of the analytic uncertainty analysis and the uncertainty analysis based on the comparisons to the benchmarks. | |||
Fluence Calculational Uncertainty. The vessel fluence (1 sigma) calculational 1, uncertainty must be demonstrated to be #20% for RTPTS and RTNDT determination. 1.4.3 In these applications, if the benchmark comparisons indicate differences greater than 20%, the calculational model must be adjusted or a correction must be applied to reduce the difference between the fluence prediction and the upper 1-sigma limit to within 20%. For other applications, the accuracy should be determined using the approach described in Regulatory Position 1.4, and an uncertainty allowance should be included in the fluence estimate as appropriate in the specific application. | |||
31 | |||
Regulatory Position FLUENCE MEASUREMENT METHODS | |||
Spectrum Coverage. The set of dosimeters should provide adequate spectrum 2.1.1 coverage. | |||
2.1. | Dosimeter Nuclear and Material Properties. Use of dosimeter materials should 2.1.1 address melting, oxidation, material purity, total and isotopic mass assay, perturbations by encapsulations and thermal shields, and accurate dosimeter positioning. Dosimeter half-life and photon yield and interference should also be evaluated. | ||
2.1. | Corrections. Dosimeter-response measurements should account for fluence rate 2.1.2 variations, isotopic burnup effects, detector perturbations, self shielding, reaction interferences, and photofission. | ||
2. | Response Uncertainty. An uncertainty analysis must be performed for the 2.1.3 response of each dosimeter. | ||
Validation. Detector-response calibrations must be carried out periodically in a 2.2 standard neutron field. | |||
2. | Fast-Neutron Fluence. The E > 1 MeV fast-neutron fluence for each 2.3 measurement location must be determined using calculated spectrum-averaged cross-sections and individual detector measurements. As an alternative, the detector responses may be used to determine reaction probabilities or average reaction rates. | ||
The M/C ratios, the standard deviation and bias between calculation and measurement, must be determined. | Measurement-to-Calculation Ratios. The M/C ratios, the standard deviation and 2.3 bias between calculation and measurement, must be determined. | ||
REPORTING PROVISIONS | |||
Neutron Fluence and Uncertainties Details of the absolute fluence calculations, associated methods qualification and 3.1 fluence adjustments (if any) should be reported. Justification and a description of any deviations from the provisions of this guide should be provided. | |||
Calculated multigroup neutron fluences and fluence rates should be reported. 3.2 The value and basis of any bias or model adjustment made to improve the 3.2 measurement-to-calculation agreement must be reported. | |||
32 | |||
Regulatory Position Calculated integral fluences and fluence rates for E > 1 MeV and their uncertainties 3.3 should be reported. | |||
3.4 | Measured and calculated integral E > 1 MeV fluences or reaction rates and 3.4 uncertainties for each measurement location should be reported. The M/C ratios and the spectrum averaged cross-section should also be reported. | ||
The results of the standard field validation of the measurement method should be 3.5 reported. | |||
3.5 | Specific Activities and Average Reaction Rates The specific activities at the end of irradiation and measured average reaction rates 3.5 with uncertainties should be reported. | ||
All corrections and adjustments to the measured quantities and their justification 3.5 should be reported. | |||
33 | |||
TABLE 2. THRESHOLD DETECTORS RECOMMENDED FOR | |||
PRESSURE VESSEL DOSIMETRY | |||
Nominal Applicable Threshold (MeV) ASTM Standards | |||
237 Np(n,f)FP* 0.69 E 705 (Ref. 83) | |||
93 Nb(n,n)93mNb 0.97 E 1297 (Ref. 84) | |||
238 U(n,f)FP 1.45 E 704 (Ref. 85) | |||
58 Ni(n,p)58Co 2.05 E 264 (Ref. 86) | |||
54 Fe(n,p)54Mn 2.32 E 263 (Ref. 87) | |||
46 Ti(n,p)46Sc 3.76 E 526 (Ref. 88) | |||
63 Cu(n,)60Co 4.65 E 523 (Ref. 89) | |||
* FP indicates fission product. | |||
34 | |||
35 Cycle-Dependent Assembly/Pin Multi-Group Powers & Library Exposures (Optional) | |||
(1.1.1) (1.1.2.2) | |||
(1.2) | |||
Isotopic Component Macroscopic Continuous Fission MCNP* Material Cross- Energy Spectra Source Compositions Sections Library | |||
(1.1.1) | |||
(1.3.2) | |||
Core/Internals/ | |||
Vessel MCNP | |||
Geometry Fluence | |||
3-D Vessel Relevant Regulatory Positions are indicated. | |||
Fluence Estimate | |||
* The MCNP code is used as an example. | |||
Figure 2. Monte Carlo Calculation Methodology | |||
36 | |||
(1.4)* | |||
Fluence Calculation Calculation Methods Methods Revision (When Necessary) | |||
Benchmark Analytic Uncertainty Uncertainty Analysis (1.4.2) | |||
Analysis | |||
(1.4.1) | |||
(1.4.2.1) (1.4.2.2) (1.4.2.3) | |||
Sensitivity Calculation of Calculation Plant Plant Calculational Calculation Analysis Operating Simulator of Simulator Model Input Model Benchmark of Benchmark (M/MX) Reactor Model Input Benchmarks Input Model Input Problem Measurements Estimated Capsule/ Comparison Calculation Uncertainty Cavity Simulator Comparison Comparison of Results Benchmark Uncertainties Determination Measurement Measurement of Results of Results Solution | |||
& Biases (M/MX)*X Data Data (X) | |||
Combine Agreement Resolution of Ca , BCa Benchmark Acceptable Differences Uncertainty/ | |||
Bias Estimates | |||
(1.4.3) Cb, BCb Combine | |||
* Relevant Regulatory Positions are indicated. | |||
Analytic & | |||
Benchmark Estimates Calculational Uncertainty Based on Ca (Cb) Analytic (Benchmark) Uncertainty Analysis BCa (BCb) Calculation Bias Based on Analytic (Benchmark) Uncertainty Analysis C, B C | |||
C (BC) Estimated Calculation Uncertainty (Bias) | |||
Figure 3. Calculation Methodology Qualification Procedure | |||
*Relevant Regulatory Positions are indicated. | |||
6 | = [Equation (6)] | ||
Figure 4. Uncertainty-Dependent Fluence Determination | |||
38 | |||
Fluence Measurement Methods | |||
(2.1.3)* (2.2) | |||
Dosimeter Dosimeter Periodic Validation Certified Fluence & | |||
Measurement Uncertainty in Standard and/or Spectrum-Averaged Uncertainties Analysis Reference Fields Cross-Section | |||
& Biases Ma , BMa vM , BMv Combine Uncertainty & | |||
Bias Estimates Measurement Uncertainty Based on Analytic (Validation) | |||
M M | |||
M | |||
, ( ) | |||
a M | |||
v - Uncertainty Analysis | |||
,B | |||
Measurement Bias Based on M | |||
B , (B ) | |||
a M | |||
v - Analytic (Validation) Uncertainty Analysis | |||
* Relevant Regulatory Positions are indicated. | |||
Estimated Measurement M, (BM ) - Uncertainty (Bias) | |||
Figure 5. Measurement Qualification Procedure | |||
REFERENCES | |||
1. J.F. Carew et al., "Application of Neutron Transport Green's Functions to the Calculation of Pressure Vessel Fluence," Nuclear Science and Engineering, Vol. 91, p. 279, 1985. | |||
2. Standard Guide for Characterizing Neutron Exposures in Iron and Low Alloy Steels in Terms of Displacements per Atom (DPA), ASTM E693, American Society for Testing and Materials, Philadelphia, 1994.1 | |||
3. Guide to the Expression of Uncertainty in Measurement, International Organization for Standardization, ISBN 92-67-10188-9, Genéve, 1993. [ISO, Care Postal 56, CH-1211 Genéve 20, Switzerland]. | |||
4. J.F. Carew et al., "Pressure Vessel Fluence Calculation Benchmark Problems and Solutions," Draft NUREG/CR-6115 (BNL NUREG-52395), USNRC, September 1999.2 | |||
5. R.E. Maerker et al., "Applications of the LEPRICON Unfolding Procedures to the Arkansas Nuclear One-Unit 1 Reactor," Nuclear Science and Engineering, Vol. 93(2), pp. 137-170, | |||
June 1986. | |||
6. R.E. Maerker, "LEPRICON Analysis of Pressure Vessel Surveillance Dosimetry Inserted into H.B. Robinson-2 During Cycle 9," Nuclear Science and Engineering, Vol. 96(4), pp. | |||
263-289, August 1987. | |||
7. R.E. Maerker, "Analysis of the NESDIP2 and NESDIP3 Radial Shield and Cavity Experiments," NUREG/CR-4886 (Oak Ridge National Laboratory, ORNL/TM-10389), | |||
USNRC, May 1987.3 | |||
8. C.Y. Fu and D.M. Hetrick, "Update of ENDF/B-V Mod-3 Iron: Neutron-Producing Reaction Cross-Sections and Energy Angle Correlations," ORNL/TM-9964, ENDF-341, Oak Ridge National Laboratory, July 1986.4 | |||
9. Standard Guide for Application of ASTM Evaluated Cross Section Data File, Matrix E | |||
706 (IIB), ASTM E1018-95, American Society for Testing and Materials, Philadelphia, | |||
1995.1 | |||
1 Available from the American Society for Testing and Materials, 1916 Race Street, Philadelphia, PA 19103-1187. | |||
2 Copies are available for inspection or copying for a fee from the NRC Public Document Room at 11555 Rockville Pike, Rockville, MD; the PDR's mailing address is USNRC PDR, Washington, DC 20555; telephone (301)415-4737; fax (301)415- | |||
3548. | |||
3 Copies may be purchased at current rates from the U.S. Government Printing Office, Post Office Box 37082, Washington, DC | |||
20013-7082 (telephone (202)512-2249 or (202)512-2171); or from the National Technical Information Service by writing NTIS | |||
at 5285 Port Royal Road, Springfield, VA 22161. Copies are available for inspection or copying for a fee from the NRC Public Document Room at 11555 Rockville Pike, Rockville, MD; the PDR's mailing address is USNRC PDR, | |||
Washington, DC 20555; telephone (301)415-4737; fax (301)415-3548. | |||
4 Available from the Radiation Safety Information Computational Center, Oak Ridge National Laboratory, Post Office Box 2008, Oak Ridge, TN 37831-6362. | |||
40 | |||
4 | 10. B. Petrovic, H.L. Hanshaw, and A. Haghighat, Evaluation of Anisotropy Effects in Pressure Vessel Fluence Calculations Using the BUGLE-93 Library, Transactions of the ANS, Vol. 71, 1994.5 | ||
11. R.W. Roussin et al., "VITAMIN-C: The CTR Processed Multigroup Cross-Section Library for Neutronics Studies, ORNL/RSIC-37 (ENDF-296), Oak Ridge National Laboratory, July 1980.4 | |||
12. R.W. Roussin et al., "VITAMIN-E: 174 Neutron, 38 Gamma-Ray Multigroup Cross- Section Library for Deriving Application-Dependent Working Libraries for Radiation Transport Calculations," DLC-113, Oak Ridge National Laboratory, November 1987.4 | |||
13. W.E. Ford, III, et al., "Modification Number One to the 100n - 21g Cross Section Library," | |||
ORNL/TM-5249 (Available as DLC-37D/EPR from RSIC), Oak Ridge National Laboratory, March 1976.4 | |||
14. "SAILOR: Coupled, Self-Shielded, 47-Neutron, 20-Gamma Ray, P3, Cross Section Library for Light Water Reactors," DLC-76-SAILOR, Oak Ridge National Laboratory, July 1987.4 (Available at <http://epicws.epm.ornl.gov/rsic.html> .) | |||
15. M.L. Williams et al., "The ELXSIR Cross-Section Library for LWR Pressure Vessel Irradiation Studies: Part of the LEPRICON Computer Code System," NP-3654, Electric Power Research Institute, September 1984.6 | |||
16. D.T. Ingersoll et al., Production and Testing of the VITAMIN-B6 Fine-Group and the BUGLE-93 Broad-Group Neutron/Photon Cross-Section Libraries Derived from ENDF/B- | |||
VI Nuclear Data, ORNL-6795, NUREG/CR-6214, January 1995; and J.E. White et al., | |||
BUGLE-96: A Revised Multigroup Cross-Section Library for LWR Applications Based on ENDF/B-VI Release 3, presented at the American Nuclear Society Radiation & Shielding Topical Meeting, April 21-25, 1996, Falmouth, MA, April 1996.3 | |||
17. R.E. Maerker, M.L. Williams, and B.L. Broadhead, "Accounting for Changing Source Distributions in Light Water Reactor Surveillance Dosimetry Analysis," Nuclear Science and Engineering, Vol. 94, pp. 291-308, 1986. | |||
18. R.E. Maerker, M.L. Williams, and B.L. Broadhead, "TIMEPATCH: A Module in the LEPRICON Computer Code System for Evaluating Effects of Time-Dependent Source Distributions in PWR Surveillance Dosimetry," EPRI Interim Report, December 1985. | |||
(Available from the Radiation Safety Information Computational Center, ORNL, as part of PSR-277/LEPRICON, or at <http://epicws.epm.ornl.gov/rsic.html> .)3 | |||
5 Available from the American Nuclear Society, 555 N. Kensington Avenue, La Grange Park, Illinois 60525. | |||
6 Available from EPRI Research Reports Center, P.O. Box 50490, Palo Alto, CA 94303. | |||
41 | |||
19. M. Todosow and J.F. Carew, "Evaluation of Selected Approximations Used in Pressure Vessel Fluence Calculations," Transactions of the American Nuclear Society, Vol. 46, p. | |||
658, June 1984.5 | |||
20. W.J. Eich, "Advanced Recycle Methodology Program," Part II, Chapter 5, Research Project | |||
118-1, Electric Power Research Institute, 1976.7 | |||
21 A. Ahlin et al., "CASMO - A Fuel Assembly Burnup Program," AE-RF-76-4158 (Rev. | |||
Ed.), Studsvik Energiteknik AB, 1978. (Available from Studsvik of America Inc., 1087 Beacon St., Newton, MA 02159). | |||
22. A.L. Aronson et al., "Evaluation of Methods for Reducing Pressure Vessel Fluence," BNL- | |||
NUREG-32876, Brookhaven National Laboratory, March 1983.2 | |||
23. G.P. Cavanaugh et al., "Reduction in Reactor Vessel Irradiation Through Fuel Management," Transactions of the American Nuclear Society, Vol. 45, p. 98, October | |||
1983.5 | |||
24. M. Todosow et al., "Pressure Vessel Fluence Reduction Through Selective Fuel Assembly Replacement," Transactions of the American Nuclear Society, Vol. 45, p. 595, October | |||
1983.5 | |||
25. D. Cokinos et al., "Pressure Vessel Damage Fluence Reduction by Low-Leakage Fuel Management," Transactions of the American Nuclear Society, Vol. 45, p. 594, October | |||
1983.5 | |||
26. M.L. Williams, "DOTSOR: A Module in the LEPRICON Computer Code System for Representing the Neutron Source Distribution in LWR Cores," EPRI Interim Report, December 1985. (Available from the Radiation Safety Information Computational Center, ORNL, as part of PSR-277/LEPRICON, or at <http://epicws.epm.ornl.gov/rsic.html> .)3 | |||
27. A. Haghighat, M. Mahgerefteh, and B. Petrovic, Evaluation of the Uncertainties in the Source Distribution for Pressure Vessel Neutron Fluence Calculations, Nuclear Technology, Vol. 109, 54-75, January 1995. | |||
28. "DOT 3.5 - A Two-Dimensional Discrete Ordinate Transport Code," CCC-276, Oak Ridge National Laboratory, 1978.4 In RSICC Computer Code Collection, CCC-650, Doors3.2, ORNL, July 1998. (Available at <http://epicws.epm.ornl.gov/rsic.html> .) | |||
29. W.A. Rhoades and R.L. Childs, "An Updated Version of the DOT 4 One-and-Two- Dimensional Neutron/Photon Transport Code," ORNL-5851, Oak Ridge National Laboratory, July 1982.4 In RSICC Computer Code Collection, CCC-650, Doors3.2, ORNL, July 1998. (Available at <http://epicws.epm.ornl.gov/rsic.html> .) | |||
7 Copies are available from the National Technical Information Service by writing NTIS at 5285 Port Royal Road, Springfield, VA 22161; telephone (703)487-4650; <http://www.ntis.gov/ordernow>. | |||
42 | |||
30. F.B. Kam et al., "Pressure Vessel Fluence Analysis and Neutron Dosimetry," | |||
NUREG/CR-5049 (Oak Ridge National Laboratory, ORNL/TM-10651), December 1987.3 | |||
31. "Standard Guide for Application of Neutron Transport Methods for Reactor Vessel Surveillance, E 706 (IID)," ASTM E482-89, American Society for Testing and Materials, Philadelphia, 1989.1 | |||
32. B.G. Petrovic and A. Haghighat, Effects of SN Method Numerics on Pressure Vessel Neutron Fluence Calculations, Nuclear Science and Engineering, 122, 167-193, February | |||
1996. | |||
33. B.G. Petrovic and Haghighat, Analysis of Inherent Oscillations in Multidimensional SN | |||
Solutions of the Neutron Transport Equation, Nuclear Science and Engineering, 124, 31- | |||
37, September 1996. | |||
34. B.G. Petrovic and A.Haghighat, New Directional Theta-Weighted (DTW) Differencing Scheme and Reduction of Estimated Pressure Vessel Fluence Uncertainty, Proceedings of the Ninth International Symposium on Reactor Dosimetry, H. Ait Abderrahim, P. Dhondt and B. Osmera, Eds., World Scientific Publ. Co., 1998. | |||
35. P. G. Laky and N. Tsoulfanidis, Neutron Fluence at the Pressure Vessel of a Pressurized Water Reactor Determined by the MCNP Code, Nuclear Science & Engineering, 121, 433, | |||
1995. | |||
36. J. C. Wagner, A. Haghighat, and B. G. Petrovic, Monte Carlo Transport Calculations and Analysis for Reactor Pressure Vessel Neutron Fluence, Nuclear Technology, 114, No. 3, | |||
373, 1996. | |||
37. J. C. Wagner, A. Haghighat, and B. G. Petrovic, Investigation of Pressure Vessel Fluence Calculation with Monte Carlo, Transactions of the American Nuclear Society, 68, Part A, | |||
446, June 1993. | |||
38. W. T. Urban et al., PCA Benchmark Solutions Using MCNP and THREEDANT, LA- | |||
UR-93-31-21, Los Alamos National Laboratory, 1993. In ASTM STP 1228, pp. 376-383, | |||
1994. | |||
39. A. Avery et al., Calculations of Pressure Vessel Fluence in PWRs Using ENDF/B-VI | |||
Data, Proceedings of the 8th International Conference on Radiation Shielding, Arlington, Texas, April 24-28, 1994, Vol. 2, p. 677, American Nuclear Society, 1994.5 | |||
40. S. Power, An Analysis of the HB Robinson Unit 2 PWR Using the Monte Carlo Code MCBEND, Proc. of the 7th ASTM-Euratom Symposium on Reactor Dosimetry, Strasbourg ASTM, 1990.1 | |||
43 | |||
41. J.C. Wagner and A. Haghighat, Automated Variance Reduction of Monte Carlo Using Discrete Ordinates Adjoint Functions, Nuclear Science & Engineering, Vol. 128, 186-208, | |||
1998. | |||
42. A. Haghighat et al., Performance of the Automated Adjoint Accelerated MCNP(A3MCNPTM1) for Simulation of a BWR Core Shroud Problem, Proceedings of the M&C 99 Meeting, Madrid, Spain, Vol. 2, 1381-1392, Senda Editorial, S.A., Madrid, Spain, September 27, 1999. | |||
43. J. F. Breismeister (Ed.), MCNP - A General Monte Carlo N-Particle Transport Code, Version 4C, LA-12625, Los Alamos National Laboratory, March 1997.3 | |||
44. G. A. Wright et al., MCBEND - A Fluence Tool from AEA Technology, Reactor Dosimetry, ASTM STP 1398, J. G. Williams et al., Eds., ASTM, West Conshohoken, PA. | |||
2000.1 | |||
45. T. E. Booth, A Sample Problem in Variance Reduction in MCNP, LA-10363-MS, Los Alamos National Laboratory, 1985. | |||
46. L. L. Carter and E.D. Cashwell, Particle Transport Simulation with the Monte Carlo Method, ERDA Critical Review Series, TID-26607, 1975.7 | |||
47. J.F. Carew, D.K. Min, and A.L. Aronson, "Spectral Effects in the Extrapolation of Pressure Vessel Surveillance Capsule Measurements," Nuclear Technology, Vol. 55, No. 3, p. 565, December 1981. | |||
48. M.L. Williams, P. Chowdhury, and B.L. Broadhead, DOTSYN: A Module for Synthesizing Three-Dimensional Fluxes in the LEPRICON Computer Code System, EPRI Interim Report, December 1985. (Available from the Radiation Safety Information Computational Center, ORNL, as part of PSR-277/LEPRICON, or at <http://epicws.epm.ornl.gov/rsic.html> .)3 | |||
49. R.E. Maerker, "Analysis of the VENUS-3 Experiments," NUREG/CR-5338 (Prepared for the NRC by Oak Ridge National Laboratory, ORNL/TM-11106), USNRC, August 1989.2 | |||
50. L. Lois and T. Collins, "Reactor Cavity Dosimetry - A Regulatory Perspective," Transactions of the American Nuclear Society, Vol. 63, p. 431, June 1991.4 | |||
51. N. Tsoulfanidis, "Neutron Energy Spectra in the Core and Cavity of the ANO-2 PWR," NP- | |||
4238, Electric Power Research Institute, September 1985.5 | |||
52. B.N. Taylor and C.E. Kuyatt, "Guidelines For Evaluating and Expressing the Uncertainty of NIST Measurement Results," NIST Technical Note-1297, National Institute For Standards and Technology, January 1993.1 | |||
53. J.F. Carew et al., "Pressure Vessel Fluence Benchmark Calculations," BNL-NUREG-34715, Brookhaven National Laboratory, February 1984.3 | |||
44 | |||
54. D.M. Cokinos et al., "Benchmarking of Pressure Vessel Fluence Calculations," Transactions of the American Nuclear Society, Vol. 46, p. 636, June 1984.4 | |||
55. B.L. Broadhead et al., "LEPRICON Adjustment Module: A Generalized Linear Least Squares Data Analysis Program with Application to PWR Surveillance Dosimetry," EPRI Interim Report, March 1985. (Available from the Radiation Safety Information Computational Center, ORNL, as part of PSR-277/LEPRICON, or at <http://epicws.epm.ornl.gov/rsic.html> .)3 | |||
. | |||
56. R.E. Maerker, B.L. Broadhead, and J.J. Wagschal, "Theory of a New Unfolding Procedure in Pressurized Water Reactor Pressure Vessel Dosimetry and Development of an Associated Benchmark Data Base," Nuclear Science and Engineering, Vol. 91(4), pp. 369-392, December | |||
1985. | |||
57. W.N. McElroy et al., "A Computer-Automated Iterative Method for Neutron Flux Spectral Determination by Foil Activation," AFL-TR-67-41, Vol. I (available as RSIC Program No. | |||
CCC-112/SAND II), 1967.3 | |||
58. C.A. Oster et al., "A Modified Monte Carlo Program for SAND-II with Solution Weighing and Error Analysis," Hanford Engineering Development Laboratory, HEDL-TME 76-60, 1976.2 | |||
59. C.R. Green, J.A. Halbleib, and J.V. Walker, "A Technique for Unfolding Neutron Spectra from Activation Measurements," Sandia Corporation, SC-RR-67-746, 1967. (Available as RSIC Program No. CCC-108/SPECTRA.)4 | |||
60. F.G. Perey, "Least-Squares Dosimetry Unfolding: The Program STAY'SL," ORNL/TM-6062, Oak Ridge National Laboratory, 1977. (Available as RSIC Program No. PSR-113.)4 | |||
61. F. Schmittroth, "FERRET Data Analysis Code," HEDL-TME 79-40, Hanford Engineering Development Laboratory, 1979. (Available as RSIC Program No. PSR-145.)4 | |||
62. F.W. Stallmann, "LSL-M2: A Computer Program for Least-Squares Logarithmic Adjustment of Neutron Spectra," NUREG/CR-4349 (Prepared for NRC by Oak Ridge National Laboratory, ORNL/TM-9933), March 1986. (Available as RSIC Program No. PSR-233.)3,4 | |||
63. "Standard Guide for Application of Neutron Spectrum Adjustment Methods in Reactor Surveillance (IIA), ASTM E944-96, American Society for Testing and Materials, Philadelphia, 1996.1 | |||
64. "Standard Guide for Benchmark Testing of Reactor Dosimetry in Standard and Reference Neutron Fields, ASTM E2005-99, American Society for Testing and Materials, Philadelphia, | |||
1999.1 | |||
65. I. Remec and F.B. Kam, "H. B. Robinson-2 Pressure Vessel Benchmark," NUREG/CR-6453 (Prepared for NRC by Oak Ridge National Laboratory, ORNL/TM-13204), February 1998.2 | |||
45 | |||
66. M.L. Williams and M. Asgari, Impact of ENDF/B-VI Cross-Section Data on H.B. Robinson Cycle 9 Dosimetry Calculations, NUREG/CR-6071, October 1993.3 | |||
67. E.P. Lippincott et al., "Evaluation of Surveillance Capsule and Reactor Cavity Dosimetry from H.B. Robinson Unit-2, Cycle 9," NUREG/CR-4576 (WCAP-11104), USNRC, February 1987.2 | |||
68. M. Austin, Sense of Direction: An Observation of Trends in Materials Dosimetry in the United Kingdom, Proceedings of the Fourth ASTM Euratom Symposium on Reactor Dosimetry, Gaithersburg, MD, March 22-26, 1982, NUREG/CP-0029, USNRC, Vol 1, August 1982.1 | |||
69. J. Butler et al., The PCA Replica Experiment Part I: Winfrith Measurements and Calibrations, AEEW-R 1736, Part I, UKAEA, Winfrith, United Kingdom, January 1984.7 | |||
70. M. D. Carter and I. J. Curl, NESTOR Shielding and Dosimetry Improvement Programme, AEEW-M 2329, 1986. | |||
71. J. Butler et al., Review of the NESTOR Shielding and Dosimetry Improvement Programme (NESDIP), Reactor Dosimetry: Methods, Applications and Standardization, Sixth ASTM | |||
Euratom Symposium on Reactor Dosimetry, Jackson Hole, Wyoming, May 31- June 6, 1987, STP 1001, ASTM, May 1989.1 | |||
72. W.N. McElroy, "LWR Pressure Vessel Surveillance Dosimetry Improvement Program: PCA | |||
Experiments and Blind Test," NUREG/CR-1861 (Prepared for NRC by Hanford Engineering Development Laboratory, HEDL-TME 80-87), July 1981.2 | |||
73. I. Remec and F.B. Kam, "Pool Critical Assembly Pressure Vessel Facility Benchmark," | |||
NUREG/CR-6454 (Prepared for NRC by Oak Ridge National Laboratory, ORNL/TM-13205), | |||
July 1997.2 | |||
74. F.W. Stallman et al., Reactor Calculation `Benchmark PCA Blind Test Results, ORNL/NUREG/TM-428, March 1981.2 | |||
75. D. K. Min, A. L. Aronson, and J. F. Carew, Analysis of the ORNL Pool Critical Assembly Pressure Vessel Dosimetry Benchmark Experiment, BNL-NUREG-29047, Brookhaven National Laboratory, February 1981.3 | |||
76. P. D'hondt et al., "Contributions of the Venus-Engineering Mock-Up Experiments to the LWR- | |||
PV Surveillance," Proceedings of the 7th ASTM-Euratom Symposium on Reactor Dosimetry, Strasbourg, France, ASTM, August 1990.1 | |||
77. L. Leenders, LWR-PVS Benchmark Experiment VENUS-3 with Partial Length Shield Assemblies, SCK-CN, MOL, Belgium, FCP/VEN/01, September 1998. (See also Protection of Neutron Embrittlement in the Reactor Pressure Vessel: Venus-1 and Venus-3 Benchmarks, Nuclear Energy Agency, 2000. (Available at WWW.NEA.FR)) | |||
46 | |||
78. G. Hehn and B.C. Na, New NEA Benchmarks Reveal Decisive Improvements in Calculating Fast Neutron Fluence for Predictions of Embrittlement in Reactor Pressure Vessels, Reactor Dosimetry, ASTM STP 1398 (J.G. Williams et al., Eds.), ASTM, West Conshohoken, PA. | |||
2000.1 | |||
79. A. Haghighat, H. Ait Abderrahim, and G.E. Sjoden, Accuracy and Performance of PENTRANTM Using the VENUS-3 Benchmark Experiment, Reactor Dosimetry, ASTM STP | |||
1398 (J.G. Williams et al., Eds.), ASTM, West Conshohoken, PA. 2000.1 | |||
80. R.E. Maerker and B.A. Worley, "Activity and Fluence Calculations for the Startup and Two- Year Irradiation Experiments Performed at the Poolside Facility," NUREG/CR-3886 (Prepared for NRC by Oak Ridge National Laboratory, ORNL/TM-4265), October 1984.2 | |||
81. W.N. McElroy, "LWR Pressure Vessel Surveillance Dosimetry Improvement Program. PSF | |||
Experiments Summary and Blind Test Results," NUREG/CR-3320, Vols. 1-4 USNRC, | |||
(Prepared for the NRC by Hanford Engineering Development Laboratory), July 1986 through July 1992.2 | |||
82. G. L. Guthrie, E. P. Lippincott, E. D. McGarry, Light Water Reactor Pressure Vessel Surveillance Dosimetry Improvement Program: PSF Blind Test Workshop Minutes, Westinghouse Hanford Company, April 1994.1 | |||
83. "Standard Test Method for Measuring Reaction Rates by Radioactivation of Neptunium-237," | |||
ASTM E 705-90, ASTM, Philadelphia, 1991.1 | |||
84. "Standard Test Method for Measuring Fast-Neutron Reaction Rates by Radioactivation of Niobium," ASTM E1297-89, ASTM, Philadelphia, 1989.1 | |||
85. "Standard Test Method for Measuring Reaction Rates by Radioactivation of Uranium-238," | |||
ASTM E 704-90, ASTM, Philadelphia, 1991.1 | |||
86. "Standard Test Method for Measuring Fast-Neutron Reaction Rates by Radioactivation of Nickel," ASTM E 264-87, ASTM, Philadelphia, 1987.1 | |||
87. "Standard Test Method for Measuring Fast-Neutron Reaction Rates by Radioactivation of Iron," | |||
ASTM E 263-88, ASTM, Philadelphia, 1988.1 | |||
88. "Standard Test Method for Measuring Fast-Neutron Reaction Rates by Radioactivation of Titanium," ASTM E 526-87, American Society for Testing and Materials, Philadelphia, 1987.1 | |||
89. "Standard Test Method for Measuring Fast-Neutron Reaction Rates by Radioactivation of Copper," ASTM E 523-87, ASTM, Philadelphia, 1987.1 | |||
90. "Standard Guide for Sensor Design and Irradiation for Reactor Surveillance, E 706 (IIC)," | |||
ASTM E 844-86, ASTM, Philadelphia, 1986.1 | |||
47 | |||
91. "Standard Test Method for Application and Analysis of Radiometric Monitors for Reactor Vessel Surveillance, E 706 (IIIA)," ASTM E 1005-84, ASTM, Philadelphia, 1984.1 | |||
92. "Standard General Methods for Detector Calibration and Analysis of Radionuclides," ASTM E | |||
181-82, ASTM, Philadelphia, 1983.1 | |||
93. R. Gold et al., "Neutron Dosimetry with Solid-State Track Recorders in the Three-Mile Island Unit-2 Reactor Cavity," Nuclear Tracks, Vol. 103, p. 447, 1985. | |||
94. F.H. Ruddy et al., "Solid-State Track Recorder Neutron Dosimetry in Light-Water Reactor Pressure Vessel Surveillance Mockups in Reactor Dosimetry,"Proceedings of the 5th ASTM- | |||
EURATOM Symposium on Reactor Dosimetry, Geesthacht, Federal Republic of Germany, Sept. | |||
24-28, 1984, EUR 9869, Commission of the European Communities, D. Reidel Publishing Co., | |||
1985.1 | |||
95. "Standard Test Method for Application and Analysis of Solid State Track Recorder (SSTR) | |||
Monitors for Reactor Surveillance, E 706 (IIIB)," ASTM E 854-90, ASTM, Philadelphia, 1990.1 | |||
96. "Standard Test Method for Application and Analysis of Helium Accumulation Fluence Monitors for Reactor Vessel Surveillance, E 706(IIIC)," ASTM E 910 95, ASTM, Philadelphia, 1995.1 | |||
97. J.A. Grundl and C.M. Eisenhauer, Compendium of Benchmark Neutron Fields for Reactor Dosimetry, NBSIR 85-3151, National Bureau of Standards, Gaithersburg, MD, January 1986.8 | |||
98. A. Hawari et al., "Materials Dosimetry Reference Facility," Proceedings of the Eighth ASTM- | |||
Euratom Symposium on Reactor Dosimetry, Vail, Colorado, August 1993.1 | |||
99. "Standard Guide for Benchmark Testing of Light Water Reactor Caalculations, ASTM E2006- | |||
99, ASTM, Philadelphia, 1999.1 | |||
100. "Standard Practice for Determining Neutron Fluence, Fluence Rate and Spectra by Radioactivation Techniques, ASTM E261-98, ASTM, Philadelphia, 1998.1 | |||
8 Available from the National Institute of Standards and Technology, NIST Publication Productions, Room A635, Gaithersburg, MD 20899. | |||
48 | |||
REGULATORY ANALYSIS | |||
1. STATEMENT OF THE PROBLEM | |||
The U.S. Nuclear Regulatory Commission (NRC) has promulgated regulations to ensure the structural integrity of the reactor pressure vessel for light water power reactors. Specific fracture toughness requirements for normal operation and for anticipated operational occurrences for power reactors are set forth in Appendix G, "Fracture Toughness Requirements," to 10 CFR Part 50, | |||
"Domestic Licensing of Production and Utilization Facilities." These requirements are imposed through 10 CFR 50.60. Additionally, in response to concerns over potential pressurized thermal shock (PTS) events in pressurized water reactors (PWRs), the NRC issued 10 CFR 50.61, | |||
"Fracture Toughness Requirements for Protection Against Pressurized Thermal Shock Events." | |||
To satisfy the requirements of Appendix G; 10 CFR 50.61; Criterion 14, Reactor Coolant Pressure Boundary, Criterion 30, Quality of Reactor Coolant Pressure Boundary, and Criterion | |||
31, Fracture Prevention of Reactor Coolant Pressure Boundary, of Apppendix A, General Design Criteria for Nuclear Power Plants, to 10 CFR Part 50, methods for determining the fast neutron fluence (E > MeV) are necessary to estimate the fracture toughness of the pressure vessel materials. Appendix H, "Reactor Vessel Material Surveillance Program Requirements," to 10 CFR | |||
Part 50 requires the installation of surveillance capsules, including material test specimens and flux dosimeters, to provide data on material damage correlations as a function of fluence. | |||
The neutron fluence is attenuated by several decades between the core and vessel. This attenuation results in a strong sensitivity of the calculated vessel fluence to the physical description of the core and vessel internals and the numerical calculation of the neutron transport, and it makes an accurate determination of the pressure vessel fluence difficult. As a result, a wide range of methods of varying reliability and accuracy have been used to determine the reactor vessel fluence. | |||
Consequently, comparisons of measured and calculated fluences have shown varying degrees of agreement, and in some cases conservatisms have been required in licensing analyses to accommodate the observed measurement-to-calculation differences. | |||
Over the past decade, substantial improvements have been made in both the calculation and measurement of the pressure vessel fluence. These improvements have stemmed from both NRC | |||
and industry programs. These include the development and improvement of computer codes, calculational models, measurement techniques and basic cross-section data, and the systematic qualification of the fluence methods by comparison to NRC-sponsored benchmark experiments. | |||
These calculation and measurement improvements provide increased accuracy in the fluence determinations that are an essential part of meeting the requirements of Appendices G to 10 CFR | |||
Part 50 and 10 CFR 50.61. This is especially important for plants seeking to renew their operating licenses. | |||
The wide variation in fluence calculation methods has resulted in lengthy plant-specific reviews and made it difficult to confirm, during the review process, that the actual fluence is adequately bounded by the various calculational methods used. This calculation and dosimetry | |||
49 | |||
guide would provide standardized methods and procedures that would allow these reviews to be greatly simplified, and would improve confidence in the calculated fluence values. | |||
===2. OBJECTIVE=== | |||
The objective of this guide is to provide state-of-the-art calculation and measurement procedures that are acceptable to the NRC staff for determining pressure vessel fluence. These procedures would yield a more accurate and reliable vessel fluence determination than the procedures that have been used in the past. The improved accuracy and realistic assessment of the uncertainty in the calculation would provide assurance that the fluence value is appropriate for use in evaluating compliance with the regulations. | |||
1 | ===3. ALTERNATIVES=== | ||
The alternatives to issuing the vessel fluence calculation and dosimetry regulatory guide are as follows: | |||
3.1 Branch Technical Position The pressure vessel fluence methods provided by the regulatory guide could be included in a branch technical position. However, this is not considered an acceptable alternative since the branch position does not provide the same high level of input and participation, especially public, industry, and the ACRS input, that the vessel fluence analysis requires. | |||
3.2 NUREG-Series Report The vessel fluence procedures could be published in a NUREG-series report. However, these reports also do not receive the required input and participation from the public, industry, and the ACRS and they are not appropriate for providing regulatory guidance. | |||
3.3 Discussions with Licensees The detailed fluence calculational and measurement methods that are considered acceptable to the NRC staff could be provided to the licensees through individual reviews and discussions on a case-by-case basis. This alternative is basically the same as the current practice and is equivalent to taking no action. Individual licensee discussions are extremely time-consuming for both the NRC staff and the licensee, they lead to highly individual analyses and reviews, and they do not result in an established standard. | |||
===4. COSTS AND BENEFITS=== | |||
4.1 Benefits The methods described in this guide may be used for all fluence determinations used in vessel fracture toughness evaluations, including the determination of the fluence used in calculating the pressure vessel material values of RTNDT specified in Appendix G to 10 CFR Part | |||
50 imposed by 10 CFR 50.60, and the values of RTPTS specified in 10 CFR 50.6 | |||
===1. The regulatory=== | |||
50 | |||
guide would improve the accuracy and reliability of these evaluations and provide consistent and reliable uncertainty estimates by incorporating state-of-the-art methods and procedures for determining the fluence and the fluence uncertainty. The guide would also ensure the completeness of licensee vessel fluence submittals and improve the efficiency of staff reviews. | |||
The improved fluence determination will provide fluence and uncertainty estimates that are more reliable and understandable. Thus, it will provide more reliable and accurate information for the PTS screening criteria of 10 CFR 50.61 and improve application of this rule. In this regard it is noted that, for a pressure vessel near the PTS screening criteria of 10 CFR 50.61, a 25% reduction in calculated end-of-license fluence, which is typical of existing uncertainties, will reduce the calculated vessel failure frequency by approximately a factor of three. | |||
4.2 NRC Costs The NRC costs for reviewing fluence-related submittals would be reduced substantially by the issuance of this guide. For estimating the costs, it is assumed that of the ~60 PWRs, half have a RTPTS within ~40-F of the PTS screening criterion (or other temperature limit) at the end of license and will require a detailed review. Assuming each submittal requires a staff week and only half of the PWR owners submit revised fluence analyses, the total NRC cost is approximately 15 staff weeks. If the licensees used the methods given in the guide, this cost could be reduced to approximately 2 staff weeks. | |||
4.3 Licensee Costs Increased costs to the licensee would result from changes in the fluence calculation and measurement procedures. The calculational costs would be one-time costs and have been estimated in Table RA-1. The licensee costs resulting from the changes in the measurement procedures have been estimated in Table RA-2. | |||
===5. DECISION RATIONALE=== | |||
It is recommended that the proposed regulatory guide be issued because (1) a high level of participation of the NRC, ACRS, industry, and the public is reflected in the guide, (2) a methodology standard would be established, and (3) inefficient use of NRC staff and licensee resources during the review process would be eliminated. | |||
The alternatives identified above for providing acceptable fluence methods to the licensees do not provide the advantages listed above associated with issuance of this regulatory guide. In particular, while these alternative approaches may result in the same or slightly increased cost to licensees, they result in a highly inefficient use of the NRC staff resources. The alternatives to a regulatory guide are therefore judged to be unacceptable. | |||
51 | |||
a | TABLE RA-1. ADDITIONAL LICENSEE CALCULATION COSTS | ||
Calculation Tasks Staff Weeks (a) Modifications to calculational models +2 (b) Additional calculation benchmarking and +8 qualification (c) Calculation uncertainty analysis +3 (d) Calculation documentation and reporting +2 (e) Reduced Licensing Activities* -4 Total Additional Licensee Cost +11 staff weeks TABLE RA-2. ADDITIONAL LICENSEE MEASUREMENT COSTS | |||
Measurement Task Staff Weeks (a) Additional Quality Control +1 (b) Dosimeter Response Corrections +1 (c) Periodic Detector Calibration +2 (d) Response Uncertainty Analysis +2 (e) Additional Measurement +1 Documentation and Reporting Total Additional Licensee Cost + 7 staff weeks | |||
* This estimate reflects the reduced licensee costs from not attending meetings with the NRC staff and responding to questions as a result of following the procedures in the regulatory guide. | |||
52 | |||
BACKFIT ANALYSIS | |||
The regulatory guide does not require a backfit analysis as described in 10 CFR 50.109(c) | |||
because it does not impose a new or amended provision in the Commission rules or a regulatory staff position interpreting the Commission rules that is either new or different from a previous applicable staff position. In addition, this regulatory guide does not require the modification or addition to systems, structures, components, or design of a facility or the procedures or organization required to design, construct, or operate a facility. Rather a licensee or applicant can select a preferred method for achieving compliance with license or the rules or the orders of the Commission as described in 10 CFR 50.109(a)(7). The regulatory guide provides the opportunity to use the methods described in the guide for all fluence determinations used in vessel fracture toughness evaluations, including the determination of the fluence used in calculating the pressure vessel materials' values of RTNDT for use in 10 CFR Part 50, Appendix G (imposed by 10 CFR | |||
50.60), and the values of RTPTS in accordance with 10 CFR 50.61.}} | |||
{{RG-Nav}} | {{RG-Nav}} |
Latest revision as of 05:55, 24 November 2019
ML010890301 | |
Person / Time | |
---|---|
Issue date: | 03/31/2001 |
From: | Office of Nuclear Regulatory Research |
To: | |
Jones W (301)415-7558 | |
References | |
-nr, DG-1025, DG-1053 RG-1.190 | |
Download: ML010890301 (56) | |
U.S. NUCLEAR REGULATORY COMMISSION March 2001 REGULATORY
GUIDE
OFFICE OF NUCLEAR REGULATORY RESEARCH
REGULATORY GUIDE 1.190
(Previous drafts were DG-1053 and DG-1025)
CALCULATIONAL AND DOSIMETRY METHODS
FOR DETERMINING
PRESSURE VESSEL NEUTRON FLUENCE
Regulatory guides are issued to describe and make available to the public such information as methods acceptable to the NRC staff for implementing specific parts of the NRCs regulations, techniques used by the staff in evaluating specific problems or postulated accidents, and data needed by the NRC staff in its review of applications for permits and licenses. Regulatory guides are not substitutes for regulations, and compliance with them is not required. Methods and solutions different from those set out in the guides will be acceptable if they provide a basis for the findings requisite to the issuance or continuance of a permit or license by the Commission.
This guide was issued after consideration of comments received from the public. Comments and suggestions for improvements in these guides are encouraged at all times, and guides will be revised, as appropriate, to accommodate comments and to reflect new information or experience. Written comments may be submitted to the Rules and Directives Branch, ADM, U.S. Nuclear Regulatory Commission, Washington, DC 20555-0001.
Regulatory guides are issued in ten broad divisions: 1, Power Reactors; 2, Research and Test Reactors; 3, Fuels and Materials Facilities; 4, Environmental and Siting;
5, Materials and Plant Protection; 6, Products; 7, Transportation; 8, Occupational Health; 9, Antitrust and Financial Review; and 10, General.
Single copies of regulatory guides (which may be reproduced) may be obtained free of charge by writing the Distribution Services Section, U.S. Nuclear Regulatory Commission, Washington, DC 20555-0001, or by fax to (301)415-2289, or by email to DISTRIBUTION@NRC.GOV. Electronic copies of this guide are available on the internet at NRCs home page at <WWW.NRC.GOV> in the Reference Library under Regulatory Guides. This guide is also in the Electronic Reading Room at NRCs home page, along with other recently issued guides, Accession Number ML010890301.
TABLE OF CONTENTS
Page
A. INTRODUCTION
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
B. DISCUSSION
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
C. REGULATORY POSITION
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1. NEUTRON FLUENCE CALCULATIONAL METHODS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.1 Input Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.2 Core Neutron Source . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.3 Fluence Calculation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4 Methodology Qualification and Uncertainty Estimates . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2. NEUTRON FLUENCE MEASUREMENT METHODS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.1 Measurement Procedures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.2 Validation in Standard and Reference Neutron Fields . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.3 Fluence Determination from Detector Measurements . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
2.4 Ex-vessel Dosimetry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
3. REPORTING . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.1 Fluence Methods . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.2 Multigroup Fluences . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.3 Integral Fluences . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.4 Comparisons of Calculation and Measurement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3.5 Specific Activities and Average Reaction Rates . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
D. IMPLEMENTATION
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29 References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
LIST OF TABLES
1. Summary of Regulatory Positions on Calculation and Dosimetry . . . . . . . . . . . . . . . . . . . . . . . . 30
2. Threshold Detectors Recommended for Pressure Vessel Dosimetry . . . . . . . . . . . . . . . . . . . . . . 34 LIST OF FIGURES
1. Discrete Ordinates Calculation Methodology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
2. Monte Carlo Calculation Methodology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
3. Calculation Methodology Qualification Procedure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
4. Fluence Determination . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
5 Measurement Qualification Procedure . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39 iii
A. INTRODUCTION
This regulatory guide has been developed to provide state-of-the-art calculations and measurement procedures that are acceptable to the NRC staff for determining pressure vessel fluence.
The U.S. Nuclear Regulatory Commission (NRC) has promulgated regulations that ensure the structural integrity of the reactor pressure vessel for light-water-cooled power reactors.
Chapter 15, Accident Analysis, of the Standard Review Plan (NUREG-0800, Standard Review Plan for the Review of Safety Analysis Reports for Nuclear Power Plants) assumes the pressure vessel does not fail. Specific fracture toughness requirements for normal operation and for anticipated operational occurrences for power reactors are set forth in Appendix G, "Fracture Toughness Requirements," of 10 CFR Part 50, "Domestic Licensing of Production and Utilization Facilities." The requirements of Appendix G are imposed by 10 CFR 50.60. Additionally, in response to concerns over potential pressurized thermal shock (PTS) events in pressurized water reactors (PWRs), the NRC issued 10 CFR 50.61, "Fracture Toughness Requirements for Protection Against Pressurized Thermal Shock Events."
To satisfy the requirements of both Appendix G and 10 CFR 50.61, methods for determining the fast neutron fluence (E > 1 MeV) are necessary to estimate the fracture toughness of the pressure vessel materials. Appendix H, "Reactor Vessel Material Surveillance Program Requirements," to 10 CFR Part 50 requires the installation of surveillance capsules, including material test specimens and flux dosimeters, to provide data for material damage correlations as a function of fluence.
The fracture toughness of pressure vessel materials is related to a parameter called the material's "reference temperature for nil-ductility transition," or simply reference temperature, and is denoted as RTNDT. The RTNDT is defined in Revision 2 of Regulatory Guide 1.99, "Radiation Embrittlement of Reactor Vessel Materials," by a correlation of the fluence (E > 1 MeV), material chemistry (concentrations of Cu and Ni), initial reference temperature, and margin to account for uncertainties in the correlation and input values. In 10 CFR 50.61, evaluation of the reference temperature based on the best estimate of the fast neutron fluence at the end of the license period is required, and the corresponding reference temperature is termed RTPTS.
This guide is intended to ensure the accuracy and reliability of the fluence determination required by General Design Criteria (GDC) 14, 30, and 31 of Appendix A, General Design Criteria for Nuclear Power Plants, to 10 CFR Part 50. The guide describes methods and assumptions acceptable to the NRC staff for determining the pressure vessel neutron fluence.
These methods are directly applicable to the determination of RTNDT and RTPTS. The licensee may propose alternative methods. Alternative methods will be considered on a plant-specific basis, especially in cases involving unusual plant characteristics or factors that require different methods and assumptions. It is recognized that, when the embrittlement of a reactor vessel material is not significant and there is a large margin to the RTNDT limits, more approximate methods for determining the fluence may be appropriate (Reference 1).
1
Compliance with this guide is not a regulatory requirement of the USNRC. However, if a licensee elects to use this guide to determine pressure vessel neutron fluence, implementation of the guide would not be satisfied unless the licensee complies with certain specific provisions identified in the Regulatory Position of the guide. The use of the following terms is explained to clarify compliance with these regulatory positions.
Must - Necessary provisions if implementation is to be satisfied.
Should - Provisions that are expected to be complied with unless it is not possible because of specific circumstances (for example, data needed to meet the position are not available).
May - Provisions that are acceptable and recommended, but are to be applied at the option of the licensee.
The information collections contained in this regulatory guide are covered by the requirements of 10 CFR Part 50, which were approved by the Office of Management and Budget, approval number 3150-0011. If a means used to impose an information collection does not display a currently valid OMB control number, the NRC may not conduct or sponsor, and a person is not required to respond to, the information collection.
B. DISCUSSION
The methods and assumptions described in this guide are for the calculation and measurement of vessel fluence for core and vessel geometrical and material configurations that are typical of current PWR and boiling water reactor (BWR) designs. This guide does not address the determination of surveillance specimen material properties or the correlation between material properties and neutron fluence. The methodology presented is intended as a best estimate, rather than a bounding or conservative fluence determination. For example, in the RTPTS correlation called for in 10 CFR 50.61, the best-estimate fluence is used to calculate the shift in RTPTS.
Uncertainty in the shift prediction (e.g., from uncertainty in the fluence, chemistry factor, or shift correlation) has been included separately in an explicit margin term. While the E > 1 MeV
fluence has been selected as the exposure parameter for the RTNDT and RTPTS correlations, the procedures described in this guide determine the damage fluence spectrum (from 0.1 to 15 MeV)
and are applicable to other exposure units, such as iron displacements per atom - dpa (Ref. 2).
The determination of the pressure vessel fluence is based on both calculations and measurements; the fluence prediction is made with a calculation, and the measurements are used to qualify the calculational methodology. Because of the importance and the difficulty of these calculations, the methods must be qualified by comparison to measurements to ensure a reliable and accurate vessel fluence determination. In this qualification, measurement-to-calculation comparisons are used to identify biases (i.e., systematic errors) in the calculations and to provide reliable estimates of the fluence uncertainties.1 When the measurement data are of sufficient quality and quantity that they allow a reliable estimate of the calculational bias (i.e., they represent a statistically significant measurement data base), the comparisons to measurement may be used to
1 For a discussion of the terms bias and uncertainty as used in this guide, see Sections 3.2 and 3.3 of Reference 3.
2
(1) determine the effect of the various modeling approximations and any calculational bias and, if appropriate, (2) modify the calculations by applying a correction to account for bias or by model adjustment or both. Typically, plant-specific data alone are not sufficient to determine a calculational bias for an individual plant. As an additional qualification, the sensitivity of the calculation to the important input and modeling parameters must be determined and combined with the uncertainties of the input and modeling parameters to provide an independent estimate of the overall calculational uncertainty.
The prediction of the vessel fluence must be made by an "absolute" fluence calculation in which the transport of the neutron flux is calculated from the core out to the vessel and cavity, rather than a simple spatial extrapolation of the fluence measurements.
The calculations of the pressure vessel fluence consist of the following steps: (1)
determination of the geometrical and material input data, (2) determination of the core neutron source, (3) propagation of the neutron fluence from the core to the vessel and into the cavity, and
(4) qualification of the calculational procedure. These steps are discussed in detail in Regulatory Positions 1.1 through 1.4. In Regulatory Position 2, the use of surveillance dosimetry as an in situ verification for the calculations is described. Reporting is discussed in Regulatory Position 3.
The major regulatory positions are summarized in Table 1, Summary of Regulatory Positions on Calculation and Dosimetry.
As an indication of current practice, selected codes and cross-section libraries are listed in the references; however, it is the responsibility of the licensee to demonstrate their acceptability in a specific application.
C. REGULATORY POSITION
1. NEUTRON FLUENCE CALCULATIONAL METHODS
This guide describes the application and qualification of a methodology acceptable to the NRC staff for determining the best-estimate neutron fluence experienced by materials in the beltline region of light water reactor (LWR) pressure vessels, as well as for determining the overall uncertainty associated with those best-estimate values. Both the deterministic discrete ordinates method and the statistical Monte Carlo transport method are considered acceptable for vessel fluence determination and are described in Regulatory Position 1.3.
Appendix B, "Quality Assurance Criteria for Nuclear Power Plants and Fuel Reprocessing Plants," to 10 CFR Part 50 requires that this methodology be properly qualified. Qualification includes determination of the uncertainty in the reactor vessel fluence as described in Regulatory Position 1.4. The uncertainty of the fluence must be 20% (1 ) or less when the fluence is used to determine RTPTS and RTNDT for complying with 10 CFR 50.61 and Revision 2 of Regulatory Guide 1.99, "Radiation Embrittlement of Reactor Vessel Materials," respectively. It should be recognized that this 20% uncertainty value has been included in the margin term for the RTPTS.
Uncertainty of fluence for other applications should be determined using Regulatory Position 1.4 and included as an uncertainty allowance in the fluence estimate, as appropriate for the specific
3
application. For example, when performing probabilistic risk assessment (PRA) evaluations such as those described in Regulatory Guide 1.154, "Format and Content of Plant-Specific Pressurized Thermal Shock Safety Analysis Reports for Pressurized Water Reactors," an explicit uncertainty term is required in addition to the best-estimate fluence values. In this case, it is acceptable to use the fluence uncertainty determined per the Regulatory Position 1.4 of this guide.
An overview of the calculational methods described in Regulatory Positions 1.1 through
1.3 is provided in Figures 1 and 2. The procedures for qualifying these methods and determining the pressure vessel fluence are summarized in Figures 3 and 4.
1.1 Input Data
1.1.1 Materials and Geometry Detailed material and geometrical input data should be used to define the physical characteristics that determine the attenuation of the neutron flux from the core to the locations of interest on the pressure vessel. These data include material compositions, regional temperatures, and geometry of the pressure vessel, core, and internals. The geometrical input data include the dimensions and locations of the fuel assemblies, reactor internals (baffle, shroud, former plates, core support barrel, thermal shield, and neutron pads), the pressure vessel (including identification and location of all welds and plates) and cladding, and surveillance capsules. For cavity dosimetry, input data should also include the width of the reactor cavity and the material compositions of the support structure and concrete (biological) shielding, including water content, rebar, and steel. The data should, to the extent possible, be based on documented and verified as- built dimensions and plant-specific materials. The isotopic compositions of important constituent nuclides within each region should be based on as-built materials data. In the absence of plant- specific information, nominal compositions and design dimensions may be used; however, in this case conservative estimates of the variations in the compositions and dimensions should be made and accounted for in the determination of the fluence uncertainty (Regulatory Position 1.4.1). The determination of the concentrations of the major isotopes responsible for the fluence attenuation (e.g., iron and water) should be emphasized. The water number densities should be based on plant full-power operating temperatures and pressures, as well as standard steam tables. The data should account for axial and radial variations in water density caused by temperature differences in the core and inside the core barrel, as well as the presence of in-channel and downcomer voids in the case of BWRs. The in-channel and bypass water and fuel channel may be combined to determine a homogenized material region. This approximation has been evaluated in Reference 4 and results in less than a ~3% additional uncertainty in the vessel fluence.
1.1.2 Cross-Sections The calculational method to estimate vessel fluence should use the neutron cross-sections over the energy range from ~0.1 MeV to ~15 MeV and should apply the latest version of the Evaluated Nuclear Data File (currently ENDF/B-VI). These data have been thoroughly reviewed and tested relative to experimental benchmarks.2 Cross-section sets based on earlier or equivalent
2 It should be noted that in many applications the ENDF/B-IV and the first three MODs of the ENDF/B-V iron cross-sections result in as much as ~20% underprediction of the vessel inner-wall fluence and ~35% underprediction of the cavity fluence (Refs.
5-7). Updated ENDF/B-VI iron cross-section data (Ref. 8) have been demonstrated to provide a more accurate determination of the flux attenuation through iron (Refs. 5, 6) and are strongly recommended. These new iron data are included in ENDF/B-VI.
4
nuclear-data sets that have been thoroughly benchmarked for a specific application may be used for that specific application.3 However, when the evaluated cross-section data change, the effect of these changes on the licensee-specific methodology must be evaluated and the fluence estimates updated when the effects are significant (see reporting requirements in Regulatory Position 3).
1.1.2.1 Multigroup Libraries. Transport theory codes that are used to determine the neutron fluence, which employ a multigroup approximation, require cross-section libraries in which the basic data contained within the ENDF files has been preprocessed into a multigroup structure. The development of a multigroup library should consider the adequacy of the group structure, the energy dependence of the flux used to average the cross-sections over the individual groups, and the order of the Legendre expansion of the scattering cross-section. Sufficient details of the energy and angular dependence of the differential cross-sections (e.g., the minima in the iron total cross-section) should be included to preserve the accuracy in attenuation characteristics.
1.1.2.2 Constructing a Multigroup Library. The construction of the multigroup library involves the selection of a problem-independent, fine-multigroup, master library containing data for all required isotopes. This master library should include a sufficiently large number of groups
(/100) that differences between the shape of the assumed flux spectrum and the true flux have a negligible effect on the multigroup data. This library typically includes ~50-100 energy groups above ~0.1 MeV. In addition, a minimum of a P-3 Legendre expansion of the scattering cross- section must be used for typical LWR configurations (Ref. 10). Several libraries satisfying these provisions are available from RSICC, the Oak Ridge National Laboratory Radiation Safety Information Computational Center (Refs. 11-13).
The number of groups may be reduced, with little loss in accuracy, by collapsing the data in the master library over spectra that more closely approximate the true spectra. The vessel fluence calculations should be performed with a job library that has been determined by collapsing the master library. This reduction may be accomplished with a one-dimensional calculation that includes the discrete regions of the core, vessel internals, by-pass and downcomer water, pressure vessel, reactor cavity, shield, and support structures. The resulting job library should consist of multigroup cross-sections based on the region-specific isotopic compositions. This job library should include ~20 energy groups above ~0.1 MeV. The adequacy of the job library must be demonstrated by comparing calculations for a representative configuration performed with both the master and job libraries. In these comparisons, the threshold-detector cross-sections and reaction rates obtained with the fine-group structure should be preserved in the multigroup calculations. In addition to the validation of the collapsing of the master library, additional qualification of the job library is provided by the benchmarking comparisons of Regulatory Position 1.4.2.
There are several ~50-group libraries available from RSICC that were generated using light-water reactor (LWR) spectra for the group collapsing from master libraries (Refs. 14-16) and may be used for LWR application. These libraries contain microscopic cross-sections as well as some premixed macroscopic cross-sections for relevant mixtures. Because these prepackaged
3 The ASTM E1018 standard (Ref. 9) provides additional guidance on the selection of dosimeter reaction rate cross-sections.
5
libraries were designed specifically for LWR pressure vessel fluence calculations, their applicability in any atypical application (for example, when an accurate thermal flux is required)
should be verified prior to use.
1.2 Core Neutron Source The determination of the neutron source for the pressure vessel fluence calculations should include the temporal, spatial, and energy dependence together with the absolute source normalization.
The spatial dependence of the source should be based on depletion calculations that simulate core operation or measured data. The depletion calculations can be performed in either two or three dimensions. Three-dimensional calculations will provide the source in both the radial and axial directions. If two-dimensional calculations are used, the axial effects should be based on measured data.
The core neutron source should be determined by the power distribution (which varies significantly with fuel burnup), the power level, and the fuel management scheme. The detailed state-point dependence should be accounted for (Refs. 17, 18); however, if this is not feasible, averaging over the operating power distribution is acceptable and may be obtained by either (1)
averaging representative power distributions within the cycle or (2) assuming the cycle-average assembly power distribution is well approximated by the accumulated-exposure distribution at the end of the cycle.
A best-estimate power distribution may be used for reactor vessel neutron fluence calculations. However, this best-estimate power distribution must be updated if changes in core loadings, surveillance measurements, or other information indicate a significant change in projected fluence values. This updating may be avoided by using a conservative generic average power distribution, provided no measured distribution yields higher power levels for the important peripheral assemblies.
The peripheral assemblies, which contribute the most to the vessel fluence, have strong radial power gradients, and these gradients should not be neglected (Ref. 19). In the case of PWRs, the fuel pins closest to the core periphery tend to have reduced relative power while for BWRs the peripheral fuel pins may have increased relative power. The pin-wise source distribution should be used for best-estimate calculations and the peripheral-assembly pin-wise source data should be obtained from core depletion calculations. The pin-wise source distribution should represent the absolute source distribution in the assembly.
The local source should be determined as the product of the fission rate and the neutron yield. The energy dependence of the source (i.e., the spectrum) and the normalization of the source to the number of neutrons per megawatt must account for the fact that changes in the isotopic fission fractions with fuel exposure (caused by Pu build-up) result in variations in the fission spectra, the number of neutrons produced per fission, and the energy released per fission.
These effects tend to increase the fast neutron source per megawatt of power for high-burnup assemblies. The variations in these physics parameters with fuel exposure may be obtained from
6
standard lattice physics depletion calculations (Refs. 20, 21). This effect is particularly important for plants that have adopted the PWR low-leakage refueling schemes (vs. out-in-in three-batch fuel management) in which once-, twice-, or thrice-burned fuel is located in the high-importance peripheral locations (Refs. 22-25). The harder spectrum in the BWR fuel regions having a high void fraction will have a similar effect on the isotopic fission fractions and on the neutron-source normalization and spectrum.
A planar-octant representation is generally acceptable for the octant-symmetric fuel- loading patterns typically employed in LWRs. However, in the case of BWRs in which the jet- pump positioning is quadrant-symmetric, a quadrant-symmetric model may be required. For determining the peak fluence, fuel-loading patterns that are not octant symmetric (e.g., as in certain low-leakage patterns) may be represented in octant geometry using the octant having the highest fluence. For evaluating dosimetry, the octant in which the dosimetry is located may be used.
When the actual planar core rectangular geometry cannot be modeled (e.g., in the case of (r,) discrete ordinates calculations), the core geometry may be described using an (r,) planar representation. To accurately represent the important peripheral assembly geometry, a -mesh of at least 40 angular intervals over an octant must be applied. The (r,) representation should reproduce the true physical assembly area to within ~0.5% and the pin-wise source gradients to within ~10%. The assignment of the (x,y) pin-wise powers to the individual (r,) mesh intervals should be made on a fractional area or equivalent basis (Ref. 26). Reference 26 is particularly useful if the radial mesh is a function of (). The overall source normalization should be performed with respect to the (r,) source so that differences between the core area in the (r,)
representation and the true core area do not bias the fluence predictions. A discussion of the uncertainties in determining the core neutron source is provided in Reference 27.
Determination of the three-dimensional (3-D) fluence at the vessel using (r,z)- and (r)-
geometry calculations may also be appropriate (see Regulatory Position 1.3.4). If these calculations are used to provide an axial correction factor, since the (r,z)- and the (r)- dependent fluences enter as a ratio (see Equation 4), the source specification may be less stringent if consistent sources are used in both the (r,z) and (r) calculations.
1.3 Fluence Calculation
1.3.1 Discrete Ordinates Transport Calculation The transport of neutrons from the core to locations of interest in the pressure vessel may be determined with a two-dimensional discrete ordinates transport program (Refs. 28-30) in (r,)
and, when appropriate, in (r,z)- and (r)- geometries.4,5 When calculating a horizontal plane of the core/vessel geometry in which the rectangular (x,y) geometry of the core boundary and the
4 Additional information concerning the application of transport methods to reactor vessel surveillance is provided in ASTM
Standard E 482-89 (Ref. 31)
5 If DOT (Ref. 29) is used, the "-weighted" option (MODE-5 in DOT 4.3) is considered to be more accurate than the "weighted"
option (MODE-3 in DOT 3.5 or 4.3) for flux extrapolation and is recommended. Also, a value of = 0.3 is generally adequate (Refs. 4, 32, 33). Reference 34 includes a description of the directional -weighted scheme which may also be used.
7
cylindrical (r) geometry of the vessel are mixed, a more accurate description is provided by the variable (r,) mesh option (Ref. 29) and may be applied.
An azimuthal () mesh using at least 40 intervals over an octant in (r,) geometry in the horizontal plane should provide an accurate representation of the spatial distribution of the material compositions and source described in Regulatory Position 1.2. The radial mesh in the core region should be ~2 intervals per inch for peripheral assemblies and may be much more coarse for assemblies more than approximately two assembly pitches removed from the core- reflector interface. In excore regions, a spatial mesh that ensures the flux in any group changes by less than a factor of ~2 between adjacent intervals should be applied, and a radial mesh of at least
~3 intervals per inch in water and ~1.5 intervals per inch in steel should be used. Because of the relatively weak axial variation of the fluence, a coarse axial mesh of ~0.5 interval per inch may be used except near material and source interfaces, where flux gradients can be large. An S8 fully symmetric angular quadrature must be used as a minimum for determining the fluence at the vessel. However, in reactor cavity fluence calculations a higher order quadrature may be needed, depending on the width of the cavity and the axial location at which the fluence is being calculated.
Where computer-storage limitations prevent the implementation of these mesh densities in a single-model representation, the calculation should be performed in two or more "bootstrap"
steps rather than compromising the spatial mesh or quadrature (the number of groups used usually does not affect the storage limitations, only the execution time). In this approach, the problem volume is divided into overlapping regions. In a two-step bootstrap calculation, for example, a transport calculation is performed for the cylinder defined by 0 < r < RN with a fictitious vacuum- boundary condition applied at RN. From this initial calculation a boundary source is determined at the radius RO = RN - and is subsequently applied as the internal-boundary condition for a second transport calculation from RO to R (the true outer boundary of the problem). The adequacy of the overlap region must be tested (e.g., by decreasing the inner radius of the outer region) to ensure that the use of the fictitious boundary condition at RN has not unduly affected the boundary source at RO or the results at the vessel.
A point-wise flux convergence criterion of .0.001 should be used, and a sufficient number of iterations should be allowed within each group to ensure convergence. To avoid negative fluxes and improve convergence, a weighted-difference model should be used.5 The adequacy of the spatial mesh and angular quadrature, as well as the convergence criterion, must be demonstrated by tightening the numerics until the resulting changes are negligible (Ref. 32). In discrete ordinates codes, the spatial mesh and the angular quadrature should be refined simultaneously. In some cases, these evaluations can be adequately performed with a one- dimensional model.
The application of the discrete ordinates methods described in this section is illustrated in detail in Reference 4.
1.3.2 Monte Carlo Transport Calculation The transport of neutrons from the core to locations of interest in the pressure vessel may be determined by using Monte Carlo Transport (Refs. 4 and 35-46). The Monte Carlo method has
8
the advantage of allowing an exact representation of the problem geometry (Refs. 43-44); i.e., the rectangular fuel assemblies, the barrel, thermal shield and vessel cylindrical regions, and the axial geometry including upper and lower core reflectors, axial dosimeter locations, and vessel circumferential weld locations. In addition, the Monte Carlo method allows a continuous (as well as a multigroup) energy description of the nuclear cross-sections and flux solution. It is noteworthy, however, that in order to obtain comparable accuracy this method typically requires substantially longer computing times than the discrete ordinates method.
The cross-section data used in the Monte Carlo calculations should satisfy the requirements of Regulatory Position 1.1.2. However, if a continuous energy cross-section library is used, the requirements concerning the Legendre expansion of the scattering cross-sections and the number of energy groups are not applicable. In addition, to ensure an accurate integration over the dosimeter cross-section, the energy bin structure used to determine dosimeter response scoring should satisfy the requirements of Regulatory Position 1.1.2 concerning the selection of the multigroup library group structure.
While the Monte Carlo programs allow a great deal of flexibility in defining problem geometry, the input is generally relatively complex, requiring the definition of a large number of surfaces and cells, and is a potential source of error. Consequently, the geometry input should be checked thoroughly (using plots, volume checks, etc.) against the model design data to validate the definition of the geometry.
In pressure vessel evaluations, the flux is required at specific dosimeter and vessel surface and internal locations. The calculated flux is edited in area and volume tally regions at these locations. Typically, the size of these tally regions is increased to increase the number of particle histories and reduce the statistical error in the flux estimate. The sensitivity of the flux estimate to the volume or area of the tally region should be determined to ensure that the spatial dependence of the flux does not introduce a bias in the flux estimate. If the size of the tally region introduces a bias into the flux estimate, the calculation should be rerun with a smaller tally region and an increased number of particle histories or the Monte Carlo prediction should be adjusted to eliminate the calculational bias. This adjustment may be made, for example, using a correction based on a deterministic calculation.
In addition to the area and volume tally regions, point detectors may be used to estimate the flux at the locations of interest. However, the use of point detectors is complicated and the following cautions should be observed. To minimize the variance in point detector edits, Monte Carlo codes may define an average-flux region surrounding the detector within which the collision estimation is determined. The neutron flux in this region is assumed to be uniform. In order to ensure that the point estimate is not biased, the flux estimate should be shown to be insensitive to the selection of the average-flux region. If the size of the average-flux region introduces a bias into the flux estimate, the calculation should be rerun with a smaller average-flux region and an increased number of particle histories. No material boundaries should be included in the average- flux region surrounding the point detector. Since the next-event particle tracking used in point detectors does not account for reflecting, periodic, or white boundary conditions, the point detector flux estimate is underpredicted when particles from these boundaries can contribute to the flux estimate. Therefore, point detectors should not be located where particles from these
9
boundaries can contribute to the flux estimate. In order to ensure that a sufficient number of particle histories are included and that the region close to the detector is adequately sampled, the statistical edits should be reviewed thoroughly, especially if the relative error is greater than 5%.
In the Monte Carlo analysis, a calculational uncertainty is introduced as a result of the finite number N of particle histories sampled. The Monte Carlo analysis provides estimates of the mean flux, relative error (standard deviation/mean) and terms related to the higher moments, such as the variance of the variance (VOV), based on the calculated particle histories. In order to ensure that the calculation has converged and the estimated (mean) flux is valid, a reliable estimate of the sampling uncertainty is required. However, because of the finite sampling of the physical phase space, important paths in the geometry and peaks in the cross-sections may not be adequately sampled, resulting in significant errors in the estimated mean and standard deviation.
Monte Carlo programs include statistical criteria to provide assurance that the calculations have converged and that the estimated mean and relative error are valid. When performing Monte Carlo pressure vessel fluence calculations, all statistical tests provided by the Monte Carlo code should be satisfied or justification for accepting the results should be reported. The Monte Carlo code should provide sufficient statistical testing (in addition to the usual tally standard deviation)
to ensure that the scoring phase space has been adequately sampled. Representative statistical tests that are used to ensure tally convergence include (e.g., Reference 43):
(1) The flux mean should not have a significant monotonic dependence (either increasing or decreasing) on N for the last half of the problem,
(2) The percent relative error of the mean flux should be less than ~10%, except for point detectors for which it should be less than ~5%,
(3) The relative error should have a ~ 1/ N dependence for the last half of the problem,
(4) The figure of merit defined as
1 FOM ' (Equation 1)
R 2T
where R is the relative error and T is the problem computing time, should not have a significant monotonic dependence on N for the last half of the problem, and
(5) The VOV should be less than ~0.10 for all tallies and decrease as ~1/N for the last half of the problem.
As a result of the strong flux attenuation between the core and the pressure vessel, very few neutrons (~less than one in a thousand) actually reach the vessel. Since the relative error (and variance) in the Monte Carlo calculation only decreases as 1/ N , the number of neutrons tracked and resulting computer times required to obtain acceptable accuracy are very large. In practice, the computing time may be reduced by using variance reduction techniques that increase the number of neutrons reaching the locations of interest (such as the dosimetry capsules and vessel inner-wall) or improve the particle tracking efficiency, or both. Typically, the number of particles that contribute to the tally is increased by increasing the number of particles in the regions of the particle phase space that are judged (based on prior knowledge) to contribute significantly to the
10
calculation. However, if this judgment is incorrect and the neutrons believed to be unimportant make a significant contribution to the tally, the Monte Carlo estimated fluence and relative error (and variance) will be erroneous. Consequently, the variance reduction methods used in performing pressure vessel fluence analyses should be qualified by comparing the Monte Carlo variance reduction predictions with estimates made without the application of the variance reduction technique. Since the calculation without variance reduction will typically involve a large number of particle histories and long computing times, this comparison can be made using a simpler representative configuration (e.g., one-dimensional geometry).
The following Monte Carlo variance reduction methods may be used in performing pressure vessel fluence analyses: (1) neutron energy cutoff, (2) source biasing, (3) geometry splitting with Russian Roulette, and (4) weight windows. Specific concerns associated with the application of these methods are described in the following.
In the application of the neutron energy cutoff technique, neutron tracking efficiency is improved by terminating tracking when the particle weight falls below a selected input energy.
This method assumes that particles with energy below the cutoff do not contribute to the tally and, consequently, may underestimate or bias the tally edit. While a cutoff energy of 1-MeV can be used to calculate the E > 1-Mev fluence without introducing any error, the cutoff energy for the dosimeter response calculation depends on the specific dosimeter reaction cross-section and should be selected low enough to ensure that the calculation provides acceptable accuracy. If this method is employed, the bias introduced by the energy cutoff should be estimated by comparison with an unbiased calculation and the Monte Carlo calculation adjusted to account for the bias.
The source biasing method increases the number of source particles in the regions of phase space that are believed (based on prior knowledge) to provide the most significant contribution to the Monte Carlo tally. For example, the neutrons born in the peripheral fuel assemblies contribute most to the vessel fluence and the number of these particles can be increased (and their weight decreased proportionally) to improve the particle tracking efficiency. The space- and energy- dependent importance of the core neutrons to a specific location on the vessel can be estimated using a precalculated adjoint flux and used to determine the source biasing (e.g., Refs. 41, 42, and
44). To ensure that the biased sampling of the source is producing reliable Monte Carlo estimates, the source biasing method should be qualified by comparison with an unbiased calculation for a representative problem.
Geometrical splitting with Russian Roulette reduces the variance of the tally edit by increasing the number of neutrons that reach the vessel. This is accomplished by assigning numerical importance values to each spatial cell such that the cell importance values increase as the neutrons approach the vessel (Ref. 46). When a neutron travels toward the vessel passing from a region of lower importance I to a region of higher importance IN, it is split into IN/I
neutrons. Conversely, when a neutron travels away from the vessel passing from a region of higher importance I to a region of lower importance IN, it undergoes Russian Roulette with a survival probability of IN/I. The reliability of this technique depends on the selection of the cell importances I. It is recommended that the importance in adjacent cells should not change by more than a factor of ~3, and the cell radial optical thickness should be less than ~2 mean free paths (Ref. 43). To ensure that the selection of the cell importance distribution provides reliable Monte
11
Carlo estimates, the geometrical splitting with Russian Roulette method should be qualified by comparison with an unbiased calculation for a representative problem.
The weight windows method may also be used to improve the efficiency of the neutron tracking in Monte Carlo pressure vessel analyses (e.g., Refs. 41, 42, and 45). In the weight windows method, the neutron weight in each space-energy cell is limited to a specific range or weight window by splitting and Russian Roulette. The window is defined so that the neutron weight is inversely proportional to the importance of the space-energy cell so that the scoring variance is minimized (Ref. 43). In practice, the specification of the weight windows is difficult and requires a detailed knowledge of the space-energy-dependent cell importance. Consequently, in order to ensure that the definition of the weight windows is providing reliable Monte Carlo estimates, it is important that the weight windows method be qualified by comparison with a calculation performed without variance reduction.
The application of Monte Carlo methods described in this section is illustrated in detail in Reference 4.
1.3.3 Fluence Determination The transport calculations (both discrete ordinate and Monte Carlo) may be performed in either the forward or adjoint modes. When several transport calculations are needed for a specific geometry, assembly importance factors may be precalculated by either performing forward calculations with a unit source (with the desired pin-wise source distribution) specified in the fuel assembly of interest or by performing adjoint calculations. The adjoint fluxes are used to determine the fluence at a specific (field) location, while the forward fluxes from the unit-source calculations determine the fluence at all locations in the problem. Once calculated, these factors contain the required information from the transport solution, and by weighting the assembly importance factors with the source distribution of interest, the vessel (or capsule) fluence may be determined without additional transport calculations, assuming the in-vessel geometry, material, and source distribution within the assembly remain the same.
In performing calculations of surveillance capsule fluence (Regulatory Position 1.4), it should be noted that the capsule fluence is extremely sensitive to the representation of the capsule geometry and internal water region (if present), and the adequacy of the capsule representation and mesh must be demonstrated using sensitivity calculations (as described in Regulatory Position
1.4.1). The capsule fluence and spectra are sensitive to the radial location of the capsule and its proximity to material interfaces (e.g., at the vessel, thermal shield, and concrete shield in the cavity), and these should be represented accurately. The core shroud and baffle former plates can result in a 5-10% reduction in the surveillance capsule dosimeter response and should be included in the model.
When fluence reduction schemes have introduced strong axial or azimuthal heterogeneities into the source (e.g., an axially zoned replacement of fuel with stainless steel for fluence reduction), these should be modeled in detail. When the transport calculation is performed using the discrete ordinates method, a finer spatial mesh and tighter convergence criteria may be appropriate to ensure an accurate solution. These schemes may also entail a 3-D flux calculation.
12
To account for the neutron spectrum dependence of RTNDT when the E > 1 -MeV fluence is extrapolated from the inside of the pressure vessel to the T/4- and 3T/4-vessel locations, a spectral lead factor (which accounts for the change in neutron spectrum between downcomer and vessel internal locations) must be applied to the fluence for the calculation of RTNDT (Ref. 47). This spectral lead factor has been included in the Equation 3 attenuation formula of Revision 2 of Regulatory Guide 1.99, and therefore is not required when this formula is used. However, when this formula is not used, the spectral lead factor must be applied to the fluence at the vessel internal locations. Displacements per atom is used as the extrapolation parameter in Regulatory Guide 1.99 and is an acceptable parameter for extrapolating the E > 1-MeV fluence to vessel internal locations and determining this lead factor.
1.3.4 Synthesis of the 3-D Fluence When 3-D effects are important and three-dimensional transport calculations are not practical, a 3-D fluence representation may be constructed by synthesizing calculations of lower dimensions using the expression Ng(r,2,z) = Ng(r,2)(Lg(r,z) (Equation 2)
where g(r,) is the group-g transport solution in (r,) geometry for a representative plane and Lg(r,z) is a group-dependent axial shape factor. Two simple methods available for determining Lg(r,z) are defined by the expressions Lg(r,z) ' P(z) (Equation 3)
where P(z) is the peripheral-assembly axial power distribution, and L g(r,z) = Ng(r,z)/Ng(r) (Equation 4)
where g(r) and g(r,z) are one- and two-dimensional group-g flux solutions, respectively, for a cylindrical representation of the geometry that preserves the important axial source and attenuation characteristics (Refs. 4 and 48). The (r,z) plane should correspond to the azimuthal location of interest (e.g., peak vessel fluence or dosimetry locations) or a conservatively -
averaged (r,z) plane. The source per unit height for both the (r,)- and (r)- models should be identical, and the true axial source density should be used in the (r,z) model.
Equation 3 is only applicable when (1) the axial source distribution for all important peripheral assemblies is approximately the same or is bounded by a conservative axial power shape and (2) the attenuation characteristics do not vary axially over the region of interest. In addition, since the axial flux distribution tends to flatten as it propagates from the core to the pressure vessel, for typical axial power shapes, use of Equation 3 will tend to overpredict axial flux maxima and underpredict minima. This underprediction is nonconservative and can be large near the top and bottom reflectors, as well as when minima are strongly localized as occurs in some fluence-reduction schemes. Because of the neglect of the former plates, this method may result in an additional 5-10% underprediction of the fluence.
13
Equation 4 is applicable when the axial source distribution and attenuation characteristics vary radially but do not vary significantly in the azimuthal () direction within a given annulus.
For example, this approximation is not appropriate when strong axial fuel-enrichment variations are present only in selected peripheral assemblies.
In summary, an (r,)-geometry fluence calculation and a knowledge of the peripheral assembly axial power distribution are needed when using Equation 3. Use of this equation may result in fluence overpredictions near the midplane at relatively large distances from the core, e.g.,
in the cavity, and underpredictions at axial locations beyond the beltline at relatively large radial distances from the core. Conservatism may be included in the latter case by using the peak axial power for all elevations.
Both radial and axial fluence calculations are needed when using Equation 4; thus, it is generally more accurate in preserving the integral properties of the three-dimensional fluence.
Both Equation 3 and Equation 4 assume separability between the axial and azimuthal fluence calculations, which is only approximately true.
When these simple synthesis techniques are not applicable, multichannel synthesis methods may be used. In the multichannel synthesis calculation, the fluence is represented as the sum N
Ng(r,2,z) = j aiNgi(r,2)Ngi(r,z)/Ngi(r) (Equation 5)
i'1 where the gi are basis flux solutions, typically representing specific regions of the core and vessel geometries, and the weighting coefficients ai are determined to provide an optimum prediction of the vessel fluence. It should be emphasized, however, that the accuracy of this method is sensitive to the selection of the basis functions, especially at region interfaces, and three-dimensional calculations should be considered where strong axial or azimuthal heterogeneities exist. This synthesis technique has been applied to a calculational benchmark in Reference 28 and to an experimental benchmark in Reference 49.
1.3.5 Cavity Fluence Calculations Accurate cavity fluence calculations are relied on to analyze dosimeters located in the reactor cavity (Ref. 50). The calculation of the neutron transport in the cavity is made difficult by (a) the strong attenuation of the E > 1 MeV fluence through the (5.5- to 10.0-inch thick) vessel and the resulting increased sensitivity to the iron inelastic-scattering cross-section and (b) the possibility of neutron streaming (i.e., strong directionally dependent) effects in the low-density materials (air and vessel insulation) in the cavity. Because of the increased sensitivity to the iron cross-sections, ENDF/B-VI cross-section data should be used for cavity fluence calculations.2 As indicated in Regulatory Position 1.1.2, properly benchmarked alternative cross-sections may also be used. However, for cavity applications, the benchmarking must include comparisons for operating reactor cavities or simulated cavity environments.
14
Typically, the width of the cavity together with the close-to-midplane locations of the dosimetry capsules result in minimal cavity streaming effects. Consequently, in discrete ordinates calculations, an S8 angular quadrature is acceptable. However, when off-midplane locations are analyzed, the adequacy of the S8 quadrature to determine the streaming component must be demonstrated with higher-order Sn calculations. In addition, since the radial mesh in the (r,z)
calculation is generally finer than the z-mesh in the cavity resulting in narrow spatial-mesh intervals, a -weighted difference model should be used.5 The cavity fluence is sensitive to both the material compositions (Regulatory Position
1.1.1) and the local geometry (e.g., the presence of detector wells) of the concrete shield, and these should be represented as accurately as possible. Benchmark measurements involving simulated reactor cavities that are recommended for methods evaluation are given in Regulatory Position 1.4.2.2. The measured energy spectrum for a typical PWR cavity is described in Reference 51. When both in-vessel and cavity dosimetry measurements are available, an additional verification of the measurements and calculations may be made by comparing the vessel inner-wall fluence determined from (1) the absolute fluence calculation, (2) the extrapolation of the in-vessel measurements, and (3) the extrapolation of the cavity measurements.
Measurements performed in reactor cavities and the associated calculations are described in References 5 and 6.
1.4 Methodology Qualification and Uncertainty Estimates The neutron transport calculational methodology must be qualified, and flux uncertainty estimates must be determined. The neutron flux undergoes several decades of attenuation before reaching the vessel, and the calculation is sensitive to the material and geometrical representation of the core and vessel internals, the neutron source, the nuclear cross-section data, and the numerical schemes used in its determination. The uncertainty estimates are used to determine the appropriate uncertainty allowance to be included in the application of the fluence estimate. While adherence to the guidelines described here will generally result in accurate fluence estimates, the overall methodology must be qualified in order to quantify uncertainties, identify any potential biases in the calculations, and provide confidence in the fluence calculations. In addition, while the methodology, including computer codes and data libraries used in the calculations, may have been found to be acceptable in previous applications, the qualification ensures that the licensee's implementation of the methodology is valid.
The methods qualification consists of three parts: (1) the analytic uncertainty analysis (Regulatory Position 1.4.1), (2) the comparison with benchmarks and operating reactor measurements (Regulatory Position 1.4.2), and (3) the estimate of uncertainty in the calculated fluence (Regulatory Position 1.4.3). These three steps in the overall qualification procedure are discussed below and are outlined in Figure 3.
1.4.1 Analytic Uncertainty Analysis To demonstrate the accuracy of the methodology, an analytic uncertainty analysis must be performed. This analysis includes identification of the important sources of uncertainty. For typical fluence calculations, these sources include:
15
C Nuclear data (transport and dosimeter reaction cross-sections and fission spectra),
C Geometry (location of internals and deviations from the nominal dimensions),
C Isotopic composition of material (density and composition of coolant water, core, core barrel, former plates, thermal shield, pressure vessel with cladding, and concrete shield),
C Neutron sources (space and energy distribution, burnup dependence),
C Methods error (mesh density, angular expansion, convergence criteria, macroscopic group cross-sections, fluence perturbation by surveillance capsules, spatial synthesis, and cavity streaming).
This list of uncertainty components is not necessarily exhaustive, and other uncertainties that are specific to a particular reactor or a particular calculational method should be considered.
In typical applications, the fluence calculational uncertainty is dominated by a few easily identified uncertainty components, such as the geometry and core neutron source, which results in a substantial simplification of the uncertainty analysis.
The sensitivity of the flux to the significant component uncertainties should be determined by a series of transport sensitivity calculations in which the calculational model input data and modeling assumptions are varied and the effect on the calculated flux is determined. (A typical sensitivity would be -10-15% decrease in vessel E > 1 MeV fluence per centimeter increase in vessel inside radius.) Estimates of the expected uncertainties in these input parameters must be made and combined with the corresponding fluence sensitivities to determine the total analytically determined calculation uncertainty ca (i.e., standard deviation) as indicated in Figure 3. The uncertainties should be combined in a statistical (root-sum-of-the-squares) fashion to determine the total fluence uncertainty, accounting for the correlation between uncertainties when necessary (see, e.g., Section 5.2 of Ref. 3). The known systematic errors (or biases) should be combined algebraically, recognizing the sign of each contribution, to determine the overall calculational bias Bca (Ref. 52).6 The component uncertainties should be based on measurement or on the acceptable deviations included in the design specifications. The sensitivity calculations may be performed in one dimension when the model sensitivity does not require a two-dimensional representation.
If the sensitivities are determined using Monte Carlo transport calculations, the magnitude of the model variations and the number of particle histories may have to be increased to provide reliable estimates of the calculation sensitivities. If the increased model variations are large, several calculations may be required to determine the non-linear dependence of the sensitivity.
Perturbation calculations using correlated sampling may also be used to determine the sensitivity to model geometry, material compositions, and cross-sections. Adjoint Monte Carlo calculations may be used to determine the sensitivity to input such as fission spectra. Alternatively, if the sensitivities can be determined with sufficient accuracy in one-dimensional geometry, the sensitivities can be determined using a one-dimensional model.
6 The systematic errors affect the fluence prediction in a specific direction (i.e., increase or decrease the estimate). Consequently, when these errors (or biases) are applied, the signs (or directions) must be recognized and they must be combined algebraically.
The independent random uncertainties have no specific direction (i.e., an increase or decrease) associated with them, and they are combined without recognizing the sign of the error.
16
A sensitivity analysis, in which the influence of each of these uncertainties on the calculated group fluences has been considered, is included in References 5 and 6. Since the uncertainties used in these analyses are common to many pressurized water reactors, the sensitivities (including correlations) may be used to determine initial uncertainty estimates. These variance estimates should be modified as the methods and cross-section data change or if the reactor of interest differs substantially from the reactors analyzed in these references. The referenced sensitivity analysis provides guidance for such modifications.
1.4.2 Comparisons with Benchmark Measurements and Calculations Calculational methods must be validated by comparison with measurement and calculational benchmarks. The fluence calculation methods must be validated against (1)
operating reactor measurements that provide in-vessel surveillance capsule dosimetry or ex-vessel cavity measurements or both, (2) a pressure vessel simulator benchmark that provides measurements at the inner surface and at the T/4 and 3T/4 positions within the vessel (see Regulatory Position 1.4.2.2 for a discussion of such benchmarks), and (3) the fluence calculation benchmark. The results of the validation should include comparisons of reaction rates, fluences, and group fluxes for the locations of interest (Refs. 53, 54).
The methods used to determine the plant-specific data and calculate the benchmarks must be consistent (to the extent possible) with those used to calculate the vessel fluence. That is, the same cross-sections, transport techniques, and transport code parameters that are to be used in the reactor licensing application must be employed in the calculation of the benchmark measurements and calculations.
Differences between measurements and calculations should be consistent with the combined uncertainty estimates for the measurements and calculations. (Note that the uncertainties in both the calculations and measurements will contribute to the observed measurement-to-calculation differences.) The calculated reaction rates (using the methods described in Regulatory Positions 1.1 through 1.3) typically agree with the measurements to within about 20% for in-vessel surveillance capsules and 30% for cavity dosimetry. Deviations greater than these values must be investigated and, when the cause of the deviation is determined to be an error in the calculation, the calculations must be modified.
If the calculation is modified using a least-squares statistical approach (Refs. 55-62), the measurement-to-calculation (M/C) ratios should be reviewed, prior to adjustment, in order to identify any systematic trends in the M/C data. Specific examples include systematic overprediction (or underprediction) of the dosimeter response for (1) dosimeters with low-energy thresholds, (2) dosimeters in the cavity (versus inner-wall and accelerated capsule locations), (3)
dosimeters with long half-life reaction products, and (4) dosimeters monitoring specific cycles. If systematic trends, interdependence, or correlations appear in the M/C data, these should be explained and removed (prior to adjustment) or accounted for in the statistical adjustment
17
procedure.7 Also, since the fluence adjustment is sensitive to the input uncertainties, realistic uncertainty values should be used and documented. The dosimeter measurement uncertainty should include the response uncertainty caused by dosimeter mislocation. Additional guidance on least-squares adjustment methods is provided in Ref. 63. It should be noted, however, that the available least-squares statistical adjustment codes differ substantially in their approach and their capabilities. Consequently, the selection of the least-squares code and the details of the analysis are highly application-specific and are outside the scope of this guide.
The simulator benchmarks provide accurate measurements but typically do not provide an accurate representation of the actual plant configuration. The operating reactor measurements, on the other hand, represent the actual as-built plant configuration but, typically, include substantial measurement uncertainties. The simulator benchmarks together with the operating reactor measurements generally provide an acceptable measurement data base. The comparisons of the calculations to the measurement benchmarks should be used to estimate the calculational bias Bcb and uncertainty cb. If the bias estimate is considered reliable, a correction equal to the bias may be applied to the calculation to determine a best-estimate fluence.8 If the calculations or measurements are adjusted to improve the M/C agreement, using either a least-squares approach or a direct multiplicative bias, the adjustment and basis must be reported.
1.4.2.1 Operating Reactor Measurements. Well documented fluence dosimetry measurements for operating power reactors may be used for methods and data qualification.9 Descriptions of these configurations include three-dimensional geometry, reactor operating conditions and, in some cases, both in-vessel and ex-vessel measurement data.
Comparisons of measurements and calculations must be performed for the specific reactor being analyzed or for reactors of similar design. Plant-specific measurements have the advantage of including the as-built materials and geometry and the actual reactor operating conditions. An especially accurate determination of the fluence can be obtained when both in-vessel and cavity dosimetry are available.
These measurements should not be used to bias or adjust the fluence calculations unless a statistically significant number of measurements is available, the various dosimeter measurements are self-consistent, and a reliable estimate of the calculational bias can be determined.8,9 Similarly, plant-specific biases should not be used unless sufficient reliable measurement data are available. As capsule and cavity measurements become available, they should be incorporated
7 For example, if the M/C data for the Np-237 dosimeters tends to be high, indicating a systematically low calculation or high measurement, the errors in the Np-237 calculated or measured responses are correlated. If this error is not removed prior to the adjustment, the important sources of correlated error should be identified (e.g., the Np-237 dosimeter reaction cross-section, reaction product half-life and Np-237 photo-fission correction) and included in the adjustment procedure to allow the proper treatment of this error.
8 The adequacy of the measurement data base for determining a bias in the calculations depends on the magnitude of the bias and is therefore problem dependent. For example, if the bias is small an accurate estimate of the bias will require either (1) a large number of reasonably accurate measurements or (2) a small number of very accurate measurements. Specifically, the uncertainty in the bias should be substantially less than the bias itself. For example, if the calculated fluence is to be increased by 10% , the uncertainty in this increase should be substantially less than 10%.
9 Additional guidance on benchmark testing is provided in the ASTM standard E2005-99 (Ref. 64) and examples of operating reactor benchmark measurements are provided in Reference 65.
18
into the operating reactor measurements data base, and the calculational biases and uncertainties should be updated as necessary.
1.4.2.2 Pressure Vessel Simulator Measurements. Several pressure vessel simulator benchmarks are available (Refs. 7, 49, 65-82) and may be used for methods qualification. These benchmark experiments were carried out by several laboratories, and dosimetry measurements using different techniques were compared to provide experimental results with well known and documented uncertainties. Three configurations were used in these experiments.
C Pressure Vessel Simulator Experiments. The Oak Ridge National Laboratory (ORNL)
Pool Critical Assembly (PCA) is a full-scale-section mockup of a pressure vessel with dosimetry measurements at the inner surface of the vessel and at the T/4, T/2 and 3T/4 locations within the vessel wall (Refs. 72-75). The PCA benchmark experiment includes a simulated surveillance capsule that allows the determination of the effect of the capsule on the dosimetry measurements. This benchmark is characterized by a relatively simple geometry with generally less uncertainty in region compositions, temperatures, and source distributions than in operating power reactor benchmarks.
C Experiments Simulating Operating Reactor Conditions. The CEN/SK Laboratory VENUS benchmark experiment (Refs. 49, 76-79) and the ORNL Pool Side Facility (PSF)
benchmark experiment (Refs. 80-82) simulate the operating reactor: (1) pin-wise power distributions in peripheral fuel assemblies, (2) three-dimensional effects caused by partial- length shield assemblies, and (3) heterogeneities caused by neutron pads attached to the core barrel.
C Pressure Vessel Cavity Measurements. The H.B. Robinson-2 cavity measurement (Refs. 65-67) and the NESDIP cavity simulator experiment (Refs. 7, 68-71) may be used as benchmarks to qualify cavity fluence calculational methods. Calculations of dosimetry measurements performed in an operating reactor cavity are described in References 5 and 6.
1.4.2.3 Calculational Benchmark. The vessel fluence benchmark problems provided by the NRC in NUREG/CR-6115 (Ref. 4) should be used for methods qualification. This report provides the detailed information defining the problems and corresponding numerical solutions for a set of PWR and BWR pressure vessel fluence benchmark problems. The geometry, materials, and space- and energy-dependent source are fixed by the problem specification. PWR
benchmark problems have been specified for (1) a standard core loading pattern, (2) a low-leakage core loading pattern, and (3) a partial-length shield assembly core design. Since BWR fuel loading patterns are presently not being designed for vessel fluence reduction, only a single BWR
problem is specified. In addition, MCNP Monte Carlo calculations have been performed for both the BWR problem and the PWR problems for the standard core and partial-length shield assembly core loadings. Comparisons of the vessel fluence determined by using both the MCNP and DORT computer codes are included.
The calculation of the benchmark problems allows a detailed assessment and verification of the numerical procedures, code implementation, and the various modeling approximations
19
relative to state-of-the-art solutions for representative operating configurations.10 If the differences between the benchmark problem calculation and the reference solution are substantially larger than what would be expected based on the differences in the methods approximations and nuclear data used in the two calculations, the agreement is considered unacceptable. In this case, the calculation should be reviewed and the differences between the two solutions explained. (Note that NUREG/CR-6115 includes sensitivity calculations and comparisons with Monte Carlo predictions which provide an indication of the accuracy of the methods approximations used in the NUREG/CR-6115 reference calculations.) When the cause of the deviation is determined to be an error in the calculation, the calculational method must be revised.
1.4.3 Estimate of Fluence Calculational Bias and Uncertainty The overall fluence calculation bias Bc and uncertainty c must be determined by an appropriate combination of (1) the analytic uncertainty analysis results of Regulatory Position
1.4.1 (Bca and ca) and (2) the results of the uncertainty analysis based on the comparisons to the operating reactor and simulator benchmark measurements of Regulatory Position 1.4.2 (Bcb and cb).
This combination may be a weighted average that accounts for the reliability of the individual estimates.11 The bias Bc may be applied as a multiplicative correction to the calculated fluence to determine the best-estimate value, and c should be used (when required) as the (1-sigma)
uncertainty in the calculated best-estimate fluence.
The fluence accuracy requirements are generally application specific; however, a vessel fluence uncertainty of 20% (1 sigma) is acceptable for RTPTS and RTNDT determination. In the determination of RTNDT and RTPTS, if the overall calculation uncertainty c is greater than 20%, the calculational model must be adjusted or a correction must be applied to reduce the difference between the fluence prediction and the upper 1-sigma limit to within 20%. For example, if c is greater than 20% but less than 30%, the fluence should be determined from the calculated fluence c by the relation N 6 1 B c [Fc(%) 20]/100 > Nc (Equation 6)
If c is greater than 30%, the methodology of this regulatory guide is not applicable and the application will be reviewed on an individual basis. The procedure used to determine the fluence is illustrated in Figure 4.
10
The calculation being used as the benchmark must be the NUREG/CR-6115 original referenced benchmark calculation, and not just a second independent calculation of the benchmark.
11 Because the weighting of the analytic and benchmark uncertainty estimates depends on the details of the specific application that can vary widely, it is not possible to specify a practical and generically valid prescription for determining the weights.
However, the following example illustrates factors that should be considered. When as-built measurements of the vessel diameter are available and reasonable estimates of the core neutron source and other input uncertainties can be determined, the analytic uncertainty estimate should be reliable and have an uncertainty of -15% (1 sigma). Assume that there are a statistically significant number of accurate ( <~ 5%) operating reactor measurements and the uncertainty estimate based on this data has an uncertainty of -20%. The uncertainty estimate based on the vessel simulator measurements is assumed to be less certain and has an uncertainty of -25%. Using a weighted mean in which the weight is inversely proportional to the square of the standard deviation of the estimate (i.e., -2), the weights are wA = 0.5, wo = 0.3 and wS = 0.2.
20
The fluence calculation methods of this section are summarized in Table 1, Summary of Regulatory Positions on Calculation and Dosimetry.
2. NEUTRON FLUENCE MEASUREMENT METHODS
Dosimetry measurements provide independent estimates of specific activities and isotopic production rates that are used for validating the neutron transport calculations. Fluence is obtained from the response of passive integral detectors placed in surveillance capsules and, more recently, in the ex-vessel cavity. Procedures for performing the measurements and assessing uncertainty, including the validation of measurement methods and detector response in standard neutron fields, is described in this section.
The fluence measurement methods of this section are summarized in Table 1, Summary of Regulatory Positions on Calculation and Dosimetry. The procedure used to qualify the measurement methods is presented in Figure 5. The guidance provided in this Regulatory Position is intended for new dosimetry and, to the extent practical, for existing dosimetry.
2.1 Measurement Procedures Power reactor dosimetry measurement methods use passive integral detectors, which are typically activation detectors and solid-state track recorders. Most of the detectors respond to neutrons with energies above a characteristic reaction threshold. These detectors should be selected with substantial nonoverlapping energy regions (i.e., with well-separated thresholds) to provide coarse spectrum information as well as an estimate of the neutron fluence.
2.1.1 Specification and Application of Dosimeters Neutron dosimetry for pressure vessel surveillance may consist of as-built packages of threshold dosimeters placed in surveillance capsules during reactor construction. The selected dosimeter set should provide adequate spectrum coverage. A common set of fast neutron integral detectors that may be employed in these packages is listed in Table 2. These neutron detectors are discussed in References 83 through 89. The 54Fe(n,p)54Mn and 58Ni(n,p)58Co activation reactions are useful for monitoring neutrons with energies above ~2 MeV. The 46Ti(n,p)46Sc reaction is useful for monitoring neutrons with energies above approximately 4 MeV, and the 63Cu(n,)60Co reaction is useful for measuring neutrons with energies above approximately 5 MeV. Fission monitors are important since they typically have reaction thresholds below 1 MeV. Activities determined by assaying one or more fission products from the 238U(n,f)FP fast-neutron fission reaction are useful for monitoring neutrons with energies from approximately 1.5 to 7 MeV.
Activities determined by assaying fission products from the 237Np(n,f)FP reaction are useful for measuring neutrons with energies from approximately 0.7 to 6 MeV; however, there is a significant gamma-ray background associated with neptunium dosimetry that is due to the large quantity of naturally occurring 233Pa existing in secular equilibrium with 237Np. Compton scattering of the protactinium decay photons generates a large background continuum, which reduces the effective signal-to-noise ratio and lowers the accuracy of the radioassay. The
93 Nb(n,n)93mNb activation reaction is useful for monitoring neutrons with energies above approximately 1 MeV, but involves the counting of a very soft photon, which can complicate the
21
corrections made to the measured activity to compensate for photon attenuation within the dosimeter.
In addition to the fast-neutron threshold monitors, a thermal monitor such as 59Co should be included among the selected dosimeters to allow the determination of the thermal neutron fluence. The thermal fluence is needed to assess the effects of interfering activities in threshold monitors produced by low-energy activation of impurity elements. In instances where the thermal flux is relatively high, it can also be used for assessing the effects of burnout of the activated species.
Taken together with a low-energy detector such as cobalt (to estimate the thermal neutron fluence for determining interference from low-energy activations), the Table 2 dosimeter set provides satisfactory neutron energy spectrum coverage for pressure vessel dosimetry. Alternative detector sets that are used should provide equivalent spectrum coverage. Detector selection criteria and related recommendations in ASTM E 844 (Ref. 90) and E 1005 (Ref. 91) should be followed (see Table 2).
The application of the dosimeter should be evaluated relative to the expected neutron flux and the total irradiation time in order to assess the effects of dosimeter burnout, and the burnout and burn-in of reaction products. For example, when using 58Ni monitors in the presence of a high thermal flux, the burnout of the nickel dosimeter as well as the 58Co and 58mCo reaction products should be considered. The resonant absorption of epithermal neutrons in 238U and the breeding of 239Pu, which can undergo thermal fission and contribute to the fission product activity, should also be considered.
Application of activation detectors involves aspects of the measurement process that must be carefully controlled and documented to obtain accurate results and establish reasonable uncertainty estimates. Where applicable, procedures in ASTM Standards E 181 (Ref. 92), E 844 (Ref. 90) and E 1005 (Ref. 91) and methods devoted to individual radiometric sensors must be used as indicated in Table 2. Specific regulatory positions associated with the dosimetry measurements are indicated in the following.
2.1.1.1 Dosimeter Nuclear and Material Proprieties. The dosimeter nuclear and material properties are important considerations in the selection and composition of the surveillance dosimetry package. The physiochemical properties of the dosimeter materials should be compatible with the prevailing service conditions. Dosimeter melting, for example, can result in the loss of monitor activity and unreliable measurements. Consequently, the melting point of the dosimeters and any associated materials should be carefully considered relative to the service environment ambient temperature and the radiation heating within the dosimeter. Although the melting point of typical radiometric monitors is sufficiently high to preclude melting, even mild heating can be problematical for certain dosimeters. For instance, diffusion losses of 54Mn activity from iron dosimeters can occur when heated to temperatures greater than 700-C. In addition, the dosimeter should be chemically stable and corrosion resistant. Fissionable materials are especially prone to oxidation and, consequently, fission monitors should be hermetically sealed to prevent oxidation and subsequent loss of material.
22
The nuclear properties to be considered in dosimeter selection include the reaction cross- section, decay half-life, gamma-ray yield, and the fission yield when employing fission dosimeters. The reaction cross-section should be of sufficient magnitude to allow the dosimeter to be used as a fluence monitor, consistent with other requirements on the dosimeter mass and geometry, as discussed in following sections.
The photon yield associated with the nuclear decay should be large enough to allow an accurate radio assay, and the energy of the particular photons counted should minimize interference from neighboring spectral lines. For example, the 24Na reaction product, created from the 27Al(n,) reaction, emits photons with energies of 1.369 and 2.754 MeV. Since there is less interference from the contaminant gamma rays associated with the higher-energy photon, it is preferable to count the high-energy line during the radioassay. In the case of fission product radioassay involving the 95Zr-95Nb fission product, three dominant spectral lines are observed at
724, 756, and 765 keV. However, the two high-energy lines are in close proximity and complicate the determination of the background subtraction. Consequently, counting the 724 keV
line provides the more accurate result.
The half-life of the activated species should be long enough, relative to the elapsed time from the end of the irradiation to the time at which the dosimeters are actually counted, to permit nuclear counting to proceed with reasonable counting statistics. This time interval is often a significant consideration. The half-life should also be long enough to provide an accurate indication of the fluence during the irradiation period.
The fission product yields associated with fission dosimeters should be of sufficient magnitude to permit an accurate radioassay. The yields (and decay constants) of the fission products 99Mo and 140Ba-140La are known with greater accuracy than most fission products;
however, their half-life is generally too short for these isotopes to be useful as neutron fluence monitors in power reactor applications. Useful fission products include 95Zr-95Nb, 103Ru, 144Ce, and especially 137Cs-137mBa, because of its 30-year half-life.
2.1.1.2 Material Composition. The dosimeter materials should be pure enough to ensure there is no significant error in the response of the dosimeter from extraneous activities.
Cursory specifications of materials regarding impurities are often unreliable. Specifically, fissile residuals in 237Np and 238U and minute amounts of cobalt in copper and nickel dosimeters (fractional parts per million) should be determined by mass spectrography or radioactivation analysis. Tracking the principal decay half-life is also a useful technique for determining the presence of extraneous activities from impurity elements.
Dosimetric materials should be handled with care to prevent cross contamination of the dosimeters and, when necessary, the dosimeters should be cleaned of surface contaminants. For example, the handling of aluminum monitors can introduce sodium (24Na) as an external surface contaminant. Neutron activation of the sodium will produce 24Na activity that will interfere with the radioassay of the 24Na activity produced from the desired dosimetric reaction.
Post-irradiation quality assurance may be accomplished at the National Institute for Standards and Technology (NIST) for dosimeters already in surveillance capsules by additional
23
irradiation in a thermal neutron field or a standard fission spectrum. These additional exposures and analyses can provide data on the dosimeter masses and impurities (for example, the original cobalt-59 impurity content in copper used for the copper-63 (n, ) reaction).
2.1.1.3 Encapsulation. The detector capsule design must take into account possible activation interference and neutron spectrum perturbation. Thermal neutron shields that eliminate interference from thermal neutron reactions in some detectors must be designed to accommodate radiation heating and should be placed apart from low-energy detectors (see ASTM Standard E
844, Ref. 90). The thickness of the thermal shield should be selected based on the expected thermal flux and duty cycle. (Note that the effective cadmium cutoff energy is a function of shield thickness.) The design of the shields and their placement relative to other dosimeters should be carefully considered in order to minimize their impact on neighboring foils. For example, the use of a cadmium shield will generate a local gamma field that could affect adjacent fission foils by producing a photofission component in the measured activity.
2.1.1.4. Isotopic Mass. The mass of the dosimeters should be selected to permit the production of sufficient activity to allow an accurate radioassay. However, a dosimeter mass that is too large can lead to high activities that result in excessive dead-time losses during nuclear counting and should be avoided. Metal alloys can be used for diluting materials that have high nuclear cross-sections (e.g., dilution of cobalt in an aluminum alloy). However, stoichiometry and isotopic analysis should be well documented for dosimeters that are not of pure natural elements.
2.1.1.5. Geometry and Location. The geometric configuration of the dosimeters can have a significant effect on the dosimeter response. Most radiometric monitors are in the shape of thin circular activation foils, although other shapes are available. Dosimeters should be small enough to allow placement at the desired location and large enough to provide sufficient activity to allow an accurate measurement. The dosimeter thickness can have a significant effect on the dosimeter self-shielding (of the neutron flux) and the dosimeter response for strongly absorbing foils. In addition, increased dosimeter thickness increases the effects of scattering and photon attenuation within the dosimeter and can lead to large response corrections. When selecting fissionable dosimeters, it should be recognized that reduced foil thickness increases fission product loss from nuclear recoil out of the foil, which reduces the measured dosimeter activity.
The placement of the dosimeters can also have a significant effect on the dosimeter measurement. The surroundings of a dosimeter (e.g., adjacent dosimeters or material interface)
can influence detector response. Strongly and weakly absorbing foils that respond over the same energy range should be sufficiently separated within the dosimetry package to minimize interference effects. The arrangement of the dosimeters within the package should minimize scattering effects (viz., thick dosimeters should not be stacked so that large scattering corrections are required for adjacent foils).
The location of individual dosimeters must be determined accurately and recorded, because fluence gradients in out-of-core positions are generally severe. In the pressure vessel cavity, establishing azimuthal position can be as important as the radial location. Specially designed mounting arrangements that ensure accurate radial, azimuthal, and vertical positions
24
should be used for cavity dosimetry. Vertical wires may be used. Comprehensive and accurate detector location information should be maintained.
2.1.1.6 Solid-State Track Recorders. Solid state track recorders (SSTRs) are integral detectors that employ fission reactions. These sensors directly record the tracks of fission fragments from a thin fissionable deposit (Refs. 93 and 94). The principal advantages of these detectors are wide sensitivity ranges and a permanent measurement record. Because the application of SSTRs employs fissionable deposits in the nanogram to picogram range, details of the measurements should be well documented, and standard neutron field calibration should be performed prior to the application. ASTM Standard E 854 (Ref. 95) provides additional information concerning the use of SSTRs.
2.1.1.7 Helium Accumulation Fluence Monitors. The helium accumulation fluence monitor (HAFM) is another type of stable-product neutron monitor. These monitors rely on the generation of helium gas via (n,) reactions that occur in a variety of materials when exposed to a neutron field. The amount of helium generated in a specific material is a function of the neutron fluence and is determined by vaporizing the sample and measuring the helium content by mass spectrometry. HAFMs are routinely used in conjunction with radiometric monitors and, in fact, several commonly used radiometric monitors also undergo (n,) reactions and can be analyzed for helium accumulation as well as for radiometric activity. ASTM Standard E910 (Ref. 96) provides additional guidance on the use of HAFMs.
2.1.2 Detector Response Measurements A measured response must be provided to allow comparison of the calculations and measurements. Typical measured responses include the specific activity at end-of-irradiation (given in disintegrations per second per nucleus), the measured isotopic production (for example, helium atoms per initial atom of material), and the total reactions (for example, fissions per initial atom of material). When comparing calculations with measurements, corrections should be included for detector response perturbations, interfering reactions, and, when applicable, burnup and photofission. In the calculation of end of irradiation activities (or in the conversion of end of irradiation activities into reaction rates), the power-history and the half-life of the dosimeter activation products must be included (Ref. 91). Photofission corrections can vary considerably (from 2-15%) depending upon the location of the dosimetry and the type of reactor. Fission yields should be those specified in the relevant ASTM Standards, the ENDF library, or the validated job library. In situ neutron field perturbations (e.g., by the surveillance capsule and detector encapsulation) must be accounted for if they are not an integral part of the neutron transport calculation.
When reporting measurements, these corrections should be described and quantified along with any other effects that have a significant impact on the measurements. This is especially important because pressure vessel surveillance dosimetry often involves comparison of measurements carried out by different organizations and over long periods of time.
2.1.3 Measurement Uncertainties Regulatory Position 1.4 states that the calculations must be validated by comparison with measured benchmarks. In order to perform this validation, the uncertainty associated with the
25
measured response must be determined for each dosimeter type. The uncertainty must be included in the documentation of the measured results.
The important uncertainties must be quantified and included in an uncertainty table that summarizes the specific components of uncertainty contributing to each detector response.
Typical sources of uncertainty include (1) background counting, counting efficiency, counting statistics, and fraction of the sample counted, (2) decay constant and fission yield uncertainty, (3)
uncertainty in detector weight, geometry, composition and isotopic purity, (4) power history and time of irradiation, (5) dosimeter location uncertainty, and (6) cross-section uncertainty when the fluence is determined. Additional sources are identified in Regulatory Positions 2.1.1 and 2.1.2 and in the dosimeter-specific ASTM Standards of Table 2. The entries in the table should be identified as standard deviations, upper bounds, or appropriate fractions of the correction. Each entry should be described, along with the method used to combine the entries and determine the total response uncertainty. The evaluation of the measurement uncertainty for each dosimeter is important for ensuring a meaningful comparison with the calculations.
2.2 Validation in Standard and Reference Neutron Fields The dosimeter measurements used to benchmark the calculations must meet the requirements in Appendix B to 10 CFR Part 50 for a quality assurance program (for surveillance measurements). To ensure long-term measurement consistency and confirm measurement uncertainties, dosimetry measurements must be performed every few years12 in well characterized neutron fields. If there are changes in dosimeter evaluation procedures, equipment, or personnel that could have a significant effect on the dosimetry, the validation must be updated. This validation may be performed using the Materials Dosimetry Reference Facility (MDRF) or using reference fission neutron sources (Refs. 97, 98). Neutron field referencing may be used as a detector response calibration (Ref. 99).
The validation is accomplished by exposing each type of detector to a certified neutron fluence in the reference neutron field and by determining the fluence using the measurement method to be validated. A calculated spectrum-averaged cross-section, generally specified along with the certified neutron fluence, must be used to derive the measurement value. If measured and certified neutron fluence agree (to within the combined uncertainty of the measured and certified fluence), the detector measurement method, including the detector cross-section, is validated. If the fluences disagree, this calculation-to-experiment (C/E) ratio represents a bias associated with the detector response measurement or the detector cross-section or both. In this case, the detector measurement methods and input parameters should be re-examined in order to eliminate the bias. If after re-examination the bias is still present, the bias may be used directly as a detector response calibration factor.
The results of the neutron field validation procedure should be reported in terms of C/E
ratios for the individual detectors and should include an uncertainty table. The standard neutron field validation may be used, as appropriate, to simplify the uncertainty table called for in
12 This measurement validation procedure does not require transport calculations.
26
Regulatory Position 2.1.3 by reducing or eliminating many uncertainties in the activity measurement, nuclear decay parameters, and detector cross-sections.
Aside from validating the measurement method, standard and reference neutron fields may be used for quality assurance of critical features of the detector response. Examples are activation interference by impurities, proper determination of fission product activities, and mass assay of fissionable deposits for track recorders.
Interlaboratory comparisons, such as those carried out under the NRC-sponsored LWR
Surveillance Dosimetry Improvement Program (Refs. 7, 49, 72, 80, 81), have been useful in validating the quality assurance programs of various laboratories. These interlaboratory comparisons have been used to identify problems in the measurement procedures and establish consistency between the various participating laboratories.
2.3 Fluence Determination from Detector Measurements As stated in Regulatory Position 1.4, the calculated fast-neutron fluence values must be validated by comparison with measurement benchmarks. These comparisons validate the calculational methodology by determining the bias and uncertainty in the calculated fluences, E >
1 MeV. A fast-neutron fluence, E > 1 MeV, must be obtained for each detector as the quotient of the measured reaction probability and the effective E > 1-MeV spectrum averaged cross-section (Ref. 100) based on a neutron transport calculation. These fluences and a suitably weighted average fluence must be reported together with the associated measurement uncertainty determined using the methods of Regulatory Positions 2.1.3 and 2.2 and accounting for the cross- section uncertainty. An alternative to deriving a neutron fluence from the detector responses is to directly compare specific end-of-irradiation activities, measured reaction probabilities, isotopic production, total reactions, average reaction, or measured reaction rates with results from the neutron transport calculation. The response selected should provide a reliable indication of the uncertainty and bias in the E > 1 MeV fluence calculation.
The measurement biases and uncertainties must be documented for the selected detector response. Based on the measurement-to-calculation (M/C) comparisons, an average M/C ratio should be determined as a suitably weighted average of the individual M/C values. In determining the weighting for this average, as a minimum the measurement and spectral average cross-section uncertainties should be considered. An estimate of the uncertainty in the neutron spectrum may be included as an additional contributor to the dosimeter weighting.
When justification can be provided, an individual detector may be declared suspect and given reduced weight or discarded. However, in this case the justification and procedure used must be documented.
2.4 Ex-vessel Dosimetry As-built in-vessel surveillance capsule dosimetry cannot easily be updated. Furthermore, the dosimetry irradiated with metallurgical specimens is only available at infrequent intervals.
However, additional and upgraded dosimetry is important for understanding and following vessel
27
exposures, especially for low-leakage core modifications. The ex-vessel cavity may be used as an alternative site for installing additional improved dosimetry. Recent pressure vessel benchmark experiments (Refs. 7, 67, 94) have demonstrated that the ex-vessel dosimetry can provide useful exposure information within the pressure vessel wall (Refs. 72, 81). When placed at appropriate circumferential locations, this dosimetry is a good monitor of the effectiveness of low-leakage core strategies.
3. REPORTING
When fluence determinations are required by the regulations, the licensee's documentation describing the determination of pressure vessel fluence must provide a complete description of the methods used to calculate and measure the neutron fluences. In applying the methodology of this guide, the details of the application and the results should be reported as described in this section.
A topical report providing the details of the fluence methodology, as outlined in Regulatory Positions 1 and 2, should be satisfactory for compliance with this guide. Plant-specific measurements and data should be reported in the pressure vessel surveillance reports.
The specific regulatory positions on reporting are given in this section and are summarized in Table 1, Summary of Regulatory Positions on Calculation and Dosimetry.
3.1 Fluence Methods The methods used to calculate the integral and multigroup fluences and fluence rates and associated methods qualification should be reported. The calculational uncertainty analysis, benchmark comparisons, and the determination of the overall fluence uncertainty should be described in detail. If the calculated fluence is adjusted, the fluence adjustment and associated uncertainty together with justification must be reported. A discussion of any deviations from the procedures provided in this regulatory guide should be included. The source of the cross-section data, the numerical methods (e.g., quadrature, mesh, and convergence criteria), and the treatment of special effects (e.g., fuel burnup, axial effects, and pin-power distributions) should be described in detail.
3.2 Multigroup Fluences The calculated absolute multigroup neutron fluences and fluence rates at the peak wall (or other limiting) location, surveillance locations, and T/4 and 3T/4 positions within the pressure vessel should be reported. The multigroup energy boundaries should be included.
3.3 Integral Fluences The calculated E > 1 MeV integral fluences and fluence rates at the vessel inner wall locations (determined as described in Regulatory Position 1.4.3), together with the uncertainties, should be reported.
28
3.4 Comparisons of Calculation and Measurement If the qualification of the calculation methods is performed using fluence comparisons, the measured and calculated E > 1 MeV integral neutron fluences should be reported. If the methods qualification is performed using reaction rate comparisons (or other responses as described in Regulatory Position 2.3), the calculated and measured reaction rates (or reaction probabilities)
should be reported. The spectrum-averaged reaction cross-section used to relate the fluence and reaction rate and the method for its determination should also be provided. In either case, the M/C ratios and the measurement uncertainty should be reported for the average fluence and for each detector at each measurement location. If these measurements are used to adjust the calculation, the adjustment should be described and the M/C ratios before and after adjustment should be reported.
3.5 Specific Activities and Average Reaction Rates The specific activities at the end of irradiation and the measured average reaction rates should be reported, together with the associated uncertainty tables and the power-time history.
For each dosimeter, provide the reaction type, dosimeter material and form (wire, foil, etc.),
weight-percent and isotopic-percent of the target material, fission yields, and half-lives. The corrections for the detector response perturbations, interfering reactions, and photofission should also be described. The results of the measurement validation should also be described. If thermal-neutron fluence rate measurements have been performed, these should be reported together with the uncertainty and the dosimeter thermal-neutron shield material.
D. IMPLEMENTATION
The purpose of this section is to provide information to applicants and licensees regarding the NRC staffs plans for using this regulatory guide.
Except in those cases in which an applicant or licensee proposes an acceptable alternative method for complying with the specified portions of the NRCs regulations, the methods described in this guide will be used in the evaluation of applications for new licenses and for evaluating compliance with 10 CFR 50.61 and Appendix G to 10 CFR Part 50, which is imposed by 10 CFR 50.60.
29
TABLE
1. SUMMARY O
F. REGULATORY POSITION
S ON
CALCULATION AND DOSIMETRY
FLUENCE CALCULATION METHODS Regulatory Position Fluence Determination. Absolute fluence calculations, rather than extrapolated 1.3 fluence measurements, must be used for the fluence determination.
Modeling Data. The calculation modeling (geometry, materials, etc.) should be 1.1.1 based on documented and verified plant-specific data.
Nuclear Data. The latest version of the Evaluated Nuclear Data File (ENDF/B) 1.1.2 should be used for determining nuclear cross-sections. Cross-section sets based on earlier or equivalent nuclear-data sets that have been thoroughly benchmarked are also acceptable. When the recommended cross-section data change, the effect of these changes on the licensee-specific methodology must be evaluated and the fluence estimates updated when the effects are significant.
Cross-Section Angular Representation. In discrete ordinates transport 1.1.2 calculations, a P3 angular decomposition of the scattering cross-sections (at a minimum) must be employed.
Cross-Section Group Collapsing. The adequacy of the collapsed job library must 1.1.2 be demonstrated by comparing calculations for a representative configuration performed with both the master library and the job library.
Neutron Source. The core neutron source should account for local fuel isotopics 1.2 and, where appropriate, the effects of moderator density. The neutron source normalization and energy dependence must account for the fuel exposure dependence of the fission spectra, the number of neutrons produced per fission, and the energy released per fission.
End-of-Life Predictions. Predictions of the vessel end-of-life fluence should be 1.2 made with a best-estimate or conservative generic power distribution. If a best estimate is used, the power distribution must be updated if changes in core loadings, surveillance measurements, or other information indicate a significant change in projected fluence values.
Spatial Representation. Discrete ordinates neutron transport calculations should 1.3.1 incorporate a detailed radial- and azimuthal-spatial mesh of ~2 intervals per inch radially. The discrete ordinates calculations must employ (at a minimum) an S8 quadrature and (at least) 40 intervals per octant.
30
Regulatory Position Multiple Transport Calculations. If the calculation is performed using two or 1.3.1 more "bootstrap" calculations, the adequacy of the overlap regions must be demonstrated.
Point Estimates. If the dimensions of the tally region or the definition of the 1.3.2 average-flux region introduce a bias in the tally edit, the Monte Carlo prediction should be adjusted to eliminate the calculational bias. The average-flux region surrounding the point location should not include material boundaries or be located near reflecting, periodic, or white boundaries.
Statistical Tests. The Monte Carlo estimated mean and relative error should be 1.3.2 tested and satisfy all statistical criteria.
Variance Reduction. All variance reduction methods should be qualified by 1.3.2 comparison with calculations performed without variance reduction.
Capsule Modeling. The capsule fluence is extremely sensitive to the geometrical 1.3.3 representation of the capsule geometry and internal water region, and the adequacy of the capsule representation and mesh must be demonstrated.
Spectral Effects on RTNDT. In order to account for the neutron spectrum 1.3.3 dependence of RTNDT, when it is extrapolated from the inside surface of the pressure vessel to the T/4 and 3T/4 vessel locations using the E > 1-MeV fluence, a spectral lead factor must be applied to the fluence for the calculation of RTNDT.
Cavity Calculations. In discrete ordinates transport calculations, the adequacy of 1.3.5 the S8 angular quadrature used in cavity transport calculations must be demonstrated.
Methods Qualification. The calculational methodology must be qualified by both 1.4.1,
(1) comparisons to measurement and calculational benchmarks and (2) an analytic 1.4.2, uncertainty analysis. The methods used to calculate the benchmarks must be 1.4.3 consistent (to the extent possible) with the methods used to calculate the vessel fluence. The overall calculational bias and uncertainty must be determined by an appropriate combination of the analytic uncertainty analysis and the uncertainty analysis based on the comparisons to the benchmarks.
Fluence Calculational Uncertainty. The vessel fluence (1 sigma) calculational 1, uncertainty must be demonstrated to be #20% for RTPTS and RTNDT determination. 1.4.3 In these applications, if the benchmark comparisons indicate differences greater than 20%, the calculational model must be adjusted or a correction must be applied to reduce the difference between the fluence prediction and the upper 1-sigma limit to within 20%. For other applications, the accuracy should be determined using the approach described in Regulatory Position 1.4, and an uncertainty allowance should be included in the fluence estimate as appropriate in the specific application.
31
Regulatory Position FLUENCE MEASUREMENT METHODS
Spectrum Coverage. The set of dosimeters should provide adequate spectrum 2.1.1 coverage.
Dosimeter Nuclear and Material Properties. Use of dosimeter materials should 2.1.1 address melting, oxidation, material purity, total and isotopic mass assay, perturbations by encapsulations and thermal shields, and accurate dosimeter positioning. Dosimeter half-life and photon yield and interference should also be evaluated.
Corrections. Dosimeter-response measurements should account for fluence rate 2.1.2 variations, isotopic burnup effects, detector perturbations, self shielding, reaction interferences, and photofission.
Response Uncertainty. An uncertainty analysis must be performed for the 2.1.3 response of each dosimeter.
Validation. Detector-response calibrations must be carried out periodically in a 2.2 standard neutron field.
Fast-Neutron Fluence. The E > 1 MeV fast-neutron fluence for each 2.3 measurement location must be determined using calculated spectrum-averaged cross-sections and individual detector measurements. As an alternative, the detector responses may be used to determine reaction probabilities or average reaction rates.
Measurement-to-Calculation Ratios. The M/C ratios, the standard deviation and 2.3 bias between calculation and measurement, must be determined.
REPORTING PROVISIONS
Neutron Fluence and Uncertainties Details of the absolute fluence calculations, associated methods qualification and 3.1 fluence adjustments (if any) should be reported. Justification and a description of any deviations from the provisions of this guide should be provided.
Calculated multigroup neutron fluences and fluence rates should be reported. 3.2 The value and basis of any bias or model adjustment made to improve the 3.2 measurement-to-calculation agreement must be reported.
32
Regulatory Position Calculated integral fluences and fluence rates for E > 1 MeV and their uncertainties 3.3 should be reported.
Measured and calculated integral E > 1 MeV fluences or reaction rates and 3.4 uncertainties for each measurement location should be reported. The M/C ratios and the spectrum averaged cross-section should also be reported.
The results of the standard field validation of the measurement method should be 3.5 reported.
Specific Activities and Average Reaction Rates The specific activities at the end of irradiation and measured average reaction rates 3.5 with uncertainties should be reported.
All corrections and adjustments to the measured quantities and their justification 3.5 should be reported.
33
TABLE 2. THRESHOLD DETECTORS RECOMMENDED FOR
PRESSURE VESSEL DOSIMETRY
Nominal Applicable Threshold (MeV) ASTM Standards
237 Np(n,f)FP* 0.69 E 705 (Ref. 83)
93 Nb(n,n)93mNb 0.97 E 1297 (Ref. 84)
238 U(n,f)FP 1.45 E 704 (Ref. 85)
58 Ni(n,p)58Co 2.05 E 264 (Ref. 86)
54 Fe(n,p)54Mn 2.32 E 263 (Ref. 87)
46 Ti(n,p)46Sc 3.76 E 526 (Ref. 88)
63 Cu(n,)60Co 4.65 E 523 (Ref. 89)
- FP indicates fission product.
34
35 Cycle-Dependent Assembly/Pin Multi-Group Powers & Library Exposures (Optional)
(1.1.1) (1.1.2.2)
(1.2)
Isotopic Component Macroscopic Continuous Fission MCNP* Material Cross- Energy Spectra Source Compositions Sections Library
(1.1.1)
(1.3.2)
Core/Internals/
Vessel MCNP
Geometry Fluence
3-D Vessel Relevant Regulatory Positions are indicated.
Fluence Estimate
- The MCNP code is used as an example.
Figure 2. Monte Carlo Calculation Methodology
36
(1.4)*
Fluence Calculation Calculation Methods Methods Revision (When Necessary)
Benchmark Analytic Uncertainty Uncertainty Analysis (1.4.2)
Analysis
(1.4.1)
(1.4.2.1) (1.4.2.2) (1.4.2.3)
Sensitivity Calculation of Calculation Plant Plant Calculational Calculation Analysis Operating Simulator of Simulator Model Input Model Benchmark of Benchmark (M/MX) Reactor Model Input Benchmarks Input Model Input Problem Measurements Estimated Capsule/ Comparison Calculation Uncertainty Cavity Simulator Comparison Comparison of Results Benchmark Uncertainties Determination Measurement Measurement of Results of Results Solution
& Biases (M/MX)*X Data Data (X)
Combine Agreement Resolution of Ca , BCa Benchmark Acceptable Differences Uncertainty/
Bias Estimates
(1.4.3) Cb, BCb Combine
- Relevant Regulatory Positions are indicated.
Analytic &
Benchmark Estimates Calculational Uncertainty Based on Ca (Cb) Analytic (Benchmark) Uncertainty Analysis BCa (BCb) Calculation Bias Based on Analytic (Benchmark) Uncertainty Analysis C, B C
C (BC) Estimated Calculation Uncertainty (Bias)
Figure 3. Calculation Methodology Qualification Procedure
- Relevant Regulatory Positions are indicated.
= [Equation (6)]
Figure 4. Uncertainty-Dependent Fluence Determination
38
Fluence Measurement Methods
(2.1.3)* (2.2)
Dosimeter Dosimeter Periodic Validation Certified Fluence &
Measurement Uncertainty in Standard and/or Spectrum-Averaged Uncertainties Analysis Reference Fields Cross-Section
& Biases Ma , BMa vM , BMv Combine Uncertainty &
Bias Estimates Measurement Uncertainty Based on Analytic (Validation)
M M
M
, ( )
a M
v - Uncertainty Analysis
,B
Measurement Bias Based on M
B , (B )
a M
v - Analytic (Validation) Uncertainty Analysis
- Relevant Regulatory Positions are indicated.
Estimated Measurement M, (BM ) - Uncertainty (Bias)
Figure 5. Measurement Qualification Procedure
REFERENCES
1. J.F. Carew et al., "Application of Neutron Transport Green's Functions to the Calculation of Pressure Vessel Fluence," Nuclear Science and Engineering, Vol. 91, p. 279, 1985.
2. Standard Guide for Characterizing Neutron Exposures in Iron and Low Alloy Steels in Terms of Displacements per Atom (DPA), ASTM E693, American Society for Testing and Materials, Philadelphia, 1994.1
3. Guide to the Expression of Uncertainty in Measurement, International Organization for Standardization, ISBN 92-67-10188-9, Genéve, 1993. [ISO, Care Postal 56, CH-1211 Genéve 20, Switzerland].
4. J.F. Carew et al., "Pressure Vessel Fluence Calculation Benchmark Problems and Solutions," Draft NUREG/CR-6115 (BNL NUREG-52395), USNRC, September 1999.2
5. R.E. Maerker et al., "Applications of the LEPRICON Unfolding Procedures to the Arkansas Nuclear One-Unit 1 Reactor," Nuclear Science and Engineering, Vol. 93(2), pp. 137-170,
June 1986.
6. R.E. Maerker, "LEPRICON Analysis of Pressure Vessel Surveillance Dosimetry Inserted into H.B. Robinson-2 During Cycle 9," Nuclear Science and Engineering, Vol. 96(4), pp.
263-289, August 1987.
7. R.E. Maerker, "Analysis of the NESDIP2 and NESDIP3 Radial Shield and Cavity Experiments," NUREG/CR-4886 (Oak Ridge National Laboratory, ORNL/TM-10389),
USNRC, May 1987.3
8. C.Y. Fu and D.M. Hetrick, "Update of ENDF/B-V Mod-3 Iron: Neutron-Producing Reaction Cross-Sections and Energy Angle Correlations," ORNL/TM-9964, ENDF-341, Oak Ridge National Laboratory, July 1986.4
9. Standard Guide for Application of ASTM Evaluated Cross Section Data File, Matrix E
706 (IIB), ASTM E1018-95, American Society for Testing and Materials, Philadelphia,
1995.1
1 Available from the American Society for Testing and Materials, 1916 Race Street, Philadelphia, PA 19103-1187.
2 Copies are available for inspection or copying for a fee from the NRC Public Document Room at 11555 Rockville Pike, Rockville, MD; the PDR's mailing address is USNRC PDR, Washington, DC 20555; telephone (301)415-4737; fax (301)415-
3548.
3 Copies may be purchased at current rates from the U.S. Government Printing Office, Post Office Box 37082, Washington, DC
20013-7082 (telephone (202)512-2249 or (202)512-2171); or from the National Technical Information Service by writing NTIS
at 5285 Port Royal Road, Springfield, VA 22161. Copies are available for inspection or copying for a fee from the NRC Public Document Room at 11555 Rockville Pike, Rockville, MD; the PDR's mailing address is USNRC PDR,
Washington, DC 20555; telephone (301)415-4737; fax (301)415-3548.
4 Available from the Radiation Safety Information Computational Center, Oak Ridge National Laboratory, Post Office Box 2008, Oak Ridge, TN 37831-6362.
40
10. B. Petrovic, H.L. Hanshaw, and A. Haghighat, Evaluation of Anisotropy Effects in Pressure Vessel Fluence Calculations Using the BUGLE-93 Library, Transactions of the ANS, Vol. 71, 1994.5
11. R.W. Roussin et al., "VITAMIN-C: The CTR Processed Multigroup Cross-Section Library for Neutronics Studies, ORNL/RSIC-37 (ENDF-296), Oak Ridge National Laboratory, July 1980.4
12. R.W. Roussin et al., "VITAMIN-E: 174 Neutron, 38 Gamma-Ray Multigroup Cross- Section Library for Deriving Application-Dependent Working Libraries for Radiation Transport Calculations," DLC-113, Oak Ridge National Laboratory, November 1987.4
13. W.E. Ford, III, et al., "Modification Number One to the 100n - 21g Cross Section Library,"
ORNL/TM-5249 (Available as DLC-37D/EPR from RSIC), Oak Ridge National Laboratory, March 1976.4
14. "SAILOR: Coupled, Self-Shielded, 47-Neutron, 20-Gamma Ray, P3, Cross Section Library for Light Water Reactors," DLC-76-SAILOR, Oak Ridge National Laboratory, July 1987.4 (Available at <http://epicws.epm.ornl.gov/rsic.html> .)
15. M.L. Williams et al., "The ELXSIR Cross-Section Library for LWR Pressure Vessel Irradiation Studies: Part of the LEPRICON Computer Code System," NP-3654, Electric Power Research Institute, September 1984.6
16. D.T. Ingersoll et al., Production and Testing of the VITAMIN-B6 Fine-Group and the BUGLE-93 Broad-Group Neutron/Photon Cross-Section Libraries Derived from ENDF/B-
VI Nuclear Data, ORNL-6795, NUREG/CR-6214, January 1995; and J.E. White et al.,
BUGLE-96: A Revised Multigroup Cross-Section Library for LWR Applications Based on ENDF/B-VI Release 3, presented at the American Nuclear Society Radiation & Shielding Topical Meeting, April 21-25, 1996, Falmouth, MA, April 1996.3
17. R.E. Maerker, M.L. Williams, and B.L. Broadhead, "Accounting for Changing Source Distributions in Light Water Reactor Surveillance Dosimetry Analysis," Nuclear Science and Engineering, Vol. 94, pp. 291-308, 1986.
18. R.E. Maerker, M.L. Williams, and B.L. Broadhead, "TIMEPATCH: A Module in the LEPRICON Computer Code System for Evaluating Effects of Time-Dependent Source Distributions in PWR Surveillance Dosimetry," EPRI Interim Report, December 1985.
(Available from the Radiation Safety Information Computational Center, ORNL, as part of PSR-277/LEPRICON, or at <http://epicws.epm.ornl.gov/rsic.html> .)3
5 Available from the American Nuclear Society, 555 N. Kensington Avenue, La Grange Park, Illinois 60525.
6 Available from EPRI Research Reports Center, P.O. Box 50490, Palo Alto, CA 94303.
41
19. M. Todosow and J.F. Carew, "Evaluation of Selected Approximations Used in Pressure Vessel Fluence Calculations," Transactions of the American Nuclear Society, Vol. 46, p.
658, June 1984.5
20. W.J. Eich, "Advanced Recycle Methodology Program," Part II, Chapter 5, Research Project
118-1, Electric Power Research Institute, 1976.7
21 A. Ahlin et al., "CASMO - A Fuel Assembly Burnup Program," AE-RF-76-4158 (Rev.
Ed.), Studsvik Energiteknik AB, 1978. (Available from Studsvik of America Inc., 1087 Beacon St., Newton, MA 02159).
22. A.L. Aronson et al., "Evaluation of Methods for Reducing Pressure Vessel Fluence," BNL-
NUREG-32876, Brookhaven National Laboratory, March 1983.2
23. G.P. Cavanaugh et al., "Reduction in Reactor Vessel Irradiation Through Fuel Management," Transactions of the American Nuclear Society, Vol. 45, p. 98, October
1983.5
24. M. Todosow et al., "Pressure Vessel Fluence Reduction Through Selective Fuel Assembly Replacement," Transactions of the American Nuclear Society, Vol. 45, p. 595, October
1983.5
25. D. Cokinos et al., "Pressure Vessel Damage Fluence Reduction by Low-Leakage Fuel Management," Transactions of the American Nuclear Society, Vol. 45, p. 594, October
1983.5
26. M.L. Williams, "DOTSOR: A Module in the LEPRICON Computer Code System for Representing the Neutron Source Distribution in LWR Cores," EPRI Interim Report, December 1985. (Available from the Radiation Safety Information Computational Center, ORNL, as part of PSR-277/LEPRICON, or at <http://epicws.epm.ornl.gov/rsic.html> .)3
27. A. Haghighat, M. Mahgerefteh, and B. Petrovic, Evaluation of the Uncertainties in the Source Distribution for Pressure Vessel Neutron Fluence Calculations, Nuclear Technology, Vol. 109, 54-75, January 1995.
28. "DOT 3.5 - A Two-Dimensional Discrete Ordinate Transport Code," CCC-276, Oak Ridge National Laboratory, 1978.4 In RSICC Computer Code Collection, CCC-650, Doors3.2, ORNL, July 1998. (Available at <http://epicws.epm.ornl.gov/rsic.html> .)
29. W.A. Rhoades and R.L. Childs, "An Updated Version of the DOT 4 One-and-Two- Dimensional Neutron/Photon Transport Code," ORNL-5851, Oak Ridge National Laboratory, July 1982.4 In RSICC Computer Code Collection, CCC-650, Doors3.2, ORNL, July 1998. (Available at <http://epicws.epm.ornl.gov/rsic.html> .)
7 Copies are available from the National Technical Information Service by writing NTIS at 5285 Port Royal Road, Springfield, VA 22161; telephone (703)487-4650; <http://www.ntis.gov/ordernow>.
42
30. F.B. Kam et al., "Pressure Vessel Fluence Analysis and Neutron Dosimetry,"
NUREG/CR-5049 (Oak Ridge National Laboratory, ORNL/TM-10651), December 1987.3
31. "Standard Guide for Application of Neutron Transport Methods for Reactor Vessel Surveillance, E 706 (IID)," ASTM E482-89, American Society for Testing and Materials, Philadelphia, 1989.1
32. B.G. Petrovic and A. Haghighat, Effects of SN Method Numerics on Pressure Vessel Neutron Fluence Calculations, Nuclear Science and Engineering, 122, 167-193, February
1996.
33. B.G. Petrovic and Haghighat, Analysis of Inherent Oscillations in Multidimensional SN
Solutions of the Neutron Transport Equation, Nuclear Science and Engineering, 124, 31-
37, September 1996.
34. B.G. Petrovic and A.Haghighat, New Directional Theta-Weighted (DTW) Differencing Scheme and Reduction of Estimated Pressure Vessel Fluence Uncertainty, Proceedings of the Ninth International Symposium on Reactor Dosimetry, H. Ait Abderrahim, P. Dhondt and B. Osmera, Eds., World Scientific Publ. Co., 1998.
35. P. G. Laky and N. Tsoulfanidis, Neutron Fluence at the Pressure Vessel of a Pressurized Water Reactor Determined by the MCNP Code, Nuclear Science & Engineering, 121, 433,
1995.
36. J. C. Wagner, A. Haghighat, and B. G. Petrovic, Monte Carlo Transport Calculations and Analysis for Reactor Pressure Vessel Neutron Fluence, Nuclear Technology, 114, No. 3,
373, 1996.
37. J. C. Wagner, A. Haghighat, and B. G. Petrovic, Investigation of Pressure Vessel Fluence Calculation with Monte Carlo, Transactions of the American Nuclear Society, 68, Part A,
446, June 1993.
38. W. T. Urban et al., PCA Benchmark Solutions Using MCNP and THREEDANT, LA-
UR-93-31-21, Los Alamos National Laboratory, 1993. In ASTM STP 1228, pp. 376-383,
1994.
39. A. Avery et al., Calculations of Pressure Vessel Fluence in PWRs Using ENDF/B-VI
Data, Proceedings of the 8th International Conference on Radiation Shielding, Arlington, Texas, April 24-28, 1994, Vol. 2, p. 677, American Nuclear Society, 1994.5
40. S. Power, An Analysis of the HB Robinson Unit 2 PWR Using the Monte Carlo Code MCBEND, Proc. of the 7th ASTM-Euratom Symposium on Reactor Dosimetry, Strasbourg ASTM, 1990.1
43
41. J.C. Wagner and A. Haghighat, Automated Variance Reduction of Monte Carlo Using Discrete Ordinates Adjoint Functions, Nuclear Science & Engineering, Vol. 128, 186-208,
1998.
42. A. Haghighat et al., Performance of the Automated Adjoint Accelerated MCNP(A3MCNPTM1) for Simulation of a BWR Core Shroud Problem, Proceedings of the M&C 99 Meeting, Madrid, Spain, Vol. 2, 1381-1392, Senda Editorial, S.A., Madrid, Spain, September 27, 1999.
43. J. F. Breismeister (Ed.), MCNP - A General Monte Carlo N-Particle Transport Code, Version 4C, LA-12625, Los Alamos National Laboratory, March 1997.3
44. G. A. Wright et al., MCBEND - A Fluence Tool from AEA Technology, Reactor Dosimetry, ASTM STP 1398, J. G. Williams et al., Eds., ASTM, West Conshohoken, PA.
2000.1
45. T. E. Booth, A Sample Problem in Variance Reduction in MCNP, LA-10363-MS, Los Alamos National Laboratory, 1985.
46. L. L. Carter and E.D. Cashwell, Particle Transport Simulation with the Monte Carlo Method, ERDA Critical Review Series, TID-26607, 1975.7
47. J.F. Carew, D.K. Min, and A.L. Aronson, "Spectral Effects in the Extrapolation of Pressure Vessel Surveillance Capsule Measurements," Nuclear Technology, Vol. 55, No. 3, p. 565, December 1981.
48. M.L. Williams, P. Chowdhury, and B.L. Broadhead, DOTSYN: A Module for Synthesizing Three-Dimensional Fluxes in the LEPRICON Computer Code System, EPRI Interim Report, December 1985. (Available from the Radiation Safety Information Computational Center, ORNL, as part of PSR-277/LEPRICON, or at <http://epicws.epm.ornl.gov/rsic.html> .)3
49. R.E. Maerker, "Analysis of the VENUS-3 Experiments," NUREG/CR-5338 (Prepared for the NRC by Oak Ridge National Laboratory, ORNL/TM-11106), USNRC, August 1989.2
50. L. Lois and T. Collins, "Reactor Cavity Dosimetry - A Regulatory Perspective," Transactions of the American Nuclear Society, Vol. 63, p. 431, June 1991.4
51. N. Tsoulfanidis, "Neutron Energy Spectra in the Core and Cavity of the ANO-2 PWR," NP-
4238, Electric Power Research Institute, September 1985.5
52. B.N. Taylor and C.E. Kuyatt, "Guidelines For Evaluating and Expressing the Uncertainty of NIST Measurement Results," NIST Technical Note-1297, National Institute For Standards and Technology, January 1993.1
53. J.F. Carew et al., "Pressure Vessel Fluence Benchmark Calculations," BNL-NUREG-34715, Brookhaven National Laboratory, February 1984.3
44
54. D.M. Cokinos et al., "Benchmarking of Pressure Vessel Fluence Calculations," Transactions of the American Nuclear Society, Vol. 46, p. 636, June 1984.4
55. B.L. Broadhead et al., "LEPRICON Adjustment Module: A Generalized Linear Least Squares Data Analysis Program with Application to PWR Surveillance Dosimetry," EPRI Interim Report, March 1985. (Available from the Radiation Safety Information Computational Center, ORNL, as part of PSR-277/LEPRICON, or at <http://epicws.epm.ornl.gov/rsic.html> .)3
.
56. R.E. Maerker, B.L. Broadhead, and J.J. Wagschal, "Theory of a New Unfolding Procedure in Pressurized Water Reactor Pressure Vessel Dosimetry and Development of an Associated Benchmark Data Base," Nuclear Science and Engineering, Vol. 91(4), pp. 369-392, December
1985.
57. W.N. McElroy et al., "A Computer-Automated Iterative Method for Neutron Flux Spectral Determination by Foil Activation," AFL-TR-67-41, Vol. I (available as RSIC Program No.
CCC-112/SAND II), 1967.3
58. C.A. Oster et al., "A Modified Monte Carlo Program for SAND-II with Solution Weighing and Error Analysis," Hanford Engineering Development Laboratory, HEDL-TME 76-60, 1976.2
59. C.R. Green, J.A. Halbleib, and J.V. Walker, "A Technique for Unfolding Neutron Spectra from Activation Measurements," Sandia Corporation, SC-RR-67-746, 1967. (Available as RSIC Program No. CCC-108/SPECTRA.)4
60. F.G. Perey, "Least-Squares Dosimetry Unfolding: The Program STAY'SL," ORNL/TM-6062, Oak Ridge National Laboratory, 1977. (Available as RSIC Program No. PSR-113.)4
61. F. Schmittroth, "FERRET Data Analysis Code," HEDL-TME 79-40, Hanford Engineering Development Laboratory, 1979. (Available as RSIC Program No. PSR-145.)4
62. F.W. Stallmann, "LSL-M2: A Computer Program for Least-Squares Logarithmic Adjustment of Neutron Spectra," NUREG/CR-4349 (Prepared for NRC by Oak Ridge National Laboratory, ORNL/TM-9933), March 1986. (Available as RSIC Program No. PSR-233.)3,4
63. "Standard Guide for Application of Neutron Spectrum Adjustment Methods in Reactor Surveillance (IIA), ASTM E944-96, American Society for Testing and Materials, Philadelphia, 1996.1
64. "Standard Guide for Benchmark Testing of Reactor Dosimetry in Standard and Reference Neutron Fields, ASTM E2005-99, American Society for Testing and Materials, Philadelphia,
1999.1
65. I. Remec and F.B. Kam, "H. B. Robinson-2 Pressure Vessel Benchmark," NUREG/CR-6453 (Prepared for NRC by Oak Ridge National Laboratory, ORNL/TM-13204), February 1998.2
45
66. M.L. Williams and M. Asgari, Impact of ENDF/B-VI Cross-Section Data on H.B. Robinson Cycle 9 Dosimetry Calculations, NUREG/CR-6071, October 1993.3
67. E.P. Lippincott et al., "Evaluation of Surveillance Capsule and Reactor Cavity Dosimetry from H.B. Robinson Unit-2, Cycle 9," NUREG/CR-4576 (WCAP-11104), USNRC, February 1987.2
68. M. Austin, Sense of Direction: An Observation of Trends in Materials Dosimetry in the United Kingdom, Proceedings of the Fourth ASTM Euratom Symposium on Reactor Dosimetry, Gaithersburg, MD, March 22-26, 1982, NUREG/CP-0029, USNRC, Vol 1, August 1982.1
69. J. Butler et al., The PCA Replica Experiment Part I: Winfrith Measurements and Calibrations, AEEW-R 1736, Part I, UKAEA, Winfrith, United Kingdom, January 1984.7
70. M. D. Carter and I. J. Curl, NESTOR Shielding and Dosimetry Improvement Programme, AEEW-M 2329, 1986.
71. J. Butler et al., Review of the NESTOR Shielding and Dosimetry Improvement Programme (NESDIP), Reactor Dosimetry: Methods, Applications and Standardization, Sixth ASTM
Euratom Symposium on Reactor Dosimetry, Jackson Hole, Wyoming, May 31- June 6, 1987, STP 1001, ASTM, May 1989.1
72. W.N. McElroy, "LWR Pressure Vessel Surveillance Dosimetry Improvement Program: PCA
Experiments and Blind Test," NUREG/CR-1861 (Prepared for NRC by Hanford Engineering Development Laboratory, HEDL-TME 80-87), July 1981.2
73. I. Remec and F.B. Kam, "Pool Critical Assembly Pressure Vessel Facility Benchmark,"
NUREG/CR-6454 (Prepared for NRC by Oak Ridge National Laboratory, ORNL/TM-13205),
July 1997.2
74. F.W. Stallman et al., Reactor Calculation `Benchmark PCA Blind Test Results, ORNL/NUREG/TM-428, March 1981.2
75. D. K. Min, A. L. Aronson, and J. F. Carew, Analysis of the ORNL Pool Critical Assembly Pressure Vessel Dosimetry Benchmark Experiment, BNL-NUREG-29047, Brookhaven National Laboratory, February 1981.3
76. P. D'hondt et al., "Contributions of the Venus-Engineering Mock-Up Experiments to the LWR-
PV Surveillance," Proceedings of the 7th ASTM-Euratom Symposium on Reactor Dosimetry, Strasbourg, France, ASTM, August 1990.1
77. L. Leenders, LWR-PVS Benchmark Experiment VENUS-3 with Partial Length Shield Assemblies, SCK-CN, MOL, Belgium, FCP/VEN/01, September 1998. (See also Protection of Neutron Embrittlement in the Reactor Pressure Vessel: Venus-1 and Venus-3 Benchmarks, Nuclear Energy Agency, 2000. (Available at WWW.NEA.FR))
46
78. G. Hehn and B.C. Na, New NEA Benchmarks Reveal Decisive Improvements in Calculating Fast Neutron Fluence for Predictions of Embrittlement in Reactor Pressure Vessels, Reactor Dosimetry, ASTM STP 1398 (J.G. Williams et al., Eds.), ASTM, West Conshohoken, PA.
2000.1
79. A. Haghighat, H. Ait Abderrahim, and G.E. Sjoden, Accuracy and Performance of PENTRANTM Using the VENUS-3 Benchmark Experiment, Reactor Dosimetry, ASTM STP
1398 (J.G. Williams et al., Eds.), ASTM, West Conshohoken, PA. 2000.1
80. R.E. Maerker and B.A. Worley, "Activity and Fluence Calculations for the Startup and Two- Year Irradiation Experiments Performed at the Poolside Facility," NUREG/CR-3886 (Prepared for NRC by Oak Ridge National Laboratory, ORNL/TM-4265), October 1984.2
81. W.N. McElroy, "LWR Pressure Vessel Surveillance Dosimetry Improvement Program. PSF
Experiments Summary and Blind Test Results," NUREG/CR-3320, Vols. 1-4 USNRC,
(Prepared for the NRC by Hanford Engineering Development Laboratory), July 1986 through July 1992.2
82. G. L. Guthrie, E. P. Lippincott, E. D. McGarry, Light Water Reactor Pressure Vessel Surveillance Dosimetry Improvement Program: PSF Blind Test Workshop Minutes, Westinghouse Hanford Company, April 1994.1
83. "Standard Test Method for Measuring Reaction Rates by Radioactivation of Neptunium-237,"
ASTM E 705-90, ASTM, Philadelphia, 1991.1
84. "Standard Test Method for Measuring Fast-Neutron Reaction Rates by Radioactivation of Niobium," ASTM E1297-89, ASTM, Philadelphia, 1989.1
85. "Standard Test Method for Measuring Reaction Rates by Radioactivation of Uranium-238,"
ASTM E 704-90, ASTM, Philadelphia, 1991.1
86. "Standard Test Method for Measuring Fast-Neutron Reaction Rates by Radioactivation of Nickel," ASTM E 264-87, ASTM, Philadelphia, 1987.1
87. "Standard Test Method for Measuring Fast-Neutron Reaction Rates by Radioactivation of Iron,"
ASTM E 263-88, ASTM, Philadelphia, 1988.1
88. "Standard Test Method for Measuring Fast-Neutron Reaction Rates by Radioactivation of Titanium," ASTM E 526-87, American Society for Testing and Materials, Philadelphia, 1987.1
89. "Standard Test Method for Measuring Fast-Neutron Reaction Rates by Radioactivation of Copper," ASTM E 523-87, ASTM, Philadelphia, 1987.1
90. "Standard Guide for Sensor Design and Irradiation for Reactor Surveillance, E 706 (IIC),"
ASTM E 844-86, ASTM, Philadelphia, 1986.1
47
91. "Standard Test Method for Application and Analysis of Radiometric Monitors for Reactor Vessel Surveillance, E 706 (IIIA)," ASTM E 1005-84, ASTM, Philadelphia, 1984.1
92. "Standard General Methods for Detector Calibration and Analysis of Radionuclides," ASTM E
181-82, ASTM, Philadelphia, 1983.1
93. R. Gold et al., "Neutron Dosimetry with Solid-State Track Recorders in the Three-Mile Island Unit-2 Reactor Cavity," Nuclear Tracks, Vol. 103, p. 447, 1985.
94. F.H. Ruddy et al., "Solid-State Track Recorder Neutron Dosimetry in Light-Water Reactor Pressure Vessel Surveillance Mockups in Reactor Dosimetry,"Proceedings of the 5th ASTM-
EURATOM Symposium on Reactor Dosimetry, Geesthacht, Federal Republic of Germany, Sept.
24-28, 1984, EUR 9869, Commission of the European Communities, D. Reidel Publishing Co.,
1985.1
95. "Standard Test Method for Application and Analysis of Solid State Track Recorder (SSTR)
Monitors for Reactor Surveillance, E 706 (IIIB)," ASTM E 854-90, ASTM, Philadelphia, 1990.1
96. "Standard Test Method for Application and Analysis of Helium Accumulation Fluence Monitors for Reactor Vessel Surveillance, E 706(IIIC)," ASTM E 910 95, ASTM, Philadelphia, 1995.1
97. J.A. Grundl and C.M. Eisenhauer, Compendium of Benchmark Neutron Fields for Reactor Dosimetry, NBSIR 85-3151, National Bureau of Standards, Gaithersburg, MD, January 1986.8
98. A. Hawari et al., "Materials Dosimetry Reference Facility," Proceedings of the Eighth ASTM-
Euratom Symposium on Reactor Dosimetry, Vail, Colorado, August 1993.1
99. "Standard Guide for Benchmark Testing of Light Water Reactor Caalculations, ASTM E2006-
99, ASTM, Philadelphia, 1999.1
100. "Standard Practice for Determining Neutron Fluence, Fluence Rate and Spectra by Radioactivation Techniques, ASTM E261-98, ASTM, Philadelphia, 1998.1
8 Available from the National Institute of Standards and Technology, NIST Publication Productions, Room A635, Gaithersburg, MD 20899.
48
REGULATORY ANALYSIS
1. STATEMENT OF THE PROBLEM
The U.S. Nuclear Regulatory Commission (NRC) has promulgated regulations to ensure the structural integrity of the reactor pressure vessel for light water power reactors. Specific fracture toughness requirements for normal operation and for anticipated operational occurrences for power reactors are set forth in Appendix G, "Fracture Toughness Requirements," to 10 CFR Part 50,
"Domestic Licensing of Production and Utilization Facilities." These requirements are imposed through 10 CFR 50.60. Additionally, in response to concerns over potential pressurized thermal shock (PTS) events in pressurized water reactors (PWRs), the NRC issued 10 CFR 50.61,
"Fracture Toughness Requirements for Protection Against Pressurized Thermal Shock Events."
To satisfy the requirements of Appendix G; 10 CFR 50.61; Criterion 14, Reactor Coolant Pressure Boundary, Criterion 30, Quality of Reactor Coolant Pressure Boundary, and Criterion
31, Fracture Prevention of Reactor Coolant Pressure Boundary, of Apppendix A, General Design Criteria for Nuclear Power Plants, to 10 CFR Part 50, methods for determining the fast neutron fluence (E > MeV) are necessary to estimate the fracture toughness of the pressure vessel materials. Appendix H, "Reactor Vessel Material Surveillance Program Requirements," to 10 CFR Part 50 requires the installation of surveillance capsules, including material test specimens and flux dosimeters, to provide data on material damage correlations as a function of fluence.
The neutron fluence is attenuated by several decades between the core and vessel. This attenuation results in a strong sensitivity of the calculated vessel fluence to the physical description of the core and vessel internals and the numerical calculation of the neutron transport, and it makes an accurate determination of the pressure vessel fluence difficult. As a result, a wide range of methods of varying reliability and accuracy have been used to determine the reactor vessel fluence.
Consequently, comparisons of measured and calculated fluences have shown varying degrees of agreement, and in some cases conservatisms have been required in licensing analyses to accommodate the observed measurement-to-calculation differences.
Over the past decade, substantial improvements have been made in both the calculation and measurement of the pressure vessel fluence. These improvements have stemmed from both NRC
and industry programs. These include the development and improvement of computer codes, calculational models, measurement techniques and basic cross-section data, and the systematic qualification of the fluence methods by comparison to NRC-sponsored benchmark experiments.
These calculation and measurement improvements provide increased accuracy in the fluence determinations that are an essential part of meeting the requirements of Appendices G to 10 CFR Part 50 and 10 CFR 50.61. This is especially important for plants seeking to renew their operating licenses.
The wide variation in fluence calculation methods has resulted in lengthy plant-specific reviews and made it difficult to confirm, during the review process, that the actual fluence is adequately bounded by the various calculational methods used. This calculation and dosimetry
49
guide would provide standardized methods and procedures that would allow these reviews to be greatly simplified, and would improve confidence in the calculated fluence values.
2. OBJECTIVE
The objective of this guide is to provide state-of-the-art calculation and measurement procedures that are acceptable to the NRC staff for determining pressure vessel fluence. These procedures would yield a more accurate and reliable vessel fluence determination than the procedures that have been used in the past. The improved accuracy and realistic assessment of the uncertainty in the calculation would provide assurance that the fluence value is appropriate for use in evaluating compliance with the regulations.
3. ALTERNATIVES
The alternatives to issuing the vessel fluence calculation and dosimetry regulatory guide are as follows:
3.1 Branch Technical Position The pressure vessel fluence methods provided by the regulatory guide could be included in a branch technical position. However, this is not considered an acceptable alternative since the branch position does not provide the same high level of input and participation, especially public, industry, and the ACRS input, that the vessel fluence analysis requires.
3.2 NUREG-Series Report The vessel fluence procedures could be published in a NUREG-series report. However, these reports also do not receive the required input and participation from the public, industry, and the ACRS and they are not appropriate for providing regulatory guidance.
3.3 Discussions with Licensees The detailed fluence calculational and measurement methods that are considered acceptable to the NRC staff could be provided to the licensees through individual reviews and discussions on a case-by-case basis. This alternative is basically the same as the current practice and is equivalent to taking no action. Individual licensee discussions are extremely time-consuming for both the NRC staff and the licensee, they lead to highly individual analyses and reviews, and they do not result in an established standard.
4. COSTS AND BENEFITS
4.1 Benefits The methods described in this guide may be used for all fluence determinations used in vessel fracture toughness evaluations, including the determination of the fluence used in calculating the pressure vessel material values of RTNDT specified in Appendix G to 10 CFR Part 50 imposed by 10 CFR 50.60, and the values of RTPTS specified in 10 CFR 50.6
1. The regulatory
50
guide would improve the accuracy and reliability of these evaluations and provide consistent and reliable uncertainty estimates by incorporating state-of-the-art methods and procedures for determining the fluence and the fluence uncertainty. The guide would also ensure the completeness of licensee vessel fluence submittals and improve the efficiency of staff reviews.
The improved fluence determination will provide fluence and uncertainty estimates that are more reliable and understandable. Thus, it will provide more reliable and accurate information for the PTS screening criteria of 10 CFR 50.61 and improve application of this rule. In this regard it is noted that, for a pressure vessel near the PTS screening criteria of 10 CFR 50.61, a 25% reduction in calculated end-of-license fluence, which is typical of existing uncertainties, will reduce the calculated vessel failure frequency by approximately a factor of three.
4.2 NRC Costs The NRC costs for reviewing fluence-related submittals would be reduced substantially by the issuance of this guide. For estimating the costs, it is assumed that of the ~60 PWRs, half have a RTPTS within ~40-F of the PTS screening criterion (or other temperature limit) at the end of license and will require a detailed review. Assuming each submittal requires a staff week and only half of the PWR owners submit revised fluence analyses, the total NRC cost is approximately 15 staff weeks. If the licensees used the methods given in the guide, this cost could be reduced to approximately 2 staff weeks.
4.3 Licensee Costs Increased costs to the licensee would result from changes in the fluence calculation and measurement procedures. The calculational costs would be one-time costs and have been estimated in Table RA-1. The licensee costs resulting from the changes in the measurement procedures have been estimated in Table RA-2.
5. DECISION RATIONALE
It is recommended that the proposed regulatory guide be issued because (1) a high level of participation of the NRC, ACRS, industry, and the public is reflected in the guide, (2) a methodology standard would be established, and (3) inefficient use of NRC staff and licensee resources during the review process would be eliminated.
The alternatives identified above for providing acceptable fluence methods to the licensees do not provide the advantages listed above associated with issuance of this regulatory guide. In particular, while these alternative approaches may result in the same or slightly increased cost to licensees, they result in a highly inefficient use of the NRC staff resources. The alternatives to a regulatory guide are therefore judged to be unacceptable.
51
TABLE RA-1. ADDITIONAL LICENSEE CALCULATION COSTS
Calculation Tasks Staff Weeks (a) Modifications to calculational models +2 (b) Additional calculation benchmarking and +8 qualification (c) Calculation uncertainty analysis +3 (d) Calculation documentation and reporting +2 (e) Reduced Licensing Activities* -4 Total Additional Licensee Cost +11 staff weeks TABLE RA-2. ADDITIONAL LICENSEE MEASUREMENT COSTS
Measurement Task Staff Weeks (a) Additional Quality Control +1 (b) Dosimeter Response Corrections +1 (c) Periodic Detector Calibration +2 (d) Response Uncertainty Analysis +2 (e) Additional Measurement +1 Documentation and Reporting Total Additional Licensee Cost + 7 staff weeks
- This estimate reflects the reduced licensee costs from not attending meetings with the NRC staff and responding to questions as a result of following the procedures in the regulatory guide.
52
BACKFIT ANALYSIS
The regulatory guide does not require a backfit analysis as described in 10 CFR 50.109(c)
because it does not impose a new or amended provision in the Commission rules or a regulatory staff position interpreting the Commission rules that is either new or different from a previous applicable staff position. In addition, this regulatory guide does not require the modification or addition to systems, structures, components, or design of a facility or the procedures or organization required to design, construct, or operate a facility. Rather a licensee or applicant can select a preferred method for achieving compliance with license or the rules or the orders of the Commission as described in 10 CFR 50.109(a)(7). The regulatory guide provides the opportunity to use the methods described in the guide for all fluence determinations used in vessel fracture toughness evaluations, including the determination of the fluence used in calculating the pressure vessel materials' values of RTNDT for use in 10 CFR Part 50, Appendix G (imposed by 10 CFR 50.60), and the values of RTPTS in accordance with 10 CFR 50.61.