Regulatory Guide 1.77: Difference between revisions

From kanterella
Jump to navigation Jump to search
(Created page by program invented by StriderTol)
(StriderTol Bot change)
 
Line 14: Line 14:
| page count = 8
| page count = 8
}}
}}
{{#Wiki_filter:May 1974 U.S. ATOMIC, ENERGY COMMISSION
{{#Wiki_filter:U.S. ATOMIC, ENERGY COMMISSION  
                              REGULATORY                                                                                     GUIDE
REGULATORY  
                              DIRECTORATE OF REGULATORY STANDARDS
DIRECTORATE OF REGULATORY STANDARDS
                                                                REGULATORY GUIDE 1.77 ASSUMPTIONS USED FOR EVALUATING A CONTROL ROD
May 1974 GUIDE
                            EJECTION ACCIDENT FOR PRESSURIZED WATER REACTORS
REGULATORY GUIDE 1.77 ASSUMPTIONS USED FOR EVALUATING A CONTROL ROD  
EJECTION ACCIDENT FOR PRESSURIZED WATER REACTORS


==A. INTRODUCTION==
==A. INTRODUCTION==
reactor is normally limited by the-design of the control rod system to a value well below that which would result Section 50.34, 'Contents of applications: technical                         in serious damage to the reactor system. However, a information," of 10 CFR Part 50, "Licensing of Pro                                   postulated failure of the control rod system provides the duction and Utilization Facilities," requires that each                             potential for a relatively high rate of reactivity insertion application for a construction permit or operating                                   which, if large enough, could cause a prompt power license provide an analysis and evaluation of the design                             burst. For U0 2 fuel, a large fraction of this generated and performance of structures, systems, and components                               nuclear energy is stored momentarily in the fuel and of the facility with the objective of assessing the                                 then released to the rest of the system. If the fuel energy potential risk to public health and safety resulting from                           densities were high enough, there would exist the operation of the facility. General Design Criterion 28,                             potential for prompt rupture of fuel pins and the
Section 50.34, 'Contents of applications: technical information," of 10 CFR Part 50, "Licensing of Pro duction and Utilization Facilities," requires that each application for a construction permit or operating license provide an analysis and evaluation of the design and performance of structures, systems, and components of the facility with the objective of assessing the potential risk to public health and safety resulting from operation of the facility. General Design Criterion 28,  
"Reactivity Limits," of Appendix A, "General Design                                 consequent rapid heat transfer to the water from finely Criteria for Nuclear Power Plants," to 10 CFR Part 50,                               dispersed molten U0 2 . Prompt fuel element rupture is requires the reactivity control system to be designed                               defined herein as a rapid increase in internal fuel rod with appropriate limits on the potential amount and rate                             pressure due to extensive fuel melting, followed by rapid of reactivity increase to assure that the effects of                                 fragmentation and dispersal of fuel cladding into the postulated reactivity accidents can neither result in                               coolant. This is accompanied by the conversion of damage to the reactor coolant pressure boundary greater                             nuclear energy, deposited as overpower heat in the fuel than limited local yielding nor sufficiently disturb the                             and in the coolant, to mechanical energy which, in core, its support structures, or other reactor pressure                             sufficient quantity, could conceivably disarrange the vessel internals to impair significantly the capability to                           reactor core or breach the primary system.
"Reactivity Limits," of Appendix A, "General Design Criteria for Nuclear Power Plants," to 10 CFR Part 50,  
requires the reactivity control system to be designed with appropriate limits on the potential amount and rate of reactivity increase to assure that the effects of postulated reactivity accidents can neither result in damage to the reactor coolant pressure boundary greater than limited local yielding nor sufficiently disturb the core, its support structures, or other reactor pressure vessel internals to impair significantly the capability to cool the core. General Design Criterion 28 also requires that these postulated reactivity accidents include con sideration of the rod ejection accident unless such an accident is prevented by positive means.


cool the core. General Design Criterion 28 also requires that these postulated reactivity accidents include con                                      The Regulatory staff has reviewed the available sideration of the rod ejection accident unless such an                              experimental information concerning fuel failure thresh accident is prevented by positive means.                                            olds. In general, failure consequences for U0 2 have been insignificant below 300 cal/g for both irradiated and This guide identifies acceptable analytical methods                         unirradiated fuel rods. Therefore, a calculated radial and assumptions that may be used in evaluating the                                   average energy density of 280 cal/g at any axial fuel consequences of a rod ejection accident in uranium                                   location in any fuel rod as a result of a postulated rod oxide-fueled pressurized water reactors (PWRs). In some                             ejection accident provides a conservative maximum limit cases, unusual site characteristics, plant design features,                           to ensure that core damage will be minimal and that or other factors may require different assumptions                                   both short-term andlong-term core cooling capability which will be considered on an individual basis. The                                 will not be impaired.
This guide identifies acceptable analytical methods and assumptions that may be used in evaluating the consequences of a rod ejection accident in uranium oxide-fueled pressurized water reactors (PWRs). In some cases, unusual site characteristics, plant design features, or other factors may require different assumptions which will be considered on an individual basis. The Advisory Committee on Reactor Safeguards has been consulted concerning this guide and has concurred in the regulatory position.


Advisory Committee on Reactor Safeguards has been consulted concerning this guide and has concurred in the regulatory position.                                                                       For the postulated control rod ejection accident, a  
==B. DISCUSSION==
The rate at which reactivity can be inserted into the core of a uranium oxide-fueled water-cooled power reactor is normally limited by the-design of the control rod system to a value well below that which would result in serious damage to the reactor system. However, a postulated failure of the control rod system provides the potential for a relatively high rate of reactivity insertion which, if large enough, could cause a prompt power burst. For U0 2 fuel, a large fraction of this generated nuclear energy is stored momentarily in the fuel and then released to the rest of the system. If the fuel energy densities were high enough, there would exist the potential for prompt rupture of fuel pins and the consequent rapid heat transfer to the water from finely dispersed molten U0 2. Prompt fuel element rupture is defined herein as a rapid increase in internal fuel rod pressure due to extensive fuel melting, followed by rapid fragmentation and dispersal of fuel cladding into the coolant. This is accompanied by the conversion of nuclear energy, deposited as overpower heat in the fuel and in the coolant, to mechanical energy which, in sufficient quantity, could conceivably disarrange the reactor core or breach the primary system.
 
The Regulatory staff has reviewed the available experimental information concerning fuel failure thresh olds. In general, failure consequences for U0 2 have been insignificant below 300 cal/g for both irradiated and unirradiated fuel rods. Therefore, a calculated radial average energy density of 280 cal/g at any axial fuel location in any fuel rod as a result of a postulated rod ejection accident provides a conservative maximum limit to ensure that core damage will be minimal and that both short-term andlong-term core cooling capability will not be impaired.
 
For the postulated control rod ejection accident, a mechanical failure of a control rod mechanism housing is assumed such that the reactor coolant system pressure would eject the control rod and drive shaft to the fully withdrawn position.


==B. DISCUSSION==
USAEC REGULATORY GUIDES
mechanical failure of a control rod mechanism housing is assumed such that the reactor coolant system pressure The rate at which reactivity can be inserted into the                         would eject the control rod and drive shaft to the fully core of a uranium oxide-fueled water-cooled power                                    withdrawn position.
Copies of published guides may be obtained by request Indicating the divisions desired to the US. Atomic Energy Commission. Washington, D.C. 20545.
 
Regulatory Guides we issued to describe and' make available to the public Attention: Director of Regulatory Standards. Comments and suggestions for methods acceptale to the AEC Regulatory staff of implemnenting specfifc parts of Improvements In thes guldes are encouraged and should be sent to the Secretary the Commission's regulations, to delineate techniques used -by the staff In of the Commission. US. Atomic Energy Commission, Washington, D.C. 20645 evaluating specific problems or postulated accidents, or to provide guidance to Attention. Chief, Public Proceedings Staff.
 
applicants. 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 era issued in the following ten broad divisions:
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.
 
1. Power Rleactors S. Products
2. Research and Test Reactors


USAEC REGULATORY GUIDES                                      Copies of published guides may be obtained by request Indicating the divisions desired to the US. Atomic Energy Commission. Washington, D.C. 20545.
===7. Transportation ===
3. Fuels and Materials Facilities
8. Occupational Health Published guides will be revised periodically, as appropriate, to accommodate
4. Environmental and Siting
9. Antitrust Review comments and to rflect new Information or experience.


Regulatory Guides we issued to describe and' make available to the public          Attention: Director of Regulatory Standards. Comments and suggestions for Improvements In thes guldes are encouraged and should be sent to the Secretary methods acceptale to the AEC Regulatory staff of implemnenting specfifc parts of the Commission's regulations, to delineate techniques used -by the        staff In  of the Commission. US. Atomic Energy Commission, Washington, D.C. 20645 evaluating specific problems or postulated accidents, or to provide guidance to    Attention. Chief, Public Proceedings Staff.
L. Materials Ond Plant Protection
1


applicants. Regulatory Guides are not substitutes for regulations and compliance with them Isnot required. Methods and solutions different from those set out In    The guides era issued in the following ten broad divisions:
===0. General===
the guides will be acceptable If they provide a basis for the findings requisite to  1. Power Rleactors                        S. Products the Issuance or continuance of a permit or license by the Commission.                2. Research and Test Reactors             


===7. Transportation===
A sufficient number of initial reactor states to completely bracket all possible operational conditions of interest should be analyzed to assure examination of upper bounds on ultimate damage. In areas of uncer tainty, the appropriate minimum or maximum para meters relative to nominal or expected values should be used to assure a conservative evaluation. The initial reactor states should include consideration of at least the following:
                                                                                      3. Fuels and Materials Facilities         
Zero power (hot standby) - Beginning of Life (BOL)
and End of Life (EOL);
Low power
- BOL and EOL;
Ful


===8. Occupational Health===
====l. power ====
                                                                                      4. Environmental and Siting                9. Antitrust Review Published guides will be revised periodically, asappropriate, to accommodate        L. Materials Ond Plant Protection        10. General comments and to rflect new Information or experience.
- BOL and EOL.


A sufficient number of initial reactor states to            consequences of this accident for a pressurized water completely bracket all possible operational conditions of       reactor.
The effects of the loss of primary system integrity as a result of the failed control rod housing should be included in the analysis. It should- also be shown that failure of one control rod housing will not lead to failure of other control rod housings.


interest should be analyzed to assure examination of upper bounds on ultimate damage. In areas of uncer                             
The approach that should be used in the radiological analysis of a control rod ejection accident is to deter mine the amount of each gaseous radionuclide released to the primary containment and, with this information in conjunction with the procedures set forth in Appendix B of this guide, to determine the radiological consequences of this accident for a pressurized water reactor.


==C. REGULATORY POSITION==
==C. REGULATORY POSITION==
tainty, the appropriate minimum or maximum para meters relative to nominal or expected values should be                Acceptable assumptions and evaluation models for used to assure a conservative evaluation. The initial              analyzing a rod. ejection accident in PWRs are presented reactor states should include consideration of at least the      in Appendices A (Physics and Thermal-Hydraulics) and following:                                                        B (Radiological Assumptions) of this guide. By use of these appendices, it should be shown that:
Acceptable assumptions and evaluation models for analyzing a rod. ejection accident in PWRs are presented in Appendices A (Physics and Thermal-Hydraulics) and B (Radiological Assumptions) of this guide. By use of these appendices, it should be shown that:  
      Zero power (hot standby) - Beginning of Life (BOL)
1. Reactivity excursions will not result in a radial average fuel enthalpy greater than 280 cal/g at any axial location in any fuel rod.
                                  and End of Life (EOL);          1. Reactivity excursions will not result in a radial Low power                  - BOL and EOL;                  average fuel enthalpy greater than 280 cal/g at any axial Full.power                - BOL and EOL.                  location in any fuel rod.


The effects of the loss of primary system integrity        2. Maximum reactor pressure during any portion of as a result of the failed control rod housing should be          the assumed transient will be less than the value that will included in the analysis. It should- also be shown that          cause-stresses to-exceed the Emergency Condition-stress failure of one control rod housing will not lead to failure      limits as defined in Section III of the ASME Boiler and of other control rod housings.                                    Pressure Vessel Code.'
2. Maximum reactor pressure during any portion of the assumed transient will be less than the value that will cause-stresses to-exceed the Emergency Condition-stress limits as defined in Section III of the ASME Boiler and Pressure Vessel Code.'  
      The approach that should be used in the radiological        3. Offsite dose consequences will be well within the analysis of a control rod ejection accident is to deter          guidelines of 10 CFR Part 100, "Reactor Site Criteria."
3.
mine the amount of each gaseous radionuclide released to the primary containment and, with this information                  1 Copies may be -obtained from the American Society of in conjunction with the procedures set forth in                  Mechanical Engineers, United Engineering Center, 345 East 47th Appendix B of this guide, to determine the radiological          Street, New York, New York 10017.


-     1.77-2
Offsite dose consequences will be well within the guidelines of 10 CFR Part 100, "Reactor Site Criteria."
1Copies may be -obtained from the American Society of Mechanical Engineers, United Engineering Center, 345 East 47th Street, New York, New York 10017.
 
-
1.77-2


APPENDIX A
APPENDIX A
                                    PHYSICS AND THERMAL-HYDRAULICS
PHYSICS AND THERMAL-HYDRAULICS
                                                                the magnitude of the transient. Pressure and temperature The assumptions described below should be applied are mainly significant with respect to their effect on the in evaluating the physics and thermal-hydraulic behavior amount of reactivity inserted if there exists a positive of the reactor system for a control rod ejection accident.
The assumptions described below should be applied in evaluating the physics and thermal-hydraulic behavior of the reactor system for a control rod ejection accident.
 
I.
 
The ejected rod worth should be calculated based on the maximum worth rod resulting from the following conditions: (a) all control banks at positions corres ponding to values for maximum allowable bank inser tions at a given power level and (b) additional fully or partially inserted misaligned or inoperable rod or rods if allowed by operating procedures. Sufficient parametric studies should be performed to determine the worth of the most reactive control rod in each rod group for different control rod configurations, both expected and unexpected. The worth of single rods in rod groups should be evaluated during startup physics tests and compared with values used in the rod ejection analysis.
 
The accident should be reanalyzed if the rod worths used in the initial analysis are found to be noncon servative. Calculated rod worths should be increased, if necessary, to account for calculational uncertainties in parameters such as neutron cross sections and power asymmetries due to xenon oscillations.


moderator coefficient.
2. The reactivity insertion rate due to an ejected rod should be determined from differential control rod worth curves and calculated transient rod position versus time curves. If differential rod worth curves are not available for the reactor state of interest, conservatism should be included in the calculation of reactivity insertion through consideration of the nonlinearity in reactivity addition as the rod passes through the active core. The rate of ejection should be calculated based on the maximum pressure differential and the weight and'
cross-sectional area of the control rod and drive shaft, assuming no pressure barrier restriction.


I. The ejected rod worth should be calculated based
3.
                                                                  5. ' The fuel thermal properties such'"as fuel-clad gap on the maximum worth rod resulting from the following heat transfer coefficient and fuel thermal conductivity conditions: (a) all control banks at positions corres
                                                                .should be conservatively chosen, depending upon the ponding to values for maximum allowable bank inser transient phenomenon being investigated. For conditions tions at a given power level and (b) additional fully or          of a zero or positive moderator coefficient (usually at partially inserted misaligned or inoperable rod or rods if        beginning of life), for example, high heat transfer allowed by operating procedures. Sufficient parametric            parameters would reduce the Doppler feedback and studies should be performed to determine the worth of            increase any positive moderator feedback effects and the most reactive control rod in each rod group for              hence tend to increase the magnitude of the reactivity different control rod configurations, both expected and          transient. For a negative moderator coefficient, high unexpected. The worth of single rods in rod groups                heat transfer parameters could cause the magnitude of should be evaluated during startup physics tests and              the transient to decrease if a given quantity of heat compared with values used in the rod ejection analysis.          produces more feedback in the moderator than in the The accident should be reanalyzed if the rod worths              fuel. In the consideration of pressure pulses which may used in the initial analysis are found to be noncon              be generated, high moderator heating rates could cause servative. Calculated rod worths should be increased, if          significant pressure gradients to develop in the moder necessary, to account for calculational uncertainties in          ator channels. In computing the average enthalpy of the parameters such as neutron cross sections and power              hottest fuel pellet during the excursion for power cases, asymmetries due to xenon oscillations.                            low heat transfer would be conservative.


2. The reactivity insertion rate due to an ejected rod            6. The specific heat of U0 2 has been determined should be determined from differential control rod                experimentally and is a deterministic factor in the worth curves and calculated transient rod position versus calculated amount of stored energy (enthalpy) in the time curves. If differential rod worth curves are not fuel. Recommended values in the range of 25 to 902fC
The calculation of effective delayed neutron frac tion (Beff) and prompt neutron lifetime (2*) should be based on the well-known definitions resulting from per turbation theory, such as those described by Henry (Ref.
available for the reactor state of interest, conservatism are the data reported by Moore and Kelly (Ref. 2). In should be included in the calculation of reactivity the range of 900 to 2842!C, the data obtained by Hein insertion through consideration of the nonlinearity in and Flagella (Ref. 3), Leibowitz, Mishler, and Chasanov reactivity addition as the rod passes through the active (Ref. 4), and Chasanov (Ref. 5) are recommended for core. The rate of ejection should be calculated based on the heat capacity of the fuel. These recommended values the maximum pressure differential and the weight and'            are for clean core conditions. Possible variation in the cross-sectional area of the control rod and drive shaft, specific heat due to bumup should be investigated and assuming no pressure barrier restriction.


appropriate values used, if necessary.
1), using available experimental delayed neutron data and averaging by the fraction of fission in the various fissionable materials. In cases where the accident is quite sensitive to Peff (where the ejected rod worth >Peff), the minimum calculated value for the given reactor state should be used. For smaller transients, conservatism in the value should include consideration of not only the initial power rise (which increases with decreasing P), but also the power reduction after the trip. Similar don siderations should also be applied to determine an appropriately conservative value of 2* to be used.


3. The calculation of effective delayed neutron frac
4. The initial reactor coolant 1 pressure, core inlet temperature, and flow rate used in the analysis should be conservatively chosen with respect to their influence on the magnitude of the transient. Pressure and temperature are mainly significant with respect to their effect on the amount of reactivity inserted if there exists a positive moderator coefficient.
                                                                    7. The moderator reactivity coefficients due to voids, tion (Beff) and prompt neutron lifetime (2*) should be coolant pressure changes, and coolant temperature based on the well-known definitions resulting from per            changes should be calculated based on the various turbation theory, such as those described by Henry (Ref.          assumed conditions of the fuel and moderator using
  1), using available experimental delayed neutron data              standard transport and diffusion theory codes. If no and averaging by the fraction of fission in the various            three-dimensional space-time kinetics calculation is per fissionable materials. In cases where the accident is quite        formed, the reactivity feedback due to these coefficients sensitive to Peff (where the ejected rod worth >Peff), the        should be conservatively weighted to account for the minimum calculated value for the given reactor state              variation in their spatial Importance in the missing should be used. For smaller transients, conservatism in          'dimension(s). If boric acid shim is used in the moder the value should include consideration of not only the             ator, the highest boron concentration corresponding to initial power rise (which increases with decreasing P), but        the initial reactor state should be assumed.


also the power reduction after the trip. Similar don siderations should also be applied to determine an appropriately conservative value of 2*to be used.                  8. The Doppler coefficient should be calculated based on the effective resonance integrals and should include
5. ' The fuel thermal properties such'"as fuel-clad gap heat transfer coefficient and fuel thermal conductivity
  4. The initial reactor coolant 1pressure, core inlet              corrections for pin shadowing (Dancoff correction).
.should be conservatively chosen, depending upon the transient phenomenon being investigated. For conditions of a zero or positive moderator coefficient (usually at beginning of life), for example, high heat transfer parameters would reduce the Doppler feedback and increase any positive moderator feedback effects and hence tend to increase the magnitude of the reactivity transient. For a negative moderator coefficient, high heat transfer parameters could cause the magnitude of the transient to decrease if a given quantity of heat produces more feedback in the moderator than in the fuel. In the consideration of pressure pulses which may be generated, high moderator heating rates could cause significant pressure gradients to develop in the moder ator channels. In computing the average enthalpy of the hottest fuel pellet during the excursion for power cases, low heat transfer would be conservative.
                                                                      Calculations of the Doppler coefficient of reactivity temperature, and flow rate used in the analysis should be should be based on and should compare conservatively conservatively chosen with respect to their influence on
                                                              1.77-3


with available experimental data such as those of Hell.          characteristics and changes in flux shapes should be strand (Ref. 6). Since the Doppler coefficient reflects the       investigated, and the conservatism of the flux shapes change in reactivity as a function of fuel temperature,          used for reactivity input and feedback, peak energy uncertainties in predicting fuel temperatures at different        deposition, total energy, and gross heat transfer to the power levels should be reflected by conservatism in the          coolant should be evaluated. Also, sensitivity studies on applied value of the Doppler coefficient. If no three            variations of the Doppler effect, power distribution, fuel dimensional space-time kinetics calculation is performed,        element heat transfer parameters, and other relevant the reactivity effect of spatially weighting the core             parameters should be included.
6. The specific heat of U0 2 has been determined experimentally and is a deterministic factor in the calculated amount of stored energy (enthalpy) in the fuel. Recommended values in the range of 25 to 902fC
are the data reported by Moore and Kelly (Ref. 2). In the range of 900 to 2842!C, the data obtained by Hein and Flagella (Ref. 3), Leibowitz, Mishler, and Chasanov (Ref. 4), and Chasanov (Ref. 5) are recommended for the heat capacity of the fuel. These recommended values are for clean core conditions. Possible variation in the specific heat due to bumup should be investigated and appropriate values used, if necessary.


average temperature rise in both the axial and radial directions should be calculated.                                  13. The pressure surge should be calculated on the basis of conventional heat transfer from the fuel,. a conserva
7.
9. Control rod reactivity insertion during trip versus            tive metal-water reaction threshold, and prompt heat time should be obtained by combining the differential            generation in the coolant to determine the variation of rod worth curve with a rod velocity curve based on                heat flux with time and the volume surge. The volume maximum design limit values for scram insertion times.            surge should then be used in the calculation of the If the rod worth curve (reactivity vs. depth of insertion)        pressure transient, taking into account fluid transport in is not obtained from a "true" representation (i.e., an x,          the system, heat transfer to the steam generators, and y, z, t or an r, z, t calculation); the conservatism of the        the action of the pressurizer relief and.safety valves. No approximate calculation should be shown. The dif                  credit should be taken for the possible pressure reduc ference in the depth of insertion at zero power and at            tion caused by the assumed failure of the control rod full power should be accounted for in calculating the              pressure housing.


available scram reactivity.
The moderator reactivity coefficients due to voids, coolant pressure changes, and coolant temperature changes should be calculated based on the various assumed conditions of the fuel and moderator using standard transport and diffusion theory codes. If no three-dimensional space-time kinetics calculation is per formed, the reactivity feedback due to these coefficients should be conservatively weighted to account for the variation in their spatial Importance in the missing
'dimension(s). If boric acid shim is used in the moder ator, the highest boron concentration corresponding to the initial reactor state should be assumed.
 
8.
 
The Doppler coefficient should be calculated based on the effective resonance integrals and should include corrections for pin shadowing (Dancoff correction). 
Calculations of the Doppler coefficient of reactivity should be based on and should compare conservatively
1.77-3
 
with available experimental data such as those of Hell.
 
strand (Ref. 6). Since the Doppler coefficient reflects the change in reactivity as a function of fuel temperature, uncertainties in predicting fuel temperatures at different power levels should be reflected by conservatism in the applied value of the Doppler coefficient. If no three dimensional space-time kinetics calculation is performed, the reactivity effect of spatially weighting the core average temperature rise in both the axial and radial directions should be calculated.
 
9. Control rod reactivity insertion during trip versus time should be obtained by combining the differential rod worth curve with a rod velocity curve based on maximum design limit values for scram insertion times.
 
If the rod worth curve (reactivity vs. depth of insertion)
is not obtained from a "true" representation (i.e., an x, y, z, t or an r, z, t calculation); the conservatism of the approximate calculation should be shown. The dif ference in the depth of insertion at zero power and at full power should be accounted for in calculating the available scram reactivity.
 
10. The reactor trip delay time, or the amount of time which elapses between the instant the sensed parameter (e.g., pressure or neutron flux) reaches the level for which protective action is required and the onset of negative reactivity insertion, should be based on maxi mum values of the following: (a) time required for instrument channel to produce a signal, (b) time for the trip breaker to open, (c) time for the coil to release the rods, and (d) time required before scram rods enter the core if the tips lie above the core-reflector interface.
 
11. The computer code used for calculating the tran sient should be a coupled thermal, hydrodynamic, and nuclear model with the following capabilities: (a) incor poration of all major reactivity feedback mechanis ms, (b) at least six delayed neutron groups, (c) both axial and radial segmentation of the fuel element, (d) coolant flow provision, and (e) control rod scram initiation on either coolant system pressure or neutron flux.
 
12. The analytical models and computer codes used should be documented and justified and the con servatism of the models and codes should be evaluated both by comparison with experiment, as available, and with more sophisticated spatial kinetics codes. In par ticular, the importance of two- or three-dimensional flux characteristics and changes in flux shapes should be investigated, and the conservatism of the flux shapes used for reactivity input and feedback, peak energy deposition, total energy, and gross heat transfer to the coolant should be evaluated. Also, sensitivity studies on variations of the Doppler effect, power distribution, fuel element heat transfer parameters, and other relevant parameters should be included.
 
13. The pressure surge should be calculated on the basis of conventional heat transfer from the fuel,. a conserva tive metal-water reaction threshold, and prompt heat generation in the coolant to determine the variation of heat flux with time and the volume surge. The volume surge should then be used in the calculation of the pressure transient, taking into account fluid transport in the system, heat transfer to the steam generators, and the action of the pressurizer relief and.safety valves. No credit should be taken for the possible pressure reduc tion caused by the assumed failure of the control rod pressure housing.


14. The number of fuel rods experiencing clad failure.
14. The number of fuel rods experiencing clad failure.


10. The reactor trip delay time, or the amount of time            should be calculated and used to obtain the amount of which elapses between the instant the sensed parameter            contained fission product inventory released to the (e.g., pressure or neutron flux) reaches the level for            reactor coolant system. It should be assumed that clad which protective action is required and the onset of              failure occurs if the heat flux equals or exceeds the value negative reactivity insertion, should be based on maxi            corresponding to the onset of the transition from mum values of the following: (a) time required for                nucleate to film boiling (DNB), or for other appropriate instrument channel to produce a signal, (b) time for the          causes.
should be calculated and used to obtain the amount of contained fission product inventory released to the reactor coolant system. It should be assumed that clad failure occurs if the heat flux equals or exceeds the value corresponding to the onset of the transition from nucleate to film boiling (DNB), or for other appropriate causes.


trip breaker to open, (c) time for the coil to release the                                                                    K.,
The margin to DNB is expressed in terms of a departure from nucleate boiling ratio (DNBR). The DNBR at any position in the hottest channel is the ratio of the DNB heat flux to the actual heat flu
rods, and (d) time required before scram rods enter the core if the tips lie above the core-reflector interface.                The margin to DNB is expressed in terms of a departure from nucleate boiling ratio (DNBR). The
11. The computer code used for calculating the tran              DNBR at any position in the hottest channel is the ratio sient should be a coupled thermal, hydrodynamic, and              of the DNB heat flux to the actual heat flu


====x. The DNB====
====x. The DNB ====
nuclear model with the following capabilities: (a) incor          heat flux should be evaluated using correlations based on poration of all major reactivity feedback mechanis ms,            recognized studies and experimental heat transfer DNB
heat flux should be evaluated using correlations based on recognized studies and experimental heat transfer DNB  
(b) at least six delayed neutron groups, (c) both axial            data. A minimum DNBR should be determined from the and radial segmentation of the fuel element, (d) coolant          evaluation of the experimental data to ensure a 95%
data. A minimum DNBR should be determined from the evaluation of the experimental data to ensure a 95%  
flow provision, and (e) control rod scram initiation on            probability with a 95% confidence level that DNB has either coolant system pressure or neutron flux.                    not occurred for the fuel element being evaluated. One example of a correlation which has been used to date is
probability with a 95% confidence level that DNB has not occurred for the fuel element being evaluated. One example of a correlation which has been used to date is given by Tong (Ref. 7). The use of this correlation'and the above probabilities and confidence level yields a mini mum DNBR of 1.30. Other DNB or clad failure correlations may be used if they are adequately justified by analytical methods and supported by sufficient experimental data.
12. The analytical models and computer codes used                  given by Tong (Ref. 7). The use of this correlation'and the should be documented and justified and the con                    above probabilities and confidence level yields a mini servatism of the models and codes should be evaluated              mum DNBR of 1.30. Other DNB or clad failure both by comparison with experiment, as available, and              correlations may be used if they are adequately justified with more sophisticated spatial kinetics codes. In par            by analytical methods and supported by sufficient ticular, the importance of two- or three-dimensional flux          experimental data.


REFERENCES
REFERENCES
I. A. F. Henry, "Computation of Parameters Appear                 3. R. A. Hein and P. N. Flagella, "Enthalpy ing in the Reactor Kinetics Equations," WAPD-142,                  Measurements of U0 2 and Tungsten to 32600 K,"
I. A. F. Henry, "Computation of Parameters Appear ing in the Reactor Kinetics Equations," WAPD-142, December 1955.
    December 1955.                                                    GEMP-578, February 1968.
 
2. G. E. Moore and K. K. Kelley, "High Temperature Heat Contents of Uranium, Uranium Dioxide and Uranium Trioxide," Journal of the American Chemical Society, Vol. 69, p. 2105, 1949.
 
3. R.
 
A.
 
Hein and P. N. Flagella, "Enthalpy Measurements of U0 2 and Tungsten to 32600K,"  
GEMP-578, February 1968.


2. G. E. Moore and K. K. Kelley, "High Temperature                4. L Leibowitz, L. W. Mishler, and M. G. Chasanov, Heat Contents of Uranium, Uranium Dioxide and                      "Enthalpy of Solid Uranium Dioxide from 25000K
4. L Leibowitz, L. W. Mishler, and M. G. Chasanov,  
    Uranium Trioxide," Journal of the American                        to its Melting Point," Journal of Nuclear Materials, Chemical Society, Vol. 69, p. 2105, 1949.                          Vol. 29, pp. 356-358, 1969.
"Enthalpy of Solid Uranium Dioxide from 25000K  
to its Melting Point," Journal of Nuclear Materials, Vol. 29, pp. 356-358, 1969.


1.774
1.774 K.,


5. M. G. Chasanov, Argonne National Laboratory,              Oxide," Nuclear Science and Engineering, Vol. 8, Reactor Development Progress Report, ANL-7618,           pp. 497-506,1960.
5. M. G. Chasanov, Argonne National Laboratory, Reactor Development Progress Report, ANL-7618, September 1969.


September 1969.
6. E. Hellstrand et al., "The Temperature Coefficient of the Resonance Integral for Uranium Metal and Oxide," Nuclear Science and Engineering, Vol. 8, pp. 497-506,1960.


7. L. S. Tong, "Prediction of Departure from Nucleate Boiling for an Axially Non-Uniform Heat Flux
7. L. S. Tong, "Prediction of Departure from Nucleate Boiling for an Axially Non-Uniform Heat Flux Distribution," Journal of Nuclear Energy, Vol. 21, pp. 241-248, 1967.
6. E. Hellstrand et al., "The Temperature Coefficient        Distribution," Journal of Nuclear Energy, Vol. 21, of the Resonance Integral for Uranium Metal and          pp. 241-248, 1967.


1.77-5
1.77-5


APPENDIX B
APPENDIX B
                                          RADIOLOGICAL ASSUMPTIONS
RADIOLOGICAL ASSUMPTIONS
      The assumptions given below should be applied in               pressure defined in the technical specifications for determining a conservative source term and subsequent               containment leak testing.
The assumptions given below should be applied in determining a conservative source term and subsequent transport of activity and resulting doses to the public for use in evaluating the radiological consequences of a control rod ejection accident.
 
1. The assumptions related to the release of radioactive material to the primary containment are as follows:
a.
 
The case resulting in the largest source term should be selected for evaluation.
 
b. The nuclide inventory in the fuel elements potentially breached should be calculated, and it should be assumed that all gaseous constituents in the fuel-clad gaps are released.


transport of activity and resulting doses to the public for use in evaluating the radiological consequences of a                    I. Release of fission products to the secondary control rod ejection accident.                                      system should be computed by assuming that all fission products released from the fuel clad are uniformly
c.
1. The assumptions related to the release of radioactive            mixed in the primary coolant volume.


material to the primary containment are as follows:
The amount of activity accumulated in the fuel-clad gap should be assumed to be 10% of the iodines and 10% of the noble gases accumulated at the end of core life, assuming continuous maximum full power operation.
                                                                          j. The primary-to-secondary leak rate limitation a. The case resulting in the largest source term              incorporated or to be incorporated as a technical should be selected for evaluation.                                  specification requirement should be assumed to exist until the primary system pressure falls below the b. The nuclide inventory in the fuel elements                  secondary system pressure.


potentially breached should be calculated, and it should be assumed that all gaseous constituents in the fuel-clad                k. The release of fission products from the gaps are released.                                                  secondary system should be evaluated with the assump.
d.
 
No allowance should be given for activity decay prior to accident initiation, regardless of the reactor status for the selected case.
 
e. The nuclide inventory of the fraction of the fuel which reaches or exceeds the initiation temperature of fuel melting (typically 2842fC) at any time during the course of the accident should be calculated, and  
100% of the noble gases and 25% of the iodine contained in this fraction should be assumed to be available for release from the containment.
 
f. The effects of radiological decay during holdup in the containment or other buildings should be taken into account.
 
& The reduction in the amount of radioactive material available for leakage to the environment by containment sprays, recirculating filter systems, or other engineered safety features may be taken into account, but the amount of reduction in concentration of radioactive materials should be evaluated on a case-by case basis.
 
h. The primary reactor containment should be assumed to leak at the leak rate incorporated or to be incorporated as a technical specification requirement at peak accident pressure for the first 24 hours, and at 50%
of this leak rate for the remaining duration of the accident.'
Peak accident pressure is the maximum t The effect on containment leakage under accident conditions of features provided to reduce the leakage of radioactive materials from the containment should be evaluated on a case-by-case basis.
 
pressure defined in the technical specifications for containment leak testing.
 
I.
 
Release of fission products to the secondary system should be computed by assuming that all fission products released from the fuel clad are uniformly mixed in the primary coolant volume.
 
j.
 
The primary-to-secondary leak rate limitation incorporated or to be incorporated as a technical specification requirement should be assumed to exist until the primary system pressure falls below the secondary system pressure.
 
k. The release of fission products from the secondary system should be evaluated with the assump.


tion of a coincident loss of offsite power.
tion of a coincident loss of offsite power.


c. The amount of activity accumulated in the fuel-clad gap should be assumed to be 10% of the iodines and 10% of the noble gases accumulated at the end of                2. Acceptable assumptions for atmospheric diffusion core life, assuming continuous maximum full power                  and dose conversion are:
2.
operation.                                                              a. The 0-to-8-hour ground-level release concentra d. No allowance should be given for activity decay            tions may be reduced by a factor ranging from one to a prior to accident initiation, regardless of the reactor              maximum of three (see Figure 1) for additional dis status for the selected case.                                        persion produced by the turbulent wake of the reactor building in calculating potential exposures. The volu e. The nuclide inventory of the fraction of the                metric building wake correction, as defined in Section fuel which reaches or exceeds the initiation temperature            3-3.5.2 of Meteorology and Atomic Energy 1968 (Ref. 1),
 
of fuel melting (typically 2842fC) at any time during                should be used only in the 0-to-8-hour period; it is used the course of the accident should be calculated, and                with a shape factor of 1/2 and the minimum cross
Acceptable assumptions for atmospheric diffusion and dose conversion are:  
100% of the noble gases and 25% of the iodine                      sectional area of the reactor buiding only.
a.
 
The 0-to-8-hour ground-level release concentra tions may be reduced by a factor ranging from one to a maximum of three (see Figure 1) for additional dis persion produced by the turbulent wake of the reactor building in calculating potential exposures. The volu metric building wake correction, as defined in Section  
3-3.5.2 of Meteorology and Atomic Energy 1968 (Ref. 1),  
should be used only in the 0-to-8-hour period; it is used with a shape factor of 1/2 and the minimum cross sectional area of the reactor buiding only.
 
b.


contained in this fraction should be assumed to be available for release from the containment.                              b. No correction should be made for depletion of the effluent plume of radioactive iodine due to deposi f. The effects of radiological decay during holdup            tion on the ground or for the radiological decay of in the containment or other buildings should be taken                iodine in transit.
No correction should be made for depletion of the effluent plume of radioactive iodine due to deposi tion on the ground or for the radiological decay of iodine in transit.


into account.
c. For the first 8 hours, the breathing rate of a person offsite should be assumed to be 3.47 x 10-4 m3 /sec. From 8 to 24 hours following the accident, the breathing rate should be assumed to be 1.75 x 10-4 m3/sec. From 24 hours until the end of the accident, the rate should be assumed to be 2.32 x 10-
m3/sec.


c. For the first 8 hours, the breathing rate of a
(These values were developed from the average daily.
      & The reduction in the amount of radioactive                  person offsite should be assumed to be 3.47 x 10-4 material available for leakage to the environment by                m3 /sec. From 8 to 24 hours following the accident, the containment sprays, recirculating filter systems, or other          breathing rate should be assumed to be 1.75 x 10-4 engineered safety features may be taken into account,                m3 /sec. From 24 hours until the end of the accident, the but the amount of reduction in concentration of                      rate should be assumed to be 2.32 x 10- m3 /sec.


radioactive materials should be evaluated on a case-by              (These values were developed from the average daily.
breathing rate [2 x 107 cm3/day] assumed in a report (Ref. 2) of ICRP.2)
d.


case basis.                                                          breathing rate [22 x 107 cm3 /day] assumed in a report (Ref. 2) of ICRP. )
The iodine dose conversion factors are also given in Reference 2.
      h. The primary reactor containment should be assumed to leak at the leak rate incorporated or to be                    d. The iodine dose conversion factors are also given incorporated as a technical specification requirement at            in Reference 2.


peak accident pressure for the first 24 hours, and at 50%
e. External whole body doses should be calculated using "infinite cloud" assumptions, i.e., the dimensions of the cloud are assumed to be large compared to the distance that the gamma rays and beta particles travel.
of this leak rate for the remaining duration of the                      e. External whole body doses should be calculated accident.' Peak accident pressure is the maximum                    using "infinite cloud" assumptions, i.e., the dimensions of the cloud are assumed to be large compared to the t
                                                                    distance that the gamma rays and beta particles travel.


The effect on containment leakage under accident            "Such a cloud would be considered an infinite cloud for conditions of features provided to reduce the leakage of radioactive materials from the containment should be evaluated            2 on a case-by-case basis.                                                    International Commission on Radiological Protection.
"Such a cloud would be considered an infinite cloud for  
2 International Commission on Radiological Protection.


1.77-6
1.77-6


a receptor at the center because any additional [gamma                 g. The atmospheric diffusion model should be as and] beta emitting material beyond the cloud dimen               follows:
a receptor at the center because any additional [gamma and] beta emitting material beyond the cloud dimen sions would not alter the flux of [gamma rays and] beta particles to the receptor." (Ref. 3) Editorial additions were made to the quotation so that gamma as well as beta emitting material could be considered. Under these conditions, the rate of energy absorption per unit volume is equal to the rate of energy released per unit volume. For an infinite, uniform cloud containing X  
sions would not alter the flux of [gamma rays and] beta particles to the receptor." (Ref. 3) Editorial additions were made to the quotation so that gamma as well as beta                     (1) The basic equabion for atmospheric dif emitting material could be considered. Under these               fusion from a ground-level point source is:
curies of beta radioactivity per cubic meter, the beta dose in air at the cloud center is:  
conditions, the rate of energy absorption per unit volume is equal to the rate of energy released per unit                           XIQ        *OyOz volume. For an infinite, uniform cloud containing X               where curies of beta radioactivity per cubic meter, the beta dose in air at the cloud center is:
=D'
                                                                    X      =   the short-term average centerline value        of
= 0.457"E9 The surface body dose rate from beta emitters in the infinite cloud can be approximated as being one-half this amount (i.e., PD` = 0.23 Tp,. For gamma emitting material, the dose rate in air at the cloud center is:  
                          = 0.457"E9
S=  
                          =D'                                                    the ground-level concentration (Ci/m 3 )
From a semi-infinite cloud, the gamma dose rate in air is:
                                                                    Q      =    amount of material released (Ci/sec)
7D' = 0.25 EY
The surface body dose rate from beta emitters in the               u            windspeed (m/sec)
7x where PD' =
infinite cloud can be approximated as being one-half this           oy      =    the horizontal standard deviation of the amount (i.e., PD` = 0.23 Tp,. For gamma emitting                                 plume (meters), [see Figure V-1, Ref. 5].
beta dose rate from an infinite cloud (rad/sec)
material, the dose rate in air at the cloud center is:               z    =   the vertical standard deviation of the plume (meters) [see Figure V-2, Ref. 51.
,yD; =
gamma dose rate from an infinite cloud (rad/sec)
=
average beta energy per disintegration (Mev/dis)
=
average gamma energy per disintegration (Mev/dis)
X
=
concentration of beta or gamma emitting isotope in the cloud (Ci/m 3)  
f.
 
The following specific assumptions are accep table with respect to the radioactive cloud dose calcula tions:
(1) The dose at any distance from the reactor should be calculated based on the maximum concentra tion in the plume at that distance, taking into account special meteorological, topographical, and other char acteristics which may affect the maximum plume con centration. These site-related characteristics must be evaluated on a case-by-case basis. In the case of beta radiation, the receptor is assumed to be exposed to an infinite cloud at the maximum ground-level concentra tion at that distance from the reactor. In the case of gamma radiation, the receptor is assumed to be exposed to only one-half the cloud owing to the presence of the ground. The maximum cloud concentration always should be assumed to be at ground level.
 
(2) The appropriate average beta and gamma energies emitted per. disintegration, as given in the Table of Isotopes (Ref. 4), should be used.
 
g.


S=
The atmospheric diffusion model should be as follows:
                                                                            (2) For time periods greater than 8 hours, the From a semi-infinite cloud, the gamma dose rate in air          plume should be assumed to meander and spread is:                                                              uniformly over a 22.50 sector. The resultant equation is:
(1) The basic equabion for atmospheric dif fusion from a ground-level point source is:  
                    7D' = 0.25 EY  7x
XIQ
                                                                                        *dQ =2.032 GzUX
*OyOz where X
                                                                  where where PD' =      beta dose rate from an infinite cloud                  x    =    distance from the point of release to the (rad/sec)                                                       receptor; other variables are as given in
=  
      ,yD; =     gamma dose rate from an infinite cloud                          paragraph g. (1), above.
the short-term average centerline value of the ground-level concentration (Ci/m 3)
Q
=
amount of material released (Ci/sec)  
u windspeed (m/sec)
oy
=
the horizontal standard deviation of the plume (meters), [see Figure V-1, Ref. 5]. 
z
=  
the vertical standard deviation of the plume (meters) [see Figure V-2, Ref. 51.


(rad/sec)
(2) For time periods greater than 8 hours, the plume should be assumed to meander and spread uniformly over a 22.50 sector. The resultant equation is:
            =   average beta energy per disintegration (Mev/dis)                                                  (3) The atmospheric diffusion model3 for
* dQ =2.032 GzUX
            =   average gamma energy per disintegration          ground-level releases is based on the information in the (Mev/dis)                                       table below.
where x
=  
distance from the point of release to the receptor; other variables are as given in paragraph g. (1), above.


X    =    concentration of beta or gamma emitting isotope in the cloud (Ci/m 3 )
(3)  
                                                                    Time f. The following specific assumptions are accep            Following table with respect to the radioactive cloud dose calcula        Accident                    Atmospheric Conditions tions:
ground-level table below.
                                                                  0-8          Pasquill Type F, wind speed I m/sec, uniform
            (1) The dose at any distance from the reactor        hours        direction should be calculated based on the maximum concentra tion in the plume at that distance, taking into account          8-24        Pasquill Type F, wind speed I m/sec, special meteorological, topographical, and other char          hours        variable direction within a 22.50 sector acteristics which may affect the maximum plume con centration. These site-related characteristics must be          1-4          (a) 40% Pasquill Type D, wind speed 3 evaluated on a case-by-case basis. In the case of beta          days                m/sec radiation, the receptor is assumed to be exposed to an                        (b) 60% Pasquill Type F, wind speed 2 infinite cloud at the maximum ground-level concentra                                m/sec tion at that distance from the reactor. In the case of                        (c) wind direction - variable within a 22.50
  gamma radiation, the receptor is assumed to be exposed                              sector.


to only one-half the cloud owing to the presence of the ground. The maximum cloud concentration always
Time Following Accident The atmospheric diffusion model3 for releases is based on the information in the Atmospheric Conditions
                                                                        3 should be assumed to be at ground level.                                This model should be used until adequate site meteorological data are obtained. In some cases, available
0-8 Pasquill Type F, wind speed I m/sec, uniform hours direction
            (2) The appropriate average beta and gamma          Information, such as meteorology, topography, and geographical energies emitted per. disintegration, as given in the Table    location, may dictate the use of a more restrictive model to of Isotopes (Ref. 4), should be used.                          insure a conservative estimate of potential offsite exposures.
8-24 hours
1-4 days Pasquill Type F, wind speed I m/sec, variable direction within a 22.50 sector (a) 40% Pasquill Type D, wind speed 3 m/sec (b) 60% Pasquill Type F, wind speed 2 m/sec (c) wind direction - variable within a 22.50
sector.
 
3 This model should be used until adequate site meteorological data are obtained. In some cases, available Information, such as meteorology, topography, and geographical location, may dictate the use of a more restrictive model to insure a conservative estimate of potential offsite exposures.


1.77-7
1.77-7


Time.                                                         Time Following                                                    Following Accident                   Atmospheric Conditions           Accident               Atomospheric Conditions
Time.
                                                            4-30       (c) 33.3% Pasquill Type F, wind speed 2 days           m/sec (d) Wind direction - 33.3% frequency in a
 
4-30        (a) 33.3% Pasquill Type C, wind speed 3                          22.50 sector.
Following Accident  
4-30
(a)
days (b)
Atmospheric Conditions  
33.3% Pasquill Type C, wind speed 3 m/sec
33.3% Pasquill Type D, wind speed 3 m/sec Time Following Accident Atomospheric Conditions
4-30  
(c) 33.3% Pasquill Type F, wind speed 2 days m/sec (d) Wind direction - 33.3% frequency in a  
22.50 sector.


days            m/sec                                                (4) Figures 2(A) and 2(B) give the ground-level (b) 33.3% Pasquill Type D, wind speed 3        release atmospheric diffusion factors based on the m/sec                                      parameters given in paragraph g&(3), above.
(4) Figures 2(A) and 2(B) give the ground-level release atmospheric diffusion factors based on the parameters given in paragraph g&(3), above.


REFERENCES
REFERENCES
1. D. H. Slade, ed., "Meteorology and Atomic Energy         4. C. M. Lederer, J. M. Hollander, and I. Perlman,
1. D. H. Slade, ed., "Meteorology and Atomic Energy  
    -  1968," TID-24190, Division of Technical                "Table of Isotopes," Sixth Edition, University of Information, U.S. Atomic Energy Commission, July            California, Berkeley, Lawrence Radiation Lab
-
    1968.                                                      oratory.
1968,"
TID-24190,
Division of Technical Information, U.S. Atomic Energy Commission, July
1968.
 
2. Report of Committee I1, "Permissible Dose for Internal Radiation," ICRP Publication 2, 1959.
 
3. D. H. Slade, op. cir., Section 7.4.1.1.
 
4. C. M. Lederer, J. M. Hollander, and I. Perlman,  
"Table of Isotopes," Sixth Edition, University of California, Berkeley, Lawrence Radiation Lab oratory.


2. Report of Committee I1, "Permissible Dose for Internal Radiation," ICRP Publication 2, 1959.            5. F. A. Gifford, Jr., "Use of Routine Meteorological Observations for Estimating Atmospheric Dis
5. F. A. Gifford, Jr., "Use of Routine Meteorological Observations for Estimating Atmospheric Dis
3. D. H. Slade, op. cir., Section 7.4.1.1.                      .persion," Nuclear Safety, Vol. 2, No. 4, June 1961.
.persion," Nuclear Safety, Vol. 2, No. 4, June 1961.


1.77-8}}
1.77-8}}


{{RG-Nav}}
{{RG-Nav}}

Latest revision as of 02:06, 17 January 2025

Assumptions Used for Evaluating a Control Rod Ejection Accident for Pressurized Water Reactors
ML003740279
Person / Time
Issue date: 05/31/1974
From:
Office of Nuclear Regulatory Research
To:
References
RG-1.77
Download: ML003740279 (8)
v • d • e


U.S. ATOMIC, ENERGY COMMISSION

REGULATORY

DIRECTORATE OF REGULATORY STANDARDS

May 1974 GUIDE

REGULATORY GUIDE 1.77 ASSUMPTIONS USED FOR EVALUATING A CONTROL ROD

EJECTION ACCIDENT FOR PRESSURIZED WATER REACTORS

A. INTRODUCTION

Section 50.34, 'Contents of applications: technical information," of 10 CFR Part 50, "Licensing of Pro duction and Utilization Facilities," requires that each application for a construction permit or operating license provide an analysis and evaluation of the design and performance of structures, systems, and components of the facility with the objective of assessing the potential risk to public health and safety resulting from operation of the facility. General Design Criterion 28,

"Reactivity Limits," of Appendix A, "General Design Criteria for Nuclear Power Plants," to 10 CFR Part 50,

requires the reactivity control system to be designed with appropriate limits on the potential amount and rate of reactivity increase to assure that the effects of postulated reactivity accidents can neither result in damage to the reactor coolant pressure boundary greater than limited local yielding nor sufficiently disturb the core, its support structures, or other reactor pressure vessel internals to impair significantly the capability to cool the core. General Design Criterion 28 also requires that these postulated reactivity accidents include con sideration of the rod ejection accident unless such an accident is prevented by positive means.

This guide identifies acceptable analytical methods and assumptions that may be used in evaluating the consequences of a rod ejection accident in uranium oxide-fueled pressurized water reactors (PWRs). In some cases, unusual site characteristics, plant design features, or other factors may require different assumptions which will be considered on an individual basis. The Advisory Committee on Reactor Safeguards has been consulted concerning this guide and has concurred in the regulatory position.

B. DISCUSSION

The rate at which reactivity can be inserted into the core of a uranium oxide-fueled water-cooled power reactor is normally limited by the-design of the control rod system to a value well below that which would result in serious damage to the reactor system. However, a postulated failure of the control rod system provides the potential for a relatively high rate of reactivity insertion which, if large enough, could cause a prompt power burst. For U0 2 fuel, a large fraction of this generated nuclear energy is stored momentarily in the fuel and then released to the rest of the system. If the fuel energy densities were high enough, there would exist the potential for prompt rupture of fuel pins and the consequent rapid heat transfer to the water from finely dispersed molten U0 2. Prompt fuel element rupture is defined herein as a rapid increase in internal fuel rod pressure due to extensive fuel melting, followed by rapid fragmentation and dispersal of fuel cladding into the coolant. This is accompanied by the conversion of nuclear energy, deposited as overpower heat in the fuel and in the coolant, to mechanical energy which, in sufficient quantity, could conceivably disarrange the reactor core or breach the primary system.

The Regulatory staff has reviewed the available experimental information concerning fuel failure thresh olds. In general, failure consequences for U0 2 have been insignificant below 300 cal/g for both irradiated and unirradiated fuel rods. Therefore, a calculated radial average energy density of 280 cal/g at any axial fuel location in any fuel rod as a result of a postulated rod ejection accident provides a conservative maximum limit to ensure that core damage will be minimal and that both short-term andlong-term core cooling capability will not be impaired.

For the postulated control rod ejection accident, a mechanical failure of a control rod mechanism housing is assumed such that the reactor coolant system pressure would eject the control rod and drive shaft to the fully withdrawn position.

USAEC REGULATORY GUIDES

Copies of published guides may be obtained by request Indicating the divisions desired to the US. Atomic Energy Commission. Washington, D.C. 20545.

Regulatory Guides we issued to describe and' make available to the public Attention: Director of Regulatory Standards. Comments and suggestions for methods acceptale to the AEC Regulatory staff of implemnenting specfifc parts of Improvements In thes guldes are encouraged and should be sent to the Secretary the Commission's regulations, to delineate techniques used -by the staff In of the Commission. US. Atomic Energy Commission, Washington, D.C. 20645 evaluating specific problems or postulated accidents, or to provide guidance to Attention. Chief, Public Proceedings Staff.

applicants. 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 era issued in the following ten broad divisions:

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.

1. Power Rleactors S. Products

2. Research and Test Reactors

7. Transportation

3. Fuels and Materials Facilities

8. Occupational Health Published guides will be revised periodically, as appropriate, to accommodate

4. Environmental and Siting

9. Antitrust Review comments and to rflect new Information or experience.

L. Materials Ond Plant Protection

1

0. General

A sufficient number of initial reactor states to completely bracket all possible operational conditions of interest should be analyzed to assure examination of upper bounds on ultimate damage. In areas of uncer tainty, the appropriate minimum or maximum para meters relative to nominal or expected values should be used to assure a conservative evaluation. The initial reactor states should include consideration of at least the following:

Zero power (hot standby) - Beginning of Life (BOL)

and End of Life (EOL);

Low power

- BOL and EOL;

Ful

l. power

- BOL and EOL.

The effects of the loss of primary system integrity as a result of the failed control rod housing should be included in the analysis. It should- also be shown that failure of one control rod housing will not lead to failure of other control rod housings.

The approach that should be used in the radiological analysis of a control rod ejection accident is to deter mine the amount of each gaseous radionuclide released to the primary containment and, with this information in conjunction with the procedures set forth in Appendix B of this guide, to determine the radiological consequences of this accident for a pressurized water reactor.

C. REGULATORY POSITION

Acceptable assumptions and evaluation models for analyzing a rod. ejection accident in PWRs are presented in Appendices A (Physics and Thermal-Hydraulics) and B (Radiological Assumptions) of this guide. By use of these appendices, it should be shown that:

1. Reactivity excursions will not result in a radial average fuel enthalpy greater than 280 cal/g at any axial location in any fuel rod.

2. Maximum reactor pressure during any portion of the assumed transient will be less than the value that will cause-stresses to-exceed the Emergency Condition-stress limits as defined in Section III of the ASME Boiler and Pressure Vessel Code.'

3.

Offsite dose consequences will be well within the guidelines of 10 CFR Part 100, "Reactor Site Criteria."

1Copies may be -obtained from the American Society of Mechanical Engineers, United Engineering Center, 345 East 47th Street, New York, New York 10017.

-

1.77-2

APPENDIX A

PHYSICS AND THERMAL-HYDRAULICS

The assumptions described below should be applied in evaluating the physics and thermal-hydraulic behavior of the reactor system for a control rod ejection accident.

I.

The ejected rod worth should be calculated based on the maximum worth rod resulting from the following conditions: (a) all control banks at positions corres ponding to values for maximum allowable bank inser tions at a given power level and (b) additional fully or partially inserted misaligned or inoperable rod or rods if allowed by operating procedures. Sufficient parametric studies should be performed to determine the worth of the most reactive control rod in each rod group for different control rod configurations, both expected and unexpected. The worth of single rods in rod groups should be evaluated during startup physics tests and compared with values used in the rod ejection analysis.

The accident should be reanalyzed if the rod worths used in the initial analysis are found to be noncon servative. Calculated rod worths should be increased, if necessary, to account for calculational uncertainties in parameters such as neutron cross sections and power asymmetries due to xenon oscillations.

2. The reactivity insertion rate due to an ejected rod should be determined from differential control rod worth curves and calculated transient rod position versus time curves. If differential rod worth curves are not available for the reactor state of interest, conservatism should be included in the calculation of reactivity insertion through consideration of the nonlinearity in reactivity addition as the rod passes through the active core. The rate of ejection should be calculated based on the maximum pressure differential and the weight and'

cross-sectional area of the control rod and drive shaft, assuming no pressure barrier restriction.

3.

The calculation of effective delayed neutron frac tion (Beff) and prompt neutron lifetime (2*) should be based on the well-known definitions resulting from per turbation theory, such as those described by Henry (Ref.

1), using available experimental delayed neutron data and averaging by the fraction of fission in the various fissionable materials. In cases where the accident is quite sensitive to Peff (where the ejected rod worth >Peff), the minimum calculated value for the given reactor state should be used. For smaller transients, conservatism in the value should include consideration of not only the initial power rise (which increases with decreasing P), but also the power reduction after the trip. Similar don siderations should also be applied to determine an appropriately conservative value of 2* to be used.

4. The initial reactor coolant 1 pressure, core inlet temperature, and flow rate used in the analysis should be conservatively chosen with respect to their influence on the magnitude of the transient. Pressure and temperature are mainly significant with respect to their effect on the amount of reactivity inserted if there exists a positive moderator coefficient.

5. ' The fuel thermal properties such'"as fuel-clad gap heat transfer coefficient and fuel thermal conductivity

.should be conservatively chosen, depending upon the transient phenomenon being investigated. For conditions of a zero or positive moderator coefficient (usually at beginning of life), for example, high heat transfer parameters would reduce the Doppler feedback and increase any positive moderator feedback effects and hence tend to increase the magnitude of the reactivity transient. For a negative moderator coefficient, high heat transfer parameters could cause the magnitude of the transient to decrease if a given quantity of heat produces more feedback in the moderator than in the fuel. In the consideration of pressure pulses which may be generated, high moderator heating rates could cause significant pressure gradients to develop in the moder ator channels. In computing the average enthalpy of the hottest fuel pellet during the excursion for power cases, low heat transfer would be conservative.

6. The specific heat of U0 2 has been determined experimentally and is a deterministic factor in the calculated amount of stored energy (enthalpy) in the fuel. Recommended values in the range of 25 to 902fC

are the data reported by Moore and Kelly (Ref. 2). In the range of 900 to 2842!C, the data obtained by Hein and Flagella (Ref. 3), Leibowitz, Mishler, and Chasanov (Ref. 4), and Chasanov (Ref. 5) are recommended for the heat capacity of the fuel. These recommended values are for clean core conditions. Possible variation in the specific heat due to bumup should be investigated and appropriate values used, if necessary.

7.

The moderator reactivity coefficients due to voids, coolant pressure changes, and coolant temperature changes should be calculated based on the various assumed conditions of the fuel and moderator using standard transport and diffusion theory codes. If no three-dimensional space-time kinetics calculation is per formed, the reactivity feedback due to these coefficients should be conservatively weighted to account for the variation in their spatial Importance in the missing

'dimension(s). If boric acid shim is used in the moder ator, the highest boron concentration corresponding to the initial reactor state should be assumed.

8.

The Doppler coefficient should be calculated based on the effective resonance integrals and should include corrections for pin shadowing (Dancoff correction).

Calculations of the Doppler coefficient of reactivity should be based on and should compare conservatively

1.77-3

with available experimental data such as those of Hell.

strand (Ref. 6). Since the Doppler coefficient reflects the change in reactivity as a function of fuel temperature, uncertainties in predicting fuel temperatures at different power levels should be reflected by conservatism in the applied value of the Doppler coefficient. If no three dimensional space-time kinetics calculation is performed, the reactivity effect of spatially weighting the core average temperature rise in both the axial and radial directions should be calculated.

9. Control rod reactivity insertion during trip versus time should be obtained by combining the differential rod worth curve with a rod velocity curve based on maximum design limit values for scram insertion times.

If the rod worth curve (reactivity vs. depth of insertion)

is not obtained from a "true" representation (i.e., an x, y, z, t or an r, z, t calculation); the conservatism of the approximate calculation should be shown. The dif ference in the depth of insertion at zero power and at full power should be accounted for in calculating the available scram reactivity.

10. The reactor trip delay time, or the amount of time which elapses between the instant the sensed parameter (e.g., pressure or neutron flux) reaches the level for which protective action is required and the onset of negative reactivity insertion, should be based on maxi mum values of the following: (a) time required for instrument channel to produce a signal, (b) time for the trip breaker to open, (c) time for the coil to release the rods, and (d) time required before scram rods enter the core if the tips lie above the core-reflector interface.

11. The computer code used for calculating the tran sient should be a coupled thermal, hydrodynamic, and nuclear model with the following capabilities: (a) incor poration of all major reactivity feedback mechanis ms, (b) at least six delayed neutron groups, (c) both axial and radial segmentation of the fuel element, (d) coolant flow provision, and (e) control rod scram initiation on either coolant system pressure or neutron flux.

12. The analytical models and computer codes used should be documented and justified and the con servatism of the models and codes should be evaluated both by comparison with experiment, as available, and with more sophisticated spatial kinetics codes. In par ticular, the importance of two- or three-dimensional flux characteristics and changes in flux shapes should be investigated, and the conservatism of the flux shapes used for reactivity input and feedback, peak energy deposition, total energy, and gross heat transfer to the coolant should be evaluated. Also, sensitivity studies on variations of the Doppler effect, power distribution, fuel element heat transfer parameters, and other relevant parameters should be included.

13. The pressure surge should be calculated on the basis of conventional heat transfer from the fuel,. a conserva tive metal-water reaction threshold, and prompt heat generation in the coolant to determine the variation of heat flux with time and the volume surge. The volume surge should then be used in the calculation of the pressure transient, taking into account fluid transport in the system, heat transfer to the steam generators, and the action of the pressurizer relief and.safety valves. No credit should be taken for the possible pressure reduc tion caused by the assumed failure of the control rod pressure housing.

14. The number of fuel rods experiencing clad failure.

should be calculated and used to obtain the amount of contained fission product inventory released to the reactor coolant system. It should be assumed that clad failure occurs if the heat flux equals or exceeds the value corresponding to the onset of the transition from nucleate to film boiling (DNB), or for other appropriate causes.

The margin to DNB is expressed in terms of a departure from nucleate boiling ratio (DNBR). The DNBR at any position in the hottest channel is the ratio of the DNB heat flux to the actual heat flu

x. The DNB

heat flux should be evaluated using correlations based on recognized studies and experimental heat transfer DNB

data. A minimum DNBR should be determined from the evaluation of the experimental data to ensure a 95%

probability with a 95% confidence level that DNB has not occurred for the fuel element being evaluated. One example of a correlation which has been used to date is given by Tong (Ref. 7). The use of this correlation'and the above probabilities and confidence level yields a mini mum DNBR of 1.30. Other DNB or clad failure correlations may be used if they are adequately justified by analytical methods and supported by sufficient experimental data.

REFERENCES

I. A. F. Henry, "Computation of Parameters Appear ing in the Reactor Kinetics Equations," WAPD-142, December 1955.

2. G. E. Moore and K. K. Kelley, "High Temperature Heat Contents of Uranium, Uranium Dioxide and Uranium Trioxide," Journal of the American Chemical Society, Vol. 69, p. 2105, 1949.

3. R.

A.

Hein and P. N. Flagella, "Enthalpy Measurements of U0 2 and Tungsten to 32600K,"

GEMP-578, February 1968.

4. L Leibowitz, L. W. Mishler, and M. G. Chasanov,

"Enthalpy of Solid Uranium Dioxide from 25000K

to its Melting Point," Journal of Nuclear Materials, Vol. 29, pp. 356-358, 1969.

1.774 K.,

5. M. G. Chasanov, Argonne National Laboratory, Reactor Development Progress Report, ANL-7618, September 1969.

6. E. Hellstrand et al., "The Temperature Coefficient of the Resonance Integral for Uranium Metal and Oxide," Nuclear Science and Engineering, Vol. 8, pp. 497-506,1960.

7. L. S. Tong, "Prediction of Departure from Nucleate Boiling for an Axially Non-Uniform Heat Flux Distribution," Journal of Nuclear Energy, Vol. 21, pp. 241-248, 1967.

1.77-5

APPENDIX B

RADIOLOGICAL ASSUMPTIONS

The assumptions given below should be applied in determining a conservative source term and subsequent transport of activity and resulting doses to the public for use in evaluating the radiological consequences of a control rod ejection accident.

1. The assumptions related to the release of radioactive material to the primary containment are as follows:

a.

The case resulting in the largest source term should be selected for evaluation.

b. The nuclide inventory in the fuel elements potentially breached should be calculated, and it should be assumed that all gaseous constituents in the fuel-clad gaps are released.

c.

The amount of activity accumulated in the fuel-clad gap should be assumed to be 10% of the iodines and 10% of the noble gases accumulated at the end of core life, assuming continuous maximum full power operation.

d.

No allowance should be given for activity decay prior to accident initiation, regardless of the reactor status for the selected case.

e. The nuclide inventory of the fraction of the fuel which reaches or exceeds the initiation temperature of fuel melting (typically 2842fC) at any time during the course of the accident should be calculated, and

100% of the noble gases and 25% of the iodine contained in this fraction should be assumed to be available for release from the containment.

f. The effects of radiological decay during holdup in the containment or other buildings should be taken into account.

& The reduction in the amount of radioactive material available for leakage to the environment by containment sprays, recirculating filter systems, or other engineered safety features may be taken into account, but the amount of reduction in concentration of radioactive materials should be evaluated on a case-by case basis.

h. The primary reactor containment should be assumed to leak at the leak rate incorporated or to be incorporated as a technical specification requirement at peak accident pressure for the first 24 hours2.777778e-4 days <br />0.00667 hours <br />3.968254e-5 weeks <br />9.132e-6 months <br />, and at 50%

of this leak rate for the remaining duration of the accident.'

Peak accident pressure is the maximum t The effect on containment leakage under accident conditions of features provided to reduce the leakage of radioactive materials from the containment should be evaluated on a case-by-case basis.

pressure defined in the technical specifications for containment leak testing.

I.

Release of fission products to the secondary system should be computed by assuming that all fission products released from the fuel clad are uniformly mixed in the primary coolant volume.

j.

The primary-to-secondary leak rate limitation incorporated or to be incorporated as a technical specification requirement should be assumed to exist until the primary system pressure falls below the secondary system pressure.

k. The release of fission products from the secondary system should be evaluated with the assump.

tion of a coincident loss of offsite power.

2.

Acceptable assumptions for atmospheric diffusion and dose conversion are:

a.

The 0-to-8-hour ground-level release concentra tions may be reduced by a factor ranging from one to a maximum of three (see Figure 1) for additional dis persion produced by the turbulent wake of the reactor building in calculating potential exposures. The volu metric building wake correction, as defined in Section

3-3.5.2 of Meteorology and Atomic Energy 1968 (Ref. 1),

should be used only in the 0-to-8-hour period; it is used with a shape factor of 1/2 and the minimum cross sectional area of the reactor buiding only.

b.

No correction should be made for depletion of the effluent plume of radioactive iodine due to deposi tion on the ground or for the radiological decay of iodine in transit.

c. For the first 8 hours9.259259e-5 days <br />0.00222 hours <br />1.322751e-5 weeks <br />3.044e-6 months <br />, the breathing rate of a person offsite should be assumed to be 3.47 x 10-4 m3 /sec. From 8 to 24 hours2.777778e-4 days <br />0.00667 hours <br />3.968254e-5 weeks <br />9.132e-6 months <br /> following the accident, the breathing rate should be assumed to be 1.75 x 10-4 m3/sec. From 24 hours2.777778e-4 days <br />0.00667 hours <br />3.968254e-5 weeks <br />9.132e-6 months <br /> until the end of the accident, the rate should be assumed to be 2.32 x 10-

m3/sec.

(These values were developed from the average daily.

breathing rate [2 x 107 cm3/day] assumed in a report (Ref. 2) of ICRP.2)

d.

The iodine dose conversion factors are also given in Reference 2.

e. External whole body doses should be calculated using "infinite cloud" assumptions, i.e., the dimensions of the cloud are assumed to be large compared to the distance that the gamma rays and beta particles travel.

"Such a cloud would be considered an infinite cloud for

2 International Commission on Radiological Protection.

1.77-6

a receptor at the center because any additional [gamma and] beta emitting material beyond the cloud dimen sions would not alter the flux of [gamma rays and] beta particles to the receptor." (Ref. 3) Editorial additions were made to the quotation so that gamma as well as beta emitting material could be considered. Under these conditions, the rate of energy absorption per unit volume is equal to the rate of energy released per unit volume. For an infinite, uniform cloud containing X

curies of beta radioactivity per cubic meter, the beta dose in air at the cloud center is:

=D'

= 0.457"E9 The surface body dose rate from beta emitters in the infinite cloud can be approximated as being one-half this amount (i.e., PD` = 0.23 Tp,. For gamma emitting material, the dose rate in air at the cloud center is:

S=

From a semi-infinite cloud, the gamma dose rate in air is:

7D' = 0.25 EY

7x where PD' =

beta dose rate from an infinite cloud (rad/sec)

,yD; =

gamma dose rate from an infinite cloud (rad/sec)

=

average beta energy per disintegration (Mev/dis)

=

average gamma energy per disintegration (Mev/dis)

X

=

concentration of beta or gamma emitting isotope in the cloud (Ci/m 3)

f.

The following specific assumptions are accep table with respect to the radioactive cloud dose calcula tions:

(1) The dose at any distance from the reactor should be calculated based on the maximum concentra tion in the plume at that distance, taking into account special meteorological, topographical, and other char acteristics which may affect the maximum plume con centration. These site-related characteristics must be evaluated on a case-by-case basis. In the case of beta radiation, the receptor is assumed to be exposed to an infinite cloud at the maximum ground-level concentra tion at that distance from the reactor. In the case of gamma radiation, the receptor is assumed to be exposed to only one-half the cloud owing to the presence of the ground. The maximum cloud concentration always should be assumed to be at ground level.

(2) The appropriate average beta and gamma energies emitted per. disintegration, as given in the Table of Isotopes (Ref. 4), should be used.

g.

The atmospheric diffusion model should be as follows:

(1) The basic equabion for atmospheric dif fusion from a ground-level point source is:

XIQ

  • OyOz where X

=

the short-term average centerline value of the ground-level concentration (Ci/m 3)

Q

=

amount of material released (Ci/sec)

u windspeed (m/sec)

oy

=

the horizontal standard deviation of the plume (meters), [see Figure V-1, Ref. 5].

z

=

the vertical standard deviation of the plume (meters) [see Figure V-2, Ref. 51.

(2) For time periods greater than 8 hours9.259259e-5 days <br />0.00222 hours <br />1.322751e-5 weeks <br />3.044e-6 months <br />, the plume should be assumed to meander and spread uniformly over a 22.50 sector. The resultant equation is:

  • dQ =2.032 GzUX

where x

=

distance from the point of release to the receptor; other variables are as given in paragraph g. (1), above.

(3)

ground-level table below.

Time Following Accident The atmospheric diffusion model3 for releases is based on the information in the Atmospheric Conditions

0-8 Pasquill Type F, wind speed I m/sec, uniform hours direction

8-24 hours

1-4 days Pasquill Type F, wind speed I m/sec, variable direction within a 22.50 sector (a) 40% Pasquill Type D, wind speed 3 m/sec (b) 60% Pasquill Type F, wind speed 2 m/sec (c) wind direction - variable within a 22.50

sector.

3 This model should be used until adequate site meteorological data are obtained. In some cases, available Information, such as meteorology, topography, and geographical location, may dictate the use of a more restrictive model to insure a conservative estimate of potential offsite exposures.

1.77-7

Time.

Following Accident

4-30

(a)

days (b)

Atmospheric Conditions

33.3% Pasquill Type C, wind speed 3 m/sec

33.3% Pasquill Type D, wind speed 3 m/sec Time Following Accident Atomospheric Conditions

4-30

(c) 33.3% Pasquill Type F, wind speed 2 days m/sec (d) Wind direction - 33.3% frequency in a

22.50 sector.

(4) Figures 2(A) and 2(B) give the ground-level release atmospheric diffusion factors based on the parameters given in paragraph g&(3), above.

REFERENCES

1. D. H. Slade, ed., "Meteorology and Atomic Energy

-

1968,"

TID-24190,

Division of Technical Information, U.S. Atomic Energy Commission, July

1968.

2. Report of Committee I1, "Permissible Dose for Internal Radiation," ICRP Publication 2, 1959.

3. D. H. Slade, op. cir., Section 7.4.1.1.

4. C. M. Lederer, J. M. Hollander, and I. Perlman,

"Table of Isotopes," Sixth Edition, University of California, Berkeley, Lawrence Radiation Lab oratory.

5. F. A. Gifford, Jr., "Use of Routine Meteorological Observations for Estimating Atmospheric Dis

.persion," Nuclear Safety, Vol. 2, No. 4, June 1961.

1.77-8


Retrieved from "https://kanterella.com/w/index.php?title=Regulatory_Guide_1.77&oldid=5446666"