Text

Redacted Attachment 5: General Atomics 30441R00022, Rev. B, Reactor-Based Molybdenum-99 Supply System Project, Source Term Analysis Design Calculation Report.
	ML17132A243
Person / Time
Site:	University of Missouri-Columbia
Issue date:	01/24/2017
From:	; General Atomics
To:	; Office of Nuclear Reactor Regulation
Shared Package
ML17132A252	List: ML17132A240; ML17132A241; ML17132A243; ML17132A244; ML17132A245; ML17132A246;
References
	30441R00022, Rev. B
	Download: ML17132A243 (44)
	v • d • e

Attachment 5 RELEASED 30441 R00022 COM Apvd 2017/01/27 Revision B REACTOR-BASED MOLYBDENUM-99 SUPPLY SYSTEM PROJECT SOURCE TERM ANALYSIS DESIGN CALCULATION REPORT Prepared by General Atomics for the U.S. Department of Energy/National Nuclear Security Administration and Nordion Canada Inc.

Cooperative Agreement DE-NA0002773 GA Project 30441 WBS 1110

~ --

nordion

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B REVISION HISTORY Revision Date Description of Changes A 240CT16 Information Release Revised to include only in-pool source terms, incorporate MURR B 24JAN17 comments, and reflect changes to temperature and gap gas POINT OF CONTACT INFORMATION PREPARED BY:

Name Position Email Phone Brian Cluggish Scientist brian.cluggish@ga .com 858-676-7162 Chris Ellis Engineer chris.ellis@ga .com 858-455-3141 APPROVED BY:

Name Position Email Phone B. Schleicher Chief Engineer Bob.Schleicher@ga.com 858-455-4 733 K. Murray Project Manager Katherine .Murray@ga.com 858-455-3272 K. Partain Quality Engineer Katherine.Partain@ga.com 858-455-3225 DESIGN CONTROL SYSTEM DESCRIPTION D R&D DISC QA LEVEL SYS t8J DV&S D DESIGN D T&E N I NIA D NA ii

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B TABLE OF CONTENTS REVISION HISTORY ....................................................................................................................ii POINT OF CONTACT INFORMATION ........................................................................................ ii DESIGN CONTROL SYSTEM DESCRIPTION ............................................................................ ii ACRONYMS ................................................................................................................................. v 1 INTRODUCTION ............................................................................................................... 1 2 APPLICABLE DOCUMENTS ........................................................................................... 2 3 SCOPE AND OBJECTIVES ............................................................................................. 3 4 SOURCE TERMS ............................................................................................................. 4 4.1 Methods ....... .. .. .......... ..... ............................. .. .... .. ...... ...................................... .. ........ 4 4.1.1 Total Nuclide Inventory at the End of Irradiation ......................... .... .. ... ..... ... ... .. 4 4.1.2 Nuclide Inventory in the Hottest Rod at the End of Irradiation ..................... ...... 6 4.1 .3 Post-Irradiation Evolution of the Nuclide Inventory ..... .......... ..... .. .. ... ............ .... . 6 4.1.4 Nuclide Inventory in the Gap Gas of the Hottest Rod ......... ......................... ...... 7 4.1.5 Decay Heat for the Hottest Rod ........... .... ............. ... .. ...................................... 11 4.1.6 Post-Irradiation Total Nuclide Inventory and Decay Heat.. ... .......... ................. 16 4.2 Inputs .. ........... .... ........................................... .... ................... .. ...... ... ... .... ..... .......... .. . 17 4.2.1 Operating Conditions ......... ........................... ............................. ...................... 17 4.2.2 Geometry ..... ....... ........ ... ............ .... ..... ............................................................. 18 4.2.3 Power Density and Temperature Profiles .. .... ........... .. ... ........ .... .... .. .... .......... .. 19 4.3 Results ........... ............. ............ ... ............. .............................. .... .. .... .... .... .... ............ 20 5 MARGINS AND UNCERTAINTIES ................................................................................ 23 6

SUMMARY

...................................................................................................................... 24 7 REFERENCES ................................................................................................................ 25 APPENDIX A - CALCULATION FILES FOR MCNP6/0RIGEN2/MCODE .............................. A-1 APPENDIX B- CALCULATION FILES FOR OTSTM .............................................................. B-1 APPENDIX C - EXAMPLE RELEASE-TO-BIRTH RATIO CALCULATION ........................... C-1 LIST OF FIGURES Figure 1. Simplified diagram of nuclide transport in the SGE facility .. .......... .. ...... .. .... ..... .. ....... ... 2 Figure 2. Power density profile in the hottest rod ........ ................. ... .. ....... .......... .......... .............. 19 Figure 3. Temperature profile in the hottest rod ...... .... ... .. ......... .......... .... ......... .... .... ..... ............. 20 Figure 4. Relative uncertainty in the decay heat versus cooling time for a 3 week irradiation ... 24 iii

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B LIST OF TABLES Table 1: MCODE input parameters .... .. .... ... .. ... ......... ... ......... ...... ....... .............. ........ ... .. ... ............. 5 Table 2: Gamma source spectra (per rod) vs. time since end of irradiation .... ...... .... ........ .. .. ..... 16 Table 3: Input parameters for calculating source terms ..... ........... .... ..... ... .. .... ... ...... .. ..... ...... ... ... 17 Table 4: Summary of activity summed over nuclides .... ........ ..... ...... ......... .... ........... ...... .. .. .. ... .... 21 Table 5: Activity of selected fission gases in the gap gas of the hottest rod ....... .. ..... ..... ....... .. ... 21 Table 6: Summary of decay heat results ............. .... ........ ...... ......... ... .... ...... .... ... ..... .. ...... ....... ... . 22 Table 7: Percent gamma energy deposition versus time after irradiation .. .. ..... .. .. ... ... ......... .... ... 22 Table 8: Decay heat absorption versus time after irradiation ... .. ..... ... .. ....... .. ........ .... ..... ..... .... .... 22 iv

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B ACRONYMS Acronym Description ANS American Nuclear Society ANSI American National Standards Institute GA General Atomics LEU Low Enriched Uranium MATLAB MATrix LABoratory (software by Mathworks, Inc.)

MCNP6 Monte .Q.arlo !i-.E.article code version §.

MCODE MCNP-QRIGEN DEpletion Program MURR University of Missouri Research Reactor NIST National Institute of Standards and Technology Isotope Generation and Depletion Code Matrix Exponential ORIGEN2 Method , version 2 OTSTM Once Through Source Term Model RB-MSS Reactor-Based Molybdenum-99 Supply System SGE Selective Gas Extraction SI Systeme International (international system of units) v

Attachment 5 Source Term Ana lysis Design Calculation Report 30441 R00022/B 1 INTRODUCTION Predicting the generation, decay, and transport of radionuclides in the Reactor-Based Molybdenum-99 Supply System (RB-MSS) is important for defining shielding , cooling , and ventilation requirements. In addition , radiation source terms are needed for analyzing the effects of potential accidents. The purpose of this document is to report the source terms (and the procedures used to calculate them) used in the design and accident analysis of the in-pool portion of the RB-MSS . It is written in accordance with the Quality Assurance Program Document for the RB-MSS project [QAPD-30441-11]. Although the RB-MSS is not a power reactor, the decay heat from the radiation sources is calculated according to ANSl/ANS-5.1-2014 [ANS 2014], and the release fractions for volatile fission products are calculated according to ANSl/ANS-5.4-2011

[ANS 2011 ], since these methods provide conservative estimates of the decay heat and the fission gas release, respectively. NUREG-1537 [NRC 1996] was used as a guide in presenting the source term information .

The RB-MSS generates molybdenum-99 (Mo-99) by irradiating two target assemblies in the MURR reactor pool. Each assembly holds up to 11 rods so that up to 22 rods are irradiated at one time. Each rod consists of a stack of uranium oxide (U02) pellets enclosed in Zircaloy cladding. The initial uranium content of the rods is uranium-238 (U-238) enriched with uranium-235 (U-235).

A simplified diagram of the RB-MSS is shown in Figure 1. The blue arrows show the path of the "uranium stream ," that is the uranium plus any other actinides and fission products that go with it.

The uranium stream moves through the system and changes physical and chemical form but all the actinides stay with the uranium stream . The uranium stream starts off in the target rods in the reactor pool where they undergo a number of irradiation periods . The rods are then allowed to cool for at least 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months before being removed and sent to the hot cells, labeled C1-A, C1-B , and C1-C in Figure 1. The processing of the irradiated rods after removal from the pool will be discussed in an update to this report.

1

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B Cl-A Cl-8 Cl-C Uranium Stream

. - . . - - Rack 1----------=--.~ De-clad Process Vessel Irradiate Cool _ ____,

Temp Storage Pool I

-*I Mo-99 shipment

Te, Rh, Ru , Cu , Zn , Cd , In Exhaust Noble Gas Trap (Kr, Xe)

Mo, Ga, Ge, As , Se, Tc, Sb Figure 1. Simplified diagram of nuclide transport in the SGE facility Radioactive material is generated in the in-pool portion of the RB-MSS when the target rods are irradiated . Most of the radionuclides remain in the U02 pellets, but during irradiation a small fraction diffuse out of the pellets and into the gap between the pellet and the cladding . This report presents calculations of the radioactive source terms used in the design and accident analysis of the in-pool portion of the RB-MSS .

This report is organized as follows. The scope and objectives of this report are given in Section 3.

The assumptions and methods for calculating the source terms , followed by their results, are presented in Section 4. The uncertainty in the source terms is discussed in Section 5. A summary is presented in Section 6. References are provided in Section 7. Description of the input and output files used in the source term calculations are presented in Appendices A and B. A detailed example of the release of iodine-131 (1-131) and krypton-88 (Kr-88) from the U02 pellets is provided in Appendix C.

2 APPLICABLE DOCUMENTS A list of applicable documents is provided below.

2

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B Do cument Number Document Title 30441R00002 MCNP6 (Version 1.0) Verification Test Report 3 0441R00003 ORIGEN2 (Version 2.2) Verification Test Report 3 0441R00006 MCODE (Version 2.2GA) Verification Test Report 3 0441 R00021 Target Assembly Thermal Analysis 3 0441R00023 OTSTM (Version 1.0) Software Verification Test Report Mo-99 Target Assembly Nuclear Design for Once-Through 30441R00031 Operation OTSTM (Version 1.0) Software Design Requirements, 30441R00034 Description , and Verification Report 3 0441S00001 Molybdenum-99 Supply System Requirements Document Quality Assurance Program Document - Phase 11 , Reactor-QAPD-30441-11 Based Molybdenum 99 Supply System (RB-MSS )

3 SCOPE AND OBJECTIVES The objective of this report is to present source terms for use in designing the RB-MSS . Source terms are needed for determining the radiation dose received by personnel working near the pool and shielding requirements for the transfer casks used to remove the rods from the pool. In addition, the heat load from radioactive decay after the end of irradiation is needed for designing cooling systems for the rods in the pool.

Source terms are also needed for performing accident analysis of the RB-MSS , such as the amount of radionuclides that would be released if a rod were to break open.

To meet the needs for design and accident analysis , the calculation of the following source terms are presented in this report:

The total radionuclide inventory and decay heat load in a full "set" of 22 target rods

The nuclide inventory and decay heat load in the hottest of the 22 rods

The nuclide inventory in the "gap gas" of the hottest rod, volatile fission products released by the U02 pellets into the pellet-cladding gap during irradiation These source terms are calculated from 0 hours0 days 0 hours 0 weeks 0 months after the end of irradiation, when the activity and decay heat is highest, to 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months after the end of irradiation, just before the rods are removed from the pool.

3

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B 4 SOURCE TERMS This section presents calculation methods of the RB-MSS source terms followed by their results.

Section 4.1 discusses the methods used, Section 4.2 discusses the inputs to the calculation, and Section 4.3 discusses the results .

For all source terms the activity in Curies of nuclide i , Acti , is defined by the following equation :

kN*

A . = 1 1 ct1 - 3.7 x 1010 Here Ni and il.i are the number of atoms and the decay constant of nuclide i , and 3.7 x 10 10 is the number of decays per second per Curie.

4.1 Methods Calculation of the total nuclide inventory for 22 rods at the end of irradiation was performed with the assistance of MCODE [Xu 2006, 30441 R00006] , an open source linkage-program for combining MCNP6 [LANL 2013 , 30441 R00002] and ORIGEN2 [Ludwig 2002 , 30441 R00003].

The MCODE results were then used as the input for the Once Through Source Term Model (OTSTM) [30441 R00023, 30441 R00034], which was used to calculate post-irradiation decay of the nuclide inventory in the 22 rod set, the hottest rod , and the gap gas, as well as the decay heats .

All equations in this section are written assuming Systeme International (SI) units, unless otherwise noted.

4.1 .1 Total Nuclide Inventory at the End of Irradiation Calculation of the nuclide inventory in a set of target rods at discharge (the end of the final irradiation period) was performed using MCODE to coordinate MCNP6 and ORIGEN2. A list of the libraries used is provided in Table 1. Tables of the calculation files used are presented in Appendix A.

At the start of this coupled simulation, the thermal power in the target during irradiation , Pset

was held constant and the neutron flux was assumed to be uniform over the entire volume of U02 target material. The following steps were alternated by MCODE to simulate a sequence of depletion (irradiation) and decay (shut down) periods:

Depletion step: MCODE, after a series of calls to MCNP6 and ORIGEN2, updates the ORIGEN2 cross section library using the MCNP6 calculated one-group cross sections.

The MCNP6 cross section library closest to the average temperature in the target rods was used. The updated cross sections apply to nuclides contributing to at least 99.9% of 4

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B all neutron absorptions as calculated by ORIGEN2 . In this method the transmutation of isotopes are consistent between MCNP6 and ORIGEN2. Additionally, MCODE iterates on the ORIGEN2 input flux in multiple calls to achieve the desired power. After the cross section library is updated and input flux is converged, MCODE calls ORIGEN2 again to perform a depletion over the time step using its full burnup chain. The depletion was for a single MCNP6 material representing the entire U02 volume within two target assemblies (22 target rods).

Decay step: Using the target material composition at the end of the previous depletion time step , MCODE calls ORIGEN2 to perform the decay calculation using its decay library.

MCODE then updates the target material composition in a new MCNP6 input for the subsequent time step.

Neutron irradiation of the target rods results in the production of additional nuclides due to fission and neutron absorption . The nuclide inventory further evolves with time due to nuclear decay, as well as additional fission and neutron absorption during the depletion steps.

See Table 1 for a list of the time steps modeled and more details on the MCODE input parameters .

Note that this method resets the total target power to its fresh value at the beginning of each new depletion step , when in reality the power would slowly decrease with burn up. From a source term perspective , this method will yield conservative results. The MCNP6 geometry modeled the MURR with maximum burnup driver fuel which yielded the maximum target assembly power.

Table 1: MCODE input parameters Parameter Value ORIGEN2 decay library decay.lib ORIGEN2 gamma library gxuo2brm.lib ORIGEN2 cross section library pwru.lib MCNP6 cross section library ENDF/B-Vll.1 neutron absorption fraction threshold 0.999 Depletion Specification step 1 (depletion) step 2 (decay) step 3 (depletion)

~

step 4 (decay) LJ step 5 (depletion) 5

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B 4.1.2 Nuclide Inventory in the Hottest Rod at the End of Irradiation The target rods are all initially identical and the number of fission events in a set of target rods is proportional to the fission power in the set. Therefore, the nuclide inventory in the hottest target rod at the end of irradiation is obtained from the nuclide inventory of the set by using the following equation:

Here t 0 and tEoI are the time at the start of the first irradiation and at the end of the last, respectively; Ni(rod, t) and Ni(set, t) are the number of atoms of nuclide i in the hottest rod and in the full set, respectively, at a time t; and Nrod is the number of rods in the set. Note that Ni(set, t 0 ) is equal to zero except for U-238, U-235, and 0-16. In this way , the nuclide inventory in the hottest rod is calculated from 0 to 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months after the end of irradiation.

4.1.3 Post-Irradiation Evolution of the Nuclide Inventory The decay of the nuclides into daughter nuclides is calculated by OTSTM. This MATLAB-based program solves a generalized Bateman equation [Bateman 191 O] successively for each nuclide in each decay chain, starting with the top level nuclide and working down through the generations of daughter nuclides. The general equation is:

dN*1 dt = -A *N* +

I I I AJ.

J)J -+1*N*]

j Here Ni is the number of atoms of nuclide i at time t, Ai is the decay rate of nuclide i , fj-+i is the fraction of nuclides j that decay into nuclide i, and the last term is the source of nuclide i due to the radioactive decay of parent nuclides.

An analytical solution to this equation is found as follows starting with the top-level nuclide in the decay chain, which has no parent species. For that nuclide, the solution is:

N 0 (t) = N0 (0) exp(-A. 0 t)

Here the subscript "O" refers to the top-level parent nuclide, and N0 (0) is its number of atoms at time 0. For each succeeding generation the solution for the parent species is inserted in the equation for the daughter species, which is then solved analytically. Writing the general solution for a parent nuclide as a sum of exponential decays of nuclides k:

Nj = Lk ajk exp(-Akt) 6

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B where the ajk are constants, then the general solution for the daughter species is:

Ni = NJO)T/i + [L' A.jfj-+i L

' ajk T/k-T/i]

A* _A.

j k l k T/i = =

exp(-A.it) , T/k exp(-A.kt)

Here j denotes the parents of nuclide i. Thus, the solution for the top-level nuclides can be used to find the solutions for their daughter nuclides, which can then in turn be used to find the solutions for their daughters, and so on for all succeeding generations.

Note that the solution does not diverge if A.i = A.k because:

This is accounted for numerically when calculating the number densities of the daughter species.

The decay rates and branching fraction data for each nuclide were obtained from the ORIGEN2 decay library. All 1307 nuclides in the ORIGEN2 decay library are tracked . The input and output files used by OTSTM for these calculations are discussed in Appendix B. OTSTM has been verified to have excellent agreement with ORIGEN2 in modeling radioactive decay [Cluggishb] .

4.1.4 Nuclide Inventory in the Gap Gas of the Hottest Rod At the high temperatures experienced by the rod during irradiation, almost all of the fission products are volatile and can diffuse out of the U02 pellets. The nuclide inventory in the gap gas of the hottest rod at the end of irradiation (t = O) is calculated as :

Here Ni(gap, t) is the number of atoms of nuclide i in the gap gas at a time t after the end of the final irradiation, and (i) i is the "release to birth ratio" for nuclide i, (i.e. the fraction of atoms of nuclide i that are released from the pellets). The release-to-birth ratios are calculated according to ANSI standard ANSl/ANS-5.4-2011 as described in Section 4 .1.4.1.

The nuclide inventories in the gap gas at times after the end of irradiation (t > 0) are calculated by solving the nuclide decay equations as described in Section 4.1 .3. It is assumed that no further nuclides are released after the end of irradiation since the rod is much cooler and fission is no longer occurring .

7

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B Since the rods are only irradiated in the pool, the gap gas can only be released into the pool or after at least 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months of cooling . Therefore , it is assumed that only the volatile halogen and noble gas fission products are released when the cladding is breached. Although the non-volatile nuclides in the gap are not released , they do contribute to the inventory of volatile nuclides through decay.

4.1.4.1 Release to Birth Ratios ANSl/ANS -5.4-2011 provides a method for calculating the release to birth ratio , R/B , for each nuclide. The release-to-birth ratio for a nuclide is the fraction of atoms of that nuclide that escape from the U02 matrix during irradiation . For fission products with half-lives less than 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months the release-to-birth ratio for the ith nuclide is calculated using equation 1 of Reference [ANS 2011 ],

reproduced here:

H ere~ is the surface area-to-volume ratio of the U02 matrix, and ai, Di, and ili are the precursor coefficient, diffusion coefficient, and decay rate, respectively, for nuclide i. The precursor coefficients are given in Table 1 of Reference [ANS 2011 ]. The diffusion coefficient in units of m2/s depends on the temperature and fission rate density in the U02 as follows [ANS 2011]:

Here T and P are the highest temperature (in Kelvin) and fission rate density (in fissions/m 3/s) reached in the U02 during irradiation. The fission rate density is related to the power density, h, by:

. h F=-

Er Here Er is the recoverable energy per fission event. Reference [ANS 2011] specifies that Er =

1

- - -- J = 200 MeV.

l .0 Xl0 10 7t The surface to volume ratio depends on the temperature as follows [ANS 2011 ]:

~ = {120 cm- 1 T :::::; Tunk V 650 cm- 1 T > Tlink 8

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B Here Tunk is temperature at which the fission gas bubbles are assumed to become interlinked on the U02 grain boundaries allowing the gas to be released. It depends on the burnup as follows

[ANS 2011]:

9800

_ I ( B ) + 273 Bu::; 18.2 MWd/kgU Tlink - n 17 6 u

{

1434 - 12.85Bu + 273 Bu> 18.2 MWd/kgU Here Bu is the burnup in megawatt days per kg of initial uranium . It is calculated as follows :

Bu (MWd) = h (~) trun(days) kgU Mu/rrR~l Here the power density is specified in MW/m 3 , trunCdays) is the total irradiation time in days, Mu is the initial mass of uranium in a rod in kg , RP is the radius of the pellets, and L is the total height of the stack of pellets in the rod.

For nuclides with half-lives between 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months and 60 days the release-to-birth ratio is calculated using equation 3 of [ANS 2011 ], reproduced here:

Here aasmKr and ilasmKr are the precursor coefficient and decay rate , respectively for metastable Kr-85 .

Longer lived nuclides are not covered by ANSl/ANS-5.4-2011 . Their release-to-birth ratios are estimated using the following formula:

(!!..) = 3 (co th ( y'µ) _ .!.) , µ=max ( ili,-1-) 9 2 Bi {µ µ trun a *D* (~)

I I V Here coth is the hyperbolic cotangent function and µ is a dimensionless constant. This is the solution to the Booth equation given in equation 3 in the Foreword of [ANS 2011] , except that for long lived nuclides ili has been replaced by l/trun* since the escape of the nuclides from the pellet is limited by the total run time rather than radioactive decay.

Since the diffusion coefficient and surface-to-volume ratio at a given point in a target rod depends on the temperature , T , and the fission rate density, P, at that point, the average release-to-birth 9

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B ratio for a rod must be calculated by integrating the local, power density weighted , release-to-birth ratio over the rod and then dividing by the total power in the rod :

R 2rr f0Rp rdr f0L dz h(r, z) (~) . (r , z)

(fj)i = R L l 2rr f0 P rdr f0 dz h(r, z)

Here r and z are the radial and axial coordinates in the rod , respectively, and Rp and L are the radius of the pellets and the length of the rod, respectively. The power density, h(r , z) , is included in the integral because the local rate of production of fission products is proportional to the local power density.

The temperature profile in a target rod is calculated as described in 30441 R00021. The heat equation is solved in steady state using the power density calculated by ORIGEN2, the temperature dependent thermal conductivity of U02, and the temperature dependent heat flux in the pellet-cladding gap.

The procedure used to calculate the release-to-birth ratios in the hottest rod at the end of irradiation is as follows:

1. A cylindrical grid is overlaid on the rod, with Nr equally spaced nodes in the radial direction and Nz equally spaced nodes in the axial direction.

2. The temperature at each grid node, T(r, z) is calculated as described in 30441 R00021 .

3. The fission rate density, F(r, z) and burnup, Bu (r , z) are calculated at each point on the grid.

4. The release-to-birth ratio for fission product i at each grid node, (~) i (r, z) is calculated for each fission product using the procedure set out in ANSl/ANS-5.4-2011. The local burnup and temperature are used to obtain the interlinkage temperature and then the surface-to-volume ratio of the U02. The local temperature and fission rate density are used to calculate the diffusion coefficient.

5. A factor of 5 conservatism (recommended in Section 8 of ANSl/ANS-5.4-2011) is applied to (R/B) at each grid node.

6. The rod averaged release fraction is then obtained by integrating over the volume of the rod. The integral is performed using trapezoidal integration.

10

Attachment 5 Source Term Analysis Design Calcu lation Report 30441 R00022/B A detailed example calculation of the net release-to-birth ratios for 1-131 and Kr-88 is presented in APPENDIX C.

4.1.5 Decay Heat for the Hottest Rod The decay heat for the hottest rod is calculated by OTSTM as set out in standard ANSl/ANS-5.1-2014 [ANS 2014). ANSl/ANS-5 .1-2014 defines 4 contributions to the decay heat:

1. Pd, the base decay heat from fission products

2. Pde* the decay heat correction due to neutron capture by fission products during irradiation

3. PdHE

the decay heat due to the heavy element nuclides U-239 and Np-239

4. PdA* the decay heat due to other actinides Items 1 and 3 are proportional to the fission power, whereas items 2 and 4 increase slightly more rapidly with increasing fission power than items 1 and 3.

The MCODE/MCNP6/0RIGEN2 calculations discussed in Section 4.1.1 indicate that over 99%

of the fission events during a 3 week irradiation are due to U-235 fissions . Therefore , for purposes of calculating the decay heat, we assume that all the fission power comes from U-235 fission . As recommended in Section C of Reference [ANS 2014) the energy re leased per fission event is set at Er = 200 MeV.

The decay heat is calculated as described in the standard ; no alternate models are used . Also, to be conservative , the expression for the neutron flux and neutron absorption cross sections recommended in the standard are used , rather than those calculated with MCNP6 for the calculations described in Section 4.1.1. Because the irradiation time for each batch of rods is fairly short and the removal of Mo-99 from the pellets proceeds quickly, the anticipated operating conditions for the RB-MSS are well within the limits given in Reference [ANS 2014). These limits include: less than 3 fissions per initial fissile atom, less than 4 years of operation , and less than 1010 seconds of cooling time after irradiation.

Combining equations 2 and 3 of Reference [ANS 2014) gives the following formula for the base fission product decay heat:

Nirr 23 Pd(t) = P~od LL~~ exp(-Ajtci ) (1 - exp(-Ajtirr,i)]

f i=1 j=l J N irr tci =t + L (tirr,k + torf.k) k=i+l 11

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B Here Prod is the fission power in the hottest rod , tci is the cooling time after the ith irradiation, Nirr is the number of irradiations, t irr,i is the length of the ith irradiation , toff,i is the time the reactor is shut down before the ith irradiations, t is the cooling time after the final irradiation, and the aj and Aj are fitting coefficients to empirical decay heat data. The values of the aj and A.j are given in Table 5 of Reference [ANS 2014). Note that there is no sum over fissionable nuclides because we assume that 100% of the fission power comes from U-235 fission . The uncertainty in Pd_ is calculated using the same formula as for Pd_ itself, except that the factors of aj are replaced with factors of pj, which are fitting coefficients given in Table 5 of Reference [ANS 2014]. The other terms in the decay heat calculation are defined such that they do not contribute additional uncertainty.

For cooling times less than 10,000 seconds after the final irradiation the fission product decay heat due to neutron capture by fission products during irradiation is given by equation 13 of Reference [ANS 2014], reproduced here:

The factor G(t), which accounts for neutron capture in fission products , is given by equation 14 of Reference [ANS 2014], reproduced here:

G(t ) = (3 .24 x 10- 6 + 5.23 x 10- 10 t)r 0 .4tjJ T =L(tirr,i + toff,i) i Here r is the total time from the beginning of the first irradiation to the end of the last and tjJ is the number of fissions per U-235 atom, which is calculated as:

Here F is the fission rate density defined in Section 4.1.4.1, trun is the total irradiation time , mu 235 is the atomic mass of U-235, and E is the initial U-235 enrichment in percent. This equation gives a value of tjJ about 1% higher than the formula provided in equation C.6 of Reference [ANS 2014].

For cooling times greater than 10,000 seconds after the final irradiation the fission product decay heat due to neutron capture in fission products is given by equation 22 of Reference [ANS 2014],

reproduced here:

12

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B Here Pdcs is the correction due to neutron capture by Cs-133 , and PdE is the contribution to the decay heat from neutron capture for the remaining fission products . Combining equations 15-20 of Reference [ANS 2014] gives the following expression for Pdcs:

Prod { 1 - a4 a3 - a4 }

Pdcs(t) = -E f Y3E4 A. 4 exp(-A.4t) A.

4 + 0'4

+ A.

0'3 - 4 - 0'4

a3 = exp(-a3 <1>Trun)

U4 =: exp(-A4Trun -a4<l>Trun)

Here Y3 is the effective cumulative fission yield of Cs-133, £ 4 is the energy per decay for Cs-134, A. 4 is the decay constant for Cs-134, a3 and a 4 are the neutron absorption cross sections for Cs-133 and Cs-134 , respectively. Values for these quantities are given in Section 4.4.2 of Reference

[ANS 2014]. Also , is the effective neutron flux during irradiation . The effective neutron flux is calculated using equation 18 of Reference [ANS 2014], reproduced here:

<1> = 104 (z.ss x 10 13 _!___)

Eef f E

Eetf =: z+ 1 3

10 ( h )

S =: 10 6 Mu/(nR~ L)

Here Eetf is the effective enrichment (as defined in Section 4.4.2 of ANSl/ANS-5 .1 -2014) , Sis the specific power density in units of megawatts per ton of initial uranium and h is the power density.

10 3 The factor of in the formula for S converts it from Watts per kilogram to megawatts per ton ,

10 6 and the factor of 10 4 in the formula for converts it from n/cm 2/s to n/m 2/s. Note that this approximation gives a value of the neutron flux 2-3 times higher than the more detailed modeling of the RB-MSS by MCODE and MCNP6. Nonetheless, the ANSl/ANS-5 .1-2014 method is used in calculating the decay heat because it is more conservative .

The contribution to the decay heat from neutron capture for the remaining fission products (other than neutron capture by Cs-133) is calculated using equation 21 of Reference [ANS 2014],

reproduced here:

Here His a constant interpolated from the values given in Table 13 of Reference [ANS 2014].

The decay heat of U-239 and Np-239 for one irradiation period is calculated using equations 23 through 25 of Reference [ANS 2014]. Combining those equations and summing over irradiation periods gives the following expression for PdHE :

13

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B Here, Eu and ENp are the average energies from the decay of one atom of U-239 and Np-239, respectively; and ilu and ii.Np are the decay constants for U-239 and Np-239, respectively. Values for these quantities are given in Section 5.1 of Reference [ANS 2014). The parameter R is the atoms of U-239 produced by U-238 capture per fission evaluated at the time of shutdown, and is calculated using the following equation :

0"233<1> Er 1 - t:/100)

R = n 238 - - .- = 10 ( 2.58 x 1013 a 238 - ----

F mu Eerr Here n 238 is the number density of U-238 atoms , mu 238 is the atomic mass of U-238, and a 238 is the U-238 atomic cross section for neutron capture. A value of a238 = 2.7 barns is given in Table C.2 of Reference [ANS 2014). This value is over twice the cross-section calculated by MCNP6 for the source term calculations discussed in Section 4.1 .1; nevertheless, the larger value from ANSl/ANS-5.1-2014 was used in order to be conservative . Using these equations with the expected values of Er= 200 MeV and£ = 19.75% gives R = 0.42. Note that both and Pare proportional to the power density, but R is independent of power density.

The decay heat of the other actinides is calculated using equation 26 of Reference [ANS 2014),

reproduced here:

, [10- 4 A(t)Bua(t) ]

Pd.A(t) = Pd(t) -----

Eerr a(t) = 1.5 exp (- 10t9 5 )

S ) trun Bu =( 1000 tday Here Bu is the uranium burnup in units of gigawatt days per ton of initial uranium (with S in units of megawatts per ton), tday is the number of seconds in a day, and A is obtained by interpolating from the values given in Table 14 of Reference [ANS 2014).

The procedure used to calculate the decay heat in one rod irradiated Nirr times at a time t after the final irradiation is as follows:

14

Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B

1. The uncorrected fission product decay heat is calculated for the each irradiation.
2. The decay heats for the individual irradiations are added up to get Pd_ , the uncorrected, fission product, decay heat for the Nirr irradiations.
3. A cylindrical grid is overlaid on the rod, with Nr equally spaced nodes in the radial direction and Nz equally spaced nodes in the axial direction.
4. The specific power density, S, neutron flux , <JJ , number of neutron captures per fission , R, number of fissions per U-235 atom, t/J, and burnup, Bu, are calculated at each grid node as shown above.
5. The neutron absorption correction to the decay heat density, hde* and the decay heat densities from U-239 and Np-239, hdHE* and other actinides, hdA* are calculated at each node on the grid .
6. The total decay heat in the rod , Pd 1 , at time tis obtained by integrating the corrected decay heat densities over the volume of the rod :
Pde= 2rr Rp rdr JL dzhde(r,z)
J0 0 Here Pde* PdHE
and PdA are the corrections to the decay heat due to neutron capture by fission products, Np-239 and U-239, and the other actinides , as defined in Reference [ANS 2014]. The integrals are performed numerically using trapezoidal integration .
4.1.5.1 Correction due to gamma ray escape During in-pool cooling, not all of the decay heat in the target rods is absorbed in the U02. Some of the gamma rays escape from the rods and are absorbed in the surrounding water.
The decay heat absorbed in a target rod, Pabs* is given by the following equation .
Pabs = [1 - fy + fyfabs]Pdl Here fr is the fraction of the decay heat due to gamma emission, fab s is the fraction of decay heat due to gamma emission that is absorbed in the rod, and Pd 1 is the decay heat in one rod.
15
Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B The fy parameter was obtained from the ORIGEN2 output files produced from the calculations described in Section 4.1.1. Obtaining fabs also requ ires data from the ORIGEN2 output files , in addition to a separate MCNP6 model. In this model , calculations were performed using the F6 (energy deposition) tally. The geometry modeled was a single target rod centered within a 1 meter radius sphere of water. Only photons (gammas) were tracked in the model. The nuclide inventory in the rod was obtained at multiple time steps by taking the ORIGEN2 output concentrations and dividing by 22 to get the value for a single rod . The gamma ray source definition (the SDEF card) was also obtained from the ORIGEN2 output files . The spectra of the gamma ray source at 5 different times after the end of irradiation are provided in Table 2. The gamma ray source includes the contribution of fission products plus actinides and daughters. The MCNP6 cross section library used the 2012 continuous-energy photoatomic data from ENDF/B-Vl.8 (library name mcplib84 ). It was assumed that 100% of the gammas originate in the U02 - target clad activation was neglected and assumed to be small.
Table 2: Gamma source spectra (per rod) vs. time since end of irradiation Mean Energy Gamma Release Rates (gammas/sec)
(MeV) 0 hour0 days 0 hours 0 weeks 0 months 0.5 hour5.787037e-5 days 0.00139 hours 8.267196e-6 weeks 1.9025e-6 months 1 hour1.157407e-5 days 2.777778e-4 hours 1.653439e-6 weeks 3.805e-7 months 3 hours3.472222e-5 days 8.333333e-4 hours 4.960317e-6 weeks 1.1415e-6 months 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months 1.000E-02 1.649E+15 5.981E+14 4.861 E+14 3.455E+14 2.719E+14 2.500E-02 3.989E+14 1.557E+14 1.296E+14 9.349E+13 7.601 E+13 3.750E-02 3.278E+14 1.447E+14 1.261E+14 1.020E+14 8.905E+13 5.750E-02 3.468E+14 1.149E+14 9.432E+13 6.657E+13 5.152E+13 8.500E-02 2.823E+14 1.113E+14 8.957E+13 6.406E+13 5.376E+13 1.250E-01 2.724E+14 1.397E+14 1.242E+14 1.014E+14 9.215E+13 2.250E-01 6.181E+14 2.394E+14 1.963E+14 1.473E+14 1.274E+14 3.750E-01 3.793E+14 1.442E+14 1.006E+14 5.835E+13 4.554E+13 5.750E-01 6.276E+14 3.606E+14 3.162E+14 2.451 E+14 2.077E+14 8.500E-01 6.514E+14 3.518E+14 2.899E+14 1.782E+14 1.311 E+14 1.250E+OO 4.476E+14 2.162E+14 1.640E+14 9.009E+13 5.718E+13 1.750E+OO 1.440E+14 8.023E+13 6.800E+13 4.886E+13 3.878E+13 2.250E+OO 8.936E+13 4.791E+13 3.543E+13 1.615E+13 7.195E+12 2.750E+OO 3.778E+13 1.771E+13 1.286E+13 5.005E+12 2.078E+12 3.500E+OO 2.325E+13 4.636E+12 3.350E+12 1.193E+12 3.304E+11 5.000E+OO 1.258E+13 7.273E+10 4.964E+10 2.455E+10 1.149E+10 7.000E+OO 1.016E+11 4.132E-02 9.586E-03 9.586E-03 9.591E-03 9.500E+OO 2.001E+07 1.091 E-03 1.091E-03 1.091 E-03 1.091 E-03 4.1.6 Post-Irradiation Total Nuclide Inventory and Decay Heat The post-irradiation nuclide inventory of the 22 rod set is calculated by scaling from that of the hottest rod by inverting the equation in Section 4.1.2, which gives:
16
Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B Here t ;:::: tEoI is the post-irradiation time at which the nuclide inventory of the set is to be calculated.
The decay heat is proportional to the fission power. Since the pre-irradiation decay heat is negligible, the total decay heat for the 22 rod target, Pdset* is obtained by multiplying the decay heat generated in the hottest rod by the ratio of the total fission power to the power in the hottest rod :
4.2 Inputs This section describes the numeric inputs used in calculating the source terms. The input parameters are summarized in Table 3 and further described in the subsequent sub-sections.
Table 3: Input parameters for calculating source terms Parameter t--~~~~~~~~~~~~~~~~~~~~~~~~~+
Value Fission ower for a set of 22 rods, P
~
Down time between irradiations, t 0
~
Number of weeks of irradiation, N*
Pellet radius, R 2.5mm Rod len th, L 60 cm Initial mass of uranium in 1 rod 0.1062 k U-235 enrichment, E 19.75wt. %
Number of radial nodes in rod rid, N 25 Number of axial nodes in rod rid, N 51 4.2.1 Operating Conditions To be conservative, the operating conditions that would produce the highest radioactivity and decay heat terms were used in calculating the source terms. The linear fission power density, Pi.
profiles in each of the 22 rods in the two target assemblies, as well as the total fission power in each rod, Prod
and in the two target assemblies , Pset* has been calculated over the complete range of expected operating conditions for the RB-MSS [30441 R00031]. This data was used in choosing the following conservative operating conditions:
17
Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B
1. The source term activity of the target rods at the end of irradiation is approximately proportional to the burnup, so the maximum value of the fission power was conservatively used . For two target assemblies containing 22 target rods , the maximum expected fission power is calculated to be -
with the power dropping slowly with time ~
[30441 R00031]. To be conservative and account for uncertainty, two standard deviations
(- ) were added to this maximum power and the drop in power with time was LJ neglected. Therefore, the source term was calculated assuming a constant fission power of Pset =- during each irradiation.
2. To maximize the gap gas activity, the linear power density profile from the operating conditions that give the highest expected linear power density was used as the linear power density profile for the hottest rod . This produces the highest temperatures and therefore the highest release-to-birth ratios. Including two standard deviations of uncertainty, the expected peak linear power is PL = - and the power in the hottest ~
rod is Prod =- [30441 R00031 ]. Note that these values of PL and Prod are higher than the values in the hottest rod at the highest maximum total power and so are doubly LJ conservative. Furthermore , the temperature profile was calculated under the assumption that the cooling water flow was at 85% of the nominal rate for the entire which gives higher temperatures in the rod .
3. Under normal operation , up to 22 target rods are expected to be irradiated for - at up to - per week. However, occasionally a batch of rods may be irradiated for
- [30441 S00001 ]. Therefore , to be conservative, the source terms were calculated ~
for a 22 rod set that undergoes Nirr = I irradiations, each t irr = . . . - in length and separated by shut down periods of duration t 011 = ~* The total irradiation time is LJ trun=****
4.2.2 Geometry All of the target rods are identical. The U02 pellets in the rods are 95% of theoretical density and have a radius RP = 2.5 mm . Each pellet is fabricated from uranium enriched to= 19.75 U-235 wt. %. The stack of pellets in each rod is of length L = 0.6 m. The total initial mass of uranium is Mu = 106.2 g per rod. The pellets are fabricated with dishes and chamfers on the ends to allow for thermal expansion, so that the mass of uranium is about 2% less than would be expected if the pellets were ideal cylinders.
The effect of variations in the density and enrichment of the pellets is accounted for in the -
~
of uncertainty in the target power. Furthermore , since the target power is fixed at -
changes in the density and enrichment have a negligible effect on the calculated source terms.
LJ 18
Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B For the gap gas and decay heat calculations a numerical grid is imposed on the stack of pellets.
The grid has Nr = 25 equally spaced radial grid points and Nz = 51 equally spaced axial grid points.
4.2.3 Power Density and Temperature Profiles The power density profile used for the hottest rod is shown in Figure 2 [30441 R00031, 30441 R00021 ]. The power density is essentially proportional to the neutron flux.
Figure 2. Power density profile in the hottest rod The power density is slightly lower along the centerline of the rod due to the neutron self-shielding effect.
The temperature profile in the target rod is needed for calculating the fission gas release. The temperature profile for the hottest rod is shown in Figure 3 [30441 R00021 ]. The temperature is highest where the linear power density is highest.
19
Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B 2000
.-.. 1500
(.)
I-1000 500 0
60 r (cm) 0 z (cm)
Figure 3. Temperature profile in the hottest rod The volumetric-average temperature in the hottest rod is 1224 K. For this reason , the MCNP6 1200 K cross section library was used in calculating the total nuclide inventory. Cross section temperature dependence becomes less significant at temperatures above 1200 K, therefore the cross sections are appropriate for source term calculations during irradiation.
4.3 Results The total activity and decay heat in a set of 22 rods and the hottest rod are shown in Table 4 at 0 to 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months after the end of irradiation. The total activity of the set decreases by a factor of 5 with 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months of cooling , and the activity of the set is 21.45 times larger than the hottest rod.
20
Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B Table 4: Summary of activity summed over nuclides Time (hours) 22 Rod Total (Ci) Hottest Rod (Ci) 0 2.542E+06 1.185E+05 1 8.106E+05 3.779E+04 3 6.291 E+05 2.933E+04 6 5.292E+05 2.467E+04 Tables of the activity per nuclide for the 22 rod set, the hottest rod , and the halogens and noble gases in the gap gas are provided in the spreadsheet that accompanies this document [Cluggish 2016c].
The activity of the most important nuclides in the gap gas of the hottest rod is shown in Table 5, with the total activity of all halogen and noble gas isotopes shown in the last row. Note that the activities of 1-132, Kr-85 , and Xe-135 are higher at 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months than at 0 hours0 days 0 hours 0 weeks 0 months due to the decay of Te-132 , Kr-85m , and 1-135, respectively. Also, the sum of the individual activities in Table 5 is slightly lower than the total Br, I, Kr, and Xe activity shown in Table 5, because some minor, short lived nuclides are not included in Table 5. The activity of the gap gas is much smaller than the whole rod, but decreases relatively slowly due to the production of 1-132 by the decay of Te-132 .
Table 5: Activity of selected fission gases in the gap gas of the hottest rod Fission Product Half-life RIB Activit Ci at 0 hrs Activit Ci at 6 hrs 1- 131 8.041 d 26.2934 25 .8531 1-132 2.3 h 144.6951 84.4133 1-133 20 .8 h 45 .7961 37 .9656 1-134 52.6 m 26 .0561 0.7470 1-135 6.611 h 30 .5511 16.2867 Kr-85 10.72 0.0316 0.0316 Kr-85m 4.481 h 5.9623 2.3670
~
Kr-87 1.272 h 6.2891 0.2400 Kr-88 2.839 h 12.0723 2.7903 Kr-89 Kr-90 3.17 m 32.32 s 2.2662 0.8907 0.0000 0.0000 Ll Xe-133 5.245 d 64.8080 64.0750 Xe-133m 2.19 d 1.5804 1.5520 Xe-135 9.089 h 8.1066 13.3210 Xe-135m 15.29 m 5.5078 2.6087 Xe-137 3.83 m 3.0510 0.0000 Xe-138 14.17 m 5.9107 0.0000 Xe-139 39.5 s 1.0359 0.0000 Total Br, I, Kr, Xe 417.6132 256.0431 21 I
__ _J
Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B The total decay heat of the set of 22 rods and that of the hottest rod are shown in Table 6 at 0 to 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months after the end of irradiation. The decay heat decreases by a factor of 10 during 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months of cooling.
Table 6: Summary of decay heat results Time (hours) 22 Rod Total (W) Hottest Rod (W) 0 34469.2 1607.0 1 6596.8 307 .5 3 4538.5 211 .6 6 3513.5 163.8 The fractions of the gamma heat absorbed in the U02, cladding , and cooling water are provided in Table 7. The table shows that of the total gamma energy generated in the U02 about half is deposited in the target rod .
Table 7: Percent gamma energy deposition versus time after irradiation Material 1 hour1.157407e-5 days 2.777778e-4 hours 1.653439e-6 weeks 3.805e-7 months 3 hours3.472222e-5 days 8.333333e-4 hours 4.960317e-6 weeks 1.1415e-6 months 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months U02 41 .7% 44.0% 46.2%
Cladding 7.0% 6.9% 6.8%
Surrounding Water 50.8% 49.0% 46.9%
The fraction of the decay heat that comes from gammas, and the fraction of the total decay heat that is absorbed in the U02 are shown in Table 8. The symbols are defined in Section 4.1.5.1 .
Table 8: Decay heat absorption versus time after irradiation Quantity 1 hour1.157407e-5 days 2.777778e-4 hours 1.653439e-6 weeks 3.805e-7 months 3 hours3.472222e-5 days 8.333333e-4 hours 4.960317e-6 weeks 1.1415e-6 months 6 hours6.944444e-5 days 0.00167 hours 9.920635e-6 weeks 2.283e-6 months Fraction gamma power absorbed by U02, fabs 0.417 0.440 0.462 Fraction decay heat due to gammas, fv 0.588 0.576 0.568 Fraction total decay heat absorbed by U02, 1 - fv + fvfabs 0.657 0.677 0.694 Decay heat generated by the hottest rod (W), Pd 1 308 212 164 Decay heat absorbed by U02 in hottest rod (W), Pabs 202 143 114 This table shows that about two-thirds of the decay heat in the rods is absorbed in the U02. The rest is absorbed in the cladding or the surrounding cooling water.
22
Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B 5 MARGINS AND UNCERTAINTIES The uncertainty in the fission power has been implicitly included by calculating the source terms at the nominal maximum fission power plus two standard deviations so the uncertainty in the fission power is not added to the uncertainties in the source term and decay heat calculations.
The remaining uncertainty is the statistical uncertainty in the MCNP6 calculation and the uncertainty in the decay heat.
There is statistical uncertainty in the calculation of the nuclide inventory at discharge (the end of irradiation) as calculated using MCNP6. The average standard deviation in the neutron flux and one-group cross section tallies of the MCNP6 calculations is about** This uncertainty carries
~
~
through to the source terms for the fission gas release.
The uncertainty in the decay heat calculated using ANSl/ANS-5.1-2014 is given by equation 7 of Reference [ANS 2014):
2
= ( b.~d) 2
+ (b.Prod) pd Prod As described in Reference [ANS 2014) the uncertainty in the fission power is implicitly accounted for in the source terms so only !!.Pd is included in the uncertainty and is calculated as described in Reference [ANS 2014) using the following equation:
Nirr tci =t + L (tirr,k + torf.k) k=i+1 The pj and Aj coefficients are obtained from Table 5 of Reference [ANS 2014), and other terms are defined in Section 4.1.5. Figure 4 shows !!.Pd/ Pd as a function of time since the end of irradiation. The uncertainty is about 3% immediately after the end of irradiation, dropping to about 2% for cooling times over 10 seconds.
23
Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B 0.023 0.022 a_"O 0.021
. a_
. "O
Tunk*
shown in Figure 7. Since A. and a are independent of position, the release-to-birth ratio profile of Kr-88 has the same shape as that of 1-131 but is smaller by a factor of:
The net release-to-birth ratio in the rod , (R/B) is calculated by performing a power density weighted integral of R/B over the volume of the U02 pellets in the rod and dividing by the total fission power in the rod as follows:
R f.L R
(!_!_) = 2n f.0 Prdr 0 dz hB = (RP dr (L dz [2nrRph R]
B Prod 10 Rp 10 L ProdL B C-7
Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B This is the equivalent of Equation 2 in Reference [ANS 2011). Here Prod = 2rr f0 Rv rdr f: dz h .
The quantity in the square brackets is called the "weighted release-to-birth ratio." Its value for 1-131 is shown in Figure 10 as a function of position in the hottest rod.
Nuclide= 1311
..c.
""O c::
co I...
[5j f
.0
..c.
C)
"Ci) ca 0:::
Figure 10. Weighted release-to-birth ratio profile for 1-131 in the hottest rod The weighted Rf B is small on the axis of the target rod because r ~ O there . Therefore, most of the 1-131 that is released into the gap gas comes from locations where 0.05 < r < 0.1 cm and 5 <
z < 30 cm. The weighted release-to-birth ratio profile for Kr-88 has the same shape but is a factor of 3.39 smaller.
Integrating the weighted R/B over the radius and length of the U02 pellets in the rod using
=-
release-to-birth ratio for Kr-88 is (R/B) =**I trapezoidal integration gives a net release-to-birth ratio for 1-131 of (R/B) The net a factor of 3.39 smaller. The release-to-birth ratio of 1-131 is larger because its half-life is much larger.
C-8
Attachment 5 Source Term Analysis Design Calculation Report 30441 R00022/B The numbers of atoms of 1-131 and Kr-88 in the gap gas at the end of irradiation are obtained by multiplying the numbers of atoms of 1-131 and Kr-88 in the rod by their respective (R/B) . The numbers of atoms of 1-131 and Kr-88 in the gap gas at later times are found by evolving the nuclide inventory in the gap as described in Section 4.1.3 of this report. For 1-131, this includes the production of 1-131 by the decay of ground state and metastable Te-131 , even though tellurium is not volatile in the rod after the end of irradiation.
C-9
+ t:ENEIUll. ATOMleS P.O. BOX 85608 SAN DIEGO, CA 92186-5608 (858) 455-3000

ML17132A243

Text

SUMMARY

Navigation menu

Search