ML16306A065
| ML16306A065 | |
| Person / Time | |
|---|---|
| Site: | Rhode Island Atomic Energy Commission |
| Issue date: | 11/01/2016 |
| From: | State of RI, Atomic Energy Comm, Nuclear Science Ctr |
| To: | Office of Nuclear Reactor Regulation |
| Shared Package | |
| ML16306A063 | List: |
| References | |
| Download: ML16306A065 (54) | |
Text
1 Fuel Failure Analysis 161031 Fuel Failure Mechanisms
- 1.
There are two basic types of fuel failure accidents that are relevant to the RINSC reactor:
A.
Mechanistic failures, in which the cladding is breached as a result of mechanical damage, and B.
Non-Mechanistic failures, in which the cladding is breached as a result of overheating.
- 2.
Of the two types of failure, the mechanistic failure is considered to be far more likely than the non-mechanistic failure at the RINSC facility. The most realistic scenario that could lead to a mechanistic fuel failure is fuel handling, which occurs on an annual basis during fuel inspections. Non-mechanistic failure is considered to be extremely unlikely. The transient analyses for the RINSC core show that over power excursions in which the over power protections operate as designed, lead to a peak fuel temperature of 88.8 C1. The loss of coolant accident analysis shows that for the maximum credible accident scenario, the peak cladding temperature for the most limiting fuel plate is still below the blister temperature of the cladding, given the predicted time that it would take for the pool level to drain down to the top of the fuel2.
- 3.
For a major non-mechanistic power plant fuel failure accident, there are five phases of failure3:
A.
Coolant Activity Release For power plants, initially there is usually a coolant system failure in which radionuclides in the coolant leak out of the system. If the cooling system is unable to keep the cladding sufficiently cool, cladding failure will occur.
B.
Gap Activity Release Inside the fuel, a gap exists between the fuel and the cladding. Radionuclides accumulate in this gap. When the cladding is breached, these nuclides are free to escape. During this phase, a few percent of the more volatile nuclides, such as Cs and I escape, while the bulk of the fission products are retained in the fuel.
2 C.
Early In-Vessel Release During this phase, the fuel begins to degrade and melt until it slumps to the bottom of the reactor vessel. During this phase, all of the noble gases, and a significant fraction of the volatile nuclides such as iodine and cesium are released.
D.
Ex-Vessel Release For a power plant, if the bottom of the reactor vessel fails, the molten core debris will interact with the concrete structural material beneath the vessel.
This will cause some of the less volatile nuclides to be release directly into the containment.
E.
Late In-Vessel Release At the same time that the ex-vessel release phase is occurring, and for a longer duration than that phase, late in-vessel release occurs. During this phase, volatile nuclides that had deposited in the cooling system during the in-vessel phase, are release into containment.
- 4.
Of the fuel failure phases that occur in nuclear power plants, only the gap and the onset of in-vessel release phases are relevant to the RINSC reactor. The RINSC reactor does not have an appreciable amount of activity buildup in the coolant due to activation, so coolant activity release is not a major concern. Gap activity release is relevant for any fuel cladding failure, regardless of whether it is due to mechanistic or non-mechanistic events, so this release consequence must be considered. It is plausible, though unlikely that the onset of early in-vessel release could occur, and that the integrity of the cladding could begin to be compromised during a fuel failure at the RINSC facility. However, the transient and loss of coolant accident analyses show that it is not plausible for this phase to progress to the point in which the core material slumps to the bottom of the pool. Ex vessel release, and late in vessel release are not considered to be credible for the RINSC reactor.
Radionuclides Released from Fuel to Coolant
- 1.
The available regulatory guidance regarding the quantity of radioactive material that gets into the coolant from the damaged fuel is geared toward nuclear power plants, rather than low power research reactors. The guidance is divided into information that pertains to low pressure boiling water reactors (BWR), versus information related to high pressure, pressurized reactors (PWR). Of the two types of reactors, the RINSC reactor is more similar to a BWR than a PWR. Consequently, the guidance for BWRs was used.
3
- 2.
The radionuclides that are expected to be released as a result of gap and early in-vessel fuel failure have been grouped on the basis of chemical similarity4. There are seven groups that are expected to get into the coolant5:
Group Elements Noble Gases Xe, Kr Halogens I, Br Alkali Metals Cs, Rb Tellurium Group Te, Sb, Se, Ba, Sr Noble Metals Ru, Rh, Pd, Mo, Tc, Co Lanthanides La, Zr, Nd, Eu, Nb, Pm, Pr, Sm, Y, Cm, Am Cerium Group Ce, Pu, Np Release Fraction from Fuel to Coolant
- 1.
Regulatory Guide 1.183 indicates that during the gap and early in-vessel release phases of fuel failure, the release fractions from the fuel to the confinement air are6:
Group Gap Release Phase Early In-Vessel Release Phase Total Noble Gases 0.05 0.95 1.0 Halogens 0.05 0.25 0.3 Alkali Metals 0.05 0.20 0.25 Tellurium Metals 0.00 0.05 0.05 Ba, Sr 0.00 0.02 0.02 Noble Metals 0.00 0.0025 0.0025 Lanthanides 0.00 0.0002 0.0002 Cerium Group 0.00 0.0005 0.0005
- 2.
This assumes an accident scenario in which there is not an appreciable amount of coolant through which these isotopes would have to pass in order to reach the confinement air.
- 3.
For the RINSC fuel failure scenario, it is anticipated that during these release phases there will still be an appreciable amount of coolant in the pool, so these release fractions can be considered to be the fractions that are released into the coolant.
Release Fraction from Coolant to Confinement Air
- 1.
As mentioned earlier, the credible fuel failure accident scenarios do not involve an instantaneous coolant loss, so the fuel failure scenario described in the guidance that best fits the facility is the fuel handling accident scenario.
4
- 2.
All of the fission products that are released from the fuel are in particulate form except for7:
A.
Elemental Iodine (I2)
B.
Organic Iodide C.
Noble Gases
- 3.
Particulates:
A.
The particulates are assumed to be retained by the water in the pool11.
Consequently, only the iodines and noble gases must be considered in the analysis.
- 4.
Iodine:
A.
If there were not an appreciable amount of coolant acting as a barrier between the fuel and the confinement air:
- 1.
The chemical form of the iodine that is released from the coolant into the confinement air would be7:
A.
Cesium Iodide (CsI) 95%
B.
Elemental Iodine (I2)4.85%
C.
Organic Iodide 0.15%
- 2.
The cesium iodide is expected to completely dissociate in the pool water, after which, the iodine is expected to re-evolve as elemental iodine (I2) due to the low pH of the water relatively instantaneously8.
B.
However, RINSC Technical Specifications require that the height of the pool water over the top of the fuel meat be at least 23 feet 7 inches for both natural convection operation, and forced convection operation9.
C.
If the depth of the water above the damaged fuel is 23 feet or greater, then 99.5% of the total iodine in the coolant remains in there, and the composition of the iodine that reaches the confinement air is10:
- 1.
Elemental Iodine (I2) 57%
- 2.
Organic Iodide 43%
- 5.
Noble Gases:
A.
The retention of the noble gases in the pool water is assumed to be negligible11.
5
- 6.
Therefore, the isotopes that are relevant due to the fact that a significant fraction of them will be transported to the confinement air are the iodines and noble gases.
A.
Given the fact that at the RINSC reactor the water depth is at least 23 feet over the fuel meat, the fraction of the halogens that are released into the pool water that will reach the confinement air is 0.5% of the content in the coolant.
B.
The noble gases will be unaffected by the coolant, so it is anticipated that 100% of the noble gases that are released into the pool water will reach the confinement air.
Overall Release Fraction from Fuel to Confinement Air
- 1.
Of the activity available to be released from the fuel, the overall release fraction to the confinement air is:
(Release Fraction from Fuel to Coolant)(Release Fraction from Coolant to Air)
- 2.
This means that the release fractions from the fuel to the confinement air are:
Group Gap Release Phase Early In-Vessel Release Phase Total Noble Gases 0.05 0.95 1.0 Halogens (0.05)(0.005) = 0.00025 (0.25)(0.005) = 0.00125 0.0015 Release Rates12
- 1.
The guidance for BWRs indicates that the durations of each of these release phases are:
Phase Duration Gap Release 0.5 hours5.787037e-5 days <br />0.00139 hours <br />8.267196e-6 weeks <br />1.9025e-6 months <br /> Early In-Vessel Release 1.5 hours5.787037e-5 days <br />0.00139 hours <br />8.267196e-6 weeks <br />1.9025e-6 months <br />
- 2.
The early in-vessel release phase begins after the gap release has ended.
- 3.
The guidance indicates that during each release phase, it is allowable to use a linear model of the activity release over the duration of the phase.
Accident Assumptions
- 1.
The highest power fuel plate is damaged under at least 23 feet of water.
- 2.
The reactor has been operated long enough for the fission product inventory to reach saturation.
6 Peak Power Fuel Plate Fission Fragment Inventory
- 1.
Highest Power Fuel Plate A.
The memo regarding the Steady-State Thermal-Hydraulic Analysis for Forced-Convective Flow in the Rhode Island Nuclear Science Center (RINSC) Reactor, indicates that power of the highest power plate during full power reactor operation at 2 MW is 9.653 kW.13
- 2.
Fission Rate A.
The energy associated with each fission event that occurs in the reactor is 200 MeV per fission.
B.
Converting from MeV to MW - Seconds:
Watt MW Joule Second Watt MeV Joule X
fission MeV 6
13 10 10 6.1 200
= 3.2 X 10-17 MW - second per fission C.
Therefore the fission rate at 1 MW power is:
1
[
3.2 1017
]
= 3.1 1016 D.
In order to generate 9.653 kW of power in the highest power fuel plate, the fission rate in that plate must be:
1 653
.9 1000 sec 10 1.3 16 kW kW MW ond MW fission X
= 2.992 X 1014 fission / second ~ 3 X 1014 fission / second
- 2.
Fission Nuclide Production Rate A.
The fission nuclide production rate for the i th fission product nuclide is the product of the fission rate and the fission product yield (i) for the i th fission product:
Fission Nuclide Production Rate = (3 X 1014 fission / second)( i)
7
- 3.
Fission Nuclide Decay Rate A.
The fission nuclide decay rate for the i th fission product nuclide is the product of the decay constant for the i th fission product nuclide (i), and the number of atoms of the nuclide that are present (Ni):
Fission Nuclide Decay Rate = (i)(Ni)
- 4.
Fission Product Saturation A.
Fission product saturation occurs when the production rate and decay rate are the same. Therefore, for the i th fission product, saturation is when:
(3 X 1014 fission / second)( i) = (i)(Ni sat)
B.
Therefore, if we wanted to estimate the number of atoms of the i th fission product in the highest power fuel plate at saturation, it would be:
i i
sat i
ond fission atoms ond fission X
N
sec sec 10 3
14 C.
However, if we want to estimate the activity of the i th fission product in the highest power fuel plate at saturation, it would be:
Activity (Bq) = (i)(Ni sat) = (3 X 1014 fission / second)( i)
D.
The activity can be converted to units of Ci by using the conversion factor:
1 Ci = 3.7 X 1010 Bq
- 5.
As an example, consider I-131 saturation in the core, which has a yield of = 0.0277 atoms per fission, and a decay constant of = 9.98 X 10-7 per second:
A.
Saturation Activity in the Highest Power Fuel Plate
fission atoms I
ond fission X
131 0277
.0 sec 10 3
14
= 8.3 X 1012 I-131 atoms per second
= 8.3 X 1012 Bq I-131
8
Bq X
Ci I
Bq X
10 12 10 7.3 1
131 10 3.8
= 2.25 X 102 Ci B.
Most of the fission products do not get out of the fuel matrix. Of the isotopes that get into the pool water, there is so much solvent in comparison to solute that the vast majority of the isotopes would stay dissolved in the pool water14.
The fission products that are both volatile and long lived enough to potentially escape from the fuel matrix to the pool, as well as the decay constants and cumulative fission yields of each of those isotopes are15:
Nuclide Fission Number Lamarsh Decay Highest Power Highest Power Yield of Seconds Half Life Constant Fuel Plate Activity Fuel Plate Activity (Atoms/Fission)
(s)
( /s)
(Bq)
(Ci)
I-131 0.0277 8.04 d
8.64E+04 6.95E+05 9.98E-07 8.31E+12 2.25E+02 I-132 0.0413 2.28 h
3.60E+03 8.21E+03 8.44E-05 1.24E+13 3.35E+02 I-133 0.0676 20.8 h
3.60E+03 7.49E+04 9.26E-06 2.03E+13 5.48E+02 I-134 0.0718 52.3 m
6.00E+01 3.14E+03 2.21E-04 2.15E+13 5.82E+02 I-135 0.0639 6.7 h
3.60E+03 2.41E+04 2.87E-05 1.92E+13 5.18E+02 Kr-85m 0.0133 4.4 h
3.60E+03 1.58E+04 4.38E-05 3.99E+12 1.08E+02 Kr-85 0.00285 10.76 y 3.15E+07 3.39E+08 2.04E-09 8.55E+11 2.31E+01 Kr-87 0.0237 76 m
6.00E+01 4.56E+03 1.52E-04 7.11E+12 1.92E+02 Kr-88 0.0364 2.79 h
3.60E+03 1.00E+04 6.90E-05 1.09E+13 2.95E+02 Xe-133m 0.00189 2.26 d
8.64E+04 1.95E+05 3.55E-06 5.67E+11 1.53E+01 Xe-133 0.0677 5.27 d
8.64E+04 4.55E+05 1.52E-06 2.03E+13 5.49E+02 Xe-135m 0.0105 15.7 m
6.00E+01 9.42E+02 7.36E-04 3.15E+12 8.51E+01 Xe-135 0.0672 9.2 h
3.60E+03 3.31E+04 2.09E-05 2.02E+13 5.45E+02 Fuel Plate Activity Lamarsh Half Life Note that the fission yields are cumulative, and include not only the yield of the nuclide, but also take into account the yields of the short lived precursors as well.
Activity Available to be Released from the Fuel16
- 1.
If a fuel plate were to be completely denuded on one side, the amount of the plate activity that would be available to be released into the coolant would depend on the temperature of the fuel, and the surface area of the fuel that is exposed. Since the fuel temperature even during full power operation is very low, diffusion from the fuel matrix into the coolant would be essentially zero. Consequently, the only fission products that would be available to get into the coolant would be from the denuded surface due to the kinetic energy associated with fission fragment recoil.
- 2.
As a result, of the total fission fragment inventory in the damaged fuel plate, the fraction of the activity that could be released to the coolant is the fraction that is in the volume defined by the surface area of the exposed fuel matrix, and the fission fragment recoil depth.
9
- 3.
Therefore, we can estimate the source term that is available to be released to the coolant by finding the ratio between the volume of the fuel meat that is within the recoil range, to the total volume of the fuel meat. Using the minimum volumes maximizes the activities in each volume, and is therefore conservative.
- 3.
NUREG / CR-2079 PNL-3691 Analysis of Credible Accidents for Argonaut Reactors indicates that the range of fission fragment recoil in aluminum is:
1.37 X 10-3 cm
- 4.
For the RINSC reactor, the minimum total volume of the fuel meat is17:
Length = (22.50 inches)(2.54 cm) / in = 57.15 cm Width = (2.32 inches)(2.54 cm) / in = 5.89 cm Thickness = (0.020 inches)(2.54 cm) / in = 0.05 cm Total Fuel Meat Volume = (57.15 cm)(5.89 cm)(0.05 cm) = 16.83 cm3
- 5.
If a fuel plate were stripped of one side of the cladding, the volume within the recoil range of the fission fragments would be:
(57.15 cm Length)(5.89 cm Width)(1.37 X 10-3 cm) = 0.46 cm3
- 6.
The ratio of fission fragments that are available to be released to the coolant, to the total fuel plate fission fragment inventory is:
(0.46 cm3)(16.83 cm3) = 0.027 = 2.7%
- 7.
Therefore, the activity that is available to be released from the fuel into the coolant is only 2.7% of the fuel plate activity. Overall:
Nuclide Highest Power Source Term Fuel Plate Activity (Ci)
(Ci)
I-131 2.25E+02 6.06E+00 I-132 3.35E+02 9.04E+00 I-133 5.48E+02 1.48E+01 I-134 5.82E+02 1.57E+01 I-135 5.18E+02 1.40E+01 Kr-85m 1.08E+02 2.91E+00 Kr-85 2.31E+01 6.24E-01 Kr-87 1.92E+02 5.19E+00 Kr-88 2.95E+02 7.97E+00 Xe-133m 1.53E+01 4.14E-01 Xe-133 5.49E+02 1.48E+01 Xe-135m 8.51E+01 2.30E+00 Xe-135 5.45E+02 1.47E+01 Source Term
10 Gap Release
- 1.
During Gap Release, the release fractions from the fuel to the coolant are:
A.
Iodine 0.05 B.
Noble Gas 0.05
- 2.
Example Calculation:
A.
For I-131 we concluded that the saturation activity in the hottest fuel plate would be 2.25 X 102 Ci.
B.
Of this, the activity in the fuel meat volume that is available to be released to the coolant is:
(2.25 X 102 Ci)(2.7%) = 6.06 Ci C.
The release fraction for both halogens and noble gases is 0.05 during gap release, so the amount that gets into the coolant during this phase of release is:
(0.05)(6.06 Ci) = 3.03 X 10-1 Ci I-131
- 3.
Overall:
Nuclide Source Coolant Term Activity (Ci)
(Ci)
I-131 6.06E+00 3.03E-01 I-132 9.04E+00 4.52E-01 I-133 1.48E+01 7.40E-01 I-134 1.57E+01 7.86E-01 I-135 1.40E+01 6.99E-01 Kr-85m 2.91E+00 1.46E-01 Kr-85 6.24E-01 3.12E-02 Kr-87 5.19E+00 2.59E-01 Kr-88 7.97E+00 3.98E-01 Xe-133m 4.14E-01 2.07E-02 Xe-133 1.48E+01 7.41E-01 Xe-135m 2.30E+00 1.15E-01 Xe-135 1.47E+01 7.36E-01 Coolant Activity From Gap Release
11
- 4.
Release Fractions from Coolant to Confinement Air A.
Iodine 0.005 B.
Noble Gas 1.0
- 5.
Example Calculations:
A.
- 1.
The amount of I - 131 the gets into the coolant during gap release is 3.03 X 10-1 Ci I-131.
- 2.
The amount that gets into the confinement air is:
(3.03 X 10-1 Ci)(0.005) = 1.52 X 10-3 Ci B.
Noble Gas
- 1.
The amount of Kr - 85m the gets into the coolant during gap release is 1.46 X 10-1 Ci.
- 2.
The amount that gets into the confinement air is:
(1.46 X 10-1 Ci)(1) = 1.46 X 10-1 Ci Kr - 85m
- 6.
Total Activity in Confinement Air Due to Gap Release Nuclide Coolant Confinement Air Activity Activity (Ci)
(Ci)
I-131 3.03E-01 1.52E-03 I-132 4.52E-01 2.26E-03 I-133 7.40E-01 3.70E-03 I-134 7.86E-01 3.93E-03 I-135 6.99E-01 3.50E-03 Kr-85m 1.46E-01 1.46E-01 Kr-85 3.12E-02 3.12E-02 Kr-87 2.59E-01 2.59E-01 Kr-88 3.98E-01 3.98E-01 Xe-133m 2.07E-02 2.07E-02 Xe-133 7.41E-01 7.41E-01 Xe-135m 1.15E-01 1.15E-01 Xe-135 7.36E-01 7.36E-01 Confinement Air Activity From Gap Release
12
- 7.
Duration of Gap Phase A.
The Gap Phase of release begins at time zero, and continues for 0.5 hours5.787037e-5 days <br />0.00139 hours <br />8.267196e-6 weeks <br />1.9025e-6 months <br />.
- 8.
Confinement Air Gap Activity Build-Up Rate A.
Rate of Activity Released = Total Activity Released / 0.5 Hours B.
For I-131 Rate of I-131 Release = (1.52 X 10-3) / (0.5 Hours) = 3.03 X 10-3 Ci / hr C.
Overall:
Nuclide Total Gap Release Gap Activity Confinement Air Build-Up Activity Rate (Ci)
(Ci / hr)
I-131 1.52E-03 3.03E-03 I-132 2.26E-03 4.52E-03 I-133 3.70E-03 7.40E-03 I-134 3.93E-03 7.86E-03 I-135 3.50E-03 6.99E-03 Kr-85m 1.46E-01 2.91E-01 Kr-85 3.12E-02 6.24E-02 Kr-87 2.59E-01 5.19E-01 Kr-88 3.98E-01 7.97E-01 Xe-133m 2.07E-02 4.14E-02 Xe-133 7.41E-01 1.48E+00 Xe-135m 1.15E-01 2.30E-01 Xe-135 7.36E-01 1.47E+00 Confinement Air Gap Activity Build-UP Rate
- 9.
Confinement Volume A negative pressure is maintained in the confinement building so that all of the air that exits the buiding will exit through a stack. If an airborne RAM release is detected, the Emergency Air Handling System is activated, and the airflow is directed through an emergency filter prior to reaching the stack.
During facility re-licensing, the volume of the confinement building was determined to be approximately 203,695 cubic feet. The volume of the pool structure and water was determined to be 21,740 cubic feet, leaving 181,955 cubic feet of open space.
The control room takes up about 3,612 cubic feet of this space. Converting the free volume of the confinement room to cubic centimeters:
13
3 9
3 3
3 3
3 10 15
.5 54
.2 12 1
181955 cm X
in cm ft in ft
- 10.
Confinement Air Concentration Build-Up Rate A.
If we make the assumption that all of the activity is dispersed uniformly in the confinement volume, the rate at which the concentrations of each isotope increase is:
(Activity Build-Up Rate) / (Confinement Volume)
B.
For I-131 Concentration Build-Up Rate = (3.03 X 10-3 Ci / hr) / (5.15 X 109 cc)
= 5.88 X 10-13 Ci / cc per hour C.
Overall Gap Release Confinement Air Confinement Air Confinement Air Gap Activity Gap Concentration Activity Build-Up Rate Build-Up Rate (Ci)
(Ci / hr)
(Ci / cc-hr)
I-131 1.52E-03 3.03E-03 5.89E-13 I-132 2.26E-03 4.52E-03 8.78E-13 I-133 3.70E-03 7.40E-03 1.44E-12 I-134 3.93E-03 7.86E-03 1.53E-12 I-135 3.50E-03 6.99E-03 1.36E-12 Kr-85m 1.46E-01 2.91E-01 5.65E-11 Kr-85 3.12E-02 6.24E-02 1.21E-11 Kr-87 2.59E-01 5.19E-01 1.01E-10 Kr-88 3.98E-01 7.97E-01 1.55E-10 Xe-133m 2.07E-02 4.14E-02 8.03E-12 Xe-133 7.41E-01 1.48E+00 2.88E-10 Xe-135m 1.15E-01 2.30E-01 4.46E-11 Xe-135 7.36E-01 1.47E+00 2.86E-10 Gap Phase Confinement Activity Concentration Rate (0 to 0.5 Hours)
- 11.
Emergency Filter Release Fraction A.
When the Emergency Air Handling System is activated, all of the air from the confinement room is exhausted through an emergency filter. The emergency filter consists of:
14
- 1.
Roughing Filter
- 2.
HEPA Filter
- 3.
Charcoal Filter
- 4.
HEPA Filter B.
The proposed Technical Specifications associated with the Emergency Filter efficiency are:
3.5.2.3 The emergency filter shall be at least 99% efficient at removing iodine.
3.5.2.4 The Emergency Filter System Exhaust Absolute Filter shall be certified by the manufacturer to have a minimum efficiency of 99.97% for removing 0.3 micron diameter particulates.
C.
There are two sets of HEPA filters that the concentration of airborne particles must go through, so th concnetration of particles that are 0.3 microns or greater in the confinement room that will reach the stack is:
(0.03%)(0.03%)(Confinement Concentration)
= (0.0009%)(Confinement Concentration)
= (0.00009)(Confinement Concentration)
= (9 X 10-5)(Confinement Concentration)
D.
Of the iodine that reaches the stack, only the non-organic iodine is adsorbed by the charcoal filter. The organic iodine is not adsorbed.
E.
The concentration of non-organic iodine in the confinement room that will reach the stack is:
(1%)(Confinement Concentration)
= (0.01)(Confinement Concentration)
F.
The noble gases are unaffected by the HEPA filters, but are slowed by the charcoal filter. We will assume that all of the noble gases are released to the stack.
G.
If we continue with our I-131 example:
- 1.
We concluded that the concentration build-up rate for I-131 in the confinement air was 5.89E-13 Ci / cm3 per hour.
15
- 2.
Therefore, if the concentration were to build up for the entire 0.5 hr duration of the gap release, without exhausting any of it through the stack, the concentration would be:
(5.89 X 10-13 Ci / cc - hour)(0.5 hour5.787037e-5 days <br />0.00139 hours <br />8.267196e-6 weeks <br />1.9025e-6 months <br />) = 2.95 X 10-13 Ci / cc
- 3.
Of this concentration of iodine, the relative amounts of organic and non-organic iodine are:
(43% Organic Iodine)( 2.95 X 10-13 Ci / cc) = 1.27 X 10-13 Ci / cc (57% Non-Organic Iodine)( 2.95 X 10-13 Ci / cc) = 1.68 X 10-13 Ci / cc
- 4.
All of the organic Iodine is assumed to get through the filter, and only 1% of the non-organic iodine is assumed to get through. Therefore, the total amount that gets through the filter is:
1.27 X 10-13 Ci / cc + (0.01)(1.68 X 10-13 Ci / cc) = 1.28 X 10-13 Ci / cc
- 5.
Overall Confinement Air Final Concentration Concentration Emergency Filter Gap Concentration Gap Activity of Organic of Non-Organic Release Build-Up Rate Concentration Iodine Iodine Concentration (Ci / cc-hr)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
I-131 5.89E-13 2.94E-13 1.27E-13 1.68E-13 1.28E-13 I-132 8.78E-13 4.39E-13 1.89E-13 2.50E-13 1.91E-13 I-133 1.44E-12 7.18E-13 3.09E-13 4.09E-13 3.13E-13 I-134 1.53E-12 7.63E-13 3.28E-13 4.35E-13 3.32E-13 I-135 1.36E-12 6.79E-13 2.92E-13 3.87E-13 2.96E-13 Kr-85m 5.65E-11 2.83E-11 2.83E-11 Kr-85 1.21E-11 6.06E-12 6.06E-12 Kr-87 1.01E-10 5.04E-11 5.04E-11 Kr-88 1.55E-10 7.74E-11 7.74E-11 Xe-133m 8.03E-12 4.02E-12 4.02E-12 Xe-133 2.88E-10 1.44E-10 1.44E-10 Xe-135m 4.46E-11 2.23E-11 2.23E-11 Xe-135 2.86E-10 1.43E-10 1.43E-10 Emergency Filter Release Concentrations
- 12.
Therefore, the concentration of the nuclides that enter into the stack, is equivalent to the concentrations that are released from the emergency filter.
- 13.
Radioactive Material Release Rate A.
The dispersion calculations use radioactive material release rates, rather than activity concentrations.
B.
The rate at which RAM is exhausted from the stack is equivalent to the rate at which it is exhausted from the emergency filter.
16 C.
The higher the filter exhaust flow rate, the higher the RAM release rate from the stack.
D.
Technical Specification 3.5.2.2 limits flow through the emergency filter to a maximum of 1500 cfm.
s cm X
s in cm ft in ft
/
10 08
.7 60 min 54
.2 12 min 1500 3
5 3
3 3
3 3
E.
Therefore, the RAM release rate from the stack is:
s Ci s
cm Flowrate Exhaust Filter Emergency cm Ci ion Concentrat Exhaust Filter Emergency
3 3
F.
Continuing with the I-131 example, the concentation of I-131 in the emergency filter exhaust was 1.28 X 10-13 Ci / cm3, so the I-131 release rate is:
s Ci X
s cm X
cm Ci X
/
10 06
.9 10 08
.7 10 28
.1 8
3 5
3 13
E.
Overall:
Initial Stack Stack RAM Release Concentration Rate (Ci / cc)
(Ci / s)
I-131 1.28E-13 9.08E-08 I-132 1.91E-13 1.35E-07 I-133 3.13E-13 2.22E-07 I-134 3.32E-13 2.35E-07 I-135 2.96E-13 2.09E-07 Kr-85m 2.83E-11 2.00E-05 Kr-85 6.06E-12 4.29E-06 Kr-87 5.04E-11 3.57E-05 Kr-88 7.74E-11 5.48E-05 Xe-133m 4.02E-12 2.84E-06 Xe-133 1.44E-10 1.02E-04 Xe-135m 2.23E-11 1.58E-05 Xe-135 1.43E-10 1.01E-04 Stack RAM Release Rate
17
- 14.
Atmospheric Dispersion A.
Air Turbulance Conditions
- 1.
Non-Fumigation Conditions Atmospheric stability is a measure of the turbulence in the plume, and it affects the rate of the dispersion of the plume. The more turbulent the air is, the greater the dispersion rate. There are six classifications of atmospheric stability, ranging from Pasquill Type A through Pasquill Type F, in which A is extremely unstable, and F is moderately stable. While this suggests that Pasquill Type F is conservative because it minimizes the dispersion rate and maximizes the airborne RAM concentrations at ground level, it is not necessarily the most conservative from the standpoint of the distance at which the plume reaches ground level. Consequently, Pasquill Types B, D, and F were considered.
- 2.
Fumigation Conditions Fumigation is a condition in which there is an air inversion above the stack release which forces the plume to stay below the inversion. As a result, under fumigation conditions, the plume reaches the ground level at distances that are much shorter than what occurs under non-fumigation conditions. NRC Regulatory Guide 1.145 Section 2.1.2.b indicates that for coastal sites that are located less than 3.2 kilometers from a large body of water, fumigation conditions should be assumed to exist. Since the stack is approximately 150 meters from the Narragansett Bay, fumigation conditions were also considered.
- 3.
Conditions Under Consideration Initially for this analysis, an air turbulance condition of Pasquill F type was considered because it was assumed that this level of air stability would result in the highest concnetration of RAM when the plume reached ground level. NRC pointed out that this assumption might not give the most conservative results because under more turbulent conditions, the plume reaches ground level much sooner than under stable conditions. They suggested that stability types Pasquill B and D be considered. Additionally, given the facilitys proximity to a large body of water, the guidance indicated that fumigation conditions must also be considered. Consequently, non-fumigation conditions Pasquill B, D, and F were considered, as well as the fumigation condition.
18 B.
Dispersion Equations Atmoshperic dispersion calculations estimate the concentration of some material in air for a given release rate, under specified atmospheric conditions, at some distance away from the source. A Gaussian Straight Line Plume Model is used.
- 1.
Non-Fumigation Equation Section 1.3.2.a of NRC Regulatory Guide 1.145 (November 1982 Rev.
- 1) indicates that the equation for stack releases under nonfumigation conditions is:
2 2
2 1
z eh h
z y
e u
Q
- 2.
Fumigation Equation Section 1.3.2.b of NRC Regulatory Guide 1.145 (November 1982 Rev.
- 1) indicates that the equation for stack releases under nonfumigation conditions is:
19 e
he y
h u
Q
2
/
1) 2
(
1
- 3.
The factors for these equations are:
A.
is the concetration of the material in the air (Activity /
Volume)
B.
Q is the material release rate (Activity / Time)
C.
y is the horizontal dispersion coefficient (Distance)
D.
z is the vertical dispersion coefficient (Distance)
E.
hs is the stack height above the plant grade (Distance)
F.
ht is the maximum terrain height above plant grade between the release point, and the point of interest (Distance)
G.
he is the effective stack height = hs - ht with he = 0 if ht> hs (Distance)
H.
uh is the average wind speed of the plume at the release height (Distance / Time)
I.
uhe is the average wind speed of the fumigation layer at the effective stack height (Distance / Time)
J.
x is the downwind distance (Distance)
K.
y is the horizontal distance at right angles to the plume centerline (Distance)
L.
z is the height above the ground (Distance)
C.
Stack Height The release is at the level of the stack (hs = 115 ft), with no significant change in terrain height between the release point and the site boundary (ht = 0), so the effective stack height is:
he = hs - ht = 35m = 0 = 115 ft
= (115 ft)(12 in / ft)(2.54 cm / in)(m / 100 cm)
= 35.052 m D.
Wind Direction Wind direction is assumed to be constant, As a result, all of the concentration of RAM will be along one line of direction, rather than dispersed across more than one direction. Consequently, this assumption is conservative.
E.
Wind Speed The wind speed is assumed to be one meter per second. Higher wind speeds increase dilution because the RAM released per unit time is added to a larger
20 volume of air passing by the release point. Consequently, this assumption is conservative.
F.
Coordinates of Interest The dispersion coefficents are quantitative measures of how much the plume spreads out in the horizontal (y) and vertical (z) directions at some given distance (x) down wind of the release point. The material concentration in the plume as a function of distance from the plume centerline is a Gaussian distribution, with maximum concentration at the centerline. The dispersion coeficients are the standard deviations in each direction. Therefore, 68% of the plume is within y distance from the centerline in the y - direction, and z distance from the centerline in the z - direction. For this analysis, the y and z values are taken to be zero so that the concentrations calculated will be associated with the centerline of the plume, and therefore will be the most conservative.
The distance between the confinement stack and the site boundary is 50 meters. RINSC has the authority to prevent the public from entering this boundary.
The facility sits on the Bay Campus of the University of Rhode Island. There are residential areas to the north and south of the campus. To the west of the campus, woods and businesses abut the campus. To the northwest along South Ferry Road, there are residences as well. Using Google Maps to estimate the distances between the stack release point, and the residences to the north, nothwest, and south, the nearest residence was determined to be 500 meters from the stack:
Nearest Residence on South Ferry Road
21 Distance to Nearest Residence in Saounderstown Nearest Residence in Bonnet Shores
22 Therefore the closest residence is 1600 feet from the stack:
(1600 ft)(m / 3.28 ft) = 487.8 m 500 m The maximum concentration point occurs whenNeed
Reference:
2(z)2 = (he)2
= [
2 2 ]
1/2 From this, the maximum concentration point is the x distance from the stack that corresponds to the z that was calculated.
In addition to calculating where the maximum concentration point occurs, numerical methods were used to verify that itt occurs where it is predicted to occur.
Consequently, the three points down wind of the stack that we are interested in knowing the RAM concentrations at ground level are:
- 1.
Site Boundary (x = 50 meters)
- 2.
Nearest Residence (x = 500 meters)
- 3.
Maximum concentration point The dispersion coefficients take these distances into account. The dispersion coefficient curves given in the US Atomic Energy Commissions Meteorology and Atomic Energy, 1968, pp. 102 - 103 are:
23
24
25 G.
Non-Fumigation Example Calculation for Site Boundary Concentrations Under Pasquill F Conditions
- 1.
The Site Boundary distance from the stack is 50 meters.
- 2.
The y dispersion coefficient at 50 meters is:
Using linear extrapolation to determine what y is at 50 meters for the Pasquill F curve on the y vs. x graph, two data points are:
x1 = 100 m y1 = 4 m x2 = 200 m y2 = 8 m The slope of the line is:
=
(2 1)
(2 1) =
(8 4 )
(200 100 ) =
4 100 = 0.04 The point - slope form of the line is:
(y2 - y1) = m(x2 - x1)
The y - intercept occurs at point (0,b):
(y2 - b) = m(x2 - 0)
(y2 - b) = mx2 y2 - mx2 = b b = y2 - mx2 b = (8 m) - (0.04)(200 m) = 0.0 m Consequently, the general equation is:
y = mx + b y = (0.04) x + 0.0 m For our case x = 50 m, so y is:
y = y = (0.016)(50 m) + 0.08 m = 2.0 m
26
- 3.
The z dispersion coefficient at 50 meters is:
Using linear extrapolation to determine what y is at 50 meters for the Pasquill F curve on the y vs. x graph, two data points are:
x1 = 100 m y1 = 2.4 m x2 = 200 m y2 = 4 m The slope of the line is:
=
(2 1)
(2 1) =
(4 2.4 )
(200 100 ) =
1.6 100 = 0.016 The point - slope form of the line is:
(y2 - y1) = m(x2 - x1)
The y - intercept occurs at point (0,b):
(y2 - b) = m(x2 - 0)
(y2 - b) = mx2 y2 - mx2 = b b = y2 - mx2 b = (4 m) - (0.016)(200 m) = 0.8 m Consequently, the general equation is:
y = mx + b y = (0.016) x + 0.8 m For our case x = 50 m, so z is:
z = y = (0.016)(48 m) + 0.08 m = 1.6 m
- 4.
Therefore, the dispersion coefficients for a distance of 50 meters from the stack, under non-fumigation conditions, with Pasquill Type F turbulance are:
y = 2.0 meters z = 1.6 meters
27
- 5.
Consequently, the values of each of the factors in the non-fumigation dispersion equation are:
he = 35.052 m uh = 1 m / s y = 2.0 m z = 1.6 m
- 6.
Plugging these values into the non-fumigation dispersion equation:
2 2
2 1
z eh h
z y
e u
Q
2 2
)
6.1
(
2
)
052 35
(
)
/
1()
6.1()
2
(
1 m
m e
s m
m m
Q
=
1 10.0531 3 / 239.97 Q
= 6.03 X 10-106 s / m3 Converting the volume to units of cm3:
Q
= (6.03 X 10-106 s / m3)(m / 100 cm)3 = 6.03 X10-112 s / cm3
- 7.
Using this to convert the stack release rate of I-131 into a concentration:
For I-131, the release rate out of the stack was 9.08 X 10-8 Ci / s.
(9.08 X 10-8 Ci / s)( 6.03 X10-112 s / cm3) = 5.4 X 10-119 Ci / cm3 This is the concentration of I-131 at ground level, at a distance of 50 meters down wind of the stack, if the air turbulance condition is Pasquill F, and the wind speed is 1 m / s.
H.
Methodology for Determining Concentrations at the Site Boundary, Nearest Residiential, and Maximum Concentration Points
- 1.
In order to get a sense of where the maximum concentration point is, dispersion coefficients were determined for distances ranging from the 50 meter site boundary point, to 10,000 meters out from the stack, including the 500 meter residential point.
28
- 2.
X/Q was calculated for each of these distances.
- 3.
Isotope concentrations were calculated for each of these distances, and the distance with the highest concentration data was highlighted.
- 4.
A second set of dispersion coefficients were determined for distances at smaller increments around the point of maximum concentration found in the first set of data, in order to get a more precise location.
- 5.
X/Q was calculated for each of these new distances.
- 6.
Isotope concentrations were calculated for each of these new distances, and the distance with the highest concentration data was highlighted.
The distance associated with these concentrations is the maximum concentration point.
- 7.
Additionally, the point of maximum concentration was also calculated using the following formulaNeed formal reference for Ronald Fjeld, Quantitative Environmental Risk Analysis :
2(z)2 = (he)2 and solving for z.
The z vs. x graph was used to determine the distance x associated with the calculated value of z. This distance is the predicted point of maximum concentration point.
- 8.
For each case, both methods for determining the maximum concentration point were found to be consistent with each other.
29 H.
Pasquill B Condition Results Stack Site Boundary Nearest Residence RAM Release 50 m 200 m 500 m 1000 m 2000 m 5000 m 10000 m
=
3.141593 Rate Concentration Concentration Concentration Concentration Concentration Concentration Concentration (Ci / s)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc) h(s)
=
115 ft
=
35.052 m I-131 9.08E-08 2.54E-20 8.90E-18 5.09E-18 1.36E-18 2.72E-19 1.45E-20 0.00E+00 I-132 1.35E-07 3.79E-20 1.33E-17 7.58E-18 2.03E-18 4.06E-19 2.17E-20 0.00E+00 h(t)
=
0 ft
=
0 m I-133 2.22E-07 6.20E-20 2.17E-17 1.24E-17 3.32E-18 6.65E-19 3.55E-20 0.00E+00 I-134 2.35E-07 6.59E-20 2.31E-17 1.32E-17 3.53E-18 7.06E-19 3.77E-20 0.00E+00 h(e)
=
115 ft
=
35.052 m I-135 2.09E-07 5.87E-20 2.05E-17 1.17E-17 3.14E-18 6.28E-19 3.35E-20 0.00E+00 Kr-85m 2.00E-05 5.60E-18 1.96E-15 1.12E-15 3.00E-16 6.00E-17 3.20E-18 0.00E+00 h(e)^2
=
=
1228.643 m^2 Kr-85 4.29E-06 1.20E-18 4.20E-16 2.40E-16 6.43E-17 1.29E-17 6.86E-19 0.00E+00 Kr-87 3.57E-05 9.99E-18 3.50E-15 2.00E-15 5.35E-16 1.07E-16 5.71E-18 0.00E+00 (y)
=
85 m Kr-88 5.48E-05 1.53E-17 5.37E-15 3.07E-15 8.22E-16 1.64E-16 8.76E-18 0.00E+00 Xe-133m 2.84E-06 7.96E-19 2.79E-16 1.59E-16 4.27E-17 8.53E-18 4.55E-19 0.00E+00 (z)
=
55 m Xe-133 1.02E-04 2.85E-17 9.98E-15 5.71E-15 1.53E-15 3.06E-16 1.63E-17 0.00E+00 Xe-135m 1.58E-05 4.42E-18 1.55E-15 8.85E-16 2.37E-16 4.74E-17 2.53E-18 0.00E+00 (z)^2
=
3025 m^2 Xe-135 1.01E-04 2.83E-17 9.91E-15 5.66E-15 1.52E-15 3.03E-16 1.62E-17 0.00E+00 U(h)
=
1 m/s Paquill 50 m 200 m 500 m 1000 m 2000 m 5000 m 10000 m
/Q
=
5.56E-05 s/m^3
=
5.6E-11 s/cm^3 B
X/Q X/Q X/Q X/Q X/Q X/Q X/Q (s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(y) 9.5 35 85 160 300 650 1250 (z) 8 20 55 130 350 3000 X/Q 2.8E-13 9.8E-11 5.6E-11 1.5E-11 3.0E-12 1.6E-13 h(e)^2
=
1228.643 m^2 (z)^2
=
614.3214 m^2 Stack RAM Release 50 m 100 m 150 m 200 m 250 m 300 m 350 m 400 m 500 m (z)
=
25 m Rate Concentration Concentration Concentration Concentration Concentration Concentration Concentration Concentration Concentration (Ci / s)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc) x[(z)=25]
250 m I-131 9.08E-08 2.54E-20 1.91E-18 4.99E-18 8.90E-18 1.18E-17 9.99E-18 8.26E-18 6.99E-18 5.09E-18 I-132 1.35E-07 3.79E-20 2.84E-18 7.45E-18 1.33E-17 1.76E-17 1.49E-17 1.23E-17 1.04E-17 7.58E-18 I-133 2.22E-07 6.20E-20 4.65E-18 1.22E-17 2.17E-17 2.88E-17 2.44E-17 2.02E-17 1.71E-17 1.24E-17 I-134 2.35E-07 6.59E-20 4.94E-18 1.29E-17 2.31E-17 3.06E-17 2.59E-17 2.14E-17 1.81E-17 1.32E-17 I-135 2.09E-07 5.87E-20 4.40E-18 1.15E-17 2.05E-17 2.72E-17 2.30E-17 1.91E-17 1.61E-17 1.17E-17 (y)
Kr-85m 2.00E-05 5.60E-18 4.20E-16 1.10E-15 1.96E-15 2.60E-15 2.20E-15 1.82E-15 1.54E-15 1.12E-15 x1
=
100 m Kr-85 4.29E-06 1.20E-18 9.01E-17 2.36E-16 4.20E-16 5.58E-16 4.72E-16 3.90E-16 3.30E-16 2.40E-16 x2
=
200 m Kr-87 3.57E-05 9.99E-18 7.49E-16 1.96E-15 3.50E-15 4.64E-15 3.92E-15 3.25E-15 2.75E-15 2.00E-15 (y1)
=
18 m Kr-88 5.48E-05 1.53E-17 1.15E-15 3.01E-15 5.37E-15 7.12E-15 6.03E-15 4.98E-15 4.22E-15 3.07E-15 (y2)
=
35 m Xe-133m 2.84E-06 7.96E-19 5.97E-17 1.56E-16 2.79E-16 3.70E-16 3.13E-16 2.59E-16 2.19E-16 1.59E-16 Xe-133 1.02E-04 2.85E-17 2.14E-15 5.60E-15 9.98E-15 1.32E-14 1.12E-14 9.27E-15 7.84E-15 5.71E-15 m
=
0.17 Xe-135m 1.58E-05 4.42E-18 3.32E-16 8.69E-16 1.55E-15 2.05E-15 1.74E-15 1.44E-15 1.22E-15 8.85E-16 b
=
1 Xe-135 1.01E-04 2.83E-17 2.12E-15 5.56E-15 9.91E-15 1.31E-14 1.11E-14 9.20E-15 7.79E-15 5.66E-15 x(50 m)
=
9.5 (z)
Pasquill 50 m 100 m 150 m 200 m 250 m 300 m 350 m 400 m 500 m x1
=
100 m B
X/Q X/Q X/Q X/Q X/Q X/Q X/Q X/Q X/Q x2
=
200 m (s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(z1)
=
12 m (z2)
=
20 m (y) 9.5 18 25 35 40 50 60 68 85 (z) 8 12 15 20 28 32 38 44 55 m
=
0.08 b
=
4 X/Q 2.8E-13 2.1E-11 5.5E-11 9.8E-11 1.3E-10 1.1E-10 9.1E-11 7.7E-11 5.6E-11 x(50 m)
=
8 Atmospheric Dispersion Calculator Calculated Maximum Concentration Point 2(z)^2 = h(e)^2 Linear Extrapolation Data for X/Q at 50 m
30 I.
Pasquill D Condition Results Stack Site Boundary Nearest Residence RAM Release 50 m 200 m 500 m 1000 m 2000 m 5000 m 10000 m
=
3.141593 Rate Concentration Concentration Concentration Concentration Concentration Concentration Concentration (Ci / s)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc) h(s)
=
115 ft
=
35.052 m I-131 9.08E-08 5.45E-45 9.99E-20 7.81E-18 6.90E-18 3.09E-18 1.18E-18 3.45E-19 I-132 1.35E-07 8.12E-45 1.49E-19 1.16E-17 1.03E-17 4.60E-18 1.76E-18 5.14E-19 h(t)
=
0 ft
=
0 m I-133 2.22E-07 1.33E-44 2.44E-19 1.91E-17 1.68E-17 7.53E-18 2.88E-18 8.42E-19 I-134 2.35E-07 1.41E-44 2.59E-19 2.02E-17 1.79E-17 8.00E-18 3.06E-18 8.94E-19 h(e)
=
115 ft
=
35.052 m I-135 2.09E-07 1.26E-44 2.30E-19 1.80E-17 1.59E-17 7.12E-18 2.72E-18 7.96E-19 Kr-85m 2.00E-05 1.20E-42 2.20E-17 1.72E-15 1.52E-15 6.80E-16 2.60E-16 7.61E-17 h(e)^2
=
=
1228.643 m^2 Kr-85 4.29E-06 2.57E-43 4.72E-18 3.69E-16 3.26E-16 1.46E-16 5.58E-17 1.63E-17 Kr-87 3.57E-05 2.14E-42 3.92E-17 3.07E-15 2.71E-15 1.21E-15 4.64E-16 1.36E-16 (y)
=
72 m Kr-88 5.48E-05 3.29E-42 6.03E-17 4.71E-15 4.16E-15 1.86E-15 7.12E-16 2.08E-16 Xe-133m 2.84E-06 1.71E-43 3.13E-18 2.45E-16 2.16E-16 9.67E-17 3.70E-17 1.08E-17 (z)
=
33 m Xe-133 1.02E-04 6.11E-42 1.12E-16 8.76E-15 7.74E-15 3.46E-15 1.32E-15 3.87E-16 Xe-135m 1.58E-05 9.48E-43 1.74E-17 1.36E-15 1.20E-15 5.37E-16 2.05E-16 6.00E-17 (z)^2
=
1089 m^2 Xe-135 1.01E-04 6.07E-42 1.11E-16 8.70E-15 7.69E-15 3.44E-15 1.31E-15 3.84E-16 U(h)
=
1 m/s Paquill 50 m 200 m 500 m 1000 m 2000 m 5000 m 10000 m
/Q
=
7.62E-05 s/m^3
=
7.6E-11 s/cm^3 D
X/Q X/Q X/Q X/Q X/Q X/Q X/Q (s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(y) 4 16 40 72 145 250 550 (z) 3 9
20 33 50 90 150 X/Q 6.00E-38 1.1E-12 8.6E-11 7.6E-11 3.4E-11 1.3E-11 3.80E-12 h(e)^2
=
1228.643 m^2 (z)^2
=
614.3214 m^2 Stack RAM Release 200 m 300 m 400 m 500 m 600 m 700 m 800 m 900 m 1000 m (z)
=
25 m Rate Concentration Concentration Concentration Concentration Concentration Concentration Concentration Concentration Concentration (Ci / s)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc) x[(z)=25]
700 m I-131 9.08E-08 9.99E-20 3.63E-18 5.09E-18 7.81E-18 7.90E-18 8.35E-18 7.90E-18 7.45E-18 6.90E-18 I-132 1.35E-07 1.49E-19 5.42E-18 7.58E-18 1.16E-17 1.18E-17 1.25E-17 1.18E-17 1.11E-17 1.03E-17 I-133 2.22E-07 2.44E-19 8.86E-18 1.24E-17 1.91E-17 1.93E-17 2.04E-17 1.93E-17 1.82E-17 1.68E-17 I-134 2.35E-07 2.59E-19 9.42E-18 1.32E-17 2.02E-17 2.05E-17 2.17E-17 2.05E-17 1.93E-17 1.79E-17 I-135 2.09E-07 2.30E-19 8.38E-18 1.17E-17 1.80E-17 1.82E-17 1.93E-17 1.82E-17 1.72E-17 1.59E-17 (y)
Kr-85m 2.00E-05 2.20E-17 8.01E-16 1.12E-15 1.72E-15 1.74E-15 1.84E-15 1.74E-15 1.64E-15 1.52E-15 x1
=
100 m Kr-85 4.29E-06 4.72E-18 1.72E-16 2.40E-16 3.69E-16 3.73E-16 3.95E-16 3.73E-16 3.52E-16 3.26E-16 x2
=
200 m Kr-87 3.57E-05 3.92E-17 1.43E-15 2.00E-15 3.07E-15 3.10E-15 3.28E-15 3.10E-15 2.92E-15 2.71E-15 (y1)
=
8 m Kr-88 5.48E-05 6.03E-17 2.19E-15 3.07E-15 4.71E-15 4.77E-15 5.04E-15 4.77E-15 4.49E-15 4.16E-15 (y2)
=
16 m Xe-133m 2.84E-06 3.13E-18 1.14E-16 1.59E-16 2.45E-16 2.47E-16 2.62E-16 2.47E-16 2.33E-16 2.16E-16 Xe-133 1.02E-04 1.12E-16 4.08E-15 5.71E-15 8.76E-15 8.86E-15 9.37E-15 8.86E-15 8.35E-15 7.74E-15 m
=
0.08 Xe-135m 1.58E-05 1.74E-17 6.32E-16 8.85E-16 1.36E-15 1.37E-15 1.45E-15 1.37E-15 1.30E-15 1.20E-15 b
=
0 Xe-135 1.01E-04 1.11E-16 4.04E-15 5.66E-15 8.70E-15 8.80E-15 9.30E-15 8.80E-15 8.29E-15 7.69E-15 x(50 m)
=
4 (z)
Pasquill 200 m 300 m 400 m 500 m 600 m 700 m 800 m 900 m 1000 m x1
=
100 m D
X/Q X/Q X/Q X/Q X/Q X/Q X/Q X/Q X/Q x2
=
200 m (s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(z1)
=
5 m (z2)
=
9 m (y) 16 25 32 40 47 52 60 65 72 (z) 9 14 16 20 22 25 28 30 33 m
=
0.04 b
=
1 X/Q 1.1E-12 4.0E-11 5.6E-11 8.6E-11 8.7E-11 9.2E-11 8.7E-11 8.2E-11 7.6E-11 x(50 m)
=
3 Atmospheric Dispersion Calculator Calculated Maximum Concentration Point 2(z)^2 = h(e)^2 Linear Extrapolation Data for X/Q at 50 m
31 J.
Pasquill F Condition Results Stack Site Boundary Nearest Residence RAM Release 50 m 200 m 500 m 1000 m 2000 m 5000 m 10000 m
=
3.141592654 Rate Concentration Concentration Concentration Concentration Concentration Concentration Concentration (Ci / s)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc) h(s)
=
115 ft
=
35.052 m I-131 9.08E-08 5.45E-119 1.91E-32 3.45E-20 3.18E-18 5.27E-18 3.09E-18 1.63E-18 I-132 1.35E-07 8.12E-119 2.84E-32 5.14E-20 4.74E-18 7.85E-18 4.60E-18 2.44E-18 h(t)
=
0 ft
=
0 m I-133 2.22E-07 1.33E-118 4.65E-32 8.42E-20 7.76E-18 1.29E-17 7.53E-18 3.99E-18 I-134 2.35E-07 1.41E-118 4.94E-32 8.94E-20 8.24E-18 1.37E-17 8.00E-18 4.24E-18 h(e)
=
115 ft
=
35.052 m I-135 2.09E-07 1.26E-118 4.40E-32 7.96E-20 7.33E-18 1.21E-17 7.12E-18 3.77E-18 Kr-85m 2.00E-05 1.20E-116 4.20E-30 7.61E-18 7.00E-16 1.16E-15 6.80E-16 3.60E-16 h(e)^2
=
=
1228.643 m^2 Kr-85 4.29E-06 2.57E-117 9.01E-31 1.63E-18 1.50E-16 2.49E-16 1.46E-16 7.72E-17 Kr-87 3.57E-05 2.14E-116 7.49E-30 1.36E-17 1.25E-15 2.07E-15 1.21E-15 6.42E-16 (y)
=
160 m Kr-88 5.48E-05 3.29E-116 1.15E-29 2.08E-17 1.92E-15 3.18E-15 1.86E-15 9.86E-16 Xe-133m 2.84E-06 1.71E-117 5.97E-31 1.08E-18 9.95E-17 1.65E-16 9.67E-17 5.12E-17 (z)
=
35 m Xe-133 1.02E-04 6.11E-116 2.14E-29 3.87E-17 3.57E-15 5.91E-15 3.46E-15 1.83E-15 Xe-135m 1.58E-05 9.48E-117 3.32E-30 6.00E-18 5.53E-16 9.16E-16 5.37E-16 2.84E-16 (z)^2
=
1225 m^2 Xe-135 1.01E-04 6.07E-116 2.12E-29 3.84E-17 3.54E-15 5.87E-15 3.44E-15 1.82E-15 U(h)
=
1 m/s Paquill 50 m 200 m 500 m 1000 m 2000 m 5000 m 10000 m
/Q
=
F X/Q X/Q X/Q X/Q X/Q X/Q X/Q
=
3.44246E-05 s/m^3
=
3.4E-11 s/cm^3 (s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(y) 2 8
20 40 70 160 280 (z) 1.6 4
8.5 15 22 35 46 X/Q 6.00E-112 2.1E-25 3.8E-13 3.5E-11 5.8E-11 3.4E-11 1.80E-11 h(e)^2
=
1228.642704 m^2 (z)^2
=
614.321352 m^2 Stack RAM Release 1000 m 1500 m 2000 2500 m 3000 m 3500 m 4000 m 4500 m 5000 m (z)
=
25 m Rate Concentration Concentration Concentration Concentration Concentration Concentration Concentration Concentration Concentration (Ci / s)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc) x[(z)=25]
2500 m I-131 9.08E-08 3.18E-18 4.09E-18 5.27E-18 5.45E-18 4.72E-18 4.09E-18 3.54E-18 3.36E-18 3.09E-18 I-132 1.35E-07 4.74E-18 6.09E-18 7.85E-18 8.12E-18 7.04E-18 6.09E-18 5.28E-18 5.01E-18 4.60E-18 I-133 2.22E-07 7.76E-18 9.97E-18 1.29E-17 1.33E-17 1.15E-17 9.97E-18 8.64E-18 8.20E-18 7.53E-18 I-134 2.35E-07 8.24E-18 1.06E-17 1.37E-17 1.41E-17 1.22E-17 1.06E-17 9.18E-18 8.71E-18 8.00E-18 I-135 2.09E-07 7.33E-18 9.43E-18 1.21E-17 1.26E-17 1.09E-17 9.43E-18 8.17E-18 7.75E-18 7.12E-18 (y)
Kr-85m 2.00E-05 7.00E-16 9.01E-16 1.16E-15 1.20E-15 1.04E-15 9.01E-16 7.81E-16 7.41E-16 6.80E-16 x1
=
100 m Kr-85 4.29E-06 1.50E-16 1.93E-16 2.49E-16 2.57E-16 2.23E-16 1.93E-16 1.67E-16 1.59E-16 1.46E-16 x2
=
200 m Kr-87 3.57E-05 1.25E-15 1.60E-15 2.07E-15 2.14E-15 1.85E-15 1.60E-15 1.39E-15 1.32E-15 1.21E-15 (y1)
=
4 m Kr-88 5.48E-05 1.92E-15 2.46E-15 3.18E-15 3.29E-15 2.85E-15 2.46E-15 2.14E-15 2.03E-15 1.86E-15 (y2)
=
8 m Xe-133m 2.84E-06 9.95E-17 1.28E-16 1.65E-16 1.71E-16 1.48E-16 1.28E-16 1.11E-16 1.05E-16 9.67E-17 Xe-133 1.02E-04 3.57E-15 4.58E-15 5.91E-15 6.11E-15 5.30E-15 4.58E-15 3.97E-15 3.77E-15 3.46E-15 m
=
0.04 Xe-135m 1.58E-05 5.53E-16 7.11E-16 9.16E-16 9.48E-16 8.22E-16 7.11E-16 6.16E-16 5.85E-16 5.37E-16 b
=
0 Xe-135 1.01E-04 3.54E-15 4.55E-15 5.87E-15 6.07E-15 5.26E-15 4.55E-15 3.94E-15 3.74E-15 3.44E-15 x(50 m)
=
2 (z)
Pasquill 1000 m 1500 m 2000 m 2500 m 3000 m 3500 m 4000 m 4500 m 5000 m x1
=
100 m F
X/Q X/Q X/Q X/Q X/Q X/Q X/Q X/Q X/Q x2
=
200 m (s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(z1)
=
2.4 m (z2)
=
4 m (y) 40 50 70 80 100 120 140 150 160 (z) 15 17 22 25 28 30 32 34 35 m
=
0.016 b
=
0.8 X/Q 3.5E-11 4.5E-11 5.8E-11 6.0E-11 5.2E-11 4.5E-11 3.9E-11 3.7E-11 3.4E-11 x(50 m)
=
1.6 Atmospheric Dispersion Calculator Calculated Maximum Concentration Point 2(z)^2 = h(e)^2 Linear Extrapolation Data for X/Q at 50 m
32 K.
Fumigation Example Calculation for Site Boundary Concentrations Under Pasquill F Conditions
- 1.
NRC Regulatory Guide 1.145 section 1.3.2.b provides the dispersion equation for fumigation conditions, and indicates that:
A.
e he y
h u
Q
2
/
1) 2
(
1
B.
A wind speed of u = 2 meters / second is considered to be conservative. Since the initial conservative assumption was made that the wind speed is u = 1 meter / second, and this is more conservative that an assuption of 2 meters / second, the analysis was done with u = 1 meter / second.
C.
The lateral dispersion coefficient y is determined under air turbulent conditions of Pasquill F.
- 2.
The Site Boundary distance from the stack is 50 meters.
- 3.
The y dispersion coefficient at 50 meters under Pasquill F conditions was found to be 2 meters.
- 4.
Consequently, the values of each of the factors in the fumigation dispersion equation are:
he = 35.052 m uhe = 1 m / s y = 2.0 m
- 6.
Plugging these values into the non-fumigation dispersion equation:
)
052 35
()
/
1()
2
(
)
2
(
1
)
2
(
1 2
/
1 2
/
1 m
s m
m h
u Q
e he y
= 5.69 X 10-3 s / m3 Converting the volume to units of cm3:
Q
= (5.69 X 10-3 s / m3)(m / 100 cm)3 = 5.69 X10-9 s / cm3
- 7.
Using this to convert the stack release rate of I-131 into a concentration:
33 For I-131, the release rate out of the stack was 9.08 X 10-8 Ci / s.
(9.08 X 10-8 Ci / s)(5.69 X10-9 s / cm3) = 5.17 X 10-16 Ci / cm3 This is the concentration of I-131 at ground level, at a distance of 50 meters down wind of the stack, if the air turbulance is in a fumigation condition, and the wind speed is 1 m / s.
34 L.
Fumigation Results Stack Site Boundary Nearest Residence RAM Release 50 m 200 m 500 m 1000 m 2000 m 5000 m 10000 m
=
3.141593 Rate Concentration Concentration Concentration Concentration Concentration Concentration Concentration (Ci / s)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc) 2
=
6.283185 I-131 9.08E-08 5.18E-16 1.27E-16 5.18E-17 2.54E-17 1.45E-17 6.45E-18 3.72E-18 I-132 1.35E-07 7.72E-16 1.90E-16 7.72E-17 3.79E-17 2.17E-17 9.61E-18 5.55E-18 (2)^(1/2) =
2.506628 I-133 2.22E-07 1.26E-15 3.10E-16 1.26E-16 6.20E-17 3.55E-17 1.57E-17 9.09E-18 I-134 2.35E-07 1.34E-15 3.30E-16 1.34E-16 6.59E-17 3.77E-17 1.67E-17 9.65E-18 h(s)
=
115 ft
=
35.052 m I-135 2.09E-07 1.19E-15 2.93E-16 1.19E-16 5.87E-17 3.35E-17 1.49E-17 8.59E-18 Kr-85m 2.00E-05 1.14E-13 2.80E-14 1.14E-14 5.60E-15 3.20E-15 1.42E-15 8.21E-16 h(t)
=
0 ft
=
0 m Kr-85 4.29E-06 2.44E-14 6.00E-15 2.44E-15 1.20E-15 6.86E-16 3.04E-16 1.76E-16 Kr-87 3.57E-05 2.03E-13 4.99E-14 2.03E-14 9.99E-15 5.71E-15 2.53E-15 1.46E-15 h(e)
=
115 ft
=
35.052 m Kr-88 5.48E-05 3.12E-13 7.67E-14 3.12E-14 1.53E-14 8.76E-15 3.89E-15 2.25E-15 Xe-133m 2.84E-06 1.62E-14 3.98E-15 1.62E-15 7.96E-16 4.55E-16 2.02E-16 1.17E-16 (y)
=
6 m Xe-133 1.02E-04 5.81E-13 1.43E-13 5.81E-14 2.85E-14 1.63E-14 7.23E-15 4.18E-15 Xe-135m 1.58E-05 9.01E-14 2.21E-14 9.01E-15 4.42E-15 2.53E-15 1.12E-15 6.48E-16 U(he)
=
1 m/s Xe-135 1.01E-04 5.76E-13 1.42E-13 5.76E-14 2.83E-14 1.62E-14 7.18E-15 4.15E-15 Paquel 50 m 200 m 500 m 1000 m 2000 m 5000 m 10000 m
/Q
=
1.90E-03 s/m^3
=
1.9E-09 s/cm^3 F
X/Q X/Q X/Q X/Q X/Q X/Q X/Q (s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(y) 2 8
20 40 70 160 280 X/Q 5.7E-09 1.4E-09 5.7E-10 2.8E-10 1.6E-10 7.1E-11 4.1E-11 h(e)^2
=
m^2 (z)^2
=
m^2 Stack RAM Release 50 m 100 m 150 m 200 m 250 m 300 m 350 m 400 m 500 m (z)
=
m Rate Concentration Concentration Concentration Concentration Concentration Concentration Concentration Concentration Concentration (Ci / s)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc) x[(z)=25]
m I-131 9.08E-08 5.18E-16 2.54E-16 1.73E-16 1.27E-16 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 I-132 1.35E-07 7.72E-16 3.79E-16 2.57E-16 1.90E-16 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 I-133 2.22E-07 1.26E-15 6.20E-16 4.21E-16 3.10E-16 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 I-134 2.35E-07 1.34E-15 6.59E-16 4.47E-16 3.30E-16 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 I-135 2.09E-07 1.19E-15 5.87E-16 3.98E-16 2.93E-16 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 (y)
Kr-85m 2.00E-05 1.14E-13 5.60E-14 3.80E-14 2.80E-14 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 x1
=
100 m Kr-85 4.29E-06 2.44E-14 1.20E-14 8.15E-15 6.00E-15 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 x2
=
200 m Kr-87 3.57E-05 2.03E-13 9.99E-14 6.78E-14 4.99E-14 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 (y1)
=
4 m Kr-88 5.48E-05 3.12E-13 1.53E-13 1.04E-13 7.67E-14 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 (y2)
=
8 m Xe-133m 2.84E-06 1.62E-14 7.96E-15 5.40E-15 3.98E-15 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 Xe-133 1.02E-04 5.81E-13 2.85E-13 1.94E-13 1.43E-13 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 m
=
0.04 Xe-135m 1.58E-05 9.01E-14 4.42E-14 3.00E-14 2.21E-14 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 b
=
0 Xe-135 1.01E-04 5.76E-13 2.83E-13 1.92E-13 1.42E-13 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 x(50 m)
=
2 (z)
Pasquel 50 m 100 m 150 m 200 m 250 m 300 m 350 m 400 m 500 m x1
=
100 m F
X/Q X/Q X/Q X/Q X/Q X/Q X/Q X/Q X/Q x2
=
200 m (s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(s / cc)
(z1)
=
2.4 m (z2)
=
4 m (y) 2 4
6 8
m
=
0.016 b
=
0.8 X/Q 5.7E-09 2.8E-09 1.9E-09 1.4E-09 x(50 m)
=
1.6 Atmospheric Dispersion Calculator Calculated Maximum Concentration Point 2(z)^2 = h(e)^2 Linear Extrapolation Data for X/Q at 50 m
35 M.
Results Summary Fumigation Fumigation Pasquill Pasquill Pasquill Pasquill Pasquill Pasquill Pasquill Pasquill B
D F
F B
D F
F 50 m 50 m 50 m 50 m 200 m 500 m 2500 m 50 m Concentration Concentration Concentration Concentration Concentration Concentration Concentration Concentration (Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
I-131 2.54E-20 5.45E-45 5.45E-119 5.18E-16 I-131 8.90E-18 7.81E-18 5.45E-18 5.18E-16 I-132 3.79E-20 8.12E-45 8.12E-119 7.72E-16 I-132 1.33E-17 1.16E-17 8.12E-18 7.72E-16 I-133 6.20E-20 1.33E-44 1.33E-118 1.26E-15 I-133 2.17E-17 1.91E-17 1.33E-17 1.26E-15 I-134 6.59E-20 1.41E-44 1.41E-118 1.34E-15 I-134 2.31E-17 2.02E-17 1.41E-17 1.34E-15 I-135 5.87E-20 1.26E-44 1.26E-118 1.19E-15 I-135 2.05E-17 1.80E-17 1.26E-17 1.19E-15 Kr-85m 5.60E-18 1.20E-42 1.20E-116 1.14E-13 Kr-85m 1.96E-15 1.72E-15 1.20E-15 1.14E-13 Kr-85 1.20E-18 2.57E-43 2.57E-117 2.44E-14 Kr-85 4.20E-16 3.69E-16 2.57E-16 2.44E-14 Kr-87 9.99E-18 2.14E-42 2.14E-116 2.03E-13 Kr-87 3.50E-15 3.07E-15 2.14E-15 2.03E-13 Kr-88 1.53E-17 3.29E-42 3.29E-116 3.12E-13 Kr-88 5.37E-15 4.71E-15 3.29E-15 3.12E-13 Xe-133m 7.96E-19 1.71E-43 1.71E-117 1.62E-14 Xe-133m 2.79E-16 2.45E-16 1.71E-16 1.62E-14 Xe-133 2.85E-17 6.11E-42 6.11E-116 5.81E-13 Xe-133 9.98E-15 8.76E-15 6.11E-15 5.81E-13 Xe-135m 4.42E-18 9.48E-43 9.48E-117 9.01E-14 Xe-135m 1.55E-15 1.36E-15 9.48E-16 9.01E-14 Xe-135 2.83E-17 6.07E-42 6.07E-116 5.76E-13 Xe-135 9.91E-15 8.70E-15 6.07E-15 5.76E-13 Fumigation Pasquill Pasquill Pasquill Pasquill B
D F
F 500 m 500 m 500 m 500 m Concentration Concentration Concentration Concentration (Ci / cc)
(Ci / cc)
(Ci / cc)
(Ci / cc)
I-131 5.09E-18 7.81E-18 3.45E-20 5.18E-17 I-132 7.58E-18 1.16E-17 5.14E-20 7.72E-17 I-133 1.24E-17 1.91E-17 8.42E-20 1.26E-16 I-134 1.32E-17 2.02E-17 8.94E-20 1.34E-16 I-135 1.17E-17 1.80E-17 7.96E-20 1.19E-16 Kr-85m 1.12E-15 1.72E-15 7.61E-18 1.14E-14 Kr-85 2.40E-16 3.69E-16 1.63E-18 2.44E-15 Kr-87 2.00E-15 3.07E-15 1.36E-17 2.03E-14 Kr-88 3.07E-15 4.71E-15 2.08E-17 3.12E-14 Xe-133m 1.59E-16 2.45E-16 1.08E-18 1.62E-15 Xe-133 5.71E-15 8.76E-15 3.87E-17 5.81E-14 Xe-135m 8.85E-16 1.36E-15 6.00E-18 9.01E-15 Xe-135 5.66E-15 8.70E-15 3.84E-17 5.76E-14 Nearest Residence Concentration Non-Fumigation Conditions Non-Fumigation Conditions Non-Fumigation Conditions Site Boundary Concentration Maximum Concentration Point
36
- 1.
A comparison of the data for non-fumigation Pasquill B, D, and F conditions, with the fumigation condition shows that in every case, the fumigation condition is the most conservative condition. Furthermore, under fumigation conditions, the maximum concentration point is the site boundary. Since the nearest residence is ten times the distance between the stack and the site boundary, the concentrations at the nearest residence are exactly one order of magnatude less than the concnetrations at the site boundary.
N.
RAM Concenterations at Points of Interest Site Boundary and Nearest Maximum Concentration Point Residence Concentration Concentration (Ci / cc)
(Ci / cc)
I-131 5.18E-16 5.18E-17 I-132 7.72E-16 7.72E-17 I-133 1.26E-15 1.26E-16 I-134 1.34E-15 1.34E-16 I-135 1.19E-15 1.19E-16 Kr-85m 1.14E-13 1.14E-14 Kr-85 2.44E-14 2.44E-15 Kr-87 2.03E-13 2.03E-14 Kr-88 3.12E-13 3.12E-14 Xe-133m 1.62E-14 1.62E-15 Xe-133 5.81E-13 5.81E-14 Xe-135m 9.01E-14 9.01E-15 Xe-135 5.76E-13 5.76E-14 RAM Concentrations at Points of Interest
- 15.
Dose Consequences A.
Dose Calculation Background Information
- 1.
Health effects of radiation dose are separated into two categories:
A.
Stochastic Effects - These effects are probabilistic, and are due to random ionization events. Consequently, there is no threshold for these effects, and the probability of occurrence is proportional to the dose received. Cancer is an example of these types of effects.
B.
Non-Stochastic Effects - These effects depend on the amount of dose received beyond a minimum threshold, and the amount
37 of damage depends on the magnatude of the dose. Skin erythmia is an example of a non-stochastic effect.
- 2.
The objective of dose limits are to minimize the risk of stochastic effects, and to prevent the occurrence of non-stochastic effects. The dose limits have been designed to be independent of whether or not the radiation dose is uniform or non-uniform. This is achieved by having effective dose limits in which the effective dose takes into consideration the risk due to the irradiation of each individual organ and equates it to the risk associated with a uniform irradiation of the whole body.
B.
Definitions
- 1.
Allowable Limit on Intake (ALI) - This is the amoung of RAM taken into the body via ingestion or inhalation that would lead to a committed effective dose equivalent of 5 Rem, or 50 Rem to any individual tissue or organ.
- 2.
Derived Air Concentration (DAC) - This is the concentration of a given radionuclide in air which if inhaled at a rate of 2 X 104 cm3 per minute for one working year (2000 hours0.0231 days <br />0.556 hours <br />0.00331 weeks <br />7.61e-4 months <br />) would result in reaching the ALI.
Therefore, as an example, I-131 has a DAC = 2 X 10-8 Ci / cm3 and an ALI = 50 Ci so the relationship between ALI and DAC is:
Ci Ci hr hr cm X
cm Ci X
50 48 1
2000 min 60 min 10 2
10 2
3 4
3 8
This means that if the reference man were to breath in the DAC for 2000 hours0.0231 days <br />0.556 hours <br />0.00331 weeks <br />7.61e-4 months <br />, at a rate of 2 x 104 Ci / minute, then they will have an uptake equivalent to the ALI.
- 3.
Derived Air Concentration - Hour (DAC - Hour) -This is the product of the concentration of RAM in the air expressed as a fraction or multiple of the DAC, and the exposure time expressed in hours.
- 4.
Absorbed Dose (D) - This is a measure of the radiation energy that is absorbed per unit mass of material of interest.
- 5.
Dose Equivalent (H) - This is the product of the absorbed dose in tissue, quality factor (Q), and all other necessary modifying factors at the location of interest. The units of dose equivalent are the rem and sievert (Sv). In general:
H = DQ
38
- 6.
Quality Factor (Q) - This is a regulatoraly defined factor to account for the fact that the type and energy of the incident radiation has an effect on the amount of biological damage that is produced per unit of absorbed energy (absorbed dose).
- 7.
Tissue Dose Equivalent (HT) - This is the dose equivalent to a specific tissue or organ due to external sources.
- 8.
Committed Dose Equivalent (HT,50) - This is the dose equivalent to organs or tissues of reference (T) that will be received from a single intake of radioactive material by an individual that will be accumulated over the 50-year period following the intake.
- 9.
Effective Dose Equivalent (HE) - This equates the risk of a non-uniform external dose, or internal dose to the risk associated with a dose that is distributed uniformly over the whole body. A regulatoraly defined weighting factor (WT) is used for each organ, and the overall effective dose equivalent is:
HE = WTHT This is the sum of the products of the dose equivalent to the organ or tissue (HT) and the weighting factors (WT) applicable to each of the body organs or tissues that are irradiated (HE = WTHT).
- 10.
Committed Effective Dose Equivalent (HE,50) - This is the effective dose equivalent accumulated over a 50 year period as a result of a single intake of radioactive material. In general:
HE,50 = WTHT,50 This is the sum of the products of the weighting factors applicable to each of the body organs or tissues that are irradiated and the committed dose equivalent to these organs or tissues (HE,50 =
WTHT,50).
- 11.
Deep Dose Equivalent (DDE) - This is the whole body dose at a depth of 1 cm due to an external exposure.
- 12.
Total Effective Dose Equivalent (TEDE) - This is the sum of the DDE due to an external dose, and the CEDE due to an internal dose from an intake of radioactive material. In general:
39 C.
Regulatory Limits10 CFR 20
- 1.
Occupational Dose Limits A.
TEDE = 5 rem / yr [10 CFR 20.1201(a)(1)(i)]
B.
DDE + CDE to any individual organ or tissue = 50 rem / yr [10 CFR 20.1201(a)(1)(ii)]
C.
The DAC and ALI may be used to determine the individuals dose and to demonstrate compliance with dose limits. [10 CFR 20.1201(d)]
D.
If the only intake of radionuclides is by inhalation, the total effective dose equivalent limit is not exceeded if the sum of the deep-dose equivalent divided by the total effective dose equivalent limit, and one of the following, does not exceed unity [10 CFR 20.1202(b)]:
(1)
The sum of the fractions of the inhalation ALI for each radionuclide, or (2)
The total number of derived air concentration-hours (DAC-hours) for all radionuclides divided by 2,000, or (3)
The sum of the calculated committed effective dose equivalents to all significantly irradiated1 organs or tissues (T) calculated from bioassay data using appropriate biological models and expressed as a fraction of the annual limit.
E.
If the identity and concentration of each radionuclide in a mixture are known, the fraction of the DAC applicable to the mixture for use in calculating DAC-hours must be either [10 CFR 20.1204(e)]:
(1)
The sum of the ratios of the concentration to the appropriate DAC value (e.g., D, W, Y) from appendix B to part 20 for each radionuclide in the mixture; or (2)
The ratio of the total concentration for all radionuclides in the mixture to the most restrictive DAC value for any radionuclide in the mixture.
F.
In order to calculate the committed effective dose equivalent, the licensee may assume that the inhalation of one ALI, or an exposure of 2,000 DAC-hours, results in a committed effective dose equivalent of 5 rems (0.05 Sv) for radionuclides that have their ALIs or DACs based on the committed effective dose equivalent. [10 CFR 20.1204(h)(1)]
40 G.
When the ALI (and the associated DAC) is determined by the nonstochastic organ dose limit of 50 rems (0.5 Sv), the intake of radionuclides that would result in a committed effective dose equivalent of 5 rems (0.05 Sv) (the stochastic ALI) is listed in parentheses in table 1 of appendix B to part 20. In this case, the licensee may, as a simplifying assumption, use the stochastic ALIs to determine committed effective dose equivalent. However, if the licensee uses the stochastic ALIs, the licensee must also demonstrate that the limit in § 20.1201(a)(1)(ii) is met. [10 CFR 20.1204(h)(2)]
- 2.
Dose Limits for Individual Members of the Public A.
TEDE = 100 mrem / yr [10 CFR 20.1301(a)(1)]
B.
A licensee shall show compliance with the annual dose limit in
§ 20.1301 by [10 CFR 20.1302(b)]:
(1)
Demonstrating by measurement or calculation that the total effective dose equivalent to the individual likely to receive the highest dose from the licensed operation does not exceed the annual dose limit; or (2)
Demonstrating that:
(i)
The annual average concentrations of radioactive material released in gaseous and liquid effluents at the boundary of the unrestricted area do not exceed the values specified in table 2 of appendix B to part 20; and (ii)
If an individual were continuously present in an unrestricted area, the dose from external sources would not exceed 0.002 rem (0.02 mSv) in an hour and 0.05 rem (0.5 mSv) in a year.
D.
External Immersion Dose vs. Internal Dose For the fuel failure accident we are concerned about the doses that individuals will recieve due to airborne radioactive materials. The airborne RAM that is released in these types of accidents consist of halogens, such as iodine, and noble gases, such as xenon and krypton.
When halogens are inhaled, part of what is inhaled is taken up and incorporated into the body. Consequently, these isotopes cause not only an external immersion dose, but also a committed internal dose.
41 Noble gases are inert, so when they are inhaled, they are not taken up and incorporated into the body. Consequently, these isotopes only cause an external immersion dose and do not contribute to an internal dose.
During a fuel failure accident, the principle halogen that is released as an airborne RAM source is iodine. When iodine is uptaken into the body, it concentrates in the thyroid. As a result, the internal dose associated with a fuel failure accident would be the Committed Dose Equivalent (HT,50) to the thyroid due to iodine.
During a fuel failure accident, the principle noble gases that are released as airborne RAM sources are krypton and xenon. These isotopes, in addition to the iodine isotopes are the are the sources for the external immersion dose. As a result, the external immersion dose that is associated with a fuel failure accident is the Deep Dose Equivalent (DDE) to the whole body due to the iodine, krypton, and xenon isotopes.
E.
What We Must Show
- 1.
Therefore, based on the regulations, we must show that:
A.
The occupational doses to individuals inside confinement are no greater than:
- 1.
- 2.
DDE + CDE to any individual organ or tissue = 50 rem
/ yr B.
The dose to the public at the site boundary is no greater than:
- 1.
TEDE = 100 mrem / yr F.
Use of the DAC to Determine Deep Dose Equivalent (DDE)
The DAC for any given nuclide can be found in 10 CFR 20 Appendix B. For Kr-85, the inhalation value of the DAC for occupational exposure are given to be:
DAC = 1 X 10-4 Ci / cm3 This means that if an individual is immersed in a concentration of 1 X 10-4 Ci / cm3 Kr-85 for 2000 hours0.0231 days <br />0.556 hours <br />0.00331 weeks <br />7.61e-4 months <br />, they would receive a DDE of 5 Rem whole body.
If and individual is immersed in an air concentration equivalent to one DAC, the average dose rate that they would be receiving would be:
42 Whole Body:
hr mrem hr rem
/
5.2
)
2000
(
5
Since this is the dose rate that is associated with an airborne RAM concentration of one DAC, we can relate the dose rate associated with being immersed in an airborne RAM concentration of one DAC by:
[2.5
] = 2.5 /
We can express the concentration of any RAM in air as a fraction or multiple of the DAC:
DAC Fraction (Multiple) = (Air Concentration) / (DAC)
Therefore, if we had a concentration of 6.06 X 10-6Ci / cm3 of Kr-85 in the confinement air (the confinement concentration after 5 minutes), the DAC fraction (multiple) if the occupational DAC were 1 X 10-4Ci / cm3 would be:
DAC Fraction (Multiple) = (Air Concentration) / (DAC)
= [6.06 106 3
] [
85 1 104 /3]
= 6.06 X 10-2 DAC If we know the concentration of RAM in air as a fraction or multiple of the DAC, then we can determine the dose rate that that an individual immersed in the air would recieve. Continuing with the Kr-85 example:
hr mrem X
DAC X
hr DAC mrem
/
10 5.1 1
10 06
.6 5.2 1
2
For a mixture of airborne radionuclides, the total dose rate can be determined either by summing the individual DAC fractions and multiplying the sum by the dose rate per DAC - hr:
)
/
(
1 5.2 hr mrem Rate Dose Total Fractions DAC hr DAC mrem
43 Or by finding the dose rate associated with each nuclide and summing the individual dose rates to get the total dose rate:
[(2.5
) (
1
)] = (/)
Consider:
- 1.
Suppose that the air has a concentration of 6.06 X 10-6Ci / cm3 of Kr-85, and 2.83 X 10-5Ci / cm3 of Kr-85m in it. The DAC fractions are:
(Air Concentration) / (DAC)
Where the DAC is defined in 10 CFR 20 for each isotope
- 2.
For Kr-85 the DAC fraction has been previously calculated to be 6.06 X 10-2 DAC.
- 3.
For Kr-85m, given that the DAC is 2 X 10-5Ci / cm3, the DAC fraction is:
(2.83 X 10-5Ci / cm3) / (2 X 10-5Ci / cm3) = 1.42 DAC
- 4.
Therefore the total DAC fraction is:
Total DAC Fraction = 6.06 X 10-2 DAC + 1.42 DAC = 1.48 DAC
- 5.
Therefore the deep dose equivalent is:
/
7.3 1
48
.1 5.2
G.
Deep Dose Equivalent to Personnel Inside Confinement for the First Five Minutes After Fuel Facility evacuation drills show that in the event of an evacuation, the confinement building can be evacuated within 5 minutes. Using the confinement concentration after five minutes, and taking into consideration that individuals inside confinement will only be exposed for five minutes, the projected dose from the gas release for individuals inside confinement when the fuel failure occurs is:
)
(
min 60 1
min 5
1 5.2 mrem DDE hr Fraction DAC hr DAC mrem
44 For the isotopes of interest:
Nuclide Confinement Occupational DAC Immersion Deep Dose Concentration DAC Fraction DAC Equivalent (microCi / cc)
(microCi / cc)
(DAC-hr) mrem I-131 2.94E-07 2.00E-08 1.47E+01 1.23E+00 3.07E+00 I-132 4.39E-07 3.00E-06 1.46E-01 1.22E-02 3.05E-02 I-133 7.18E-07 1.00E-07 7.18E+00 5.99E-01 1.50E+00 I-134 7.63E-07 2.00E-05 3.82E-02 3.18E-03 7.95E-03 I-135 6.79E-07 7.00E-07 9.70E-01 8.08E-02 2.02E-01 Kr-85m 2.83E-05 2.00E-05 1.41E+00 1.18E-01 2.94E-01 Kr-85 6.06E-06 1.00E-04 6.06E-02 5.05E-03 1.26E-02 Kr-87 5.04E-05 5.00E-06 1.01E+01 8.40E-01 2.10E+00 Kr-88 7.74E-05 2.00E-06 3.87E+01 3.22E+00 8.06E+00 Xe-133m 4.02E-06 1.00E-04 4.02E-02 3.35E-03 8.37E-03 Xe-133 1.44E-04 1.00E-04 1.44E+00 1.20E-01 3.00E-01 Xe-135m 2.23E-05 9.00E-06 2.48E+00 2.07E-01 5.17E-01 Xe-135 1.43E-04 1.00E-05 1.43E+01 1.19E+00 2.98E+00 Confinement Air Immersion Dose Therefore, if someone remains in confinement for five minutes, they will receive a deep dose equivalent of:
Individual DDE = 19.1 mrem G.
Deep Dose Equivalent to General Public at the Site Boundary Section 1.3 of NRC Regulatory Guide 1.145 (November 1982 Rev. 1) indicates that a conservative estimation of site boundary doses can be determined by it assuming that an individual spends two hours at the site boundary immediately following an accident. The deep dose equivalent would be the sum of the dose equivalents from each of the isotopes of interest.
For the isotopes of interest:
45 Nuclide Site Boundary Site Boundary Occupational DAC Immersion Deep Dose Concentration Concentration DAC Fraction DAC Equivalent (Ci / cc)
(microCi / cc)
(microCi / cc)
(DAC-hr) mrem I-131 5.18E-16 5.18E-10 2.00E-08 2.59E-02 5.18E-02 1.29E-01 I-132 7.72E-16 7.72E-10 3.00E-06 2.57E-04 5.14E-04 1.29E-03 I-133 1.26E-15 1.26E-09 1.00E-07 1.26E-02 2.53E-02 6.32E-02 I-134 1.34E-15 1.34E-09 2.00E-05 6.71E-05 1.34E-04 3.35E-04 I-135 1.19E-15 1.19E-09 7.00E-07 1.71E-03 3.41E-03 8.53E-03 Kr-85m 1.14E-13 1.14E-07 2.00E-05 5.70E-03 1.14E-02 2.85E-02 Kr-85 2.44E-14 2.44E-08 1.00E-04 2.44E-04 4.89E-04 1.22E-03 Kr-87 2.03E-13 2.03E-07 5.00E-06 4.07E-02 8.13E-02 2.03E-01 Kr-88 3.12E-13 3.12E-07 2.00E-06 1.56E-01 3.12E-01 7.81E-01 Xe-133m 1.62E-14 1.62E-08 1.00E-04 1.62E-04 3.24E-04 8.11E-04 Xe-133 5.81E-13 5.81E-07 1.00E-04 5.81E-03 1.16E-02 2.90E-02 Xe-135m 9.01E-14 9.01E-08 9.00E-06 1.00E-02 2.00E-02 5.00E-02 Xe-135 5.76E-13 5.76E-07 1.00E-05 5.76E-02 1.15E-01 2.88E-01 Site Boundary Air Immersion Dose Therefore, if a member of the general public remains at the site boundary for two hours, they will receive a dose of:
Individual DDE = 1.58 mrem H.
Deep Dose Equivalent to General Public at the Nearest Residence Again, assuming a two hour occupancy time, the deep dose equivalent at the nearest residency is:
Nuclide Nearest Residence Nearest Residence Occupational DAC Immersion Deep Dose Concentration Concentration DAC Fraction DAC Equivalent (Ci / cc)
(microCi / cc)
(microCi / cc)
(DAC-hr) mrem I-131 5.18E-17 5.18E-11 2.00E-08 2.59E-03 5.18E-03 1.29E-02 I-132 7.72E-17 7.72E-11 3.00E-06 2.57E-05 5.14E-05 1.29E-04 I-133 1.26E-16 1.26E-10 1.00E-07 1.26E-03 2.53E-03 6.32E-03 I-134 1.34E-16 1.34E-10 2.00E-05 6.71E-06 1.34E-05 3.35E-05 I-135 1.19E-16 1.19E-10 7.00E-07 1.71E-04 3.41E-04 8.53E-04 Kr-85m 1.14E-14 1.14E-08 2.00E-05 5.70E-04 1.14E-03 2.85E-03 Kr-85 2.44E-15 2.44E-09 1.00E-04 2.44E-05 4.89E-05 1.22E-04 Kr-87 2.03E-14 2.03E-08 5.00E-06 4.07E-03 8.13E-03 2.03E-02 Kr-88 3.12E-14 3.12E-08 2.00E-06 1.56E-02 3.12E-02 7.81E-02 Xe-133m 1.62E-15 1.62E-09 1.00E-04 1.62E-05 3.24E-05 8.11E-05 Xe-133 5.81E-14 5.81E-08 1.00E-04 5.81E-04 1.16E-03 2.90E-03 Xe-135m 9.01E-15 9.01E-09 9.00E-06 1.00E-03 2.00E-03 5.00E-03 Xe-135 5.76E-14 5.76E-08 1.00E-05 5.76E-03 1.15E-02 2.88E-02 Nearest Residence Air Immersion Dose Therefore, if a member of the general public remains at the nearest residence for two hours, they will receive a dose of:
Individual DDE = 1.58 X 10-1 mrem
46 I.
Use of the DAC to Determine the Committed Dose Equivalent to the Thyroid (CDE)
The ALI and DAC for any given nuclide can be found in 10 CFR 20 Appendix B. For I-131, the inhalation values for occupational exposure are given to be:
- 1.
ALI = 50 Ci This means that an intake of 50 Ci of I-131 will lead to a CEDE of 5 Rem, or 50 Rem to any individual tissue or organ. Since iodine concentrates in the thyroid, the ALI is based on a dose of 50 Rem to the thyroid.
- 2.
DAC = 2 X 10-8 Ci / cm3 This means that if an individual inhales concentration of 2 X 10-8 Ci /
cm3 I-131 at a rate of 2 X 104 cm3 per minute for 2000 hours0.0231 days <br />0.556 hours <br />0.00331 weeks <br />7.61e-4 months <br />, they would intake enough of the radionuclide to receive a CEDE of 5 Rem whole body, or 50 Rem to any individual tissue or organ:
If an individual is immersed in an air concentration equivalent to one DAC, the average dose rate that they would be receiving would be:
- 1.
Individual Tissue or Organ (in this case Thyroid):
/
25 2000 50 Since this is the dose rate that is associated with an airborne RAM concentration of one DAC, we can relate the dose rate associated with being immersed in an airborne RAM concentration of one DAC by:
[25
] = 25 /
We can express the concentration of any RAM in air as a fraction or multiple of the DAC:
DAC Fraction (Multiple) = (Air Concentration) / (DAC)
Therefore, if we had a concentration of 2.94 X 10-7Ci / cm3 of I-131 in the confinement air, the DAC fraction (multiple) if the occupational DAC were 2 X 10-8Ci / cm3 would be:
47 DAC Fraction (Multiple) = (Air Concentration) / (DAC)
= [2.94 107 3
] [
131 2 108 /3]
= 14.7 DAC If we know the concentration of RAM in air as a fraction or multiple of the DAC, then we can determine the committed dose rate that that an individual immersed in the air would recieve. Continuing with the I-131 example:
/
5.
367 1
7.
14 25
If there is more than one nuclide in the air with the same dose rate associated with exposure (either whole body or individual organ), then the DAC fractions can be added together before determining the dose rate. Consider:
- 1.
Suppose that the air has a concentration of 2.94 X 10-7Ci / cm3 of I-131, and 4.39 X 10-7Ci / cm3 of I-132 in it. The DAC fractions are:
(Air Concentration) / (DAC)
Where the DAC is defined in 10 CFR 20 for each isotope
- 2.
For I-131 the DAC fraction has been previously calculated to be 14.7 DAC.
- 3.
For I-132, given that the DAC is 3 X 10-6Ci / cm3, the DAC fraction is:
(4.39 X 10-7Ci / cm3) / (3 X 10-6Ci / cm3) = 1.46 X 10-1 DAC
- 4.
Therefore the total DAC fraction is:
Total DAC Fraction = 14.7 DAC + 1.46 X 10-1 DAC = 14.8 DAC
- 5.
For the iodines, the committed dose to the thyroid is also dependent on the amount of time that the individual is immersed. If an individual were only in the concentration of iodine for 5 minutes (0.083 hr), then the DAC fraction can be reduced:
Immersion DAC = (DAC)(Immersion Time)
Immersion DAC = (14.8 DAC)(0.083 hr) = 1.23 DAC - hr
48
- 6.
Both of these DACs are based on a committed thyroid dose of 50 rem per year, which means that the dose rate associated with an air concentration of one DAC is 25 mrem / DAC - hr.
- 7.
Therefore the committed dose to the thyroid for an individual that is immersed for 5 minutes in air with a concentration of 2.94 X 10-7Ci /
cm3 of I-131 and 4.39 X 10-7Ci / cm3 of I-132 would be:
23
.1 25
J.
Confinement Committed Dose to the Thyroid (CDE)
Halogens are an inhalation hazard because they are absorbed into the body.
The halogen of interest in the case of a fuel failure is iodine. Iodine concentrates in the thyroid. Consequently, the DAC for each isotope of Iodine is based on the amount of isotope that will result in a 50 rem dose to the thyroid over a 2000 hour0.0231 days <br />0.556 hours <br />0.00331 weeks <br />7.61e-4 months <br /> year. We have calculated that an individual immersed in air with a concentration of RAM in it equivalent to one DAC would lead to an internal dose rate of 25 mrem / hr to the thyroid.
Facility evacuation drills show that in the event of an evacuation, the confinement building can be evacuated within 5 minutes. Thererfore the projected dose from the gas release for individuals inside confinement when the fuel failure occurs is:
)
(
min 60 1
min 5
1 25 mrem CDE hr Fraction DAC hr DAC mrem
For the isotopes of interest:
Nuclide Confinement Occupational DAC Immersion Committed Dose Concentration DAC Fraction DAC Equivalent (microCi / cc)
(microCi / cc)
(DAC-hr)
(mrem)
I-131 2.94E-07 2.00E-08 1.47E+01 1.23E+00 3.07E+01 I-132 4.39E-07 3.00E-06 1.46E-01 1.22E-02 3.05E-01 I-133 7.18E-07 1.00E-07 7.18E+00 5.99E-01 1.50E+01 I-134 7.63E-07 2.00E-05 3.82E-02 3.18E-03 7.95E-02 I-135 6.79E-07 7.00E-07 9.70E-01 8.08E-02 2.02E+00 Confinement Internal Dose
49 Therefore, if someone remains in confinement for five minutes, they will receive a committed does to the thyroide of:
Individual CDE = 48 mrem K.
Site Boundary Committed Dose to the Thyroid (CDE)
Section 1.3 of NRC Regulatory Guide 1.145 (November 1982 Rev. 1) indicates that a conservative estimation of site boundary doses can be determined by it assuming that an individual spends two hours at the site boundary immediately following an accident. The deep dose equivalent would be the sum of the dose equivalents from each of the isotopes of interest.
For the isotopes of interest:
Nuclide Confinement Occupational DAC Immersion Committed Dose Concentration DAC Fraction DAC Equivalent (microCi / cc)
(microCi / cc)
(DAC-hr)
(mrem)
I-131 5.18E-10 2.00E-08 2.59E-02 5.18E-02 1.29E+00 I-132 7.72E-10 3.00E-06 2.57E-04 5.14E-04 1.29E-02 I-133 1.26E-09 1.00E-07 1.26E-02 2.53E-02 6.32E-01 I-134 1.34E-09 2.00E-05 6.71E-05 1.34E-04 3.35E-03 I-135 1.19E-09 7.00E-07 1.71E-03 3.41E-03 8.53E-02 Site Boundary Internal Dose Therefore, if someone remains at the site boundary for two hours, they will receive a committed dose to the thyroid of:
Individual CDE = 2.03 mrem L.
Nearest Residence Committed Dose to the Thyroid (CDE)
Again, assuming a two hour occupancy time, the committed dose equivalent to the thyroid at the nearest residency is:
Nuclide Confinement Occupational DAC Immersion Committed Dose Concentration DAC Fraction DAC Equivalent (microCi / cc)
(microCi / cc)
(DAC-hr)
(mrem)
I-131 5.18E-11 2.00E-08 2.59E-03 5.18E-03 1.29E-01 I-132 7.72E-11 3.00E-06 2.57E-05 5.14E-05 1.29E-03 I-133 1.26E-10 1.00E-07 1.26E-03 2.53E-03 6.32E-02 I-134 1.34E-10 2.00E-05 6.71E-06 1.34E-05 3.35E-04 I-135 1.19E-10 7.00E-07 1.71E-04 3.41E-04 8.53E-03 Nearest Residence Internal Dose
50 Therefore, if someone remains at the nearest residence for two hours, they will receive a committed dose to the thyroid of:
Individual CDE = 2.03 X 10-1 mrem M.
Determination of the Committed Effective Dose Equivalent (CEDE)
The committed dose equivalent (CDE) is the cumulative dose that an individual organ in the body would receive due to the uptake of a radioisotope, over a 50 year period. An effective dose equivalent (HE) uses a tissue weighting factor (WT) to equate the risk associated with a non-uniform dose, to the risk associated with a uniformly distributed whole body dose. If the weighting factor is applied to the CDE, then we get the committed effective dose equivalent (CEDE), which provides a measure of what the risk associated with the dose to the individual organ would be if it were evenly distributed in the whole body:
CEDE = (WT)(CDE) 10 CFR 20.2003 Defines the weighting factor (WT) for the thyroid to be:
WThyroid = 0.03 N.
Confinement Committed Effective Dose Equivalent (CEDE)
The committed dose equivalent (CDE) for an individual that is in confinement for five minutes was found to be 48 mrem. Consequently, the CEDE for this individual is:
CEDE = (48 mrem)(0.03) = 1.44 mrem O.
Site boundary Committed Effective Dose Equivalent (CEDE)
The committed dose equivalent (CDE) for an individual that is at the site boundary for two hours was found to be 2.03 mrem. Consequently, the CEDE for this individual is:
CEDE = (2.03 mrem)(0.03) = 6.09 X 10-2 mrem P.
Nearest Residence Committed Effective Dose Equivalent (CEDE)
The committed dose equivalent (CDE) for an individual that is at the nearest residence for two hours was found to be 2.03 X 10-1 mrem. Consequently, the CEDE for this individual is:
CEDE = (2.03 X 10-1 mrem)(0.03) = 6.09 X 10-3 mrem
51 Q.
Determination of the Total Effective Dose Equivalent (TEDE)
The total effective dose equivalent (TEDE) is the sum of the dose due to external sources (DDE) and internal sources (CEDE):
Confinement Total Effective Dose Equivalent (TEDE)
TEDE = 19.1 mrem + 1.44 mrem = 20.5 mrem S.
Site Boundary Total Effective Dose Equivalent (TEDE)
TEDE = 1.58 mrem + 6.09 X 10-2 mrem = 1.64 mrem T.
Nearest Residence Total Effective Dose Equivalent (TEDE)
TEDE = 1.58 X 10-1 mrem + 6.09 X 10-3 mrem = 1.64 X 10-1 mrem
- 16.
Conclusion 10 CFR 20 provides radiation dose limits to radiation workers, and to the general public. For radiation worker, the limits are:
50 rem / yr to an individual organ (CDE) 5 rem / yr whole body (TEDE)
For members of the general public, the limits are:
100 mrem / yr The doses that an individual is predicted to receive due to a fuel failure event in which the core has reached saturation activity is based on the following assumptions:
A.
Individuals inside confinement recognize the problem and evacuate within five minutes.
B.
Individuals that are exposed outside confinement remain at the site boundary for two hours.
C.
Individuals that are exposed outside confinement remain at the nearest residence for two hours.
In all cases, the predicted doses are well below the regulatory limit. A summary of the predicted doses, and any regulatory limit associated with the dose is:
52 Confinement Dose Committed Dose to Thyroid 4.80E+01 mrem 50 rem / yr CEDE 1.44E+00 mrem Immersion 1.91E+01 mrem TEDE 2.05E+01 mrem 5 rem / yr 100 mrem / yr Site Boundary Dose Committed Dose to Thyroid 2.03E+00 mrem 50 rem / yr CEDE 6.08E-02 mrem Immersion 1.58E+00 mrem TEDE 1.65E+00 mrem 5 rem / yr 100 mrem / yr Committed Dose to Thyroid 2.03E-01 mrem CEDE 6.08E-03 mrem Immersion 1.58E-01 mrem TEDE 1.65E-01 mrem 100 mrem / yr Dose Summary Dose Limits Occupational Limits General Public Limits Nearest Residence Dose
53 References
- 1.
Olson, Arne P., Rhode Island Nuclear Science Center Transient Analyses, ANL, Revised Jan. 20, 2016, Page 4.
- 2.
Feldman, Earl E., Loss of Coolant Accident (LOCA) of Safety Analysis Report of the Rhode Island Nuclear Science Center Reactor, August 2010, Page 30.
- 3.
NUREG 1465 Accident Source Terms for Light Water Nuclear Power Plants, USNRC, 1995, Page 2.
- 4.
NUREG 1465 Accident Source Terms for Light Water Nuclear Power Plants, USNRC, 1995, Section 3.2, Page 9.
- 5.
Regulatory Guide 1.183 Alternative Radiological Source Terms for Evaluating Design Basis Accidents at Nuclear Power Reactors, USNRC, July 2000, Section 3.4, Table 5, Page 15.
- 6.
Regulatory Guide 1.183 Alternative Radiological Source Terms for Evaluating Design Basis Accidents at Nuclear Power Reactors, USNRC, July 2000, Section 3.2, Table 1, Page 13.
- 7.
Regulatory Guide 1.183 Alternative Radiological Source Terms for Evaluating Design Basis Accidents at Nuclear Power Reactors, USNRC, July 2000, Section 3.5, Page 15.
- 8.
Regulatory Guide 1.183 Alternative Radiological Source Terms for Evaluating Design Basis Accidents at Nuclear Power Reactors, USNRC, July 2000, Appendix B, Section 1.3, Page B-1.
- 9.
Rhode Island Nuclear Science Center Technical Specifications, February 26, 2016, Specifications 2.2.1.2 and 2.2.2.2, Pages 12 and 14.
- 10.
Regulatory Guide 1.183 Alternative Radiological Source Terms for Evaluating Design Basis Accidents at Nuclear Power Reactors, USNRC, July 2000, Appendix B, Section 2, Page B-2.
- 11.
Regulatory Guide 1.183 Alternative Radiological Source Terms for Evaluating Design Basis Accidents at Nuclear Power Reactors, USNRC, July 2000, Appendix B, Section 3, Page B-1.
- 12.
Regulatory Guide 1.183 Alternative Radiological Source Terms for Evaluating Design Basis Accidents at Nuclear Power Reactors, USNRC, July 2000, Section 3.3, Pages 13 and 14.
54
- 13.
Steady-State Thermal-Hydraulic Analysis for Forced-Convective Flow in the Rhode Island Nuclear Science Center (RINSC) Reactor, Memo to James E. Matos from Earl E. Feldman, Dated September 3, 2010, Section 4.6.8, Page 21.
- 14.
Regulatory Guide 1.183 Alternative Radiological Source Terms for Evaluating Design Basis Accidents at Nuclear Power Reactors, USNRC, July 2000, Appendix B Section 3, Page B-2.
- 15.
Lamarsh, John R., Introduction to Nuclear Engineering, Addison-Wesley Publishing Company, 1977, Page 535.
- 16.
NUREG / CR-2079 PNL-3691 Analysis of Credible Accidents for Argonaut Reactors, USNRC, 1981, Pages 46 -47.
- 17.
EG&G Drawing Number 422873 Rhode Island Nuclear Science Center Test Research Training Reactor Fuel Plate.